for i=251:-1:1;imagesc(u(:,:,i));pause;end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
imagesc(u(:,:,1));
for i=251:-1:1;imagesc(u(:,:,i));pause;end;
for i=251:-1:1;imagesc(u(:,:,i));pause(0.1);end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=251:-1:1;imagesc(u(:,:,i));pause(0.1);end;
imagesc(sism);
for i=251:-1:1;imagesc(u(:,:,i));disppause(0.1);end;
for i=251:-1:1;imagesc(u(:,:,i));disp(i)pause(0.01);end;
for i=251:-1:1;imagesc(u(:,:,i));disp(i);pause(0.01);end;
for i=251:-1:1;imagesc(u(:,:,i));disp(i);pause(0.1);end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=251:-1:1;imagesc(u(:,:,i));disp(i);pause(0.1);end;
x=(0:5:400);
t=(0:0.5:250);
sism=zeros(501,81);for i=1:81;j=floor(sqrt((500/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
imagesc(sism);
t=(0:0.5:100);
sism=zeros(201,81);for i=1:81;j=floor(sqrt((500/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
imagesc(sism);
t=(0:0.5:50);
sism=zeros(201,81);for i=1:81;j=floor(sqrt((500/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
imagesc(sism);
x(81)
for i=1:81;sismo(:,i)=conv(w,sism(:,i));end;
clear sismo;for i=1:81;sismo(:,i)=conv(w,sism(:,i));end;
imagesc(sismo);
clear sism;sism=sismo(17:217,:);
imagesc(sism);
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=251:-1:1;imagesc(u(:,:,i));disp(i);pause(0.1);end;
for i=101:-1:1;imagesc(u(:,:,i));disp(i);pause(0.1);end;
%-- 10/04/07 12:18 AM --%
t=(0:0.5:50);
x=(0:5:400);
sism=zeros(101,81);for i=1:81;j=floor(sqrt((500/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
for i=1:81;sismo(:,i)=conv(w,sism(:,i));end;
w=[1 1 1]
w=conv(w,w);
for i=1:81;sismo(:,i)=conv(w,sism(:,i));end;
clear sism;sism=sismo(17:117,:);
imagesc(sism);
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=1:71;for j=1:81;for k=1:101;a(j,k,i)=u(i,j,k);end;end;end;
for i=1:71;imagesc(a(:,:,i));disp(i);pause(0.1);end;
for i=1:71;imagesc(a(:,:,i));disp(i);pause;end;
for i=1:71;imagesc(a(:,:,i)');disp(i);pause;end;
u = rec_down(x,0:5:1000,t,ones(71,81)*4000,sism);
u = rec_down(x,0:5:1000,t,ones(201,81)*4000,sism);
for i=1:71;imagesc(a(:,:,i)');disp(i);pause;end;
clear u sismo a i j k
save rtm01
%-- 10/04/07  7:58 PM --%
cd ../RTM/
ls
load rtm02
cd ../Anali
pwd
cd ../../AnaliseNumerica/
load rtm01
cd ../Doutorado/RTM/
%-- 10/05/07 12:58 PM --%
rho=2500;g=10;z=-1000;
sigma=[2*rho*g*z/5 rho*g*z/50 0;rho*g*z/50 2*rho*g*z/5 0;0 0 rho*g*z]
help eig
[V,D]=eig(sigma)
[V,D]=eig(sigma/1000000000)
sigma/1000000000
10^9
1000000000
[V,D]=eig(sigma/1000000000)
[V,D]=eig(-sigma/1000000000)
clear
t=(0:0.0005:.5);
t1=(-.25:0.0005:.25);
w1=(1-2*pi*(pi*30*t1).^2).*exp(-pi*(pi*30*t1).^2);
sism=zeros(1001,161);for i=1:161;j=floor(((sqrt((250/4000)^2+(((i-1)*5-400)/4000)^2)))/0.0005);sism(j,i)=1;end;
clear sism1;for i=1:161;sism1(:,i)=conv(w1,sism(:,i));end;
clear sism2;sism2=sism1(501:1501,:);
clear sism3;sism3=sism2(1:300,:);
imagesc(sism3)
125*0.0005*2000
u=rec_down(0:5:800,0:5:1000,(0:0.0005:.1495),ones(201,161)*4000,sism3);
pwd
u=rec_down(0:5:800,0:5:1000,(0:0.0005:.1495),ones(201,161)*4000,sism3);
for i=300:-1:1;imagesc(u(:,:,i));pause(0.1);end;
pwd
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
u=rec_down(0:5:800,0:5:1000,(0:0.0005:.1495),ones(201,161)*2000,sism3);
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
u=rec_down(0:5:800,0:5:1000,(0:0.0005:.1495),ones(201,161)*4000,sism3);
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
diff
help diff
u=rec_down(0:5:800,0:5:500,(0:0.0005:.1495),ones(201,161)*4000,sism3);
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
plot(u(:,81,1))
u=rec_down(0:5:800,0:5:500,(0:0.0005:.1495),ones(201,161)*4000,sism3);
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
plot(u(:,81,1))
u=rec_down(0:5:800,0:5:500,(0:0.0005:.1495),ones(201,161)*4000,sism3);
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
u=rec_down(0:5:800,0:5:500,(0:0.0005:.1495),ones(201,161)*4000,sism3);
for i=297:-4:1;imagesc(u(:,:,i));pause(0.1);end;
pwd
ls
help we_df_O4
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.1495,ones(201,161)*4000,wav
source=sism3(:,81);
ww=zeros(300,1);
ww(1:100)=source(126:225);
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.1495,ones(201,161)*4000,ww);
for i=1:4:300;imagesc(us(:,:,i));pause(0.1);end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.1495,ones(101,161)*4000,ww);
for i=1:4:300;imagesc(us(:,:,i));pause(0.1);end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(101,161)*4000,ww);
clear us;us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(101,161)*4000,ww);
clear us
clear us;us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(101,161)*4000,ww);
for i=1:4:300;imagesc(us(:,:,i));pause(0.1);end;
clear us;us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(101,161)*4000,ww(1:150));
for i=1:4:300;imagesc(us(:,:,i));pause(0.1);end;
for i=1:4:150;imagesc(us(:,:,i));pause(0.1);end;
for i=1:4:149;imagesc(us(:,:,i));pause(0.1);end;
for i=1:4:141;imagesc(us(:,:,i));pause(0.1);end;
us=we_df_O4(0:10:800,0:10:500,0:0.0005:.07475,ones(101,161)*4000,ww);
for i=1:4:141;imagesc(us(:,:,i));pause(0.1);end;
us=we_df_O4(0:2:800,0:2:500,0:0.0005:.07475,ones(101,161)*4000,ww);
us=we_df_O4(0:2:800,0:2:500,0:0.0005:.07475,ones(301,461)*4000,ww);
for i=1:4:141;imagesc(us(:,:,i));pause(0.1);end;
clear us;us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(101,161)*4000,ww(1:150));
for i=1:4:141;imagesc(us(:,:,i));pause(0.1);end;
plot(us(:,81,144))
save rtm03
us=we_df_O4(0:2:800,0:2:500,0:0.0005:.07475,ones(301,461)*4000,ww);
save rtm03
%-- 10/06/07 12:53 PM --%
ls
ls rtm*
load rtm03
us=we_df_O4(0:10:800,0:10:500,0:0.0005:.07475,ones(301,461)*4000,ww);
for i=1:4:141;imagesc(us(:,:,i));pause(0.1);end;
for i=1:4:141;imagesc(us(:,:,i));pause;end;
us=we_df_O4(0:1:800,0:1:500,0:0.0005:.07475,ones(601,921)*4000,ww);
clear u us
clear sism sism2 sism1
us=we_df_O4(0:1:800,0:1:500,0:0.0005:.07475,ones(601,921)*4000,ww);
us=we_df_O4(0:2:800,0:2:500,0:0.0005:.07475,ones(601,921)*4000,ww);
for i=1:4:300;imagesc(us(:,:,i));pause;end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww);
for i=1:4:300;imagesc(us(:,:,i));pause;end;
us=we_df_O4(0:6:800,0:6:500,0:0.0005:.07475,ones(601,921)*4000,ww);
for i=1:4:300;imagesc(us(:,:,i));pause;end;
0.07475/0.0005
us=we_df_O4(0:6:800,0:6:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:4:151;imagesc(us(:,:,i));pause;end;
us=we_df_O4(0:6:800,0:6:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:4:151;imagesc(us(:,:,i));pause;end;
plot(us(:,63,151))
help zeros
help green
src=zeros(84,134,151);
src(1,67,2:151)=diff(ww);
src(1,67,2:151)=diff(ww');
src(1,67,2:151)=diff(ww(1:151));
g=green(ones(601,921)*4000,[400 0],src,5,0.0005,150,3);
g=green(ones(151,231)*4000,[400 0],src,5,0.0005,150,3);
g=green(ones(151,231)*4000,[400 134],src,5,0.0005,150,3);
g=green(ones(84,134)*4000,[400 134],src,5,0.0005,150,3);
g=green(ones(84,134)*4000,[400 68],src,5,0.0005,150,3);
for i=1:150;imagesc(g);pause;end;
for i=1:150;imagesc(g(:,:,i));pause;end;
for i=1:150;imagesc(us(:,:,i));pause;end;
for i=1:150;imagesc(us(:,:,i));colormap('gray')pause;end;
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us1=we_df_O4(0:6:800,0:6:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us1(:,:,i));colormap('gray');pause;end;
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
for i=1:150;imagesc(us(:,:,50));colormap('gray');pause;end;
imagesc(us(:,:,80))
fftshift
hrelp fftshift
help fftshift
help fft2
imagesc(abs(fft2(us(:,:,80))))
imagesc(abs(fftshift(fft2(us(:,:,80)))))
imagesc(abs(fftshift(fft2([1 1 1;1 4 1;1 1 1]))))
imagesc(real(ifft2(fft2([1 1 1;1 4 1;1 1 1])*fft2(us(:,:,80)))))
imagesc(real(ifft2(fft2([1 1 1;1 4 1;1 1 1]).*fft2(us(:,:,80)))))
filt=zeros(84,134);
filt(41,66)=1;filt(42,66)=1;filt(43,66)=1;filt(41,67)=1;filt(42,67)=4;filt(43,67)=1;filt(41,68)=1;filt(42,68)=1;filt(43,68)=1;
imagesc(real(ifft2(fft2(filt).*fft2(us(:,:,80)))))
imagesc(real(ifft2(fftshift(fft2(filt).*fft2(us(:,:,80))))))
imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
figure;imagesc(us(:,:,80));
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;imagesc(us(:,:,80));
filt(50,80)=1;filt(51,80)=1;filt(52,80)=1;filt(50,81)=1;filt(51,81)=4;filt(51,81)=1;filt(50,82)=1;filt(51,82)=1;filt(52,82)=1;
filt=zeros(101,161);
filt(50,80)=1;filt(51,80)=1;filt(52,80)=1;filt(50,81)=1;filt(51,81)=4;filt(51,81)=1;filt(50,82)=1;filt(51,82)=1;filt(52,82)=1;
imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
figure;imagesc(us(:,:,80));
figure;imagesc(us(:,:,80),real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))));
imagesc(filt)
filt(50,80)=1;filt(51,80)=1;filt(52,80)=1;filt(50,81)=1;filt(51,81)=4;filt(52,81)=1;filt(50,82)=1;filt(51,82)=1;filt(52,82)=1;
imagesc(filt)
figure;imagesc(us(:,:,80),real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))));
figure;imagesc(us(:,:,80));
figure;imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
imagesc(filt)
imagesc(real(fftshift(fft2(filt))))
filt(50,80)=0;filt(51,80)=1;filt(52,80)=0;filt(50,81)=1;filt(51,81)=4;filt(52,81)=1;filt(50,82)=0;filt(51,82)=1;filt(52,82)=0;
imagesc(real(fftshift(fft2(filt))))
figure;imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
figure;imagesc(us(:,:,80));
filt1=zeros(101,161);
filt(51,81)=1;filt(52,81)=1;filt(53,81)=1;filt(51,82)=1;filt(52,82)=4;filt(53,82)=1;filt(51,83)=1;filt(52,83)=1;filt(53,83)=1;
figure;imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
imagesc(filt)
filt(51,80)=0;filt(50,81)=0;
figure;imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
figure;imagesc(us(:,:,80));
imagesc(filt)
filt(51,81)=0;filt(53,81)=0;filt(51,83)=0;filt(53,83)=0;
figure;imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
figure;imagesc(us(:,:,80));
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
imagesc(us)
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
imagesc(us)
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
clear us1
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;imagesc(us(:,:,80));
figure;imagesc(real(fftshift(ifft2(fft2(filt).*fft2(us(:,:,80))))))
figure;imagesc(filt);
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
filt=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
figure;imagesc(filt);
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:10:800,0:10:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:10:800,0:10:500,0:0.001:.07475,ones(601,921)*4000,ww(1:150));
for i=1:150;imagesc(us(:,:,i));colormap('gray');pause;end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:151;imagesc(us(:,:,i));colormap('gray');pause(0.1);end;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:151;imagesc(us(:,:,i));colormap('gray');pause(0.1);end;
imagesc(us(75,:,:))
for i=1:161;for j=1:151;sss(j,i)=us(75,i,j);end;end;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(25,i,j);end;end;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(15,i,j);end;end;figure;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(5,i,j);end;end;figure;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(1,i,j);end;end;figure;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(2,i,j);end;end;figure;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(3,i,j);end;end;figure;imagesc(sss)
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:161;for j=1:151;sss(j,i)=us(3,i,j);end;end;figure;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(13,i,j);end;end;figure;imagesc(sss)
plot(sss:,81)
plot(sss,81)
plot(sss(:,81))
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for i=1:161;for j=1:151;sss(j,i)=us(13,i,j);end;end;figure;plot(sss(:,81))
clear srd
clear src
clear g
clear filt*
help rec_down
ur=rec_down(0:5:800,0:5:500,0:0.0005:.1495,ones(601,921)*4000,sism3);
for i=1:4:300;imagesc(ur(:,:,i));colormap('gray');pause(0.1);end;
for i=1:161;for j=1:151;sss(j,i)=ur(13,i,j);end;end;figure;plot(sss(:,81))
for i=1:161;for j=1:151;sss(j,i)=ur(13,i,j);end;end;figure;imagesc(sss)
for i=1:161;for j=1:151;sss(j,i)=us(13,i,j);end;end;figure;imagesc(sss)
for k=1:151;for i=1:161;for j=1:101;sss(j,i)=us(j,i,k);end;end;imagesc(sss);end;
for k=1:151;for i=1:161;for j=1:101;sss(j,i)=us(j,i,k);end;end;imagesc(sss);pause;end;
for k=1:151;imagesc(us(:,:,k);end;
for k=1:151;imagesc(us(:,:,k));end;
for k=1:151;figure;imagesc(us(:,:,k));figure;imagesc(ur(:,:,151-k));pauseend;
for k=1:151;figure;imagesc(us(:,:,k));figure;imagesc(ur(:,:,151-k));pause;end;
for k=1:4:151;figure;imagesc(us(:,:,k));figure;imagesc(ur(:,:,300-k));pause;end;
help close
close all
for k=1:4:151;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(ur(:,:,300-k));pause;end;
for k=1:4:151;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(ur(:,:,151-k));pause;end;
for k=1:4:151;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(ur(:,:,k));pause;end;
for i=1:161;for j=1:101; ws(j,i,:)=fftshift(fft(us(j,i,:)));end;end;
for i=1:161;for j=1:101; wr(j,i,:)=fftshift(fft(ur(j,i,1:151)));end;end;
for i=1:151;img(:,:)=img(:,:)+conj(ws(:,:,i)).*wr(:,:,i);end;
img=zeros(101,161);for i=1:151;img(:,:)=img(:,:)+conj(ws(:,:,i)).*wr(:,:,i);end;
imgsc(img)
imagesc(img)
img=zeros(101,161);for i=1:151;img(:,:)=img(:,:)+conj(ws(:,:,i)).*wr(:,:,i);end;img=real(img);
imagesc(img)
img1=zeros(101,321);for i=1:161;img1(:,i:i+160)=img1(:,i:i+160):img;end;
imagesc(img1)
img1=zeros(101,321);for i=1:161;img1(:,i:i+160)=img1(:,i:i+160)+img;end;
imagesc(img1)
plot(img1(:,161))
imagesc(img1)
imagesc(img1);colormap('gray')
image=rtm(0:5:800,0:5:500,0:0.0005:.1495,ones(601,921)*4000,sism3,ww);
imagesc(image)
sism=zeros(251,81);for i=1:81;j=floor(sqrt((500/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
clear sism
sismdip30zeros(251,161);for i=1:161;j=floor(sqrt((500/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
x=(-400:5:400);
sismdip30=zeros(251,161);for i=1:161;j=floor(sqrt((x(i)^2+(2*500*sind(30))^2+4*500*x(i)*sind(30)^2)/4000)/0.004);sism(j,i)=1;end;
sind=sin(pi*30/180)
sismdip30=zeros(251,161);for i=1:161;j=floor(sqrt((x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.004);sism(j,i)=1;end;
imagesc(sismdip30)
sismdip30=zeros(251,161);for i=1:161;j=floor(sqrt((x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);sism(j,i)=1;end;
sismdip30=zeros(251,161);for i=1:161;j=floor(sqrt((x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sism(j,i)=1;end;
sismdip30=zeros(251,161);for i=1:161;j=floor(sqrt((x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000));disp(num2str(j));sism(j,i)=1;end;
sismdip30=zeros(251,161);for i=1:161;j=floor((sqrt((x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2))/4000)/0.0005);disp(num2str(j));sism(j,i)=1;end;
sismdip30=zeros(251,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sism(j,i)=1;end;
sismdip30=zeros(251,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000));disp(num2str(j));sism(j,i)=1;end;
sismdip30=zeros(251,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sism(j,i)=1;end;
sismdip30=zeros(400,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sism(j,i)=1;end;
imagesc(sismdip30)
sismdip30=zeros(400,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sismdip30(j,i)=1;end;
imagesc(sismdip30)
plot(ww)
plot(w1)
for i=1:161;sism1(:,i)=conv(w1,sismdip30(:,i));end;
sism1=zeros(2000,161);for i=1:161;sism1(:,i)=conv(w1,sismdip30(:,i));end;
sism1=zeros(1400,161);for i=1:161;sism1(:,i)=conv(w1,sismdip30(:,i));end;
a=conv(w1,sismdip30(:,1));
plot(a)
sism1=zeros(1400,161);for i=1:161;sism1(:,i)=conv(w1',sismdip30(:,i));end;
imagesc(sism1)
sims2=sism1(501:900,:);
imagesc(sism2)
imagesc(sims2)
figure;imagesc(sismdip30)
sismdip30=zeros(400,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sismdip30(j,i)=1/j;end;
sism1=zeros(1400,161);for i=1:161;sism1(:,i)=conv(w1',sismdip30(:,i));end;
sims2=sism1(501:900,:);
imagesc(sims2)
imagesc(sims2);colormap('gray')
sismdip30=zeros(450,161);for i=1:161;j=floor((sqrt(x(i)^2+(2*500*sind)^2+4*500*x(i)*sind^2)/4000)/0.0005);disp(num2str(j));sismdip30(j,i)=1/j;end;
sism1=zeros(1450,161);for i=1:161;sism1(:,i)=conv(w1',sismdip30(:,i));end;
sims2=sism1(501:950,:);
imagesc(sims2);colormap('gray')
image=rtm(0:5:800,0:5:700,0:0.0005:.2195,ones(601,921)*4000,sism2,ww);
image=rtm(0:5:800,0:5:700,0:0.0005:.2195,ones(601,921)*4000,sims2,ww);
imagesc(image);colormap('gray')
15/25
tan(30*pi/180)
15/26
image=rtm(0:5:800,0:5:700,0:0.0005:.2195,ones(601,921)*4000,sims2,ww);
imagesc(image);colormap('gray')
gaussiana=zeros(101,101);
gaussiana(50,50)=1;
gaussiana(50,51)=1;
gaussiana(50,52)=1;
gaussiana(51,52)=1;
gaussiana(51,51)=1;
gaussiana(51,50)=1;
gaussiana(52,50)=1;
gaussiana(52,51)=1;
gaussiana(52,52)=1;
imagesc(gaussiana)
gauss2=real(fftshift(ifft2(fft2(gaussiana).^2)));
imagesc(gauss2)
gauss2=real(fftshift(ifft2(fft2(gauss2).^2)));
imagesc(gauss2)
gauss2=real(fftshift(ifft2(fft2(gauss2).^2)));
imagesc(gauss2)
gauss2=real(fftshift(ifft2(fft2(gaussiana).^2)));
gauss2=real(fftshift(ifft2(fft2(gauss2).^2)));imagesc(gauss2);
surf(gauss2)
maxloc(gauss2)
imagesc(gauss2)
gauss3=zeros(101,101);gauss3(32:70,32:70)=gauss2(1:39,1:39);
imagesc(gauss3)
gauss3=zeros(101,101);gauss3(32:70,32:70)=gauss2(1:39,1:39);gauss3=gauss3/max(max(gauss3));
imagesc(gauss3)
velgauss=5000-1000*gauss3;
imagesc(velgauss)
velgauss1=ones(1001,1001)*5000;velgauss1(450:550,450:550)=velgauss;
imagesc(velgauss1)
velgauss1=ones(1001,1001)*5000;velgauss1(1:101,450:550)=velgauss;
imagesc(velgauss1)
velgauss1(100,:)=4000;
us=we_df_O4(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150)
us1 = sht_down_mod(0:5:800,0:5:700,0:0.0005:.2195,velgauss1,ww(1:150));
for i=1:4:440;imagesc(us1(:,:,i));colormap('gray');pause;end;
us1 = sht_down_mod(0:5:800,0:5:700,0:0.0005:.2195,velgauss1,ww(1:150));
imagesc(velgauss1)
us1 = sht_down_mod(0:5:800,0:5:700,0:0.0005:.2195,velgauss1,ww(1:150));
for i=1:4:440;imagesc(us1(:,:,i));colormap('gray');pause;end;
for i=1:4:440;imagesc(us1(:,:,i));colormap('gray');pause(0.01);end;
us1 = sht_down_mod(0:5:800,0:5:700,0:0.0005:.2195,velgauss1,ww(1:150));
imagesc(velgauss(:,423:579))
imagesc(velgauss1(:,423:579))
us1 = sht_down_mod(0:5:800,0:5:700,0:0.0005:.2195,velgauss1,ww(1:150));
for i=1:4:440;imagesc(us1(:,:,i));colormap('gray');pause(0.01);end;
for i=1:4:440;imagesc(us1(:,:,i));colormap('gray');pause(0.1);end;
clear j k sind sism sism1 sismdip30 sss ur us velgauss ws wr
clear a ans gauss* i img image img1
clear us1
save rtm03
prefdir
%-- 10/08/07  8:15 PM --%
ls
load rtm03
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150)
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));end;
for k=1:150;imagesc(us(:,:,k));pause(0.01);end;
max(us)
max(max(max(us)))
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.01);end;
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.01);end;
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.01);end;
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.01);end;
for k=1:150;imagesc(us(:,:,k));pause(0.1);end;
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.1);end;
us=sht_down_2t10x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.1);end;
us=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us(:,:,k));pause(0.1);end;
for k=1:150;imagesc(us(:,:,k));pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
us=sht_down_2t10x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us=sht_down(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
for k=1:150;surf(us1(:,:,k),linestyle,'none');pause(0.1);colormap('gray');end;
for k=1:150;surf(us1(:,:,k),'linestyle','none');pause(0.1);colormap('gray');end;
for k=1:150;surf(us1(:,:,k),'linestyle','none');view(141,66);pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:5:800,0:5:500,0:0.001:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;surf(us1(:,:,k),'linestyle','none');view(141,66);pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:10:800,0:10:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;surf(us1(:,:,k),'linestyle','none');view(141,66);pause(0.1);colormap('gray');end;
for k=1:150;imagesc(us1(:,:,k);pause(0.1);colormap('gray');end;
for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:10:800,0:10:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
figure;for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:2:800,0:2:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:5:800,0:5:500,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us1=sht_down_2t8x(0:5:1600,0:5:1000,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;for k=1:150;imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
us=sht_down_2t8x(0:5:1600,0:5:1000,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
for k=1:150;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
figure;for k=1:150;subplot(2,1,1);imagesc(us(:,:,k));subplot(2,1,2);imagesc(us1(:,:,k));pause(0.1);colormap('gray');end;
figure;imagesc(us(:,:,150)/max(max(us(:,:,150)))-us1(:,:,150)/max(max(us(:,:,150))));pause(0.1);colormap('gray');end;
figure;imagesc(us(:,:,150)/max(max(us(:,:,150)))-us1(:,:,150)/max(max(us(:,:,150))));pause(0.1);colormap('gray');
figure;a=(us(:,:,150)/max(max(us(:,:,150)))-us1(:,:,150)/max(max(us(:,:,150))));imagesc(a);pause(0.1);colormap('gray');
figure;imagesc(a)
figure;a=(us(:,:,150)/max(max(us(:,:,150)))-us(:,:,150)/max(max(us(:,:,150))));imagesc(a);pause(0.1);colormap('gray');
max(max(us(:,:,150)))
max(max(us1(:,:,150)))
figure;a=(us(:,:,150)/max(max(us(:,:,150)))-us1(:,:,150)/max(max(us1(:,:,150))));imagesc(a);pause(0.1);colormap('gray');
figure;imagesc(us)
figure;imagesc(us(:,:,150))
figure;imagesc(us1(:,:,150))
us=sht_down_2t8x(0:5:1600,0:5:1000,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;imagesc(us(:,:,150))
figure;imagesc(us(:,:,150));colormap('gray')
[A C D] = solv(2);
[A,C,D] = solv(2);
[A,C,D] = solv(2)
[A,C,D] = solv(1)
c/2.5
C/2.5
C/5
D/5
B/5
A/5
A+1+C+D
(A+1+C+D)/5
(A+1+C+D)/2.5
C/5
D/5
A/5
(A+1+C+D)/2.5
us=sht_down_2t8x(0:5:1600,0:5:1000,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;imagesc(us(:,:,150));colormap('gray')
us=sht_down_2t8x(0:5:1600,0:5:1000,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;imagesc(us(:,:,150));colormap('gray')
us1=sht_down_2t8x(0:5:1600,0:5:1000,0:0.0005:.07475,ones(601,921)*4000,ww(1:150));
figure;imagesc(us1(:,:,150));colormap('gray')
