第1-2章 线性规划 整数规划(2)

2019-05-18 18:26

∑x

i=1

n

ij

=1

xij=0 或 1

上述指派问题的可行解可以用一个矩阵表示,其每行每列均有且只有一个元素为1,其余元素均为0;可以用1,L,n中的一个置换表示。

问题中的变量只能取0或1,从而是一个0-1规划问题。一般的0-1规划问题求解极为困难。但指派问题并不难解,其约束方程组的系数矩阵十分特殊(被称为全单位模矩阵,其各阶非零子式均为±1),其非负可行解的分量只能取0或1,故约束xij=0或1可改写为xij≥0而不改变其解。此时,指派问题被转化为一个特殊的运输问题,其中

m=n,ai=bj=1。

3.2 求解指派问题的匈牙利算法

由于指派问题的特殊性,又存在着由匈牙利数学家Konig提出的更为简便的解法—匈牙利算法。算法主要依据以下事实:如果系数矩阵C=(cij)一行(或一列)中每一元素都加上或减去同一个数,得到一个新矩阵B=(bij) ,则以C或B为系数矩阵的指派问题具有相同的最优指派。

例8 求解指派问题,其系数矩阵为

?16?17

C=?

?24??17?1?0

B1=?

?7??1

1521221919191822

22?18?? 17??16?

解 将第一行元素减去此行中的最小元素15,同样,第二行元素减去17,第三行元素减去17,最后一行的元素减去16,得

04534216

7?1?? 0??0?

再将第3列元素各减去1,得

?10*37??*?0411? B2=?*

?7500??*?1350??

以B2为系数矩阵的指派问题有最优指派

?1234?

??2134??

??

由等价性,它也是例7的最优指派。

有时问题会稍复杂一些。

例9 求解系数矩阵C的指派问题

-6-

?127979??89666?

??

C=?717121412?

??

15146610?????4107106?

解:先作等价变换如下

容易看出,从变换后的矩阵中只能选出四个位于不同行不同列的零元素,但n=5,

最优指派还无法看出。此时等价变换还可进行下去。步骤如下:

(1) 对未选出0元素的行打∨; (2) 对∨行中0元素所在列打∨;

(3) 对∨列中选中的0元素所在行打∨; 重复(2)、(3)直到无法再打∨为止。

可以证明,若用直线划没有打∨的行与打∨的列,就得到了能够覆盖住矩阵中所有零元素的最少条数的直线集合,找出未覆盖的元素中的最小者,令∨行元素减去此数,∨列元素加上此数,则原先选中的0元素不变,而未覆盖元素中至少有一个已转变为0,且新矩阵的指派问题与原问题也等价。上述过程可反复采用,直到能选取出足够的0元素为止。例如,对例5变换后的矩阵再变换,第三行、第五行元素减去2,第一列元素加上2,得

0?127979??50*2

?89666??2300*???

?7?717121412?→?0*1057

???

?6?15146610?80*0?9

??4?636?4107106???0

?7?6

2?0??5?∨? 4?2??∨

?7

?4? ?0??11??0

0

38842030100504

2?0??3? ?4?0??

现在已可看出,最优指派为??

?12345?

。 ??

?24135?

§4 对偶理论与灵敏度分析

4.1 原始问题和对偶问题

考虑下列一对线性规划模型:

max和

cTx s.t. Ax≤b,x≥0 (P)

minbTy s.t. ATy≥c,y≥0 (D)

-7-

称(P)为原始问题,(D)为它的对偶问题。

不太严谨地说,对偶问题可被看作是原始问题的“行列转置”:

(1) 原始问题中的第j列系数与其对偶问题中的第j行的系数相同; (2) 原始目标函数的各个系数行与其对偶问题右侧的各常数列相同; (3) 原始问题右侧的各常数列与其对偶目标函数的各个系数行相同; (4) 在这一对问题中,不等式方向和优化方向相反。 考虑线性规划:

mincTx

min

s.t.Ax=b,x≥0

把其中的等式约束变成不等式约束,可得

?A??b?

cTx s.t. ??x≥??,x≥0

??A???b?

它的对偶问题是

?y?

?AT?1?≤c

?y2?

令y=y1?y2,其中y1和y2分别表示对应于约束Ax≥b和?Ax≥?b的对偶变量组。

max

[b

T

?y?

?bT?1?

?y2?

]s.t. AT

[]则上式又可写成

maxbTy s.t. ATy≤c

原问题和对偶的对偶约束之间的关系:

min max

变量

行约束

?≥0?

?≤0 行约束?无限制??≥0?

?≤0 变量?=0?

?≤0??≥0 ?=0??≥0?

?≤0 ?无限制?

4.2 对偶问题的基本性质

1o 对称性:对偶问题的对偶是原问题。

2o 弱对偶性:若x是原问题的可行解,y是对偶问题的可行解。则存在

