人力资源调度的优化模型

1、人力资源调度的优化模型 标准化文件发布号：（9456EUATWKMWUB-WUNN-INNULDDQTY-KII 人力资源调度的优化模型 摘要 本文主要硏究人力资源调度的最优化问题。人力资源调度问题中所要处理的数据 之间的关系是比较繁琐的，所以如何有效地设置决策变量，找出相互关系是我们建立 模型的突破口。上述模型属于多元函数的条件极值问题的范围，然而许多实际问题归 结出的这种形式的优化模型，起决策变量个数n和约束条件m般比较大，并且最优 解往往在可行域的边界上取到，这样就不能简单地用微分法求解，数学规划是解决这 类问题的有效方法。 根据所给的“PE公司”技术人员结构及工资情况表、不同项目和各种

2、人员的收费 标准表格，为了在满足客户对专业技术人员结构要求的前提下，使“PE公司”每天的 直接收益最大，我们首先对不同项口的不同技术人员的分配个数进行假设，从而得到 了 “PE公司”每天总收入I和每天总支出C,所以每天的直接收益U = I-C.这就是 公司每天直接收益的U标函数。在此基础上我们建立了基于Matlab软件上的线性规划 方法一和基于Lindo6.0软件上的整数线性规划方法二来求解这个模型。首先我们 Matlab软件运行这个函数，得到求得的值恰好是整数，满足题意，在题LI的约束条件 下得到的最大公司效益是27150元，此时的人员分布如下表所示： 技术人员 A B C D 高级工程师

3、1 5 2 1 工程师 6 3 6 2 助理工程师 2 5 2 1 技术员 1 3 1 0 因为对题中的数据稍做改动时得出的答案就会出现小数的现象.为了更好的解决 该问题，我们又引入了一个很好地能处理整数的软件Lindo6. 0，得到了各个有效的数 据n并在模型扩展中运用已建立的程序对所得的结果进行灵敏度分析，即讨论在收费 标准不变的情况下技术人员结构对公司收益的影响以及在技术人员结构不变的情况下 收费标准对公司收益的影响，并且进一步分析在怎样的范围内最优解保持不变，并联 系社会实际进行了一定的分析。最后在适当简化模型的同时对模型进行了改进和推 广预示了高素质人才在现代社会中将发挥着越来越重要

4、的作用。 关键词:人力资源调度；决策变量；可行域；灵敏度分析；博弈论 问题重述: “PE公司”是一家从事电力工程技术的中美合资公司，现有41个专业技术人员, 其结构和相应的工资水平分布如表1所示。 表1公司的人员结构及工资情况 高级工程师 工程师 助理工程师 技术员 人数 9 17 10 5 日工资（元） 250 200 170 110 目前，公司承接有4个工程项目，其中2项是现场施工监理.分别在A地和B 地，主要工作在现场完成；另外2项是工程设计，分别在C地和D地，主要工作在办 公室完成。由于4个项目来源于不同客户，并且工作的难易程度不一，因此，各项目 的合同对有关技术人员的收费标准不同，具

5、体情况如表2所示。 表2不同项目和各种人员的收费标准 高级工程师 工程师 助理工程师 技术员 A 1000 800 600 500 收费 B 1500 800 700 600 （元/ C 1300 900 700 400 天） D 1000 800 700 500 为了保证工程质量.各项目中必须保证专业人员结构符合客户的要求，具体情况 如表3所7F : 表3 ：各项目对专业技术人员结构的要求 A B C D 高级工程师 1 3 25 2 12 工程师 22 N2 M2 28 助理工程师 22 N2 22 Ml 技术员 21 N3 21 总计 W10 W16 W11 W18 说明： 表中T3表示

6、“大于等于1,小于等于3,其他有“符号的同理； 项目D,由于技术要求较高，人员配备必须是肋理工程师以上，技术员不能参 加； 高级工程师相对稀缺，而且是质量保证的关键，因此，各项目客户对高级工程 师的配备有不能少于一定数目的限制。各项目对其他专业人员也有不同的限制 或要求； 各项目客户对总人数都有限制； 由于C、D两项目是在办公室完成，所以每人每天有50元的管理费开支。 由于收费是按人工计算的，而且4个项目总共同时最多需要的人数是 10+16+11+18二55,多于公司现有人数41。因此需解决的问题是：如何合理的分配 现有的技术力量，使公司每天的直接收益最大？并写出相应的论证报告。 二模型假设和

