非线性规划之分式规划
文章目录1. 分式规划例子2. 分式规划例题1. 分式规划例子在数学规划问题中,若目标函数为分式函数,且约束条件中的函数是线性的,则称线性规划分式规划,简称分式规划。分式规划通常可表示为如下形式:minpTx+αqTx+βs.t.{Ax≤b,x≥0.\min \frac{\bm{p}^\text{T}\bm{x} + \alpha}{\bm{q}^\text{T}\bm{x} + \beta}
1. 分式规划例子
在数学规划问题中,若目标函数为分式函数,且约束条件中的函数是线性的,则称线性规划分式规划,简称分式规划。
分式规划通常可表示为如下形式:
min
p
T
x
+
α
q
T
x
+
β
s.t.
{
A
x
≤
b
,
x
≥
0
.
\min \frac{\bm{p}^\text{T}\bm{x} + \alpha}{\bm{q}^\text{T}\bm{x} + \beta} \\ \text{s.t.} \begin{cases} \bm{Ax} \leq \bm{b}, \\ \bm{x} \geq \bm{0}. \\ \end{cases}
minqTx+βpTx+αs.t.{Ax≤b,x≥0.
其中,
α
\alpha
α,
β
\beta
β 为常数;
p
\bm{p}
p,
q
\bm{q}
q 为
n
n
n 维列向量;
A
\bm{A}
A 为
m
×
n
m \times n
m×n 阶矩阵。
这一类问题有类似于线性规划问题的极好的性质。
- 若分式规划问题存在最优解,则最优解可在可行域顶点上达到。
- 任一局部极小值即全局极小值。
由 Charnes 和 Cooper 于 1962 年提出的用单纯刑法求解分式规划问题的方法:
设集合上
S
=
{
x
∈
R
n
∣
A
x
≤
b
,
b
≥
0
}
S = \{ x \in \mathbb{R}^n \mid \bm{Ax} \leq \bm{b}, \bm{b} \geq \bm{0} \}
S={x∈Rn∣Ax≤b,b≥0} 是有界闭集,且对
∀
x
∈
S
\forall \bm{x} \in S
∀x∈S,有
q
T
x
+
β
>
0
\bm{q}^\text{T}\bm{x} + \beta > 0
qTx+β>0。引入新变量
z
z
z,令
z
=
1
q
T
x
+
β
z = \frac{1}{\bm{q}^\text{T}\bm{x} + \beta}
z=qTx+β1,
y
=
z
x
\bm{y} = z \bm{x}
y=zx,则以上模型可转化为线性规划模型:
min
p
T
y
+
α
z
s.t.
{
A
y
−
b
z
≤
0
,
q
T
y
+
β
z
=
1
,
y
≥
0
,
z
≥
0.
\min \bm{p}^\text{T}\bm{y} + \alpha z \\ \text{s.t.} \begin{cases} \bm{Ay} - \bm{b}z \leq 0, \\ \bm{q}^\text{T}\bm{y} + \beta z = 1, \\ \bm{y} \geq 0, z \geq 0. \\ \end{cases}
minpTy+αzs.t.⎩⎪⎨⎪⎧Ay−bz≤0,qTy+βz=1,y≥0,z≥0.
至此,可用单纯形法来求解此规划,并最终得到原分式规划的最优解。
原分式
p
T
x
+
α
q
T
x
+
β
\frac{\bm{p}^\text{T}\bm{x} + \alpha}{\bm{q}^\text{T}\bm{x} + \beta}
qTx+βpTx+α 转化为
p
T
x
+
α
q
T
x
+
β
=
p
T
x
q
T
x
+
β
+
α
q
T
x
+
β
=
p
T
x
z
+
α
z
=
p
T
y
+
α
z
.
\begin{aligned} \frac{\bm{p}^\text{T}\bm{x} + \alpha}{\bm{q}^\text{T}\bm{x} + \beta} &= \frac{\bm{p}^\text{T}\bm{x}}{\bm{q}^\text{T}\bm{x} + \beta} + \frac{\alpha}{\bm{q}^\text{T}\bm{x} + \beta} \\ & = \bm{p}^\text{T} \bm{x} z + \alpha z \\ & = \bm{p}^\text{T} \bm{y} + \alpha z. \\ \end{aligned}
qTx+βpTx+α=qTx+βpTx+qTx+βα=pTxz+αz=pTy+αz.
原约束
A
x
≤
b
\bm{Ax} \leq \bm{b}
Ax≤b 转化为:
A
x
≤
b
A
y
z
≤
b
A
y
≤
b
z
A
y
−
b
z
≤
0
\begin{aligned} \bm{Ax} & \leq \bm{b} \\ \bm{A}\frac{\bm{y}}{z} & \leq \bm{b} \\ \bm{Ay} &\leq \bm{b} z \\ \bm{Ay} - \bm{b}z & \leq 0 \\ \end{aligned}
AxAzyAyAy−bz≤b≤b≤bz≤0
z
=
1
q
T
x
+
β
z = \frac{1}{\bm{q}^\text{T}\bm{x} + \beta}
z=qTx+β1 转化为:
z
(
q
T
x
+
β
)
=
1
q
T
z
x
+
β
z
=
1
q
T
y
+
β
z
=
1
\begin{aligned} z (\bm{q}^\text{T}\bm{x} + \beta) & = 1 \\ \bm{q}^\text{T} z \bm{x} + \beta z & = 1 \\ \bm{q}^\text{T} \bm{y} + \beta z & = 1 \\ \end{aligned}
z(qTx+β)qTzx+βzqTy+βz=1=1=1
2. 分式规划例题
例 求解下列分式规划:
min
−
2
x
1
+
x
2
+
2
x
1
+
3
x
2
+
4
s.t.
{
−
x
1
+
x
2
≤
4
,
2
x
1
+
x
2
≤
14
,
x
2
≤
6
,
x
1
≥
0
,
x
2
≥
0.
\min \frac{-2x_1 + x_2 + 2}{x_1 + 3x_2 + 4} \\ \text{s.t.} \begin{cases} -x_1 + x_2 \leq 4, \\ 2x_1 + x_2 \leq 14, \\ x_2 \leq 6, \\ x_1 \geq 0, x_2 \geq 0. \\ \end{cases}
minx1+3x2+4−2x1+x2+2s.t.⎩⎪⎪⎪⎨⎪⎪⎪⎧−x1+x2≤4,2x1+x2≤14,x2≤6,x1≥0,x2≥0.
解 令 z = 1 x 1 + 3 x 2 + 4 z = \frac{1}{x_1 + 3x_2 + 4} z=x1+3x2+41, y = z x \bm{y} = z\bm{x} y=zx,则原分式规划问题可转化为如下等价的线性规划模型:
A = [ − 1 1 2 1 0 1 ] , b = [ 4 14 6 ] , y = [ y 1 y 2 ] , p = [ − 2 1 ] , q = [ 1 3 ] , α = 2 , β = 4. \bm{A} = \begin{bmatrix} -1 & 1 \\ 2 & 1 \\ 0 & 1 \\ \end{bmatrix}, \bm{b} = \begin{bmatrix} 4 \\ 14 \\ 6 \\ \end{bmatrix}, \bm{y} = \begin{bmatrix} y_1 \\ y_2 \\ \end{bmatrix}, \bm{p} = \begin{bmatrix} -2 \\ 1 \\ \end{bmatrix}, \bm{q} = \begin{bmatrix} 1 \\ 3 \\ \end{bmatrix}, \\ \alpha = 2, \beta = 4. A=⎣⎡−120111⎦⎤,b=⎣⎡4146⎦⎤,y=[y1y2],p=[−21],q=[13],α=2,β=4.
p T y + α z = − 2 y 1 + y 2 + 2 z A y − b z ≤ 0 → − y 1 + y 2 − 4 z ≤ 0 2 y 1 + y 2 − 14 z ≤ 0 y 2 − 6 z ≤ 0 q T y + β z = 1 → y 1 + 3 y 2 + 4 z = 1 \bm{p}^\text{T} \bm{y} + \alpha z = -2 y_1 + y_2 + 2z \\ \bm{Ay} - \bm{b} z \leq 0 \rightarrow \\ -y_1 + y_2 - 4 z \leq 0 \\ 2y_1 + y_2 - 14 z \leq 0 \\ y_2 - 6 z \leq 0 \\ \bm{q}^\text{T} \bm{y} + \beta z = 1 \rightarrow \\ y_1 + 3y_2 + 4z = 1 \\ pTy+αz=−2y1+y2+2zAy−bz≤0→−y1+y2−4z≤02y1+y2−14z≤0y2−6z≤0qTy+βz=1→y1+3y2+4z=1
min − 2 y 1 + y 2 + 2 z s.t. { − y 1 + y 2 − 4 z ≤ 0 2 y 1 + y 2 − 14 z ≤ 0 y 2 − 6 z ≤ 0 y 1 + 3 y 2 + 4 z = 1 y 1 ≥ 0 , y 2 ≥ 0 , z ≥ 0. \min -2y_1 + y_2 + 2z \\ \text{s.t.} \begin{cases} -y_1 + y_2 - 4z \leq 0 \\ 2y_1 + y_2- 14z \leq 0 \\ y_2 - 6z \leq 0 \\ y_1 + 3y_2 + 4z = 1 \\ y_1 \geq 0, y_2 \geq 0, z \geq 0. \\ \end{cases} min−2y1+y2+2zs.t.⎩⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎧−y1+y2−4z≤02y1+y2−14z≤0y2−6z≤0y1+3y2+4z=1y1≥0,y2≥0,z≥0.
首先化成标准型:
min
w
=
−
2
y
1
+
y
2
+
2
z
max
(
−
w
)
=
2
y
1
−
y
2
−
2
z
−
M
y
6
s.t.
{
−
y
1
+
y
2
−
4
z
+
y
3
=
0
,
2
y
1
+
y
2
−
14
z
+
y
4
=
0
,
y
2
−
6
z
+
y
5
=
0
,
y
1
+
3
y
2
+
4
z
+
y
6
=
1
,
y
j
≥
0
,
z
≥
0.
\min w = -2y_1 + y_2 + 2z \\ \max (-w) = 2y_1 - y_2 - 2z - M y_6 \\ \text{s.t.} \begin{cases} -y_1 + y_2 - 4z + y_3 = 0, \\ 2y_1 + y_2 - 14z + y_4 = 0, \\ y_2 - 6z + y_5 = 0, \\ y_1 + 3y_2 + 4z + y_6 = 1, \\ y_j \geq 0, z \geq 0. \\ \end{cases}
minw=−2y1+y2+2zmax(−w)=2y1−y2−2z−My6s.t.⎩⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎧−y1+y2−4z+y3=0,2y1+y2−14z+y4=0,y2−6z+y5=0,y1+3y2+4z+y6=1,yj≥0,z≥0.
初始单纯形表:
. | y 1 y_1 y1 | y 2 y_2 y2 | y 3 y_3 y3 | y 4 y_4 y4 | y 5 y_5 y5 | y 6 y_6 y6 | z | 右边 |
---|---|---|---|---|---|---|---|---|
. | -2-M | 1-3M | 0 | 0 | 0 | 0 | 2-4M | -M |
y 3 y_3 y3 | -1 | 1 | 1 | 0 | 0 | 0 | -4 | 0 |
y 4 y_4 y4 | 2 | 1 | 0 | 1 | 0 | 0 | -14 | 0 |
y 5 y_5 y5 | 0 | 1 | 0 | 0 | 1 | 0 | -6 | 0 |
y 6 y_6 y6 | 1 | 3 | 0 | 0 | 0 | 1 | 4 | 1 |
最优单纯形表:
. | y 1 y_1 y1 | y 2 y_2 y2 | y 3 y_3 y3 | y 4 y_4 y4 | y 5 y_5 y5 | y 6 y_6 y6 | z | 右边 |
---|---|---|---|---|---|---|---|---|
. | 0 | 52 11 \frac{52}{11} 1152 | 0 | 5 11 \frac{5}{11} 115 | 0 | M + 12 11 M+\frac{12}{11} M+1112 | 0 | 12 11 \frac{12}{11} 1112 |
y 3 y_3 y3 | 0 | 4 | 1 | 0 | 0 | 1 | 0 | 1 |
y 1 y_1 y1 | 1 | 23 11 \frac{23}{11} 1123 | 0 | 2 11 \frac{2}{11} 112 | 0 | 7 11 \frac{7}{11} 117 | 0 | 7 11 \frac{7}{11} 117 |
y 5 y_5 y5 | 0 | 15 11 \frac{15}{11} 1115 | 0 | − 3 11 -\frac{3}{11} −113 | 1 | 6 11 \frac{6}{11} 116 | 0 | 6 11 \frac{6}{11} 116 |
z z z | 0 | 5 22 \frac{5}{22} 225 | 0 | 1 22 \frac{1}{22} 221 | 0 | 1 11 \frac{1}{11} 111 | 1 | 1 11 \frac{1}{11} 111 |
故:
y
1
=
7
11
y_1 = \frac{7}{11}
y1=117,
y
2
=
0
y_2 = 0
y2=0,
z
=
1
11
z = \frac{1}{11}
z=111 是以上线性规划模型的最优解。
故原分式规划的最优解为
x
1
=
y
1
z
=
7
x_1 = \frac{y_1}{z} = 7
x1=zy1=7,
x
2
=
y
2
z
=
0
x_2 = \frac{y_2}{z} = 0
x2=zy2=0。
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
所有评论(0)