计划的合理性、适应性及对计划的管理效率，直接决定着型号研制项目计划的生命力。网络计划技术（PERT/CPM、GERT）在实际工程项目的管理中虽已取得较好的效果，但是也存在一些问题，如建立网络模型困难、繁琐，模型与实际情况差异大；在分析方法上，常用的参数计算方法虽然简便，但是由于型号研制工作的不确定性，结果往往难以反映实际情况，用处很少；同时对重大活动不能以灵活的调控方式给予充裕的时间资源，很难进行计划优化。而理想的计划应当能够反映并代表模型的随机性的统计规律性，具有一般性，应能兼顾和协调计划活动的随机性、重要程度等具体要求，合理均衡地分配时间资源，否则就会出现计划不得不在研制工作中不断修改，失去计划的指导作用。因此有必要研究一种有效的计划分析和优化方法。

 

网络计划技术是以网络图或条线图来描述系统元素的多边联系的系统表达和分析模型，尽管可以反映各项活动之间的次序关系，但是很难与型号研制项目的产品系统（或项目系统）的树状分枝（并列）、分层（嵌套）结构形式相对应，而工作分解结构（Work Breakdown Structure, WBS）确定了型号研制项目产品内在的分支分层结构体系， 反映了要完成产品所需进行的所有工作及其相互关系，适合于作为研制部门以产品为中心编制研制计划的基准模型。因此采用WBS活动树作为计划总体模型，可以使计划模型与工程模型结构一致，建模容易，便于理解，使用方便。在分析方法中，应用Monte Carlo随机模拟法，对时间参数进行动态、预测性的模拟，特别符合型号研制项目工作的不确定性和随机性。之后，用相对偏差和最小二乘法建立寻优判据，随机投点直接搜索寻优，逐步优选，最终实现型号研制项目计划的优化。

 

（1）自下向上计算：

设起始活动i最早开始时间ESi已经确定，且已知持续时间为t_i, 则起始活动的最早完成时间

EF_i=ES_i+ t_i             1)

对非起始活动j，最早开始时间

ES_j=max ( EF_i ),   i∈I_j   2)

I_j为活动j的紧前活动集合。而非起始活动j的最早完成时间

EF_j=ES_j+ t_j              3)

对顶层活动k按3)式计算出的结果，即为完成计划全部活动的最早时间，即工期T。

（2）自上向下计算：

对顶层活动k，令最迟完成时间

LK_k=EK_k=T                    4)

最迟开始时间

LS_k= LF_k—t_k= Es_k                         5)

对非顶层活动j

LF_j= LS_{f ,}                                   6)

 LS_j= LF_j—t_j                              7)

f为活动j的紧后活动节点。

（3）时差计算

对非顶层活动j，时差

FF_j= ES_f  —EF_j                           8)

f仍为活动j的紧后活动节点。

对各活动i, 总时差

TF_i= LF_i—EF_i= LS_i—ES_i               9)

当TF_i=0时，则说明活动i是关键活动。相邻关键活动构成的自底向顶的路径就是关键路线。

 

（4）完成概率计算

当活动节点层次数目较多（大于10层）时，按中心极限定理，工期T近似服从正态分布，其数学期望

T_e=∑t_j       j∈J                          10)

方差

S_e²=∑σ_j²       j∈J                        11)

J为某关键路线活动集合。若指定工期T_zd，则

服从正态分布，利用标准正态分布表可以得到对于指定工期T_zd的指定完成概率P_zd。

在工程实际工作中，WBS活动节点层次数目很少大于10层，且由于研制工作的不确定性和随机性，所以仅仅采用上述式10)、11)计算的T_e及S_e²有较大的误差，必须进行改进。

 

（1）Monte Carlo随机模拟分析

采用Monte Carlo随机模拟法，可以准确地得出模型总体的随机特性。用MC方法模拟运行第k次（k=1，2，3，……   N）时，随机抽样得出第k个计划

{ t_ik }, i =1，2，3，……   n,

n为活动树节点数，t_ik 为活动i 的第k次时间抽样值。

按前述方法计算活动时间参数，寻找关键路线，计算研制工期，得出工期这一随机变量的一次取值

T_k= EF_ok         13)

EF_ok 为顶层活动在第k次模拟运行中最早完成时间。对于N次模拟运行结果

{T_i },  i=1，2，3，……   N          14)

计算样本均值

 

 

T_s 与S_s即为工期T的数学期望与方差。对模拟运行结果14)进行统计分析，对分布参数做各种点估计、区间估计和峰度、偏度分析，判别T的分布类型，有利于全面了解工期T的规律性，为决策提供重要参考。

（2）工期分布函数

对模拟运行结果14) 适当划分时间区间，按区间统计频数，可得工期T的近似分布密度函数：

