power_2.m

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

Fs=1000;
t=0:1/Fs:1;
xn=cos(20*pi*t)+cos(40*pi*t)+randn(size(t));
pxx1=(abs(fft(xn(1:256))).^2+abs(fft(xn(257:512))).^2+abs(fft(xn(513:768))).^2)/(256*3);
plot((0:255)/256*Fs,10*log10(pxx1));