levenberg_Marquardt算法

%% 函数功能:利用Levenberg_Marquardt算法求解无约束的非线性最小二乘问题  
function [x_curr iter] = Levenberg_Marquardt（Difff var x0 lamda beita eps）  %eps为精度，
    % 相关参数  
    max_iter    = 100;  
    flag        = 1;  
    x_curr      = x0;                                                 
  
    for iter = 1 : max_iter  
        if flag  
            J   = double（subs（jacobian（Diff var） var x_curr））;      % 雅可比矩阵  
            H   = transpose（J） * J;                                    % 赫森矩阵      
              
            if iter == 1  
                f_curr = double（subs（f var x_curr））;  
            end  
        end  
  
        % 计算新Hessen矩阵  
        H_new = H + （lamda * eye（8 8））;  
  
        % 计算步长d  
        diff = double（subs（Diff var x_curr））;