%========================开始提取加噪信号的各类特征值================================
for n=1:1:50;
m=n*Ns;
x=(n-1)*Ns;
for i=x+1:m; %提取加噪信号'signal_with_noise=y+noise'的前256个元素，抽取50次
y0(i)=signal_with_noise(i);
end
Y=fft(y0); %对调制信号进行快速傅里叶算法（离散）
y1=hilbert(y0) ; %调制信号实部的解析式
factor=0; %开始求零中心归一化瞬时幅度谱密度的最大值gamma_max
for i=x+1:m;
factor=factor+y0(i);
end
ms=factor/(m-x);
an_i=y0./ms;
acn_i=an_i-1;
end
gamma_max=max(fft(acn_i.*acn_i))/Ns
