跳到主要内容

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

Author

Casablanca

Description

日本語版

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

(i) 次の非線形計画問題を考える。

(P) Maximize  θ(x)subject to  xX\begin{aligned} \text{(P) } &\text{Maximize } \ \theta(\boldsymbol{x}) \\ &\text{subject to } \ \boldsymbol{x} \in X \end{aligned}

ただし、(P) の決定変数は xRnx \in \mathbb{R}^n であり、θ:RnR\theta : \mathbb{R}^n \rightarrow \mathbb{R}XRnX \subseteq \mathbb{R}^n は以下のように定義された目的関数と実行可能領域である。

θ(x)=(i=1nxi)1n,X={xRn|i=1nxi=1,xi0 (i=1,,n)}\theta(\boldsymbol{x}) = \left( \prod_{i=1}^{n} x_i \right)^{\frac{1}{n}}, \quad X = \left\{ \boldsymbol{x} \in \mathbb{R}^n \middle| \sum_{i=1}^{n} x_i = 1, \, x_i \geqq 0 \ (i = 1, \ldots, n) \right\}

問題 (P) は唯一の最適解 x\boldsymbol{x}^* を持ち、関数 θ\thetaR+n\mathbb{R}_{+}^n 上で凹関数(すなわち、θ-\theta は凸関数)であることが知られている。 ただし、R+n={xRnxi>0 (i=1,,n)}\mathbb{R}_{+}^n = \{ \boldsymbol{x} \in \mathbb{R}^n \mid x_i > 0 \ (i = 1, \ldots, n) \} である。

以下の (a), (b), (c)(c) に答えよ。

(a) 問題 (P) のカルーシュ・キューン・タッカー条件 (Karush-Kuhn-Tucker 条件) を書け。(問題 (P) が最大化問題であることに注意すること。)

(b) 問題 (P) の最適解 x\boldsymbol{x}^* を求めよ。

(c)(c) γiR,γi0 (i=1,,n)\gamma_i \in \mathbb{R}, \, \gamma_i \geqq 0 \ (i = 1, \ldots, n) とする。問題 (P) の最適解 x\boldsymbol{x}^* を利用して、以下の算術幾何平均の不等式が成り立つことを示せ。

1ni=1nγi(i=1nγi)1n\frac{1}{n} \sum_{i=1}^{n} \gamma_i \geqq \left( \prod_{i=1}^{n} \gamma_i \right)^{\frac{1}{n}}

(ii) 正の整数 nn に対して、Fn\mathcal{F}_nRn\mathbb{R}^n から R\mathbb{R} への非負の凸関数の集合とする。以下の (A), (B) に答えよ。

(A) fFnf \in \mathcal{F}_n が与えられたとき、関数 gf:RnRg_f : \mathbb{R}^n \rightarrow \mathbb{R}gf(x)=f(x)2 (xRn)g_f(\boldsymbol{x}) = f(\boldsymbol{x})^2 \ (\boldsymbol{x} \in \mathbb{R}^n) と定義する。そのとき、任意の fn=1Fnf \in \bigcup_{n=1}^{\infty} \mathcal{F}_n に対して、gfg_f が凸関数であることを示せ。

(B) 正の数 αR\alpha \in \mathbb{R}fFnf \in \mathcal{F}_n が与えられたとき、関数 hf,α:RnRh_{f,\alpha} : \mathbb{R}^n \rightarrow \mathbb{R}hf,α(x)=f(x)α (xRn)h_{f,\alpha}(\boldsymbol{x}) = f(\boldsymbol{x})^{\alpha} \ (\boldsymbol{x} \in \mathbb{R}^n) と定義する。 そのとき、すべての αα\alpha \geqq \alpha^*fn=1Fnf \in \bigcup_{n=1}^{\infty} \mathcal{F}_n に対して、hf,αh_{f,\alpha} が凸関数であるような最小な αR\alpha^* \in \mathbb{R} を求めよ。その際、α\alpha^* が最小であることを示せ。

English Version

Kai

(i)

(a)

