利用二分法求方程在区间[-2,-1]的根3*X*X*X-4*X*X-5*X+13=0在区间[-2,-1]的根 20
二分法求解的中间结果如图所示:二分法求根区间:[-2,-1]第一次迭代:-1.500000第二次迭代:-1.750000第三次迭代:-1.625000.......第十六...
二分法求解的中间结果如图所示:
二分法求根区间:[-2,-1]
第一次迭代:-1.500000
第二次迭代:-1.750000
第三次迭代:-1.625000
.... ...
第十六次迭代:-1.548904
第十七次迭代:-1.548912
要求:编程实现解方程的二分法,并把求解的中间结果输出到一个文本文件中。 展开
二分法求根区间:[-2,-1]
第一次迭代:-1.500000
第二次迭代:-1.750000
第三次迭代:-1.625000
.... ...
第十六次迭代:-1.548904
第十七次迭代:-1.548912
要求:编程实现解方程的二分法,并把求解的中间结果输出到一个文本文件中。 展开
2个回答
展开全部
给你个例子,希望对你有帮助!
将下面的代码复制到M文件,保存成bisection.m文件
function [x,k]=demimethod(a,b,f,emg)
% a,b :求解区间的两个端点
% f :所求方程的函数名
% emg :精度指标
% x:所求近似解
% k: 循环次数
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=0;
while abs(b-a)>emg
if fab==0
x=(a+b)/2;
return;
elseif fa*fab<0
b=(a+b)/2;
else
a=(a+b)/2;
end
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=k+1;
end
x=(a+b)/2;
结果:
>> f=@(x)x.^3+2*x.^2+x-5;
>> demimethod(-2,2,f,10^-4)
ans =
1.1164
参考文章:
MATLAB--二分法求解方程程序(完整版)
http://hi.baidu.com/2008zhenghui/blog/item/4793594bf68b682809f7ef89.html
将下面的代码复制到M文件,保存成bisection.m文件
function [x,k]=demimethod(a,b,f,emg)
% a,b :求解区间的两个端点
% f :所求方程的函数名
% emg :精度指标
% x:所求近似解
% k: 循环次数
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=0;
while abs(b-a)>emg
if fab==0
x=(a+b)/2;
return;
elseif fa*fab<0
b=(a+b)/2;
else
a=(a+b)/2;
end
fa=feval(f,a);
fab=feval(f,(a+b)/2);
k=k+1;
end
x=(a+b)/2;
结果:
>> f=@(x)x.^3+2*x.^2+x-5;
>> demimethod(-2,2,f,10^-4)
ans =
1.1164
参考文章:
MATLAB--二分法求解方程程序(完整版)
http://hi.baidu.com/2008zhenghui/blog/item/4793594bf68b682809f7ef89.html
2009-12-16
展开全部
12.2..1.
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询