计算高阶累积量

程序是计算高阶累积量的。最高到6阶，其他可以以此类推，十分好用

function [f1f2h]=cumn1（x）
%此函数功能与cumn类似，可求输入信号x的高阶累积量，并通过分别设定门限值t_f1和t_f2来识别信号的调制方式，并根据识别结果的不同，分别赋予h=2、4、8
%被rate_bpsk，rate_qpskrate_8psk程序调用
sig=x;
m20=mean（sig.^2）;
m21=mean（sig.*conj（sig））;
m40=mean（sig.^4）;
m41=mean（sig.^3.*conj（sig））;
m42=mean（（sig.^2）.*（conj（sig）.^2））;
m63=mean（（sig.^3）.*（conj（sig）.^3））;
c21=m21;
c40=m40-3*m20^2;
c41=m41-3*m21.*m20;
c42=m42-（abs（m20））^2-2*m21^2;
c63=m63-9*c42*c21-6*c21.^3;
f1=（abs（c40））/（abs（c42））;
f2=（abs（c41））/（abs（c42））;
t_f1=0.5;
t_f2=0.5;
if f1>t_f1&f2>t_f2
    h=2;
else if f1>t_f1&f2        h=4;
    else h=8;
    end
end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         678  2014-04-27 12:31  cumn1.m