2楼: Originally posted by 月只蓝 at 2015-06-24 21:14:11
问题4~8如果调用MATLAB现有函数可以很快解决，不过现在题目对算法都规定死了，只能根据这一章学习的内容，依葫芦画瓢写代码。

2楼: Originally posted by 月只蓝 at 2015-06-24 21:14:11
问题4~8如果调用MATLAB现有函数可以很快解决，不过现在题目对算法都规定死了，只能根据这一章学习的内容，依葫芦画瓢写代码。

function erfenfa_solve
clear all;clc
i=0; %二分次数记数
a=1; %求根区间左端
b=2; %求根区间右端
fa=f(a); %计算f(a)的值
fb=f(b); %计算f(b)的值
c=(a+b)/2; %计算区间中点
fc=f(c); %计算区间中点f(c)
while abs(fc)>=1e-2; %判断f(c)是否为零点
if fa*fc>=0; %判断左侧区间是否有根
fa=fc;
a=c;
else fb=fc;
b=c;
end
c=(a+b)/2;
fc=f(c);
i=i+1;
end

fprintf('\n\t%s%.2f\n\t%s%d\n\t%s%.3f','方程根x=',c,'迭代次数i=',i,'方程残差fval=',f(c)) %计算结果输出

function y=f(x)
y=2*x^3-5*x-1; %函数f(t)的表达式

方程根x=1.67
        迭代次数i=8
        方程残差fval=0.010

fsolve(@(x) 2*x^3-5*x-1,1)

ans =

    1.6730