matlab中的循环平稳工具箱

其中包括：循环平稳的自相关函数，循环谱分析等内容。

cyclostationary_toolbox/cyclic_3rd_order_cumulant.m010064400316540031654000000030340643656516400233250ustar00acmcacmc00003040001726function C3=cyclic_3rd_order_cumulant（x1x2x3alphamax_tau）
%
% CYCLIC_3RD_ORDER_CUMULANT
%              
%              calculates the cyclic third order cumulant of 
%              three signals x1x2x3 at frequency alpha
%             
%              C3（k*alphatau1tau2）=E{（x1（t）-E{x1（t）}） *
%                                      （x2（t+tau1）-E{x2（t+tau1）} *
%             			       （x3（t+tau2）-E{x3（t+tau2）} *
%                                      exp（-jk（alpha）t）  }
%              for k=0 ... 1/alpha
%             
%
% USAGE
%              C3=cyclic_3rd_order_cumulant（xyalphamax_tau）
%

% File: cyclic_3rd_order_cumulant.m
% Last Revised: 25/11/97
% Created: 25/11/97
% Author: Andrew C. McCormick
% （C） University of Strathclyde

% Simple error checks
if nargin~=5
  error（‘Incorrect number of arguments for function cyclic_3rd_order_cumulant‘）;
end
if alpha>2*pi
  error（‘Cyclic frequency must be less than 2 pi in function cyclic_3rd_order_cumulant‘）;
end


% Remove cyclic mean from signals
cmx1=cyclic_mean（x1alpha）;
cmx2=cyclic_mean（x2alpha）;
cmx3=cyclic_mean（x3alpha）;
lx=length（x1）;
t=0:lx-1;
T=ceil（2*pi/alpha）-1;
for k=1:lx
  x1（k）=x1（k）-1/（2*pi）*sum（cmx1.*exp（j*alpha*（0:T）*（k-1）））;
  x2（k）=x2（k）-1/（2*pi）*sum（cmx2.*exp（j*alpha*（0:T）*（k-1）））;
  x3（k）=x3（k）-1/（2*pi）*sum（cmx3.*exp（j*alpha*（0:T）*（k-1）））;
end

C3=zeros（max_taumax_tauT+1）;

ix=1:lx-max_tau-1;

for tau1=0:max_tau
  for tau2=0:max_tau
    for k=0:T
      C3（tau1+1tau2+1k+1）=mean（x1（ix）.*x2（tau1+ix） ...
	  .*x3（tau2+ix）.*exp（j*k*alpha*t（ix）））;
    end
  end
end



cyclostationary_toolbox/cyclic_3rd_order_cumulant_fast.m010064400316540031654000000060740643655537000243470ustar00acmcacmc00003040001726

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

     文件      58880  2000-01-31 19:11  Contents.m

     文件       1564  1997-11-25 23:10  cyclic_3rd_order_cumulant.m

     文件       3132  1997-11-25 22:04  cyclic_3rd_order_cumulant_fast.m

     文件       2171  1997-11-25 23:11  cyclic_4th_order_cumulant.m

     文件       3983  1997-11-25 22:50  cyclic_4th_order_cumulant_fast.m

     文件        638  1997-11-24 20:41  cyclic_autocorrelation.m

     文件       1113  1997-11-24 20:48  cyclic_autocorrelation_fast.m

     文件        686  1997-11-24 20:44  cyclic_autocovariance.m

     文件       1215  1998-03-26 19:45  cyclic_autocovariance_fast.m

     文件       1047  1997-11-25 19:47  cyclic_correlation_coefficient.m

     文件       1381  1998-04-23 18:48  cyclic_cross_correlation.m

     文件       2767  1998-06-08 20:41  cyclic_cross_correlation_fast.m

     文件       1668  1998-04-23 18:51  cyclic_cross_covariance.m

     文件       3090  1998-06-24 20:09  cyclic_cross_covariance_fast.m

     文件       1435  1998-04-01 19:11  cyclic_cross_periodogram.m

     文件       1597  1997-11-25 18:51  cyclic_cumulants_fast.m

     文件        683  1997-11-25 23:12  cyclic_mean.m

     文件        671  1997-11-25 19:59  degree_of_cyclostationarity.m

     文件       1895  2000-01-31 19:09  fft_accumulation_method.m

     文件       2094  1998-05-05 19:12  frac_delay_cyclic_ac.m

     文件        972  1997-11-25 18:59  get_impulses.m

     文件        921  1998-04-23 18:19  scd.m

     文件       1368  2000-01-31 19:10  strip_spectral_correlation.m

     文件       1397  1998-06-24 20:03  synchronous_average.m

     文件        624  1997-11-25 18:14  wvd.m

----------- ---------  ---------- -----  ----

                96992                    25