欢迎光临散文网 会员登陆 & 注册

挺有用的常微分方程(五)

2023-07-16 20:48 作者:不能吃的大鱼  | 我要投稿

呀呼!一些烦心的事终于结束了,现在可以静下心来写一些东西喽!

近半年都没有更新……属实是因为这半年来太忙了,学习上的事太多,都要处理,心力交瘁,就没再多写。现在好了,终于放假了,所以可以更新啦!

上次把之前积攒的一篇专栏补全了点东西之后就发了出去,介绍了有关常系数齐次线性微分方程的一些相关内容。也许有一些和大家之前接触过的部分微分方程的概念和内容有比较大的出入,不知道各位能接受多少。希望大家都能理解吧!

现在,我们就要继续将一维线性微分方程继续推及,研究非齐次线性微分方程。


Chapter  Three  常系数线性方程 

3.3  稳定多项式

我们已经掌握了基本的解决线性微分方程的一般方法,还对解的形式和性质有了一些初步的认知。从我们之前的讨论当中,不难得出,常系数线性微分方程的解一般都与指数函数关系密切。因此,指数函数的性质对于解出来的解函数而言,有着重要的作用。

我们都知道,对于指数函数e%5E%7Bax%7D而言,如果a为正数,则x趋于正无穷时,函数无界;趋于负无穷,函数趋于0。反之则相反。一般而言,从微分方程的角度来说,需要应用的很多场合函数的变量都是非负的,这说明函数在非负半轴的性质十分重要。因此,指数上的系数a对于我们解出来的微分方程的性质就起到了关键作用。

对于常系数线性微分方程而言,其解一般具有形式:

x%3D%5Csum_%7B%E5%8D%95%E6%A0%B9%7D%20c_ie%5E%7B%5Clambda%20_i%20t%7D%2B%5Csum_%7B%E9%87%8D%E6%A0%B9%7D%20%5Cbigg%20(%20%5Csum_%7Bj%3D1%7D%5Em%20%20c_j%20t%5E%7Bj-1%7De%5E%7B%5Clambda%20_k%20t%7D%5Cbigg)

因此,各个特征值%5Clambda%20_k就是影响解函数的性质的关键。

我们知道,对于任意的解函数,都应有:

%5Clambda%20_k%3D%5Calpha%20_k%2Bi%5Cbeta%20_k%20

的形式,这是一个复数。而由Euler公式,我们能够得到:

e%5E%7B%5Clambda%20_k%20t%7D%3De%5E%7B%5Calpha%20_%20k%20t%2Bi%5Cbeta%20_k%20t%7D%3De%5E%7B%5Calpha%20_k%20t%7D(%5Ccos%20%5Cbeta%20_k%20t%2Bi%5Csin%20%5Cbeta%20_k%20t)

这说明,对于解函数是否有界,起到决定作用的是特征值的实部——%5Calpha%20_k

如果%5Cforall%20%5Calpha%20_k%EF%BC%9C%200,则称方程:

p(%5Ctext%20D)x%3D0

中的多项式p为稳定多项式

稳定多项式的概念对于计算而言是比较重要的,但是我们只关注它的一些性质即可,如果有需要的小伙伴可以自行去阅读相关的计算著作。

由于我们知道,指数函数的增长速度是相当快的,以至于一般的任何函数都不能与之比较。因此,也不难理解,即使是对于重根部分,存在幂函数和指数函数相乘的形式,只要满足稳定条件,函数最终也是趋近于0的。

因此,不难得出以下结论:

%5Cexists%20M%2C%5Cgamma%20%EF%BC%9E0%EF%BC%8C%5Cforall%20%5Cvarphi%20(t)%EF%BC%8C%7C%5Cvarphi%20(t)%7C%EF%BC%9CMe%5E%7B-%5Cgamma%20%20t%7D%5Cquad%20(p(%5Ctext%20D)%5Cvarphi%20%3D0%2Ct%5Cge%200)

证明十分容易。取:

%5Czeta%20%20%3D%5Csup%20%5C%7B%20-%5Calpha%20_k%5C%7D%2B%5Cdelta%2B%5Cvarepsilon%20%3D%5Ceta%2B%20%5Cdelta%20%2B%5Cvarepsilon%20%20%20%5Cquad%20(%5Cdelta%2C%5Cvarepsilon%20%20%20%EF%BC%9E0)

则显然有:

%5Cforall%20%5Calpha%20_k%EF%BC%8Ce%5E%7B%5Calpha%20_k%20t%7D%EF%BC%9Ce%5E%7B-%5Ceta%20%20t%7D%5Cquad%20(t%5Cge%200)

故而可以推知:

%5Cforall%20%5Calpha%20_k%EF%BC%8Ct%5Eme%5E%7B%5Calpha%20_k%20t%7D%EF%BC%9Ct%5Eme%5E%7B-%5Ceta%20%20t%7D%5Cquad%20(t%5Cge%200%2Cm%5Cle%20n)

由于:

t%5Eme%5E%7B-%5Cdelta%20%20t%7D%5Crightarrow%200%5Cquad%20(t%5Crightarrow%20%2B%E2%88%9E)

因此在t≥0时,解函数中的任意一部分都有界M_s,于是就有:

进一步,我们得到:

t%5Eme%5E%7B-%5Czeta%20%20t%7D%3D%20t%5Eme%5E%7B-%5Cdelta%20t%7D%5Ccdot%20e%5E%7B-(%5Ceta%2B%5Cvarepsilon%20)%20t%7D%EF%BC%9CM_se%5E%7B-(%5Ceta%20%2B%5Cvarepsilon%20)%20t%7D

