问答文章1 问答文章501 问答文章1001 问答文章1501 问答文章2001 问答文章2501 问答文章3001 问答文章3501 问答文章4001 问答文章4501 问答文章5001 问答文章5501 问答文章6001 问答文章6501 问答文章7001 问答文章7501 问答文章8001 问答文章8501 问答文章9001 问答文章9501

二阶常系数齐次线性差分方程的一般解法

发布网友 发布时间:2024-10-01 09:11

我来回答

1个回答

热心网友 时间:2024-11-10 01:51

前言:这篇文章写于2019年10月,当时在学习概率论中的贝叶斯公式时,遇到一些需要用差分模型来解决的问题,于是记录下了这篇文章。

解释术语

我们称形式为(我们在中学学过的符号)的方程为差分方程。

一个[n]阶差分方程可以写成

其中最大指数和最小指数的差是[n]。

因此,一个二阶差分方程是以下形式

如果方程可以写成

其中[a]可能依赖于[a],则称为线性方程。

如果[a]我们称方程为齐次方程。否则,它被称为非齐次方程。

如果系数[a]是常数,这意味着[a]与指数无关,我们说方程具有常数系数。

因此,本文要讨论的是一种差分方程,其形式为

或者我们将其转换为形式(你将在以下段落中了解为什么我要这样做)

其中[a](否则它将不是二阶差分方程)

如果函数有一个附加条件如[a],使得我们可以求解唯一的函数,我们称这些条件为初始条件。

现在我们可以开始求解它了。

方法1:使用几何序列

显然,[a]是满足方程的一个解。但如果我们只关注零解,那可能有点无聊,甚至无益。我们主要关心的是通解,这意味着所有可能的解都可以用它来表示。

我们首先考虑一个形式为[a]且满足方程的几何序列。假设公比为[a],那么我们有



比较系数,我们得到

所以根据维达公式,[a]和[a]是方程的解

我们可以轻松地解这个方程,得到

设[a]([a]可以是[a]或[a])。那么[a]是一个几何序列。为了不失一般性,我们令[a],那么[a]因此

由于[a],[a],然后我们有

因此

如果[a],即[a],我们可以得到

注意到[a]。因此

其中[a]

如果[a],那么[a]和[a]可以写成[a]。所以通过使用德莫弗公式,我们可以得到

其中[a]。等式意味着解是[a]的形式。同样,我们可以通过初始条件求解[a]。

如果[a],即[a],那么

因此

这意味着解是形式[a]的形式,其中[a]可以通过使用初始条件求解。

更好的计算方法

当我们得到公式[a]时,我们可以得到另一个类似的等式

因此,当[a]时,通过克莱姆法则,我们得到

因此,我们得到了相同的答案,但计算量大大减少!

方法2:使用特征方程和根

在了解了方法1之后,你可能觉得计算很复杂,容易出错(也许作者本人也有一些),最好找到一个更简单的解,以减少计算量并降低错误率。因此,这种方法应运而生。

我们首先猜测(我不知道提出这个想法的人的想法,但我们只能遵循他的步骤)解的形式为[a]。我们将其代入方程,将得到

我们假设[a],然后我们可以得到[a]。这个方程是我们在方法1中看到的方程,但现在这个方程有了新的意义,即其解直接是初始差分方程的解。我们称[a]的方程[a]的特征方程。特征方程的解称为特征根。

现在我们有一个问题:在大多数情况下,方程[a]有两个不同的根,这两个根都包含在通解中吗?你可能可以在下面的引理中找到答案:

引理:齐次方程的叠加

引理:如果[a]是差分方程[a]的解,那么[a]也是方程的解。

证明:基于条件,我们有

然后

Q.E.D.

现在又出现了一个新问题:什么是通解的形式。以方程[a]为例,它是它的通解[a]、[a]还是其他?定理1告诉我们答案。

定理1:不同特征根的解

定理:如果[a]是特征方程[a]的两个不同根(这意味着[a]),那么差分方程的通解是

其中[a]是常数,[a]

证明:基于引理,[a]是差分方程[a]的解

因此,我们需要证明的是,对于任何初始值[a]和[a],我们可以得到差分方程的唯一特解(这个陈述等价于任何方程的解都可以用通解表示)。

给定初始值[a]和[a],我们可以得到线性系统:

由于[a],系数行列式的

因此,线性系统有唯一解[a]。因此,方程的通解是[a]。

Q.E.D.

我们还有一个问题要做:如果[a]怎么办?仍然,有一个定理为我们提供了答案:

定理2:重复特征根的解

定理:如果特征根的两个根相同([a]),它们都可以表示为[a],差分函数的通解是

其中[a]

证明不难,你可以自己试试。(*^_^*)

方法3:使用生成函数

