非负矩阵在一类代数方程的应用毕业论文.docVIP

非负矩阵在一类代数方程中的应用 摘 要：　本文运用非负矩阵理论中的Perron- Frobenius定理,讨论了在计算数学中大量出现的一类代数方程的根分布情况. 　Perron -Frobenius定理代数方程根分布定理 形如, , j = 1, 2, …, n, ( 1) 的代数方程作为许多差分方程的特征方程,在计算数学中大量出现因而,对(1)的根分布情况,特别是其最大模的根进行讨论是很有意义的.本文利用非负矩阵理论中著名的Perron -Frobenius定理,得出了(1)的根分布定理.该定理较全面地反映了其根的分布情况,它将关于(1)的结果统一在一起,并且有所推广此外,本文还给出了(1)的最大模正根的二个估计式最后还运用本文的结果,讨论了常微分方程数值解法中的一类含参数的线性多步法的零稳定性 . （2） 其中,, j = 1, 2, …, n,则称为方程(1)的友阵显然这里的是非负矩阵由(2)所定义的n阶矩阵A的特征方程是(1) . 通过行列式计算，整理后结果显然是：即为(1) . 定义2　 设为n阶不可约非负矩阵,又令h为的 指导教师：杨 芳 作者简介：李小鹏(1985—)，男，陕西陇县人. 的特征值的个数若h=1,则称为本原阵;若h 1,则称为具有指标h的循环阵我们不加证明地叙述如下的Perron-Frobenius定理(证明见[2]或[3])　设为n阶不可约非负矩阵,则1）的谱半径是的特征值2）对应于的特征向量x 0,3）是的简单特征值）当的任一元素增加时,也增加若为具有指标h( h 1)的n阶不可约非负矩阵则的模为的h个特征是，（）,j=0,1, … ,h-1.也就是说，它们均匀地分布在以原点为圆心为半径的圆周上（证明从略）. 定理3 设++…+，，…，0，1 …=n ,为不可约非负矩阵的特征多项式，设A的循环指标为h,则h是，，…，的最大公因子，即. 2 主要结果 定理4 对方程(1)设其n个根分别为记最大根模为p ,即p = m ax则p是(1)唯一的正单根,所有其它根的模均不超过（1）不为0，1…=n，即（1）++…+=0, （4） h=gcd，则当h=1时，p是（1）p. 沿用的记号，当h1时，p仍是（1）（1）h个模等于p的根：，j=0,1, …,h-1. （）n-h个根的模均小于h. 方程(1)的系数增加时,则其最大正单根p也增加 　首先证明(1)的友阵是不可约非负矩阵.…时，的直接图显然是强连接的，故为不可约非负矩阵.若在中有些不为0，则在A的直接图中增加了连通路径，这显然不消弱直接图的强连接性，所以是不可约非负矩阵. 下面证明定理中的各个结论: 由引理1知，方程（1）n个根等于（2）的n个特征值，所以.于是，由定理1中的4）可得出，又由定理1中的1）可知p是（1）p的唯一性证明如下：将（1）式中的n次多项式除以定义为上的函数，即 : （6） R上，与的零点相同，又因 ， （7）在单调增，故是唯一使=0，亦即使（1）式成立的正实数，于是得证，而和则可由定理2，定理3直接得出，其中所涉及到的p的唯一性可同上推证. 定理5 对于方程（1）的最大正单根，有如下估计式: ， （）. （9），易证 ， （10） ， 将此结果应用与（1）的友阵A便可得证. 推论1在（1）中，当时，有p=1. 3 一个应用 在常微分方程的数值解法中, 有一类含参数的线性多步法 （11）其中sR是一个可调整的实参数 [6]中曾且不会破坏零稳本文指出:在零稳定性条件的限制下方法(11) 的绝对稳定区域有时不可能无限增大(例如,当k = 4时)　 含实参数s的线性多步法(11)为零稳定的,其充分必要条件是s满足 （1） （13）的所有根均在单位圆内或圆上，而在单位圆上的根只能是单根，仍记p= ,下面按照s的值进行讨论： 1）若s1,则有1,知p1. 2）若0s1,则 （i）当s=0，，根为0（k-1重根）和1（单重根）； （ii）当 s = 1,,根为,j=0,1, …k-1,均为单根；