• 拉格朗日对偶性


    0 前言

    本文承接上一篇博文拉格朗日乘子法和KKT条件http://www.cnblogs.com/liaohuiqiang/p/7805954.html讲讲拉格朗日对偶性的问题。

    在约束优化问题中,常常用拉格朗日对偶性来将原始问题转为对偶问题,通过解对偶问题的解来得到原始问题的解。

    1 为什么要利用对偶?

    首先要明确,对偶问题的解不一定直接等于原问题的解(弱对偶),但是,对偶问题有两点性质。

    1.1 满足某些条件时,对偶问题直接等于原问题的解(强对偶)

    1.2 无论原始问题是否是凸的,对偶问题都是凸优化问题

    显然,在某些情况下,直接对对偶问题求解可以得到原问题的解,而且对偶问题是凸优化,易于求解。所以利用对偶来求解是很有用的。

    1 原始问题

    考虑原始的优化问题如下

    先定义拉格朗日函数$ L(x, alpha , eta ) = f(x) + sum _{i=1}^q alpha_i g_i(x) + sum _{j=q+1}^m eta_j h_j(x) $

    考虑函数$$Theta_P(x) = max_{alpha, eta: alpha_igeqslant 0}L(x, alpha , eta ) $$

    显然满足约束时函数等于f(x),不满足约束时函数等于正无穷,容易得到$ Theta_P(x) = egin{Bmatrix}f(x), x satisfy constraints\ {+ infty, others }end{Bmatrix} $

    极小化该函数得到的$$min_x  Theta_P(x) =min_x max_{alpha, eta: alpha_igeqslant 0}L(x, alpha , eta ) $$和原始的优化问题是等价。

    同时,我们定义原始问题的最优值为$$ p^*=min_x  Theta_P(x) $$

    2 对偶问题

    考虑函数$$Theta_D(alpha , eta) = min_x L(x, alpha , eta ) $$

    极大化该函数得到$$max_{alpha, eta: alpha_igeqslant 0} Theta_D(alpha ,eta) =max_{alpha, eta:alpha_igeqslant 0}min_x L(x, alpha , eta ) $$我们把这个极大化问题视为对偶问题。

    同时,我们定义对偶问题的最优值为$$ d^*=max_{alpha, eta: alpha_igeqslant 0} Theta_D(alpha ,eta) $$

    3 原始问题和对偶问题的关系

    可以得知$ d^*leqslant p^* $,因为p是先求最大的一块区域然后在这块区域求最小,d是先求最小的一块区域然后在这块区域求最大,最大里面的最小,总会比最小里面的最大要大(或等于)。

    $ d^*leqslant p^* $,称为弱对偶,对于所有优化问题都成立,这个时候我们可以得到原始问题的一个下界。

    如果$  d^*= p^* $,称为强对偶,满足某些条件才成立,这时可以用解对偶问题替代原始问题。那么满足什么样的条件可以得到强对偶呢?

    如果原问题是一个凸优化,并且不等式约束g(x)是严格可行的,即存在x对所有i有$ g_i(x)<0 $,那么强对偶成立。这里需要注意的是,这里的条件只是强对偶成立的一种情况,对于非凸的问题也有可能是强对偶。

    强对偶成立时,将拉格朗日函数分别对原变量x和对偶变量α和β分别求导,令导数等于零(还需要满足KKT条件),即可求解对偶问题的解,也就求得了原问题的解。

    4 对偶求解和拉格朗日乘子法求解

    对于拉格朗日乘子法求解,当原问题是凸优化的时候,考虑满足KKT的条件,对于拉格朗日函数求导等于0即可得到解。

    对于对偶问题求解,转换为对偶函数求解,为了得到原始问题的解,要求强对偶,而在强对偶(不一定满足原问题凸优化)的情况下,考虑满足KKT的条件,对于拉格朗日函数求导等于0即可得到解。

    当原问题是凸优化的时候,强对偶和KKT条件是互为充要条件的。

  • 相关阅读:
    取出某个月有几天
    重建索引时,一些数值
    java代理概念
    java内部类和静态内部类
    Lamdba表达式的代码使用讲解
    java 中Vector的使用详解
    mysql 安装失败 start service执行不下去
    Mysql中ENGINE=InnoDB AUTO_INCREMENT=2 DEFAULT CHARSET=utf8怎么转换为sql sever2008的代码
    通过命令行发送邮件
    Servlet开发总结(一)
  • 原文地址:https://www.cnblogs.com/liaohuiqiang/p/7818448.html
Copyright © 2020-2023  润新知