借助 Matlab 使用二分法求解方程的根

借助 Matlab 使用二分法求解方程的根第一次使用Matlab,遂将过程详细记录之。图中标注①是工作目录,即代码存放的目录;标注②是编辑器,即我们写代码的地方;标注③是命令行,是我们执行语句的地方。本次实验我们是在这里执行二分法的函数。例题:应用二分法求解方程\(x^3-x-1=0\)在区间\([1,

大家好,欢迎来到IT知识分享网。

第一次使用 Matlab,遂将过程详细记录之。

20210317185838

图中标注 ① 是工作目录,即代码存放的目录;标注 ② 是编辑器,即我们写代码的地方;标注 ③ 是命令行,是我们执行语句的地方。本次实验我们是在这里执行二分法的函数。

例题: 应用二分法求解方程 \(x^3 – x – 1 = 0\) 在区间 \([1, 1.5]\) 内的数值解 \(x_k\),要求绝对误差小于 \(10^{-8}\).

解答 如下。

代码:

half.m 脚本:

function x = half(a, b, tol)% tol 是 tolerance 的缩写,表示绝对误差
c = (a + b) / 2; k = 1;
m = 1 + round((log(b - a) - log(2 * tol)) / log(2)); % <1>
while k <= m
    if f(c) == 0
        c = c;
    return;
    elseif f(a) * f(c) < 0
        b = (a + b) / 2;
    else
        a = (a + b) / 2;
    end
    c = (a + b) / 2; k = k + 1;
end
x = c; % 这里加分号是为了不再命令行中输出
k % 不加分号就会在控制台输出
c

f.m 脚本,这是 half.m 中调用的 f() 函数。

function y = f(x)
y = x^3 - x -1;

然后我们在命令行执行:

20210317190826

可以看出,最后求解得到的 \(x = 1.3247\),即输出的 ans,迭代次数 \(k = 27\).

关于代码 half.m 中的标注 <1>,有如下解释:

20210317191251

20210317191337

注意,在 Matlab 中,log() 函数的底是 \(e\).

补充例题(感兴趣的朋友可以自行测试):

20210317191627

免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://yundeesoft.com/33758.html

(0)

相关推荐

发表回复

您的电子邮箱地址不会被公开。 必填项已用 * 标注

关注微信