太原理工大学计算机数值方法实验报告.docxVIP

本科实验报告 课程名称： 计算机数值方法 实验地点： 专业班级：_学号： 学生姓名： 指导教师： 成 绩： 年 月曰 太原理工大学学生实验报告 学院名称 计算机科学与技术学院 专业班级 学号 学生姓名 实验日期 成绩 课程名称 计算机数值方法 实验题目 实验一方程求根 、实验目的和要求: 实验目的： 熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程 进彳亍根的求解。求方程：f(x)二x*x*x+4*x*x?10=0在［1,2］内的一个实根，且要求满足精度 |x -xn|0.001 实验要求： 1.应用C,C++或JAVA编出通用程序，源程序要有详细的注释和说明； 2?比较计算结果，对不同方法进行比较分析； 3.实验完成，提交实验结果并写出报告，分析计算结果是否符合问题的耍求，找出计 算成功的原因或计算失败的教训。 二、 实验内容和原理： 增值寻根法：基本思想为：从初始值xO开始，按规定的一个初始步长h來增 值。令x (n+1) =x (n) +h,(n=0,l,2...),同时计算f (x (n+1)).在增值过程屮会遇到三种 情况:1. f (x(n+l)) =0,此时x(n+l)即为方程根。 f (x(n))和f (x (n+1))同号，说明区间内无根。 f (x(n))和f (x (n+1))同号，说明区间内冇根，则把步长缩小，直至满足精度要求为 止，x(n)或x(n+l)就是满足精度的近似根。 二分法：基本思想为：设f(x)在［a,b］内连续，.ft f(a)*f(b)0,则方程f(x)=O在(a,b) 内有实根x*.然后逐步对分区间［a,b］,通过判断两端点函数值乘积的符号，进一步缩小有 根区间从而求出满足精度要求的近似值。 牛顿迭代法：基本思想为给定一个初始值由牛顿法的迭代公式： x(n+1 )=x(n)-f(x(n))/f9(x(n)) (n=0,1,2...)逐步求出 x (n),直至(x (n+1) -x(n))的绝对值小 于给定精度，则x (n+1)即可作为近似值。 双点割线法：由给出的两个初始近似值，再根据双点割线法迭代公式： x(n+1 )=x(n)-(f(x(n))/(f(x(n))-f(x(n-1))) *(x(n)-x(n-1 )),(n= 1,2,3...)逐 步求出 x ( n ),直至 x(n+l)-x(n)的绝对值满足精度，则x (n+1)即口J作为近似值。 单点割线法：由给出的两个初始近似值，再根据双点割线法迭代公式： x(n+l)=x(n)-(f(x(n))/(f(x(n))-f(x(0)))*(x(n)-x(0)),(n=l,2,3...)逐步求岀 x (n),直至 x(n+l)-x(n) 的绝对值满足精度，则x (n+1)即可作为近似值。 三、 主要仪器设备：笔记本电脑 四、操作方法： 源代码： 增值寻根法： # includestdio.h double fun(double x){ 〃原函数 retum(x*x*x+4*x*xJ0);//求解方程 f(x)=x*x*x+4*x*x?10=0 的根，精度为 10-3. } int main(){ double a=1.25,h=l,x=a; printff初始近似值为:%lf\n”,a); do{ if(fun(x)==O){printf(H根为：%f\x);return 0;} /*如果初始值函数值为0,则初始 值即为根勺 else if(fun(x)*fun(x-t-h)0) /*如果 fun(x)和 fun(x+h)同号则使 X 加 h 并跳岀本 次循环执行下一次勺 {x=x+h;continue;} else if(fun(x)*fun(x+h)vO) 〃若异号则说明在X和X+h Z间存在函数根，则缩 〃短步长继续寻根 {h=h/10.0;} }while(h0.001);//当不满足精度要求吋继续执行循环体 printfC根为：％i\n”,x);〃跳出循环说明满足精度耍求则x可近似作为方程根 return 0; } 二分法： #includestdio.h # includemath.h #define esp le-3 〃精度 double f(double x) 〃原函数 { return (x*x*x+4*x*x?10); } double root(double (*fun)(double),double left,double right,double deviation)//用二分法求 方程根 { 〃其中形参*fun为指向原函数的指针 double x,y; while(fabs(right-left)deviation)//当不满足精度要求继续执行循环体 { x=(right+left)