信号与系统作者张延华等第6章-z变换信号与系统书稿-6-10课案.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * 中国联合网络有限公司北京市分公司委托研究开发项目 中国联合网络有限公司北京市分公司委托研究开发项目 中国联合网络有限公司北京市分公司委托研究开发项目 中国联合网络有限公司北京市分公司委托研究开发项目 ThemeGallery PowerTemplate §6-10 应用示例及MATLAB实验 国家“十二五”规划教材——《信号与系统》 重点 难点 MATLAB仿真实验的设计 MATLAB仿真实验的设计 内容安排 6-10-1 正弦输入的响应 6-10-2 噪声序列的滤波 6-10-3 MACD决策 6-10-1 正弦输入的响应 如前所述，线性时不变系统的系统函数 一般可以用 一个关于z的有理函数形式来描述： (6-10-1) 设系统是因果、稳定的，即 的全部极点均位于单位圆内， 则系统的频率响应特性可以通过正弦输入的响应来确定。 令系统输入序列 为： 式中C和 均为实数。 的z变换为： 如果初始条件为零，则系统输出 的z变换为： (6-10-2) 该响应的极点是 的极点加上 的根。 解出一对共轭复数根 。现对 进行部分分式展开： 式中 是z的多项式，其阶次低于N。系数c是待定系数， 是c的共轭复数，根据系数公式有： 6-10-1 正弦输入的响应 (6-10-3) (6-10-4) 因此，系统输出 的z变换 为： (6-10-5) 6-10-1 正弦输入的响应 对上式取逆z变换： (6-10-6) 其中等式右端第一项 是 的逆z变换，因为系统已知是稳定的， 的根全部位于单位圆内，故当 时， ，所以 就是在正弦输入序列 作用下的暂态响应。等式右端第二项 是由一对共轭复数极点 形成的一对共轭复数项，它的逆z变换 可写成通式： (6-10-7) 6-10-1 正弦输入的响应 式(6-10-7)中的系数和极点可分别写成它的极坐标的形式： 和 将上述关系式代入式(7-10-7)，并利用欧拉公式得到： 6-10-1 正弦输入的响应 就是在正弦序列 作用下，系统的稳态响应。 (6-10-8) 显然，当 时，响应 并不趋于零，所以 因此有： 6-10-1 正弦输入的响应 6-10-2 噪声序列的滤波 指数加权移动平均(EWMA)滤波器的单位样值响应 已知为: 因为 为有限值，所以EWMA是一个FIR滤波器。 令 ，则可以将L点EWMA滤波器扩展成一个IIR滤波器。 因此，IIR EWMA滤波器的单位样值响应 就为： (6-10-9) 对上式求z变换，就可得到IIR EWMA滤波器的系统函数 (6-10-10) 式中 当 时，由式(7-10-10)可得到IIR EWMA滤波器的 递归方程如下： (6-10-12) 式(6-10-12)在金融领域又被称之为EMA(Exponential Moving Average)滤波器或指数移动平均线，一般用下面的标准形式给出： 6-10-2 噪声序列的滤波 (6-10-13) 对上式重新整理，得其等价形式为： (6-10-14) 如果 ，则 ，这将意味着系统没有滤波作用。 如果 从0开始增大，且 从0到增加 ，则系统的频率响应曲线有一个陡降。 6-10-2 噪声序列的滤波 下面程序给出了IIR EWMA或EMA滤波器的频率响应： a=0.2; W=0:.01:1; omega=W*pi; z=exp(j*omega); H=a*z./(z-(1-a)); magH = abs(H); angH = 180/pi*unwrap(angle(H)); subplot(211); plot(W,magH); grid on; axis([0,1,0,1]) xlabel(frequency in pi units); ylabel(|H|); subplot(212); plot(W,angH); grid on xlabel(frequency in pi units); ylabel(Phase in pi Radian