6.1范数逼近

  1. 基本的范数逼近问题
  2. 罚函数逼近

基本的范数逼近问题

凸优化第六章逼近与拟合 6.1范数逼近

其中凸优化第六章逼近与拟合 6.1范数逼近,且凸优化第六章逼近与拟合 6.1范数逼近凸优化第六章逼近与拟合 6.1范数逼近一种范数。

范数逼近问题的解有时又被称为凸优化第六章逼近与拟合 6.1范数逼近在范数凸优化第六章逼近与拟合 6.1范数逼近的近似解。

凸优化第六章逼近与拟合 6.1范数逼近表示问题的残差。

解释:

(1)几何解释:在A的列空间上找到一个在范数凸优化第六章逼近与拟合 6.1范数逼近下离b最近的点。

(2)估计的解释:假设y=Ax+v,y是测量值,v是噪声,x是待估计的参数向量。给定y=b,找到最好的估计值凸优化第六章逼近与拟合 6.1范数逼近,使得凸优化第六章逼近与拟合 6.1范数逼近最小。

(3)设计的解释:x是设计变量,b是期望得到的最好的结果,而Ax是实际的结果,找到最好的设计值凸优化第六章逼近与拟合 6.1范数逼近,使得凸优化第六章逼近与拟合 6.1范数逼近最小。

例子

(1)最小二乘逼近:

凸优化第六章逼近与拟合 6.1范数逼近

对x极小化目标函数,凸优化第六章逼近与拟合 6.1范数逼近,令其为0,得到凸优化第六章逼近与拟合 6.1范数逼近,解得凸优化第六章逼近与拟合 6.1范数逼近

(2)Chebyshev或极小极大逼近

凸优化第六章逼近与拟合 6.1范数逼近

问题试图极小化最大残差,故称为极小极大逼近问题。该问题可以描述为一个线性规划:

凸优化第六章逼近与拟合 6.1范数逼近

其中凸优化第六章逼近与拟合 6.1范数逼近,约束的意思就是Ax-b的任意分量都在凸优化第六章逼近与拟合 6.1范数逼近之间。这里凸优化第六章逼近与拟合 6.1范数逼近

(3)残差绝对值之和逼近

凸优化第六章逼近与拟合 6.1范数逼近

在估计领域,称为一种鲁棒估计器。

凸优化第六章逼近与拟合 6.1范数逼近

其中凸优化第六章逼近与拟合 6.1范数逼近,不同于极小极大逼近问题,这里凸优化第六章逼近与拟合 6.1范数逼近的分量不一定是相等的。

罚函数逼近

罚函数逼近问题具有形式:

凸优化第六章逼近与拟合 6.1范数逼近

其中r是残差,凸优化第六章逼近与拟合 6.1范数逼近是残差的第i个分量,凸优化第六章逼近与拟合 6.1范数逼近称为残差罚函数,设凸优化第六章逼近与拟合 6.1范数逼近是凸函数,凸优化第六章逼近与拟合 6.1范数逼近表示对每个残差的每个分量的费用或惩罚。在多数情况下,罚函数是对称,非负的,且满足凸优化第六章逼近与拟合 6.1范数逼近.

例子:

(1)带有四区的线性罚函数:凸优化第六章逼近与拟合 6.1范数逼近,即对小于a的残差不惩罚。

(2)对数障碍罚函数:凸优化第六章逼近与拟合 6.1范数逼近,对大于a的残差给与无穷大惩罚。

(3)二次罚函数:凸优化第六章逼近与拟合 6.1范数逼近

凸优化第六章逼近与拟合 6.1范数逼近

 

下图给出了选取矩阵凸优化第六章逼近与拟合 6.1范数逼近,然后比较凸优化第六章逼近与拟合 6.1范数逼近的一范数,二范数逼近结果,以及带有死区的罚函数(a=0.5)和对数障碍罚函数(a=1)的结果:

凸优化第六章逼近与拟合 6.1范数逼近

 

Huber罚函数

凸优化第六章逼近与拟合 6.1范数逼近

凸优化第六章逼近与拟合 6.1范数逼近时,是二次罚函数,而在凸优化第六章逼近与拟合 6.1范数逼近时,类似线性罚函数。

M=1时Huber罚函数的图像如下:

凸优化第六章逼近与拟合 6.1范数逼近

Huber罚函数在实际问题中效果非常好。

凸优化第六章逼近与拟合 6.1范数逼近

上图实线是Huber罚函数的拟合结果,而虚线是直线凸优化第六章逼近与拟合 6.1范数逼近对42个点的最小二乘拟合结果。可以看出最小二乘拟合结果偏离了实际的位置。这是因为在图中有两个野值凸优化第六章逼近与拟合 6.1范数逼近,且野值大于M,此时最小二乘对野值的惩罚是凸优化第六章逼近与拟合 6.1范数逼近,较大,所以会使原来拟合的直线想野值方向旋转,即从实线向虚线变化,而Huber罚函数对野值的惩罚较小,所以影响也较小,从而会得到很好的效果。

相关文章: