本專題將講述以多種方式求解線性方程組,也作為本人在《數值計算與優化》課程中學到知識的總結與具體**實現。主要用到資料為《數值計算方法(第3版)》。
解線性代數組是科學研究與工程計算中經常遇到的問題。 此專題討論以下n階線性方程組⎧⎩
⎨⎪⎪⎪
⎪⎪⎪a
11x1+
a12x2
+⋯+a
1nxn
=b1a
21x1+
a22x2
+⋯+a
2nxn
=b2⋯
an1x
1+an
2x2+
⋯+an
nxn=
bn的解法,這裡假定係數行列式不為零。若用矩陣表示,則上述方程組可表示為1a
x=b
目前,計算機上常用的解線性方程組的方法大致可分為「直接法」與「迭代法」兩大類。「直接法」是指那些在沒有捨入誤差影響的條件下經有限步四則運算可求得準確解的方法,而「迭代法」則是一種逐次逼近的方法
2。 本專題主要介紹如下幾類求解方法:
《數值計算方法(第3版)》page 34, para 1. ↩
《數值計算方法(第3版)》page 34, para 3. ↩
數值計算解線性方程組
高斯消元法 function x solvegauss a,b n size a,1 x zeros n,1 for j 1 n 1 for i j 1 n mul a i,j a j,j a i,a i,mul a j,b i b i mul b j endendfor i n 1 1 sum 0...
MATLAB 線性方程組求解
clc,clear all close all 高斯消去法 a 2 3 4 3 5 2 4 3 30 線性方程組的係數矩陣 b 6 5 32 線性方程組的右端列向量 m,n size a 測量係數矩陣的維數 if m n fprint 線性方程組的係數矩陣非方陣 break end fprintf ...
MATLAB線性方程組求解
對於一般的,有唯一解的線性方程組,我們可以轉換成矩陣的形式 a x bax b ax b 則可以用矩陣運算求解x,即x a b 求解齊次線性方程組基礎解系的函式是null z null a 表示返回矩陣a的基礎解系組成的矩陣。z還滿足ztz i z null a,r 得出的z不滿足ztz i,但得出...