楊 鳳

(溫州大學(xué)數(shù)理與電子信息工程學(xué)院，浙江溫州 325035)

本文考慮下列非奇異復(fù)對(duì)稱線性系統(tǒng)

其中 W ,T ∈ Rn×n均為對(duì)稱矩陣，i為虛數(shù)單位．這類線性系統(tǒng)廣泛存在于科學(xué)計(jì)算和工程應(yīng)用中，例如渦流問(wèn)題[1]、光散射成像[2]、分子動(dòng)力學(xué)和流體動(dòng)力學(xué)[3]等問(wèn)題．

對(duì)于問(wèn)題（1），當(dāng) W ,T ∈ Rn×n，且其中一個(gè)矩陣是對(duì)稱正定的，另一個(gè)矩陣是對(duì)稱半正定時(shí)，前人已經(jīng)提出了許多的迭代法．基于線性系統(tǒng)系數(shù)矩陣的 Hermite和斜 Hermite分裂，Bai等[4]提出了HSS（Hermitian and Skew-Hermitian Splitting）迭代法，此后基于此類的迭代法層出不窮．例如為了避免求解系數(shù)矩陣為斜Hermite線性方程組，Bai等[5]提出了修正的HSS迭代法；為了加快MHSS（Modified Hermitian and Skew-Hermitian Splitting）迭代法的收斂速度，Bai等[6]又提出了PMHSS（Preconditioned Modified Hermitian and Skew-Hermitian Splitting）迭代法；Wang等[7]提出了CRI（Combination Method of Real Part and Imaginary Part）迭代法，并證明了CRI迭代法迭代矩陣譜半徑的上界比PMHSS迭代法迭代矩陣譜半徑的上界更?。?/p>

1 CRI迭代法

2 CRI變型迭代法

CRI迭代法[7]的迭代格式如下：

其中α為任意大于 0的實(shí)數(shù)．在復(fù)對(duì)稱線性系統(tǒng)（1）的兩邊同時(shí)乘以(1 - i )有：由CRI迭代法迭代格式，得到了CRI變型迭代法的迭代格式

從而對(duì)?α＞0，CRI變型迭代法對(duì)于線性系統(tǒng)（1）都無(wú)條件收斂．

證明：對(duì)CRI變型迭代法迭代矩陣做適當(dāng)?shù)淖冃斡校?/p>

類似定理1的證明，設(shè)λ為矩陣φ(α)的特征值，由u表示T-1W的特征值知

當(dāng)0＜α＜1時(shí)，h(α)單調(diào)遞減；當(dāng)α＞1時(shí)，h(α)單調(diào)遞增．則α*=1時(shí)，h(α)取得極小值，將α*=1代入ρ(φ(α) ) 中，可得ρ(φ( 1 ))=，在定理 2的條件下，由推論 1知

3 數(shù)值實(shí)驗(yàn)

使用了一個(gè)復(fù)對(duì)稱線性系統(tǒng)例子，通過(guò)CRI變型迭代法和CRI迭代法進(jìn)行比較，驗(yàn)證了CRI變型迭代法求解復(fù)對(duì)稱線性系統(tǒng)（1）的有效性和可行性．在實(shí)驗(yàn)中，分別從迭代次數(shù)（IT）、計(jì)算時(shí)間（CPU）、迭代誤差（ERR）這三個(gè)角度進(jìn)行比較，選擇初始向量 x0為零向量，迭代誤，其中 xk為當(dāng)前的迭代殘量．在CRI變型迭代法和CRI迭代法的比較中，令參數(shù)α=1．

例1 復(fù)對(duì)稱不定線性系統(tǒng)（1.1）[5]的形式如下：

其中矩陣M和K是慣性和剛度矩陣，CV和τK是粘滯阻尼矩陣，ω是駕駛圓頻率，τ是阻尼系數(shù)．

令 M = I,CV= 1 0I ，K是5點(diǎn)中心差分，在一個(gè)均勻網(wǎng)格上近似于具有齊次狄利克雷邊界條件的負(fù)拉普拉斯算子的矩陣，而均勻網(wǎng)格[0,1]×[0,1]的尺寸大小矩陣K也是一個(gè)n階塊三對(duì)角矩陣，具有如下的張量積形式 ( n = m2)：

分別令ω =4,τ =1選擇m=16, 32, 64, 128時(shí)，滿足W是對(duì)稱半正定矩陣，T是對(duì)稱正定矩陣，以及本文的假設(shè)CRI變型和CRI這兩種迭代法求解例1的實(shí)驗(yàn)結(jié)果見(jiàn)表1，表 1實(shí)驗(yàn)數(shù)據(jù)表明從迭代次數(shù)（IT），計(jì)算時(shí)間（CPU）和迭代誤差（ERR）這三個(gè)角度來(lái)看，CRI變型迭代法都優(yōu)于CRI迭代法，并且對(duì)于CRI變型迭代法求解的線性系統(tǒng)維數(shù)越大，迭代次數(shù)越?。?/p>

表1 例1的數(shù)值實(shí)驗(yàn)結(jié)果

4 結(jié) 論

本文分析了CRI迭代法收斂半徑以及收斂半徑的最優(yōu)參數(shù)，在此基礎(chǔ)上進(jìn)一步提出了求解復(fù)對(duì)稱線性系統(tǒng)的CRI變型迭代法，并給出了CRI變型迭代法的收斂分析以及最優(yōu)收斂因子．在取最優(yōu)參數(shù)的條件下，驗(yàn)證了在滿足一定條件時(shí)，CRI變型迭代法比CRI迭代法的收斂半徑更小，最后通過(guò)數(shù)值實(shí)驗(yàn)結(jié)果體現(xiàn)了此結(jié)論，并驗(yàn)證了CRI變型迭代法的優(yōu)越性．