7、符号说明： 1. 模型假设 根据问题的要求，并为了达到将问题进一步明确抽象的目的，在我们的模型中有如 下的假设： 1) PE公司是在同一时间接手A、B、C、D四个工程项目。 2) PE公司的各种技术人员分工相当明确，如高级工程师不能兼职工程师的工 作。 3) PE公司的专业技术人员在接手工程期间不存在着请假，缺席的现象。 4) 在PE公司接手这四个工程的间段，市场物价稳定。各级技术人员的收费标准 分别在如下的范围内： 高级工程师的收费最低不低于800元，最高不超过1600元； 工程师的收费最低不低于500元，最高不超过1200元； 助理工程师的收费最低不低于300元，最高不超过900元； 技术

8、员的收费最低不低于150元，最高不超过600元。 5) 假设A, B. C, X是如下四个矩阵： A 1 .】2，13 9 】4 9 219 22 9 d23。“24 31，。32，“33，“34 仃42 43 9 44 =孙 i 上2 上13，勺 4 “21 9 ”22 ”23 ”24 “31、32 33 “34 编,”42，43 “44 c = C I, C2 , C3 9 C|4 9 c2i、C22，C23，C24，C3，32、33，C34，C41，C42C43，C44 I T X -|, X2，13，14，-21 “22 *23, -24，*31，“32 “33 *34，人41 *42

9、 *43，人 44 2. 符号说明 xH,xI2,x13,x14分别代表的是在A、B、C、D项目中高级工程师的人数安排 心宀2，心3宀4分别代表的是在A、B、C、D项目中工程师的人数安排 切宀2宀，心分别代表的是在A、B、C、D项目中助理工程师的人数安排 心宀，心宀分别代表的是在A、B、C、D项目中技术员的人数安排 如如如，分别代表的是高级工程师在A、B、C、D项目中的每天收费标准 “21 ,“22, “23, “24分别代表的是工程师在A、B、C、D项目中的每天收费标准 5，也2卫33，34分别代表的是助理工程师在A、B、C、D项目中的每天收费标准 “41，“42, “43, “44分别代表

10、的是技术员在人、B、C、D项目中的每天收费标准 垢九九九分别代表的是PE公司为A、B、C、D项目中的高级工程师每天所支付 的费用爼血血上24分别代表的是PE公司在A、B、C、D项目中的工程师每天所支付 的费用编,血厶3,爼分别代表的是PE公司为A、B、C、D项目中的肋理工程师每天所 支付的费用 尽兀分别代表的是PE公司为A、B、C、D项目中的技术员每天所支付的费 用 注：PE公司支出费用包括技术人员的工资和C、D项目中每个人员每天的50元管 理费。 三模型的建立: 模型I :基于Matlab的线性规划方法 根据题意以及上面的符号说明可以得到下列A, B, C的值 A二1000 1500 130

11、0 1000 800 800 900 800 600 700 700 700 500 600 400 500 B二250 250 300 300 200 200 250 250 170 170 220 220 110 110 160 160 C二AB二750 1250 1000 700 600 600 650 550 430 530 480 480 390 490 240 340 于是得到目标函数： 4 MaxZ = CX =工（知-bjx 七 S.t. xH+xI2+xI3 + x14 9 x2l + x22 + x23 + x24 17 x31 + x32 + X33 + x34 10 些

12、1 + X42 + x43 +x44 5 AH+x21+x3l+A41 10 x12 + x22 + x32 +x42 16 A13 + X23 + “33 + 勺3I X14 + X24 + X34 +%44 18 Ai 1 3 召2 S 5 A13 - 2 m S 2 X24 8 1 XIPX|4, X34, X4l, X43 2 S X2 9 册3 , X21,兀22，“23，“24，“319 “32，“33 3 “2 ()=m 我们首先来观察表1和表3因为A、B、C、D四个工程需要的技术员最低限分别 是1、3、1、0,而PE公司的技术员恰好只有5人，所以关于技术员的调度就已经确 定A工

13、程1人，B程3人工程1人工程0人。即： 兀41 =1，心2 =3,心3 =1，44 = 又有C工程中高级工程师的数量已定,XJ3 = 2 ,因此其 实只有11个决策变量在影响最终公司效益。 用 Mat lab 6. 5 软件中的函数x, fval=l inprog (C, A, b, Aeq, beq, vlb, vub)解，具体 程序看附录 1 :得出数据 x=l, 5, 2, 1, 6, 3, 6, 2, 2, 5, 2, 1, 1, 3,1, 0T Max Z二CX二27150 通过对高级工程师的人数变化时（其他因素全都不变）,分析其对最大公司效益的影 响,分别取高级工程师的人数是9,

