线性时不变离散时间系统的频域分析技术总结.docxVIP

云南大学信息学院数字信号处理实验报告 PAGE  PAGE 13 实验四 线性时不变离散时间系统的频域分析 实验室名称：格物楼 1204 实验时间：2015年11月6日 姓名： 成绩  教师签名： 年 月 日 一、实验目的 （1）学习MATLAB软件及其在信号处理中的应用，加深对常用连续时间信号的理解。 （2）连续时间信号在时域和频域中的抽样效果。 （3）熟悉MATLAB模拟滤波器的。 二、实验内容 Q4.1 修改程序P3.1中，取三个不同的M值，当时计算并画出式(2.13)所示滑动平均滤波器的幅度和相位谱。证明由幅度和相位谱表现出的对称类型。它表示了那种类型的滤波器？你现在能解释习题Q2.1的结果吗？修改后的程序P3.1也可以用于计算并画出线性时不变离散时间系统的频率响应，该系统用形如式(4.12)的传输函数描述。 Q4.2 使用修改后的程序P3.1，计算并画出当时传输函数 的因果线性时不变离散时间系统的频率响应。它表示那种类型的滤波器？ Q4.3 对下面的传输函数重做习题Q4.2 4.2和4.3给出的两个滤波器之间的区别是什么？你将选择哪一个滤波器来滤波，为什么？ Q4.6 使用zplane分别生成4.2和4.2所确定的两个滤波器的零极点图。讨论你的结果。 Q4.7 用程序P4.1计算并画出近似理想低通滤波器的冲激响应。低通有限冲激响应滤波器的长度是多少？在程序P4.1中，那个语句确定滤波器的长度？那个参数控制截止频率？ Q4.8 修改程序P4.1 ，计算并画出式(4.39)所示的长度为20，截止角频率为的有限冲激响应低通滤波器的冲激响应。 Q4.9 修改程序P4.1 ，计算并画出式(4.39)所示的长度为15，截止角频率为的有限冲激响应低通滤波器的冲激响应。 Q4.10 编写一个MATLAB程序，计算并画出式(4.39)所示有限冲激响应低通滤波器的振幅响应。使用这个程序，选取几个不同的N值，画出振幅响应并讨论你的结果。 Q4.11运行程序P4.2，计算并画出一个长度为2的滑动平均滤波器的增益响应。从图中验证3dB截止频率在π/2处。 Q4.23 用MATLAB产生如下两个因果系统传输函数的零极点图 研究生成的零极点图，你可以推断它们的稳定性吗？ 三、实验器材及软件 1. 微型计算机1台 MATLAB 7.0软件 四、实验原理 熟悉MATLAB中产生信号和绘制信号的基本命令，MATLAB已被开发成能对数据向量或矩阵运算的工具。序列以向量的形式储存，并且所有的信号被限定为因果的和有限长的，采用何种步骤执行程序，要根据MATLAB所运行的平台来决定。 若{h[n]}表示一个线性时不变离散时间系统的冲激响应，对{h[n]}做离散时间傅里叶变换得到其频率响应，即。 通常，是一个周期为2π的ω的复值函数，可以根据实部，虚部或者幅度相位来表示该函数。因此， 其中和分别是的实部和虚部，并且又|| 称为幅度响应，而称为线性时不变离散时间系统的相位响应。 线性时不变系统的增益函数g(ω)定义为 dB，增益函数的相反数a(ω)=-g(ω),称为衰减或损益函数。 对于用实冲激响应h[n]描述的离散时间系统，幅度函数是ω的偶函数，即|| =|| ；而相位函数是ω的奇函数，即。同样，是ω的偶函数，是ω的奇函数。 线性时不变离散时间系统的频率响应可以由输出序列y[n]的傅里叶变换与输入序列x[n]的傅里叶变换相比得到，即/。 线性时不变离散时间系统的冲激响应{h[n]}的z变换H(z)，称为传输函数或系统函数。H(z)可由输出序列y[n]的z变换Y(z)与输入序列x[n]的z变换X(z)相比得到，即H(z)=Y(z)/X(z)。 稳定因果系统函数的传输函数H(z)的所有极点都必须严格在单位圆内。 线性常系数差分方程描述的线性时不变系统，传输函数H(z)可以表示为 。 五、实验步骤 (1)打开MATLAB (2)新建M文件 (3)编写代码 (4)运行代码 (5)得到并分析结果 六、实验记录（数据、图表、波形、程序等） Q4.1 n = 0:100; s1 = cos(2*pi*0.05*n); % A low-frequency sinusoid s2 = cos(2*pi*0.47*n); % A high frequency sinusoid x = s1+s2; % Implementation of the moving average filter M = input(Desired length of the filter = ); num = ones(1,M); den = filter(num,1,x)