原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
迭代法(二分、弦截、牛顿)
（1）牛顿迭代法
实验代码如下：
#include stdafx.h
#include iostream
#include math.h
using namespace std;//求y=x*x-3*x-1的根
float f(float x)
{return (pow(x,2)-3*x-1);
}
float f1(float x) //f1(x)是f(x)的导数
{ return (2*x-3);
}
void main()
{ float x,x1,y1,y2;
cout请输入一个初始值:X=endl;
cinx;
do{
x1=x;
y1=f(x);
y2=f1(x1);
x=x1-y1/y2;
}
while (fabs(x-x1)=1e-5);
cout方程的根是：x1=x1endl;
}
弦截迭代法
实验代码如下：
#include stdafx.h
#includestdio.h
#includemath.h
#define e 1e-5
double f(double x)
{
return(pow(x,2)-3.0*x-1.0);
}
double f1(double X1, double X2)
{
double y;
y = (double)(X1*f(X2) - X2*f(X1))/(f(X2) - f(X1));
return(y);
}
double Root(double X1, double X2)
{
double X,y,y1;
y1 = f(X1);
do
{
X = f1(X1, X2);
y = f(X);
if (y*y1 0)
{
y1 = y;
X1 = X;
}
else
{
X2 = X;
}
}
while (fabs(y) = e);
return(X);
}
void main()
{
double X1,X2,f1,f2,X;
printf(======== 求 x2-3x-1= 0 的全部实根（弦截法） =======\n\n);
do
{
printf(请输入X的最小值：X1 = );
scanf(%lf,X1);
printf(请输入X的最大值：X2 = );
scanf(%lf,X2);
f1 = f(X1);
f2 = f(X2);
}
while (f1*f2 = 0);
X = Root(X1,X2);
printf(\n该方程的近似根为:%lg\n,X);
getchar();
}
二分迭代法
实验源代码如下：
#include stdafx.h
#includeiostream
#define f(x) (x*x-3*x-1)
using namespace std;
void main()
{
float a,b,c;
cout请输入两端点a，b的值：endl;
couta=;
cina;
coutb=;
cinb;
while ((b-a)1e-5)
{
c=(a+b)/2;
if(f(c)==0) break;
else if(f(a)*f(c)0) b=c;
else a=c;
}
cout该方程的根为 x=cendl;
}
文档评论（0）