差分方程 matlab,差分方程的解法分析及MATLAB实现(程序)

差分方程的解法分析及MATLAB 实现(程序)

摘自：张登奇,彭仕玉.差分方程的解法分析及其MATLAB 实现[J]. 湖南理工学院学报.2014(03) 引言

线性常系数差分方程是描述线性时不变离散时间系统的数学模型,求解差分方程是分析离散时间系统的重要内容.在《信号与系统》课程中介绍的求解方法主要有迭代法、时域经典法、双零法和变换域法[1].

1 迭代法

例1 已知离散系统的差分方程为)1(3

1)()2(81)1(43)(-+=-+--n x n x n y n y n y ,激励信号为)()4

3()(n u n x n =,初始状态为21)2(4)1(=-=-y y ，.求系统响应. 根据激励信号和初始状态,手工依次迭代可算出24

59)1(,25)0(==y y . 利用MATLAB 中的filter 函数实现迭代过程的m 程序如下:

clc;clear;format compact;

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)求齐次解.特征方程为081432=+-αα,可算出4

1 , 2121==αα.高阶特征根可用MATLAB 的roots 函数计算.齐次解为. 0 , )4

1()21()(21≥+=n C C n y n n h (2)求方程的特解.将)()4

3()(n u n x n =代入差分方程右端得自由项为 ?????≥?==-?+-1,)4

3(9130 ,1)1()43(31)()43(1n n n u n u n n n 当1≥n 时,特解可设为n p D n y )4

3()(=,代入差分方程求得213=D . (3)利用边界条件求完全解.当n =0时迭代求出25)0(=y ,当n ≥1时,完全解的形式为 ,)4

3(213 )41()21()(21n n n C C n y ?++=选择求完全解系数的边界条件可参考文[4]选)1(),0(-y y .根据边界条件求得35,31721=-=C C .注意完全解的表达式只适于特解成立的n 取值范围,其他点要用

)(n δ及其延迟表示,如果其值符合表达式则可合并处理.差分方程的完全解为