高阶累积量 matlab

clc ;
clear all;
close all;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x1=cell2mat（struct2cell（load（‘data.mat‘）））;
            %Abar = A（RecInd:）;
% 二阶距 M20 = E[X（k）X（k）] M21 = E[X（k）X‘（k）]
m20（:） = mean（x1（:）.*x1（:））;
m21（:） = mean（x1（:）.*conj（x1（:）））;

% 四阶距 M40 = E[X（k）X（k）X（k）X（k）]  M41 = E[X（k）X（k）X（k）X‘（k）]  M42 = E[X‘（k）X（k）X（k）X‘（k）]
m42（:） = mean（x1（:）.*x1（:）.*conj（x1（:））.*conj（x1（:）））;
m40（:） = mean（x1（:）.*x1（:）.*x1（:）.*x1（:））;
%=========计算高阶累计量===========      
c21（:）=m21（:）;
c42（:）= m42（:）-abs（m20（:））.^2-2*m21（:）.^2;
c40（:） = m40（:）-3*m20（:）.^2;
result1=c42/c40;%自定义的识别特征
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%C60%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 二阶距 M20 = E[X（k）X（k）] M21 = E[X（k）X‘（k）]
x1（:） = x1（:）-mean