C语言用二分法求方程的近似解的方法
来源:互联网 发布:免费永久域名 编辑:程序博客网 时间:2024/05/21 08:45
C语言用二分法求方程的近似解的方法
基础知识以下地址:http://jingyan.baidu.com/article/597a06438def54312a524376.html
我的应用如下:求x^5+x^3+7=0,求取x的值是多少?
编译环境:Window8 +C Free
源程序如下:
#include "stdio.h"
#include "math.h"
#include "stdlib.h"
double func (double x)//这个函数主要实现求y=x^3+x^5+7 的值
{
double y;
y=pow(x,3)+pow(x,5)+7; //函数表达式可自己定义
return y;
}
我的应用如下:求x^5+x^3+7=0,求取x的值是多少?
编译环境:Window8 +C Free
源程序如下:
#include "stdio.h"
#include "math.h"
#include "stdlib.h"
double func (double x)//这个函数主要实现求y=x^3+x^5+7 的值
{
double y;
y=pow(x,3)+pow(x,5)+7; //函数表达式可自己定义
return y;
}
double CalRoot(double m,double n)//这里用的迭代法和二分法的结合吧!
{//这个自定义函数就是不断的取y(m)和y(n)的值,和他们的中间值y(x)的比较
//不断取中间值y=func(x),当y接近于0时,那么结束循环!此时x的值便是我们求到
//的近似值!
double x,y,y1,y2,a,b;
a=m; b=n;
while (1)
{
x=(a+b)/2;
y=func(x);
if(fabs(y)<=1e-5) //相当于|y|<0.00001,当y为浮点型数据时,用这样的方式来判断与0的大小!
{
return x;//这个if语句也是判断循环结束的标志!
}
y1=func(a);
y2=func(b);
if(y1*y<0) b=x;//如果y1和y为异号的话就将a,b中间值x传递给b
if(y2*y<0) a=x;//如果y1和y为异号的话就将a,b中间值x传递给a
}
}
{//这个自定义函数就是不断的取y(m)和y(n)的值,和他们的中间值y(x)的比较
//不断取中间值y=func(x),当y接近于0时,那么结束循环!此时x的值便是我们求到
//的近似值!
double x,y,y1,y2,a,b;
a=m; b=n;
while (1)
{
x=(a+b)/2;
y=func(x);
if(fabs(y)<=1e-5) //相当于|y|<0.00001,当y为浮点型数据时,用这样的方式来判断与0的大小!
{
return x;//这个if语句也是判断循环结束的标志!
}
y1=func(a);
y2=func(b);
if(y1*y<0) b=x;//如果y1和y为异号的话就将a,b中间值x传递给b
if(y2*y<0) a=x;//如果y1和y为异号的话就将a,b中间值x传递给a
}
}
void main ( )//原来程序中过多的部分我已经全部删除掉了,这版本比较简洁!
{
double a,b,a1,b1,x;
a=-3; b=-1; //a,b的值是搜索根的区间范围,可以适当调整,范围不可过大这里!
a1=a; b1=b;//还要注意func(a)和func(b)那么他们的值必须为异号才行!这是二分法使用的前提!
x=CalRoot(a1,b1);//也就是说范围a,b的取值要注意很多这里我取得是a=-3;b=-1
printf("X=%f\n",x);
}
0 0
- C语言用二分法求方程的近似解的方法
- 二分法求方程的近似解
- 二分法——求近似方程的解
- 分治法--二分法求方程近似解
- 二分法求方程的解
- 二分法求方程近似根
- C语言之基本算法23—二分法求方程近似根
- 二分法求解方程的近似解(sicily 1017)
- 用牛顿迭代法和二分法求方程的根【C语言】
- C代码:二分法求三次方程近似根
- 二分法求方程2x³-4x²+3x-6=0的近似解
- 用二分法求方程的根
- 二分法求方程根--C语言
- 方程求近似解方法Java实现
- 二分法求方程的根
- 二分法求方程的根
- 二分法求方程的根
- 二分法求方程的根
- TLD源码学习——cvCalcOpticalFlowPyrLK源代码
- 第二周 统计成绩
- 几篇负载均衡与热备的文章
- PAT (Advanced) 1024. Palindromic Number (25)
- 有什么好计较 一转眼 生命到尽头
- C语言用二分法求方程的近似解的方法
- 九度:1445<并查集>
- 2013年人气最高的JavaScript框架排名
- 白话经典算法系列之六 快速排序 快速搞定
- KMP算法
- Eclipse中设置编码的方式
- iOS7判断获取相册和拍照的权限
- 【marks】rails资料
- Vertex Buffer Sharing between DirectX11 and Directx9(转自GAMEDEV)