人力资源调度的优化模型

1、人力资源调度的优化模型摘要本文主要研究人力资源调度的最优化问题。人力资源调度问题中所要处理的数据之间的关系是 比较繁琐的，所以如何有效地设置决策变量，找出相互关系是我们建立模型的突破口。上述模型属 于多元函数的条件极值问题的围，然而许多实际问题归结出的这种形式的优化模型，起决策变量个 数 n 和约束条件 m 一般比较大，并且最优解往往在可行域的边界上取到，这样就不能简单地用微分 法求解，数学规划是解决这类问题的有效方法。根据所给的“ PE 公司”技术人员结构及工资情况表、不同项目和各种人员的收费标准表格，为 了在满足客户对专业技术人员结构要求的前提下，使“PE公司”每天的直接收益最大，我们首先

2、对不同项目的不同技术人员的分配个数进行假设，从而得到了“ PE 公司”每天总收入 I 和每天总支出 C ，所以每天的直接收益 U I C ，这就是公司每天直接收益的目标函数。在此基础上我们建立 了基于 Matlab 软件上的线性规划方法一和基于 Lindo6.0 软件上的整数线性规划方法二来求解这个 模型。首先我们 Matlab 软件运行这个函数，得到求得的值恰好是整数，满足题意，在题目的约束条 件下得到的最大公司效益是 27150 元，此时的人员分布如下表所示：项目技术人员ABCD高级工程师1521工程师6362助理工程师2521技术员1310因为对题中的数据稍做改动时得出的答案就会出现小数

3、的现象，为了更好的解决该问题，我们 又引入了一个很好地能处理整数的软件 Lindo6.0 ，得到了各个有效的数据。并在模型扩展中运用已 建立的程序对所得的结果进行灵敏度分析，即讨论在收费标准不变的情况下技术人员结构对公司收 益的影响以及在技术人员结构不变的情况下收费标准对公司收益的影响，并且进一步分析在怎样的 围最优解保持不变，并联系社会实际进行了一定的分析。最后在适当简化模型的同时，对模型进行 了改进和推广，预示了高素质人才在现代社会中将发挥着越来越重要的作用。关键词 ：人力资源调度；决策变量；可行域；灵敏度分析；博弈论. 问题重述：“PE公司”是一家从事电力工程技术的中美合资公司， 现有

4、41 个专业技术人员， 其结构和相应的工资水平分布如表 1 所示。表 1 公司的人员结构及工资情况高级工程师工程师助理工程师技术员人数917105日工资（元）250200170110目前，公司承接有 4个工程项目，其中 2项是现场施工监理，分别在 A地和 B地， 主要工作在现场完成；另外 2项是工程设计，分别在 C地和 D地，主要工作在办公室完 成。由于 4 个项目来源于不同客户，并且工作的难易程度不一，因此，各项目的合同对 有关技术人员的收费标准不同，具体情况如表 2 所示。表 2 不同项目和各种人员的收费标准高级工程师工程师助理工程师技术员A1000800600500收费B15008007

5、00600（元/ 天）C1300900700400D1000800700500为了保证工程质量，各项目中必须保证专业人员结构符合客户的要求，具体情况如 表 3 所示：表 3 ：各项目对专业技术人员结构的要求ABCD高级工程师 工程师 助理工程师 技术员 总计13221 1025223162221 111228 118说明：表中“ 13”表示“大于等于 1，小于等于 3”，其他有“”符号的同理；项目 D，由于技术要求较高， 人员配备必须是助理工程师以上， 技术员不能参加； 高级工程师相对稀缺， 而且是质量保证的关键， 因此，各项目客户对高级工程师 的配备有不能少于一定数目的限制。 各项目对其他专

6、业人员也有不同的限制或要 求；各项目客户对总人数都有限制；由于 C、D两项目是在办公室完成，所以每人每天有 50 元的管理费开支。由 于 收 费 是 按 人 工 计 算 的 ， 而 且 4 个 项 目 总 共 同 时 最 多 需 要 的 人 数 是 10+16+11+18=55，多于公司现有人数 41。因此需解决的问题是： 如何合理的分配现 有的技术力量，使公司每天的直接收益最大？并写出相应的论证报告。二 . 模型假设和符号说明 :1. 模型假设 根据问题的要求，并为了达到将问题进一步明确抽象的目的，在我们的模型中有如 下的假设：1) PE公司是在同一时间接手 A、B、C、D 四个工程项目。2

