数字信号处理实验报告材料43484.pdf

文档 数字信号处理 实验报告 实验一 序列的傅立叶变换 一、实验目的 1.进一步加深理解 DFS,DFT 算法的原理 ; 2.研究补零问题 ; 3.快速傅立叶变换（ FFT）的应用。 二、 实验步骤 1.复习 DFS 和 DFT 的定义，性质和应用； 2 熟悉 MATLAB 语言的命令窗口、编程窗口和图形窗口的使用； 3 利用提供的 程序例子编写实验用程序； 4.按实验容上机实验，并进行实验结果分析； 5.写出 完整的实验报告，并将程序附在后面 。 三、 实验容 1.周期方波序列的频谱 试画出下面四种情况下的的幅度频谱 , 并分析补零后，对信号频谱的影响。 x (n) cos(0.48 n) cos(0.52 n) 2.有限长序列 x(n) 的 DFT （1）取 x(n)(n=0:10) 时，画出 x(n) 的频谱 X(k) 的幅度； （2）将（ 1）中的 x(n) 以补零的方式，使 x(n)加长到 (n:0~100)时，画出 x(n) 的频 谱 X(k) 的幅度； （3）取 x(n)(n:0~100) 时，画出 x(n) 的频谱 X(k) 的幅度。利用 FFT 进行谱分析 x(t ) 2 sin(4 t ) 5 cos(8 t ) 3.已知：模拟信号 以 t=0.01n(n=0:N-1)进行采样，求 N 点 DFT 的幅值谱。 请分别画出 N=45; N=50;N=55;N=60 时的幅值曲线。 四、 实验数据分析 文档 1.周期方波序列的频谱分析 首先定义一个功能函数 dfs function[Xk]=dfs(xn,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^nk; Xk=xn*WNnk; (1)L=5 ，N=20; %题 1.(1) L=5;N=20;% 对于（ 2 ），（3），（4 ）问，只要修改 L,N 的数值就好。 n=1:N; xn=[ones(1,L),zeros(1,N-L)]; Xk=dfs(xn,N); magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]); k=[-N/2:N/2]; figure(1) subplot(2,1,1); stem(n,xn);xlabel('n');ylabel('xtide(n)'); title('DFS of SQ.wave:L=5,N=20'); subplot(2,1,2); stem(k,magXk); axis([-N/2,N/2,0,16]); xlabel('k');ylabel('Xtide(k)'); 文档 (2)L=5 ，N=40; (3).L=5 ，N=60 文档 （4 ）L=7 ，N=60; 结果分析 ：虽然周期序列不存在 FT,但是一个周期序列可以利用其 DFS系数X(k) 表示它的频谱分布规律， 从以上各频谱图可以看出， 随着补零点数的增加， 周期 文档 序列的谐波次数越来越多， 其频谱的包络线越来越平滑连续， 更能反