Layer With Value density or velocity

Program ini berisi menampilkan lapisan 2 dimensi yang mempunyai nilai densitas atau kecepatan // this program contain showing a 2-D layer that have a value density or velocity

Program ini prototipe untuk tomografi // This is Prototipe for tomography

Pada penjalanan program kita diminta untuk menuliskan koordinat yang akan di chek atau dicari kecepatannya, hasilnya ditampilkan di layar

script :

% PROGRAM LAYER TOMOGRAFI

clear;

clc;

a=0:100;b=0:100;

[x,y]=meshgrid(a,b);

plot(x,y,’.white’),hold on % horee..coba coba kata “white” eh bisa

plot(0,10,’*r’);plot(0,30,’*r’);plot(0,60,’*r’);plot(0,90,’*r’),hold on;

% gambar template

a=0:100;b=80:100;

[x,y]=meshgrid(a,b);

plot(x,y,’.b’),hold on

a=50:100;b=60:80;

[x,y]=meshgrid(a,b);

plot(x,y,’.b’),hold on

a=0:50;b=60:80;

[x,y]=meshgrid(a,b);

plot(x,y,’.g’)

a=0:100;b=40:60;

[x,y]=meshgrid(a,b);

plot(x,y,’.g’)

a=0:100;b=20:40;

[x,y]=meshgrid(a,b);

plot(x,y,’.c’)

a=0:100;b=0:20;

[x,y]=meshgrid(a,b);

plot(x,y,’.m’)

a=input(‘masukkan nilai a ‘)

b=input(‘masukkan nilai b ‘)

% inisiasi kecepatan tiap lapisan

if (b<=20)

v=3500;

elseif b<=40

v=3000;

elseif ((b<=80)&(a<=50))||((b<=60)&(a>50))

v=2500;

else

v=2000;

end

disp(['maka kecepatan adalah =' num2str(v)])