跳到主要内容

京都大学 情報学研究科 数理工学専攻 2023年8月実施 凸最適化

Author

Casablanca

Description

日本語版

English Version

Let c=(c1,c2,,cn)c=(c_1,c_2, \ldots , c_n) , where the superscript \top denotes transposition. Consider the following linear programming problem P:

P:Minimizecysubject toi=1nyi1 y0\begin{aligned} \text{P} : &\text{Minimize} &\boldsymbol{c}^{\top} \boldsymbol{y}\\ &\text{subject to} &\sum_{i=1}^{n} y_{i} \leq 1 \\ &\text{ } &\boldsymbol{y} \geq 0 \end{aligned}

where the decision variable of problem P is the vector y=(y1,y2,yn)R\boldsymbol{y} = (y_1, y_2, \ldots y_n)^\top \in \mathbb{R}^\top.

Answer the following questions (i) and (ii)

(i) write out a dual problem of problem P.

(ii) Show that problem P has an optimal solution.

Let YY be the set of optimal solutions of problem P.

Answer the following quetions (iii) and (iv).

(iii) Show that YY is a convex set.

(iv) Suppose that c1=c2==cnc_{1} = c_2 = \cdots = c_n and c1<0c_1 < 0. Consider the following optimization problem Q:

Q:Minimize12xxcxsubject toxY\begin{aligned} \text{Q} : &\text{Minimize} &\frac{1}{2} {\boldsymbol{x}^{\top} \boldsymbol{x} - \boldsymbol{c}^{\top} \boldsymbol{x}}\\ &\text{subject to} &\boldsymbol{x} \in Y \end{aligned}

where the decision variable of problem Q is the vector xRn\boldsymbol{x} \in \mathbb{R}^n. Obtain an optimal solution of problem Q by using Karush-Kuhn-Tucker conditions.

题目描述

c=(c1,c2,,cn)\boldsymbol c=(c_1,c_2,\ldots,c_n)^\top,其中 \top 表示转置。考虑以 y=(y1,,yn)\boldsymbol y=(y_1,\ldots,y_n)^\top 为变量的线性规划

P:最小化cy满足i=1nyi1,y0.\begin{aligned} \mathrm P:\quad &\text{最小化}\quad \boldsymbol c^\top\boldsymbol y\\ &\text{满足}\quad\sum_{i=1}^n y_i\le1,\qquad \boldsymbol y\ge0. \end{aligned}

回答:

  1. 写出 P 的对偶问题。
  2. 证明 P 有最优解。

YY 为 P 的最优解集合,继续回答:

  1. 证明 YY 是凸集。
  2. 假设 c1=c2==cnc_1=c_2=\cdots=c_nc1<0c_1<0。考虑
    Q:最小化12(xxcx)满足xY,\begin{aligned} \mathrm Q:\quad &\text{最小化}\quad \frac12\bigl(\boldsymbol x^\top\boldsymbol x- \boldsymbol c^\top\boldsymbol x\bigr)\\ &\text{满足}\quad\boldsymbol x\in Y, \end{aligned}
    其中 xRn\boldsymbol x\in\mathbb R^n。利用 KKT 条件求 Q 的一个最优解。

考点

  • 线性规划对偶与解的存在性:为单个资源不等式和非负变量构造对偶,并证明紧致可行域上的线性目标可达最优值。
  • 最优解集合的凸性:利用线性目标与凸可行域证明任意两个最优解的凸组合仍最优。
  • KKT 条件:在对称目标系数下先刻画最优面 YY,再求其上的严格凸二次函数最优解。

Kai

(i)

We have Lagrangian: L(y,λ,ν)=cy+λ(1y1)νyL(y,\lambda,\nu) = c^\top y + \lambda (\mathbf{1}^\top y - \mathbf{1} ) - \nu^\top y.

Obtain Lagrange dual function: d(λ,ν)=infy((c+λ1ν)yλ)=λd(\lambda, \nu) = \inf_{y} ((c^\top + \lambda \mathbf{1}^\top - \nu^\top)y - \lambda) = - \lambda.

Then we write the dual problem D:

D:Maximizeλsubject toc+λ1ν=0 0,ν0\begin{aligned} \text{D} : &\text{Maximize} &-\lambda\\ &\text{subject to} &c + \lambda \mathbf{1} - \nu = 0\\ &\text{ } &\geq 0, \nu \succeq 0 \end{aligned}

(ii)

Since λ1c-\lambda \mathbf{1} \preceq c, the optimal value vv of dual problem D satisfies: vmin{c1,c2,cn}v \leq \min \{ c_1, c_2, \ldots c_n \} then the linear programming PP is bounded, thus has an optimal solution.

(iii)

For any y1,y2Yy_1,y_2 \in Y , we know that for any y~\widetilde{y} which satisfies the constraints of P,

cy1cy~,cy2cy~c^\top y_1 \leq c^\top \widetilde{y}, \quad c^\top y_2 \leq c^\top \widetilde{y}

for θ[0,1]\theta \in [0,1],

θy1+(1θ)y20,\theta y_1 + (1 - \theta) y_2 \succeq \mathbf{0}, \quad
1(θy1+(1θ)y2)1\mathbf{1}^\top (\theta y_1 + (1-\theta) y_2) \leq 1
c(θy1+(1θ)y2)=θcy1+(1θ)cy2cy~c^\top (\theta y_1 + (1-\theta)y_2) = \theta c^\top y_1 + (1-\theta)c^\top y_2 \leq c^\top \widetilde{y}

then θy1+(1θ)y2Y\theta y_1 + (1-\theta)y_2 \in Y, YY is a convex set.

(iv)

Since c1=c2=cn<0c_1 = c_2 = \ldots c_n < 0, Y={yy0,1y=1}Y = \{ y | y \succeq \mathbf{0}, \mathbf{1}^\top y = 1 \}, we can rewrite Q as

Q:Minimize12xxcxsubject tox0,1x=1\begin{aligned} Q:&\text{Minimize} &\frac{1}{2} x^\top x - c^\top x \\ &\text{subject to} &x \succeq \mathbf{0}, \mathbf{1}^\top x = 1 \\ \end{aligned}

thus we get the Lagrangian: L(x,λ,μ)=12xxcxλx+μ(11x)L(x,\lambda, \mu) = \frac{1}{2} x^\top x - c^\top x - \lambda^\top x + \mu(1-\mathbf{1}^\top x). Then the KKT-condition:

KKT-conditions: {xcλμ1=0λ0,λx=0x01x=1\text{KKT-conditions: } \left\{ \begin{aligned} x - c - \lambda - \mu \mathbf{1} &= 0 \\ \lambda \succeq 0, -\lambda^\top x &=0 \\ x &\succeq 0\\ \mathbf{1}^\top x &= 1 \end{aligned} \right.

it is obviously that

x=[1n,1n,,1n],λ=0,μ=[1nc1,1nc1,,1nc1]x^* = [\frac{1}{n} , \frac{1}{n}, \ldots , \frac{1}{n}]^\top , \lambda = \mathbf{0}, \mu = [\frac{1}{n} - c_1, \frac{1}{n} - c_1, \ldots, \frac{1}{n} - c_1]^\top

satiesfies the KKT-condtions. Hence xx^* is an optimal solution.