erfenfa.m

来自「共轭梯度法」· M 代码 · 共 29 行

function [b]=erfenfa(x0,y0,p)
xmax=6;ymax=6;
%ymin=-6;ymin=-6;
esp=1*10^(-4);
x1=p(1,1)*(x0+xmax)/4;
y1=p(2,1)*(y0+ymax)/4;
x2=3*p(1,1)*(x0+xmax)/4;
y2=3*p(1,1)*(y0+ymax)/4;
% k=1;
while (forward(x1,y1)-forward(x2,y2)>esp)
if(forward(x1,y1)<forward(x2,y2))
    xmax=p(1,1)*(x0+xmax)/2;
    ymax=p(2,1)*(x0+xmax)/2;
else
    x0=p(1,1)*(x0+xmax)/2;
    y0=p(2,1)*(x0+xmax)/2;
end
  x1=p(1,1)*(x0+xmax)/4;
  y1=p(2,1)*(y0+ymax)/4;
  x2=3*p(1,1)*(x0+xmax)/4;
  y2=3*p(1,1)*(y0+ymax)/4;
%   k=k+1;
end
b=zeros(2,1);
% b(1,1)=x1;
% b(2,1)=y1;
 b(1,1)=p(1,1)*(x0+xmax)/2;
 b(2,1)=p(2,1)*(x0+xmax)/2;