第三讲经济与金融中的优化模型

1、第三讲第三讲 经济与金融中的优化问题经济与金融中的优化问题 本章主要介绍用本章主要介绍用LINDO/LINGO软件求解经济、软件求解经济、金融和市场营销方面的几个优化问题的案例。金融和市场营销方面的几个优化问题的案例。3.1 经济均衡问题极其应用经济均衡问题极其应用 在市场经济活动中，当市场上某种产品的价格越高时，生在市场经济活动中，当市场上某种产品的价格越高时，生产商越是愿意扩大生产能力产商越是愿意扩大生产能力(供应能力供应能力)，提高更多的产品，提高更多的产品满足市场需求；但市场价格太高时，消费者的消费欲望满足市场需求；但市场价格太高时，消费者的消费欲望(需求能力需求能力)会下降。反之，当

2、市场上某种商品的价格越低会下降。反之，当市场上某种商品的价格越低时，消费者的消费欲望时，消费者的消费欲望(需求能力需求能力)会上升，但生产商的供会上升，但生产商的供应能力会下降。如果生产商的供应能力和消费者的需求能应能力会下降。如果生产商的供应能力和消费者的需求能力长期不匹配，就会导致经济不稳定。在完全市场竞争的力长期不匹配，就会导致经济不稳定。在完全市场竞争的环境中，我们总是认为经济活动应当达到均衡环境中，我们总是认为经济活动应当达到均衡(equilibrium)，即生产和消费，即生产和消费(供应能力和需求能力供应能力和需求能力)达到达到平衡，不再发生变化，这时该产品的价格就是市场的清算平衡

3、，不再发生变化，这时该产品的价格就是市场的清算价格。价格。 下面考虑两个简单的单一市场及双边市场的具体实例。并下面考虑两个简单的单一市场及双边市场的具体实例。并介绍经济均衡思想在拍卖与投标问题、交通流分配问题中介绍经济均衡思想在拍卖与投标问题、交通流分配问题中的应用案例。的应用案例。生产商（甲）生产商（甲）消费者（乙）消费者（乙）单价单价/（万元（万元/t）供应能力供应能力/t单价单价/（元（元/t）需求能力需求能力/t1292244.543636482.2583.1.1 单一生产商、单一消费者的情形单一生产商、单一消费者的情形例例3.1 假设市场上只有一个生产商假设市场上只有一个生产商(记为

4、甲记为甲)和一个消费和一个消费者者(记为乙记为乙)。对某种产品，他们在不同价格下的供应能。对某种产品，他们在不同价格下的供应能力和需求能力如表力和需求能力如表3-1所示。举例来说，表中数据的含所示。举例来说，表中数据的含义是：当单价低于义是：当单价低于2万元但大于或等于万元但大于或等于1万元时，甲愿万元时，甲愿意生产意生产2t产品，乙愿意购买产品，乙愿意购买8t产品；当价格低于产品；当价格低于9万元万元但大于或等于但大于或等于4.5万元时，乙愿意购买万元时，乙愿意购买4t产品，甲愿意产品，甲愿意生产生产8t产品；依次类推。那么市场的清算价格应该是多产品；依次类推。那么市场的清算价格应该是多少？

5、少？ 问题分析问题分析 仔细观察一下表仔细观察一下表3-1就可以看出来，这个具体问题的解是就可以看出来，这个具体问题的解是一目了然的：清算价格显然应该是一目了然的：清算价格显然应该是3万元万元/t，因为此时供求，因为此时供求平衡平衡(都是都是6t)。为了能够处理一般情况，下面通过建立优。为了能够处理一般情况，下面通过建立优化模型来解这个问题。化模型来解这个问题。 这个问题给人的第一印象似乎没有明确的目标函数，不太这个问题给人的第一印象似乎没有明确的目标函数，不太像是一个优化问题。不过，我们可以换一个角度来想问题：像是一个优化问题。不过，我们可以换一个角度来想问题：假设市场上还有一个虚拟的经销商

6、，他是甲乙进行交易的假设市场上还有一个虚拟的经销商，他是甲乙进行交易的中介。那么，为了使自己获得的利益最大，他将总是以可中介。那么，为了使自己获得的利益最大，他将总是以可能的最低价格从甲购买产品，在以可能的最高价格卖给乙，能的最低价格从甲购买产品，在以可能的最高价格卖给乙，直到进一步的交易无利可图为止。例如，最开始的直到进一步的交易无利可图为止。例如，最开始的2t产品产品他会以他会以1万元的单价从甲购买，以万元的单价从甲购买，以9万元的单价卖给乙；接万元的单价卖给乙；接下来的下来的2t产品会以产品会以2万元的单价从甲购买，再以万元的单价从甲购买，再以4.5万元的万元的单价卖给乙；再接下来的单价

7、卖给乙；再接下来的2t产品他只能以产品他只能以3万元的单价从万元的单价从甲购买，再以甲购买，再以3万元的单价卖给乙万元的单价卖给乙(其实这次交易他已经只其实这次交易他已经只是保本，但我们仍然是假设这笔交易会发生，例如他为了是保本，但我们仍然是假设这笔交易会发生，例如他为了使自己的营业额尽量大使自己的营业额尽量大)；最后，如果他继续购买甲的产；最后，如果他继续购买甲的产品卖给乙，他一定会亏本，所以他肯定不会交易。因此，品卖给乙，他一定会亏本，所以他肯定不会交易。因此，市场清算价格就是市场清算价格就是3万元。根据这个想法，我们就可以建万元。根据这个想法，我们就可以建立这个问题的线性规划模型。立这个

