实验4解线性方程组的直接法.doc

《实验4解线性方程组的直接法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西华数学与计算机学院上机实践报告课程名称：数值计算方法B年级：2009级上机实践成绩：指导教师：严常龙姓名：上机实践名称：解线性方程组的直接法学号：上机实践日期：2011.12.2上机实践编号：4上机实践时间：8:00~9:40一、目的1．通过本实验，加深对解线性方程组的直接法——高斯列主消元法和LU分解法的构造过程的理解；2．能将各种方法的算法描述正确地改编为程序并上机实现；3．熟悉并掌握各种方法的适用对象及优缺点，学会针对不同问题选择不同方法。二、内容与设计思想自选线性方程组，编制一个程序，分别用高斯列主元消元法和LU分解法求解。然后完成作业：分

2、别用高斯列主元消元法和LU分解法求解下面的方程组（以下方程组供选择）。三、使用环境操作系统：windowsxp软件平台：vc++6.0四、核心代码及调试过程#include#include#includevoidmain(){inti,j,k,p,N;floata[50][51],t,x[50];floatmax,temp;printf("请输入方程未知数的个数N:");scanf("%d",&N);printf("请输入方程的增广矩阵的数据,数据之间按回车键隔开!");for(i=0

3、;i=0){max=a[k][k];}e

4、lse{max=a[k][k];}for(i=k;i

5、

6、(max<-a[i][k])){max=0;if(a[i][k]>=0){max=a[i][k];}else{max=-a[i][k];}}}for(i=k;i

7、

8、(max=-a[i][k])){p=i;break;}}for(j=0;j

9、选主列的元。for(i=k+1;i=0;k--){for(i=k-1;i>=0;i--){t=(a[i][k]/a[k][k]);for(j=0;j

10、i][j]=a[i][j]-t*a[k][j];}}}printf("增广矩阵的系数对角线化得到的矩阵：");for(i=0;i#include

11、ath.h>#includefloata[50][51],x[50],u[50][51];voidmain(){inti,j,k,p,N;floatsum,t;printf("请输入方程未知数的个数N:");scanf("%d",&N);//*************************************输入矩阵a的数据。printf("请输入方程系数矩阵a的数据,数据之间按回车键隔开!");for(i=0;i

12、