吳宇航，閻少宏，彭美葉

（華北理工大學理學院，河北唐山063009）

一類特殊反對角方程組的追趕法及其實現(xiàn)

吳宇航，閻少宏，彭美葉

（華北理工大學理學院，河北唐山063009）

反三對角方程組；非奇異矩陣；YH分解；追趕法

研究了反三對角方程組的求解問題。首先給出了反三對角矩陣A的定義，其次證明了滿足嚴格反對角占優(yōu)的反三角矩陣為非奇異矩陣，然后通過利用YH矩陣分解的方法，推導得出了反三對角方程組的追趕法，最后運用算例進行演示。

0 引言

隨著現(xiàn)代工業(yè)和科學技術的發(fā)展，線性方程組的應用出現(xiàn)在經(jīng)濟管理、工程計算等各個領域，許多應用會導出一些具有特殊結構的稀疏線性方程組的計算問題［1，2］。伴隨著這些方程組的出現(xiàn)，尋找簡便而且準確的求解方法就顯得十分重要而且具有現(xiàn)實意義。

反三對角方程組是反對角方程組中比較常見的一類，在力學、流體力學、工程學等領域有很重要的應用。事實上，因為許多實際問題通常會采用微分方程模型來描述，然后用有限元方法、差分方法等來求其數(shù)值解，這類問題的最終解決又歸結為解大型線性方程組的問題。在對角線性方程組的解法中，追趕法因其計算公式簡單，運算量和存儲量小，在科學領域中被廣泛運用。因此，在借鑒文獻［3－6］的思想后，建立求解反三對角方程組的追趕法，然后算例進行演示。

1 預備知識

定義1若方陣A＝（aij）的元素當1≤i≤n－2，1≤j≤n－i－1且3≤i≤n，n－i＋3≤j≤n時，均有aij＝0，則稱此矩陣為反三對角矩陣。

定義2若反三對角矩陣A滿足條件（1），則A為嚴格反對角占優(yōu)矩陣。

引理1對于任意階數(shù)不小于2的反三對角矩陣A，一般記：

若滿足條件（1），則矩陣A為非奇異矩陣。

證明：利用行列式初等變換把矩陣A轉(zhuǎn)化成三對角矩陣。即，

當n＝4wn＝4w＋1（w為任意大于零的整數(shù)）時，

當n＝4w＋2n＝4w＋3（w為任意大于零的整數(shù)）時，

通過三對角矩陣是非奇異矩陣的證明過程可以引證得滿足條件（1）的反三對角矩陣A為非奇異矩陣［6］（可逆矩陣）。

定理1對于非奇異的反三對角矩陣A都可以分解為如公式（2）形式的一個反三角矩陣與一個三角矩陣的乘積［7］。

證明：假設反三對角矩陣A可實現(xiàn)矩陣YH分解。利用矩陣乘法，可得

根據(jù)公式（3）可得一系列關系式如下：

由不等式（5）及條件（1）可知

即

同理可證

由a1＝ln≠0，c2＝p2ln，由公式（8）得0，同理可求得

從而由公式（3）可求出pi。

這就是說，由A的假設條件（1），完全確定了｛li｝，｛pi｝，｛mi｝，實現(xiàn)了A的YH分解。

2 反三對角線性方程組的追趕法

在實際問題中，會經(jīng)常遇到如下形式的線性方程組

這種方程組稱為反三對角方程組，簡記為Ax＝f。

在引理1的條件下有l(wèi)i≠0，即追趕法可以進行計算，故將求解方程組Ax＝f化為依次求解：

算法1

第1步：解方程組Yy＝f，即“追”過程，算法如下：

第2步：解方程組Hx＝y(tǒng)，即“趕”過程，算法如下：

其中l(wèi)i，pi，mi的計算見公式（3）。

上述算法1就是求解反三對角方程組的追趕法。又由不等式（8）的估計，即追趕法計算過程的中間變量有界，不會產(chǎn)生大的變化，可以有效地算出結果。

3 算例演示

用追趕法求解下列反三對角方程組：

根據(jù)反三對角方程組有

利用公式（3）可推導出

把以上結果代入公式（14），就可以得到關于xi（i＝1，2，…，7）的解：

4 結論

針對反三對角方程組的特點，沿用解三對角方程組時LU分解和追趕法的基本思想，證得滿足嚴格反對角占優(yōu)的反三對角矩陣是可逆的，可以分解成YH形式并推導出了反三對角方程組的追趕法。計算結果表明，此種追趕法求解n階反三對角方程組只需要O（5n－4）的運算量，可在線性時間內(nèi)完成求解工作。

［1］ 張明望．一類非單調(diào)線性互補問題的高階Dikin型仿射尺度算法［J］．數(shù)學雜志．2004，24（05）：585－590．

［2］ 邱茂路．矩陣約當標準化的一個新方法［J］．數(shù)學雜志．2001，21（02）：237－240．

［3］ 李文強，馬民．求解循環(huán)三對角方程組的追趕法［J］．科技導報．2009，27（14）：69－72．

［4］ 倪有義，蔡靜．反五對角與擬反五對角方程組的追趕法［J］．數(shù)學雜志．2014，34（01）：137－144．

［5］ 楊愛民，閻少宏，夏國坤，等．求解三對角方程組的并行追趕算法［J］．河北理工大學學報：自然科學版．2008，30（01）：107－109．

［6］ 薛正林，吳開騰．求擬反三對角線性方程組的一種數(shù)值方法［J］．內(nèi)江師范學院學報．2016，31（02）：4－7．

［7］ 關治，陸金甫．數(shù)值分析基礎［M］．北京：高等教育出版社，1998：254－256．

［8］ 李慶揚，王能超，易大義．數(shù)值分析（第5版）［M］．北京：清華大學出版社，2008：159－160．

Pursuant Method of A Kind of Special Anti－diagonal Equations and Its Realization

WU Yu－h(huán)ang，YAN Shao－h(huán)ong，PENG Mei－ye
（College of Science，North China University of Science and Technology，Tangshan Hebei 063009，China）

anti－tridiagonal equation；nonsingular matrix；YHdecomposition；pursuant method

The problem of solving the anti－tridiagonal equation is discussed．The definition of the antitridiagonal matrix Ais firstly given．Secondly，it is proved that meet strictly diagonally dominant antitriangular matrix is nonsingular，and then by taking advantage of the YHdecomposition of matrix method，it is concluded that the pursuant method of anti－tridiagonal equations．Finally，a numerical example is used to demonstrate．

O241．6

A

2095－2716（2016）04－0027－05

2016－05－03

2016－09－22