用于求解休假排队网络稳态分布的预处理GMRES法.pdf

第32卷 第3期 工 程 数 学 学 报 Vo1．32No．3 2015年06fl CHINESEJOURNALOFENGINEERING MATHEMATICS June2015 doi：10．3969／j．issn．1005—3085．2015．03．008 文章编号：1005—3085(2015)03—0391—06 用于求解休假排队网络稳态分布的预处理GMRES法术 杨淑伶 (广东工业大学应用数学学院，广州 510090) 摘 要：休假排队网络的稳态分布的计算在众多应用领域具有非常重要的意义 ，稳态分布向量 可通过求解一个奇异的线性系统来得到．但由于该线性系统维数 巨大且结构复杂， 不易直接求解，本文提出了一种带分块下三角预处理算子的GMRES迭代法求解该系 统．该预处理GMRES法具有易于构造且快速收敛的优点．数值试验 的结果验证了该 算法的优越性． 关键词：排队网络；休假；QBD过程；预处理GMRES法 分类号：AMS(2000)60K25；65F10 中图分类号：0226；0242 文献标识码：A 1 引言 休假排队在工程系统、管理系统和计算机系统中具有广泛的应用．已有的主要研究工 作是采用矩阵几何解[]等理论方法研究单个的排队系统并得到稳态分布各项指标的解析表 达式，但对于某些无穷小生成元非常复杂的排队模型，无法采用理论分析方法来研究其稳 态分布的各项指标 ．因而，从数值求解角度研究排 队模型具有极其重要的意义[2I3】．后一 领域 已有的工作主要局限于研究无休假Markovian排队网络4【15】，但很少涉及休假排队网 络的研究．前者的无穷小生成元通常包含三对角矩阵，而后者的则包含分块三对角矩阵， 因而数值计算也更为困难．文献f61研究了单个M ／M ／2队伍的休假排队系统的瞬时状态 分布．在本文中，我们将建立一个新的包含两个M／M／2队伍的溢流休假排队网络模型， 其中单个队伍的休假规则和文献 6【，71相同．和文献f4，61研究瞬时状态分布不同，本文感 兴趣的是该排队网络模型的稳定状态分布．需要特别指出的是，由于溢流排队网络中允许 顾客在不同队伍之间的转移，导致排队网络的Kolmogorov平衡方程 中系数矩阵更大型且 结构更复杂，从而快速有效的算法研究显得更为必要．本文基于系数矩阵的结构，提出了 一 种带分块下三角预处理算子的GMRES算法求解休假排队网络的稳态分布 向量，且通过 几个数值实例论证 了该算法的快速收敛性． 2 休假排队网络的稳定状态模型 本节先简要地描述本文研究的休假排队网络模型．假设两个队伍的状态空间都是有 收稿 日期：2013—10—24． 作者简介：杨淑伶 (1976年llfl生)，女，硕士，讲师．研究方向：数值计算 基金项目：国家 自然科学基金青年基金． 392 工 程 数 学 学 报 第32卷 限的且相互独立．在第 个队伍中(= 1，2)，系统空间容量为他 ，顾客到达率遵循参数 为入k的泊松分布 ．在每个队伍 中各有两名服务员，休假规则如下 ：当队伍为空时，服务 员立即同时开始一次休假，休假时问长度 服从参数为0的指数分布 ．休假结束时，如 果队伍中至少有一名顾客在等待服务，则服务员开始工作，否则继续下一次同样长度的 休假 ．第 个队伍中，服务员的休假次数不超过 次．在第 次休假结束后，无论队伍 中是否有顾客，服务员都处于工作状态．两名服务员按照先到先服务的原则给顾客提供 服务，服务时间服从指数分布，参数分别假设为 1和 ％2．当队伍 中新到一名顾客而两 名服务员都处于空闲状态时，由服务员1提供服务 ．假设两个队伍之间遵循如下溢流原 NI8】：当第 一个队伍排满而第二个队伍未满时，顾客可以从第一个队伍转到第二个队伍 中，但不允许相反方 向的转移 ． 设 (t)=i表示在时刻t第k个队伍中的顾客数， (t)=J表示该系统服务员在时 刻t的状态．当i=0，1，… ，礼％且J=0，1，… ， 一1时，状态 (i，J)表示系统中有i名