数值分析里面经常会涉及到用MATLAB程序实现用列主元消去法分别解方程组Ax=b
具体的方法和代码以如下方程(3x3矩阵)为例进行说明:
用列主元消去法分别解方程组Ax=b,用MATLAB程序实现:
(1)
1、 实现该方程的解的MATLAB代码可以分为两种,一种是入门级别的,只是简单地计算出这道题即可,第二种是一种通用的代码,可以实现很多3x3矩阵的方程解,写好以后只需要改不同矩阵里的元素即可算出相应的解,需要建立在对MATLAB比较熟悉的基础上,具体如下:
第一种代码实现—入门级:
A=[3.01,6.03,1.99,;1.27,4.16,-1.23,;0.987,-4.81,9.34]
A1=[3.01,6.03,1.99,1;1.27,4.16,-1.23,1;0.987,-4.81,9.34,1]
B1=A1(1,1:4)
C1=A1(2,1:4)
D1=A1(3,1:4)
E1=-1.27/3.01*B1+C1
F1=-0.987/3.01*B1+D1
p1=E1(1,2)
q1=F1(1,2)
if (abs(p1)>=abs(q1))
a1=p1
a2=q1
FF1=E1
EE1=F1
else
a1=q1
a2=p1
FF1=F1
EE1=E1
end
G1=-a2/a1*FF1+EE1
H1=[E11;FF1;G1]
J1=H1(1:3,1:3)
b1=H1(1:3,4)
x1=J1\b1
第二种代码实现如下—熟练通用级:
A=[3.01,6.03,1.99,;1.27,4.16,-1.23,;0.987,-4.81,9.34]
A1=[3.01,6.03,1.99,1;1.27,4.16,-1.23,1;0.987,-4.81,9.34,1]
B1=A1(1,1:4)
C1=A1(2,1:4)
D1=A1(3,1:4)
f1=A1(1,1)
f2=A1(2,1)
f3=A1(3,1)<