工程常用算法.doc

哈尔滨工程大学机电工程学院 综合实践大作业 学年学期： 2015-2016学年春季学期 课程代号： 0907322[1] 课程名称： 工程常用算法 授课教师： 郭 健 副教授 小组成员 班级 学号 姓名 成绩 座机电话号码 座机电话号码05 郭俊杰 座机电话号码 座机电话号码10 凌泽宇 座机电话号码 座机电话号码27 张明旺 2016年 6 月 15 日 《工程常用算法》综合实践作业一 作业名称：非线性方程求根 完成日期：2015年3月 26日 组号：8 班级 学号 姓名 主要工作说明 自评成绩 座机电话号码 座机电话号码05 郭俊杰 程序流程图 A 座机电话号码 座机电话号码10 凌泽宇 C语言算法编程 座机电话号码 座机电话号码27 张明旺 总结排版 知方程，请分别用二分法、牛顿法和割线法求此方程的根。误差限取：。注意先确定出方程的有根区间。 程序流程图 二分法 牛顿迭代法 割线法 三、完整的程序及简要的注释 1.牛顿法 #include #include double f double x return x*x*x + 2*x*x - exp x ; double fd double x return 3*x*x + 4*x - exp x ; void NewtonRoot double x0, double ex double x1, fx1, fdx0; int k 1; printf 牛顿法求根\n ; while 1 fdx0 fd x0 ; if fdx0 0 printf 重新输入初始值！ ; break; else x1 x0 - f x0 /fdx0; fx1 f x1 ; printf k %-3d x %.12f f x %.12f\n, k,x1,fx1 ; if fabs x1-x0 ex break; else x0 x1; k++; int main void // 牛顿法 NewtonRoot 0,1e-12 ; return 0; 2.二分法 #include #include double f double x return x*x*x + 2*x*x - exp x ; double fd double x return 3*x*x + 4*x - exp x ; void BiPartition double a, double b, double ex printf 二分法求根\n ; if f a *f b 0 printf 区间内无根! ; return; double x, fx; int k 1; while 1 x a+b /2; fx f x ; printf k %-3d x %.12f f x %.12f\n, k,x,fx ; if b-a ex break; else if f a *fx 0 b x; else a x; k++; int main void // 二分法计算 BiPartition -2, -1, 1e-12 ; return 0; a x; k++; int main void // 二分法计算 BiPartition -2, -1, 1e-12 ; return 0; 3.割线法 #include #include double f double x return x*x*x+2*x*x-exp x ; double fd double x return 3*x*x+4*x-exp x ; void Newton double x0,double x1,double ex printf 割线法求根\n ; double x2,fx; int k 1; while 1 x2 x1-f x1 * x1-x0 / f x1 -f x0 ; fx f x2 ; printf k %-3d x2 %.12f f x %.12f\n,k,x2,fx ; if fabs x2-x1 ex break; else x0 x1; x1 x2; k++; int main void // 割线法计算 Newton 0.1,0.5, 1e-12 ; return 0; 程序运行结果 对不同实现方法的运行结果进行比较 问题与总结 一、计算公式及计算方法