7、) PE公司的各种技术人员分工相当明确，如高级工程师不能兼职工程师的工作。3) PE公司的专业技术人员在接手工程期间不存在着请假，缺席的现象。4) 在 PE公司接手这四个工程的间段， 市场物价稳定。 各级技术人员的收费标准分 别在如下的围：高级工程师的收费最低不低于 800 元，最高不超过 1600 元； 工程师的收费最低不低于 500 元，最高不超过 1200 元； 助理工程师的收费最低不低于 300 元，最高不超过 900 元； 技术员的收费最低不低于 150 元，最高不超过 600 元。5) 假设 A，B，C， X 是如下四个矩阵： a11, a12 , a13 , a14 , a21

8、, a22 , a 23 , a 24 ,a31 , a32 , a33 , a34 , a41 , a42 , a43 , a44 b11 ,b12 ,b13 ,b14 ,b21 ,b22 , b23 , b24 , b31 , b32 , b33 , b34 ,b41,b42,b43,b44 C c11 ,c12 , c13 , c14 , c21 , c22 ,c23,c24,c31,c32,c33,c34,c41,c42 ,c43, c44 X x11 , x12 , x13 , x14 , x21, x22,x23, x24, x31, x32 ,x33, x34 , x41, x4

9、2, x43, x44T2. 符号说明x11, x12 , x13 , x14 分别代表的是在 A、B、C、D项目中高级工程师的人数安排x21,x22,x23,x24分别代表的是在 A、B、C、D项目中工程师的人数安排x31, x32, x33 , x34分别代表的是在 A、B、C、D项目中助理工程师的人数安排x41,x42,x43,x44分别代表的是在 A、B、C、D项目术员的人数安排a11 , a12 , a13 , a14分别代表的是高级工程师在 A、B、C、D项目中的每天收费标准a21,a22,a23 ,a24分别代表的是工程师在 A、B、C、D项目中的每天收费标准a31 , a32

10、, a33 ,a34 分别代表的是助理工程师在 A、B、C、D项目中的每天收费标准a41,a42,a43 ,a44分别代表的是技术员在 A、B、C、D项目中的每天收费标准b11 , b12 ,b13 ,b14 分别代表的是 PE公司为 A、B、C、D项目中的高级工程师每天所支付 的费用b21,b22,b23,b24分别代表的是 PE公司在 A、B、C、D项目中的工程师每天所支付的 费用 b31,b32,b33,b34分别代表的是 PE公司为 A、B、C、D项目中的助理工程师每天所支付 的费用b41,b42,b43,b44 分别代表的是 PE公司为 A、B、C、D项目中的技术员每天所支付的费 用

11、注： PE 公司支出费用包括技术人员的工资和 C、D项目中每个人员每天的 50 元管 理费。三. 模型的建立：模型：基于 Matlab 的线性规划方法根据题意以及上面的符号说明可以得到下列 A， B， C的值A=1000 1500 1300 1000 800 800 900 800 600 700 700 700 500 600 400 500B=250 250 300 300 200 200 250 250 170 170 220 220 110 110 160 160C=AB=750 1250 1000 700 600 600 650 550 430 530 480 480 390 490

12、 240 340 于是得到目标函数：4MaxZ CX(aij bij ) xiji,j 1s.t.x11x12x13x149x21x22x23x2417x31x32x33x3410x41x42x43x445x11x21x31x4110x12x22x32x4216x13x23x33x4311x14x24x34x4418x113x125x132x142x24812x11 , x14 , x34 , x41 , x 43x12 , x13 , x 21, x 22 , x 23 , x24 , x31 , x32 , x 33x42 x44我们首先来观察表 1 和表 3，因为 A、B、C、D四个工程

13、需要的技术员最低限分别是 1、3、1、0，而 PE公司的技术员恰好只有 5 人，所以关于技术员的调度就已经确定， A工程 1人，B程 3人，工程 1人，工程 0人。即： x41 1,x42 3,x43 1,x44 0。又有 C工程中高级工程师的数量已定， x13 2 ，因此其实只有 11 个决策变量在影响最终公司 效益。用 Matlab 6.5 软件中的函数 x,fval=linprog(C,A,b,Aeq,beq,vlb,vub)解 , 具体程序看附录 1：得出数据 X=1,5,2,1,6,3,6,2,2,5,2,1,1,3,1,0TMax Z=CX=27150通过对高级工程师的人数变化时

