差分方程的变换解.pptVIP

差分方程的变换解*第1页，共16页，星期日，2025年，2月5日学习使用Matlab的符号运算Z变换和反Z变换方法。以及反Z变换中的部分分式展开法。加深对Z变换的理解。学习用Matlab计算差分方程的方法。加深对离散系统Z变换分析的理解，对零输入响应、零状态响应的理解。实验目的*第2页，共16页，星期日，2025年，2月5日实验原理与说明1、Z变换和反Z变换的符号运算：MATLAB的符号运算工具箱中，专门提供了Z变换和反Z变换的函数。正变换的调用格式为F=ztrans(f)式中，f为时间函数的符号表达式，F为Z变换式，也是符号表达式。反变换的调用格式为f=iztrans(F)式中，F为Z变换式的符号表达式，f为时间函数，是符号形式。为了改善公式的可读性，MATLAB提供了pretty函数，调用格式为Pretty(f)式中，f为符号表达式。*第3页，共16页，星期日，2025年，2月5日实验原理与说明2、求反Z变换的部分分式法若为有理式，则可表达为MATLAB提供了一个对residue()，其调用形式为[r,p,k]=residue(N,D)进行部分分式展开的函数式中，N和D分别为的分子多项式和分母多项式的系数向量，r为部分分式的系数向量，p为极点向量，k为多项式的系数向量。*第4页，共16页，星期日，2025年，2月5日实验原理与说明3、差分方程的Z变换解若线性常系数差分方程描述的系统为：（1）已知零输入初始值和对上式两边取变换有：上式的第一项为零输入响应，第二项为零状态响应。*第5页，共16页，星期日，2025年，2月5日实验原理与说明（2）已知系统初始值和对原方程式两边取变换有*第6页，共16页，星期日，2025年，2月5日试求下列序列的Z变换。计算示例1以1例说明(e)解：用Matlab计算的命令如下：F=ztrans(sym(k-3))%计算(e)F=z/(z-1)^2-3*z/(z-1)F=subs(F,z^-1)%根据反折性质，变量代换z换成1/zF=1/z/(1/z-1)^2-3/z/(1/z-1)*第7页，共16页，星期日，2025年，2月5日计算示例2已知，收敛域，求。解部分分式展开式为或*第8页，共16页，星期日，2025年，2月5日计算示例3描述某离散系统的差分方程为激励信号，若初始条件，试分别求其零输入响应、零状态响应和全响应。*第9页，共16页，星期日，2025年，2月5日计算示例3程序运行后在命令窗口显示的结果：零状态响应nn-(-1)+1/3(-2)+2/3零输入响应nn-5(-1)+2(-2)全响应nn-6(-1)+7/3(-2)+2/3*第10页，共16页，星期日，2025年，2月5日实验内容1求下列序列的变换，并注明收敛域。(a)(b)(c)(d)*第11页，共16页，星期日，2025年，2月5日实验内容2求下列的逆变换。(a)(b)(e)(d)(c)(f)*第12页，共16页，星期日，2025年，2月5日