频域特征值MATLAB程序

function [MNFMDFAP] = FrenquencyFeatures（data）

%能量
FFT=abs（fft（data））;
AP=sum（FFT.^2/256）;

%平均功率
fs=1000;
[LW]=size（data）;
for i=1:W
a=xcorr（data（:i）‘unbiased‘）;
AFFT=abs（fft（a））;
PSD=10*log10（AFFT）;
df=fs/L;
p=（sum（AFFT（1:L/2-1））+sum（AFFT（1:L/2）））/2.*df;
pf=（sum（AFFT（1:L/2-1）.*（1:L/2-1）‘.*df）+sum（AFFT（1