《二分法牛顿迭代法》由会员分享，可在线阅读，更多相关《二分法牛顿迭代法（10页珍藏版）》请在金锄头文库上搜索。
1、哈尔滨工业大学（威海）实验报告纸- 1 -2014 级硕士研究生数值分析上机实习级硕士研究生数值分析上机实习 (第一次第一次)姓名：乔永亮姓名：乔永亮 学号：学号：14S030125 学院：船舶与海洋工程学院学院：船舶与海洋工程学院实习题目：实习题目：分别用二分法和 Newton 迭代法求方程的02010223xxx根.实习目的：实习目的：掌握两种解法，体会两种解法的收敛速度.实习要求：实习要求：用 C 程序语言编程上机进行计算，精确到 8 位有效数字.报告内容：报告内容：1. 确定实根的个数以及所在区间.解：对函数求导，得。32( )21020f xxxx2( )34100fxxx易知恒成立
2、，所以函数没有极值，只有一个实根。又可以知道,( )0fx(x)f(1)0f方程在区间有一个实根，且为奇数重根，可以二分法和 Newton 求解(2)0f(1,2)2. 将最后两次计算结果填入下表（保留 8 位数字）： 方法kkx1k1kx二分法131.3687744141.3688354Newton 法21.368814831.36880823. 实习过程中遇到哪些问题？如何解决？有何心得体会？在编程的过程中由于对基本计算原理的理解有一定不足，同时对编程语言的不熟悉，导致在编程过程中错误百出，耗费了大量时间。但是通过课本以及网络对所需知识的不断学习，通过尝试不同的方法，最终还是得到了几种不同
3、的思路与方法。通过这次编程，深深的感受到自己的不足，同时也明白了数学与计算机编程的紧密结合，不努力提高自己在当今社会就要被淘汰。哈尔滨工业大学（威海）实验报告纸- 2 -4. 两种解法的计算程序（此页写不下时可以加页）：二分法（二分法（FortranFortran 语言）语言）programprogram Analysis1 realreal:a,b,c,m realreal:fa,fca=1.b=2.m=0.0001 !- dodo whilewhile(abs(b-a)=m)c=(a+b)/2fa=a*3+2.*a*a+10.*a-20fc=c*3+2.*c*c+10.*c-20ifif(
4、fa*fc=m)x1=x0-f(x0)/df(x0)x0=x1i=i+1writewrite(*,“(i4,f10.7)“)i,x0哈尔滨工业大学（威海）实验报告纸- 3 -endend dodo pausepause endend programprogram Analysis2 !定义原函数- functionfunction f(x) implicitimplicit nonenone realreal f,x f=x*3+2.*x*x+10.*x-20. returnreturn endend functionfunction !定义原函数的导函数- functionfunction
5、df(x) implicitimplicit nonenone realreal df,x df=3.*x*x+4.*x+10. returnreturn endend functionfunction哈尔滨工业大学（威海）实验报告纸- 4 -2014 级硕士研究生数值分析上机实习级硕士研究生数值分析上机实习 (第二次第二次)姓名：姓名： 学号：学号： 学院：学院：实习题目：实习题目：计算 8 阶三对角矩阵的行列式.)235. 0,274. 1,235. 0(tridiagA实习目的：实习目的：掌握计算行列式的方法.实习要求：实习要求：首先选择一种算法，然后用 C 程序语言编程上机进行计算.报
6、告内容：报告内容：1. 简单描述所采用的算法：2. 计算结果： A3. 实习过程中遇到哪些问题？如何解决？有何心得体会？哈尔滨工业大学（威海）实验报告纸- 5 -4. 写出 C 语言计算程序（此页写不下时可以加页）：哈尔滨工业大学（威海）实验报告纸- 6 -2014 级硕士研究生数值分析上机实习级硕士研究生数值分析上机实习 (第三次第三次)姓名：姓名： 学号：学号： 学院：学院：实习题目：实习题目：分别用 Jacobi 迭代法和 Gauss-Seidel 迭代法求解线性方程组 9 . 14 . 35 . 16 . 84 . 22 . 78 . 17 . 2 7 . 64 . 38 . 91 .
7、 2zyxzyxzyx实习目的：实习目的：感受两种迭代法的收敛速度.实习要求：实习要求：首先构造收敛的 Jacobi 迭代法和 Gauss-Seidel 迭代法，然后用 C 程序语言编程上机进行求解，初始值均取为 0，精确到 4 位小数.报告内容：报告内容：1. 写出收敛的 Jacobi 迭代法和 Gauss-Seidel 迭代法：哈尔滨工业大学（威海）实验报告纸- 7 -2. 将最后一次迭代次数 与相应的迭代结果填入下表 k迭代法k)(kx)(ky)(kzJacobi Gauss-Seidel3. 实习过程中遇到哪些问题？如何解决？有何心得体会？4. C 语言计算程序：哈尔滨工业大学（威海）
8、实验报告纸- 8 -哈尔滨工业大学（威海）实验报告纸- 9 -2014 级硕士研究生数值分析上机实习级硕士研究生数值分析上机实习 (第四次第四次)姓名：姓名： 学号：学号： 学院：学院：实习题目：实习题目：分别用复化梯形公式和复化 Simpson 公式计算积分的近似值和. dxxxIx20e2sinnImS实习目的：实习目的：体会两种复化求积公式的收敛性与收敛速度.实习要求：实习要求：用 C 程序语言编程上机进行计算，结果要有八位有效数字.报告内容：报告内容：1. 写出求和的复化求积公式：nImS2. 写出具有八位有效数字的计算结果： ；=nnI；=mmS3. 实习过程中遇到哪些问题？如何解决？有何心得体会？哈尔滨工业大学（威海）实验报告纸- 10 -4. C 语言计算程序（此页写不下时可以加页）：