14、( 其他因素全都不变 ), 分析其对最大公司效益的影 响, 分别取高级工程师的人数是 9,10,11,12 的情况 :得出结果如下表所示高级工程师人数x11x12x13x14.91.00005.00002.00001.0000101.70205.00002.00001.2980112.56485.00002.00001.4352123.00005.00002.00002.0000x21x22x23x24,96.00003.00006.00002.0000105.29803.10866.00002.5933114.43523.85066.00002.7143124.00004.01876.000

15、02.9813高级工程师人数x31x32x33x3492.00005.00002.00001.0000102.00004.89142.00001.1086112.00004.14942.00001.8506122.00003.98132.00002.0187高级工程师人数x41x42x43x4491.00003.00001.00000.0000101.00003.00001.00000.0000111.00003.00001.00000.0000121.00003.00001.00000.0000上述表格缺陷在于人员分配个数出现了小数 , 这跟实际问题相违背。 分析其原因主要 在于：我们用这个

16、 Matlab6.5 软件做出来的就是基于用单纯形法引入松弛变量而得出来 的。因为松弛问题是作为一个线形规划问题，其可行解的集合是一个凸集，任意两个可 行解的凸集组合仍为可行解。由于整数规划问题的可行解一定也是松弛问题的可行解 (反之则不一定)，所以前者最优解的目标函数值不会优于后者最优解的目标函数值。在一般情况下，松弛问题的最优解不会刚好满足变量的整数约束条件，因而不是整 数规划的可 行解，自然就 不是整 数规划的最优 解。我们用 Matlab 6.5 中函数 x,fval=linprog(C,A,b,Aeq,beq,vlb,vub) 求解出来的当高级工程师人数变化时出现 了小数现象，就是上

17、述所述的问题。此时，若对松弛问题的这个最优解中不符合整数要 求的分量简单的取整，所得到的解不一定是整数规划问题的最优解，甚至也不一定是整 数规划问题的可行解。基于这个问题我们引入了解这个模型的第二种方法。 模型 :基于 LinDo6.0 的整数规划方法：对于该问题我们有同于 4.1 的目标函数及约束条件 , 即MaxZCXi, j(aij1bijs.t.x11x12x13x149x21x22x23x2417x31x32x33x3410x41x42x43x445x11x21x31x4110x12x22x32x4216x13x23x33x4311x14x24x34x4418x113x125x132

18、x142x24841x11 , x14 , x34 , x41 , x43x12 , x13 , x 21 , x22 , x23 , x24 , x31 , x32 , x33) xijx42x44用 Lindo6.0 求解问题 , 程序具体见附录 2:得出数据 X=1,5,2,1,6,3,6,2,2,5,2,1,1,3,1,0 TMaxZ=CX=27150 结果跟方法一的相同，从而也验证了结果的正确性。 下面着重通过人数变化（其他因素都不变）对公司最大效益的影响进行分析，得出下列 数据：高级工程师人数变化总人数公司效益最大值 （单位为 元）9412715010422785011432855

19、0124429250134529250工程师人数变化总人数公司效益最大值 （ 单位为 元）174127150184227700194328250204428800214529350224629900234730450244831000254931550265032100275132150285232150助理工程师的人数变化总人数公司效益最大值 （ 单位为 元）1041271501142276301243281101344285901445290701546295501647300301748305101849309901950314702051319502152324302253329102

20、35433390245533870技术员人数的变化总人数公司效益最大值 （ 单位为 元）5412715064227590743280308442847094528910104629250114729590124829930134930270145030410155130410四. 模型分析：主要采用灵敏度分析法 :上面表格中除了告诉我们问题的最优解和最优值以外， 还能挖掘出许多隐含着的有用 信息：1. （在允许围）每增加一个高级工程师使得公司效益增加 700 元，如高级工程师的人 数从 9 个增加到 10个公司效益就增大了 700 元。增加一个工程师使得公司效益增加 550 元，增加一个助理工