cTx≤bTy。

3o 无界性:若原问题(对偶问题)为无界解,则其对偶问题(原问题)无可行解。

?是原问题的可行解,y?是对偶问题的可行解,4o 可行解是最优解时的性质:设x

?=by?时,x?,y?是最优解。 当cx

5o 对偶定理:若原问题有最优解,那么对偶问题也有最优解;且目标函数值相同。

?,y?分别是原问题和对偶问题的最优解,则 6o 互补松弛性:若x

?(Ax??b)=0,x?(Ay??c)=0 y

例10 已知线性规划问题

minω=2x1+3x2+5x3+2x4+3x5

s.t. x1+x2+2x3+x4+3x5≥4

-8-

T

T

T

TT

2x1?x2+3x3+x4+x5≥3 xj≥0,已知其对偶问题的最优解为y1=解。

解 先写出它的对偶问题

maxz=4y1+3y2

y1+2y2≤2 ①

y1?y2≤3 ②

2y1+3y3≤5 ③ y1+y2≤2 ④ 3y1+y2≤3 ⑤ y1,y2≥0

将y1,y2的值代入约束条件,得②,③,④为严格不等式;由互补松弛性得

*****x2=x3=x4=0。因 y1,y2>0;原问题的两个约束条件应取等式,故有

**

x1+3x5=4 **2x1+x5=3

*

*

*

j=1,2,L,5

4*3

,y2=;z=5。试用对偶理论找出原问题的最优

55

求解后得到x1=1,x5=1;故原问题的最优解为

X=[10001]';ω=5 。

4.3 灵敏度分析

但实际上这些系数往往是估在以前讨论线性规划问题时,假定aij,bi,cj都是常数。计值和预测值。如市场条件一变,cj值就会变化;aij往往是因工艺条件的改变而改变;

*

*

**

bi是根据资源投入后的经济效果决定的一种决策选择。因此提出这样两个问题:当这

些系数有一个或几个发生变化时,已求得的线性规划问题的最优解会有什么变化;或者

这些系数在什么范围内变化时,线性规划问题的最优解或最优基不变。这里我们就不讨论了。

4.4 参数线性规划

参数线性规划是研究aij,bi,cj这些参数中某一参数连续变化时,使最优解发生变化的各临界点的值。即把某一参数作为参变量,而目标函数在某区间内是这参变量的线性函数,含这参变量的约束条件是线性等式或不等式。因此仍可用单纯形法和对偶单纯形法进行分析参数线性规划问题。

§5 投资的收益和风险

5.1 问题提出

市场上有n种资产si(i=1,2,L,n)可以选择,现用数额为M的相当大的资金风险损失率为作一个时期的投资。这n种资产在这一时期内购买si的平均收益率为ri,

qi,投资越分散,总的风险越少,总体风险可用投资的si中最大的一个风险来度量。

-9-

购买si时要付交易费,(费率pi),当购买额不超过给定值ui时,交易费按购买ui

计算。另外,假定同期银行存款利率是r0,既无交易费又无风险。(r0=5%)

已知n=4时相关数据如表1。

表1

qi

si

s1 s2 s3 s4

ri(%) pi(%) ui(元)

28 2.5 1 103 21 1.5 2 198 23 5.5 4.5 52 25 2.6 6.5 40

试给该公司设计一种投资组合方案,即用给定资金M,有选择地购买若干种资产或存银行生息,使净收益尽可能大,使总体风险尽可能小。

5.2 符号规定和基本假设 符号规定:

si:第i种投资项目,如股票,债券

ri,pi,qi:分别为si的平均收益率,交易费率,风险损失率 ui:si的交易定额 r0:同期银行利率

xi:投资项目si的资金 a:投资风险度 Q:总体收益

基本假设:

1. 投资数额M相当大,为了便于计算,假设M=1; 2. 投资越分散,总的风险越小;

3. 总体风险用投资项目si中最大的一个风险来度量; 4. n种资产si之间是相互独立的;

5. 在投资的这一时期内,ri,pi,qi,r0为定值,不受意外因素影响; 6. 净收益和总体风险只受ri,pi,qi影响,不受其它因素干扰。 5.3 模型的分析与建立

1. 总体风险用所投资的si中最大的一个风险来衡量,即 max{qixi|i=1,2,L,n} 2.购买si所付交易费是一个分段函数,即 交易费=?

?pixi,?piui,

xi>uixi≤ui

而题目所给定的定值ui(单位:元)相对总投资M很少,piui更小,可以忽略不计,这样购买si的净收益为(ri?pi)xi。

3. 要使净收益尽可能大,总体风险尽可能小,这是一个多目标规划模型:

-10-


第1-2章 线性规划 整数规划(2).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:关于做好春节期间走访慰问老干部工作的通知

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: