非线性方程求根数值试验非线性方程求根数值试验
试验一：二分法求根
function [x k]=bisect1(fun1,a,b,tol) if nargin tol
t=a+(b-a)/2 ;
ft=feval(fun1,t) ;
if sign(ft)==sign(fb)
fb=ft;
b=t;
elseif sign(ft)==sign(fa) fa=ft;
a=t;
end
k=k+1;
end
x=a+(b-a)/2;
>> fun1=inline('x^3-x-1');
[x k]=bisect1(fun1,1,2,0....
非线性方程求根数值试验
试验一：二分法求根
function [x k]=bisect1(fun1,a,b,tol) if nargin <4 tol=1e-5; end
fa=feval(fun1,a);
fb=feval(fun1,b);
k=1;
while abs(b-a)/2>tol
t=a+(b-a)/2 ;
ft=feval(fun1,t) ;
if sign(ft)==sign(fb)
fb=ft;
b=t;
elseif sign(ft)==sign(fa) fa=ft;
a=t;
end
k=k+1;
end
x=a+(b-a)/2;
>> fun1=inline('x^3-x-1');
[x k]=bisect1(fun1,1,2,0.005) x =
1.32421875
k =8
>> [x k]=bisect1(fun1,1,2,0.00005) x =
1.32473754882813
k =15
本文档为【非线性方程求根数值试验】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483