ch3example16prg5.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 11 行

% ch3example16prg5.m
fs=500;            % 采样率
Df=1;              % 频率分辨率
N=floor(fs/Df)+1;  % 计算的序列点数
t=0:1/fs:(N-1)/fs; % 截取信号的时间段
F=0:Df:fs;         % 功率谱估计的频率分辨率和范围
xk=sin(2*pi*50*t)+2*sin(2*pi*130*t)+randn(1,length(t));
                   % 截取时间段上的离散信号样点序列
[Pxx,F] = psd(xk,512,500,hamming(256),128); 
plot(F,10*log10(Pxx/(512/2))); xlabel('freq Hz');ylabel('PSD dB')
Pav=sum(Pxx/(512/2)) % 通过功率谱计算信号功率