循环-08. 二分法求多项式单根(20)
2015-06-09 23:42
274 查看
二分法求函数根的原理为:如果连续函数f(x)在区间[a, b]的两个端点取值异号,即f(a)f(b)<0,则它在这个区间内至少存在1个根r,即f(r)=0。
二分法的步骤为:
检查区间长度,如果小于给定阈值,则停止,输出区间中点(a+b)/2;否则
如果f(a)f(b)<0,则计算中点的值f((a+b)/2);
如果f((a+b)/2)正好为0,则(a+b)/2就是要求的根;否则
如果f((a+b)/2)与f(a)同号,则说明根在区间[(a+b)/2, b],令a=(a+b)/2,重复循环;
如果f((a+b)/2)与f(b)同号,则说明根在区间[a, (a+b)/2],令b=(a+b)/2,重复循环;
本题目要求编写程序,计算给定3阶多项式f(x)=a3x3+a2x2+a1x+a0在给定区间[a, b]内的根。
输入格式:
输入在第1行中顺序给出多项式的4个系数a3、a2、a1、a0,在第2行中顺序给出区间端点a和b。题目保证多项式在给定区间内存在唯一单根。
输出格式:
在一行中输出该多项式在该区间内的根,精确到小数点后2位。
输入样例:
3 -1 -3 1
-0.5 0.5
输出样例:
0.33
二分法的步骤为:
检查区间长度,如果小于给定阈值,则停止,输出区间中点(a+b)/2;否则
如果f(a)f(b)<0,则计算中点的值f((a+b)/2);
如果f((a+b)/2)正好为0,则(a+b)/2就是要求的根;否则
如果f((a+b)/2)与f(a)同号,则说明根在区间[(a+b)/2, b],令a=(a+b)/2,重复循环;
如果f((a+b)/2)与f(b)同号,则说明根在区间[a, (a+b)/2],令b=(a+b)/2,重复循环;
本题目要求编写程序,计算给定3阶多项式f(x)=a3x3+a2x2+a1x+a0在给定区间[a, b]内的根。
输入格式:
输入在第1行中顺序给出多项式的4个系数a3、a2、a1、a0,在第2行中顺序给出区间端点a和b。题目保证多项式在给定区间内存在唯一单根。
输出格式:
在一行中输出该多项式在该区间内的根,精确到小数点后2位。
输入样例:
3 -1 -3 1
-0.5 0.5
输出样例:
0.33
#include <stdio.h> #define CriticalValue 1e-5 double fun(double a3, double a2, double a1, double a0, double x) { return x * (x * (a3 * x + a2) + a1) + a0; } int main(int argc, char const *argv[]) { double a3, a2, a1, a0; double a, b, fa, fb, mid; scanf("%lf%lf%lf%lf", &a3, &a2, &a1, &a0); scanf("%lf%lf", &a, &b); while(b - a> CriticalValue) { mid = (a + b) / 2; if (fun(a3, a2, a1, a0, mid) == 0) { printf("%.2f\n", mid); break; } else if(fun(a3, a2, a1, a0, mid) * fa> 0) { a = (a + b) / 2; fa = fun(a3, a2, a1, a0, a); } else if( fun(a3, a2, a1, a0, mid)* fb> 0) { b = (a + b) / 2; fb = fun(a3, a2, a1, a0, b); } } if(fun(a3, a2, a1, a0, (a + b) / 2) != 0) printf("%.2f\n", (a + b) / 2); return 0; }
相关文章推荐
- 2015060908 - 为什么中国的程序员总被称为码农?
- 树莓派学习笔记之登录树莓派
- Asp的两个内置对象request和response
- Android系统启动流程
- PHP环境安装
- 黑马程序员--Java基础Day07
- maven 学习
- shell学习日志
- zookeeper节点数与watch的性能测试
- oracle 11g 忘记了sys,system,scott密码
- UVA - 455 Periodic Strings
- 循环-07. 爬动的蠕虫(15)
- ZOOKEEPER解惑
- 循环-06. 统计一行文本的单词个数(15)
- 2015060907 - 为什么有的工程师早晨醒来头发成为沟壑,角色无所谓就去上班呢?
- 中断处理
- [R语言绘图]条状图barplot
- pointer指向链表的某结点(非尾结点), 且不知道该结点的前驱结点(也不知道更前驱的结点), 请给出删除该结点的方法
- Word Break II
- 初识docker(一)