国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

解線性方程組的預(yù)條件AOR迭代法分析

2017-01-03 18:19:53趙秋霞
課程教育研究·下 2016年11期
關(guān)鍵詞:線性方程組

趙秋霞

【摘要】自然科學(xué)的諸多領(lǐng)域的許多問(wèn)題最終都轉(zhuǎn)化為大型線性方程組的求解,而這些方程組的求解一般采用迭代法。 對(duì)迭代法而言,當(dāng)?shù)仃嚨淖V半徑小于1時(shí),譜半徑越小其收斂速度越快,有效降低迭代矩陣譜半徑的方法就是對(duì)線性方程組本身進(jìn)行預(yù)處理。因此預(yù)條件方法成為一個(gè)熱點(diǎn)問(wèn)題。本文對(duì)幾個(gè)預(yù)條件AOR迭代法進(jìn)行程序?qū)崿F(xiàn), 并對(duì)結(jié)果進(jìn)行分析。

【關(guān)鍵詞】線性方程組 ?迭代解法 ?預(yù)條件方法 ?AOR方法 ? 譜半徑

【中圖分類(lèi)號(hào)】G642 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 【文獻(xiàn)標(biāo)識(shí)碼】A ? ? ?【文章編號(hào)】2095-3089(2016)11-0194-01

Analysis of Preconditioned AOR Iterative Methods in Solving Linear Systems

ZHAO Qiu-Xia.

(Department of Applied Mathematics, Yuncheng University, Yuncheng, 044000)

【Abstrct】The solutions of many problems in the field of natural science are eventually turned into the solutions of large linear systems. Generally, the linear systems are solved by iterative methods. The smaller it is, the faster the method converges for iterative method when the spectral radius of the iterative matrix is smaller than 1. The effective method to decrease the spectral radius of iterative matrices is to precondition linear systems. Therefore, the study of the preconditioned methods is a hot topic. In this thesis, some preconditioned AOR iterative methods are implemented by computer programs and the results are compared with basic iterative methods.

【Keywords】 Linear systems; Iterative solution; Preconditioned methods; AOR method; ?Spectral radius.

1.引言

數(shù)學(xué)、物理、流體力學(xué)和工程技術(shù)等領(lǐng)域中許多問(wèn)題的解決, 最終都轉(zhuǎn)化為大型線性方程組

(1.1)

的求解[1,2,3], 其中非奇異, 是已知的, ?是未知的. 對(duì)大型線性方程組的求解通常采用迭代解法,常用的包括AOR迭代法,SOR迭代法等[2]. 用迭代法求解線性方程組時(shí)要考慮迭代法的收斂性, 通常用迭代矩陣的譜半徑來(lái)判斷[2,3,4,5]. 當(dāng)?shù)ㄊ諗繒r(shí), 還要考察其收斂速度, 其收斂速度越快越好,為了減小迭代矩陣的譜半徑, 常對(duì)原方程組進(jìn)行預(yù)條件處理, 此類(lèi)方法稱為預(yù)條件方法.

預(yù)條件方法的主要思想是在方程組(1.1)兩端同時(shí)左乘一個(gè)非奇異矩陣 (稱為預(yù)條件因子), 將方程組(1.1)轉(zhuǎn)化為同解方程組

, ? ? ? ? ? ? ? ? ? ? ? ?(1.2)

使得將上述迭代法應(yīng)用于方程組(1.2)時(shí)迭代矩陣有較小的譜半徑.

文[3]的作者提出了下面兩種預(yù)條件因子,并在理論上證明了預(yù)條件的有效性,本文對(duì)在這兩種預(yù)條件因子作用下的預(yù)條件AOR迭代法進(jìn)行程序?qū)崿F(xiàn)和算例分析.

1),相應(yīng)的預(yù)條件AOR迭代法的迭代矩陣記為.

2), 其中為正實(shí)數(shù), 相應(yīng)的預(yù)條件AOR迭代法的迭代矩陣記為.

2. 預(yù)條件AOR迭代法的程序?qū)崿F(xiàn)與算例分析

對(duì)算例2.1[2]進(jìn)行MATLAB編程實(shí)現(xiàn). 取初值, 用AOR迭代法和預(yù)條件AOR迭代法均得到方程組(2.1)的精確解

