DSP实验2006.docVIP

数 字 信 号 处 理 实 验 指 导 书 电子信息工程专业 2003年9月 目 录 实验一 离散序列的傅氏变换 试验二 离散系统的频域分析 试验三 离散傅氏变换及其应用 试验四 快速傅氏变换 试验五 IIR数字滤波器的设计IR数字滤波器的设计 熟悉利用Matlab实现离散序列的傅氏变换，研究序列傅氏变换的幅频特性与相频特性。 二、试验内容及要求： 编制Matlab程序，实现以下内容 1研究序列x(n)=0.8nu(n)的离散时间傅氏变换。 （绘制图形，示例s101.m） 2研究序列x(n)=RN(n)的离散序列傅氏变换。 （绘制图形，示例s102.m） 三、 试验报告应包含以下内容 试验目的。 编制Matlab程序，实现以下内容 1 一个LTI系统，其差分方程为y(n)-0.5y(n-1)=x(n)+0.5x(n-1) ⑴求系统的脉冲响应及频率响应。 ⑵如果系统的输入为x(n)=cos，求系统响应y(n)。 （绘制图形，示例s201.m） 2令xa(t)=e-1000|t|，求出并绘制其傅氏变换Xa(jΩ)。用三种不同的采样频率对其进行采样，分别求出并画出其离散序列傅氏变换X(ejω)。⑴fs=5kHz ⑵fs=2kHz ⑶fs=1kHz （绘制图形，示例s202.m） 三、 试验报告应包含以下内容 试验目的。 编制Matlab程序，实现以下内容 1 x(n)=R5(n)，求N分别取10、20、50时的X(k)，与X(ejω)比较。 （绘制图形，示例s301.m） 2利用DFT计算两序列x1(n)=[1,2,2,3]，x2(n)=[1,2,3,4,3,2]的线性卷积。 （绘制图形，示例s302.m） 三、 试验报告应包含以下内容 试验目的。 1 已知一信号x(n)=R10(n)，用CZT算法求其前30点的复频谱X(zk)，给定A0=0.8，W0=1.2，?0=?/3，?0=?/10。 （绘制图形，示例s401.m） 2 利用FFT分析信号频率成分。一被噪声污染的信号很难看出它所包含的频率分量，如一个由50Hz和120Hz正弦信号构成的信号，受到均指随机噪声的干扰，数据采样频率为1000Hz，通过FFT来分析其信号频率成分。 （示例s402.m） 三、 试验报告应包含以下内容 试验目的。 试验五 IIR数字滤波器的设计 一、实验目的： 1、掌握脉冲响应不变法与双线性变换法设计IIR数字滤波器的设计理论。 2、了解脉冲响应不变法与双线性变换法设计IIR滤波器的应用。 实验内容：编制Matlab程序，实现以下内容 1、设计低通数字滤波器，通带截止频率Ωp = 0.2πrad，δ1＝1dB，阻带截止频率Ωs = 0.3πrad，δ2＝15dB，用脉冲响应不变法与双线性变换法设计，T=1。（模拟滤波器采用巴特沃斯滤波器原型，绘制图形） 2、一个含有高频噪声的信号保存在数据文件noisbloc中，试用脉冲响应不变法或双线性变换法设计的低通滤波器对其进行除噪，观察比较输入与输出信号。 三、实验原理： 设计原理：利用模拟滤波器来设计数字滤波器要先根据给出的性能指标设计出相应的模拟滤波器系统函数Ha(s)，然后由Ha(s)经变换得到所需的数字滤波器的系统函数H(z)。 1、脉冲响应不变法： 模拟滤波器系统函数，可化简为部分分式合形式，由冲激响应不变准则h(n)=Tha(nT)，经过推导得出所设计的数字滤波器系统函数，比较Ha(s)与H(z)，可看到二者各部分分式系数Ai相同，H(z)各极点eSiT 与Ha(s)各极点si对应，因此只要将模拟滤波器的系统函数Ha(s)分解成部分分式和形式，就能立即写出相应的数字滤波器的系统函数H(z)。 对于同一传输函数，存在着多种不同的算法结构，即对同一个系统有多种不同的实现方案。不同的方案，数学运算复杂度不同，所用加法器、乘法器、延迟器数目不同，系统特性对于乘法器系数变化的灵敏度不同。选择乘法器少的结构，提高速度选择延迟器少的结构，减少存储电路选择对乘法器系数变化灵敏度低的结构，提高系统稳定性，算法结构的选择对于一个系统的实现非常重要。设计出H(z)后用并联型系统结构实现，与级联型结构比较，并联型结构的误差更小。 2、双线性变换法：（为克服混叠失真这一缺点而提出） 首先把整个s平面压缩到某一中介s1平面的一条宽度为2π/T（-π/T～π/T ）的横带里，压缩关系为，利用欧拉公式展开，扩展到整个s平面，令s =jΩ， s1=jΩ1，，再由标准变换关系z=eS1T，将此横带变换到整个z平面上，变换关系为，，这样使s平面与z平面是一一对应关系，可消除混叠失真。设计出H(z)后用并联型系统结构实现。 、实验：试验目的。 试验IR数字滤波器的设计 一、实验目的： 1、掌握窗口