提取地震子波的完整代码（Wave Extract for Earthquake HOSP matlab）

提取地震子波的完整代码，不仅包含了子波提取代码，还包含了高阶统计量工具箱的完整代码。

function window_4_order_temp=window_4_order（n）
%%% 窗函数 
for tao1=-n:n
    
 if（ abs（tao1）<=n/2 ） d1=1-6*（abs（tao1）/n）^2+6*（abs（tao1）/n）^3;
 elseif （ n/2<=abs（tao1）<=n ）  d1=2*（1-abs（tao1）/n）^3;
 else d1=0.0;
 end
  d_draw（tao1+n+1）=d1;
    for tao2=-n:n
     if（ abs（tao2）<=n/2 ） d2=1-6*（abs（tao2）/n）^2+6*（abs（tao2）/n）^3;
     elseif （ n/2<=abs（tao2）<=n ）  d2=2*（1-abs（tao2）/n）^3;
     else d2=0.0;
     end
        for tao3=-n:n
         if（ abs（tao3）<=n/2 ） d3=1-6*（abs（tao3）/n）^2+6*（abs（tao3）/n）^3;
         elseif （ n/2<=abs（tao3）<=n ）  d3=2*（1-abs（tao3）/n）^3;
         else d3=0.0;
         end
         
         tao12=tao2-tao1;
         if（ abs（tao12）<=n/2 ） d12=1-6*（abs（tao12）/n）^2+6*（abs（tao12）/n）^3;
         elseif （ n/2<=abs（tao12）<=n ）  d12=2*（1-abs（tao12）/n）^3;
         else d12=0.0;
         end
          tao23=tao3-tao2;
         if（ abs（tao23）<=n/2 ） d23=1-6*（abs（tao23）/n）^2+6*（abs（tao23）/n）^3;
         elseif （ n/2<=abs（tao23）<=n ）  d23=2*（1-abs（tao23）/n）^3;
         else d23=0.0;
         end
          tao13=tao3-tao1;
         if（ abs（tao13）<=n/2 ） d13=1-6*（abs（tao13）/n）^2+6*（abs（tao13）/n）^3;
         elseif （ n/2<=abs（tao13）<=n ）  d13=2*（1-abs（tao13）/n）^3;
         else d13=0.0;
         end
            window_4_order_temp（tao1+n+1tao2+n+1tao3+n+1）=（（d1*d2*d3*d12*d23*d13））;
        end
    end
end
tao1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4452  2011-06-09 21:49  高阶统计量\hosa\ARMAQS.M

     文件       2862  2011-06-09 21:49  高阶统计量\hosa\ARMARTS.M

     文件       4044  2011-06-09 21:49  高阶统计量\hosa\ARMASYN.M

     文件       4165  2011-06-09 21:49  高阶统计量\hosa\ARORDER.M

     文件       4034  2011-06-09 21:49  高阶统计量\hosa\ARRCEST.M

     文件       5449  2011-06-09 21:49  高阶统计量\hosa\BICEPS.M

     文件       5139  2011-06-09 21:49  高阶统计量\hosa\BICEPSF.M

     文件       4206  2011-06-09 21:49  高阶统计量\hosa\BICOHER.M

     文件       4717  2011-06-09 21:49  高阶统计量\hosa\BICOHERX.M

     文件       5854  2011-06-09 21:49  高阶统计量\hosa\BISPECD.M

     文件       6508  2011-06-09 21:49  高阶统计量\hosa\BISPECDX.M

     文件       5243  2011-06-09 21:49  高阶统计量\hosa\BISPECI.M

     文件       1974  2011-06-09 21:49  高阶统计量\hosa\BISPECT.M

     文件       5867  2011-06-09 21:49  高阶统计量\hosa\bispest.asv

     文件       4939  2011-06-09 21:49  高阶统计量\hosa\CONTENTS.M

     文件       2069  2011-06-09 21:49  高阶统计量\hosa\CUM2EST.M

     文件       3133  2011-06-09 21:49  高阶统计量\hosa\CUM2X.M

     文件       2902  2011-06-09 21:49  高阶统计量\hosa\CUM3EST.M

     文件       3611  2011-06-09 21:49  高阶统计量\hosa\CUM3X.M

     文件       4463  2011-06-09 21:49  高阶统计量\hosa\CUM4EST.M

     文件       5704  2011-06-09 21:49  高阶统计量\hosa\CUM4X.M

     文件       2764  2011-06-09 21:49  高阶统计量\hosa\CUMEST.M

     文件       4516  2011-06-09 21:49  高阶统计量\hosa\CUMTRUE.M

     文件       5955  2011-06-09 21:49  高阶统计量\hosa\DOA.M

     文件       5605  2011-06-09 21:49  高阶统计量\hosa\DOAGEN.M

     文件       8677  2011-06-09 21:49  高阶统计量\hosa\GLSTAT.M

     文件       6512  2011-06-09 21:49  高阶统计量\hosa\HARMEST.M

     文件       3757  2011-06-09 21:49  高阶统计量\hosa\HARMGEN.M

     文件       4260  2011-06-09 21:50  高阶统计量\hosa\HOSAHELP.M

     文件        706  2011-06-09 21:50  高阶统计量\hosa\HOSAVER.M

............此处省略148个文件信息