数值方法方程求根实验报告.docxVIP

PAGE 3 PAGE 一、实验概述： 1 实验目的 通过采用迭代法的程序设计，系统的理解和掌握c语言的编程技巧。学会利用c语言进行科学计算，学会查资料和工具书。 2 实验原理 利用构造数列，利用数列的收敛性逐步逼近方程的根来求解 3实验环境 Visual C++6.0 WIN10操作系统 二、实验内容： 1实验方案 1.1用简单迭代法求ex+10x-2=0的方程的根，当满足|Xn+1-Xn|=10-6时结束迭代，并说明迭代收敛的理由。 构造等价方程x=0.2-0.1ex,则φ(x)=0.2-0.1ex。 给定初值x0=0,构造迭代公式Xn+1=0.2-0.1exn。 在VC++6.0中编写程序如下并运行： #includestdio.h #includemath.h #define ESP 1e-6 void main() { double x0,x1,b; printf(INPUT x0:); scanf(%lf,x0); do { x1=0.2-0.1*exp(x0); b=x0; x0=x1; } while(fabs(x0-b)=ESP); printf(%f\n,x0); } 1.2 用牛顿法求方程x=tgx的最小正根和在x=100附近的根，当满足|Xn+1-Xn|=10-6时结束迭代。并讨论初值的变化对收敛的影响。 （1）对原方程变形为x=arctanX+π,令f(x)=x-arctanx-π,f’(x)= （2）选定初值x0=4.5,构造牛顿迭代公式 （3）在VC++6.0中编写程序如下并运行： #includestdio.h #includemath.h #define ESP 1e-6 #define PI 3.1415926 void main() { double x0,x1,f,f1,a,b; printf(INPUT x0:); scanf(%lf,x0); a=floor(x0/PI); do { f=x0-atan(x0)-a*PI; f1=x0*x0/(1+x0*x0); x1=x0-f/f1; b=x0; x0=x1; } while(fabs(x0-b)=ESP); printf(%f\n,x0); } 2实验步骤 1.1 a)程序编译出错： b)找到出错位置： c)修改为“x1=0.2-0.1*exp(x0)”后调试运行为： 调试成功，运行程序得出结果。 1.2 a)编译运行出错： b)找到出错位置 c)修改为“f=x0-atan(x0)-a*PI;”，调试运行。运行得出结果。 3实验结论（数据及分析结果） 对于方程ex+10x-2=0，输入初值为0时结果如下： 因为x∈(0.1),φ’(x)=-0.1ex1,所以该迭代格式收敛。 对于方程x=tgx，输入初值4.5时运行结果如下: 当输入x为100时： 牛顿迭代是局部收敛的，故迭代在方程的根的附近是收敛的，所以初值的选择对牛顿迭代的收敛性有影响，若初值选在根的附近则迭代收敛，若初值选择离根远则发散。 4实验小结（收获体会） 通过这次实验，基本掌握了利用C语言解决数值计算中的方程求根问题，从最初的编写算法到调试再到得出结果，虽然有困难但是通过翻阅资料，查工具书等，最终顺利完成了任务，同时也加深了对于迭代法的认识。 三、指导教师评语及成绩 评 语 评语等级 优 良 中 及格 不及格 1.实验报告按时完成,字迹清楚,文字叙述流畅,逻辑性强 2.实验方案设计合理 3.实验过程（实验步骤详细,记录完整,数据合理,分析透彻） 4实验结论正确. 成 绩： 指导教师签名： 批阅日期： 附录1：源 程 序 1.1 #includestdio.h #includemath.h #define ESP 1e-6 void main() { double x0,x1,b; printf(INPUT x0:); scanf(%lf,x0); do { x1=0.2-0.1*exp(x0); b=x0; x0=x1; } while(fabs(x0-b)=ESP); printf(%f\n,x0); } 1.2 #includestdio.h #includemath.h #define ESP 1e-6 #define PI 3.1415926 void main() { double x0,x1,f,f1,a,b; printf(INPUT x0:); scanf(%lf,x0); a=floor(x0/PI); do { f=x0-atan(x0)-a*PI;