| 编程中国 | 业界新闻 | 技术文章 | 视频教程 | 下载频道 | 程序源码 | 个人空间 | 编程论坛
全能ASP/PHP/ASP.NET主机,支持月付专业 MSSQL 数据库空间,支持月付专业 MySQL 数据库空间,支持月付学习型 ASP/PHP/ASP.NET 主机 30元/年
高端软件开发 = 年薪十万不是梦赛孚耐:软件保护加密专家身份认证令牌USB KEY 
共有 116 人关注过本帖
标题:[求助]关于求解的算法
收藏  订阅  推荐  打印 
HaPpY随心
Rank: 2
等级:注册会员
帖子:119
积分:1298
注册:2007-9-28
[求助]关于求解的算法

求多项式的解时
有利用牛顿迭代法求方程 的根
利用二分法求方程 的根。
例如2X3-4X2+3X-6=0
用两种方法区别在哪里
搜索更多相关主题的帖子: 牛顿  算法  求解  方程  
2007-10-3 15:57
HaPpY随心
Rank: 2
等级:注册会员
帖子:119
积分:1298
注册:2007-9-28

#include <stdio.h>
#include <math.h>
float f(float x)
{
float y;
y=((2*x-4.0)*x+3)*x-6.0;
return y;
}

float xpiont(float x1,float x2)
{
float y;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
return y;
}

float root (float x1,float x2)
{
float x,y,y1;
y1=f(x1);
do
{
x=xpiont(x1,x2);
y=f(x);
if(y*y1>0)
{
y1=y;
x1=x;
}
else
x2=x;
}while(fabs(y)>=0.0001);
return x;
}

void main()
{
float x1,x2,f1=1.0,f2=1.0,x;
while(f1*f2>=0)
{
printf("input x1,x2:\n");
fflush(stdin);
scanf("%f,%f",&x1,&x2);
f1=f(x1);
f2=f(x2);
}
x=root(x1,x2);
printf("A root of equation is %8.4f\n",x);
}

2007-10-15 12:34
HaPpY随心
Rank: 2
等级:注册会员
帖子:119
积分:1298
注册:2007-9-28

上篇为2分法
这篇为牛顿迭代法

#include <cmath>
#include <iostream>
using namespace std;
void main()
{
double x=1,x0,f,f1;
do
{
x0=x;
f=((2.0*x0-3.0)*x0+4.0)*x0-6.0;
f1=(3*2.0*x0+2*(-3.0))*x0+4.0;
x=x0-f/f1;
}while(fabs(x-x0)>=1e-5);
cout <<"x="<<x<<endl;
}
2007-11-5 16:18
icelake
Rank: 1
等级:新手上路
帖子:16
积分:264
注册:2007-1-17

高手

2007-11-5 16:45
关于我们 | 广告合作 | 编程中国 | 清除Cookies | Archiver | WAP | TOP

编程中国 版权所有,并保留所有权利。鲁ICP备08000592号
Powered by Discuz, Processed in 0.073033 second(s), 9 queries.
Copyright©2004-2008, BCCN.NET, All Rights Reserved