实验二_DFT_和FFT实验报告.docxVIP

实验二 DFT 和FFT 谢稳文 1、实验目的 认真复习周期序列 DFS、有限长序列DFT 的概念、旋转因子的定义、以及DFS 和DFT 的性质等有关内容；复习基2-FFT 的基本算法，混合基-FFT 的基本算法、Chirp-Z 变换的算 法等快速傅立叶变换的方法。 掌握有限长序列的循环移位、循环卷积的方法，对序列共轭对称性的含义和相关内容加 深理解和掌握，掌握利用DFT 分析序列的频谱特性的基本方法。 掌握 FFT 算法的基本原理和方法、Chirp-Z 变换的基本原理和方法，掌握利用FFT 分析 序列的频谱特性的方法。 熟悉利用 MATLAB 进行序列的DFT、FFT 的分析方法。 2、实验内容 a. 设周期序列( ) { ( n) = …,0,1,2,3,0,1,2,3,0,1,2,3,….}，求该序列的离散傅立叶级数 X?? (k) = DFS[(n)]，并画出DFS 的幅度特性。(离散傅里叶变换和快速傅里叶变换). b. 设周期方波序列为 x(n)= (m=0,) 其中 N 为基波周期，L/N 是占空比。 (1) 用L 和N求| X?? (k) |的表达式； (2) 当L 和N 分别为：L=5，N=20；L=5，N=40；L=5，N=60 以及L=7，N=60 时画出 DFS 的幅度谱； (3) 对以上结果进行讨论，总结其特点和规律 L=5; N=20; k=[-N/2:N/2]; xn=[ones(1,L),zeros(1,N-L)]; Xk=dfs(xn,N); y=abs([Xk(N/2+1:N) Xk(1:N/2+1)]); stem(k,y); c. 设有限长序列x(n) = {0,1,2,3}，计算DTFT[x(n)]=X(ejω)，并画出它的幅度谱；然后 利用kw1=k,k=0,1,2,3对X(ejω)进行采样，并证明它等于实验a 中的(k)。 X(ejω)==e-jw+2e-j2w+3e-j3w X(ej0)=1+2+3=6; X(ej2/4)=-2+2j; X(ej4/4)=2; X(ej6/4)=-2-2j; d. 序列x(n)=R4(n)，计算DTFT[x(n)]=X(ejω)，并绘制其幅度和相位谱。 (1) 计算x(n)的4 点DFT，并绘制DFT 的幅度与相位谱； (2) 将x(n)补零形成8 点序列，计算8 点DFT，并绘制幅度与相位谱，求频率分辨率； (3) 将x(n)补零形成16 点序列，计算16 点DFT，并绘制幅度与相位谱，求频率分辨率； （1） （2）实验代码: （3） e. 序列x(n)=cos(0.48πn)+cos(0.52πn) (1) 求x(n)的10 点DFT，并画出它幅度与相位谱； (2) 求x(n)的100 点DFT，并画出它幅度与相位谱； 根据实验结果，讨论 DFT 进行谱分析的条件。 f. 序列x(n)=5(0.9)nR11(n) (1) 求循环反转序列x((-n))11，并绘制x(n)和x((-n))11 的波形；求两序列的DFT，验证 DFT 的循环反转性质。 (2) 把序列x(n)分解成圆周共轭奇分量xoc(n)和圆周共轭偶分量xec(n)，并求出对应的 DFT，验证DFT 的圆周共轭对称性质。 (3) 绘制x((n+4))11R11(n)、x((n-3))15R15(n)和x((n-6))15 的波形，验证序列圆周移位性质。 n=0:10; x=5*(0.9).^n; y=x(mod(-n,11)+1); X=dft(x,11); Y=dft(y,11); subplot(2,2,1); stem(n,real(X)); subplot(2,2,2); stem(n,imag(X)); subplot(2,2,3); stem(n,real(Y)); subplot(2,2,4); stem(n,imag(Y)); 2） g. 序列x1(n)={1,2,1}，x2(n)={1,2,3,2}， (1) 编制程序在时域中直接计算4 点循环卷积1 4 2 y(n) = x (n)?x (n)，并绘制波形图； (2) 编制程序先计算x1(n)和x2(n)的4 点DFT，再利用IDFT 计算1 2 4 x (n)?x (n)，比较 结果以上两种计算结果是否一致？ (3) 计算1 2 5 x (n)?x (n)和1 2 6 x (n)?x (n)，分析循环卷积点数N 对循环卷积的影响，并 比较哪种循环卷积结果与线性卷积是一致的？总结出用循环卷积计算线性卷积的条件。 function y=circonvl(x1,x2,N) if length(x1)N error(N must be = the length of x1) end if len