交叉小波，相干小波

交叉小波，相干小波代码及示例，通俗易懂，操作简单，

function [meanthetaRsigmaconfanglekappa]=anglemean（theta）;
% calculates the mean of angles
%
%    [meanthetaanglestrengthsigmaconfanglekappa]=anglemean（theta）; 
%    
% anglestrength: can be thought of as the inverse variance. [varies between 0 and one]
% sigma: circular standard deviation
% confangle: a 95% confidence angle （confidence of the mean value）
% kappa: an estimate of kappa in the Von Mises distribution
%
% check: http://www.cosy.sbg.ac.at/~reini/huber_dutra_freitas_igarss_01.pdf
%
% Aslak Grinsted 2002
theta=mod（theta（:）2*pi）;
n=length（theta）;
S=sum（sin（theta））;
C=sum（cos（theta））;
meantheta=atan2（SC）;

if nargout<2
    return
end
%if （（S>C）&（C>0））
%    meantheta=atan（S/C）;
%elseif （C<0）
%    meantheta=atan（S/C）+pi;
%else
%    meantheta=atan（S/C）+2*pi;
%end
%meantheta=mod（meantheta2*pi）;
Rsum=sqrt（S^2+C^2）;
R=Rsum/n;

if （R<.53）
    kappa=2*R+R^3+5*R^5/6;
elseif （R<.85）
    kappa=-0.4+1.39*R+0.43/（1-R）;
else
    kappa=1/（R^3-4

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-07-16 06:29  wtc-r16\
     文件        1612  2011-06-03 06:48  wtc-r16\anglemean.m
     文件        4003  2011-06-03 06:48  wtc-r16\ar1.m
     文件         780  2011-06-03 06:48  wtc-r16\ar1noise.m
     文件         408  2011-06-03 06:48  wtc-r16\ar1nv.m
     文件         676  2011-06-03 06:48  wtc-r16\ar1spectrum.m
     文件       54586  2011-06-03 06:48  wtc-r16\arrow.m
     文件        1270  2011-06-03 06:48  wtc-r16\boxpdf.m
     文件        1387  2011-06-03 06:48  wtc-r16\chisquare_inv.m
     文件         686  2011-06-03 06:48  wtc-r16\chisquare_solve.m
     文件       12041  2011-06-03 06:48  wtc-r16\colorbarf.m
     目录           0  2013-07-16 06:29  wtc-r16\example\
     文件       23165  2011-06-03 06:48  wtc-r16\example\faq.html
     文件        8812  2011-06-03 06:48  wtc-r16\example\faq.png
     文件        6454  2011-06-03 06:48  wtc-r16\example\faq_01.png
     文件       80940  2011-06-03 06:48  wtc-r16\example\faq_02.png
     文件       17324  2011-06-03 06:48  wtc-r16\example\faq_03.png
     文件       41485  2011-06-03 06:48  wtc-r16\example\faq_04.png
     文件       40353  2011-06-03 06:48  wtc-r16\example\faq_05.png
     文件       35024  2011-06-03 06:48  wtc-r16\example\faq_06.png
     文件        6988  2011-06-03 06:48  wtc-r16\example\wtcdemo.html
     文件        8533  2011-06-03 06:48  wtc-r16\example\wtcdemo.png
     文件       37488  2011-06-03 06:48  wtc-r16\example\wtcdemo_01.png
     文件       54028  2011-06-03 06:48  wtc-r16\example\wtcdemo_02.png
     文件       37729  2011-06-03 06:48  wtc-r16\example\wtcdemo_03.png
     文件        1081  2011-06-03 06:48  wtc-r16\formatts.m
     文件        3109  2011-06-03 06:48  wtc-r16\jao.txt
     文件        2748  2011-06-03 06:48  wtc-r16\jbaltic.txt
     文件        2025  2011-06-03 06:48  wtc-r16\loadbnm.m
     文件         644  2011-06-03 06:48  wtc-r16\normalizepdf.m
     文件        5889  2011-06-03 06:48  wtc-r16\parseArgs.m
............此处省略17个文件信息