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（‘信號的自功率譜‘）