二分法求函数方程根的matlab实现(内附例题和代码)
二分法求函数方程根的matlab实现(内附例题和代码)摘要:本文主要给出二分法求函数方程根的matlab代码,并给出例题对代码进行验证。并且代码具有很好的交互性,按照提示操作即可。对于二分法求方程的根的思想我们并不陌生,这里便不在赘述(看下图或百度即可):重点来啦,下面给出matlab代码:%%%%% 二分法求函数的实根%%%%%%%%%%Liu Deping%%%%%%%%%2020.06.09
·
二分法求函数方程根的matlab实现(内附例题和代码)
摘要:本文主要给出二分法求函数方程根的matlab代码,并给出例题对代码进行验证。并且代码具有很好的交互性,按照提示操作即可。
对于二分法求方程的根的思想我们并不陌生,这里便不在赘述(看下图或百度即可):
重点来啦,下面给出matlab代码:
%%%% 二分法求函数的实根%%%%%
%%%%% Liu Deping %%%%%
%%%%% 2020.06.09 %%%%%
clc;
clear all
format short
s=input('请输入函数表达式:f = ','s');
f=inline(s);
a = input('请输入区间左端点:a=');
b = input('请输入区间右端点:b=');
eps=input('请输入停止精度要求:eps='); %%“|b-x|<=eps”
k=1;
x=(a+b)/2;
fprintf(' k a f(a) b f(b) xk f(xk)\n ');
T=[k,a,f(a),b,f(b),x,f(x)];
while abs(T(k,4)-T(k,6))>eps/2
k=k+1;
if f(x)*f(a)==0
a=a;
b=x;
x=(a+b)/2;
T=[T;k,a,f(a),b,f(b),x,f(x)];
break
elseif f(x)*f(a)>0
a=x;
b=b;
x=(a+b)/2;
T=[T;k,a,f(a),b,f(b),x,f(x)];
elseif f(x)*f(a)<0
a=a;
b=x;
x=(a+b)/2;
T=[T;k,a,f(a),b,f(b),x,f(x)];
end
end
disp(T);
fprintf('经过%d次迭代,函数方程根的近似解为:x=%.8f\n',k-1,T(k-1,6))
下面给出例题
例题1:
例题2:
希望对大家有所帮助,喜欢的朋友点个赞,加个关注,谢谢。
更多推荐
已为社区贡献1条内容
所有评论(0)