循环平稳信号处理完整的Matlab工具箱

该工具箱包括了循环平稳信号处理的各种常用函数,很全面，可以适用于通讯,雷达和声纳等多个领域。

function C3=cyclic_3rd_order_cumulant（x1x2x3alphamax_tau）
%cyclostationary_toolbox/cyclic_3rd_order_cumulant.m
%function 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



%function C3=cyclic_3rd_order_cumulant_fast（x1x2x3Tmax_tau）
%
% CYCLIC_3RD_ORDER_CUMULANT_FAST
%              
%              calculates the cyclic third order cumulant of 
%              three signals x1x2x3 at frequency alpha using a fast
%              approximation based on the synchronous average of the
%              time varying autocorrelation.  Fundamental signal
%              period can be defined as a single period or as a sequence
%              of once per period pulse times.
%             
%              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_fast（x1x2x3alphamax_tau）
%

% File: cyclic_3rd_order_cumulant_fast.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_third_order_cumulant_fast‘）;
end
if T（1）<1
  error（‘Synchronous period must be larger than 1 in function cyclic_third_order_cumulant_fast‘）;
end

[rowscols]=size（x1）;
if rows>cols
  x1=x1

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

     目录          0  2008-09-06 15:59  Cyclostationary toolbox

     文件      38689  2008-09-06 16:13  Cyclostationary toolbox\Contents.m

     文件       1564  2007-10-08 19:39  Cyclostationary toolbox\cyclic_3rd_order_cumulant.m

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

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

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

     文件         98  2007-10-11 20:27  Cyclostationary toolbox\cyclic_autocorrelation.m

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

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

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

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

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

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

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

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

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

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

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

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

     文件        349  2008-09-09 16:58  Cyclostationary toolbox\ex1.m

     文件        289  2007-10-13 08:45  Cyclostationary toolbox\ex2.m

     文件        337  2007-10-13 06:22  Cyclostationary toolbox\ex3.m

     文件        471  2007-10-17 18:07  Cyclostationary toolbox\ex4.m

     文件        452  2007-10-17 16:49  Cyclostationary toolbox\ex5.m

     文件        650  2008-09-09 16:24  Cyclostationary toolbox\ex6.m

     文件        314  2007-10-15 18:45  Cyclostationary toolbox\ex7.m

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

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

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

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

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