9512.net
甜梦文库
当前位置:首页 >> 数学 >>

Gauss列主元消去法实验


<数值计算方法>实验报告
1.实验名称 实验 2 2.实验题目 用 Gauss 列主元消去法求解线性方程组。 Gauss 列主元消去法

?0.001x1 ? 2.000 x 2 ? 3.000 x3 ? 1.000 ? ??1.000 x1 ? 3.712 x 2 ? 4.623x3 ? 2.000 ??2.000 x1 ? 1.072 x 2 ? 5.643 x3 ? 3.000 ?

3.实验目的 加深自己对 Gauss 列主元消去法的理解和认识,并且通过做实验 或做练习来加强自己 Gauss 列主元消去法的掌握,学会并灵活运用 Gauss 列主元消去法来求解方程组。
4.基础理论

-------Gauss 列主元消去法
1.Gauss 列主元消去法的基本思想是:在进行第 k(k=1,2,...,n-1)步 消元时,从第 k 列的

akk 及以下的各元素中选取绝对值最大的元素,然 akk 的位置上,再进行消元。 akk (k=1,2,...,n-1)的绝对值很

后通过行变换将它交换到主元素

2.Gauss 列主元消去法的优点:当

小时,用 Gauss 列主元消去法来求解方程组时,可以避免所的数值结果

产生较大误差或失真。

5.实验环境 实验系统:Win 7 实验平台:VisualC++语言 6.实验过程 写出算法→编写程序→计算结果
Gauss 列元消去法的算法 Input: 方 程 组 未 知 量 的 个 数 n; 增 广 矩 阵 A ? (aij ) ? ? A1, A2,..., An ? , 其 中
T

i=1,2,…,n; j=1,2,…,n+1 Output:方程组的解 x1,x2,…,xn,或失败信息。 1. for i ←1ton-1 do;

2. temp←| aii |; 3. 4. 5. 6. 8. 9. 10. 11. 12. 13. p←I; for j←i+1 to n do if

| a ji | >temp

then

p←j; end end if temp=0 |return end if p≠I then then False;

14. 15. end

Ap

? Ai
← a ji

;//i,p 两行交换

//列选主元 16. 17. 18. 19. for j←i+1 to n do

Aj ? m ji * Ai m ji

/ aii

;

Aj ← Aj ? mji * Ai ;//消元
end

7.实验结果
原方程组的解为:

X1=-0.490396 ,

x2=-0.051035 ,

x3=0.367520

8.附录

程序清单

#include<iostream.h> #include"stdio.h" #include"math.h" void main ( )

{

int n=3,i,j,k,p; double

A[10][10]={{0.001,2.000,3.000,1.000},{-1.000,3.712,4.623,2.000},{-2.0 00,1.072,5.643,3.000}},temp,m,x[100]; for(i=0;i<n;i++) { //选主元 temp=fabs(A[i][i]); p=i; for(k=i+1;k<n;k++) if(fabs(A[k][i])>temp) { temp=fabs(A[k][i]); p=k;

} if(temp==0) { printf("\n 无法求解:"); return; } if(p!=i) for(j=0;j<n+1;j++) { temp=A[i][j]; A[i][j]=A[p][j]; A[p][j]=temp;

}

//消元 for(k=i+1;k<n;k++) { m=A[k][i]/A[i][i];

for(j=i+1;j<=n;j++) A[k][j]=A[k][j]-m*A[i][j]; }

} //回代 for(i=n-1;i>=0;i--) { x[i]=A[i][n]; for(j=i+1;j<n;j++) x[i]=x[i]-A[i][j]*x[j]; x[i]=x[i]/A[i][i]; } printf("\nx=\n"); for(i=0;i<n;i++) printf("%f \n",x[i]); }


赞助商链接

更多相关文章:
试验四 Gauuss列主元消去法解线性方程组实验报告
课程实验报告课程名称 姓名 实验名称 实验目的及要求 班级 学号 实验日期 实验成绩 试验四 Gauuss 列主元消去法解线性方程组 1、熟悉 Gauss 列主元消去法,编出...
Matlab程序Gauss列主元消去法
Matlab程序Gauss列主元消去法_IT/计算机_专业资料。.Gauss 列主元消去法求解线性...数学实验“线性方程组高... 5页 免费 Gauss列主元消去法 3页 免费©...
数值分析实验,用程序实现列主元消去法解方程组
数值分析实验报告 《数值分析》实验报告实验序号:实验二 学号: 姓名: 任课教师: 专业班级:) 实验名称: 列主元消去法解方程组 1、 实验目的:用列主元 Gauss ...
Gauss列主元消去法实验
<数值计算方法>实验报告 1.实验名称 实验 2 2.实验题目 用 Gauss 列主元消去法求解线性方程组。 Gauss 列主元消去法 ?0.001x1 ? 2.000 x 2 ? 3.000 ...
计算方法实验报告 高斯列主元素消去法
再设计程序代码的过程中,充 5 分考验了我们的独立思考能力与探究精神, 并且在这次实验过程中, 我更加了解 Gauss 列主 元素消去法的实质,为以后的学习留下一个...
高斯列主元消去法实验
《数值分析》实验报告实验 2 Gauss 列主元消去法实验题目:用 Gauss 列主消去法求解线性方程组 0.001*X1+2.000*X2+3.000*X3=1.000 -1.000*X1+3.217*...
相对Gauss列主元消去法
实验题目 5 相对 Gauss 列主元消去法摘要由一般线性方程组在使用 Gauss 消去法求解时,从求解过程中可以清楚地看到,若 (k (k akk ?1) = 0 ,必须施以行...
2.3高斯列主元消去法
2.3 高斯列主元消去法解线性方程组一:问题的提出我们都知道,高斯列主元素消去...高斯列主元消去法实验 2页 1下载券 Gauss列主元消去法 4页 免费 列主元...
Gauss列主元消去法
2. 编制列主元 Gauss 消去法函数(子程序),并用于解 取初始向量 x ( 0) ...六、调试过程及实验结果 >> A=[12,-3,3;-18,3,-1;1,1,1]; >> b...
实验二:Gauss列主元消去法
列主元消去法 暂无评价 5页 免费如要投诉违规内容,请到百度文库投诉中心;如要提出功能问题或意见建议,请点击此处进行反馈。 实验二:Gauss列主元消去法 隐藏>> ...
更多相关标签:

All rights reserved Powered by 甜梦文库 9512.net

copyright ©right 2010-2021。
甜梦文库内容来自网络,如有侵犯请联系客服。zhit325@126.com|网站地图