samp9_8.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 14 行

%Samp9_8
clf;Fs=1000;   %采样频率
N=1024;Nfft=256;n=0:N-1;t=n/Fs; %数据长度、时间序列
window=hanning(256);   %选用的窗口
noverlap=128;   %分段序列重叠的采样点数（长度）
dflag='none';  %不做趋势处理
randn('state',0);   %设置产生随机数的初始状态
xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);   %带噪声原始信号
Pxx=psd(xn,Nfft,Fs,window,noverlap,dflag);   %功率谱估计
f=(0:Nfft/2)*Fs/Nfft;   %求得对应的频率向量
plot(f,10*log10(Pxx));  %绘制功率谱
xlabel('频率/Hz');ylabel('功率谱/dB');
title('PSD—Welch方法'); grid on