stack3.m

用于小波变换,调试过很好.可以求出任意多个小波系数.

v(k,1)=A0931; v(k,2)=A1129; v(k,3)=A1327;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,19)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,19)=20*pianyi;         %%%%
tmatr(k,19)=k*detat;%%%%


x=[23*pianyi 19*pianyi 15*pianyi]';%%%%
v(k,1)=A0932; v(k,2)=A1130; v(k,3)=A1328;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,20)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,20)=20.5*pianyi;         %%%%
tmatr(k,20)=k*detat;%%%%


x=[20*pianyi 16*pianyi 12*pianyi]';%%%%
v(k,1)=A1131; v(k,2)=A1329; v(k,3)=A1527;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,21)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,21)=21*pianyi;         %%%%
tmatr(k,21)=k*detat;%%%%

x=[21*pianyi 17*pianyi 13*pianyi]';%%%%
v(k,1)=A1132; v(k,2)=A1330; v(k,3)=A1528;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,22)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,22)=21.5*pianyi;         %%%%
tmatr(k,22)=k*detat;%%%%



x=[22*pianyi 18*pianyi 14*pianyi]';%%%%
v(k,1)=A1133; v(k,2)=A1331; v(k,3)=A1529;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,23)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,23)=22*pianyi;         %%%%
tmatr(k,23)=k*detat;%%%%


x=[23*pianyi 19*pianyi 15*pianyi]';%%%%
v(k,1)=A1134; v(k,2)=A1332; v(k,3)=A1530;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,24)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,24)=22.5*pianyi;         %%%%
tmatr(k,24)=k*detat;%%%%



x=[20*pianyi 16*pianyi 12*pianyi]';%%%%
v(k,1)=A1333; v(k,2)=A1531; v(k,3)=A1729;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,25)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,25)=23*pianyi;         %%%%
tmatr(k,25)=k*detat;%%%%




x=[21*pianyi 17*pianyi 13*pianyi]';%%%%
v(k,1)=A1334; v(k,2)=A1532; v(k,3)=A1730;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,26)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,26)=23.5*pianyi;         %%%%
tmatr(k,26)=k*detat;%%%%



x=[22*pianyi 18*pianyi 14*pianyi]';%%%%
v(k,1)=A1335; 
v(k,2)=A1533; 
v(k,3)=A1731;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,27)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,27)=24*pianyi;         %%%%
tmatr(k,27)=k*detat;%%%%



x=[23*pianyi 19*pianyi 15*pianyi]';%%%%
v(k,1)=A1336; v(k,2)=A1534; v(k,3)=A1732;  %%%%
tmin=x(1)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(1)/2)^2)/(speed*detat));
   v(nnormal,1)=v(k,1);
end 
for k=(nnormal+1):1024
    v(k,1)=0;
end
tmin=x(2)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(2)/2)^2)/(speed*detat));
   v(nnormal,2)=v(k,2);
end 
for k=(nnormal+1):1024
    v(k,2)=0;
end
tmin=x(3)/speed; 
ntmin=round(tmin/detat)+1;
for k=ntmin:1024
   nnormal=round(sqrt((k*detat*speed)^2-(x(3)/2)^2)/(speed*detat));
   v(nnormal,3)=v(k,3);
end 
for k=(nnormal+1):1024
    v(k,3)=0;
end
k=1:1024;
vmatr(k,28)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,28)=24.5*pianyi;         %%%%
tmatr(k,28)=k*detat;%%%%


for m=1:size(vmatr,2)
    for k=1:1024
   if vmatr(k,m)+1>=1
      vmatr(k,m)=0;
   end
    end
end
tmatr=tmatr*speed/2;
pcolor(xmatr,tmatr,vmatr);
axis([0.05 0.12 0 0.24]);
colormap gray  ,shading interp;
set(gca,'XAxislocation','top');
set(gca,'YDir','reverse');