gongpaodianfenxi.asv

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

%function [ output_args ] = gongpaodianfenxi( input_args )

detat=4.0e-7;speed=1800;pianyi=0.005;lth=1024;
k1=1:lth;
t=k1*detat;
%C25的强度，波速一般是3800M/S——4170M/S 
%C30的强度，波速一般是4200M/S以上
%load alldata.mat;
number=1024;
[sigmatr,sigmatr2]=dataget(number,detat,pianyi,speed);
 for k=1:number
   n=1;
   m=1;  A0121(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0122(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0123(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A0124(k,1)=sigmatr2(k,m);                                                        
   n=3;
   m=m+1;A0319(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0320(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A0321(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0322(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0323(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0324(k,1)=sigmatr2(k,m);                                                     
   m=m+1;A0325(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0326(k,1)=sigmatr2(k,m);                                                       
   n=5;
   m=m+1;A0517(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0518(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0519(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0520(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0521(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0522(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0523(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0524(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0525(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0526(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0527(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0528(k,1)=sigmatr2(k,m);                                                     
   n=7;
   m=m+1;A0719(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0720(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0721(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0722(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0723(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0724(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0725(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0726(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0727(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0728(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0729(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0730(k,1)=sigmatr2(k,m);                                                        
   n=9;
   m=m+1;A0921(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0922(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0923(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A0924(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0925(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0926(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0927(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0928(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0929(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0930(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A0931(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A0932(k,1)=sigmatr2(k,m);                                                      
   n=11; 
   m=m+1;A1123(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1124(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1125(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1126(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1127(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1128(k,1)=sigmatr2(k,m);                                                     
   m=m+1;A1129(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1130(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1131(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1132(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1133(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A1134(k,1)=sigmatr2(k,m);                                                       
   n=13;
   m=m+1;A1325(k,1)=sigmatr2(k,m);                                                     
   m=m+1;A1326(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1327(k,1)=sigmatr2(k,m);                                                     
   m=m+1;A1328(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1329(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1330(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1331(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1332(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1333(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1334(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1335(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1336(k,1)=sigmatr2(k,m);                                                     
   
   n=15; 
   m=m+1;A1527(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1528(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1529(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1530(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1531(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A1532(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1533(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1534(k,1)=sigmatr2(k,m);                                                       
   
   n=17;
   m=m+1;A1729(k,1)=sigmatr2(k,m);                                                       
   m=m+1;A1730(k,1)=sigmatr2(k,m);                                                        
   m=m+1;A1731(k,1)=sigmatr2(k,m);                                                      
   m=m+1;A1732(k,1)=sigmatr2(k,m); 
 end     
 figure;
   plot(A1729/4+pianyi*12,t);hold on;
   plot(A1730/4+pianyi*13,t);hold on;
   plot(A1731/4+pianyi*14,t);hold on;
   plot(A1732/4+pianyi*15,t);
   
   figure;
   plot(A1527/4+pianyi*12,t);hold on;
   plot(A1528/4+pianyi*13,t);hold on;
   plot(A1529/4+pianyi*14,t);hold on;
   plot(A1530/4+pianyi*15,t);hold on;
   plot(A1531/4+pianyi*15,t);hold on;
   plot(A1532/4+pianyi*15,t);hold on;
   plot(A1533/4+pianyi*15,t);hold on;