用弦截法求方程的根matlab,matlab 语言 用弦截法任意实数方程求实根

满意答案

dwgg2n0das8

2013.04.26

采纳率：45%    等级：12

已帮助：8202人

弦截法这个方法一般用作学习，实际用的很少。这里我提供一个较完整的弦截法求根的函数。function root=Secant(f,a,b,eps)

if(nargin==3)

eps=1.0e-6;

endf1=subs(sym(f),findsym(sym(f)),a);

f2=subs(sym(f),findsym(sym(f)),b);

if(f1==0)

root=a;

end

if(f2==0)

root=b;

endif(f1*f2>0)

disp('两端点函数值乘积大于0!');

return;

else

tol=1;

fa=subs(sym(f),findsym(sym(f)),a);

fb=subs(sym(f),findsym(sym(f)),b);

root=a-(b-a)*fa/(fb-fa);

while(tol>eps)

r1=root;

fx=subs(sym(f),findsym(sym(f)),r1);

s=fx*fa;

if(s==0)

root=r1;

else

if(s>0)

root=b-(r1-b)*fb/(fx-fb);

else

root=a-(r1-a)*fa/(fx-fa);

end

end

tol=abs(root-r1);

end

end

对于本题，输入：>> y=secant('0.15/y-y*sin(0.15^y)-1',0.1,5)结果：y = 0.1369

00分享举报