- 1、本文档共35页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
1.编写一个函数,用牛顿迭代法求方程f(x)=ax^3+bx^2+cx+d=0(^为乘方号)在1附近
的根。a、b、c、d分别为1、2、3、4。 0
提示:牛顿迭代公式为:
x1=x0-f(x0)/f(x0)=x0-(a*x0*x0*x0+b*x0*x0+c*x0+d)/(3*a*x0*x0+2*b*x0+c)
#includeiostream
#includecmath
using namespace std;
double niudun(double a,double b,double c,double d){
double x1,x0=1.0;
x1=x0-(a*x0*x0+b*x0*x0+c*x0+d)/(3*a*x0*x0+2*b*x0+c);
while(fabs(x0-x1)1e-6)
{
x0=x1;
x1=x0-(a*x0*x0+b*x0*x0+c*x0+d)/(3*a*x0*x0+2*b*x0+c);
}
return x1;
}
int main ()
{
coutax^3+bx^2+cx+d=0(^为乘方号)在一附近的根为: \n;
coutniudun(1,2,3,4)endl;
system(pause);
return 0;
}
2.编写程序,用普通迭代法求方程f(x)=x+lnx-1.7=0的近似实根, 0
迭代初值为1.5,精确到0.0001。
提示:必须把原方程同解变形为迭代形式:x=1.7-lnx
其中:lnx表示x的自然??数(头文件math.h中有对应的函数log(x))
#includeiostream
#includeiomanip
#includecmath
using namespace std;
int main()
{
double x1,x0=1.5;
x1=1.7-log(x0);
while(fabs(x0-x1)1e-4)
{
x0=x1;
x1=1.7-log(x0);
}
coutx+lnx-1.7=0的近似根为: \t;
coutsetprecision(5)x1endl;
system(pause);
return 0;
}
3.求方程5x^3+2x^2-12x-16=0在区间[1,2]上的近似根。
#includeiostream
#includeiomanip
#includecmath
using namespace std;
int main()
{
double x1,x0=1.0;
x1=x0-(5*pow(x0,3.0)+2*pow(x0,2.0)-12*x0-16)/(15*pow(x0,2.0)+4*x0-12);
while(fabs(x0-x1)1e-4)
{
x0=x1;
x1=x0-(5*pow(x0,3.0)+2*pow(x0,2.0)-12*x0-16)/(15*pow(x0,2.0)+4*x0-12);
}
cout5x^3+2x^2-12x-16=0在区间[1,2]上的近似跟为:;
coutsetprecision(5)x1endl;
system(pause);
return 0;
}
4.用迭代法求150的算术平方根,迭代公式为:x1=(x0+150/x0)/2,迭代初值x0=1,
相邻两次求出的x的差的绝对值小于10的-4次方,停止迭代。
#includeiostream
#includecmath
using namespace std;
int main()
{
double x1,x=1.0;
x1=1/2.0*(x+150.0/x);
while(fabs(x-x1)1e-4)
{
x=x1;
x1=1/2.0*(x+150.0/x);
}
cout150的算术平方根为: \n;
coutxendl;
system(pause);
return 0;
}
5.编写程序,用普通迭代法求方程f(x)=x-e^x+2.5=0的近似实根,
迭代初值为-2.5,精确到0.0001。
提示:必须把原方程化成其等价的迭代形式:x=e^x-2.5
其中:e^x表示e的x次幂(头文件math.h中有对应的函数exp(x))
#inc
文档评论(0)