论文部分内容阅读
求解非线性矩阵方程的问题主要是通过分析所给方程参数的性质来得到方程的解.由于Hermite正定解在实际中应用较多,所以我们只讨论此类解的情况.在现实生活中,方程的来源相当广泛,包括控制理论,动态规划,统计和椭圆型偏微分方程的差分方法求解等多个领域.关于此类方程的求解通常涉及到三个问题:(1)可解性问题,即方程有解的充分和必要条件;(2)数值求解问题,即有效的数值方法;(3)解的扰动分析.
首先,本文讨论了方程X+A*X-2A=Q(1)的可解性及解的性质,主要结论如下:
定理1矩阵方程(1)有解的充要条件是存在非奇异的矩阵W,Z,使得A=(W*W)Z,其中矩阵WQ1/2 ZQ1/2是列酉正交的,此时方程(1)有解X=W*W.
定理2矩阵方程(1)有解的充要条件是存在酉矩阵U1,V1和对角矩阵Γ>0,φ>0,使得(公式略)
其中Γ2+φ2=I.此时X=Q1/2*Γ2Q1/2是方程(1)的解.
定理3若方程(1)有解为X,A为非奇异矩阵,则(公式略)
定理4设方程(1)有解为X,A为非奇异矩阵,则(公式略)其中(~a)是数量方程(公式略)在(0,2/3)内的解.
定理5设方程(1)有解为X,A为非奇异矩阵,则(公式略)其中(^a)是数量方程(公式略)在[2/3,1]内的解.
定理6若A和Q满足条件(2.7),X是矩阵方程(1)的解,则(公式略)其中,α1,β1是方程(2.5)的两个正实根;α2,β2是方程(2.6)的两个正实根其次,通过两种不同的迭代方式,利用不动点定理分别求出了方程(1)的迭代解,并讨论了迭代解的收敛性,同时给出了迭代解的扰动界.
定理7若A和Q满足条件(2.7),则矩阵方程(1)在[β1I,β2I]上存在唯一解XL,且不可能有比XL更大的解.XL可由以下迭代得到(公式略)其中X0∈[β1I,β2I].
定理8若A和Q满足条件(2.7),则矩阵方程(1)在X∈[α2I,α1I]]上有极大解(^X)和极小解(-X),且(^X)=limn→∞(^X)n,(-X)=limn→∞(-X)n,其中(公式略)且(-X)0≤(-X)1≤…≤(-X)n≤…≤(^X)≤…(^X)n≤…≤(^X)1≤(^X)0.对方程(1)的任意解X,都有X≥(-X),即(-X)是方程(1)的所有解中的最小解.
定理9设A,(~A)=A+△A为非奇异矩阵,Q,(~Q)=Q+△Q∈Cn×n为Hermite正定矩阵.若A和Q满足条件(公式略)
最后,给出了方程的牛顿迭代算法并证明了其收敛阶。
定理11由牛顿迭代(4.2)产生的矩阵序列{Xk}二次收敛到XL,即存在C>0,使得(公式略)。