erfenfa.asv

来自「共轭梯度法」· ASV 代码 · 共 26 行

function [b]=erfenfa(x0,y0,p)
xmax=6;ymax=6;
%ymin=-6;ymin=-6;
esp=1*10^(-6);
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&k<10)
if(forward(x1,y1)<forward(x2,y2))
    xmax=(x0+xmax)/2;
    ymax=(x0+xmax)/2;
else
    x0=(x0+xmax)/2;
    y0=(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;