自相关功率谱的线谱提取

N=1000;
n=0:N-1;
Fs=500;
t=n/Fs;
y=4*randn（1N）;
% y=sin（100*t）;
% y2=3*randn（11000）;
 b=[1];  
 a=[0.0004 0.6 1]; 
 H=freqs（bay）;
subplot 321
plot（tH）
title（‘纯噪声信号‘）
xlabel（‘t/s‘）;
[clags]=xcorr（yN‘unbiased‘）;
subplot 323
plot（lagsc）
axis（[10 500 -3 3]）
title（‘信号的自相关函数‘）
xlabel（‘τ/ms‘）;
a=fft（cN）;
f=n*Fs/N; %频率序列
subplot 325
plot（fa）
title（‘信号的自功率谱‘）