|
,bn}是M的特征值,v=(v,v,…,vn)是相应的特征向量。定义Mv=cv转化为以下方程组:bkvk+2λakpbkXlvl=cvk,k=1,2,n、 PLVL6=0一定是真的。否则,bkvk=CVK对于所有k=1,2,n、 自c起/∈ {b,b,…,bn}(参见步骤2.1),这意味着v=0,这与特征向量的定义相矛盾。因此,我们可以设置1>v=Plvl=2λ,而不丧失一般性。我们获得v=一√卑诗省- b、 a√卑诗省- b一√bnc公司- bn公司.让cn>cn-1> ····>c>0是M的特征值。定义c:=诊断(ci)i=1,。。。,n、 Q:=cj公司- bi公司i、 j=1,2,。。。,nand Q:=AB1/2Q.(23)2.3 M=QCQ-1、在步骤2.2中,Q的列是与特征值c、c、…、,中国。对应于不同特征值的特征向量是线性独立的,因此Q是非奇异的。我们得到特征分解M=QCQ-1.2.4 1>Q=2λ>。这源自步骤2.2,其中我们假设Q中包含的每个特征向量和为2λ。3、定义:=1+2λXkak√bk-1> 0.我们让M1/2:=Q诊断(√ci)i=1,。。。,nQ公司-并用eM1表示/2T=Q diag(e√ciT)i=1,。。。,nQ公司-M1/2T的矩阵指数。定义:=A-1.M1/2eM1/2T- 我+ B1/2层eM1/2T+I,其中I表示单位矩阵。n个常微分方程组f=Mf的通解- 2AB1/21 isf(t)=eM1/2tx+eM1/2(t-t) x+2dAB-1/21,t∈ [0,T],对于x,x∈ 注册护士。要看到这一点,让我们∈ [0,T]和x,x∈ 注册护士。写入d=1/(1+2λ1>AB-1/21)showsd MAB-1/21=dAB1/21+2λAB1/21 1>AB-1/2= d1+d- 1.AB1/2=AB1/21。因此,f(t)=MeM1/2tx+eM1/2(T-t) x个= Mf(t)- 2d MAB-1/2=Mf(t)- 2AB1/21。仍需选择X和xin,以满足步骤1.2和1.3中的边界条件。首先,f(0)- f(T)=(eM1/2T- 一) (十)- x) 。步骤2.3,eM1/2T- I=Q诊断e√ciT公司i=1,。。。,nQ公司-1.- I=Q诊断e√ciT公司- 1.i=1,。。。,nQ公司-1是非奇异的。因此,f(0)=f(T)当且仅当x=x。设置x=x。
|