plot_psd.m

来自「基于matlab的Wimax数字中频DUC的仿真」· M 代码 · 共 15 行

function plot_psd(xin, Fs, Nfft)

if (nargin < 3)
  Nfft = 1024;
end

figure
[Pxx, f] = psd(xin, Nfft, 1, Nfft, Nfft/2);
Pxx = Pxx./max(Pxx);
plot((f-1/2)*Fs, 10*log10(fftshift(Pxx)),'b');
hold on; grid on; zoom on;
xlabel('Frequency (MHz)');
ylabel('Watts/Hz (dB)');
title('Power Spectral Density');
xlim([-Fs/2,Fs/2]);