8、问题的线性规划模型。 模型建立模型建立 决策变量：社甲以决策变量：社甲以1万元，万元，2万元，万元，3万元，万元，4万元的单价万元的单价售出的产品数量售出的产品数量(单位：单位：t)分别是分别是A1，A2，A3，A4，乙以，乙以9万元，万元，4.5万元，万元，3万元，万元，2.25万元的单价购买的产品数万元的单价购买的产品数量量(单位：单位：t)分别是分别是x1,x2,x3,x4。 目标函数：就是虚拟经销商的总利润，即目标函数：就是虚拟经销商的总利润，即 9x1+4.5x2+3x3+2.5x4-A1-2A2-3A3-4A4. 约束条件：约束有约束条件：约束有 供需平衡供需平衡 A1+A2+A3

9、+A4=x1+x2+x3+x4； 供应限制供应限制 A1,A2,A3,A4=2; 消费限制消费限制 x1,x2,x3,x4=0。 模型求解：（模型求解：（LINGO程序）程序） MAX 9X1 + 4.5X2 + 3X3 + 2.25X4 - A1 - 2A2 - 3A3 - 4A4 SUBJECT TO A1 + A2 + A3 + A4 - X1 - X2 - X3 - X4 = 0 A1 = 2 A2 = 2 A3 = 2 A4 = 2 X1 = 2 X2 = 2 X3 = 2 X4 = 2 END 运行分析其结果。运行分析其结果。生产商（丙）生产商（丙）消费者（丁）消费者（丁）价格价格

10、/万元万元供应能力供应能力/t价格价格/元元需求能力需求能力/t21151448368568123103.1.2 两个生产商、两个消费者的情形两个生产商、两个消费者的情形 例例3.2 假设市场上除了假设市场上除了3.1中的甲和乙外，还有另中的甲和乙外，还有另一个生产商一个生产商(记为丙记为丙)和另一个消费者和另一个消费者(记为丁记为丁)，他，他们在不同价格下的供应能力如表们在不同价格下的供应能力如表3-2所示。此外，所示。此外，从甲销售到丁的每吨产品的运输成本是从甲销售到丁的每吨产品的运输成本是1.5万元，万元，从丙销售到乙的每吨产品的运输成本是从丙销售到乙的每吨产品的运输成本是2万元，而万元

11、，而甲、乙之间没有运输成本，丙、丁之间没有运输甲、乙之间没有运输成本，丙、丁之间没有运输成本。这时，市场的清算价格应该是多少？甲和成本。这时，市场的清算价格应该是多少？甲和丙分别生产多少？乙和丁分别购买多少？丙分别生产多少？乙和丁分别购买多少？ 表表3-2 不同价格下的供应能力和消费能力不同价格下的供应能力和消费能力 问题分析问题分析 首先，我们看看为什么要考虑从甲销售到丁的产首先，我们看看为什么要考虑从甲销售到丁的产品的运输成本和从丙销售到乙的产品的运输成本。品的运输成本和从丙销售到乙的产品的运输成本。如果不考虑这些运输成本，我们就可以认为甲乙如果不考虑这些运输成本，我们就可以认为甲乙丙丁处

12、于同一个市场上，因此可以将两个生产商丙丁处于同一个市场上，因此可以将两个生产商(甲和丙甲和丙)的供应函数合并成一个供应函数，合并的供应函数合并成一个供应函数，合并后就可以认为市场上仍然只有一个供应商。类似后就可以认为市场上仍然只有一个供应商。类似地，乙和丁的需求函数也可以合并成一个需求函地，乙和丁的需求函数也可以合并成一个需求函数，合并后就可以认为市场上仍然是只有一个消数，合并后就可以认为市场上仍然是只有一个消费者。这样就回到了例费者。这样就回到了例3.1的情形。的情形。 也就是说，考虑运输成本在经济学上的含义，应也就是说，考虑运输成本在经济学上的含义，应当是认为甲乙是一个市场当是认为甲乙是一

13、个市场(地区或国家地区或国家)，而丙丁，而丙丁是另一个市场是另一个市场(地区或国家地区或国家)，运输成本也可能还，运输成本也可能还包括关税等成本，由于这个成本的存在，两个市包括关税等成本，由于这个成本的存在，两个市场的清算价可能是不同的。场的清算价可能是不同的。 仍然按照仍然按照3.1的思路，可以建立这个问题的线性的思路，可以建立这个问题的线性规划模型。规划模型。 模型建立和求解模型建立和求解 设甲以设甲以1，2，3，4(万元万元)的单价售出的产品数量的单价售出的产品数量(单单位：位：t)分别是分别是 A1,A2,A3,A4，乙以，乙以9，4.5，3，2.25(万元万元)的单价购买的产品数量的

14、单价购买的产品数量(单位：单位：t)分别是分别是X1,X2,X3,X4；丙以；丙以2，4，6，8(万元万元)的单价售出的的单价售出的产品数量产品数量(单位：单位：t)分别是分别是B1,B2,B3,B4,丁以丁以15，8，5，3(万元万元)的单价购买的产品数量的单价购买的产品数量(单位：单位：t)分别是分别是Y1,Y2,Y3,Y4,此外，假设此外，假设AX和和AY分别是甲向乙和丁分别是甲向乙和丁的供货量，的供货量，BX和和BY分别是丙向丁的供货量。这些决分别是丙向丁的供货量。这些决策变量之间的关系参见示意图策变量之间的关系参见示意图3-1。 目标函数仍然是虚拟经销商的总利润，约束条件仍然目标函数

