数字信号处理留的小作业

已知信号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估计误差