计算QBD衰减率的数值方法-计算数学专业论文.docxVIP

中文摘要摘要 中文摘要 摘要 我们在文中针对不可约非周期正常返且有限相位的QBD，提出了几种求 解其衰减率77的数值算法。其中改进的二分法是在原始的二分法基础上，根 据QBD中随机非负矩阵的基本性质设计出来的，牛顿法和割线法的提出和改进 则是建立在原始方程在所求区间内的根具有良好特性的基础之上，这些算法在 速度上比原始二分法具有很大优势。而二次特征值法则从另一个角度即求解矩 阵特征值问题来计算衰减率77。最后我们给出了一些数值实验结果来说明这些结 论。 关键词：改进的二分法，牛顿法，割线法，二次特征值法 中图法分类号：0241．6 一II— 英文摘要Abstract 英文摘要 Abstract In the thesis．we design some numerical methods fbr the caudal characteristic叩of the QBD，which is irreducible，aperiodic，positive recurrent and with finite phase．These methods include the improved bisection method，the Newton method，the cutting plane method and the quadratic eigenvalue method．Numerical experiments show that these methods ale superior to the original bisection method in convergence and computational time． Key Words： improved bisection，Newton method，cutting plane method，quadratic eigenvalue method Chinese Library Classification Number：0241．6 引言引言 引言 引言 自从上个世纪60年代，R．V Evans【2】和V Wallace【18】提出了不可约非 周期正常返QBD的稳态分布服从矩阵几何分布71n=丌0彤以来，这套理论在 过去几十年里经历了蓬勃地发展。为了揭示不可约非周期正常返QBD的稳态 分布，人们开始对随机矩阵兄进行了一系列的研究，随着研究的深入，人们 发现翮拘谱半径比R本身在刻画QBD的行为特征上起着更重要的作用，后来人 们称之为QBD的衰减率r#(Caudal Characteristic)。但无论是V Wallace【1 8】和M． E Neuts【13】设计的计算翮拘线性收敛算法还是G．Latouche和V Ramaswami 【9】设计的计算R的二阶收敛算法，都因为迭代步数过多或计算量过大等原 因而变得不够实用，于是如何绕开计算R而直接求得衰减率卵成为研究QBD领 域一个重要的课题。直到后来M．E Neuts【12】证明了衰减率叩满足代数方 程x(z)=z，0z1，这为快速求解77提供了一个重要的理论依据，于是G． Latouche和v．Ramaswami【8】8提出使用二分法来求解该方程。 本文主要的创新之处在于改进了原始的二分法，并设计了收敛速度更快的 牛顿法和割线法，以及从另外一个角度通过求解二次特征值问题来计算?7。通过 数值例子可以看出这些算法在速度上比原始的二分法更有优势。 本文的第一章介绍-]QBD暑11其衰减率77的基本概念和性质，第二章和第三 章介绍了方程)((z)=z根的基本性质和已有的求解该方程的二分法，而第四、 五、六章则是本文设计的算法的核心部分，第七章为数值例子，对各个算法进 行了一个总体的评价和比较。 一IV— 第一章预备知识第一章 第一章预备知识 第一章 预备知识 1．1 QBD简介 当我们研究一个二维状态空间{(佗，i)：n≥0，1≤i冬m)上的离散时 f司Markov链{xt，t∈Ⅳ)时，常常对其状态空间划分为U礼o t(／2)，其中f(n)= {(佗，1)，(钆，2)， ，(／2，m))，这里n≥0。我们称第一项坐标n为“层”，有时 我们也用“层”表示整个子集“n)，每一层中的状态数m可以有限也可以是无 穷。 如果一个Markov链的一步状态转移只能限制在同一层或相邻的两层，即一 步状态转移／口N台13曼1己发生在(几，歹)到(佗7，J，)，这里n7=n，佗+1．Sign一1，／2≥1。我们 称该Markov链为拟生灭过程，简称QBD。 特殊地，当该QBD的转移概率P【x1=(佗7，歹)Ix0=(佗，i)】仅依赖i，歹和n7一 佗，而不依赖于特定的层n和n7时，我们称其为时齐的拟生灭过程。此时它的转 移矩阵为分块三对角阵