14、10, 11, 12的情况： 得出结果如下表所示: 高级工程师人数 Xu Xi： X13 X14. 9 1. 0000 5. 0000 2. 0000 1. 0000 10 1. 7020 5. 0000 2. 0000 1. 2980 11 2. 5648 5. 0000 2. 0000 1. 4352 12 3. 0000 5. 0000 2. 0000 2. 0000 X：1 X22 X23 X21. 9 6. 0000 3. 0000 6. 0000 2. 0000 10 5. 2980 3. 1086 6. 0000 2.5933 11 4. 4352 3. 8506 6. 0000

15、 2.7143 12 4.0000 4. 0187 6. 0000 2. 9813 高级工程师人数 X31 X3： XS3 Xm 9 2. 0000 5. 0000 2. 0000 1. 0000 10 2. 0000 4. 8914 2. 0000 1. 1086 11 2. 0000 4. 1494 2. 0000 1. 8506 12 2. 0000 3. 9813 2. 0000 2.0187 高级工程师人数 X.I1 X12 X43 X44 9 1. 0000 3. 0000 1. 0000 0. 0000 10 1. 0000 3. 0000 1. 0000 0. 0000 11

16、1. 0000 3. 0000 1. 0000 0. 0000 12 1. 0000 3. 0000 1. 0000 0. 0000 上述表格缺陷在于人员分配个数出现了小数，这跟实际问题相违背。分析其原因主 要在于： 我们用这个Matlab6. 5软件做出来的就是基于用单纯形法引入松弛变量而得出来 的。因为松弛问题是作为一个线形规划问题，其可行解的集合是一个凸集.任意两个 可行解的凸集组合仍为可行解，由于整数规划问题的可行解一定也是松弛问题的可行 解(反之则不一定)，所以前者最优解的目标函数值不会优于后者最优解的目标函数 值 在一般情况下，松弛问题的最优解不会刚好满足变量的整数约束条件，因而不

17、是 整数规划的可行解.自然就不是整数规划的最优解n我们用Mat lab 6. 5中函数 x, fval=linprog(C, A, b, Aeq, beq, vlb, vub)求解出来的当高级工程师人数变化时出现 了小数现象.就是上述所述的问题。此时，若对松弛问题的这个最优解中不符合整数 要求的分量简单的取整.所得到的解不一定是整数规划问题的最优解.甚至也不一定 是整数规划问题的可行解。基于这个问题我们引入了解这个模型的第二种方法。 模型II :基于LinDo6. 0的整数规划方法: 对于该问题我们有同于4. 1的目标函数及约束条件，即 MaxZ = CX =工(a- _%)xx- + X2

18、+13 + X4 9 *21 + 22+ 兀23 + 兀24 17 勺1 + X32 + 兀33 + X34 10 41 + 尤42 + “3 + “4 5 + x21 + X3| + x41 10 + V22 + X32 + X42 16 和 + X23 + 勺3 + 些3 11 州4 +兀4 +心4 +心 18 3 5 和 2 和 2 花4 8 1 a-1pa-14,x34,x41,j % 2X2l，X22i V23 3 X42 0 = *44 24 9勺1 *32, “33 用Lindo6. 0求解问题，程序具体见附录2: 得出数据 X二1, 5, 2, 1, 6, 3, 6, 2, 2

19、, 5, 2, 1, 1, 3, 1, 0T MaxZ=CX=27150 结果跟方法一的相同从而也验证了结果的正确性。 F面着重通过人数变化（其他因素都不变）对公司最大效益的影响进行分析得出下 列数据: 高级工程师人数变化 总人数 公司效益最大值（单位为 元） 9 41 27150 10 42 27850 11 43 28550 12 44 29250 13 45 29250 工程师人数变化 总人数 公司效益最大值（单位为 元） 17 41 27150 18 42 27700 19 43 28250 20 44 28800 21 45 29350 22 46 29900 23 47 30450

20、 24 48 31000 25 49 31550 26 50 32100 27 51 32150 28 52 32150 助理工程师的人数变化 总人数 公司效益最大值（单位为 元） 10 41 27150 11 42 27630 12 43 28110 13 44 28590 14 45 29070 15 46 29550 16 47 30030 17 48 30510 18 49 30990 19 50 31470 20 51 31950 21 52 32430 22 53 32910 23 54 33390 24 55 33870 技术员人数的变化 总人数 公司效益最大值（单位为 元） 5

