先描述下我的问题，我想求我抽样的信号(横坐标时间，纵坐标电压)的噪声功率谱密度。方法一：由维纳辛钦定理可以知道，我需要先求时域信号的自相关函数，然后做DFT得到功率谱，在除以Fs可以得到功率谱密度

方法二：对我抽样的信号作DFT得到电压信号的频谱，然后幅值平方除以Nfft，再除以Fs可以得到功率谱密度

以上两种办法，很多大神都给出了详细的解释。

然后我的问题是我对同一个抽样信号求噪声功率谱密度结果不一样，求大神指点！

现贴下我的代码

先给出   NEW_Cal_dF35Hz_Match_4936.799MHz_LO16dBm_RF-24.6dBm.out   部分数据

-0.001049   0.000054

-0.001042   0.000061

-0.001043   0.000061

-0.001050   0.000060

-0.001062   0.000039

-0.001052   0.000055

-0.001045   0.000076

-0.001053   0.000057

-0.001061   0.000045

-0.001060   0.000050

-0.001057   0.000060

-0.001061   0.000061

-0.001070   0.000044

-0.001078   0.000037

-0.001066   0.000047

-0.001071   0.000050

-0.001062   0.000047

-0.001077   0.000032

-0.001065   0.000050

-0.001078   0.000033

-0.001070   0.000050

-0.001084   0.000032

-0.001082   0.000034

-0.001083   0.000037

-0.001085   0.000030

-0.001090   0.000026

-0.001086   0.000028

-0.001087   0.000036

-0.001086   0.000033

-0.001097   0.000029

-0.001097   0.000019

-0.001097   0.000028

-0.001096   0.000039

-0.001098   0.000027

-0.001108   0.000025

-0.001099   0.000026

-0.001107   0.000036

-0.001102   0.000026

-0.001110   0.000014

-0.001112   0.000009

-0.001119   0.000014

-0.001102   0.000038

-0.001121   0.000010

-0.001123   0.000006

-0.001125   0.000010

-0.001119   0.000026

-0.001128   0.000010

-0.001133   0.000004

-0.001140   -0.000000

-0.001128   0.000017

-0.001130   0.000006

-0.001136   -0.000012

-0.001144   -0.000004

-0.001148   -0.000004

-0.001135   0.000005

-0.001151   -0.000019

-0.001147   -0.000009

-0.001153   -0.000011

-0.001141   0.000021

-0.001147   0.000007

-0.001156   -0.000015

-0.001160   -0.000019

-0.001158   -0.000013

-0.001158   -0.000007

-0.001161   -0.000012

-0.001164   -0.000004

-0.001154   -0.000010

-0.001159   -0.000015

-0.001165   -0.000030

-0.001176   -0.000027

-0.001168   -0.000021

-0.001173   -0.000018

-0.001185   -0.000035

-0.001184   -0.000037

-0.001179   -0.000027

-0.001177   -0.000026

-0.001182   -0.000027

-0.001182   -0.000037

-0.001191   -0.000036

-0.001183   -0.000024

-0.001193   -0.000047

-0.001191   -0.000036

-0.001197   -0.000049

-0.001193   -0.000045

-0.001196   -0.000040

-0.001197   -0.000039

-0.001204   -0.000039

-0.001208   -0.000057

-0.001204   -0.000044

-0.001208   -0.000050

-0.001209   -0.000052

-0.001209   -0.000054

-0.001206   -0.000043

-0.001206   -0.000040

-0.001220   -0.000063

-0.001221   -0.000069

-0.001219   -0.000055

-0.001223   -0.000059

-0.001227   -0.000063

-0.001227   -0.000067

方法一：

clear;

clc;

data1 = load('NEW_Cal_dF35Hz_Match_4936.799MHz_LO16dBm_RF-24.6dBm.out');

amp = sqrt(data1(:,1).^2+data1(:,2).^2);

d_amp = amp - mean(amp);

a = xcorr(d_amp,'unbiased');

Fs=2e5;

N=2e4;

Cxk = fft(a,N);

Pxx_1=abs(Cxk)/Fs;

Pxx = Pxx_1;

t=0:round(N/2-1);

k=t*Fs/N;

P=10*log10(Pxx(t+1));

figure(2)

hold on;

plot(k,P);

set(gca,'XScale','log','Fontsize', 20);

untitled.png (25.29 KB, 下载次数: 6)

2018-9-19 10:43 上传

方法二：

data1 = load('NEW_Cal_dF35Hz_Match_4936.799MHz_LO16dBm_RF-24.6dBm.out');

amp = sqrt(data1(:,1).^2+data1(:,2).^2);

d_amp = amp -mean(amp);

Fs =2e5;

N = 20000;

X = fft(d_amp,N);

Px1 = X.*conj(X)/Fs/Fs/(N/Fs);

t=0:round(N/2-1);%PS=(FFT*FFT')/NFFT      PSD=(FFT*FFT')/(NFFT *Fs)

k=t*Fs/N;

figure(2)

plot(k,10*log10(Px1(t+1)*2));

2.png (41.26 KB, 下载次数: 6)

方法二

2018-9-19 10:46 上传

我想知道为什么方法一得到的噪声功率谱与方法二有这么大的区别？