跳到主要内容

東京工業大学 情報理工学院 数理・計算科学系 2017年8月実施 午前 問4

Author

peter8rabit

Description

実数 θ\theta をパラメータとする次の線形計画問題 P(θ)\mathcal{P}(\theta) を考える:

P(θ) : {max : 6x1+(3+θ)x2+(4+θ)x3constraints :3x1+2x24x352x1+2x2+x324x1+2x2+3x35x1,x2,x30\mathcal{P}(\theta) \ : \ \begin{cases} &\max \ : \ 6x_1 + (-3 + \theta)x_2 + (4+\theta)x_3 \\ &\text{constraints} \ : \\ &-3x_1 + 2x_2 - 4x_3 \leq 5 \\ &2x_1 + 2x_2 + x_3 \leq 2 \\ &4x_1 + 2x_2 + 3x_3 \leq 5 \\ & x_1, x_2, x_3 \geq 0 \end{cases}

P(θ)\mathcal{P}(\theta) の最適解と最適値をそれぞれ (x1(θ),x2(θ),x3(θ))(x_1^*(\theta), x_2^*(\theta), x_3^*(\theta))z(θ)z^*(\theta) とする.以下の問に答えよ.

(1) θ=0\theta = 0 のときの最適解 (x1(0),x2(0),x3(0))(x_1^*(0), x_2^*(0), x_3^*(0))z(0)z^*(0) を求めよ.

(2) (x1(0),x2(0),x3(0))(x_1^*(0), x_2^*(0), x_3^*(0))P(θ)\mathcal{P}(\theta) の最適解でもあるような θ\theta の範囲を求めよ.

(3) (2) で求めた θ\theta の範囲において,z(θ)z^∗(\theta)θ\theta の関数として表せ.

Kai

(1)

初期辞書

z=6x13x2+4x3x4=5+3x12x2+4x3x5=22x12x2x3x6=54x12x23x3\begin{aligned} &z = 6x_1 - 3x_2 + 4x_3 \\ &x_4 = 5 + 3x_1 - 2x_2 + 4x_3 \\ &x_5 = 2 - 2x_1 -2x_2 -x_3 \\ &x_6 = 5 - 4x_1 -2x_2 -3x_3 \end{aligned}

最適辞書

z=77x2x5x6x1=122x232x5+12x6x3=1+2x2+2x5x6x4=212+72x552x6\begin{aligned} &z = 7 - 7x_2 -x_5 - x_6 \\ &x_1 = \frac{1}{2} - 2x_2 - \frac{3}{2}x_5 + \frac{1}{2}x_6 \\ &x_3 = 1 + 2x_2 + 2x_5 - x_6 \\ &x_4 = \frac{21}{2} + \frac{7}{2}x_5 - \frac{5}{2}x_6 \end{aligned}

よって、

z=7,(x1(0),x2(0),x3(0))=(12,0,1)z^* = 7, (x_1^*(0), x_2^*(0), x_3^*(0)) = (\frac{1}{2}, 0, 1)

(2)

b=(525),xB=(x1x3x4),xN=(x2x5x6),cB=(64+θ0),cN=(3+θ00)b = \begin{pmatrix} 5 \\ 2 \\ 5 \end{pmatrix}, x_B = \begin{pmatrix} x_1 \\ x_3 \\ x_4 \end{pmatrix}, x_N = \begin{pmatrix} x_2 \\ x_5 \\ x_6 \end{pmatrix}, c_B = \begin{pmatrix} 6 \\ 4 + \theta \\ 0 \end{pmatrix}, c_N = \begin{pmatrix} -3 + \theta \\ 0 \\ 0 \end{pmatrix}
B=(341210430),N=(200210201)B = \begin{pmatrix} -3 & -4 & 1 \\ 2 & 1 & 0 \\ 4 & 3 & 0 \end{pmatrix}, N = \begin{pmatrix} 2 & 0 & 0 \\ 2 & 1 & 0 \\ 2 & 0 & 1 \end{pmatrix}

とすると、最適辞書

z=cBTB1b+(cNTcBTB1N)xNxB=B1bB1NxN\begin{aligned} z &= c_B^T B^{-1}b + (c_N^T - c_B^T B^{-1} N) x_N \\ x_B &= B^{-1}b - B^{-1}N x_N \end{aligned}

題意より、B1b0B^{-1}b \geq 0 かつ cNTcBTB1N0c_N^T - c_B^TB^{-1}N \leq 0。 よって、1θ1/2-1 \leq \theta \leq 1/2

(3)

z(θ)=6x1(0)+(3+θ)x2(0)+(4+θ)x3(0)=7+θ\begin{aligned} z^*(\theta) &= 6x_1^*(0) + (-3 + \theta)x_2^*(0) + (4+\theta)x_3^*(0) \\ &= 7 + \theta \end{aligned}