直接法解线性方程组.doc

《直接法解线性方程组.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数学与计算科学学院实验报告实验项目名称直接法解线性方程组所属课程名称数值方法A实验类型验证型实验日期2014.11.28班级信计12-学号201253100姓名成绩22一、实验概述：【实验目的】1.掌握用C语言编程实现追赶法求解三对角线性方程组；2.掌握运用高斯列主元消去法解线性方程组；3.加深对解线性方程组的直接法——高斯列主消元法和LU分解法的构造过程的理解；4.熟悉并掌握各种方法的适用对象及优缺点，学会针对不同问题选择不同方法；5.培养使用电子计算机进行科学计算和解决问题的能力。【实验原理】1.追赶法原理222.高斯列主元消去法【实验环境】1.硬件环

2、境2.软件环境（1）（2）VC++6.0二、实验内容：【实验过程】（实验步骤）1.实验步骤1）深入了解解题过程并依次写出解题算法；222）依照算法用C语言编写解题程序；3）上机时将写好的程序输入到VC++中并调试运行得出方程的解；4）比较几种方法之间的联系与区别。2.1追赶法根据以上的实验原理，在VC++编辑框中输入源程序：22由原理可知：b[1]=2,c[1]=-1,f[1]=1a[2]=-1,b[2]=2,c[2]=-1,f[2]=0a[3]=-1,b[3]=2,c[3]=-1,f[3]=0a[4]=-1,b[4]=2,c[4]=-1,f[4]=0a[

3、5]=-1,b[5]=2,f[5]=0将上述系数逐个输入运行框，并经过多次调试运行，最终运行出结果如下：222.1Gauss消元法1）数据输入main(void){floatA[4][4]={{0.4096,0.1234,0.3678,0.2943},{0.2246,0.3872,0.4015,0.1129},{0.3645,0.1920,0.3781,0.0643},{0.1784,0.4002,0.2786,0.3927}};floatb[4]={0.4043,0.1550,0.4240,-0.2557};floatx[4]={0};floatAik,

4、S;inti,j,k;intsize=4;printf("A[][]");for(i=0;i

5、or(j=k;j

6、e-1];for(k=size-2;k>=0;k--){S=b[k];for(j=k+1;j

7、15,0.1129},{0.3645,0.1920,0.3781,0.0643},{0.1784,0.4002,0.2786,0.3927}};floatb[4]={0.4043,0.1550,0.4240,-0.2557};floatx[4]={0};floatAik,S,temp;inti,j,k;floatmax;//列主元的绝对值22intcol;//列主元所在的行intsize=4;printf("A[][]");for(i=0;i

8、f("");}printf("b[]");for(i=0;