二分法matlab_帮你学Matlab二分法方程求根

今天的Matlab帮你学专栏

带你回顾计算方法课上

经典的二分法方程求根问题

函假设数  在区间  上连续，并且  ，此时就可以用二分法求解。

• 刷新伪代码：

1. a1 = a;

2. b1 = b;

3. 计算中点 

4. 如果  ，那么方程的解为  ，终止

5. 如果 

• • • 如果  ， 

    • 如果  ， 

重复上述步骤2到4，直到满足误差，停止迭代。

收敛性

 ，  是第Ñ次迭代的结果，对真为实解。

• 体现

方程  在区间【0,2】上的近似解。

 ;

• 所以区间缩放到【1,2】，

• 利用matlab计算如下：

% -------------- inputs -------------------f = @(x) 3*x^2-x-3;a = 0;b = 2;% tolerance / max iterTOL = 1e-4; NI = 50;% -------------------------------------------------------% STEP 1: initializationi = 1;fa = f(a);converge = false; % convergence flag% STEP 2: iterationwhile i<=NI% STEP 3: compute p at the i's stepp = a+(b-a)/2;fp = f(p);% STEP 4: check if meets the stopping criteriaif (abs(fp)converge = true; % bisection method converged!break; % exit out of while loopelse% STEP 5i = i+1;% STEP 6if fa*fp > 0a = p; fa = fp;elseb = p;endendend

转自matlab入门到精通微信公众号