ch3example16prg2.m

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

% ch3example16prg2.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(1:167))).^2+...
     abs(fft(xk(168:334))).^2+...
     abs(fft(xk(335:501))).^2)/3/((N/3)^2);
Pav_timedomaim=sum(xk.^2)/N  % 在时域计算信号功率
Pav_freqdomain=sum(Pxx)      % 通过功率谱计算信号功率
plot(0:3:fs,10*log10(Pxx));xlabel('freq Hz');ylabel('PSD dB');
%作出功率谱密度图