ex6_36.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 12 行

Fs=2000;
NFFT=1024;
n=0:1/Fs:1;
x=sin(2*pi*100*n)+4*sin(2*pi*500*n)+randn(size(n));  % 噪声序列
Cx=xcorr(x,’unbiased’);        % 计算序列的自相关函数
Cxk=fft(Cx,NFFT);            % 计算FFT
pxx=abs(Cxk);                % 求解PSD
t=0:round(NFFT/2-1);
k=t*Fs/NFFT;
P=10*log(pxx(t+1));           % 纵坐标的单位为dB
plot(k,P);