王靜

摘要：指出線性方程組原始矩陣的條件數(shù)與算法無關(guān)，因而用該條件數(shù)衡量直接法穩(wěn)定性不科學(xué)，提出算法的求解矩陣與算法直接相關(guān)，因此用求解矩陣的條件數(shù)反映算法穩(wěn)定性更為合理的新見解。

關(guān)鍵詞：直接法；算法；條件數(shù)；穩(wěn)定性

中圖分類號(hào)：TP301文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1009-3044(2012)16-3940-02

The New Opinion of Stability Analysis for Direct Methods

WANG Jing

(Teaching Research Office of Computer, North-Sichuan Medical Institute, Nanchong 637000, China)

Abstract: This paper points out the condition number of original coefficient matrix of linear equations has nothing to do with the algo? rithm itself, so it is not scientific by using this condition number to measure the stability of direct methods. Because the solving matrix is as? sociated with the algorithm itself directly, we present a new opinion that using the condition number of the solving matrix to reflect the sta? bility of algorithm is more reasonable.

Key words: direct methods; algorithm; condition number; stability

求解線性方程組有兩類方法：直接法、迭代法。其中直接法是指在不考慮計(jì)算誤差的前提下，通過確定的計(jì)算次數(shù)就能算出理論解的所有算法。包括各種Gauss消元法、改進(jìn)平方根法、追趕法和正交變換法。直接法不存在收斂問題，而僅考慮方程組的穩(wěn)定性問題?，F(xiàn)行的教材文獻(xiàn)分析算法穩(wěn)定性的方法是通過計(jì)算原始矩陣的條件數(shù)cond(A)來分析算法的穩(wěn)定性，公式如下：

1）分析方程組右端項(xiàng)擾動(dòng)對(duì)解的影響

假設(shè)δA=0，設(shè)δb引起δx，有

算法的穩(wěn)定性是選擇算法的關(guān)鍵，在所有軟件包中找不到Gauss消元法，就是因?yàn)樵撍惴ú环€(wěn)定。分析算法的穩(wěn)定性一定要與算法相關(guān)，用與算法無關(guān)的原始系數(shù)矩陣的條件數(shù)反映算法的穩(wěn)定性是不科學(xué)的，只有用與算法直接相關(guān)的求解矩陣的條件數(shù)衡量算法的穩(wěn)定性才更合理。改進(jìn)平方根法會(huì)擴(kuò)大矩陣條件數(shù)，因此穩(wěn)定性不好，而正交變換法更穩(wěn)定，如果原始矩陣的條件數(shù)是病態(tài)的，通常采用正交變換法。

[1]法捷耶夫,法捷耶娃.線代數(shù)計(jì)算方法[M].上海:上?？萍汲霭嫔?1965.

[2]馮康.數(shù)值計(jì)算方法[M].北京:國防工業(yè)出版社,1978.

[3]張池平.計(jì)算方法[M].北京:科學(xué)出版社,2006.

[4]孫文瑜,杜其奎,陳金如.計(jì)算方法[M].北京:科學(xué)出版社,2007.

[5]凌永祥,陳明逵.計(jì)算方法教程[M].西安:西安交通大學(xué)出版社,2005.