【MQ翻译】非线性最小二乘问题的求解方法:1. Introduction and Definitions

一直想看非线性最小二乘问题,但是每次都看得“不求甚解”、“零零散散”,刚好在网上下载到了一个英文版的教程:METHODS FOR NON-LINEAR LEAST SQUARES PROBLEMS(点此下载.pdf),感觉介绍的还不错,内容也没有特别长,所以就自己翻译了一版,供大家交流,如有错误还请指正。


在本册内容中,我们讨论以下问题:

定义1.1 最小二乘法

找到一个局部极小值 \mathbf{x^*},使得对于 \mathbf{R^n} \mapsto \mathbf{R} 的函数 f_i  :

其中,f_i (i=1,2,...,m) 是给定的函数,m \geqslant n

例1:数据拟合是一类非常重要的最小二乘问题。如图1.1所示的数据点(t_1,y_1),...,(t_m,y_m),我们可以进一步地给出它的拟合模型:

这个模型包含四个参数:\mathbf{x}=\left [ x_1,x_2,x_3,x_4 \right ]^T 。我们假设存在一个\mathbf{x^\dagger }使得:

其中, \left \{e_i}{ \right \} 表示纵坐标误差,类似“白噪音”。

 

对任何一组\mathbf{x},我们都可以计算出残差(residuals):

f_i(\mathbf{x})=y_i-M(\mathbf{x},t_i)=y_i-x_3e^{x_1t_i}-x_4e^{x_2t_i}, i=1,2,...,m

对于最小二乘拟合,我们通过使得残差平方之和最小来确定\mathbf{x^*}M(\mathbf{x^*},t) 的图像如图1.1中的实线所示。

最小二乘问题是一般问题的一个特殊变体:给定一个函数F(\mathbf{R^n} \mapsto \mathbf{R}),找到F的一组参数,使得目标函数(objective function)或成本函数(cost function)的值最小。

定义1.2 全局最小值(Global Minimizer)

总体来说,这个问题很难解决,我们只能将其简化为一个寻找F的局部最小值的问题,也就是说,在一个给定大小为 \delta 的区域内找到一个参数向量,使得F 在该区域内有最小值,其中 \delta 是一个非常小的正数。

定义1.3 局部最小值(Local Minimizer)

在接下来的介绍中,我们将讨论优化中的一些基本概念,第2章简要回顾了为一般成本函数寻找局部最小值的方法。更多详情请参考Frandsen等人(2004)的研究论文。在第三章中,我们给出了专门针对最小二乘问题的方法。

我们假设目标函数F是光滑可微的,那么就可以使用下面的泰勒展开式(Taylor expansion):

         (1.4a)

其中,g表示梯度:

         (1.4b)

H 代表海森矩阵(Hessian):

     (1.4c)

如果 \mathbf{x^*} 是一个局部极小值点且 \left \| \mathbf{h} \right \| 非常小,那么点\mathbf{​{x^*}+h} 处的F值不会比现在的F值更小。把这个结果与(1.4a)结合起来,我们可以得到:

定理1.5 局部最小值的必要条件

如果 \mathbf{x^*} 是一个局部极小值点,那么有:

我们用一个特别的名字来命名满足上述条件的点:

定义1.6 驻点

如果点 \mathbf{​{x_s}} 处满足

那么 \mathbf{​{x_s}} 就被称为F 的驻点。

因此,局部极小化也是一个驻点,但局部极大化也是一个驻点。既不是局部极大值也不是局部极小值的驻点称为鞍点(saddle point)。为了确定给定的驻点是否是局部极小值点,我们需要在泰勒级数(1.4a)中包含二阶项。

 (1.7)

根据海森矩阵的定义(1.4c),H 是对称矩阵。如果我们要求H_s 是正定(positive definite)的,那么它的特征值大于某个大于0的数\delta,并且有:

这表明,在\left \| \mathbf{h} \right \|足够小的情况下,式(1.7)右边第三项将由第二项决定。所以我们得到:

定理1.8 局部最小值的充分条件

如果驻点 \mathbf{​{x_s}} 满足 F^{''} { (\mathbf{​{x_s}})} 是正定的,那么\mathbf{​{x_s}} 是一个局部极小值点。

如果H_s是负定的(negative definite),那么\mathbf{x_s} 是一个局部极大值。如果H_s是不定矩阵(特征值既有正数又有复数),那么\mathbf{x_s} 是一个鞍点。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值