stack3.m

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

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,8)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,8)=14.5*pianyi;         %%%%
tmatr(k,8)=k*detat;%%%%

x=[20*pianyi 16*pianyi 12*pianyi]';%%%%
v(k,1)=A0525; v(k,2)=A0723; v(k,3)=A0921;  %%%%
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,9)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,9)=15*pianyi;         %%%%
tmatr(k,9)=k*detat;%%%%


x=[21*pianyi 17*pianyi 13*pianyi]';%%%%
v(k,1)=A0526; v(k,2)=A0724; v(k,3)=A0922;  %%%%
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,10)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,10)=15.5*pianyi;         %%%%
tmatr(k,10)=k*detat;%%%%

x=[22*pianyi 18*pianyi 14*pianyi]';%%%%
v(k,1)=A0527; v(k,2)=A0725; v(k,3)=A0923;  %%%%
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,11)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,11)=16*pianyi;         %%%%
tmatr(k,11)=k*detat;%%%%


x=[23*pianyi 19*pianyi 15*pianyi]';%%%%
v(k,1)=A0528; v(k,2)=A0726; v(k,3)=A0924;  %%%%
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,12)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,12)=16.5*pianyi;         %%%%
tmatr(k,12)=k*detat;%%%%



x=[20*pianyi 16*pianyi 12*pianyi]';%%%%
v(k,1)=A0727; v(k,2)=A0925; v(k,3)=A1123;  %%%%
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,13)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,13)=17*pianyi;         %%%%
tmatr(k,13)=k*detat;%%%%



x=[21*pianyi 17*pianyi 13*pianyi]';%%%%
v(k,1)=A0728; v(k,2)=A0926; v(k,3)=A1124;  %%%%
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,14)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,14)=17.5*pianyi;         %%%%
tmatr(k,14)=k*detat;%%%%

x=[22*pianyi 18*pianyi 14*pianyi]';%%%%
v(k,1)=A0729; v(k,2)=A0927; v(k,3)=A1125;  %%%%
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,15)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,15)=18*pianyi;         %%%%
tmatr(k,15)=k*detat;%%%%



x=[23*pianyi 19*pianyi 15*pianyi]';%%%%
v(k,1)=A0730; v(k,2)=A0928; v(k,3)=A1126;  %%%%
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,16)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,16)=18.5*pianyi;         %%%%
tmatr(k,16)=k*detat;%%%%



x=[20*pianyi 16*pianyi 12*pianyi]';%%%%
v(k,1)=A0929; v(k,2)=A1127; v(k,3)=A1325;  %%%%
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,17)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,17)=19*pianyi;         %%%%
tmatr(k,17)=k*detat;%%%%



x=[21*pianyi 17*pianyi 13*pianyi]';%%%%
v(k,1)=A0930; v(k,2)=A1128; v(k,3)=A1326;  %%%%
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,18)=v(k,1)+v(k,2)+v(k,3);%%%%
xmatr(k,18)=19.5*pianyi;         %%%%
tmatr(k,18)=k*detat;%%%%


x=[22*pianyi 18*pianyi 14*pianyi]';%%%%