网站大量收购独家精品文档,联系QQ:2885784924

求高阶方程的根.pdf

  1. 1、本文档共11页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多

求一元高阶方程的根

一元非线性方程求根

(1)牛顿迭代法

牛顿迭代法又称牛顿切线法:先任设定一个与真实的根

接近的值x作为第一次近似根,由x求出f(x),过(x,f(x))

00000

点做f(x)的切线,交x轴于x,把它作为第二次近似根,再由

1

x求出f(x),过(x,f(x))点做f(x)的切线,交x轴于x,……

11112

如此继续下去,直到足够接近(比如|x-x|1e-6时或

10

*

|f(x)|1e-6)真正的根x为止。

1

而f(x)=f(x)/(x-x)所以x=x-f(x)/f(x)

00101000

例如,用牛顿迭代法求下列方程在1.5附近的根

32

2x-4x+3x-6=0。

法一|x-x|1e-6

10

#includemath.h

main()

{floatx0,x1,fx,f1x;

x0=1.5;

fx=2*x0*x0*x0-4*x0*x0+3*x0-6;

f1x=6*x0*x0-8*x0+3;

x1=x0-fx/f1x;

while(fabs(x1-x0)=1e-6)

{x0=x1;

fx=2*x0*x0*x0-4*x0*x0+3*x0-6;

f1x=6*x0*x0-8*x0+3;

x1=x0-fx/f1x;

}

printf(%f\n,x1);

}

法二:|f(x)|1e-6

1

#includemath.h

main()

{floatx0,x1,fx,f1x;

x0=1.5;

fx=2*x0*x0*x0-4*x0*x0+3*x0-6;

while(fabs(fx)=1e-6)

{f1x=6*x0*x0-8*x0+3;

x1=x0-fx/f1x;

x0=x1;

fx=2*x0*x0*x0-4*x0*x0+3*x0-6;

}

printf(%f\n,x1);

}

法三:

#includemath.h

main()

{floatx1,x0,f,f1;

x1=1.5;

do{x0=x1;

f=2*x0*x0*x0-4*x0*x0+3*x0-6;

f1=6*x0*x0-8*x0+3;

x1=x0-f/f1;}while(fabs(x1-x0)=1e-5);

printf(%f\n,x1);}

(2)二分法

算法要领是:先指定一个区间[x,x],如果函数f(x)在此

12

区间是单调变化的,则可以根据f(x)和f(x)是否同号来确定

12

方程f(x)=0在区间[x,x]内是否有一个实根;如果f(x)和f(x)

1212

同号,则f(x)在区间[x,x]内无实根,要重新改变x

文档评论(0)

火龙果的春天 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档