已知信号x(n)=sin(2pf1n/Fs)+sin(2pf2n/Fs) (0£n£239),其中f1=40kHz,f2=40.5kHz,Fs=160kHz,试回答:
(a)要想分辨出两个频率点,FFT的点数至少应该为多少,并利用Matlab观察该信号的幅频特性。
(b)基于问题(a)的结果,根据幅频特性确定两个频率f1和f2的估计误差。
分享一下我的解法:
%分析:要想分辨两个频率点,则频谱的最小分辨率小于两频率差,最小分辨率
%f=fs/N<=|f1-f2|=0.5khz,得到N=320,同时N要满足2的整数倍才能使用fft,所以N=512;
clear;
n=0:239;
f1=40e+3;f2=40.5e+3;fs=160e+3;
x=sin(2pif1n/fs)+sin(2pif2n/fs);
N=512;
N1=0:N-1;
f=N1*fs/N;%横坐标换成对应的频率
y=fft(x,N);
plot(f,abs(y));
xlabel(‘khz’);ylabel(‘幅度’);title(‘N=512的幅频特性’);
[a1,b1]=max(abs(y(1:N/2)));%最大值索引
c1=(b1-1)*fs/N;%一个频率点
[a2,b2]=max(abs(y(1:b1-1)));
c2=(b2-1)*fs/N;%另一个频率点
f0=abs(c1-c2)%测量频率差值
fq=abs(f1-f2)%真实频率差值
ff=abs(f0-fq)/fq%估计误差
%结果分析:
估计误差为ff=0.25,其中f1估计值准确,f2估计值存在误差,可能原因:
%1.存在栅栏效应,漏掉了f2,导致f2测量不准确 2,存在泄露,由于N>M,泄露使得频谱分辨率减小,导致f2估计误差