Broyden方法求解非线性方程组的Matlab实现.doc

Broyden方法求解非线性方程组的Matlab实现 注：matlab代码来自网络，仅供学习参考。 把以下代码复制在一个.m文件上 function [sol, it_hist, ierr] = brsola(x,f,tol, parms) % Broydens Method solver, globally convergent % solver for f(x) = 0, Armijo rule, one vector storage % % This code comes with no guarantee or warranty of any kind. % % function [sol, it_hist, ierr] = brsola(x,f,tol,parms) % % inputs: % initial iterate = x % function = f % tol = [atol, rtol] relative/absolute % error tolerances for the nonlinear iteration % parms = [maxit, maxdim] % maxit = maxmium number of nonlinear iterations % default = 40 % maxdim = maximum number of Broyden iterations % before restart, so maxdim-1 vectors are % stored % default = 40 % % output: % sol = solution % it_hist(maxit,3) = scaled l2 norms of nonlinear residuals % for the iteration, number function evaluations, % and number of steplength reductions % ierr = 0 upon successful termination % ierr = 1 if after maxit iterations % the termination criterion is not satsified. % ierr = 2 failure in the line search. The iteration % is terminated if too many steplength reductions % are taken. % % % internal parameter: % debug = turns on/off iteration statistics display as % the iteration progresses % % alpha = 1.d-4, parameter to measure sufficient decrease % % maxarm = 10, maximum number of steplength reductions before % failure is reported % % set the debug parameter, 1 turns display on, otherwise off % debug=1; % % initialize it_hist, ierr, and set the iteration parameters % ierr = 0; maxit=40; maxdim=39; it_histx=zeros(maxit,3); maxarm=10; % if nargin == 4 maxit=parms(1); maxdim=parms(2)-1; end rtol=tol(2); atol=tol(1); n = length(x); fnrm=1; itc=0; nbroy=0; % % evaluate f at the initial iterate % compu