Jacobi迭代法解線性方程組

2021-02-18 數值分析與有限元編程

當線性方程組的規模比較大時,採用高斯消元法需要太多時間。這時就要採用迭代法求解方程組了。高斯消元法是一個O(n^3)的浮點運算的有限序列,在經過有限步計算之後理論上得到的是精確解(無捨入誤差時)。而迭代法在經過有限步迭代之後一般不產生精確解,迭代法在計算過程中逐漸減小誤差,當誤差小於容許值時停止迭代計算。方程組的係數矩陣是嚴格對角佔優矩陣時,迭代總是收斂的。

●Jacobi迭代法

對於方程組3u+v=5,u+2v=5,將其改寫為如下的形式

由於方程組的係數矩陣是嚴格對角佔優矩陣時,迭代一定收斂。使用初值[u0,v0]=[0,0]開始迭代,以下是迭代過程:

繼續迭代過程最終會收斂到解[1,2].這個迭代過程就是Jacobi迭代。

對於方程組u+2v=5,3u+v=5,由於方程組的係數矩陣不是嚴格對角佔優矩陣時,因此迭代不收斂。來看迭代過程:

設D表示係數矩陣A 的主對角部分,L表示A的主對角線下方部分,U表示A的主對角線上方部分。則A=D+L+U,AX=b可改寫為

對於上面的方程組3u+v=5,u+2v=5,寫成矩陣形式

迭代格式為

這與之前的迭代格式是一致的。

Fortran原始碼

☆☆☆  往期相關  ☆☆☆

嚴格對角佔優矩陣

高斯消去法解線性方程組及MATLAB實現

高斯消去法的算法改進

