京都大学 情報学研究科 数理工学専攻 2010年8月実施 オペレーションズ・リサーチ
Author
思齐塾, 祭音Myyura
Description
R+n={x∈Rn∣xi≥0 (i=1,…,n)},R++n={x∈Rn∣xi>0 (i=1,…,n)}.
関数 ψ:R+n→R を
ψ(x)=i=1∑nxilnxi
と定義する。ただし、 ln は自然対数を表し、 0ln0=0 とする。さらに、関数 Bψ:R+n×R++n→R を次のように定義する。
Bψ(x,y)=ψ(x)−ψ(y)−∇ψ(y)T(x−y)
ただし、 T は転置記号である。
次に、パラメータ t∈R を含む非線形計画問題 P(t) を考える。
P(t)minimizetcTx+Bψ(x,y)
subject toi=1∑nxi=1
xi≥0(i=1,...,n)
ここで、決定変数は x であり、 c∈Rn と y∈R++n は定数ベクトルである。問題 P(t) には唯一の解 x(t) が存在し、 xi(t)>0(i=1,...,n) が成り立つことが知られている。以下の問 (i)-(iv) に答えよ。
(i) 任意の x,y∈R++n に対して、 Bψ(x,y)≥0 となることを示せ。
(ii) 問題 P(t) のカルーシュ・キューン・タッカー条件 (Karush-Kuhn-Tucker 条件) を書け。
(iii) x(t) を求めよ。
(iv) ベクトル c の成分 c1,c2,...,cn に対して c1>c2>...>cn が成り立つとする。
t→∞limx(t)=(0,...,0,1)T
となることを示せ。
题目描述
记
R+n={x∈Rn∣xi≥0 (i=1,…,n)},R++n={x∈Rn∣xi>0 (i=1,…,n)}.
在 R+n 上定义
ψ(x)=i=1∑nxilnxi,
其中 ln 为自然对数,并约定 0ln0=0。对 x∈R+n、y∈R++n,定义
Bψ(x,y)=ψ(x)−ψ(y)−∇ψ(y)T(x−y),
其中上标 T 表示转置。
给定参数 t∈R,考虑以 x 为决策变量的非线性规划问题
P(t):xmins.t.tcTx+Bψ(x,y)i=1∑nxi=1,xi≥0(i=1,…,n),
其中 c∈Rn、y∈R++n 为给定常向量。已知 P(t) 存在唯一解 x(t),且其每个分量都严格为正,即 xi(t)>0 (i=1,…,n)。
完成以下各问:
-
证明对任意 x,y∈R++n,都有
Bψ(x,y)≥0.
-
写出问题 P(t) 的 Karush–Kuhn–Tucker(KKT)条件。
-
求出唯一最优解 x(t)。
-
若 c1>c2>⋯>cn,证明
t→∞limx(t)=(0,…,0,1)T.
Kai
(i) ∇ψ(y)=(logy1+1,…,logyn+1)T なので
Bψ(x,y)=i=1∑n[xilogyixi−xi+yi].
r>0 に対して ϕ(r)=rlogr−r+1 とおくと
ϕ′(r)=logr,ϕ′′(r)=r1>0.
従って ϕ は r=1 で唯一の最小値 ϕ(1)=0 を取る。各項は
xilogyixi−xi+yi=yiϕ(yixi)≥0
であるから
Bψ(x,y)≥0.
この証明は、一般の x,y∈R++n に対して成り立ち、成分和が1であることを仮定しない。
(ii) Lagrangian:
L(x,λ,μ)=tcTx+ψ(x)−ψ(y)−∇ψ(y)T(x−y)−λ(i=1∑nxi−1)−i=1∑nμixi
KKT Conditions:
∂xi∂L=tci+lnxi+1−(lnyi+1)−λ−μi=0
i=1∑nxi=1
xi≥0,μi≥0,μixi=0
⟹tci+lnxi−lnyi−λ−μi=0
(iii) From the KKT conditions:
lnxi=lnyi−tci+λ+μi
xi=yie−tci+λ+μi
If xi>0 , then μi=0 .
xi=yie−tci+λ
i=1∑nxi=i=1∑nyie−tci+λ=1
eλ=∑i=1nyie−tci1
xi=∑j=1nyje−tcjyie−tci
(iv) Since c1>c2>...>cn ,
t→∞limxi(t)=t→∞lim∑j=1nyje−tcjyie−tci
=t→∞lim∑j=1nyje−t(cj−cn)yie−t(ci−cn)
For i<n,ci−cn>0 , so limt→∞e−t(ci−cn)=0 .
For i=n,cn−cn=0 , so limt→∞e−t(cn−cn)=1 .
t→∞limxi(t)=yn0=0
t→∞limxn(t)=∑j=1nyje−t(cj−cn)yn=ynyn=1
Thus, limt→∞x(t)=(0,...,0,1)T .