資源簡介
計算信號的高階累積量，帶注釋，包括高階矩，親測可用，代碼是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評論
共有 條評論