plot_psd2.m

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

function plot_psd2(x, y, Fs, Nfft)

figure
[Pxx, f] = psd(x, 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;
[Pxx, f] = psd(y, Nfft, 1, Nfft, Nfft/2);
Pxx = Pxx./max(Pxx);
plot((f-1/2)*Fs, 10*log10(fftshift(Pxx)),'m');
hold on; grid on; zoom on;
xlabel('Frequency (MHz)');
ylabel('Watts/Hz (dB)');
title('Power Spectral Density');
xlim([-Fs/2,Fs/2]);