太原理工大学软件工程数值实验报告

《太原理工大学软件工程数值实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、本科实验报告课程名称：计算机数值方法实验地点：行勉楼A205专业班级：软件6666班学号：2014006666学生姓名：666666指导教师：66662016年3月实验一：方程求根一、实验目的和要求在许多的实际问题屮，常常会遇到一些方程求根的问题，对于非线性方程求根的研宄是很重耍的，因此我们有必耍学会一些求根方法。实验要求是我们熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。二、实验内容和原理熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在［1,2］内的一个实根，且要求满足精度

2、

3、Axn

4、

5、x;if(f(x)*f(a)==0)break;x=(a+b)/2;printf(nx(%i)=%lf,i,x);}while(fabs(x-xi)>0.5*0.00001);getch();2.牛顿法；#include"Stdio.h”#include"math.h”doublef(doublex){/*定义函数*/returnx*x*x+4*x*x-10;}doublefl(doublex){/*由函数定义出对应的导数*/main(){doublex=1.5,xi;inti=0;do{return3*x*x+8*x;/*定义初值和变量xi*//*计算X的次数*//*计算x*/

6、i++;xi=x;x=x-f(x)/fl(x);printf(’’x(%d)=%lf’’，i,x);}while(fabs(x-xi)>0.5*0.00001);getch();五、实验结果与分析(1)二分法输出结果为:(2)牛顿法输出结果为:六、讨论与心得：从两种方法可以看出，牛顿法比二分法收敛速度快，但用牛顿法得求函数的导数，计算机无法求得，得人为输进。通过编写程序，是我们对C语言也能更进一步的了解和熟练。实验二：线性方程组的直接解法实验一：Gauss消元法实验、实验目的和要求在自然科学和工程技师屮很多问题需要归结为解线性方程组的问题，如何解线性方程组就具有实际意义。我们需耍掌

7、握Gauss消元法，并且用该方法求解方程组。二、实验内容和原理合理利用Gauss消元法、LU分解法或追赶法求解下列方程组:"123_V14'012x2—8241x313三、主要仪器设备笔记本电脑；软件vc++6.0;操作方法与实验步骤（代码）Guass消元法；#include"Stdio.h”#defineN3voidGauss_eliminant(doublea[N][N+l]){/*建立系数及常数项的二维数组*/doubleb[N];/*建立解的一位数组*/doubles;/*计算解系时的变量*/intij,k;for(k=0;k

8、N-l;i++){fora=N;j>=0;j-){a[i+l][j]=a[i+l][j]-a[k]U]*a[i+l][k]/a[k][k];}}}printfC•消元后得到上三角矩阵:”)；printf("n);for(i=0;i=0;i-){/*回代过程*/s=0;fora=N-l;j>i;j-){s=a[i]Lj]*b[j]+s;

9、}b[i]=(a[i][N]-s)/a[i][i];for(i=0;i