太原理工大学数值计算实验报告.doc

《太原理工大学数值计算实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验由本人当时亲自完成仅供参考，希望可以帮助大家本科实验报告课程名称：计算机数值方法B实验项目：Gauss消元法，迭代法，方程求根，拉格朗日插值法，幂法，最小二乘法，追赶法实验地点：中区ZSA301专业班级：0000学号：学生姓名：氣宇軒昂指导教师：2012年13月13日实验一：Gauss消元法实验一、实验目的《数值计算方法》是一门实践性和实用性都很强的课程，学习的目的在于应用已有的理论基础和技术方法，实现对实际计算问题的计算机求解。本实验配合《计算数值方法》课堂教学，通过上机实践实验环节，巩固课堂和书本上所学的知识、加强实践

2、能力、提高解决实际计算问题的水平、启发创新思想。其基本目的是：(1)培养理论联系实际的思路，训练综合运用已经学过的理论和实际知识去分析和解决实际问题的能力。(2)帮助学生全面消化已学的相关课程内容，深刻理解计算数值方法课程的内涵，培养使用电子计算机进行科学计算和解决问题的能力。(3)进行基本技能训练和巩固。使学生得到选择算法、编写程序、分析数值结果、写数值试验报告、课堂讨论等环节的综合训练。二、实验内容合理利用Gauss消元法、LU分解法或追赶法求解下列方程组：①②③④（n=5,10,100,…）三、主要仪器设备1.笔记本电脑

3、；2.软件vc++6.0;四、操作方法与实验步骤(代码)1.Guass消元法；#include"Stdio.h"#defineN3voidGauss_eliminant(doublea[N][N+1])/*建立系数及常数项的二维数组*/{doubleb[N];/*建立解的一位数组*/doubles;/*计算解系时的变量*/inti,j,k;for(k=0;k=0;j--){a[i+1][j]=a[i+1][j]-a[k][j]*a[i+1

4、][k]/a[k][k];}}}printf("消元后得到上三角矩阵:");printf("");for(i=0;i=0;i--)/*回代过程*/{s=0;for(j=N-1;j>i;j--){s=a[i][j]*b[j]+s;}b[i]=(a[i][N]-s)/a[i][i

5、];}for(i=0;i

6、行如下：2.输入矩阵后执行如下：再给出N=4，解决第三个方程；1.第一步执行如下：2，输入矩阵执行如下：2.Gauss_列主元消元法；#include"Stdio.h"#defineN2voidGauss_eliminant(doublea[N][N+1])/*建立系数及常数项的二维数组*/{doubleb[N];/*建立解的一位数组*/doubles;/*计算解系时的变量*/inti,j,k,m,n;for(k=0;k0;m--){for(n=N-m;n

7、++){if(a[n][k]=0;j--){a[i+1][j]=a[i+1][j]-a[k][j]*a[i+1][k]/a[k][k];}}}printf("消元后得到上三角矩阵:");printf("");for(i=0;i

8、ntf("");}printf("回代求得原方程组的解为:");printf("");for(i=N-1;i>=0;i--)/*回代过程*/{s=0;for(j=N-1;j>i;j--){s=a[i][j]*b[j]+s;}b[i]=(a[i][N]-s)/a[i][