相關焦點

  • Scipy入門(求解方程組2)
    一.非線性方程組     在前面我們討論的大多數是關於未知數的線性方程組
  • 基於MATLAB求解三類方程組
    基於MATLAB求解三類方程組第1類 線性方程組線性方程組是指各個方程變量的最高次冪為1次的方程組。線性方程組主要有定解方程組、不定解方程組、超定方程組和奇異方程組。線性方程組的通解為齊次線性方程組的基礎解系加上線性方程組的一個特解。
  • 中國突破億階線性方程組高性能求解
    對於與浪潮聯合開發的億階線性方程組MIC加速方案,太原理工大學博士生導師,中國計算數學學會常務理事、副秘書長李明給出了上述評價。「毫不誇張的講,大部分科學與工程問題都要歸結為一個線性方程組的求解問題。」李明教授一開頭說道。他介紹說,線性方程組的求解問題是一個十分古老的問題,早在中國古代數學專著《九章算術》第八章就詳細記載了一次線性方程組的求解方法。
  • 線性方程組練習題
    線性代數的一個核心問題就是求解線性方程組:線性方程組分為兩類,齊次線性方程組,非齊次線性方程組。
  • 齊次線性方程組求解
    求解齊次線性方程組是線性代數的重要內容,也是線性代數的核心問題之一,要求同學們必需掌握。
  • 《二元一次方程組》教案
    《二元一次方程組》教案一、教學目標【知識與技能】掌握二元一次方程與二元一次方程組的概念,並了解它們的解,能正確地找出二元一次方程組的解。【過程與方法】通過類比學習、自主探究、合作交流的過程,提升類比學習的能力,樹立探究的意識。
  • 馬千子: 二元一次方程組
    迎來的第一課就是二元一次方程組。          二元一次方程是含有兩個未知數,並且所含未知數的項的次數都是一的方程叫做二元一次方程。         二元一次方程組是共含有兩個未知數的兩個一次方程所組成的方程。其中一次方程是指包含一元一次方程和二元一次方。         適合一個二元一次方程的一組未知數的值,叫做這個兩元一次方程的一個解。
  • 三元一次方程組的解法-3
    1那麼,這個方程組裡含有三個未知數,每個方程中含未知數的項的數都是1,並且一共有三個方程,像這樣的方程組叫做三元一次方程組。那麼怎麼去求解三元一次方程組呢?在求解二元一次方程組時,可以利用代入法或可以利用代入法或加減法消去一個末知末數,化成一元一次方程求解,那麼三元一次方程組呢?我們依然利用消元法,將三元一次方程組轉化為二元一次方程組,再利用消元法將二元一次方程組轉化為一元一次方程解。
  • 巧解二元一次方程組:整體代入法在二元一次方程組中的應用
    該方法在減少多項式項數,降低多項式結構複雜程度或簡化方程組等方面能起到獨到作用。下面我們就來看幾道應用整體代入法解方程組的例子。例1、解方程組:x-y-1=0,①4(x-y)-y=5。②分析:如果我們用常規的方法,把方程②的括號去掉,合併整理之後用代入消元或加減消元法也不是不可以,但略顯麻煩。
  • 量子計算機可解方程組
    量子計算機,可解線性方程組?這已不是神話。中科大6月8日發布消息,該校科學家日前在國際上首次成功實現用量子計算機求解線性方程組的實驗。  何為線性方程組?數學家認為,各個方程關於未知量均為一次的方程組,如2元1次方程組,即為線性方程組。對線性方程組的研究,中國比歐洲至少早1500年。如今,線性方程組已廣泛應用於數值計算、信號處理、經濟學、計算機科學等科學、工程領域。與日常生活緊密相關的氣象預報,就需要建立並求解百萬變量的線性方程組,來實現對溫度、氣壓、溼度等參數的模擬和預測。
  • 二元二次方程組(一)
    把兩個或者兩個以上的方程放在一起組合成方程組,一般用大括號括起來,如把 x²+2x=0和x+y=0這兩個方程組合起來便成了一個方程組。現在,我們來看一下這個方程組的構成,方程(1)是一個一元二次方程,方程(2)是一個二元一次方程,兩個方程共有兩個未知數且未知數的最高次數為2,我們把這樣的由兩個未知數的一個二次方程和一個次數不超過二次的方程所組成的方程組稱為二元二次方程組。例  判斷下列方程組是否為二元二次方程組
  • 三元一次方程組的基礎解法
    三元一次方程組的基本解法跟二元一次方程組相同,最本質的就是需要消元。
  • 《二元一次方程組》說課稿
    一、說教材首先談談我對教材的理解,《二元一次方程組》是人教版初中數學七年級下冊第八章第一節的內容,本節課的內容是二元一次方程組的概念以及二元一次方程組的解。在此之前學習了一元一次方程和解方程的步驟,為本節課打下了良好的基礎。學了本節課為後面的解二元一次方程的方法做下鋪墊。因此本節課有著承上啟下的作用。二、說學情接下來談談學生的實際情況。
  • 高等代數之線性方程組
    極大線性無關組的概念就是後面線性空間的基,也是齊次線性方程組解空間的基礎解系。從向量組的角度來看就是向量組的秩是一樣的。線性方程組解的理論跟幾何上平面的位置關係有密切的聯繫:二個平面的重合、平行、相交。
  • 二元二次方程組解法(一)
    上一期,我們知道了長什麼樣的方程組是二元二次方程組以及如何判定。從這一期開始,我們要探討如何解二元二次方程組。
  • 數學之美:牛頓-拉夫遜迭代法原理及其實現
    牛頓迭代如何迭代?直接看數學公式描述如何迭代不直觀,先來看動圖就很容易理解牛頓迭代法為什麼叫迭代法以及怎樣迭代的:牛頓迭代法是原理是根據一個初始點由前面描述知道,牛頓迭代法是用來近似求解方程的,這裡有兩個點需要說明:迭代法非常適合計算機編程實現,實際上計算機編程對於牛頓迭代法廣為應用
  • 方程組和行列式的淵源
    01理論二元線性方程組和二階行列式方法一:兩個未知數,兩個方程組,消元即可解出未知數方法二:行列式得出方程組的解
  • 初一數學,二元一次方程組的解法
    二元一次方程,是指有兩個未知數,並且未知數的指數是一次的方程,由兩個二元一次方程組成的,就是二元一次方程組。解二元一次方程組的思路,主要是消元,就是把未知數變為一個,其中,代入消元法和加減消元法是最常用的解題方法。
  • 數學之美:牛頓-拉夫遜迭代法原理及其實現
    牛頓迭代如何迭代?直接看數學公式描述如何迭代不直觀,先來看動圖就很容易理解牛頓迭代法為什麼叫迭代法以及怎樣迭代的:牛頓迭代法是原理是根據一個初始點在該點做切線,切線與X軸相交得出下一個迭代點的坐標,再在處做切線,依次類推,直到求得滿足精度的近似解為止。