記分別為用基本AOR迭代法, 預(yù)條件因子為的預(yù)條件AOR迭代法所用的迭代次數(shù), 見(jiàn)表2.1及2.2.

算例2.1[2]

考慮線性方程組:

(2.1)

其中系數(shù)矩陣A是非奇異且不可約M-矩陣[3].

表2.1:預(yù)條件因子為的情形

加速因子 松弛因子

0.3 0.5 0.3 0.6819 54 0.6727 50

0.7 0.9 1 0.3575 24 0.3037 17

0.8 1 1 0.2578 19 0.2054 14

1 1 1 0.1636 16 0.1378 10

1.2 1.2 1 0.2000 11 0.3115 12

1.7 1.7 1 0.7000 43 0.8848 116

1.9 1.9 1 0.9000 145 1.1120 溢出

表2.2:預(yù)條件因子為 的情形

加速因子 松弛因子

0.3 0.5 0.4 0.6819 54 0.6549 46

0.7 0.9 0.9 0.3575 24 0.2104 15

0.8 1 0.9 0.2578 19 0.1064 ? 12

1 1 0.9 0.1636 16 0.0483 9

1.2 1.2 0.9 0.2000 11 0.2208 12

1.7 1.7 0.9 0.7000 43 0.7075 127

1.8 1.8 0.9 0.8000 69 0.8053 溢出

分析表中數(shù)據(jù),當(dāng) 在 取值時(shí),AOR迭代法和兩種預(yù)條件AOR迭代法在 時(shí)效果最好;當(dāng) 取值大于1時(shí),且在1附近時(shí)效果更好,但隨著 越來(lái)越大,迭代法的收斂性降低, 甚至計(jì)算過(guò)程中發(fā)生溢出。

參考文獻(xiàn)

[1]A. Hadjimos. Accelerated overelaxation method. Math. Comput. 32 (1978) 149 – 157.

[2] 金一慶, 陳越, 王冬梅. 數(shù)值方法(第2版). 北京:機(jī)械工業(yè)出版社. 2000年2月.

[3] M.J. Wu, L. Wang, Y.Z. Song. Preconditioned AOR iterative methods for linear systems. Appl. Numer. Math. 57 (2007) 672–685.

[4]L. Wang, Y.Z. Song. Preconditioned AOR iterative methods for M-matrices. J. Comput. Appl. Math. 226 (2008) 114–124.

[5] 程光輝,黃廷祝,成孝予. 解線性方程組的預(yù)條件Guass-Seidel型迭代法. 應(yīng)用數(shù)學(xué)和力學(xué),第27卷,第9期,2006年9月15日出版. 1000-0887 (2006) 09-1117-05.

猜你喜歡
線性方程組
一類(lèi)整系數(shù)齊次線性方程組的整數(shù)解存在性問(wèn)題
求解非線性方程組的Newton迭代與Newton-Kazcmarz迭代的吸引域
H-矩陣線性方程組的一類(lèi)預(yù)條件并行多分裂SOR迭代法
Cramer法則推論的幾個(gè)應(yīng)用
求解單調(diào)非線性方程組的非精確正則化牛頓法及其局部收斂性
線性方程組解的判別
線性方程組解的逆向問(wèn)題的一種解法分析
保護(hù)私有信息的一般線性方程組計(jì)算協(xié)議
關(guān)于兩個(gè)線性方程組同解條件的再思考
基于Matlab實(shí)現(xiàn)線性方程組的迭代解法
东城区| 纳雍县| 锡林浩特市| 新沂市| 龙口市| 万盛区| 青岛市| 祁连县| 四子王旗| 邳州市| 金门县| 塔城市| 长葛市| 新晃| 崇明县| 蒙阴县| 普宁市| 宁乡县| 郁南县| 岑溪市| 招远市| 寻乌县| 浦江县| 肇州县| 大庆市| 略阳县| 湘潭市| 托里县| 获嘉县| 宁蒗| 开封县| 祁连县| 大埔县| 门源| 商都县| 郴州市| 承德县| 舟山市| 额尔古纳市| 德格县| 德钦县|