二分法求解非线性方程及MATLAB代码实现

二分法求解非线性方程及MATLAB代码实现

本文将介绍使用二分法求解非线性方程，并提供MATLAB代码实现。

问题描述

使用二分法求解方程 f(x)=(x-1)^3-3x+2=0 在区间[2,4]上的根，要求计算精度为10^-5，并绘制函数图像以验证计算结果的正确性。

实验目的

• 理解二分法的思想- 能够编程实现二分法，求解非线性方程

MATLAB代码实现

以下是使用二分法求解给定方程的MATLAB代码：matlab% 定义方程函数f = @(x) (x-1)^3 - 3*x + 2;

a = 2; % 区间左端点b = 4; % 区间右端点tolerance = 1e-5; % 计算精度

% 绘制函数图像x_vals = linspace(a, b, 1000);y_vals = arrayfun(f, x_vals);plot(x_vals, y_vals);xlabel('x');ylabel('f(x)');title('函数 f(x) 的图像');grid on;

% 二分法求解方程while (b - a) > tolerance c = (a + b) / 2; if f(a) * f(c) < 0 b = c; elseif f(c) * f(b) < 0 a = c; else break; % 找到精确解 endend

disp('方程的根为：');disp(c);

代码解释

1. 定义方程函数: 使用函数句柄 f = @(x) (x-1)^3 - 3*x + 2; 定义了需要求解的方程。2. 设置参数: 设置区间左右端点 a 和 b，以及计算精度 tolerance。3. 绘制函数图像: 使用 linspace 生成x坐标，使用 arrayfun 计算对应的y坐标，最后使用 plot 函数绘制函数图像，并添加标签和标题。4. 二分法求解: 使用 while 循环，不断迭代缩小区间，直到满足计算精度。在每次迭代中，计算中点 c 的函数值，根据函数值符号判断根的区间，并更新区间端点。5. 输出结果: 最后输出找到的方程根。

总结

本文介绍了使用二分法求解非线性方程的方法，并使用MATLAB代码实现了该算法。通过绘制函数图像，可以直观地验证计算结果的正确性。二分法是一种简单有效的求解非线性方程的方法，但需要注意其适用条件。