跳到主要内容

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

Author

Casablanca, 祭音Myyura

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

题目描述

回答以下两大题。

  1. 考虑以 xRn\boldsymbol{x}\in\mathbb{R}^n 为决策变量的非线性规划
(P):最大化θ(x)=(i=1nxi)1/n约束于xX={xRn | i=1nxi=1,xi0 (i=1,,n)}.\begin{aligned} (P):\quad &\text{最大化}\quad \theta(\boldsymbol{x}) =\left(\prod_{i=1}^n x_i\right)^{1/n}\\ &\text{约束于}\quad \boldsymbol{x}\in X =\left\{\boldsymbol{x}\in\mathbb{R}^n\ \middle|\ \sum_{i=1}^n x_i=1,\quad x_i\geqq0\ (i=1,\ldots,n)\right\}. \end{aligned}

已知 (P)(P) 有唯一最优解 x\boldsymbol{x}^*,且 θ\theta

R+n={xRnxi>0 (i=1,,n)}\mathbb{R}_+^n =\left\{\boldsymbol{x}\in\mathbb{R}^n\mid x_i>0\ (i=1,\ldots,n)\right\}

上为凹函数(等价地,θ-\theta 为凸函数)。完成下列各问:

  1. 写出 (P)(P) 的 KKT 条件;注意这是最大化问题。
  2. (P)(P) 的最优解 x\boldsymbol{x}^*
  3. 给定 γiR\gamma_i\in\mathbb{R}γi0\gamma_i\geqq0i=1,,ni=1,\ldots,n),利用上述 x\boldsymbol{x}^* 证明算术—几何平均不等式
1ni=1nγi(i=1nγi)1/n.\frac1n\sum_{i=1}^n\gamma_i \geqq \left(\prod_{i=1}^n\gamma_i\right)^{1/n}.
  1. 对每个正整数 nn,令 Fn\mathcal{F}_n 表示所有从 Rn\mathbb{R}^nR\mathbb{R} 的非负凸函数所成的集合。

    1. fFnf\in\mathcal{F}_n 定义 gf:RnRg_f:\mathbb{R}^n\to\mathbb{R}gf(x)=f(x)2g_f(\boldsymbol{x})=f(\boldsymbol{x})^2。证明对任意 fn=1Fnf\in\bigcup_{n=1}^{\infty}\mathcal{F}_n,函数 gfg_f 都是凸函数。
    2. 给定正数 α\alphafFnf\in\mathcal{F}_n,定义 hf,α:RnRh_{f,\alpha}:\mathbb{R}^n\to\mathbb{R}hf,α(x)=f(x)αh_{f,\alpha}(\boldsymbol{x})=f(\boldsymbol{x})^\alpha。 求最小实数 α\alpha^*,使得对每个 αα\alpha\geqq\alpha^* 以及每个 fn=1Fnf\in\bigcup_{n=1}^{\infty}\mathcal{F}_nhf,αh_{f,\alpha} 都是凸函数;还须证明所求 α\alpha^* 确实最小。

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)λxL(x, \mu,\lambda) = -\theta (x) + \mu (\boldsymbol{1}^\top x - 1)-\lambda^\top x

Since the uniform feasible point has positive objective value, every optimizer is in R+n\mathbb R_+^n. The KKT conditions are therefore

KKT-conditions {θ(x)nxi+μλi=0,i=1,,n,1x=1,x0,λ0,λixi=0,i=1,,n.\text{KKT-conditions } \left\{ \begin{aligned} &-\frac{\theta(x)}{n x_i}+\mu-\lambda_i=0,\qquad i=1,\ldots,n,\\ &\boldsymbol{1}^\top x=1,\qquad x\succeq0,\\ &\lambda\succeq0,\qquad \lambda_i x_i=0,\quad i=1,\ldots,n. \end{aligned} \right.

(b)

xx^*, μ\mu^* and λ\lambda^* satisfy the KKT conditions for

x=[1n,1n,,1n],μ=1n,λ=0.x^* = \left[\frac 1n, \frac 1n, \ldots , \frac 1n\right]^\top,\qquad \mu^*=\frac1n,\qquad \lambda^*=\boldsymbol0.

Since θ-\theta is convex, the KKT conditions are sufficient; the stated uniqueness then gives this as the unique optimal solution.

(c)(c)

If iγi=0\sum_i\gamma_i=0, the claim is immediate. Otherwise, put xi=γi/jγjx_i=\gamma_i/\sum_j\gamma_j. Since θ(x)θ(x)=1/n\theta(x)\leq\theta(x^*)=1/n,

(i=1nγi)1n=(i=1nγi)θ(x)1ni=1nγi.\left(\prod_{i=1}^{n} \gamma_i\right)^{\frac 1n} =\left(\sum_{i=1}^{n}\gamma_i\right)\theta(x) \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

For a=f(x1)0a=f(x_1)\geq0 and b=f(x2)0b=f(x_2)\geq0,

θa2+(1θ)b2(θa+(1θ)b)2=θ(1θ)(ab)20.\theta a^2+(1-\theta)b^2-(\theta a+(1-\theta)b)^2 =\theta(1-\theta)(a-b)^2\geq0.

Therefore,

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 .

For 0<α<10<\alpha<1, let f(x)=xF1f(x)=|x|\in\mathcal F_1. Then

hf,α(1)=1>hf,α(0)+hf,α(2)2=2α1,h_{f,\alpha}(1)=1> \frac{h_{f,\alpha}(0)+h_{f,\alpha}(2)}2=2^{\alpha-1},

so hf,αh_{f,\alpha} is not convex. Hence α=1\alpha^*=1.