跳转至

对偶理论

考虑一个基本的优化问题:

minxRnf(x)s.t.ci(x)=0,iEci(x)0,iI \begin{aligned} &\min_{x \in \mathbb{R}^n} && f(x) \\ &\text{s.t.} && c_i(x) = 0, & i \in \mathcal{E}\\ &&& c_i(x) \leqslant 0, & i \in \mathcal{I} \end{aligned}

下面我们将讨论该类型优化问题的拉格朗日对偶问题。

拉格朗日函数

考虑其拉格朗日函数,其基本思想是为每个约束指定一个拉格朗日乘子,以乘子为加权系数将约束增加到目标函数中:

L(x,λ,μ)=f(x)+iIλici(x)+iEμici(x),μi0 \begin{align*} L(x, \lambda, \mu) = f(x) + \sum_{i \in \mathcal{I}}\lambda_i c_i(x) + \sum_{i \in \mathcal{E}} \mu_i c_i(x), \quad \mu_i \geqslant 0 \end{align*}

上式约束条件的目的是为了保证 f(x)f(x) 在原优化问题定义下的可行点处的取值大于或等于拉格朗日函数的取值。即拉格朗日函数是原优化问题的一个下界

对偶函数

对拉格朗日函数 L(x,λ,μ)L(x, \lambda, \mu) 中的 xx 取下确界可以定义拉格朗日对偶函数,这一函数在对偶理论中起到关键性作用。

拉格朗日对偶函数

拉格朗日对偶函数 g:R+m×Rp[,+) g: \mathbb{R}^m_+ \times \mathbb{R}^p \to [-\infty, +\infty ) 是拉格朗日函数 L(x,λ,μ)L(x, \lambda, \mu) 对于 λR+m,μRp\lambda \in \mathbb{R}^m_+, \mu \in \mathbb{R}^p 关于 xx 的下确界:

g(λ,μ)=infxRnL(x,λ,μ) g(\lambda, \mu) = \inf_{x \in \mathbb{R}^n} L(x, \lambda, \mu)

弱对偶定理

根据上面的定义,可以很容易得出一个结论即

弱对偶定理

对于任意的 μ0\mu \geqslant 0λ\lambda,拉格朗日对偶函数给出了优化问题

minxRnf(x)s.t.ci(x)=0,iEci(x)0,iI \begin{aligned} &\min_{x \in \mathbb{R}^n} && f(x) \\ &\text{s.t.} && c_i(x) = 0, & i \in \mathcal{E}\\ &&& c_i(x) \leqslant 0, & i \in \mathcal{I} \end{aligned}

最优值的一个下界,即

g(λ,μ)p,μ0 g(\lambda, \mu) \leqslant p^*, \mu \geqslant 0

对偶问题

那么一个自然的问题就是,从拉格朗日对偶函数获得的下界中,哪个是最优的?为了求解该最优的下界,便有如下拉格朗日对偶问题:

maxλ,μ0g(λ,μ)=maxλ,μ0infxRnL(x,λ,μ) \max_{\lambda, \mu \geqslant 0} g(\lambda, \mu) = \max_{\lambda, \mu \geqslant 0} \inf_{x \in \mathbb{R}^n} L(x, \lambda, \mu)

即对偶问题是最大化拉格朗日函数的下确界

强对偶原理

当固定 (λ,μ)(\lambda, \mu) 时,拉格朗日函数关于 xx 可能无界,那么对偶函数在 (λ,μ)(\lambda, \mu) 处的取值为 -\infty,此时对偶函数提供的下界时无意义的。因此我们规定拉格朗日对偶函数的定义域为

domg={(λ,μ)μ0,g(λ,μ)>} \mathbf{dom } g = \{ (\lambda, \mu) \mid \mu \geqslant 0, g(\lambda, \mu) > -\infty \}

(λ,μ)domg(\lambda, \mu) \in \mathbf{dom } g 时,称其为对偶可行解。记对偶问题的最优值为 qq^*,如果有 p=qp^* = q^*,则称强对偶原理成立。

线性规划的对偶问题

下面我们具体地看看线性规划中的对偶问题。考虑线性规划

minxcxs.t.Ax=bx0 \begin{aligned} &\min_{x} && c^\top x \\ &\text{s.t.} && Ax = b\\ &&& x \geqslant 0 \end{aligned}

写出拉格朗日函数

L(x,λ,μ)=cx+λ(Axb)μx,μ0 L(x, \lambda, \mu) = c^\top x + \lambda^\top (Ax - b) - \mu^\top x, \mu \geqslant 0

固定 (λ,μ)(\lambda, \mu),则 L(x,λ,μ)L(x, \lambda, \mu) 是一个关于 xx 的一次方程,求导有

xL(x,λ,μ)=c+Aλμ \nabla_x L(x, \lambda, \mu) = c + A^\top \lambda - \mu

因此原问题对偶函数即拉格朗日函数下确界为

g(λ,μ)={λb,Aλ+cμ=0,μ0,others. \begin{align*} g(\lambda, \mu) = \begin{cases} -\lambda^\top b, & A^\top \lambda + c - \mu = 0, \mu \geqslant 0 \\ -\infty, & \text{others.} \end{cases} \end{align*}

得到对偶问题

maxλ,μλbs.t.Aλ+c=μμ0 \begin{aligned} &\max_{\lambda, \mu} && -\lambda^\top b \\ &\text{s.t.} && A^\top \lambda + c = \mu\\ &&& \mu \geqslant 0 \end{aligned}