北邮DSP Matlab 仿真实验
《数字信号处理》Matlab 实验
一. 离散信号的 FFT 分析
知识点:利用FFT 对信号频谱进行分析,用DFT 进行信号分析时基本参数的选择,以及信号经过离散时间傅立叶变换(DTFT)与有限长度离散傅立叶变换(DFT)后信号频谱上的区别。 实验教学内容:
1、用Matlab 编程上机练习。已知: N=25。这里Q=0、9+j0、3。可以推导出 ,
首先根据这个式子计算X(k)的理论值,然后计算输入序列x(n)的32个值,再利用基2时
间抽选的FFT 算法,计算x(n)的DFT X(k),与X(k)的理论值比较(要求计算结果最少6位有效数字)。
解:
format long
Q=0、9+0、3i;
WN=exp(-2*pi*1i/32);
Xk=(1-Q^32)、/(1-Q*WN 、^[0:24]);
xn=Q 、^[0:24];
Xkfft=fft(xn,32);
for (k0=1:1:25)
difference=Xk(k0)-Xkfft(k0);
end;
subplot(3,1,1);stem(abs(Xk(1:1:24)),'、');title('DFT
x(n)');xlabel('k');axis([0,35,0,15]);
subplot(3,1,2);stem(abs(Xkfft(1:1:32)),'g 、');title('FFT
x(n)');xlabel('k');axis([0,35,0,15]);
subplot(3,1,3);stem(abs(difference(1:1:25)),'r 、
');title('Xk-Xkfft');xlabel('k');axis([0,35,0,15]); 0n N-1()0 n 0, n N
n Q x n ⎧≤≤=⎨
0-=--===∑∑-=-=N k QW Q QW W n x k X N N n N N n N
N n Λ,