二分法的MATLAB代码实现

pulbieup

二分法是计算机求解方程的常用算法，很容易在matlab中实现。

%计算f（x）= 0的近似解，容差为tol
, |! ]  J. t9 o' }
% L, }' o4 t8 s/ }% Ffuction xc = bisec(f,a,b,tol)
8 E( r" i3 t: ?- }' vif sign(f(a))*sign(f(b)) >= 0
) A; q+ f5 P% g9 n    error('f(a)f(b)<0 not satisfied!')
end
& F! `, m% i# Z$ Ufa = f(a);
3 |3 O9 |; V6 I& o5 o. rfb = f(b);
; o% c  W6 v' m) G- N) cwhile (b-a)/2 > tol
2 c9 j5 o. d$ k   c = (a+b)/2;
   fc = f(c);
8 x' b. z& P* X/ Q   if fc == 0
- n: q3 D6 O/ B/ M     break
   end
   if sign(fa)*sign(fb)<0
1 O; t. s; T+ U7 D: `8 ~. L) K       b=c;fb = fc;
   else
4 s5 `/ [+ I; i       a = c;fa = fc;
   end
5 q0 u+ r6 @3 rend
3 O9 i$ {& Z6 n# j! txc = (a+b)/2

; G6 e( |, S5 O' M5 B: T5 |

7 I' j8 C& F6 z! r  n在matlab中首先要在命令行中键入：
" k9 N: U2 U4 a+ W* A& H
( m2 Z, o0 U4 t2 J# [6 j f = @(x) x^3+x-1
: n# ], m. n  ~  B& H/ a: o' l6 T# @: P
$ j2 \0 U  U6 y2 o

yin123

厉害的人儿