跳到主要内容

京都大学 情報学研究科 数理工学専攻 2007年8月実施 オペレーションズ・リサーチ

Author

思齐塾, 祭音Myyura

Description

以下の問(i), (ii)に答えよ。

(i) f:RnRf: \mathbb{R}^n \to \mathbb{R}f(0)=0f(0) = 0 である凸関数とする。以下の(a), (b)に答えよ。

(a) すべての t1t \geq 1xRnx \in \mathbb{R}^n に対して次の不等式が成り立つことを示せ。

f(tx)tf(x)f(tx) \geq tf(x)

(b) すべての xRnx \in \mathbb{R}^n に対して f(x)Mf(x) \leq M となる定数 MM が存在するとき、すべての xRnx \in \mathbb{R}^n に対して f(x)=0f(x) = 0 であることを示せ。

(ii) uunn 次元ベクトル、 QQn×nn \times n 正定値対称行列とする。 xRnx \in \mathbb{R}^n を決定変数とする次の非線形計画問題を考える。

(P):

minimizeuTxsubject toxTQxuTQu\begin{aligned} &\text{minimize} & u^T x \\ &\text{subject to} & x^T Q x \leq u^T Q u \end{aligned}

ここで TT は転置記号を表す。次の(A), (B)に答えよ。

(A) カルーシュ・キューン・タッカー条件 (Karush-Khun-Tucker 条件) を用いて、問題 (P) の最適解と目的関数の最小値を求めよ。

(B) 上の結果を用いて、次の不等式が成り立つことを示せ。

(uTu)2(uTQu)(uTQ1u)(u^T u)^2 \leq (u^T Q u)(u^T Q^{-1} u)

题目描述

回答以下两部分问题。

  1. f:RnRf:\mathbb R^n\to\mathbb R 为凸函数且 f(0)=0f(0)=0

    1. 证明对任意 t1t\geq1xRnx\in\mathbb R^n,都有

      f(tx)tf(x).f(tx)\geq t f(x).
    2. 若存在常数 MM,使得所有 xRnx\in\mathbb R^n 均满足 f(x)Mf(x)\leq M,证明 ff 在整个 Rn\mathbb R^n 上恒等于零。

  2. uunn 维向量,QQn×nn\times n 正定对称矩阵,以 xRnx\in\mathbb R^n 为决策变量考虑非线性规划

    minimizeuTx,subject toxTQxuTQu,\begin{aligned} \text{minimize}\quad &u^Tx,\\ \text{subject to}\quad &x^TQx\leq u^TQu, \end{aligned}

    其中 TT 表示转置。

    1. 使用 Karush–Kuhn–Tucker 条件求该问题的最优解和目标函数最小值。

    2. 利用所得结果证明

      (uTu)2(uTQu)(uTQ1u).(u^Tu)^2\leq(u^TQu)(u^TQ^{-1}u).

Kai

(i) (a) Since ff is a convex function, for t1t \geq 1 , let λ=1/t\lambda = 1/t , then 0<λ10 < \lambda \leq 1 . Since f(0)=0f(0) = 0 , we can write:

f(x)=f(λ(tx)+(1λ)0)λf(tx)+(1λ)f(0)=λf(tx)f(x) = f(\lambda (tx) + (1-\lambda)0) \leq \lambda f(tx) + (1-\lambda) f(0) = \lambda f(tx) .

Therefore, f(x)1tf(tx)f(x) \leq \frac{1}{t} f(tx) , which means tf(x)f(tx)tf(x) \leq f(tx) .

(i) (b) (a) より、もしある xx について f(x)>0f(x)>0 なら

f(tx)tf(x)(t1)f(tx)\geq t f(x)\qquad(t\geq1)

である。十分大きな tt を取れば右辺は MM を超え、上界 f(tx)Mf(tx)\leq M に矛盾する。したがって、すべての xx について f(x)0f(x)\leq0 である。

一方、凸性と f(0)=0f(0)=0 より

0=f(0)=f(x+(x)2)f(x)+f(x)2.0=f(0)=f\left(\frac{x+(-x)}2\right) \leq\frac{f(x)+f(-x)}2.

