power_1.m

来自「理解随机信号功率谱分析原理和方法。 (1) 生成信号」· M 代码 · 共 5 行

Fs=1000;
t=0:1/Fs:1;
xn=cos(20*pi*t)+cos(40*pi*t)+randn(size(t));
pxx=abs(fft(xn,1024)).^2/1001;
plot((0:1023)/1024*Fs,10*log10(pxx));