15、仍然是虚拟经销商的总利润，约束条件仍然是四类是四类(供求平衡、供应限制、需求限制和非负限制供求平衡、供应限制、需求限制和非负限制)，不过这时应该注意供需平衡约束应该是包括图不过这时应该注意供需平衡约束应该是包括图3-1所示所示的决策变量之间的关系：的决策变量之间的关系： AX+AY=A1+A2+A3+A4， (10) BX+BY=B1+B2+B3+B4， (11) AX+BX=x1+x2+x3+x4， (12) AY+BY=Y1+Y2+Y3+Y4。 (13) 此外的其他约束实际上只是一个简单的变量上界约束，此外的其他约束实际上只是一个简单的变量上界约束，可以用可以用“SUB”命令表示。命令表示

16、。 程序：程序：MAX 9X1 + 4.5X2 + 3X3 + 2.25X4 + 15Y1 + 8Y2 + 5Y3 + 3Y4 - 2BX - 1.5AY - A1 - 2A2 - 3A3 - 4A4 - 2B1 - 4B2 - 6B3 - 8B4 SUBJECT TO 2) - AY + A1 + A2 + A3 + A4 - AX = 0 3) - BX + B1 + B2 + B3 + B4 - BY = 0 4) - X1 - X2 - X3 - X4 + BX + AX = 0 5) - Y1 - Y2 - Y3 - Y4 + AY + BY = 0 ENDSUB A1 2 SUB

17、A2 2SUB A3 2SUB A4 2 SUB X1 2SUB X2 2SUB X3 2SUB X4 2SUB B1 1 SUB B2 3SUB B3 4SUB B4 4 SUB Y1 1SUB Y2 2SUB Y3 3SUB Y4 4运行分析其结果。运行分析其结果。3.1.3 拍卖和投标问题拍卖和投标问题 例例3.33.3 假设一家拍卖行对委托的假设一家拍卖行对委托的5类艺术品对外类艺术品对外拍卖，采用在规定日期前投标人提交投标书的方拍卖，采用在规定日期前投标人提交投标书的方式进行，最后收到了来自式进行，最后收到了来自4个投标人的投标书。每个投标人的投标书。每类项目的数量、投标人对每个项目

18、的投标价格如类项目的数量、投标人对每个项目的投标价格如表表3-3中所示。例如，有中所示。例如，有3件第件第4类艺术品；对每件类艺术品；对每件第第4类艺术品；对每件第类艺术品；对每件第4类艺术品，投标人类艺术品，投标人1，2，3，4愿意出的最高价分别为愿意出的最高价分别为6，1，3，2(货币单位，货币单位，如万元如万元)。此外，假设每个投标对人每类艺术品最。此外，假设每个投标对人每类艺术品最多只能购买一件，并且每个投标人购买的艺术品多只能购买一件，并且每个投标人购买的艺术品总数不能超过总数不能超过3件，那么，哪些艺术品能够卖出去？件，那么，哪些艺术品能够卖出去？卖给谁？这个拍卖和投标问题中每类物

19、品的清算卖给谁？这个拍卖和投标问题中每类物品的清算价应该是多少？价应该是多少？招标项目类型招标项目类型12345招标项目的数量招标项目的数量12334投标价格投标价格投标人投标人192863投标人投标人267915投标人投标人378634投标人投标人454321表表3-3 拍卖与投标信息拍卖与投标信息 问题分析问题分析 这个具体问题在实际中可能以通过对所这个具体问题在实际中可能以通过对所有投标的报价进行排序来解决，例如可有投标的报价进行排序来解决，例如可以总是将艺术品优先卖给出价最高的投以总是将艺术品优先卖给出价最高的投标人。但这种方法不太好确定每类艺术标人。但这种方法不太好确定每类艺术品的清

20、算价，所以我们这里还是借用前品的清算价，所以我们这里还是借用前面两个例子中的方法，即使设有一个中面两个例子中的方法，即使设有一个中间商希望最大化自己的利润，从而建立间商希望最大化自己的利润，从而建立这个问题的线性规划模型。这个问题的线性规划模型。 问题的一般提法和假设问题的一般提法和假设 先建立一般的模型，然后求解本例的具先建立一般的模型，然后求解本例的具体问题。体问题。 设有设有N类物品需要拍卖，第类物品需要拍卖，第j类物品的数量为类物品的数量为Sj(j=1，2，.，N)；有；有M个投标者，投标者个投标者，投标者i(i=1，2，.，M)对对第第j类物品的投标价格为类物品的投标价格为bij(假

21、设非负假设非负)。投标者。投标者i对每类对每类物品最多购买一件，切总件数不能超过物品最多购买一件，切总件数不能超过ci。我们的目标。我们的目标之一是要确定第之一是要确定第j类物品的清算价格类物品的清算价格pj，他应当满足下，他应当满足下列列假设条件假设条件： (1)成交的第成交的第j类物品的数量不超过类物品的数量不超过Sj(j=1，2，.，N)； (2)对第对第j类物品的报价低于类物品的报价低于pj的投标人将不能获得第的投标人将不能获得第j类类物品；物品； (3)如果成交的第如果成交的第j类物品的数量少于类物品的数量少于Sj(j=1，2，.，N)，可以认为可以认为pj=0(除非拍卖方另外指定一

22、个最低的保护价除非拍卖方另外指定一个最低的保护价)； (4)对对j类物品的报价高于类物品的报价高于pj的投标人有权获得第的投标人有权获得第j类物品，类物品，但如果他有权获得的物品超过但如果他有权获得的物品超过3件，那么我们假设他总件，那么我们假设他总是希望使自己的满意度最大是希望使自己的满意度最大(满意度可以用他的报价与满意度可以用他的报价与市场清算价之差来衡量市场清算价之差来衡量)。 优化模型优化模型 用用0-1变量变量xij表示是否分配一件第表示是否分配一件第j类物品给投标者类物品给投标者i，即即xij=1表示分配，而表示分配，而xij=0表示不分配。目标函数表示不分配。目标函数仍然是虚拟