P = f ( T )    17)

按区间统计频数, 可得工期T的累积分布函数：

P = F ( T )     18)

还可以得到经验分布函数

P = F_N ( T )    19)

这三个分布函数可用计算机处理成表或曲线，通过查阅，可以观察了解工期T的分布形态，并得到工期T与完工概率P的对应关系。

 

（1）按计划优化的目标要求确定指定工期T_zd或指定概率P_zd，利于上述模拟分析得到的19）式，由T_zd确定P_zd，或由P_zd确定T_zd 。

（2）给出计划工期T_jh的允许偏差极限E。

（3）进行寻优模拟。

首先按研制工作要求确定活动持续时间t和起始活动最早开始时间t_q的随机分布概率取值区间 [ p_a，p_b ]，按t和t_q的累积分布函数确定t和t_q的抽样区间，此区间小于t和t_q的随机分布区间。

其次，对t和t_q按其概率分布在上述限定的抽样区间抽样，得到t和t_q的抽样值，模拟活动按概率要求随机实现。

最后计算活动节点的各时间参数ES、EF、LS、LF、FF、TF等，统计活动的关键度D。

（4）按工期允许的偏差要求优选计划。

在N₁个模拟所得工期中，选出接近指定工期T_zd的N₂（N₂ < N₁）个较优计划工期。

 

（5）按最小二乘法进一步优选计划。

先依据累积分布函数 P = F ( T )和指定概率P_zd，得出活动持续时间指定值t_zd和起始活动最早开始时间指定值t_qz。对N₂个较优计划工期，分别用对应的模拟抽样起始点再次模拟工程运行，产生模拟数据，随机抽取出较优计划的活动持续时间t和起始活动最早开始时间t_q 。再按最小二乘法，以指定值为标准值，以抽样模拟值为计划值，建立寻优目标

 

 

Q为起始活动集合。寻优目标的意义是使活动时间计划值对指定值的相对偏差平方和达到最小。计算N₂个较优计划的δ，选出偏差较小的N₃（N₃ < N₂）个较优计划，对其记录。

最后将优选的N₃个较优计划按偏差值做排列，分别用对应的模拟抽样点再次进行抽样，再次模拟工程运行，产生模拟数据，抽样出t和t_q，计算活动节点的各时间参数ES、EF、LS、LF、FF、TF等，得到较优计划的全套数据，依此输出能满足总体和局部要求的N₃个较优计划

在上述过程中，把计划目标范围由N₁缩小到N₂再缩小到N₃，还可以依此类模拟推到更小更具体的范围，最终选出在本世纪内允许偏差极限内的相对最优的计划。

 

为更好地说明上述模拟、优化方法过程，下面给出一个简单的例子，对工期进行计算和优化。某项目的WBS活动信息数据如表1所示。

表1  WBS活动信息数据

编码	活动	a	m	b	Pz	ao	mo	bo	Pqz
1000	惯导系统	7.0	10.0	12.0	0.80
1100	惯导部件	5.0	7.0	9.0	0.80
1110	惯性平台	5.0	8.0	10.0	0.80
1111	陀螺	14.0	18.0	22.0	0.70	1.0	3.0	5.0	0.60
1112	加速度计	12.0	15.0	20.0	0.70	1.0	2.0	3.0	0.70
1113	导电环	11.0	14.0	19.0	0.80	2.0	4.0	7.0	0.60
1114	力矩电机	6.0	8.0	10.0	0.70
1120	电源系统	4.0	5.0	6.0	0.75
1121	滤波系统	2.0	3.0	5.0	0.75	3.0	4.0	6.0	0.75
1122	启动电路	2.0	3.0	4.0	0.70	4.0	7.0	9.0	0.70
1123	AC/DC电源	4.0	5.0	7.0	0.70
1200	惯导托架	5.0	6.0	8.0	0.70	5.0	6.0	7.0	0.60
1210	结构体	1.0	2.0	3.0	0.70	4.0	5.0	6.0	0.60
1211	托架底座	2.0	3.0	4.0	0.80
1220	密封胶圈	1.0	2.0	4.0	0.80
1300	卫星接收系统	4.0	8.0	12.0	0.55
1310	接收天线	5.0	7.0	10.0	0.70	2.0	4.0	5.0	0.55
1320	专用馈线	6.0	8.0	12.0	0.70	2.0	4.0	5.0	0.65

 

给出t和t_q的抽样区间对应的概率区间下限指定值P_z和P_qz，上限值为1。活动时间t和起始活动最早开始时间设定为三角分布，给出3个参数值a(乐观), m（可能）, b（悲观）。当指定工程完工率P_zd=0.8或指定工期T_zd=47.5月时，经过上述随机模拟（100次）和搜索寻优，得到一个较优进度计划如表2所示。

表2  进度计划表