21、 41 27150 6 42 27590 7 43 28030 8 44 28470 9 45 28910 10 46 29250 11 17 29590 12 48 29930 13 19 30270 14 50 30410 15 51 30410 四模型分析: 主要采用灵敏度分析法： 上面表格中除了告诉我们问题的最优解和最优值以外，还能挖掘出许多隐含着的有 用信息： 1. （在允许范围内）每增加一个高级工程师使得公司效益增加700元，如高级工程 师的人数从9个增加到10个公司效益就增大了 700元。增加一个工程师使得公司效益 増加550元，增加一个助理工程师使得公司效益增加480元，増加一

22、个技术员使得公 司效益增加440元为“影子价格，闫 高级工程师的影子价格是700元，工程师的影子价格是550元，助理工程师的影子价 格是480元，技术员的影子价格是440元。 2. 当目标函数的系数发生变化时（假定约束条件不变），最优解会改变吗？带着这个 问题.我们又用Lindo6. 0软件进行了单个系数变化的处理即在最初的最优解不变的 情况下，求出各系数允许的变化范围（以其他系数不变作为前提.求出一个目标系数 变化的范围）D公司不需要技术员，C公司对高级工程师的人数是常数2,所以他们 不会对最终公司效益产生影响，这里就只分析其他的自变量前的系数的变化。 列出表格如下： 变前系坟 此量的娄 n

23、 1 n n 薜围 n n n n n n n 公司 效益 变化 幅度 n 此变量 前的系 数 X： X：3 X： X： X.： X-13 允许的 变化范 围 0-529 480- +oo 0-580 0-530 0- +8 0- +8 0 - x 公司效 益变化 幅度 2 5 2 1 1 3 1 说明： 表格的第一行表示列出的各变量X前系数（表格中就用变量代替表示）。 第二行表示在最优解不变的情况下其系数的变化范围。如0-1249是指其前面的系数从 0变化到1249元，最优解都是不会变化的。 第三行表示当变量前面的系数在这个范围变化时，虽然最优解不发生变化，但是最 优值将发生变化，目标系数变

24、化一个单位时，最优值以表格中列的数值为幅度变化。 而还可以发现这个幅度就是得出的最优解里的相应项目中具体人员分布的人数，我们 得出的最优解是x二1, 5, 2, 1, 6, 3, 6, 2, 2, 5, 2, 1, 1, 3, 1。由上述表格还可以画出最大公 司效益与目标系数的关系图，这里只列举xn, x12, x13（ X22前面目标系数变化时，最 大公司效益和目标系数的变化关系。同理其他的都可以归结为这四类中的一类。如下 图：红线表示的就是上面表格中列的最优解保持不变时，系数变动对最优值的影响区 域。 对上面图形说明：各个图中的黑点（即折线的交点）表示最优解发生变化的临界 占 八、O 用这

25、个分析结果很容易看到，若某一个工程项目只増加其中一种技术人员的收费 标准时(该工程中其他技术人员的收费标准和其它项目中人员的收费标准，人数需求 均保持不变)，可以使公司的最大效益增加，但是人力资源调度不变。这样公司可以 根据这个标准跟对方商谈价格.在一定的人数范围内，公司的领导阶层可以尽可能地 提高收费标准，为公司获得最大的效益；同时对需求方来说他则要出低一点的收费标 准来得到项目需求的人数，这样需求方可以使支出大大减小。这个尽可能低的工资从 理论上来说有的甚至可以达到零值，但是这是不符合实际情况的。主要原因是我们算 出来的是理论值，从理论上来说是可以的。我们不能说只根据理论而去采取使高级工

26、程师的工资为250元(固定工资)，而技术员的工资还是500元这样的收费标准，即 使是人员分配还是符合项目的要求。这就要求公司和需求方共同商议来达到一个较好 的值，所以模型要跟实际联系起来才能得到更好的实际意义。 3. 对上面的第1点的人员分配做进一步的分析： (1) 对于高级工程师来说，随着高级工程师单位人数的增加，最大公司效益以700 元的幅度递增但是当人数増加到12人时，最大公司效益却不再增加。此时在达到这 个最大公司效益的上限值时，对应着一个总人数是44 (55)。为什么会出现这个现象 呢？按照正常的想法应该是高级工程师越多越好。但在此题中对此很好解释，因为高 级工程师的人数都有一个上限

27、.即他不能像我们想象的那样可以无限增加直到55。又 由于实际中高级工程师相对稀缺，而且是质量保证的关键n因此各个项目客户对高级 工程师的配备有不能少于一定数目的限制是需要的。 (2) .对于工程师来说，同样随着工程师单位人数的増加，最大公司效益以550元的 幅度递增，当人数达到27人时，最大公司效益也不再增加，而是保持在27150元这个 值不变。此时的总人数是51 (55),又出现了这样的现象。显然这里并不是纯粹的如 上面的原因.仔细观察上面的数据图可以看到此时A、B、C三个项目都已经达到了最 大的约束项限制，而在D地却没有达到，从这里就可以得出主要原因是在D地，同样 的在D地对工程师的人数需

