broydeniterate.m

来自「基于Matlab的计算方法中常用的算法程序」· M 代码 · 共 27 行

function s=BroydenIterate(x,eps)
  % 秩1的拟Newton法求非线性方程组
  % x为迭代初值，eps为允许误差值
  if nargin==1
      eps=1.0e-6;
  elseif nargin<1
      error
      return
  end
  % 第一次迭代
  x1=fx2(x);
  b1=inv(dfx2(x));
  p=-b1*x1';
  x3=x'+p;
  d=fx2(x3');
  q=(d-x1)';
  b1=b1+(q-b1*p)*p'/norm(p);
  while(norm(p)>=eps) %循环迭代
     % b1=b1+b3;
       x1=fx2(x3');
       p=-inv(b1)*(x1)';
       x3=x3+p;
       q=((fx2(x3'))-x1)';
       b1=b1+(q-b1*p)*p'/norm(p);
  end
  s=x3;
  return