到目前为止,我们没有使用微积分方法。为了证明微积分确实非常有用,我们决定使用级数理论来解方程。

我们首先考虑一个级数(我们假设[a]存在)

其中[a]称为生成函数。

由于[a]

因此

如果[a],[a]可以分解为两个部分分数

然后我们使用泰勒展开

再次我们得到了通解的一般形式[a]

如果[a],[a]可以写成

然后使用泰勒展开

仍然,[a]类似于形式[a]。

方法4:使用线性代数

线性代数在解差分方程中也很有用。不仅我们可以用这种方法解决问题,还可以用这种方法解一些差分方程的系统。有些人可能会想,如果我们想使用线性代数,我们必须找到向量与向量之间的关系,但我们如何得到它。

现在我们考虑方程[a],如果我们考虑方程[a],我们可以将其重写为[a]。

然后我们可以用矩阵形式重写方程,即

因此,我们得到了两个向量之间的关系,如果我们表示[a],那么通过计算[a],我们可以得到方程的解。

(在理论 上,一种方法是使用归纳法找到[a]的元素规律,另一种方法是使用[a]的特征分解,即[a]([a]是非奇异的)来计算。而在实践中,当[a]足够大时,我们可以使用快速幂算法来计算结果。)

一个例子:斐波那契数列

斐波那契数列是一个数列[a],给定差分方程和初始值如下:

该方程的特征函数是

通过解这个方程,我们可以得到特征根

因此,通解是

然后我们将[a]代入方程,得到

因此

所以,我们得到特解

总结

我们介绍了四种解二阶线性齐次常数系数差分方程的方法。有趣的是,所有这些方法都可以用于解[a]阶线性差分方程(仍然,上述属性是必要的)。如果你对此感兴趣,可以尝试一下。
声明声明:本网页内容为用户发布,旨在传播知识,不代表本网认同其观点,若有侵权等问题请及时与本网联系,我们将在第一时间删除处理。E-MAIL:11247931@qq.com
锅炉切圆直径调整 为什么我打开excel后是灰色的? 有没有什么值得推荐的投资港股的QDII基金? - 知乎 win10电脑屏幕亮度调不了怎么办 考焊工证具体怎么操作 广汽传祺是丰田技术吗 广汽传祺与广汽丰田有关系吗 怎样才能让不愿写作业的孩子听话? 上初中不写作业不想学习的孩子应该如何管教? 出生在2020年农历七月十二日的宝宝怎么取名有内涵? 多哈回合谈判各方呼吁 已知函数y= e^(2x),则x=? 在WORD如何插入图片作为背景 和女朋友一起自驾游,晚上怎么住 我跟女朋友去,肯定不会住大房间啊。。。那边便宜一点的酒店房间价格怎么... 两千块钱够我和我女朋友玩一天吗就住个宾馆吃个饭 和女友出去... 有必要去很棒的宾馆或是酒店吗 广西分类垃圾房售价 我电脑win10玩cf200的fps,怎么回到win7只有100了 win10玩联盟从开局fps100多打着打着突然掉到30 总之很低又秒回100多老... 能知道05和06年坐牢的人名字吗?和他什么时候出来, 无心法师中的无心的心事被谁掏空了 怎么通过QQ号查找快手? 球球大作战内吐怎么解锁_球球大作战内吐解锁方法 脾气很倔的女孩子该怎么教育 《球球大作战》电脑版怎么向中间吐球介绍_《球球大作战》电脑版怎么向... 如何根据已有内容快速制作Word表格如何根据已有内容快速制作word表格文 ... 为什么在没有输错的情况下扣扣添加好友搜索显示无相关结果 旅游景区环境保护管理制度是什么,旅游景点责任 景区或酒店有停车场吗 CPU步进1好还是步进2好 为什么我很讨厌用那些非主流头像的人?如果有人加我,他的头像是很恶心... 键盘上的sleep关机后,怎样一键开机? ...哪个好心人给几张男生非主流头像!~网上那些一般的就不要发了!要珍 ... 能否让普通键盘实现"一键关机"? 按键盘上的哪个键能叫电脑一键开机 什么是汉密顿抑郁量表 14+7+7隔离 夏雨连连,缤越置身其中,看那一抹娇艳的红 cf官网里的道具商城的“杀敌/死亡清零卡”是把杀敌总数或死亡总数清零... 铜矿粉的活性指数是什么? CF商城里卖的是所有战绩的清零卡还是负战绩清零卡? 西递宏村二日游攻略 西递宏村电话多少 怎么给异地恋的女友送花? 异地送花怎样送? 神之墓地丧尸龙魂刀怎么合成神之墓地26d咋合丧魂龙刀 风油精弄到衣服上如何清洗? 风油精弄到衣服上怎么洗掉? 如何查看手机短信聊天记录? 怎么在手机查短信聊天?