【Matlab算法】L-M法求解非线性最小二乘优化问题（附L-M法MATLAB代码）

前言

博主一头小山猪目前已开放所有文章：小山猪——经典算法专栏

活动地址：CSDN21天学习挑战赛

正文

L-M法求解非线性最小二乘优化问题

L-M法 (Levenberg-Marquardt法)原理
当矩阵 ${\left(J_k\right)}^T{J}_k$ 为病态矩阵时，用G-N算法可能得不到正确的解，甚至当 ${\left(J_k\right)}^T{J}_k$ 不可逆时， 这时G-N算法就无法计算下去。L-M算法通过采用系数矩阵阻尼的方法改造矩阵 ${\left(J_k\right)}^T{J}_k$ 的性 态，使算法能够进行下去。
L-M算法中有两个主要的步骤: 一是解方程 $[Q+\mu I]\Delta x=-\nabla S\left(x^{(k)}\right)$ (其中 $\mu$ 为阻尼系 数）求出自变量的增量，另一个是阻尼系数 $\mu$ 的调整算法。
算法步骤
用L-M法求解非线性最小二乘优化问题 $\min S(x)$ 的算法过程如下:
【1】给定初始点 $x^{(0)}$ ，选取参数 $\beta \in (0,1)，\mu >1$ 及精度 $\epsilon >0$ ，置 $k=0$ ；
【2】计算 $f\left(x^{(k)}\right),S\left(x^{(k)}\right)$ ；
【3】计算 $\nabla f\left(x^{(k)}\right)$ ；
【4】计算 $\nabla S\left(x^{(k)}\right)={\left(\nabla f\left(x^{(k)}\right)\right)}^T\ast f\left(x^{(k)}\right)$ ；
【5】令 $Q={\left(\nabla f\left(x^{(k)}\right)\right)}^T\nabla f\left(x^{(k)}\right)$ 解方程 $[Q+\mu I]\Delta x=-\nabla S\left(x^{(k)}\right)$ ；
【6】令 $x^{(k+1)}=x^{(k)}+\Delta x$ ，计算终止条件是否满足，不满足转【7】；
【7】若 $S\left(x^{(k+1)}\right)<S\left(x^{(k)}\right)+\beta {\left(\nabla S\left(x^{(k)}\right)\right)}^T\Delta x$ ，令 $\mu =\mu /v$ ，转【8】，否则令 $\mu =\mu \ast v$ ，转【5】;
【8】令 $k=k+1$ ，转【2】。

代码实现

L-M法函数如下

function [x,minf] = minLM(f,x0,beta,u,v,var,eps)
format long;
if nargin == 6
    eps = 1.0e-6;
end
S = transpose(f)*f;
k = length(f);
n = length(x0);
x0 = transpose(x0);
A = jacobian(f,var);
tol = 1;

while tol>eps
    Fx = zeros(k,1);
    for i=1:k
        Fx(i,1) = Funval(f(i),var,x0);
    end
    Sx = Funval(S,var,x0);
    Ax = Funval(A,var,x0);
    gSx = transpose(Ax)*Fx;
    Q = transpose(Ax)*Ax;
    
    while 1
        dx = -(Q+u*eye(size(Q)))\gSx;

        x1 = x0 + dx;
        for i=1:k
            Fx1(i,1) = Funval(f(i),var,x1);
        end
        Sx1 = Funval(S,var,x1);
        tol = norm(dx);
        if tol<=eps
            break;
        end

        if Sx1 >= Sx+beta*transpose(gSx)*dx
            u = v*u;
            continue;
        else
            u = u/v;
            break;
        end
    end
    x0 = x1;
end
x = x0;
minf = Funval(S,var,x);
format short;
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48

参考资料

1. 小山猪——经典算法专栏
   
2. 数据结构教程 第5版