资源描述：
《差分方程的解法分析及matlab实现(程序)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、差分方程的解法分析及MATLAB实现（程序）摘自：张登奇,彭仕玉.差分方程的解法分析及其MATLAB实现[J].湖南理工学院学报.2014(03)引言线性常系数差分方程是描述线性时不变离散时间系统的数学模型,求解差分方程是分析离散时间系统的重要内容.在《信号与系统》课程中介绍的求解方法主要有迭代法、时域经典法、双零法和变换域法[1].1迭代法例1已知离散系统的差分方程为,激励信号为,初始状态为.求系统响应.根据激励信号和初始状态,手工依次迭代可算出.利用MATLAB中的filter函数实现迭代过程的m程序如下:clc;clear;formatcom
2、pact;a=[1,-3/4,1/8],b=[1,1/3,0],%输入差分方程系数向量,不足补0对齐n=0:10;xn=(3/4).^n,%输入激励信号zx=[0,0],zy=[4,12],%输入初始状态zi=filtic(b,a,zy,zx),%计算等效初始条件[yn,zf]=filter(b,a,xn,zi),%迭代计算输出和后段等效初始条件2时域经典法用时域经典法求解差分方程:先求齐次解；再将激励信号代入方程右端化简得自由项,根据自由项形式求特解；然后根据边界条件求完全解[3].用时域经典法求解例1的基本步骤如下.（1）求齐次解.特征方程为,
3、可算出.高阶特征根可用MATLAB的roots函数计算.齐次解为（2）求方程的特解.将代入差分方程右端得自由项为当时,特解可设为,代入差分方程求得.（3）利用边界条件求完全解.当n=0时迭代求出,当n≥1时,完全解的形式为选择求完全解系数的边界条件可参考文[4]选.根据边界条件求得.注意完全解的表达式只适于特解成立的n取值范围,其他点要用及其延迟表示,如果其值符合表达式则可合并处理.差分方程的完全解为4MATLAB没有专用的差分方程求解函数,但可调用maple符号运算工具箱中的rsolve函数实现[5],格式为y=maple('rsolve({eq
4、us,inis},y(n))'),其中:equs为差分方程表达式,inis为边界条件,y(n)为差分方程中的输出函数式.rsolve的其他格式可通过mhelprsolve命令了解.在MATLAB中用时域经典法求解例1中的全响应和单位样值响应的程序如下.clc;clear;formatcompact;yn=maple('rsolve({y(n)-3/4*y(n-1)+1/8*y(n-2)=(3/4)^n+1/3*(3/4)^(n-1),y(0)=5/2,y(-1)=4},y(n))'),hn=maple('rsolve({y(n)-3/4*y(n-1
5、)+1/8*y(n-2)=0,y(0)=1,y(1)=13/12},y(n))'),3双零法根据双零响应的定义,按时域经典法的求解步骤可分别求出零输入响应和零状态响应.理解了双零法的求解原理和步骤,实际计算可调用rsolve函数实现.yzi=maple('rsolve({y(n)-3/4*y(n-1)+1/8*y(n-2)=0,y(-1)=4,y(-2)=12},y(n))'),yzs=maple('rsolve({y(n)-3/4*y(n-1)+1/8*y(n-2)=(3/4)^n+1/3*(3/4)^(n-1),y(0)=1,y(-1)=0},
6、y(n))'),4变换域法设差分方程的一般形式为.对差分方程两边取单边z变换,并利用z变换的位移公式得整理成形式有可以看出,由差分方程可直接写出和,系统函数,将系统函数进行逆z变换可得单位样值响应.由差分方程的初始状态可算出,由激励信号的初始状态可算出,将激励信号进行z变换可得,求解z域代数方程可得输出信号的象函数对输出象函数进行逆z变换可得输出信号的原函数.利用z变换求解差分方程各响应的步骤可归纳如下:（1）根据差分方程直接写出、和,的逆变换即为单位样值响应；（2）根据激励信号算出,如激励不是因果序列则还要算出前M个初始状态值；（3）根据差分方程
7、的初始状态和激励信号的初始状态算出和；（4）在z域求解代数方程得输出象函数,的逆变换即为全响应；（5）分析响应象函数的极点来源及在z平面中的位置,确定自由响应与强迫响应,或瞬态响应与4稳态响应；（6）根据零输入响应和零状态响应的定义,在z域求解双零响应的象函数,对双零响应的象函数进行逆z变换,得零输入响应和零状态响应.用变换域法求解例1的基本过程如下.根据差分方程直接写出,.系统函数的极点为.对激励信号进行z变换得.激励象函数的极点为3/4.根据差分方程的初始状态算出.根据激励信号的初始状态算出.对z域代数方程求解,得全响应的象函数.进行逆z变换得
8、全响应为其中,与系统函数的极点对应的是自由响应;与激励象函数的极点对应的是强迫响应.的极点都在z平面的单位圆内故都是瞬态响