23、的中间商的总利润仍然是虚拟的中间商的总利润(认为这些利润全部是认为这些利润全部是拍卖的利润也可以拍卖的利润也可以)，即，即 bij*xij。 (14) 除变量取值为除变量取值为0或或1的约束外，问题的约束条件只要的约束外，问题的约束条件只要是两类；每类物品的数量限制和每个投标人所能分是两类；每类物品的数量限制和每个投标人所能分到的物品的数量限制，即到的物品的数量限制，即 xij=Sj,j=1,2,.,N; (15) xij=cj,j=1,2,.,M; (16) 模型就是在约束模型就是在约束(15)、(16)下最大化目标函数下最大化目标函数(14)。 程序：程序： MODEL: TITLE 拍卖

24、与投标拍卖与投标; SETS:AUCTION: S;BIDDER : C;LINK(BIDDER,AUCTION): B, X; ENDSETS DATA:AUCTION=FILE(AUCTION.TXT);BIDDER =FILE(AUCTION.TXT);S=FILE(AUCTION.TXT);C=FILE(AUCTION.TXT);B=FILE(AUCTION.TXT); ENDDATA MAX=SUM(LINK: B*X); FOR(AUCTION(J): AUC_LIM SUM(BIDDER(I): X(I,J) S(J) ); FOR(BIDDER(I): BID_LIM SUM(

25、AUCTION(J): X(I,J) C(I) ); FOR(LINK: BND(0,X,1); END Global optimal solution found at iteration: 7 Objective value: 65.00000 运行分析其结果。运行分析其结果。 3.1.4 交通流均衡问题交通流均衡问题 例例3.43.4 某地有如图某地有如图3-2所示的一个公路网，每所示的一个公路网，每天上班时间有天上班时间有6千辆小汽车要从居民区千辆小汽车要从居民区A前往工前往工作区作区D。经过长期观察，我们得到了图中。经过长期观察，我们得到了图中5条道条道路上每辆汽车的平均行驶时间和汽

26、车流量之间路上每辆汽车的平均行驶时间和汽车流量之间的关系，如表的关系，如表3-4所示。那么，长期来看，这些所示。那么，长期来看，这些汽车将如何在每条道路上分布？汽车将如何在每条道路上分布？BADC道路道路ABACBCBDCD行驶时间行驶时间/min流量流量220521252202流量流量330531353303流量流量44054145440居民区居民区 工作区工作区 问题分析问题分析 这个问题看起来似乎与前面几个例子中的完全不同，这个问题看起来似乎与前面几个例子中的完全不同，但实际上交通流与市场经济活动类似，也存在着均衡。但实际上交通流与市场经济活动类似，也存在着均衡。 我们可以想象有一个协调

27、者，正如前面几个例子中的我们可以想象有一个协调者，正如前面几个例子中的所谓中间商可以理解为市场规律一样，实际上这里的所谓中间商可以理解为市场规律一样，实际上这里的所谓协调者也可以认为交通流的规律。交通流的规律所谓协调者也可以认为交通流的规律。交通流的规律就是每辆汽车都将选择自己使自己从就是每辆汽车都将选择自己使自己从A到到D运行时间最运行时间最少的路线，其必然的结果是无论走哪条路线从少的路线，其必然的结果是无论走哪条路线从A到到D，最终花费的时间应该是一样的最终花费的时间应该是一样的(否则，花费时间较长的否则，花费时间较长的那条线路上的部分汽车就会改变自己的路线，以缩短那条线路上的部分汽车就会

28、改变自己的路线，以缩短自己的行驶时间自己的行驶时间)。 也就是说，长期来看，这些汽车在每条道路上的分布也就是说，长期来看，这些汽车在每条道路上的分布将达到均衡状态将达到均衡状态(所谓均衡，这里的含义就是每辆汽车所谓均衡，这里的含义就是每辆汽车都不能仅仅通过自身独自改变道路节省其行使时间都不能仅仅通过自身独自改变道路节省其行使时间)。在这种想法下，我们来建立线性规划模型。在这种想法下，我们来建立线性规划模型。 优化模型优化模型 交通流的规律要求所有道路上的流量达到均衡，我们仍交通流的规律要求所有道路上的流量达到均衡，我们仍然类似例然类似例1.1和例和例1.2来考虑问题。如果车流量是一辆一辆来考虑

29、问题。如果车流量是一辆一辆增加的，那么在每条道路上车流量少于增加的，那么在每条道路上车流量少于2时，车流量会有时，车流量会有一个分布规律；但某条道路上的车流量正好超过一个分布规律；但某条道路上的车流量正好超过2时，新时，新加入的一辆车需要选择使自己堵塞时间最短的道路，这就加入的一辆车需要选择使自己堵塞时间最短的道路，这就提示我们把同一条道路上的流量分布分解成不同性质的三提示我们把同一条道路上的流量分布分解成不同性质的三个部分，也就是说，我们用个部分，也就是说，我们用Y(AB)表示道路表示道路AB上的总的上的总的流量，并进一步把他分解成三个部分：流量，并进一步把他分解成三个部分： (1)道路道路

30、AB上的流量不超过上的流量不超过2时的流量，用时的流量，用X(2，AB)表示；表示； (2)道路道路AB上的流量超过上的流量超过2但不超过但不超过3时，超过时，超过2的流量部的流量部分用分用X(3，AB)表示；表示； (3)道路道路AB上的流量超过上的流量超过3但不超过但不超过4时，超过时，超过3的流量部的流量部分用分用X(4，AB)表示。表示。 以此类推，对道路以此类推，对道路 AC，BC，BD，CD上同理可以定义类上同理可以定义类似的决策变量，因此，问题中总共有似的决策变量，因此，问题中总共有20个决策变量个决策变量Y(j)和和X(i，j)(i=2，3，4；j=AB，AC，BC，BD，CD

31、)。 问题的目标应当是使总的堵塞时间最小。用问题的目标应当是使总的堵塞时间最小。用T(i，j)表示流表示流量量X(i，j)对应的堵塞时间对应的堵塞时间 (即表即表3-3中的数据，是对每辆车中的数据，是对每辆车而言的而言的)，我们看看用，我们看看用T(i，j)X(i，j)作为总堵塞时间是作为总堵塞时间是否合适。很容易理解：后面加入道路的车辆可能又会造成否合适。很容易理解：后面加入道路的车辆可能又会造成前面进入道路的车辆的进一步堵塞，如流量为前面进入道路的车辆的进一步堵塞，如流量为3时，原先时，原先流量为流量为2的车辆实际上也只能按的车辆实际上也只能按T(3，j)的时间通过，而不的时间通过，而不是

32、是T(2，j)。也即使说，。也即使说，T(i，j)X(i，j)并不是总堵塞时间，并不是总堵塞时间，但是我们也可以发现，但是我们也可以发现，T(i，j)关于关于i是单调增加的，即不断是单调增加的，即不断增加的车流只会使以前的堵塞加剧而不可能使以前的堵塞增加的车流只会使以前的堵塞加剧而不可能使以前的堵塞减缓。所以关于决策变量减缓。所以关于决策变量X(i，j)而言，而言，T(i，j)X(i，j)与与我们希望优化的目标的单调性是一致的。因此，可以用我们希望优化的目标的单调性是一致的。因此，可以用T(i，j)X(i，j)作为目标函数进行优化。作为目标函数进行优化。 约束条件有三类：约束条件有三类： (1

33、)每条道路上的总流量每条道路上的总流量Y等于该道路上的分流量等于该道路上的分流量X的和；的和； (2)道路交汇处道路交汇处A，B，C，D(一般称为节点一般称为节点)的流量守恒的流量守恒(即即流入量等于流出量流入量等于流出量)； (3)决策变量的上限限制的，如决策变量的上限限制的，如X(2，AB)=2,X(3，AB)=1，X(4，AB)= 0Fx2) .02481442 x1 + .11678340 x2 + .11085278 x3 + PBUD - 1.213667 PRET = 0Fx3) .02615026 x1 + .11085278 x2 + .18845362 x3 + PBUD

34、- 1.234583 PRET = 0BUD) x1 + x2 + x3 = 1RET) 1.089083 x1 + 1.213667 x2 + 1.234583 x3 = 1.2345END! 告诉告诉LINDO真正的约束的起始行真正的约束的起始行QCP 5 LINGO程序：程序： MODEL: Title 简单的投资组合模型简单的投资组合模型; SETS: YEAR/1.12/; STOCKS/ A, B, C/: Mean,X; link(YEAR, STOCKS): R; STST(Stocks,stocks): COV; ENDSETS DATA: TARGET = 1.15; !