(P):Minimize  θ(x)subject to  1x=1x0\begin{aligned} \text{(P)}: \text{Minimize } \ &-\theta (x) \\ \text{subject to } \ &\boldsymbol{1}^\top x = 1 \\ &x \succeq \boldsymbol{0} \end{aligned}

Lagrangian:

L(x,μ)=θ(x)+μ(1x1)L(x, \mu) = -\theta (x) + \mu (1^\top x - 1)
KKT-conditions {1n(Πjinxj)1n1μ=0,i=1,2,,n1x=1,x0\text{KKT-conditions } \left\{ \begin{aligned} &-\frac 1n (\Pi_{j\neq i}^{n} x_j)^{\frac 1n - 1} - \mu = 0, i = 1, 2, \ldots, n \\ &1^\top x = 1, x\succeq 0 \\ \end{aligned} \right.

(b)

xx^*, μ\mu ^* satisfied KKT-conditions if x=[1n,1n,,1n]x^* = [\frac 1n, \frac 1n, \ldots , \frac 1n]^\top, μ=1n\mu = -\frac 1n

(c)(c)

(i=1nγi)1n=(i=1nγi)1nγiγi=(i=1nγi(i=1nγi)n)1n(i=1nγi)1ni=1nγi(\prod_{i=1}^{n} \gamma_i )^{\frac 1n} = (\prod_{i=1}^{n} \gamma_i )^{\frac 1n} \frac{\sum \gamma_i}{\sum \gamma_i} = (\frac{\prod_{i=1}^{n} \gamma_i}{(\sum_{i=1}^{n}\gamma_i)^n})^{\frac 1n} (\sum_{i=1}^{n} \gamma_i) \leq \frac 1n \sum_{i=1}^{n}\gamma_i

(ii)

(A)

For any fn=1Fnf \in \bigcup_{n=1}^{\infty} \mathcal{F}_n, w.l.o.g, let f:RkRf : \mathbb{R}^k \rightarrow \mathbb{R} be an nonnegative function. Then

θgf(x1)+(1θ)gf(x2)=θf(x1)2+(1θ)f(x2)2\theta g_f( x_1) + (1-\theta)g_f(x_2) = \theta f( x_1)^2 + (1-\theta) f( x_2)^2
gf(θx1+(1θ)x2)=f(θx1+(1θ)x2)2(θf(x1)+(1θ)f(x2))2g_f(\theta x_1 + (1-\theta)x_2) = f(\theta x_1 + (1-\theta)x_2) ^ 2 \leq (\theta f( x_1) + (1-\theta) f( x_2)) ^2

and consider ϕ(θ)=gf(θx1+(1θ)x2)θgf(x1)(1θ)gf(x2)\phi(\theta) =g_f(\theta x_1 + (1-\theta)x_2) - \theta g_f( x_1) - (1-\theta)g_f(x_2), by calculating Δ\Delta , easily we see:

gf(θx1+(1θ)x2)θgf(x1)+(1θ)gf(x2)g_f(\theta x_1 + (1-\theta)x_2) \leq \theta g_f( x_1) + (1-\theta)g_f(x_2)

(B)

α=1\alpha ^* = 1

for α1\alpha \geq 1: θf(x1)α+(1θ)f(x2)α(θf(x1)+(1θ)f(x2))α\theta f(x_1)^{\alpha} + (1-\theta)f(x_2) ^{\alpha} \geq (\theta f(x_1) + (1-\theta)f(x_2))^{\alpha} since θf(x1)+(1θ)f(x2)f(θx1+(1θ)x2)\theta f(x_1) + (1-\theta)f(x_2) \geq f(\theta x_1 + (1-\theta)x_2), and tαt^{\alpha} increases for t>0t>0 then (θf(x1)+(1θ)f(x2))α(f(θx1+(1θ)x2))α=h(θx1+(1θ)x2)(\theta f(x_1) + (1-\theta)f(x_2))^{\alpha} \geq (f(\theta x_1 + (1-\theta)x_2))^{\alpha } = h(\theta x_1 + (1-\theta)x_2) thus hh is convex for α1\alpha \geq 1.

If α<1\alpha < 1, let f(x)=x1αf(x) = x_1^{\alpha}, easy to see hh is not convex. hence α=1\alpha^* = 1