多目标规划方法讲义
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
目标规划的图解法
一、目标规划概述
目标规划是在线性规划的基础上,为适应经济管理 中多目标决策的需要而逐步发展起来的一个分支。
(一)、目标规划与线性规划的比较
1、线性规划只讨论一个线性目标函数在一组线性约束 条件下的极值问题;而目标规划是多个目标决策,可求 得更切合实际的解。
2、线性规划求最优解;目标规划是找到一个满意解。
对于线性多目标规划问题,可以进一步用矩阵表示:
max(min)Z CX s.t. AX b
式中:
X 为n 维决策变量向量; C 为k×n 矩阵,即目标函数系数矩阵; A 为m×n 矩阵,即约束方程系数矩阵; b 为m 维的向量,即约束向量。
多目标规划的非劣解
max(min)Z F ( X )
显然,这是一个多目标规划问题,用线性规划方法 很难找到最优解。
1、目标值和偏差变量
目标规划通过引入目标值和偏差变量,可以将目标函 数转化为目标约束。
目标值:是指预先给定的某个目标的一个期望值。
实现值或决策值:是指当决策变量xj 选定以后,目标
函数的对应值。 偏差变量(事先无法确定的未知数):是指实现值和
一 多目标规划及其非劣解
多目标规划模型
(一)任何多目标规划问题,都由两个基本部分组成: (1)两个以上的目标函数; (2)若干个约束条件。
(二)对于多目标规划问题,可以将其数学模型一般地描 写为如下形式:
max(min)
f1
(
X
)
Z F ( X ) max(min) f2 ( X )
max(min) fk ( X )
1( X )
g1
s.t.
ห้องสมุดไป่ตู้
(
X
)
2(X
)
G
g2
m ( X )
gm
式中: X [ x1, x2 ,, xn ]T 为决策变量向量。
缩写形式: max(min)Z F ( X ) s.t. ( X ) G
(1) (2)
有n个决策变量,k个目标函数, m个约束方程, 则:
Z=F(X) 是k维函数向量, (X)是m维函数向量; G是m维常数向量;
max(min)Z f1( x1, x2,, xn )
i ( x1, x2,, xn ) gi (i 1,2,, m)
f min j
fj
f
max j
(
j
2,3,, k)
方法四 目标达到法 首先将多目标规划模型化为如下标准形式:
f1( X )
min
F
(
x
)
min
f2
(X
)
fk
(
X
)
若规定3600的钢材必须用完,原式9 x1 +4 x2 ≤3600
则变为
9 x1
4x2
d
4
d
4
3600
d
4
,
d
4
0
3、达成函数(即目标规划中的目标函数)
达成函数是一个使总偏差量为最小的目标函数,记为
minZ = f(d+、d-)。
一般说来,有以下三种情况,但只能出现其中之一: ⑴.要求恰好达到规定的目标值,即正、负偏差变量要
尽可能小,则minZ = f(d++ d-)。
⑵.要求不超过目标值,即允许达不到目标值,也就是 正偏差变量尽可能小,则minZ = f(d+)。
⑶.要求超过目标值,即超过量不限,但不低于目标值, 也就是负偏差变量尽可能小,则minZ = f(d-)。
对于由绝对约束转化而来的目标函数,也照上述处理即 可。
目标值之间的差异,记为 d 。 正偏差变量:表示实现值超过目标值的部分,记为 d
+。 负偏差变量:表示实现值未达到目标值的部分,记为
d-。
在一次决策中,实现值不可能既超过目标值又未达 到目标值,故有 d+× d- =0,并规定d+≥0, d-≥0
当完成或超额完成规定的指标则表示:d+≥0, d-=0 当未完成规定的指标则表示: d+=0, d-≥0 当恰好完成指标时则表示: d+=0, d-=0
( X ) G
式中,i 是与第i个目标函数相关的权重; A是由 i (i=1,2,…,k )组成的m×m对角矩阵。
方法三 约束模型(极大极小法)
理论依据 :若规划问题的某一目标可以给出一个可供选 择的范围,则该目标就可以作为约束条件而被排除出目标 组,进入约束条件组中。
假如,除第一个目标外,其余目标都可以提出一个可供选 择的范围,则该多目标规划问题就可以转化为单目标规划 问题:
s.t. ( X ) G
多目标规划问题的求解不能只追求一个目标的最优化 (最大或最小),而不顾其它目标。 对于上述多目标规划问题,求解就意味着需要做出如下 的复合选择: ▲ 每一个目标函数取什么值,原问题可以得到最满意 的解决? ▲ 每一个决策变量取什么值,原问题可以得到最满意 的解决 ?
非劣解可以用图1说明。
1( X ) 0
(
X
)
2
(X
)
0
m ( X ) 0
在求解之前,先设计与目标函数相应的一组目标值理想
化的期望目标 fi* ( i=1,2,…,k ) ,
每一个目标对应的权重系数为 i* ( i=1,2,…,k ) , 再设 为一松弛因子。
那么,多目标规划问题就转化为:
f1( X )
当目标函数处于冲突状态时,就不会存在使所有目 标函数同时达到最大或最小值的最优解,于是我们只能 寻求非劣解(又称非支配解或帕累托解)。
二 多目标规划求解
为了求得多目标规划问题的非劣解,常常需要将 多目标规划问题转化为单目标规划问题去处理。实现 这种转化,有如下几种建模方法。
✓ 效用最优化模型 ✓ 理想点模型 ✓ 约束模型 ✓ 目标达到法 ✓ 目标规划模型
min
F
(
x)
min
f2
(X
)
fk
(
X
)
1( X ) 0
(
X
)
2
(X
)
0
m ( X ) 0
min X ,
i(X) 0
(i 1,2,, m)
fi ( X ) i fi* , (i 1,2,, k)
目标规划
(Goal programming)
目标规划概述 目标规划的数学模型
∴ d+× d- =0 成立。
2、目标约束和绝对约束
引入了目标值和正、负偏差变量后,就对某一问题 有了新的限制,既目标约束。
目标约束即可对原目标函数起作用,也可对原约束起 作用。目标约束是目标规划中特有的,是软约束。
绝对约束(系统约束)是指必须严格满足的等式或 不等式约束。如线性规划中的所有约束条件都是绝对 约束,否则无可行解。所以,绝对约束是硬约束。
4 x2
d
4
d
4
3600
图1 多目标规划的劣解与非劣解
在图1中,max(f1, f2) .就 方案①和②来说,①的 f2 目标值比②大,但其目 标值 f1 比②小,因此无 法确定这两个方案的优 与劣。
在各个方案之间, 显然:④比①好,⑤比 ④好, ⑥比②好, ⑦比 ③好……。
而对于方案⑤、 ⑥、⑦之间则无法确 定优劣,而且又没有 比它们更好的其他方 案,所以它们就被称 为多目标规划问题的 非劣解或有效解, 其余方案都称为劣解。 所有非劣解构成的集 合称为非劣解集。
方法一 效用最优化模型(线性加权法)
思想:规划问题的各个目标函数可以通过一定的方式 进行求和运算。这种方法将一系列的目标函数与效用 函数建立相关关系,各目标之间通过效用函数协调, 使多目标规划问题转化为传统的单目标规划问题:
max Z ( X )
(1)
s.t. ( X ) G (2)
是与各目标函数相关的效用函数的和函数。
(二)、目标规划的基本概念
例一、某厂计划在下一个生产周期内生产甲、乙两 种产品,已知资料如表所示。试制定生产计划,使获 得的利润最大?同时,根据市场预测,甲的销路不是 太好,应尽可能少生产;乙的销路较好,可以扩大生 产。试建立此问题的数学模型。
单位 产品 甲
资源 消耗
钢材
9
煤炭
4
设备台时 3
单件利润 70
4、优先因子(优先等级)与优先权系数
优先因子Pk 是将决策目标按其重要程度排序并表示 出来。P1>>P2>>…>>Pk>>Pk+1>>…>>PK , k=1.2…K。
权系数ωk 区别具有相同优先因子的两个目标的差别, 决策者可视具体情况而定。
5、满意解(具有层次意义的解)
对于这种解来说,前面的目标可以保证实现或部分 实现,而后面的目标就不一定能保证实现或部分实现, 有些可能就不能实现。
例二、若在例一中提出下列要求: 1、完成或超额完成利润指标 50000元; 2、产品甲不超过 200件,产品乙不低于 250件; 3、现有钢材 3600吨必须用完。
试建立目标规划模型。
分析:题目有三个目标层次,包含四个目标值。
第一目标: P1d1
第二目标:有两个要求即甲
d 2,乙
d3,但两
个具有相同的优先因子,因此需要确定权系数。本题
在用效用函数作为规划目标时,需要确定一组权值 i
来反映原问题中各目标函数在总体目标中的权重,即:
k
max i i
i1
i ( x1, x2,xn ) gi (i 1,2,, m)
式中, i 应满足:
k
i 1
i 1
向量形式: max T
s.t. ( X ) G
方法二 理想点模型(罚款模型)
例如:在例一中,规定Z1 的目标值为 50000,正、负
偏差为d+、d- ,则目标函数可以转换为目标约束,既
70 x1 + 120 x2+ d1 d1=50000,
同样,若规定 Z2=200, Z3=250 则有
x1
d
2
d
2
200
x2 d3 d3 250
d
j
,
d
j
0
( j 1.2.3)
求解多目标规划的方法大体上有以下几种:
一种是化多为少的方法 , 即把多目标化为比较容易求解 的单目标或双目标,如主要目标法、线性加权法、理想点 法等;
另一种叫分层序列法,即把目标按其重要性给出一个序列, 每次都在前一目标最优解集内求下一个目标最优解,直到 求出共同的最优解。
对多目标的线性规划除以上方法外还可以适当修正单纯形 法来求解;还有一种称为层次分析法,是由美国运筹学家 沙旦于70年代提出的,这是一种定性与定量相结合的多目 标决策与分析方法,对于目标结构复杂且缺乏必要的数据 的情况更为实用。
思想: 规划决策者对每一个目标函数都能提出所期望的值 (或称满意值); 通过比较实际值 fi 与期望值 fi* 之间的偏差来选择问题的 解,其数学表达式如下:
k
min Z i ( fi fi )2
i ( x1, x2,i1 , xn ) gi (i 1,2,, m)
或写成矩阵形式: min Z (F F )T A(F F )
可用单件利润比作为权系数即 70 :120,化简为7:12。
P2
(7
d
2
12d
3
)
第三目标:
P3
(d
4
d
4
)
目标规划模型为:
min Z
P1d1
P2
(7d
2
12d
3
)
P3
(d
4
d
4
)
70
x1 x1
120
x2
d1
d
2
d1
d
2
50000 200
9x1
x2
d
3
d
3
250
3、线性规划中的约束条件是同等重要的,是硬约束; 而目标规划中有轻重缓急和主次之分,即有优先权。
4、线性规划的最优解是绝对意义下的最优,但需花 去大量的人力、物力、财力才能得到;实际过程中, 只要求得满意解,就能满足需要(或更能满足需要)。
目前,已经在经济计划、生产管理、经营管理、市场 分析、财务管理等方面得到了广泛的应用。
多目标规划方法
多目标规划是数学规划的一个分支。 研究多于一个的目标函数在给定区域上的最优化。又称多 目标最优化。通常记为
MOP(multi-objective programming)。 在很多实际问题中,例如经济、管理、军事、科学和工程 设计等领域,衡量一个方案的好坏往往难以用一个指标来 判断,而需要用多个目标来比较,而这些目标有时不甚协 调,甚至是矛盾的。因此有许多学者致力于这方面的研究。 1896年法国经济学家 V. 帕雷托最早研究不可比较目标的优 化问题,之后,J.冯·诺伊曼、H.W.库恩、A.W.塔克、A.M. 日夫里翁等数学家做了深入的探讨,但是尚未有一个完全 令人满意的定义。
乙 资源限制
4
3600
5
2000
10 3000
120
设:甲产品 x1 ,乙产品 x2
一般有:
同时:
maxZ=70 x1 + 120 x2 9 x1 +4 x2 ≤3600 4 x1 +5 x2 ≤ 2000 3 x1 +10 x2 ≤3000 x1 , x2 ≥0
maxZ1=70 x1 + 120x2 minZ2= x1 maxZ3= x2 9 x1 +4 x2 ≤3600 4 x1 +5 x2 ≤ 2000 3 x1 +10 x2 ≤3000 x1 , x2 ≥0
一、目标规划概述
目标规划是在线性规划的基础上,为适应经济管理 中多目标决策的需要而逐步发展起来的一个分支。
(一)、目标规划与线性规划的比较
1、线性规划只讨论一个线性目标函数在一组线性约束 条件下的极值问题;而目标规划是多个目标决策,可求 得更切合实际的解。
2、线性规划求最优解;目标规划是找到一个满意解。
对于线性多目标规划问题,可以进一步用矩阵表示:
max(min)Z CX s.t. AX b
式中:
X 为n 维决策变量向量; C 为k×n 矩阵,即目标函数系数矩阵; A 为m×n 矩阵,即约束方程系数矩阵; b 为m 维的向量,即约束向量。
多目标规划的非劣解
max(min)Z F ( X )
显然,这是一个多目标规划问题,用线性规划方法 很难找到最优解。
1、目标值和偏差变量
目标规划通过引入目标值和偏差变量,可以将目标函 数转化为目标约束。
目标值:是指预先给定的某个目标的一个期望值。
实现值或决策值:是指当决策变量xj 选定以后,目标
函数的对应值。 偏差变量(事先无法确定的未知数):是指实现值和
一 多目标规划及其非劣解
多目标规划模型
(一)任何多目标规划问题,都由两个基本部分组成: (1)两个以上的目标函数; (2)若干个约束条件。
(二)对于多目标规划问题,可以将其数学模型一般地描 写为如下形式:
max(min)
f1
(
X
)
Z F ( X ) max(min) f2 ( X )
max(min) fk ( X )
1( X )
g1
s.t.
ห้องสมุดไป่ตู้
(
X
)
2(X
)
G
g2
m ( X )
gm
式中: X [ x1, x2 ,, xn ]T 为决策变量向量。
缩写形式: max(min)Z F ( X ) s.t. ( X ) G
(1) (2)
有n个决策变量,k个目标函数, m个约束方程, 则:
Z=F(X) 是k维函数向量, (X)是m维函数向量; G是m维常数向量;
max(min)Z f1( x1, x2,, xn )
i ( x1, x2,, xn ) gi (i 1,2,, m)
f min j
fj
f
max j
(
j
2,3,, k)
方法四 目标达到法 首先将多目标规划模型化为如下标准形式:
f1( X )
min
F
(
x
)
min
f2
(X
)
fk
(
X
)
若规定3600的钢材必须用完,原式9 x1 +4 x2 ≤3600
则变为
9 x1
4x2
d
4
d
4
3600
d
4
,
d
4
0
3、达成函数(即目标规划中的目标函数)
达成函数是一个使总偏差量为最小的目标函数,记为
minZ = f(d+、d-)。
一般说来,有以下三种情况,但只能出现其中之一: ⑴.要求恰好达到规定的目标值,即正、负偏差变量要
尽可能小,则minZ = f(d++ d-)。
⑵.要求不超过目标值,即允许达不到目标值,也就是 正偏差变量尽可能小,则minZ = f(d+)。
⑶.要求超过目标值,即超过量不限,但不低于目标值, 也就是负偏差变量尽可能小,则minZ = f(d-)。
对于由绝对约束转化而来的目标函数,也照上述处理即 可。
目标值之间的差异,记为 d 。 正偏差变量:表示实现值超过目标值的部分,记为 d
+。 负偏差变量:表示实现值未达到目标值的部分,记为
d-。
在一次决策中,实现值不可能既超过目标值又未达 到目标值,故有 d+× d- =0,并规定d+≥0, d-≥0
当完成或超额完成规定的指标则表示:d+≥0, d-=0 当未完成规定的指标则表示: d+=0, d-≥0 当恰好完成指标时则表示: d+=0, d-=0
( X ) G
式中,i 是与第i个目标函数相关的权重; A是由 i (i=1,2,…,k )组成的m×m对角矩阵。
方法三 约束模型(极大极小法)
理论依据 :若规划问题的某一目标可以给出一个可供选 择的范围,则该目标就可以作为约束条件而被排除出目标 组,进入约束条件组中。
假如,除第一个目标外,其余目标都可以提出一个可供选 择的范围,则该多目标规划问题就可以转化为单目标规划 问题:
s.t. ( X ) G
多目标规划问题的求解不能只追求一个目标的最优化 (最大或最小),而不顾其它目标。 对于上述多目标规划问题,求解就意味着需要做出如下 的复合选择: ▲ 每一个目标函数取什么值,原问题可以得到最满意 的解决? ▲ 每一个决策变量取什么值,原问题可以得到最满意 的解决 ?
非劣解可以用图1说明。
1( X ) 0
(
X
)
2
(X
)
0
m ( X ) 0
在求解之前,先设计与目标函数相应的一组目标值理想
化的期望目标 fi* ( i=1,2,…,k ) ,
每一个目标对应的权重系数为 i* ( i=1,2,…,k ) , 再设 为一松弛因子。
那么,多目标规划问题就转化为:
f1( X )
当目标函数处于冲突状态时,就不会存在使所有目 标函数同时达到最大或最小值的最优解,于是我们只能 寻求非劣解(又称非支配解或帕累托解)。
二 多目标规划求解
为了求得多目标规划问题的非劣解,常常需要将 多目标规划问题转化为单目标规划问题去处理。实现 这种转化,有如下几种建模方法。
✓ 效用最优化模型 ✓ 理想点模型 ✓ 约束模型 ✓ 目标达到法 ✓ 目标规划模型
min
F
(
x)
min
f2
(X
)
fk
(
X
)
1( X ) 0
(
X
)
2
(X
)
0
m ( X ) 0
min X ,
i(X) 0
(i 1,2,, m)
fi ( X ) i fi* , (i 1,2,, k)
目标规划
(Goal programming)
目标规划概述 目标规划的数学模型
∴ d+× d- =0 成立。
2、目标约束和绝对约束
引入了目标值和正、负偏差变量后,就对某一问题 有了新的限制,既目标约束。
目标约束即可对原目标函数起作用,也可对原约束起 作用。目标约束是目标规划中特有的,是软约束。
绝对约束(系统约束)是指必须严格满足的等式或 不等式约束。如线性规划中的所有约束条件都是绝对 约束,否则无可行解。所以,绝对约束是硬约束。
4 x2
d
4
d
4
3600
图1 多目标规划的劣解与非劣解
在图1中,max(f1, f2) .就 方案①和②来说,①的 f2 目标值比②大,但其目 标值 f1 比②小,因此无 法确定这两个方案的优 与劣。
在各个方案之间, 显然:④比①好,⑤比 ④好, ⑥比②好, ⑦比 ③好……。
而对于方案⑤、 ⑥、⑦之间则无法确 定优劣,而且又没有 比它们更好的其他方 案,所以它们就被称 为多目标规划问题的 非劣解或有效解, 其余方案都称为劣解。 所有非劣解构成的集 合称为非劣解集。
方法一 效用最优化模型(线性加权法)
思想:规划问题的各个目标函数可以通过一定的方式 进行求和运算。这种方法将一系列的目标函数与效用 函数建立相关关系,各目标之间通过效用函数协调, 使多目标规划问题转化为传统的单目标规划问题:
max Z ( X )
(1)
s.t. ( X ) G (2)
是与各目标函数相关的效用函数的和函数。
(二)、目标规划的基本概念
例一、某厂计划在下一个生产周期内生产甲、乙两 种产品,已知资料如表所示。试制定生产计划,使获 得的利润最大?同时,根据市场预测,甲的销路不是 太好,应尽可能少生产;乙的销路较好,可以扩大生 产。试建立此问题的数学模型。
单位 产品 甲
资源 消耗
钢材
9
煤炭
4
设备台时 3
单件利润 70
4、优先因子(优先等级)与优先权系数
优先因子Pk 是将决策目标按其重要程度排序并表示 出来。P1>>P2>>…>>Pk>>Pk+1>>…>>PK , k=1.2…K。
权系数ωk 区别具有相同优先因子的两个目标的差别, 决策者可视具体情况而定。
5、满意解(具有层次意义的解)
对于这种解来说,前面的目标可以保证实现或部分 实现,而后面的目标就不一定能保证实现或部分实现, 有些可能就不能实现。
例二、若在例一中提出下列要求: 1、完成或超额完成利润指标 50000元; 2、产品甲不超过 200件,产品乙不低于 250件; 3、现有钢材 3600吨必须用完。
试建立目标规划模型。
分析:题目有三个目标层次,包含四个目标值。
第一目标: P1d1
第二目标:有两个要求即甲
d 2,乙
d3,但两
个具有相同的优先因子,因此需要确定权系数。本题
在用效用函数作为规划目标时,需要确定一组权值 i
来反映原问题中各目标函数在总体目标中的权重,即:
k
max i i
i1
i ( x1, x2,xn ) gi (i 1,2,, m)
式中, i 应满足:
k
i 1
i 1
向量形式: max T
s.t. ( X ) G
方法二 理想点模型(罚款模型)
例如:在例一中,规定Z1 的目标值为 50000,正、负
偏差为d+、d- ,则目标函数可以转换为目标约束,既
70 x1 + 120 x2+ d1 d1=50000,
同样,若规定 Z2=200, Z3=250 则有
x1
d
2
d
2
200
x2 d3 d3 250
d
j
,
d
j
0
( j 1.2.3)
求解多目标规划的方法大体上有以下几种:
一种是化多为少的方法 , 即把多目标化为比较容易求解 的单目标或双目标,如主要目标法、线性加权法、理想点 法等;
另一种叫分层序列法,即把目标按其重要性给出一个序列, 每次都在前一目标最优解集内求下一个目标最优解,直到 求出共同的最优解。
对多目标的线性规划除以上方法外还可以适当修正单纯形 法来求解;还有一种称为层次分析法,是由美国运筹学家 沙旦于70年代提出的,这是一种定性与定量相结合的多目 标决策与分析方法,对于目标结构复杂且缺乏必要的数据 的情况更为实用。
思想: 规划决策者对每一个目标函数都能提出所期望的值 (或称满意值); 通过比较实际值 fi 与期望值 fi* 之间的偏差来选择问题的 解,其数学表达式如下:
k
min Z i ( fi fi )2
i ( x1, x2,i1 , xn ) gi (i 1,2,, m)
或写成矩阵形式: min Z (F F )T A(F F )
可用单件利润比作为权系数即 70 :120,化简为7:12。
P2
(7
d
2
12d
3
)
第三目标:
P3
(d
4
d
4
)
目标规划模型为:
min Z
P1d1
P2
(7d
2
12d
3
)
P3
(d
4
d
4
)
70
x1 x1
120
x2
d1
d
2
d1
d
2
50000 200
9x1
x2
d
3
d
3
250
3、线性规划中的约束条件是同等重要的,是硬约束; 而目标规划中有轻重缓急和主次之分,即有优先权。
4、线性规划的最优解是绝对意义下的最优,但需花 去大量的人力、物力、财力才能得到;实际过程中, 只要求得满意解,就能满足需要(或更能满足需要)。
目前,已经在经济计划、生产管理、经营管理、市场 分析、财务管理等方面得到了广泛的应用。
多目标规划方法
多目标规划是数学规划的一个分支。 研究多于一个的目标函数在给定区域上的最优化。又称多 目标最优化。通常记为
MOP(multi-objective programming)。 在很多实际问题中,例如经济、管理、军事、科学和工程 设计等领域,衡量一个方案的好坏往往难以用一个指标来 判断,而需要用多个目标来比较,而这些目标有时不甚协 调,甚至是矛盾的。因此有许多学者致力于这方面的研究。 1896年法国经济学家 V. 帕雷托最早研究不可比较目标的优 化问题,之后,J.冯·诺伊曼、H.W.库恩、A.W.塔克、A.M. 日夫里翁等数学家做了深入的探讨,但是尚未有一个完全 令人满意的定义。
乙 资源限制
4
3600
5
2000
10 3000
120
设:甲产品 x1 ,乙产品 x2
一般有:
同时:
maxZ=70 x1 + 120 x2 9 x1 +4 x2 ≤3600 4 x1 +5 x2 ≤ 2000 3 x1 +10 x2 ≤3000 x1 , x2 ≥0
maxZ1=70 x1 + 120x2 minZ2= x1 maxZ3= x2 9 x1 +4 x2 ≤3600 4 x1 +5 x2 ≤ 2000 3 x1 +10 x2 ≤3000 x1 , x2 ≥0