21、程师使得公司效益增加 480 元，增加一个技术员使得公司效益增加 440 元，上面公司效益的增加可以看作人数的潜在价值称为 “影子价格”，即高级工程师 的影子价格是 700 元，工程师的影子价格是 550 元，助理工程师的影子价格是 480 元， 技术员的影子价格是 440 元。2. 当目标函数的系数发生变化时（假定约束条件不变） ，最优解会改变吗？带着这个问 题，我们又用 Lindo6.0 软件进行了单个系数变化的处理 . 即在最初的最优解不变的情况 下，求出各系数允许的变化围 （以其他系数不变作为前提， 求出一个目标系数变化的围） 。 D公司不需要技术员， C 公司对高级工程师的人数是常数

22、 2，所以他们不会对最终公司效 益产生影响，这里就只分析其他的自变量前的系数的变化。列出表格如下：此变 量前 的系 数X11X12X13X14X21X22X23X24允许 的变 化围0-1249751- +0- +0-1249551- +551-650551- +0-599公司 效益 变化 幅度15216362此变量X31X32X33X34X41X42X43前的系 数允许的 变化围0-529480- +0-5800-5300- +0- +0- +公司效 益变化 幅度2521131说明：表格的第一行表示列出的各变量 xij 前系数（表格中就用变量代替表示） 。 第二行表示在最优解不变的情况下其系

23、数的变化围。如 0-1249 是指其前面的系数从 0 变化到 1249 元，最优解都是不会变化的。第三行表示当变量前面的系数在这个围变化时，虽然最优解不发生变化，但是最优 值将发生变化，目标系数变化一个单位时，最优值以表格中列的数值为幅度变化。而还 可以发现这个幅度就是得出的最优解里的相应项目中具体人员分布的人数， 我们得出的 最优解是 X=1,5,2,1,6,3,6,2,2,5,2,1,1,3,1 。由上述表格还可以画出最大公司效益 与目标系数的关系图，这里只列举 x11，x12,x 13，x22 前面目标系数变化时，最大公司效 益和目标系数的变化关系。同理其他的都可以归结为这四类中的一类。

24、如下图：红线表 示的就是上面表格中列的最优解保持不变时，系数变动对最优值的影响区域。对上面图形说明： 各个图中的黑点 （即折线的交点） 表示最优解发生变化的临界点。 用这个分析结果很容易看到， 若某一个工程项目只增加其中一种技术人员的收费标 准时（该工程中其他技术人员的收费标准和其它项目中人员的收费标准，人数需求均保 持不变），可以使公司的最大效益增加，但是人力资源调度不变。这样公司可以根据这 个标准跟对方商谈价格， 在一定的人数围， 公司的领导阶层可以尽可能地提高收费标准， 为公司获得最大的效益； 同时对需求方来说他则要出低一点的收费标准来得到项目需求 的人数，这样需求方可以使支出大大减小。

25、这个尽可能低的工资从理论上来说有的甚至 可以达到零值，但是这是不符合实际情况的。主要原因是我们算出来的是理论值，从理 论上来说是可以的。我们不能说只根据理论而去采取使高级工程师的工资为250 元（固定工资），而技术员的工资还是 500 元这样的收费标准，即使是人员分配还是符合项目 的要求。这就要求公司和需求方共同商议来达到一个较好的值，所以模型要跟实际联系 起来才能得到更好的实际意义。3. 对上面的第 1点的人员分配做进一步的分析：（1）. 对于高级工程师来说，随着高级工程师单位人数的增加，最大公司效益以 700 元的幅度递增。但是当人数增加到 12 人时，最大公司效益却不再增加。此时在达到这

26、 个最大公司效益的上限值时， 对应着一个总人数是 4（4 55）。为什么会出现这个现象呢？ 按照正常的想法应该是高级工程师越多越好。但在此题中对此很好解释，因为高级工程 师的人数都有一个上限， 即他不能像我们想象的那样可以无限增加直到 55。又由于实际 中高级工程师相对稀缺，而且是质量保证的关键。因此各个项目客户对高级工程师的配 备有不能少于一定数目的限制是需要的。（2）. 对于工程师来说，同样随着工程师单位人数的增加，最大公司效益以 550 元的 幅度递增。当人数达到 27人时，最大公司效益也不再增加，而是保持在 27150 元这个 值不变。此时的总人数是 51（55），又出现了这样的现象。

