理论教育 线性微分系统的反射函数理论及应用

线性微分系统的反射函数理论及应用

时间:2023-11-25 理论教育 版权反馈
【摘要】:考虑线性微分系统这里P(t)为连续矩阵函数.设X(t)为式(3.1.1)的基解矩阵,则其通解为φ(t;τ,x)=X(t)X-1(τ)x,则其反射函数为F(t,x)=φ(-t;t,x)=X(-t)X-1(t)x.定义3.1.1 称F(t)=X(-t)X-1(t)为系统(3.1.1)的反射矩阵.由反射函数的性质可推得反射矩阵F(t)具有如下性质:1) F(-t)F(t)≡F(0)=E. (3.1.2)

线性微分系统的反射函数理论及应用

考虑线性微分系统

这里Pt)为连续矩阵函数.

Xt)为式(3.1.1)的基解矩阵,则其通解为φtτx)=XtX-1τx,则其反射函数为F(t,x)=φ(-tt,x)=X(-t)X-1t)x.

定义3.1.1Ft)=X(-tX-1t)为系统(3.1.1)的反射矩阵.

由反射函数的性质可推得反射矩阵Ft)具有如下性质:

1) F(-tFt)≡F(0)=E. (3.1.2)

2) Ft)为系统(3.1.1)的反射矩阵的充要条件为Ft)为Cauchy问题

的解.

由Liouville公式得

3)对任一非奇异且满足式(3.1.2)的矩阵Ft),存在一个线性系统

Ft)为反射矩阵.

显然Ftx+x=c为式(3.1.4)的解,从而式(3.1.4)任一解x(t)有xt)+x(-t)=2x(0)成立.系统(3.1.4)称为以F(t)为反射矩阵的简单线性系统.系统(3.1.1)为简单线性系统的必要条件为Pt)≡P(-tFt.

4)以Ft)为反射矩阵的线性系统均可表示为

这里Rt)为任意连续矩阵函数,tR.

定理3.1.1 对于微分系统(3.1.1),Pt)为连续矩阵,对充分小t,存在唯一微分系统

与系统(3.1.1)具有相同的反射矩阵F(t),且Q(t)=Q(-t.若系统(3.1.1)与系统(3.1.5)皆为2ω-周期系统,则它们具有同一Poincaré映射.

Ft)为系统(3.1.1)的反射矩阵.

下面来求解

F′t)+FtQt)+QtFt)=0. (3.1.6)

由于方程Ft)-λE=0的特征根,当t=0时为1,则对t<<1,矩阵Ft)与-Ft)没有公共根,由参考文献[27]得方程(3.1.6)存在唯一解Qt.下证Qt)=Q(-t.事实上,在式(3.1.6)中令t→-t,则有

F′(-t)+F(-tQ(-t)+Q(-tF(-t)=0. (3.1.7)

FtF(-t)=E.F′tF(-t)=FtF′(-t.将式(3.1.7)两端左乘F(t)得

FtF′(-t)+FtF(-tQ(-t)+FtQ(-tF(-t)=0,

F′tF(-t)+Q(-t)+FtQ(-tF(-t)=0.

两边右乘F-1(-t)得

F′t)+Q(-tF-1(-t)+FtQ(-t)=0,

F′t)+Q(-tFt)+FtQ(-t)=0.(www.daowen.com)

由于方程(3.1.6)的解唯一,从而Qt)≡Q(-t.由此及式(3.1.6)知,微分系统(3.1.5)以Ft)为反射矩阵.

定理3.1.2Pt)为R上的连续矩阵函数,并满足

则系统(3.1.1)具有反射矩阵

显然

F(0)=EFtF(-t)=E.

下面验证Ft)满足式(3.1.3)即可.

事实上,由条件Pt)与978-7-111-47659-7-Chapter03-9.jpg可交换,则P(-t)与978-7-111-47659-7-Chapter03-10.jpg也可交换,从而978-7-111-47659-7-Chapter03-11.jpg可交换,则

F′t)=-(Pt)+P(-t))Ft),FtPt)=PtFt),从而有

Ft)为系统(3.1.1)的反射矩阵.(这里P:=P(-t))

定理3.1.3Pt+2ω)=Pt),且Pt)在R上连续,并满足关系式(3.1.8),则系统(3.1.1)的Poincaré映射为

并且当特征方程

的所有特征根λk,有Reλk<0时,系统(3.1.1)渐近稳定,若存在一个特征根λk,有Reλk>0时,该系统不稳定.

证 由定理3.1.2得

又由于

Ft)=X(-tX-1t),

F(-ω)=XωX-1(-ω)=X(-ωX(2ωX-1(-ω),

F(-ω)相似于X(2ω.由此即可得定理的结论成立.

注3.1.1 由参考文献[116]知,系统(3.1.1)的反射矩阵可表示为Ft)=eStx,其中St)=ln(X(-tX-1t))且S(-t)+St)=0.

定理3.1.4 微分系统(3.1.1)具有形如Ft,x)=eAtx的反射函数的充要

条件为

-2P(0)exp(-2P(0)t)+exp(-2P(0)tPt)+P(-t)exp(-2P(0)t)≡0,且此时Ftx)=e-2P(0)tx.

该结论直接由反射函数的定义可证.

推论3.1.1 在定理3.1.4的条件下,若Pt+2ω)=Pt),则微分系统

(3.1.1)的任一解xt)为2ω-周期解的充要条件为x(-ω)=ζ为代数方程

[e2P(0)ω-E]ζ=0

的解.若det(e2P(0)ω-E)≠0,则x=0为系统(3.1.1)的唯一2ω-周期解,且当特征方程2P(0)ω-λE=0的特征根λk满足Reλk<0时,渐近稳定.若存在λk使得Reλk>0,则该周期解不稳定.

该推论可由定理3.1.3和定理3.1.4直接推得.

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