|
发表于 2016-12-28 08:45:28
|
显示全部楼层
本帖最后由 chafenfa 于 2016-12-28 08:47 编辑
回复 20# langgege
你可以看一下下面的程序,仅供参考
%利用FFT计算SDM的PSD
fsig=1000;
M=256;
fs=2*M*fsig;
N=16384;
Ts=1/fs;
Ain=1;
ts=Ts;
Doutw=y.*blackman(N); %加窗
Dout_spect=fft(Doutw);
Dout_dB=20*log10(abs(Dout_spect));
maxdB=max(Dout_dB(1:N/2));
semilogx([0:N/2-1].*fs/N,Dout_dB(1:N/2)-maxdB);
title('Power Spectral Density');
xlabel('Frequency(Hz)');
ylabel('PSD(dB)');
接下来计算SDM的SNR,程序如下,
spectP=(abs(Dout_spect)).*(abs(Dout_spect));
Dout_dB=10*log10(spectP);
[maxdB,fin]=max(Dout_dB(1:N/2));
nb=5;
signal_bins=fin+[-(nb-1)/2:(nb-1)/2];
inband_bins=0:N/(2*M);
noise_bins=setdiff(inband_bins,signal_bins);
Ps=sum(spectP(signal_bins));
noise_bins(:,1)=[];
Pn=sum(spectP(noise_bins));
SNR=10*log10(Ps/Pn) |
|