标题:
关于差分方程
[打印本页]
作者:
yuyang911220
时间:
2016-7-4 16:05
标题:
关于差分方程
在
数学
上,
递推关系(recurrencerelation)
,也就是
差分
方程(difference equation)
,是一种
递推地
定义一个序列的方程式:序列的每一项目是定义为前一项的函数。
举个例子(
户口调查映射
(logisticmap)):
某些简单定义的递推关系式可能拥有非常复杂的(混乱的)行为,且有时候在数学的领域称为
非线性分析
中被物理学家和数学家所研究。
解一个递推关系式的意思,就是得到一种
n
的非递回函数33
[
编辑
]常系数线性齐次递推关系式
线性
字眼的意思是序列的每一项目是被定义为前一项的一种线性函数。系数和常数可能视 n 而定,甚至是非线性地。
一种特别的情况是当系数并不依照 n 而定。
齐次
意思为关系的常数项为零。
为了要得到线性递回唯一的解,必须有一些起始条件,就是序列的第一个数字无法依照该序列的其他数字而定时,且必须设定为某些数值。
[
编辑
]解线性递推关系式递推关系式的解通常是由系统的方法中找出来,通常借由使用
生成函数
(generatingfunction) (
形式幂级数
(formalpowerseries))或借由观察
r
n
是一种对
r
的特定数值之解的事实。
二阶递推关系式的形式:
我们拥有解为
r
n
:
两边除以
r
n
− 2
我们可以得到:
这就是递推关系式的
特征方程
。解出
r
可获得两个根(roots)
λ1,λ2
,且如果两个根是不同的,我们可得到解为
而如果两个根是相同的(当
A
2+4
B
=0),我们得到
C
和
D
都是常数。
换句话说,将这种
a
n
=
A
a
n
− 1 +
B
形式的方程式,用 2 代入 n 后,就得到上述的
r
2 =
A
r
+
B
。常数 "C" 和 "D" 可以从"边界条件(side conditions)" 中得到,通常会像是“已知
a
0 =
c
1
,
a
1 =
c
2
”。
[
编辑
]范例:斐波那契数(Fibonacci Number)
斐波那契数
是使用一种线性递推关系式来定义:
设若:
当n趋于无限大之极限值存在,则其值为
= Φ
恰为
黄金分割
之ㄧ值,1.618....,另一值则为0.618....,两值互为倒数,也就是说1.618....分之1=0.618....,反之亦然。
起始条件为:
因此,斐波那契数的序列为:
0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89 ...[
编辑
]常系数非齐次线性递推关系对于常系数非齐次线性递推关系,我们可以用
待定系数法
来求出它的一个特解,而它的通解就是这个特解与对应的齐次递推关系的通解的和。也可以使用
迭代法
求解,但只能得到确切的数值解,不能直接以解析式作答,该方法可利用计算机求解。
[
编辑
]时域经典法求解一般情况下,常系数线性差分方程可以写作:
则对应的齐次方程形式为:
则特征方程为:
当特征根非重根时,齐次解为:
当特征根为重根时,若
α1
为特征方程的
K
重根,齐次解为:
特解
y
p
(
n
) =
D
(
n
)
的形式由激励函数
x
(
n
)
的形式决定。
一般情况,当激励函数x(n)代入方程。
方程右方出现
n
k
的形式,则特解选择
当方程右方出现
a
n
的形式,则特解选择
当a不是特征根时
y
p
(
n
) =
A
a
n
当a是特征根时
y
p
(
n
) =(
A
1
n
+
A
0)
a
n
当a为r重根时
将特解带入原方程,求出待定系数。根据边界条件,可求出齐次节待定系数。
[
编辑
]例子我们用待定系数法来解以下的常系数非齐次线性递推关系:
对应的齐次递推关系
的齐次解是:
我们猜测特解的形式为:
代入原递推关系中,我们便得到:
比较等式两端的
3
n
项的系数,可得:
比较等式两端的
n
项的系数,可得:
比较等式两端的常数项,可得:
因此原递推关系的通解为:
[
编辑
]与差分方程的关系
数值
求解
常微分方程
时,经常会遇到递归关系。例如,求解如下
初值问题
时
如采用
欧拉法
和步长
h
,可以通过如下递归关系计算
y
0=
y
(
t
0)
,
y
1 =
y
(
t
0 +
h
),
y
2 =
y
(
t
0 + 2
h
),...
y
n
+ 1 =
y
n
+
h
f
(
t
n
,
y
n
).
线性一阶微分方程组可以用
离散化
条目中介绍的方法解析地精确离散化。
欢迎光临 电子技术论坛_中国专业的电子工程师学习交流社区-中电网技术论坛 (http://bbs.eccn.com/)
Powered by Discuz! 7.0.0