《数字信号处理》实验2 离散信号与系统的频域及复频域分析.pptxVIP

《数字信号处理实验与课程设计教程--面向工程教育》;实验目的：本实验包含3个实验子项目，通过matlab编程仿真，掌握离散信号与系统的频域及复频域分析方法，包括： 1. 利用傅里叶变换（FT）对离散信号与系统进行频域分析。 2. 利用Z变换对离散信号与系统进行复频域分析，掌握系统零极点绘制方式。 3. 熟悉Z变换的应用。;实验二 离散信号与系统的频域及复频域分析;2.1 离散时间信号与系统的频域分析---傅里叶变换(FT);3.本实验相关Matlab函数 (1)求信号x(n)的FT: X=fft(x,N) 其中：x—输入信号x(n);N---进行DFT的点数；X—x(n)的FT:X(ejω)。 (2)画信号x(n)的幅频谱和相频谱：plot(w,abs(fftshift(X))); plot(w,angle(fftshift(X))); 其中：w=2*pi*k/N---数字频率ω；(其中：k=-N/2:N/2-1) abs(fftshift(X))---x(n)的幅频谱|X(ejω)|(将频谱移动到以0为中心)。 angle(fftshift(X))---x(n)的相频谱φ(ω) (将频谱移动到以0为中心)。 二、实验环境 1. 计算机1台。 2. Windows 7 以上操作系统。 3. Matlab 7.0 以上版本软件。;三、实验参考和实验内容 1. 序列的傅里叶变换(FT) 矩形序列为：x(n)=R4(n) , 求：x(n)的傅里叶变换X(ejω)，画出其幅频谱|X(ejω)|和相频谱φ(ω)。 实验参考程序：参见实验项目参考源码\ch2\ch2prog1.m 实验内容： 1)运行该程序，画出信号 x(n)的幅频谱|X(ejω)|和相频谱φ(ω)，问： 它们各自有什么特点？ 2)编程实现：修改信号x(n)=δ(n), 画出信号 x(n)的幅频谱|X(ejω)|和相频谱φ(ω)。 3)编程实现：修改信号x(n)=anu(n) (a=0.5), 画出信号 x(n)的幅频谱|X(ejω)|和相频谱φ(ω)，问：a取什么范围时，X(ejω)才是存在的？;2. 傅里叶变换周期性、对称性和时移性 设序列为：x(n)= anu(n) (a=0.9) , 求：x(n)的傅里叶变换X(ejω)，画出其幅频谱|X(ejω)|和相频谱φ(ω)。为了观察傅里叶变换的周期性,设:做DFT点数为N=64，k=-2N~2N,即： X(ejω)≈ ω=2*pi*k/N（k=-2N~2N） 实验参考程序：参见实验项目参考源码\ch2\ch2prog2.m ch2prog2.m (主程序) 子函数:DFT2.m 实验内容： 1)运行该程序，画出信号 x(n)的幅频谱|X(ejω)|和相频谱φ(ω)。 2)观察幅频谱|X(ejω)|和相频谱φ(ω)，它们是否是周期信号？若是，则周期=? 3)幅频谱|X(ejω)|和相频谱φ(ω)是否具有对称性？哪个???偶对称？哪个是奇对称？ 4)编程实现：设序列x1(n)=x(n-n0) (n0=10), 画出信号 x(n-n0)的幅频谱和相频谱,问： 与x(n)的幅频谱和相频谱相比，有何不同？x(n-n0)的傅里叶变换理论公式是？实际波形是否验证了理论公式？;3.时域卷积定理的验证 实验内容： (1)编程实现y(n)=xa(n)*hb(n),其中：xa(n)= Ae-anTsin(Ω0nT)u(n)(0≤n50), A=1,a=0.4,Ω0=2.0734,T=1，hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)， 及Y(ejω)= FT[y(n)](取做DFT的点数N=128)，绘出|Y(ejω)|波形。 (2)编程实现Xa(ejω)=FT[Xa(n)](N=128)及Hb(ejωk)=FT[hb(n)](N=128)； 计算Y(ejω)= Xa(ejω) Hb(ejω)；绘出|Y(ejω)|波形。 问：(1)和(2)中的|Y(ejω)| 波形一致吗？为什么？ (3) 设该程序名为:ch2prog3.m,列出程序清单。（实验参考程序：参见实验项目参考源码\ch2\ch2prog3.m ) 四、实验报告要求 1. 简述实验目的。 2. 预习实验原理。 3. 实验结果及分析：包括注明程序注释、画出实验运行结果波形、回答实验中提出的问题、如果有程序设计要求的，列出所编写的程序清单并简要叙述程序调试过程。;2.2 离散时间信号与系统的复频域分析---Z变换;输出： “residuez”可将H(z)分解为简单部分分式之和： r,p,k的含义：p(1),p(2)…p(n)为列向量p;r(1),r(2)…r(n)为列向量r; k(1),k(