数值分析计算方法实验报告.pdf

..

课程计算法实验工程实验工程类型

名称名称线性程组的数值解法

验证演示综合设计其他

指导小兵成绩√

教师

1．实验目的：

〔1〕高斯列主元消去法求解线性程组的过程

〔2〕熟悉用迭代法求解线性程组的过程

〔3〕设计出相应的算法，编制相应的函数子程序

2．实验容

分别用高斯列主元消去法,Jacobi迭代法,Gauss--Saidel迭代法,超松弛迭代法求解线性程组

21003x110





341213x5

2





1234x2

3



414913x47

3．实验环境及实验文件存档名

写出实验环境及实验文件存档名

4.实验结果及分析

输出计算结果，结果分析和小结等。

解：1.高斯列主元消去法：

%用高斯列主元消去法解实验1

%高斯列主元消元法求解线性程组Ax=b

%A为输入矩阵系数,b为程组右端系数

%程组的解保存在x变量中

.专业资料.

..

formatlong;

A=[2,10,0,-3;-3,-4,-12,13;1,2,3,-4;4,14,9,-13]

b=[10,5,-2,7]

[m,n]=size(A);

ifm~=n

error(矩阵A的行数和列数必须一样);

return;

end

ifm~=size(b)

error(b的大小必须和A的行数或A的列数一样);

return;

end

ifrank(A)~=rank([A,b])

error(A矩阵的秩和增广矩阵的秩不一样,程不存在唯一解);

return;

end

c=n+1;

A(:,c)=b;

fork=1:n-1

[r,m]=max(abs(A(k:n,k)));

m=m+k-1;

if(A(m,k)~=0)

.v.

..

if(m~=k)

A([km],:)=A([mk],:);

end

A