二分法matlab编程代码

二分法的Matlab编程代码如下：

x=linspace(0,15);
y=(x-1).*(x-3).*(x-9).*(x-12);
plot(x,y);%画函数图像
%将坐标轴平移到原点
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
%将坐标轴平移到原点
xlabel('x');
ylabel('y');
syms X Y;
Y=@(X) (X-1).*(X-3).*(X-9).*(X-12);
a=2;b=6;e=0.01;%确定区间及收敛精度
sign=1;%循环进行标志
tic
while(sign==1)
c=(a+b)/2;
if(Y(c)==0)%零点即为c
fprintf('零点精确解为%f\n',c);
else
if(Y(a)*Y(c)<0)%零点在[a,c]之间
b=c;
else %零点在[c,b]之间
a=c;
end
end
if(abs(b-a)<=e|Y(c)==0)
sign=0;
fprintf('二分法求得零点解为%f\n',c);
else
sign=1;
end
end
toc

该代码实现了对函数 (x-1)*(x-3)*(x-9)*(x-12) 的零点的求解，其中 a 和 b 分别表示区间的左右端点，e 表示收敛精度。程序通过不断进行区间压缩来获取函数零点，直到区间长度小于等于精度要求。