35、R是原始数据是原始数据; R = 1.300 1.225 1.149 1.103 1.290 1.260 1.216 1.216 1.419 0.954 0.728 0.922 0.929 1.144 1.169 1.056 1.107 0.965 1.038 1.321 1.133 1.089 1.305 1.732 1.090 1.195 1.021 1.083 1.390 1.131 1.035 0.928 1.006 1.176 1.715 1.908; ENDDATA CALC:!计算均值向量计算均值向量Mean与协方差矩阵与协方差矩阵COV; for(stocks(i): Mean

36、(i) =sum(year(j): R(j,i) / size(year) ); for(stst(i,j): COV(i,j) = sum(year(k): (R(k,i)-mean(i)*(R(k,j)-mean(j) / (size(year)-1) ); ENDCALC OBJ MIN = sum(STST(i,j): COV(i,j)*x(i)*x(j); ONE SUM(STOCKS: X) = 1; TWO SUM(stocks: mean*x) = TARGET; END 运行后，可以看到投资组合的决策结果与运行后，可以看到投资组合的决策结果与LINDO模型的输模型的输出相同（

37、只有很小的计算误差）。出相同（只有很小的计算误差）。 模型也输出了均值向量模型也输出了均值向量Mean和协方差矩阵和协方差矩阵COV，结果与，结果与我们前面给出的值是一致的。此外，请注意模型中计算协我们前面给出的值是一致的。此外，请注意模型中计算协方差矩阵方差矩阵COV时，分母是样本数减去时，分母是样本数减去1（即（即“size(year)-1”）而不是样本数，这是常用的计算方法，）而不是样本数，这是常用的计算方法，主要是为了保持这个估计的无偏性（当然，样本数较大时主要是为了保持这个估计的无偏性（当然，样本数较大时两者差别不大）。两者差别不大）。 用用LINDO软件对模型进行参数分析软件对模型

38、进行参数分析 对实际投资人来说，可能不仅希望知道制定的期望投资回对实际投资人来说，可能不仅希望知道制定的期望投资回报率下的风险（回报率的方差），可能更希望知道风险随报率下的风险（回报率的方差），可能更希望知道风险随着不同的投资回报率是如何变化的，然后作出最后的投资着不同的投资回报率是如何变化的，然后作出最后的投资决策。决策。 这当然可以通过在上面的模型中不断修改约束中的参这当然可以通过在上面的模型中不断修改约束中的参数（目前为数（目前为1.15）来实现，如将）来实现，如将1.15改为改为1.2345，则，则表示投资回报率希望达到表示投资回报率希望达到23.45%（这几乎是可能达（这几乎是可能达

39、到的最大值了，因为这几乎是三种股票中最大的投资到的最大值了，因为这几乎是三种股票中最大的投资回报率，即股票回报率，即股票C的回报率）。可以想到，这时应主的回报率）。可以想到，这时应主要投资在股票要投资在股票C上。实际求解一下，可以知道最优解上。实际求解一下，可以知道最优解中投资股票中投资股票C的份额大约是的份额大约是99.6%（剩余的大约（剩余的大约0.4%投资在股票投资在股票B上）。上）。 实际上，实际上，LINDO软件可以直接完成这种参数分析过程（目软件可以直接完成这种参数分析过程（目前前LINGO软件似乎还没有这个功能）。假设利用上面的软件似乎还没有这个功能）。假设利用上面的LINDO模

40、型，对于投资回报率希望为模型，对于投资回报率希望为23.45%的情形我们的情形我们已经求解得到了结果。此时，执行菜单命令已经求解得到了结果。此时，执行菜单命令“Reports|Parametrics（参数分析）（参数分析）”，将会看到图，将会看到图3-3所示的对话框。选择约束行所示的对话框。选择约束行“RET”；然后输入新的右；然后输入新的右端项（端项（New RHS Value），我们输入），我们输入1（表示收益率为）；（表示收益率为）；选择报告类行为二维图形（选择报告类行为二维图形（“Graphics”+“2D”）；最后）；最后按按“OK”按钮即可。屏幕上将显示如图按钮即可。屏幕上将显示如

41、图3-4所示的参数分所示的参数分析的结果，图中横坐标表示的是希望达到的回报率，纵坐析的结果，图中横坐标表示的是希望达到的回报率，纵坐标（目标函数）表示的是对应的方差，这种图形在经济学标（目标函数）表示的是对应的方差，这种图形在经济学上一般被称为有效前沿面（上一般被称为有效前沿面（efficient frontier）。可以看）。可以看出，图中曲线有两个明显的转折点，当要求的回报率越过出，图中曲线有两个明显的转折点，当要求的回报率越过这两个点以后，风险（方差）增长越来越快。这两个点以后，风险（方差）增长越来越快。 那么，投资组合的决策是如何变化的呢？如果在图那么，投资组合的决策是如何变化的呢？如

42、果在图3-3所所示的参数分析对话框中选择示的参数分析对话框中选择“Text（文本）（文本）”选项，则参选项，则参数分析的结果也会显示在结果报告窗口中：数分析的结果也会显示在结果报告窗口中： 这个结果的中间这个结果的中间4行的第行的第1行说明当希望的回报率从行说明当希望的回报率从23.45%下降到下降到21.894%时，变量时，变量X1进基，即还需要购买进基，即还需要购买股票股票A；第；第2行说明当希望的回报率继续下降到行说明当希望的回报率继续下降到9.357%时，时，变量变量X3出基，即不再需要购买股票出基，即不再需要购买股票C；第；第34行说明当希行说明当希望的回报率继续下降到望的回报率继续

43、下降到8.908%或以下时，变量或以下时，变量X2出基，出基，即只需要购买股票即只需要购买股票A（注意我们前面还假设了所有资金必（注意我们前面还假设了所有资金必须全部投资到这三种股票上，没有其他投资方式或让自己须全部投资到这三种股票上，没有其他投资方式或让自己闲置）。这几个关键点正是我们在图闲置）。这几个关键点正是我们在图3-4中观察到的曲线中观察到的曲线转折点。转折点。3.2.2 存在无风险资产时的投资组合模型存在无风险资产时的投资组合模型 例例3.6 假设除了例假设除了例3.5中的三种股票外，投资人还中的三种股票外，投资人还有一种无风险的投资方式，如购买国库券，假设有一种无风险的投资方式，

44、如购买国库券，假设国库券的年收益率为国库券的年收益率为5%，如何考虑例，如何考虑例3.5中的问中的问题？题？ 问题分析问题分析 其实，无风险的投资方式（如国库券、银行存款其实，无风险的投资方式（如国库券、银行存款等）是有风险的投资方式（如股票）的一种特例，等）是有风险的投资方式（如股票）的一种特例，所以这就意味着例所以这就意味着例3.5中的模型仍然是适用的，只中的模型仍然是适用的，只不过无风险的投资方式的收益是固定的，所以方不过无风险的投资方式的收益是固定的，所以方差（包括它与其它投资方式的收益协方差）都是差（包括它与其它投资方式的收益协方差）都是0。3.2.3 考虑交易成本的投资组合模型考虑

45、交易成本的投资组合模型 例例3.73.7 继续考虑例继续考虑例3.5（期望收益率仍定为（期望收益率仍定为15%）。）。假设你目前持有的股票比例为：股票假设你目前持有的股票比例为：股票A占占50%，B占占35%，C占占15%。这个比例与例。这个比例与例3.5中得到的最中得到的最优解有所不同，但实际股票市场上每次股票买卖优解有所不同，但实际股票市场上每次股票买卖通常总有交易费，例如按交易额的通常总有交易费，例如按交易额的1%收取交易费，收取交易费，这时你是否需要对所持有的股票进行买卖（换这时你是否需要对所持有的股票进行买卖（换手），以便满足手），以便满足“最优解最优解”的要求？的要求？ 建立模型建

46、立模型 仍用决策变量仍用决策变量x1，x2和和x3分别表示投资人应当投资股票分别表示投资人应当投资股票A、B、C的比例，进一步假设购买股票的比例，进一步假设购买股票A、B、C的比例为的比例为y1，y2和和y3，卖出股票，卖出股票A、B、C的比例为的比例为z1，z2和和z3。其中，。其中，yi与与zi（i=1，2，3）中显然最多只能有一个严格取正数，）中显然最多只能有一个严格取正数，且且 x1，x2，x30， y1，y2，y30， z1，z2，z30。 由于交易费用的存在，这是约束由于交易费用的存在，这是约束x1+x2+x3=1不一定还成不一定还成立（只有不进行股票买卖，即立（只有不进行股票买卖

47、，即y1=y2=y3=z1=z2=z3=0时，时，这个约束才成立）。其实，这个关系式的本质是：当前持这个约束才成立）。其实，这个关系式的本质是：当前持有的总资金是守恒的（假设为有的总资金是守恒的（假设为“1个单位个单位”），再有交易），再有交易成本（成本（1%）的情况下，应当表示成如下形式：）的情况下，应当表示成如下形式： x1+x2+x3+0.01（y1+y2+y3+z1+z2+z3）=1. 另外，考虑到当前持有的各只股票的份额另外，考虑到当前持有的各只股票的份额ci，xi，yi与与zi（i=1，2，3）之间也应该满足守恒关系式）之间也应该满足守恒关系式 Xi=ci+yi-zi, i=1,2

48、,3. 这就是新问题的约束条件，模型的其他部分不用改变。这就是新问题的约束条件，模型的其他部分不用改变。 程序：程序： MODEL: Title 考虑交易费的投资组合模型考虑交易费的投资组合模型; SETS: STOCKS/ A, B, C/: C,Mean,X,Y,Z; STST(Stocks,stocks): COV; ENDSETS DATA: TARGET = 1.15; ! 股票的初始份额股票的初始份额; c=0.5 0.35 0.15; ! Mean是收益均值是收益均值,COV是协方差矩阵是协方差矩阵; mean=1.089083 1.213667 1.234583; COV=0.

49、01080754 0.01240721 0.01307513 0.01240721 0.05839170 0.05542639 0.01307513 0.05542639 0.09422681; ENDDATA OBJ MIN = sum(STST(i,j): COV(i,j)*x(i)*x(j); ONE SUM(STOCKS: X+0.01*Y+0.01*Z) = 1; TWO SUM(stocks: mean*x) = TARGET; FOR(stocks: ADD x = c - y + z); END 运行并分析其结果。运行并分析其结果。 3.2.4 利用股票指数简化投资组合模型利用

50、股票指数简化投资组合模型 例例3.83.8 继续考虑例继续考虑例3.5（期望收益率仍定为（期望收益率仍定为15%）。在实际）。在实际的股票市场上，一般存在成千上万的股票，这时计算两两的股票市场上，一般存在成千上万的股票，这时计算两两之间的相关性（协方差矩阵）将是一件非常费时甚至不可之间的相关性（协方差矩阵）将是一件非常费时甚至不可能的事情。例如，能的事情。例如，1000只股票就需要计算（只股票就需要计算（1000 2）=499500个协方差。能否通过一定方式避免协方差的计算，个协方差。能否通过一定方式避免协方差的计算，对模型进行简化呢？例如，例对模型进行简化呢？例如，例3.5中还给出了当时股票

51、指中还给出了当时股票指数的信息，但我们到此为止一直没有利用，我们这一节就数的信息，但我们到此为止一直没有利用，我们这一节就考虑利用股票指数队前面的模型进行修改和简化。考虑利用股票指数队前面的模型进行修改和简化。 问题分析问题分析 可以认为股票指数反映的是股票市场的大势信息，对具体可以认为股票指数反映的是股票市场的大势信息，对具体每只股票的涨跌通常是有显著影响的。我们这里最简单化每只股票的涨跌通常是有显著影响的。我们这里最简单化地假设每只股票的收益与股票指数呈线性关系，从而可以地假设每只股票的收益与股票指数呈线性关系，从而可以通过先行回归方法找出这个线性关系。通过先行回归方法找出这个线性关系。

52、3.2.5 其他目标下的投资组合模型其他目标下的投资组合模型 目前介绍的模型中都是在可能获得的收益的数学期望满足目前介绍的模型中都是在可能获得的收益的数学期望满足一定最低要求的前提下，用可能获得的收益的方差来衡量一定最低要求的前提下，用可能获得的收益的方差来衡量投资风险，将其作为最小化目标。这种做法的合理性通常投资风险，将其作为最小化目标。这种做法的合理性通常至少要有两个基本假设：至少要有两个基本假设： 可能获得的收益的分布是对称的（如正态分布）。因为这可能获得的收益的分布是对称的（如正态分布）。因为这是未来收益高于设定的最低要求的机会和地域设定的最低是未来收益高于设定的最低要求的机会和地域设

53、定的最低要求的数量（高多少、低多少）和概率是一样的。可惜的要求的数量（高多少、低多少）和概率是一样的。可惜的是，实际中这个假设往往难以验证。是，实际中这个假设往往难以验证。 投资者对风险（或偏好）的效用函数是二次的。否则为什投资者对风险（或偏好）的效用函数是二次的。否则为什么只选择收益（随机变量）的二阶矩（方差）来衡量风险么只选择收益（随机变量）的二阶矩（方差）来衡量风险是之最小化，而不采取其他阶数的矩？是之最小化，而不采取其他阶数的矩？ 一般来说，投资者实际关心的通常是未来收益地一般来说，投资者实际关心的通常是未来收益地域设定的最低要求的数量（即低多少）和概率，域设定的最低要求的数量（即低多

54、少）和概率，也就是说关心的是下侧风险（也就是说关心的是下侧风险（downside risk）。）。所以，如果分布不是对陈的，则采用收益的方差所以，如果分布不是对陈的，则采用收益的方差来衡量投资风险就不一定合适。为了克服这个缺来衡量投资风险就不一定合适。为了克服这个缺陷，可以用收益低于最低要求的数量的均值（一陷，可以用收益低于最低要求的数量的均值（一阶矩）作为下侧风险的衡量依据，即作为最小化阶矩）作为下侧风险的衡量依据，即作为最小化的目标。此外，也可以采用收益低于最低要求的的目标。此外，也可以采用收益低于最低要求的数量的二阶矩（即收益的半方差，数量的二阶矩（即收益的半方差，semivarianc

55、e）作为衡量投资风险的依据。其实，半方差计算与作为衡量投资风险的依据。其实，半方差计算与方差计算类似，只是只有当收益低于最低要求的方差计算类似，只是只有当收益低于最低要求的收益率时，才把两者之差的平方计入风险，而对收益率时，才把两者之差的平方计入风险，而对收益高于最低要求的收益率时的数据忽略不计。收益高于最低要求的收益率时的数据忽略不计。这方面的具体模型这里就不再详细介绍了。这方面的具体模型这里就不再详细介绍了。 下面介绍一个与上面这些优化目标完全不同的投资组合模下面介绍一个与上面这些优化目标完全不同的投资组合模型，这个模型虽然很简单，但却会产生一些非常有趣的现型，这个模型虽然很简单，但却会产

56、生一些非常有趣的现象。象。 例例39 假设市场上只有两只股票假设市场上只有两只股票A、B可供某个投资者购可供某个投资者购买，且该投资者对未来一年的股票市场进行了仔细分析，买，且该投资者对未来一年的股票市场进行了仔细分析，认为市场只能出现两种可能的情况（认为市场只能出现两种可能的情况（1和和2）。此外，该投）。此外，该投资者对每种情况出现的概率、每种情况出现时两只股票的资者对每种情况出现的概率、每种情况出现时两只股票的增值情况都进行了预测和分析（见表增值情况都进行了预测和分析（见表3-7，可以看出股票，可以看出股票A、B的均值和方差都是一样的）。该投资者是一位非常保守的均值和方差都是一样的）。该

57、投资者是一位非常保守的投资人，其投资目标是使两种情况下最小的收益最大化的投资人，其投资目标是使两种情况下最小的收益最大化（也就是说，不管未来发生哪种情况，他都能至少获得这（也就是说，不管未来发生哪种情况，他都能至少获得这个收益）。如何建立模型和求解？个收益）。如何建立模型和求解？情形情形发生概率发生概率股票股票A A股票股票B B1 10.80.81.01.01.21.22 20.20.21.51.50.70.73.3 市场营销问题市场营销问题 3.3.1 新产品的市场预测新产品的市场预测 例例3.10 某公司开发了一种新产品，打算与目前市场上已某公司开发了一种新产品，打算与目前市场上已有的三

58、种同类产品竞争，为了了解这种新产品在市场上的有的三种同类产品竞争，为了了解这种新产品在市场上的竞争力，在大规模投放市场前，公司营销部门进行了广泛竞争力，在大规模投放市场前，公司营销部门进行了广泛的市场调查，得到了表的市场调查，得到了表3-8。四种产品分别记为。四种产品分别记为A、B、C、D，其中，其中A为新产品，表中的数据的含义是：最近购买某为新产品，表中的数据的含义是：最近购买某种产品（用行表示）的顾客下次购买四种产品的机会（概种产品（用行表示）的顾客下次购买四种产品的机会（概率）。例如：表中第一行数据表示当前购买产品率）。例如：表中第一行数据表示当前购买产品A的顾客，的顾客，下次购买产品下

59、次购买产品A、B、C、D的概率分别为的概率分别为75%、10%、5%、10%。请你根据这个调查结果，分析新产品。请你根据这个调查结果，分析新产品A未来的未来的市场份额大概是多少？市场份额大概是多少？产品产品A AB BC CD DA A0.750.750.10.10.050.050.10.1B B0.40.40.20.20.10.10.30.3C C0.10.10.20.20.40.40.30.3D D0.20.20.20.20.30.30.30.3 问题分析问题分析 新产品进入市场后，初期的市场份额将会不断发生变化，新产品进入市场后，初期的市场份额将会不断发生变化，因此，本例中的问题是一个离

60、散动态随机过程，也就是马因此，本例中的问题是一个离散动态随机过程，也就是马氏链（氏链（Markov chain）。很明显，上面给出的表实际上）。很明显，上面给出的表实际上是转移概率矩阵（注意每行元素的和肯定为是转移概率矩阵（注意每行元素的和肯定为1）。要分析）。要分析新产品新产品A未来的市场份额，就是要计算稳定状态下每种产未来的市场份额，就是要计算稳定状态下每种产品的概率。品的概率。 模型建立模型建立 记记N为产品种数，产品编号为为产品种数，产品编号为i（i=1，2，N），转移概），转移概率矩阵的元素记为率矩阵的元素记为Tij，稳定状态下产品，稳定状态下产品i的市场份额记为的市场份额记为pi。