令:

%5Cgamma%20%3D%5Ceta%20%2B%5Cvarepsilon%20

于是得到:

%7C%5Cvarphi%20(t)%7C%EF%BC%9C%5Cbigg(%5Csum_%7B%E5%8D%95%E6%A0%B9%7D%20%7Cc_i%7CM_s%2B%5Csum_%7B%E9%87%8D%E6%A0%B9%7D%20%7Cc_j%7CM_s%5Cbigg)e%5E%7B-%5Cgamma%20t%7D%3DMe%5E%7B-%5Cgamma%20t%7D

这可以用来估计函数被控制的程度。

对于稳定性的判别法,我们有一些简单的结论。比如说:

实二次多项式:

p(x)%3Dx%5E2%2Bax%2Bb

是稳定的,当且仅当a%2Cb%EF%BC%9E0

进而可知:

实多项式:

p(x)%3Dx%5En%2Ba_1x%5E%7Bn-1%7D%2B%5Ccdots%2Ba_n

是稳定的,则可以得到%5Cforall%20a_i%EF%BC%9E0

命题1,可以用代数基本定理来考虑;注意这是个必要条件。)

我们可以简要研究一下对于稍微高次的多项式,其稳定的充分必要条件是什么。

以三次实多项式:

p(x)%3Dx%5E3%2Ba_1x%5E2%2Ba_2x%2Ba_3

为例,按代数基本定理,我们可以在复数域范围内将其拆分为:

p(x)%3D(x%2B%5Calpha%20_1)(x%2B%5Calpha_2)(x%2B%5Calpha_3)%5Cquad%20(%5Calpha_1%2C%5Calpha%20%20_2%2C%5Calpha_3%5Cin%20%5Cmathbf%20C)

则有:

%5Cbegin%20%7Bequation%7D%0A%5Cbegin%20%7Baligned%7D%0A%26a_1%3D%5Calpha_1%2B%5Calpha_2%2B%5Calpha_3%5C%5C%0A%26a_2%3D%5Calpha_1%5Calpha_2%2B%5Calpha_2%5Calpha_3%2B%5Calpha_3%5Calpha_1%5C%5C%0A%26a_3%3D%5Calpha_1%5Calpha_2%5Calpha_3%0A%5Cend%20%7Baligned%7D%0A%5Cend%20%7Bequation%7D

根据代数基本定理的结论,我们很容易知道,这三个根中,要么是只有一个实数根,要么就全是实数根。不妨设:

%5Cbegin%20%7Bequation%7D%0A%5Cbegin%20%7Baligned%7D%0A%5Calpha_2%26%3D%5Cmu%2Bi%20%5Clambda%20%5C%5C%0A%5Calpha_3%26%3D%5Cmu-i%20%5Clambda%20%0A%5Cend%20%7Baligned%7D%0A%5Cend%20%7Bequation%7D%0A

则有:

%5Cbegin%20%7Bequation%7D%0A%5Cbegin%20%7Baligned%7D%0A%26a_1%3D%5Calpha_1%2B2%5Cmu%20%5C%5C%0A%26a_2%3D2%5Calpha_1%5Cmu%20%2B%5Cmu%5E2%2B%20%5Clambda%5E2%20%5C%5C%0A%26a_3%3D%5Calpha_1(%5Cmu%20%5E2%2B%5Clambda%20%5E2)%20%0A%5Cend%20%7Baligned%7D%0A%5Cend%20%7Bequation%7D

简单计算,可以得到:

%5Cbegin%20%7Bequation%7D%0A%5Cbegin%20%7Baligned%7D%0Aa_1a_2-a_3%26%3D%5Calpha_1(%5Cmu%5E2%2B%5Clambda%5E2)%2B4%5Calpha_1%5Cmu%5E2%2B2%5Calpha_1%5E2%5Cmu%2B2%5Cmu(%5Cmu%5E2%2B%5Clambda%5E2)-%5Calpha_1(%5Cmu%5E2%2B%5Clambda%5E2)%5C%5C%0A%26%3D4%5Calpha_1%5Cmu%5E2%2B2%5Calpha_1%5E2%5Cmu%2B2%5Cmu(%5Cmu%5E2%2B%5Clambda%5E2)%5C%5C%0A%26%3D%5Cmu%20%5B4%5Calpha_1%5Cmu%2B2%5Calpha_1%5E2%2B2(%5Cmu%5E2%2B%5Clambda%5E2)%5D%0A%5Cend%20%7Baligned%7D%0A%5Cend%20%7Bequation%7D

所以,基于稳定性的要求,我们又能推知:

a_1a_2%EF%BC%9Ea_3

现在我们来看,添加了这个条件之后,我们能否得到稳定性结论。

事实上,由第二个条件,我们能直接得到:

%5Cmu%20%EF%BC%9E0

这说明两个复数根的实部已经为复数。同时,又因为第一个必要条件,我们能够得到:

a_3%EF%BC%9E0

而考虑到系数的表达式,不难发现,%5Calpha_1应为复数,这就得到了稳定性。

也就是说,对于三次实多项式而言,这两个条件合并,就是其稳定的充分必要条件了。

更高阶的结论已经超出了目前需要掌握的范围,感兴趣的小伙伴们可以去阅读其他书籍或者去网络搜索相关的资料来了解。(省流:我不会~)


(你在等什么?不会是在等思考题吧~没什么好思考的,理解就好啦!常微分方程差不多就是这样子了nia~)


最後の最後に、ありがとうございました!

挺有用的常微分方程(五)的评论 (共 条)

分享到微博请遵守国家法律