mian.m

罚函数法matlab程序.希望对大家有用.

%主函数main.m
clc
clear
a=input('请输入根的上限：');
b=input('请输入根的下限：');
e=input('请输入控制精度：');

x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
count=1;
if abs(b-a)>e
    while abs(b-a)>e
         count=count+1;
         if f(x1)>f(x2)
          a=x1;
          b=b;
          x1=x2;
          x2=a+0.618*(b-a);
         else
          a=a;
          b=x2;
          x2=x1;
          x1=a+0.382*(b-a);
         end
    end
    x=1/2*(a+b);
else
    x==1/2*(a+b);
end

fprintf('黄金分割法的解：\n')
x
fprintf('迭代次数：\n')
count