EEPW论坛

第一个图画出来100hz幅值约在20dB左右，而第二为0，为什么？

Fs=4096;
nfft=1024;         %2^n
N=1024;
n=0:1/Fs:1;
x=sin(2*pi*100*n);
X=fft(x,nfft);
mag1=abs(X)*2/N;
F1=(0:N-1)*Fs/N;   %真实频率
figure,plot(F1(1:N/2),mag1(1:N/2)); %绘出Nyquist频率之前的振幅
Pxx=abs(X).^2/length(n);   %求解PSD
t=0:round(nfft/2-1);
f=t*Fs/nfft;
P=10*log10(Pxx(t+1));   %纵坐标的单位为dB
figure;plot(f,P);

clear;
Fs=4096; %采样频率
n=0:1/Fs:1;
% 产生含有噪声的序列
xn=sin(2*pi*100*n);
window=boxcar(length(xn)); %矩形窗
nfft=1024;
[Pxx,f]=periodogram(xn,window,nfft,Fs); %直接法
plot(f,10*log10(Pxx));