27、显然这里并不是纯粹的如上 面的原因，仔细观察上面的数据图可以看到此时 A、B、C三个项目都已经达到了最大的 约束项限制，而在 D地却没有达到，从这里就可以得出主要原因是在 D 地，同样的在 D 地对工程师的人数需 28，又是上面提到的存在上限的问题。（3）. 对于助理工程师来说，随着助理工程师单位人数的增加，最大公司效益以小一 点的值 480 元的幅度递增。当人数达到 24 人时，已经达到了使最大公司效益不变的值， 此时最大值 33870元，此时总人数刚好达到题目限制的最大值 55 人，得到了较好的人 数分配，使得每个工程项目都能有足够的人力资源来进行工作。（4）. 对于技术员来说， 相对于技

28、术员单位人数的增加， 最大公司效益也是以幅度 440 元来递增。虽然技术员并没有像上面所说有上限的限制，但是在最大公司效益不再改变 时，技术员有 14 人，此时总人数是 50 人，并不能达到 55 人。此时最大公司效益已经 达到了 30410 元。这里虽然没有像上面所说的上限约束的限制，可是人数还是不能达到 最大分配。这里还是要从题目表格里仔细观察发现： D 地由于技术要求较高，技术员不 能参加。这里就给了人数一个很大的限制，所以技术员在是A，B，C三地达到最大满足后不能再增加了。这个分析结果对实际应用很有价值。像这里 A，B，C，D四个项目共需要 55 个人， 当公司有足够的人员来分配给他们

29、时，如果它不加考虑的就分配这样会出现人员浪费。 比如就拿技术员来说，他达到 14 人时公司效益不再增加。公司如果分配了 18人给需求 方，这样这四个技术员就好似没有得到任何利润，就存在着这种浪费现象。实际中公司 接手的项目往往有很多个，这里人员分配就显得更加重要，否则很可能会出现分配了这 里而在那里得不到满足。可见考虑这个问题是非常必须的。用这个模型可以很方便的解 决这些问题。五. 模型的改进与推广：（1）以上的模型还没有讲到各级技术人员分别对公司效益的产生影响。实际情况 中，不同级别的技术人员对公司的效益都是不同的，下面我们就以问题为代表来讨论高 级工程师、工程师、助理工程师、 技术员在整体

30、上和人均上每天为公司所创的效益。 （注： 在这里我们假设在 C、D 两个项目中每人每天不需要 50 元的管理开支费。）那么，在整 体上：9 个高级工程师每天为公司所创的纯收入为：1000 1 1500 5 1300 2 1000 1 250 9 9850（元）17 个工程师每天为公司所创的纯收入为：800 6 800 3 900 6 800 2 200 17 10800（元）10 个助理工程师每天为公司所创的纯收入为：600 2 700 5 700 2 700 1 170 10 510（0 元）5 个技术员每天为公司所创的纯收入为：500 1 600 3 400 1 500 0 110 5 2

31、150（元）各级技术人员整体上每天为公司所创的纯收入（表1）在人均上：每个高级工程师每天为公司所创的纯收入： 9850 9 1094.（4 元） 每个工程师每天为公司所创的纯收入： 10800 17 635.（3 元） 每个助理工程师每天为公司所创的纯收入： 5100 10 51（0 元） 每个技术员每天为公司所创的纯收入： 2150 5 430（元）各级技术人员人均上每天为公司所创的纯收入（表 2）从上面表 1 中可以看出高级工程师和工程师在公司中所起的作用是举足轻重的，由 此可见一个公司若没有高技术人员的加盟，它的效益就不会高，在激烈的竞争中就很难 立足，有可能就会面临破产的悲剧。从表 2

32、 中可以看出高级工程师每天人均为公司所创的纯收入是十分可观的， 正因为如此 现在无论各个行业各个部门都希望高素质人才加盟自己的队伍， 而国际上也把科学技术 作为衡量一个国家综合国力的重要标志，这也印证了“科学技术是生产力”的道理，符 合当今社会现实潮流。六 . 模型评价 ：模型的优点如下：1）模型的主体采取 L 软件处理数据和对其进行灵敏度分析，准确性高，容量大，逻辑 性严格，计算速度快，具有较强的说服力和适应能力。2）动态的分析了各种人员人数变化对公司效益的影响和各种人员收费标准变化对公司 效益的影响。3）从单纯的问题分析中，预见到了现今社会对高技术人才的需求程度。 模型的缺点如下：1） 我们在灵敏度分析中，对模型中最优值的影响因素只是从单个方面的变化考虑。不是十分的全面参考文献1. 胡运权，郭耀煌 . 运筹学教程 . 清华大学 . 19982. 启源 , 金星, 叶