focal=real(ifft2(fft2(sis)*conj(fft2(sis1))));
focal=real(ifft2(fft2(sis)*conj(fft2(sis1'))));
surf(focal,'linestyle','none');axis ij;view([0 90]);axis([1 479 1 999]);colormap('gray')
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
focal=real(ifftshift(ifft2(fftshift(fft2(sis))*conj(fftshift(fft2(sis1))))));
focal=real(ifftshift(ifft2(fftshift(fft2(sis)).*conj(fftshift(fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
focal=real((ifft2((fft2(sis)).*conj((fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
focal=real(ifftshift(ifft2(fftshift(fft2(sis)).*conj(fftshift(fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
focal=real((ifft2(fftshift(fft2(sis)).*conj(fftshift(fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
focal=real((ifft2((fft2(sis)).*conj((fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
focal=real(ifftshift(ifft2((fft2(sis)).*conj((fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
unfocal=real(ifftshift(ifft2((fft2(focal)).*conj((fft2(sis1))))));
surf(focal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
figure;surf(unfocal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
unfocal=real((ifft2((fft2(focal)).*conj((fft2(sis1))))));
figure;surf(unfocal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
unfocal=real((ifft2(fftshift(fft2(focal)).*conj(fftshift(fft2(sis1))))));
figure;surf(unfocal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
unfocal=real(ifftshift(ifft2(fftshift(fft2(focal)).*conj(fftshift(fft2(sis1))))));
figure;surf(unfocal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
unfocal=real(ifftshift(ifft2((fft2(focal)).*conj((fft2(sis1))))));
figure;surf(unfocal,'linestyle','none');axis ij;view([0 90]);colormap('gray')
figure;surf(flipdim(unfocal,1),'linestyle','none');axis ij;view([0 90]);colormap('gray')
figure;surf(sis,'linestyle','none');axis ij;view([0 90]);colormap('gray')
figure;surf(sis1,'linestyle','none');axis ij;view([0 90]);colormap('gray')
ls
save focal
%-- 8/19/07 10:09 PM --%
ls
load day03
figure;surf(h,(0:4:4000),ODCIGlin,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
figure;surf(h,(0:4:4000),ODCIG,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
z(1:1000,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5);z(iz,j)=1.;end;end;
z(1:1000,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5));z(iz,j)=1.;end;end;
figure;surf((-80:1:80),(0:4:4000),temp,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-80 80 0 4000])
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
clear z;z(1:1000,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5));z(iz,j)=1.;end;end;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
clear z;z(1:1000,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5));if iz < 1001;z(iz,j)=1.;end;end;end;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
clear z;z(1:1000,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5));if iz =< 1001;z(iz,j)=1.;end;end;end;
clear z;z(1:1000,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5));if iz <= 1001;z(iz,j)=1.;end;end;end;
clear z;z(1:1001,1:481)=0.;for i=1:161;for j=1:481;iz=floor(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5));if iz <= 1001;z(iz,j)=1.;end;end;end;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
clear z;z(1:1001,1:481)=0.;for i=1:161;for j=1:481;iz=floor(real(sqrt(zm^2-h(j)^2+h(j)*(-4000+(i-1)*50)/4+1.5)));if iz <= 1001;z(iz,j)=1.;end;end;end;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
max(max(z))
zeta
z = zeta;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
z = zeta;
for i=1:161;plot(z);pause;end;
for i=1:161;plot(z(i,:));pause;end;
for i=1:161;plot(z(i,:));axis ij;pause;end;
z = zeta;
zeta=zm^2-h.^2+h*(-4000:50:4000)';
zeta=zm^2-h.^2+h*(-4000:50:4000);
zeta=zm^2-h.^2+h'*(-4000:50:4000);
zeta=zm^2-h.^2'*ones(1,161)+h'*(-4000:50:4000);
for i=1:161;plot(zeta(i,:));axis ij;pause;end;
zeta=sqrt(zm^2-h.^2'*ones(1,161)+h'*(-4000:50:4000));
for i=1:161;plot(zeta(i,:));axis ij;pause;end;
for i=1:161;plot(zeta(i,:));axis ij;axis([-1200 1200 0 4000]);pause;end;
for i=1:161;plot(h,zeta(i,:));axis ij;axis([-1200 1200 0 4000]);pause;end;
for i=1:161;plot(h,real(zeta(i,:)));axis ij;axis([-1200 1200 0 4000]);pause;end;
for i=1:161;plot(h,real(zeta(i,:)'));axis ij;axis([-1200 1200 0 4000]);pause;end;
for i=1:161;plot(h',real(zeta(i,:)'));axis ij;axis([-1200 1200 0 4000]);pause;end;
for i=1:161;plot(h,real(zeta(:,i)'));axis ij;axis([-1200 1200 0 4000]);pause;end;
z = zeta;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
z = zeta;
figure;surf((-1200:5:1200),(0:4:4000),z,'linestyle','none');axis ij;view([0 90]);colormap('gray');axis([-1200 1200 0 4000])
ODCIGenv=conv2(w',z);
figure;surf(ODCIGenv,'linestyle','none');axis ij;view([0 90]);colormap('gray')
z = zeta;
ODCIGenv=conv2(w',z);
figure;surf(ODCIGenv,'linestyle','none');axis ij;view([0 90]);colormap('gray')
z = zeta;
ODCIGenv=conv2(w',z);
figure;surf(ODCIGenv,'linestyle','none');axis ij;view([0 90]);colormap('gray')
figure;surf(ODCIG+ODCIGlin+ODCIGpoint,'linestyle','none');axis ij;view([0 90]);colormap('gray')
%-- 8/20/07  8:26 PM --%
ls
pwd
cd ..
cd Matlab_Lucio\Model\
ls
help raio
help raios
ls
S = Raio(1,1,zeros(1,401),(0:10:4000),ones(1,401)*1000,(0:10:4000),zeros(1,401),(0:10:4000),zeros(1,401),4000,0.01,(0:4::2000),120);
S = Raio(1,1,zeros(1,401),(0:10:4000),ones(1,401)*1000,(0:10:4000),zeros(1,401),(0:10:4000),zeros(1,401),4000,0.01,(0:4:2000),120);
pwd
cd ray
ls
S = Raio(1,1,zeros(1,401),(0:10:4000),ones(1,401)*1000,(0:10:4000),zeros(1,401),(0:10:4000),zeros(1,401),4000,0.01,(0:4:2000),120);
halp raio
halp Raio
help Raio
S = Raio(1,1,zeros(1,401),(0:10:10000),ones(1,401)*1000,(1000:10:4000),zeros(1,401),(0:10:4000),zeros(1,401),4000,0.01,(0:4:2000),120);
S = Raio(1,1,zeros(1,401),(0:10:10000),ones(1,401)*1000,(1000:10:5000),zeros(1,401),(0:10:4000),zeros(1,401),4000,0.01,(0:4:2000),120);
S = Raio(1,1,zeros(1,1001),(0:10:10000),ones(1,1001)*1000,(1000:10:5000),zeros(1,401),(0:10:4000),zeros(1,401),4000,0.01,(0:4:2000),120);
surf(S);view(0 90);axis ij
surf(S);view([0 90]);axis ij
%-- 8/20/07  9:23 PM --%
function sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
surf(sism,'linetype','none');view([0 90]);axis ij
surf(sism,'linestyle','none');view([0 90]);axis ij
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
surf(sism,'linestyle','none');view([0 90]);axis ij
max(max(sis))
max(max(sism))
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
max(max(sism))
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
max(max(sism))
surf(sism,'linestyle','none');view([0 90]);axis ij
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
surf(sism,'linestyle','none');view([0 90]);axis ij
max(max(sism))
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
max(max(sism))
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
surf(sism,'linestyle','none');view([0 90]);axis ij
max(max(sism))
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000));
surf(sism,'linestyle','none');view([0 90]);axis ij
w=[1 1 1];
w=conv(w,w);
w=diff(diff(w));
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
surf(sism,'linestyle','none');view([0 90]);axis ij
for i=1:801;sss=conv(sism(:,i),w);end;
for i=1:801;sss=conv(sism(:,i),w');end;
for i=1:801;sss(:,i)=conv(sism(:,i),w');end;
surf(sss,'linestyle','none');view([0 90]);axis ij
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
clear sss
surf(sism,'linestyle','none');view([0 90]);axis ij
surf(sism,'linestyle','none');view([0 90]);axis ij;colormap('gray')
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
surf(sism,'linestyle','none');view([0 90]);axis ij;colormap('gray')
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
surf(sism,'linestyle','none');view([0 90]);axis ij;colormap('gray')
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
surf(sism,'linestyle','none');view([0 90]);axis ij;colormap('gray')
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
surf(sism,'linestyle','none');view([0 90]);axis ij;colormap('gray')
sism = refcoef(1000,2000,(-2000:10:6000),4000,0.001,(0:4:2000),w);
surf(sism,'linestyle','none');view([0 90]);axis ij;colormap('gray')
%-- 8/20/07 11:06 PM --%
ls
load day03
clear ODCI* adci* impul* temp spike
clear z z01 z02 zm
gama0=26.6;
z=sqrt(2000*cos(pi/180*gama0)-h.^2+2*h*2000*tan(pi/180*(0:1:60)));
z=sqrt(2000*cos(pi/180*gama0)-h.^2+2*h'*2000*tan(pi/180*(0:1:60)));
z=sqrt(2000*cos(pi/180*gama0)-(h').^2*ones(1,61)+2*h'*2000*tan(pi/180*(0:1:60)));
for i=1:61;plot(z(:,i));pause;end;
z=real(sqrt(2000*cos(pi/180*gama0)-(h').^2*ones(1,61)+2*h'*2000*tan(pi/180*(0:1:60))));
for i=1:61;plot(z(:,i));pause;end;
for i=1:61;plot(z(:,i));axis ij;axis([-1200 1200 0 4000]);pause;end;
plot(z)
plot(z);axis ij;axis([-1200 1200 0 4000])
z=real(sqrt(2000+h'*ones(1,61)*tan(pi/180*gama0)-(h').^2*ones(1,61)+2*h'*2000*tan(pi/180*(0:1:60))));
plot(z)
plot(z);axis ij;axis([-1200 1200 0 4000])
%-- 8/26/07  8:46 PM --%
dsr=(-4000:10:4000);
h=(-400:10:400);
zm=1000;
z=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr;
z=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr);
plot(real(z));
dsr=(-4000:50:4000);
clear z
z=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr);
plot(real(z));
plot(real(z));axis ij
plot(h,real(z));axis ij
dsr=(-4000:200:4000);
clear z
z=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr);
plot(h,real(z));axis ij
zcmp=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr);
zshot=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr-ones(size(h)')*dsr.^2);
zshot=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr-ones(size(h'))*dsr.^2);
plot(h,real(zshot));axis ij
dsr=(-4000:50:4000);
zshot=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr-ones(size(h'))*dsr.^2);
plot(h,real(zshot));axis ij
dsr=(-4000:10:4000);
zshot=sqrt(zm^2-(h.^2)'*ones(size(dsr))+h'*dsr-ones(size(h'))*dsr.^2);
plot(h,real(zshot));axis ij
%-- 8/27/07 11:19 PM --%
xs=0
x6=2000;
h=(-600:10:600);
xi(=0:20:4000);
xi=(0:20:4000);
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*x;
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*x)/sqrt(2);
zm=2000;
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*x)/sqrt(2);
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*xi)/sqrt(2);
plot(z)
plot(real(z))
plot(h,real(z));axis ij
xi=(0:100:4000);
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*xi)/sqrt(2);
plot(h,real(z));axis ij
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-(ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*xi)*2)/sqrt(2);
plot(h,real(z));axis ij
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-(ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*xi)*.5)/sqrt(2);
plot(h,real(z));axis ij
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-(ones(size(h'))*xi.^2+2*h'*xi+x6*ones(size(h'))*xi+xs*ones(size(h'))*xi))/sqrt(2);
plot(h,real(z));axis ij
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-2*ones(size(h'))*xi.^2+4*h'*xi+2*x6*ones(size(h'))*xi+2*xs*ones(size(h'))*xi)/sqrt(2);
plot(h,real(z));axis ij
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-2*ones(size(h'))*xi.^2+4*h'*xi+2*x6*ones(size(h'))*xi+2*xs*ones(size(h'))*xi)/sqrt(2)-250;
plot(h,real(z));axis ij
atan(540/1200)*180/pi
540/1200
tan(0.45)*1000
tan(0.45)*2000
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-0.5*ones(size(h'))*xi.^2+h'*xi+0.5*x6*ones(size(h'))*xi+0.5*xs*ones(size(h'))*xi)/sqrt(2);
plot(h,real(z));axis ij
z=sqrt(2*zm^2-xs*x6-2*xs*h'*ones(size(xi))-2*(h.^2)'*ones(size(xi))-4*ones(size(h'))*xi.^2+8*h'*xi+4*x6*ones(size(h'))*xi+4*xs*ones(size(h'))*xi)/sqrt(2);
plot(h,real(z));axis ij
%-- 8/28/07 12:28 AM --%
xi=(0:100:4000);
h=(-600:10:600);
xm=1000;zm=2000;xr=2000;xs=0;
z=sqrt(zm^2-2*xm*xr-xm^2-ones(size(h'))*xi.^2+(sqrt(2*zm^2+(ones(size(h'))*xi-xs).^2)-(xm-xr)).^2);
plot(real(z))
z=sqrt(-(h.^2)'*ones(size(xi))+2*(xm-h')*xs*ones(size(xi))+2*(xm+h')*xi+zm^2-2*xm*xs-xm^2-ones(size(h'))*xi.^2+(sqrt(2*zm^2+(ones(size(h'))*xi-xs).^2)-(xm-xr)).^2);
plot(real(z))
z=sqrt(-(h.^2)'*ones(size(xi))+2*(xm-h')*xs*ones(size(xi))+2*(xm+h')*xi+zm^2-2*xm*xs-xm^2-ones(size(h'))*xi.^2+(sqrt(2*zm^2+(ones(size(h'))*xi-xs).^2)-(xm-xr)).^2)/sqrt(2);
plot(real(z))
z=sqrt(-2*(h.^2)'*ones(size(xi))+2*(xm-h')*xs*ones(size(xi))+2*(xm+h')*xi+zm^2-2*xm*xs-xm^2-ones(size(h'))*xi.^2+(sqrt(2*zm^2+(ones(size(h'))*xi-xs).^2)-(xm-xr)).^2)/sqrt(2);
plot(real(z))
z=sqrt(-2*(h.^2)'*ones(size(xi))-2*h'*xs*ones(size(xi))+2*(xm+h')*xi+zm^2-xm^2+(sqrt(2*zm^2+ones(size(h'))*(xi-xs).^2)-(xm-xr)).^2)/sqrt(2);
plot(real(z))
z=sqrt(-4*(h.^2)'*ones(size(xi))-2*h'*xs*ones(size(xi))+2*(xm+h')*xi+zm^2-xm^2+(sqrt(2*zm^2+ones(size(h'))*(xi-xs).^2)-(xm-xr)).^2)/sqrt(2);
plot(real(z))
z=sqrt(-4*(h.^2)'*ones(size(xi))-2*h'*xs*ones(size(xi))+2*(xm+h')*xi+zm^2-xm^2+(sqrt(4*zm^2+ones(size(h'))*(xi-xs).^2)-sqrt(zm^2+(xm-xr))).^2)/sqrt(2);
plot(real(z))
%-- 10/03/07  9:29 PM --%
ls
cd Chapter08
ls
help britas
help meteroros
help meteoros
help pedras
load modelos
clean
clear
x=(0:25:4000);
x=(0:50:4000);
t=(0:4:1000);
sism=zeros(251,81);for i=1:81;j=floor(((250/4000)^2+(x(i)/4000)^2)/4);sism(j,i)=1;end;
sism=zeros(251,81);for i=1:81;j=floor(sqrt((250/4000)^2+(x(i)/4000)^2)/4);sism(j,i)=1;end;
sism=zeros(251,81);for i=1:81;j=floor(sqrt((250/4000)^2+(x(i)/4000)^2)/4);disp(j);end;
sism=zeros(251,81);for i=1:81;j=floor(sqrt((250/4000)^2+(x(i)/4000)^2)/0.004);disp(j);end;
sism=zeros(251,81);for i=1:81;j=floor(sqrt((250/4000)^2+(x(i)/4000)^2)/0.004);sism(j,i)=1;end;
imagesc(sism)
sism=zeros(251,81);for i=1:81;j=floor(sqrt((250/8000)^2+(x(i)/8000)^2)/0.004);sism(j,i)=1;end;
imagesc(sism)
sism=zeros(251,81);for i=1:81;j=floor(sqrt((500/8000)^2+(x(i)/8000)^2)/0.004);sism(j,i)=1;end;
imagesc(sism)
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;
imagesc(sism)
x=(0:25:2000);
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;
x=(0:25:2000);
imagesc(sism)
w=[1 1 1]
w=conv(w,w);
imagesc(sism)
for i=1:81;sismo(:,i)=conv(w,sism(:,i));end;
imagesc(sismo)
sism=sismo(16:266,:);
imagesc(sism)
sism=sismo(17:267,:);
imagesc(sism)
help rec_down
u = rec_down(x,0:10:700,t,ones(81,71)*4000,sism)
u = rec_down(x,0:10:700,t,ones(81,71)*4000,sism);
imagesc(u(:,:,1))
u = rec_down(x,0:10:700,t,ones(81,71)*4000,sism);
imagesc(u(:,:,1));
imagesc(sism);
max(u)
max(max(u))
max(max(max(u)))
imagesc(u(:,:,11));
for i=1:251;imagesc(u(:,:,i));pause(0.01);end;
imagesc(sism);
u = rec_down(x,0:10:700,t,ones(81,71)*4000,sism);
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
imagesc(u(:,:,1));
for i=1:251;imagesc(u(:,:,i));pause(0.01);end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
clear u
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
imagesc(u(:,:,1));
for i=1:251;imagesc(u(:,:,i));pause(0.01);end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=1:251;imagesc(u(:,:,i));pause(0.01);end;
for i=1:251;imagesc(u(i,:,:));pause(0.01);end;
for i=1:251;a(i,:,:)=u(:,:,i);end;imagesc(a(:,:,1))
for i=1:251;for j=1:81;a(i,j)=u(1,j,i);end;imagesc(a)
for i=1:251;for j=1:81;a(i,j)=u(1,j,i);end;end;imagesc(a)
clear a
for i=1:251;for j=1:81;a(i,j)=u(1,j,i);end;end;imagesc(a)
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for a=11:-1:1;disp(a);end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
imagesc(u(:,:,1));
imagesc(u(:,:,2));
imagesc(u(:,:,20));
for i=1:251;imagesc(u(:,:,i));pause(0.01);end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=1:251;imagesc(u(:,:,i));pause(0.01);end;
for i=1:251;imagesc(u(:,:,i));pause;end;
transpose
help transpose
for i=1:251;imagesc(u(:,:,i)');pause;end;
for i=1:251;imagesc(u(:,:,i));pause;end;
u = rec_down(x,0:10:700,t,ones(71,81)*4000,sism);
for i=1:251;imagesc(u(:,:,i));pause;end;
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
