陳佳宏，鄧 勇

(喀什大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，新疆 喀什 844006)

設(shè)K是一個(gè)域，用Μn,m(K)(Μn(K))表示域K上的全體n×m(n×n)矩陣環(huán)。令A(yù)∈Μn(K)，B∈Μm(K)，C∈Μn,m(K)。方程

AX-XB=C

(1)

稱為Sylvester矩陣方程，其中X∈Μn,m(K)。在式(1)中，當(dāng)C=O時(shí)，即

AX-XB=O

(2)

稱為對(duì)應(yīng)于式(1)的齊次矩陣方程;特別地，當(dāng)B=-AT時(shí)，即

AX+XAT=C

(3)

稱為L(zhǎng)yapunov矩陣方程。

在控制論、信號(hào)處理、神經(jīng)網(wǎng)絡(luò)、模型降階、圖像恢復(fù)等領(lǐng)域經(jīng)常會(huì)涉及到Sylvester矩陣方程的數(shù)值求解問題。關(guān)于Sylvester矩陣方程的求解和數(shù)值計(jì)算，目前已有很多討論[1-9]。例如，文獻(xiàn)[2]討論了Lyapunov矩陣方程的公共解，給出了其無公共解的一個(gè)充要條件；文獻(xiàn)[3]利用Jordan標(biāo)準(zhǔn)形和最小多項(xiàng)式理論，討論了Lyapunov矩陣方程有非零解的充要條件，并得到其非零解空間的維數(shù)定理；文獻(xiàn)[4]利用四元數(shù)矩陣的實(shí)分解和循環(huán)矩陣的特定結(jié)構(gòu)，借助Kronecker積，得到了四元數(shù)體上Sylvester矩陣方程循環(huán)解的存在條件及其通解形式；文獻(xiàn)[5]給出了Sylvester矩陣方程的一種基于梯度的迭代算法；文獻(xiàn)[6]改進(jìn)了傳統(tǒng)的梯度迭代法，給出了求解Sylvester矩陣方程的松弛梯度迭代法，有效地提高了收斂速度。

本文利用Sylvester算子的性質(zhì)，給出了Sylvester矩陣方程(1)在任意域K上的可解性判定，并介紹求其多項(xiàng)式解的一般形式的一種新方法。

1 Kronecker乘積

1.1 Vec算子

顯然，Vec(X)是一個(gè)nm維列向量，并且對(duì)?X,Y∈Μn,m(K)，有Vec(X+Y)=Vec(X)+Vec(Y)。

1.2 Kronecker乘積

定義2[10]設(shè)X∈Μn,m(K)，Y∈Μp,q(K)。定義矩陣X和Y的Kronecker乘積為

式中，xij是矩陣X的(i,j)元。顯然，X?Y是(np)×(mq)矩陣。

定理1[10]設(shè)矩陣X,Y,Z的乘積XYZ有定義。于是，Vec(XYZ)=(ZT?X)Vec(Y)。

1.3 Kronecker乘積的應(yīng)用

設(shè)In表示n階單位矩陣。利用Kronecker乘積，可將Sylvester矩陣方程改寫為

(Im?A-BT?In)Vec(X)=Vec(C)

(4)

的形式。事實(shí)上，由AX-XB=C，有Vec(AX-XB)=Vec(C)或Vec(AX)-Vec(XB)=Vec(C)。因此，Vec(AXIm)-Vec(InXB)=Vec(C)。再由定理1可得

(Im?A)Vec(X)-(BT?In)Vec(X)=Vec(C)或(Im?A-BT?In)Vec(X)=Vec(C)。

2 Sylvester算子

設(shè)K是任意一個(gè)域，多項(xiàng)式f,g∈K[x]。用gcd(f,g)和res(f,g)分別表示f和g的最大公因式和結(jié)式；用pA表示方陣A的特征多項(xiàng)式；mA表示A的最小多項(xiàng)式。

定義3 設(shè)A∈Μn(K)，B∈Μm(K)。線性變換

ψ(A,B):Μn,m(K)→Μn,m(K),TAT-TB

稱為關(guān)于矩陣A和B的Sylvester算子。

記φA:Μn,m(K)→Μn,m(K),TAT和TB

定義4 稱C(A,B)={T∈Μn,m(K)|AT=TB}為矩陣A和B的Sylvester(或中心化子)空間[11]。

(i)f(φA)=φf(A);

證明對(duì)自然數(shù)k和矩陣T∈Μn,m(K)，分別有

(iii)ψ(T·f(B))=ψ(T)f(B)

成立，所以定理2正確。

成立。

3 引理

記號(hào)同前，引入以下引理：

引理1 設(shè)A∈Μn(K)，B∈Μm(K)。下列陳述等價(jià):

①算子ψ(A,B)=ψ是同構(gòu)的;

②C(A,B)={O};

③在K[x]中，pA和pB無公共的素因式;

④pA和pB的結(jié)式res(pA,pB)≠0;

⑤mA和mB的結(jié)式res(mA,mB)≠0。

證明參見文獻(xiàn)[12]，在此從略。

引理2 式(1)的通解可表示為其特解和式(2)的通解之和。

證明設(shè)X和X0是式(1)的通解和特解，即

于是，有A(X-X0)-(X-X0)B=O，即(X-X0)是式(2)的通解。故式(1)的通解可表示為其特解和式(2)的通解之和。

引理3[13]設(shè)f∈K[x]是一個(gè)n次多項(xiàng)式，x是域K上的未定元。于是

1)對(duì)所有0≤k≤n，都有Δkf(x)∈K[x];

2)k!Δkf(x)=f(k)(x)，其中f(k)(x)是f(x)的k階導(dǎo)數(shù)。

4 Sylvester矩陣方程的可解性

4.1 可解性判定

設(shè)K是一個(gè)代數(shù)閉域，A∈Μn(K)，B∈Μm(K)。用σ(A)表示矩陣A的所有特征值集合(譜)。文獻(xiàn)[3]指出，C(A,B)={O}?σ(A)∩σ(B)=?。下面，將此結(jié)論推廣到任意域K上，即

定理3 Sylvester矩陣方程(1)有唯一解?pA和pB無公共的素因式。

證明首先，式(1)有唯一解?關(guān)于A和B的Sylvester算子同構(gòu)[11-12]。然后，利用引理1即可得證。

推論2pA和pB無公共的素因式?mA和mB無公共的素因式。

推論3 式(1)的解空間是一個(gè)仿射空間，并且維數(shù)等于dimKC(A,B)。

推論4 式(1)有唯一解?式(2)有唯一解并且為零。

4.2 多項(xiàng)式解

證明因μ=η-ξ，故有η=ξ+μ。在上述f(x+y)的表達(dá)式中，用μ代替x，ξ代替y后即可推出結(jié)果。

定理5 設(shè)A∈Μn(K)，B∈Μm(K)，ψ是關(guān)于A和B的Sylvester算子。若矩陣A和B的特征多項(xiàng)式無公共的素因式，則

(a)ψ是一個(gè)同構(gòu)。

證明(a)由引理1即可得證。

于是，有