matlab用列主元高斯消去法解方程组,matlab用列主元高斯消去法解方程组 写出理论知识、设计思路、算法步骤(或流程图);写出程序清单(加上必要的注释);写出程序运行操作过程与输

来源:学生作业帮助网 编辑:作业帮 时间:2024/11/06 04:32:20
matlab用列主元高斯消去法解方程组,matlab用列主元高斯消去法解方程组     写出理论知识、设计思路、算法步骤(或流程图);写出程序清单(加上必要的注释);写出程序运行操作过程与输
xXkSW+;$dX.*8m?r+VpCoFiEƚDc#8Un)ݳ|/= =B̷2ɻ猃Cc\21l*Z-m ZJѡYy'nJ XsiP(?pKSZFlJ̀-XdZ6J,-팼3/uN9XvEXKkA~c cO`K$O0@%&#D#Д-8>fuDq_F8{BxxB3;23a gX;dx2`y_( ٠aDža粍N+|q//P(8N3r;96~.8xУ8hG<ITsio(ep|kb.02ii/;貀WhiBK#Z>>Gu4L\ɧK:s8Wm˵<[ƪ奲VOi7%UX~^u aMJ!HaV]<_Q=+pPL'=E,>.|a~rMmG[+Pf@bi<VTY1*pJ B9/.8:vT sz!Zvh)D7Î6%eZ(-h.?B Qo%ʼ*b[C@HCCpB(ED `0h_FآQ>0'; QZYg1Ө@PոۧW'쨘@!O2uδ |`16 ‹JRmYp /U-npDiE(-*3Vbnظ'x"=ilH<.D ե.m`4unJrk6̺& /9m}Y63zLmw s/#ڶ)G@ ]$LѺ.yTX2966.`?a[ëΕ4S;Sb. ډ'4j(sv@Iwdٜ#=*:K\UHa7(;^(x?4 =D?C^JUJ}eWcV_=l.uPito.꩝A]4` #l:iR@)]k A+-ZM*(f^ծGqLѨ)QL A;q^PI{Ad*L jp"ěX1sxbprQ'GC@Hx&xQ1>5O#ު<>0&>$:BO-Q:j-<5h,0l@U^OUo=րJxK 0 C_4=ugz ,ޚH$>:0FX?h t5ՎQ>ޅ 2Vߢ?"BCQF?m=V3L:fluuuw톶sٙ4w/1mz/Ә]x w0ۖ'&4

matlab用列主元高斯消去法解方程组,matlab用列主元高斯消去法解方程组 写出理论知识、设计思路、算法步骤(或流程图);写出程序清单(加上必要的注释);写出程序运行操作过程与输
matlab用列主元高斯消去法解方程组,
matlab用列主元高斯消去法解方程组     
写出理论知识、设计思路、算法步骤(或流程图);
写出程序清单(加上必要的注释);
写出程序运行操作过程与输出结果;
写出对计算过程与结果的分析(如误差分析,收敛性,稳定性,计算量,方法比较等);各设计的优缺点(如特色,自己最满意之处,需改进的地方等)


要求有上面那些内容,

matlab用列主元高斯消去法解方程组,matlab用列主元高斯消去法解方程组 写出理论知识、设计思路、算法步骤(或流程图);写出程序清单(加上必要的注释);写出程序运行操作过程与输
一、
实验名称:
线性方程组高斯消去法.

二、
实验目的:
进一步熟悉理解
Guass
消元法解法思路,提高
matlab
编程能力.

三、
实验要求:
已知线性方程矩阵,
利用软件求解线性方程组的解.

四、
实验原理


消元过程:


0
)
0
(
11

a
,令乘数
)
0
(
11
)
0
(
1
1
/
a
a
m
i
i


,做(消去第
i
个方程组的
i
x

操作
1
i
m
×第
1
个方程
+

i
个方程(
i=2
,
3
,
.n


则第
i
个方程变为
1
)
1
(
2
)
1
(
2
...
i
n
in
i
b
x
a
x
a




这样消去第
2
,
3
,
.
.
.
,
n
个方程的变元
i
x
后.原线性方程组变
为:



















)
1
(
)
1
(
2
)
1
(
2
)
1
(
2
)
1
(
2
2
)
1
(
22
)
0
(
1
)
0
(
1
1
)
0
(
11
...

.


.


...

...
n
n
nn
n
n
n
n
n
b
x
a
x
a
b
x
a
x
a
b
x
a
x
a

这样就完成了第
1
步消元.

回代过程:

在最后的一方程中解出
n
x
,得:
)
1
(
)
1
(
/



n
nn
n
n
n
a
b
x

再将
n
x
的值代入倒数第二个方程,
解出
1

n
x
,
依次往上反推,

可求出方程组的

其通项为
3,...1
-
n
2,
-
n
k

/
)
(
)
1
(
1
)
1
(
)
1
(









k
kk
n
k
j
j
k
kj
k
k
k
a
x
a
b
x

五、实验内容:

function
maintest2
clc
clear
all

A=[1 3 4;2 4 5;1 4 6];
%
系数矩阵

b=[1 7 6]'
%
常数项

num=length(b)
for
k=1:num-1




for
i=k+1:num








if
A(k,k)~=0












l=A(i,k)/A(k,k);







A(i,:)=A(i,:)-A(k,:).*l;







b(i)=b(i)-b(k)*l;








end














end

end

A
b
%
回代求
x
x(num)=b(num)/A(num,num);
for
i=num-1:-1:1




sum=0;




for
j=i+1:num








sum=sum+A(i,j)*x(j);




end





x(i)=(b(i)-sum)/A(i,i);
end

x
End
六、实验结果:

A =




1.0000



3.0000



4.0000









0


-2.0000


-3.0000









0








0



0.5000
b =




1.0000




5.0000




7.5000
x =




16


-25



15