📅  最后修改于: 2020-12-22 04:56:48             🧑  作者: Mango
通过将初始条件的值放入齐次解中,可以找到差分方程为齐次线性型时的特殊解。
示例1:求解差分方程2a r -5a r-1 + 2a r-2 = 0,并找到特定的解,使得0 = 0和a 1 = 1。
解决方案:特征方程为2s 2 -5s + 2 = 0
(2s-1)(s-2)= 0
s = 和2。
因此,方程的齐次解为
a r(h) = C 1 + C 2 .2 r ……….方程(i)
将r = 0和r = 1放在等式(i)中,我们得到
a 0 = C 1 + C 2 = 0 ………..方程(a)
1 = C 1 + 2C 2 = 1 ………….方程。(b)
求解方程(a)和(b),我们有
C 1 =- 和C 2 =
因此,特定的解决方案是
<="" alt="特殊解决方案" p="" src="https://static.javatpoint.com/tutorial/dms/images/particular-solution5.png">
示例2:求解差分方程a r -4a r-1 + 4a r-2 = 0并找到特定的解,使得0 = 0和1 = 6。
解决方案:特征方程为
s 2 -4s + 4 = 0或(s-2) 2 = 0 s = 2,2
因此,方程的齐次解为
a r(n) =(C 1 + C 2 r).2 r …………..等式(i)
将r = 0和r = 1放在等式(i)中,我们得到
0 =(C 1 + 0)0.2 0 = 1∴C1 = 1
1 =(C 1 + C 2)0.2 = 6∴C1 + C 2 =3⇒C2 = 2
因此,特定的解决方案是
a r(P) =(1 + 2r).2 r 。
示例3:求解满足条件a 0 = 0和a 1 = 2的差分方程9a r -6a r-1 + a r-2 = 0。
解决方案:特征方程为
因此,方程的齐次解为
a r(h) =(C 1 + C 2 r)。 ………方程(i)
将r = 0和r = 1放在等式(i)中,我们得到
0 = C_1 = 0
a 1 =(C 1 + C 2 )。 = 2。 ∴C1 + C 2 =6⇒C2 = 6
因此,特定的解决方案是
r(P) = 6r。 。
有两种方法可以找到非齐次线性差分方程的特定解。这些如下:
1.不确定系数法:该方法用于查找非齐次线性差分方程的特定解,其RHS项R(n)由特殊形式的项组成。
在这种方法中,首先,我们根据包含一些未知常数的R(n)的类型,假定特定解的一般形式,必须确定这些常数。然后根据差分方程,我们将确定精确解。
表中显示了为特定形式的R(n)假定的特定解的一般形式,以找到确切的解。
Form of R (n) | General form to be assumed |
---|---|
Z, here z is constant | A |
Zr, here z is constant | Zr |
P (r), a polynomial of degree n | A0 rn+A1 rn-1+⋯..An |
Zr. P (r), here P(r) is a polynomial of the nth degree in r. Z is a constant. | [A0 rn+A1 rn-1+⋯..An].Zr |
例1:找到差分方程a r + 2 -3a r + 1 + 2a r = Z r ……..方程(i)的特定解
Z是一些常数。
解:解的一般形式是= A. Z r
现在将这个解决方案放在等式(i)的LHS上,我们得到
= AZ r + 2 -3AZ r + 1 + 2AZ r =(Z 2 -3Z + 2)AZ r ………方程(ii)
将等式(ii)与等式(i)的RHS相等,我们得到
(Z 2 -3Z + 2)A = 1
A = (Z≠1,Z≠2)
因此,特定的解决方案是
例2:找到差分方程a r + 2 -5a r + 1 + 6a r = 5 r的特定解………….方程(i)
解:让我们假设解的一般形式为A. 5 r 。
现在找到A的值,将此解放在等式(i)的LHS上,则变为
=A。5 r + 2 -5.A5 r + 1 + 6.A5 r
= 25A。 5 r -25A.5 r + 6A.5 r
= 6A.5 r …………方程(ii)
将公式(ii)等同于公式(i)的RHS,我们得到
a ><=""
因此,差分方程的特定解为= 0.5河
例3:求出差分方程a r + 2 + a r + 1 + a r = r.2 r ………方程(i)的特定解
解:让我们假设解的一般形式=(A 0 + A 1 r)。 2 ^ r
现在,将这些解放在等式(i)的LHS中,我们得到
= 2 r + 2 [A 0 + A 1 (r + 2)] + 2 r + 1 [A 0 + A 1 (r + 1)] + 2 r (A 0 + A 1 r)
= 4. 2 r (A 0 + A 1 r + 2A 1 )+2.2 r (A 0 + A 1 r + A 1 )+2 r (A 0 + A 1 r)
= r。 2 r (7A 1 )+2 r (7A 0 + 10A 1 )…………(ii)
将等式(ii)与等式(i)的RHS相等,我们得到
7A 1 = 1∴A 1 =
7A 0 + 10A 1 = 0∴A 0 =
因此,特定的解决方案是
2. E和∆运算符方法:
运算符E的定义: f(x)上E的运算符表示给函数x的值增加一个值。 E的运算是将(x + h)放在函数的任何x位置。这里,h是增量量。所以Ef(x)= f(x + h)
此处,E在f(x)上进行运算,因此,E是称为移位运算符的符号。
运算符∆的定义:运算∆是两步运算。
首先,在函数x由恒定递增,然后前被从后面即减去
δf(x)= +="" f(x="" h)-f(x)<="" p="">
定理1:证明E≅1+ ∆。
证明: Δ在f(x)上的运算分为两个步骤。首先,在函数增加x的值。因此,只要在f(x)中存在x,就将x + h(这里h是恒定增量)放进去,这意味着E在f(x)上的运算,即,
f(x + h)= Ef(x)。
第二,从第一步获得的值中减去原始函数,因此
δf(x)= ef(x)-1f(x)="(E-1)f(x)
因此,∆在f(x)上的运算等效于(E-1)在f(x)上的运算。
因此,我们有
e≅1+ p="" ∆。<="">
定理2:证明E n f(x)= f(x + nh)。
证明:我们知道E f(x)= f(x + h)
现在E n f(x)= EEEE …….. n次f(x)
= E n-1 [E f(x)] = E n-1 f(x + h)
= E n-2 [E f(x + h)] = E n-2 f(x + 2h)
………………….
………………….
= Ef [x +(n-1)h] = f(x + nh)。
定理3:证明E Cf(x)= CE f(x)
证明:我们知道EC f(x)= C f(x + h)= CE f(x + h)。因此证明。
E的运算对任何常数都没有影响。因此,E对任何常数的运算将等于常数本身。
通过E和∆运算符方法,我们将找到
C 0 y n + r + C 1 y n + r-1 + C 2 y n + r-2 +⋯+ C n y n = R(n)……..等式(i)
等式(i)可以写成
C 0 E r y n + C 1 E r-1 y n + C 2 E r-2 y n +⋯+ C n y n = R(n)
(C 0 E r + C 1 E r-1 + C 2 E r-2 +⋯+ C n )y n = R(n)
将C 0 E r + C 1 E r-1 + C 2 E r-2 +⋯+ C n = P(E)
所以P(E)y n = R(n)
y y n = …….方程式(ii)
为了找到不同形式的R(n)的(ii)特定解,我们有以下几种情况。
情况1:当R(n)为常数A时。
我们知道,E在任何常数上的运算将等于常数本身,即
EA = A
因此,P(E)A =(C 0 E r + C 1 E r-1 + C 2 E r-2 +⋯+ C n )A
=(C 0 + C 1 + C 2 +⋯+ C n )A
= P(1)A
因此,使用等式(ii),(i)的特定解为
y n = ,P(1)≠0
通过将E = 1放在P(E)中来获得P(1)。
情况2:当R(n)的形式为A. Z n时,其中A和Z为常数
我们有P(E)(A. Z n )= {C 0 E r + C 1 E r-1 +⋯+ C n }(AZ n )
= A {C 0 Z r + n + C 1 Z r + n-1 +⋯+ C n Z n }
= A {C 0 Z r + C 1 Z r-1 +⋯+ C n }。 Z N
= AP(Z).Z n
要获得,P(Z)将E = Z放入P(E)
因此, ,前提是P(Z)≠0
因此,y n = ,P(Z)≠0
如果A = 1,则y n =
当P(Z)= 0时,则为等式
(i)(EZ)y n = A.Z n
为此,特定的解决方案变为A。 Z n = A。 n Z n -1
(ii)(EZ) 2 y n = A.Z n
为此,特定的解决方案成为
(iii)(EZ) 3 y n = A.Z n
为此,特定的解决方案成为等等。
情况3:当R(n)是次数为m的多项式时,n为n。
我们知道E≅1+ ∆
因此,p(e)= +="" p="" p(1="" ∆)<="">
可以将∆的上升能力扩展到∆ m
⇒ =(B 0 + B 1Δ+ B 2 +Δ⋯。+ B MΔM +⋯)
⇒ .R(N)=(B 0 + B 1Δ+ B 2 +Δ⋯。+ B MΔM +⋯).R(n)的
= b 0 R(n)+ b 1 ∆ R(n)+⋯+ b m ∆ m R(n)
所有其他较高项将为零,因为R(n)是次数为m的多项式。
因此,在这种情况下,方程式(i)的特定解为
y n = b 0 R(n)+ b 1 ∆ R(n)+⋯。+ b m ∆ m R(n)。
情况4:当R(n)的形式为R(n).Z n时,其中R(n)是次数为m的多项式,Z为某个常数
我们有E r [Z n R(n)] = Z r + n R(n + r)= Z r .Z n .E r .R(n)= Z n (ZE) r R(n)
同样,我们有
[Z n R(n)] = Z n 。(R(n))= Z n [P(Z + Z∆)] -1 .R(n)
因此,在这种情况下,方程式(i)的特定解为
y n = Z n [P(Z + Z∆)] -1 .R(n)
示例1:找到差分方程的特定解
2a r + 1 -a r = 12。
解决方案:上面的等式可以写成
(2E-1)a r = 12
具体解决方案由
r = .12
将E = 1放在等式中。特定的解决方案是r = 12
示例2:找到差分方程a r -4a r-1 + 4a r-2 = 2 r的特定解。
解决方案:上面的等式可以写成
(E 2 -4E + 4)a r = 2 r
因此,P(E)= E 2 -4E + 4 =(E-2) 2
因此,特定的解决方案由
因此,p(e)=>e≅1+>δf(x)=>δf(x)=>a>