28、求是2 8,又是上面提到的存在上限的问题。 (3) .对于肋理工程师来说，随着肋理工程师单位人数的增加，最大公司效益以小一 点的值480元的幅度递增。当人数达到24人时，已经达到了使最大公司效益不变的 值.此时最大值33870元，此时总人数刚好达到题目限制的最大值55人，得到了较好 的人数分配，使得每个工程项目都能有足够的人力资源来进行工作。 (4) .对于技术员来说，相对于技术员单位人数的增加，最大公司效益也是以幅度 440元来递増。虽然技术员并没有像上面所说有上限的限制，但是在最大公司效益不再 改变时，技术员有14人，此时总人数是50人，并不能达到55人。此时最大公司效益 已经达到了 30

29、410元。这里虽然没有像上面所说的上限约束的限制，可是人数还是不 能达到最大分配。这里还是要从题目表格里仔细观察发现：D地由于技术要求较高，技 术员不能参加。这里就给了人数一个很大的限制，所以技术员在是A, B, C三地达到 最大满足后不能再增加了。 这个分析结果对实际应用很有价值像这里A, B, C, D四个项目共需要55个 人，当公司有足够的人员来分配给他们时，如果它不加考虑的就分配这样会出现人员 浪费比如就拿技术员来说，他达到14人时公司效益不再增加。公司如果分配了 18 人给需求方.这样这四个技术员就好似没有得到任何利润，就存在着这种浪费现象。 实际中公司接手的项目往往有很多个，这里人

30、员分配就显得更加重要.否则很可能会 出现分配了这里而在那里得不到满足。可见考虑这个问题是非常必须的。用这个模型 可以很方便的解决这些问题， 五模型的改进与推广: （1）以上的模型还没有讲到各级技术人员分别对公司效益的产生影响。实际情况 中.不同级别的技术人员对公司的效益都是不同的，下面我们就以问题为代表来讨论 高级工程师、工程师、助理工程师、技术员在整体上和人均上每天为公司所创的效 益。（注：在这里我们假设在c、D两个项目中每人每天不需要50元的管理开支 费n ）那么，在整体上： 9个高级工程师每天为公司所创的纯收入为： 1000 x1 + 1500 x5 + 1300 x2 + 1000 x

31、1-250 x9=9850（元） 17个工程师每天为公司所创的纯收入为： 800 x6 + 800 x3 + 900 x6 + 800 x2-200 x17 = 10800（76） 10个助理工程师每天为公司所创的纯收入为： 600 x2 + 700 x5 +700 x2 + 700 x1 170 x10 = 5100（元） 5个技术员每天为公司所创的纯收入为： 500 x1+ 600 x3+ 400 x1+500 x0 110 x5 = 2150（元） 各级技术人员整体上每天为公司所创的纯收入 （表1） 15000f 高级工程师助理工程师 在人均上： 每个高级工程师每天为公司所创的纯收入：9

32、8504-91094.4（元） 每个工程师每天为公司所创的纯收入：10800+17 35.3（元） 每个肋理工程师每天为公司所创的纯收入：5100+10 = 510（元） 每个技术员每天为公司所创的纯收入：2150一5 = 430（元） 各级技术人员人均上每天为公司所创的纯收入（表2） 从上面表1中可以看出高级工程师和工程师在公司中所起的作用是举足轻重的，由 此可见个公司若没有高技术人员的加盟 它的效益就不会高，在激烈的竞争中就很 难立足，有可能就会面临破产的悲剧。 从表2中可以看出高级工程师每天人均为公司所创的纯收入是十分可观的，正因为如 此现在无论各个行业各个部门都希望高素质人才加盟自己的队伍.而国际上也把科学 技术作为衡量一个国家综合国力的重要标志，这也印证了“科学技术是生产力”的道理， 符合当今社会现实潮流。 六模型评价： 模型的优点如下： 1）模型的主体采取L软件处理数据和对其进行灵敏度分析，准确性高，容量大，逻辑 性严格，计算速度快，具有较强的说服力和适应能力。 2）动态的分析了各种人员人数变化对公司效益的影响和各种人员收费标准变化对公司 效益的影响。 3）从单纯的问题分析中，预见到了现今社会对高技术人才的需求程度。 模型的缺点如下： 1）我们在灵敏度