f(x),f(x)0f(x),f(-x)\leq0 なので、これは f(x)=f(x)=0f(x)=f(-x)=0 を意味する。よってすべての xx に対して f(x)=0\boxed{f(x)=0} である。

(ii) (A) u=0u=0 のとき、制約は xTQx0x^TQx\leq0 であり、 QQ が正定値なので実行可能解は x=0x=0 だけである。従って最適解は x=0x^*=0 、最小値は 00 である。

以下、 u0u\ne0 とする。The Lagrangian of (P) is

L(x,λ)=uTx+λ(xTQxuTQu).L(x, \lambda) = u^T x + \lambda(x^T Q x - u^T Q u).

The KKT conditions are:

L(x,λ)=u+2λQx=0λ(xTQxuTQu)=0λ0xTQxuTQu\begin{aligned} \nabla L(x, \lambda) &= u + 2\lambda Q x = 0 \\ \lambda(x^T Q x - u^T Q u) &= 0 \\ \lambda &\geq 0 \\ x^T Q x &\leq u^T Q u \end{aligned}

停留条件から λ=0\lambda=0 は不可能なので λ>0\lambda>0 であり、相補性により制約は等号で成立する。From u+2λQx=0u + 2\lambda Q x = 0 , we have x=12λQ1ux = -\frac{1}{2\lambda} Q^{-1} u . Substituting this into xTQx=uTQux^T Q x = u^T Q u , we get

(12λQ1u)TQ(12λQ1u)=uTQu\left(-\frac{1}{2\lambda} Q^{-1} u\right)^T Q \left(-\frac{1}{2\lambda} Q^{-1} u\right) = u^T Q u
14λ2uTQ1QQ1u=uTQu\frac{1}{4\lambda^2} u^T Q^{-1} Q Q^{-1} u = u^T Q u
14λ2uTQ1u=uTQu\frac{1}{4\lambda^2} u^T Q^{-1} u = u^T Q u
λ2=uTQ1u4uTQuλ=12uTQ1uuTQu\lambda^2 = \frac{u^T Q^{-1} u}{4 u^T Q u} \Rightarrow \lambda = \frac{1}{2} \sqrt{\frac{u^T Q^{-1} u}{u^T Q u}}

Substituting x=12λQ1ux = -\frac{1}{2\lambda} Q^{-1} u into the objective function uTxu^T x , we get

uTx=uT(12λQ1u)=12λuTQ1u=uTQuuTQ1uuTQ1u=(uTQu)(uTQ1u)u^T x = u^T \left(-\frac{1}{2\lambda} Q^{-1} u\right) = -\frac{1}{2\lambda} u^T Q^{-1} u = - \sqrt{\frac{u^T Q u}{u^T Q^{-1} u}} u^T Q^{-1} u = -\sqrt{(u^T Q u)(u^T Q^{-1} u)}

Thus the optimal solution is

x=uTQuuTQ1uQ1u\boxed{x^*=-\sqrt{\frac{u^TQu}{u^TQ^{-1}u}}\,Q^{-1}u}

and the minimum value is

(uTQu)(uTQ1u).\boxed{-\sqrt{(u^TQu)(u^TQ^{-1}u)}}.

目的関数は線形、制約関数は凸であり、 x=0x=0 は制約を狭義に満たすため、KKT 条件は十分でもある。

(ii) (B) u=0u=0 なら不等式は明らかである。 u0u\ne0 のとき、 x=ux=-u

(u)TQ(u)=uTQu(-u)^TQ(-u)=u^TQu

を満たす実行可能解である。従って (A) の最小値はこの点の目的値以下であり、

(uTQu)(uTQ1u)uT(u)=uTu.-\sqrt{(u^TQu)(u^TQ^{-1}u)} \leq u^T(-u)=-u^Tu.

よって

uTu(uTQu)(uTQ1u).u^Tu\leq\sqrt{(u^TQu)(u^TQ^{-1}u)}.

両辺は非負なので二乗して

(uTu)2(uTQu)(uTQ1u)\boxed{(u^Tu)^2\leq(u^TQu)(u^TQ^{-1}u)}

を得る。