数值计算方法实验报告(例).docVIP

云南大学数学与统计学实验教学中心 实验报告 课程：—2016学年第一学期 成绩： 指导教师： 李耀堂 学生姓名：名称： ： 专业： 年级： 一、目的 二、内容三、环境四.实验方法 五、实验过程 1实验步骤 2 关键代码及其解释 3 调试过程 六、总结 七、程序源代码： 八、教师评语 云南大学数学与统计学实验教学中心 实验报告 课程：—2016学年第一学期 成绩： 指导教师： 李耀堂 学生姓名：名称： ： 专业： 年级： 一.试验目的：练习用数值方法求解给定非线性方程。 二.实验内容：求解人口方程: 要求误差小于。 三.实验：PC计算机，FORTRAN、C、C++、VB任选一种。 四.实验方法：牛顿法 牛顿法简述：牛顿法是一种特殊的迭代法，其迭代公式为： ， 当数列收敛时，其极限值即为方程的解。 定理：给定方程 1）设； 2）在上不变号，且； 3）选取，满足； 则牛顿法产生的序列收敛于在内的唯一解。 C语言编出牛顿法的源程序。 2. 开机, 打开C语言编译程序，键入所编程序源代码. 3. 调试程序, 修改错误至能正确运行. 4. 运行程序并输出计算结果. 次数 初值 eps N xx f(xx) 1 0.45 0.000001 10 0.152886 7.120609 2 0.153 0.000001 10 0.102254 0.168226 3 0.1023 0.000001 10 0.100999 0.000107 六．总结:（称为下山条件），则有可能收敛。把新的近似值看作初值的话会比原来的取得好，有可能落入局部收敛的邻域。 （2）牛顿法要求在附近不为零。亦即只能是单根, 不能求重根。可用重根加速收敛法求重根。 （3）牛顿法的每一步迭代中，都要计算一次导数值，若计算比计算函数的近似值要麻烦的多。为了避免求导数，可用差商近似代替微商 此时牛顿迭代法改为 . (4) 由于人口方程来源于实际问题, 代表人口增长率, 其真实值不会太大, 初值不应取得过大.否则会得到该方程的另外一个解 七、程序源代码： #includemath.h #define ep 1e-4 float f (float x) { float y; y=100*exp(x)+43.5*(exp(x)-1)/x-156.4; return(y); } float df (float x) { float y; y=100*exp(x)+43.5*( x*exp(x)-exp(x)+1)/(x*x); return(y); } float root(float x) { float y; if (fabs)fep) y=x; else y=root(x-f(x)/df(x)); while(nN) return(y); } main() {float y,x0; printf(“enter x0:”); scanf(“%f”,x0); y=root(x); printf(“%f%f’,y,f(y)); 八、教师评语 云南大学数学与统计学实验教学中心 实验报告 课程：—2016学年第一学期 成绩： 指导教师： 李耀堂 学生姓名：名称： 实验编号： 专业： 年级： 一.试验目的：。二.实验内容：. 三.环境：PC计算机， C。四.实验方法： ，改用行标准化，把位置上的元素先变为1。（2）高斯消元法只是把这一列中下面行的元素消为0，而高斯—约当消去法则是把这一列元素除以外全部消为0。（3）高斯—约当消去法进行n次消元，把第n列也消为只剩一个元素为1，其余均为0。因此，对于消元也是必要条件。 高斯—约当消去法算法为： 五.实验过程: 1．编程: 用C语言编出全主元高斯—约当消去法的源程序。 2. 开机, 打开C语言编译程序，键入所编程序源代码. 3. 调试程序, 修改错误至能正确运行. 4. 运行程序并输出计算结果. 计算结果：MAT A-1 IS: . 六. 实验总结： 高斯——约当消去法算法具有以下优点和缺点： 算法优点：不用换行，换列，不用回代，精度高。 算法缺点：循环语句比较难组织，已选过主元素所在行所在列的元素不能再被选作主元素，解向量的分量也不一定按次序排列。 七、程序源代码： #include “stdlib.h” #include “math.h” #include “stdio.h” int rinv(n,a) int n; double a[]; {int*is,*js,i