ch3example16prg1.m

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

% ch3example16prg1.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=abs(fft(xk)).^2/(N^2);   % 功率谱估计
Pav_timedomaim=sum(xk.^2)/N  % 在时域计算信号功率
Pav_freqdomain=sum(Pxx)      % 通过功率谱计算信号功率
plot(F,10*log10(Pxx));xlabel('freq Hz');ylabel('PSD dB') % 作出功率谱密度图