有理矩阵有理相似对角化计算机实现.docVIP

有理矩阵有理相似对角化计算机实现 　　摘要： 人工计算有理矩阵能否在有理数域上相似对角化是非常困难的，所以需要计算机来辅助实现，然而已有的数学软件对此问题的计算结果却存在着误差，于是需要研究有理矩阵在有理数域上相似对角化的算法及程序.因此在直接进行分数运算的基础上，首先使用矩阵的幂与来计算有理矩阵的特征多项式，克服了直接计算行列式的方法所存在的算法设计上的困难，其次根据有理多项式有理根的求法计算出有理矩阵的有理特征根，进而精确地计算出相应的有理特征向量，从而成功设计出判断及实现有理矩阵在有理数域上对角化的算法及相应的语言程序，使用该程序能够精确地解决有理矩阵在有理数域上相似对角化的问题。 关键词： 有理矩阵；有理相似对角化；算法；程序 中图分类号：O151.21 文献标识码：A 文章编号：1006-4311（2013）16-0194-06 0 引言 矩阵的相似对角化是重要的矩阵方法，然而人工计算却力不能及，所以需要计算机辅助实现。有理数域是常用的数域，有理矩阵更是常用的矩阵，所以判定及求解有理矩阵在有理数域上相似对角化问题非常必要，然而遗憾的是尽管已有的数学软件（如Matlab、Mathematica、Maple）能够解决矩阵对角化问题，但却只是在实数域、复数域上进行，而对于有理矩阵在有理数域上的相似对角化问题的解答存在着误差。 与相应的对角矩阵为 而在Maple中计算，却输出另一结果：可逆矩阵为 与相应的对角矩阵为 显然Mathematica与Maple的输出结果均存在着误差。 除此之外，这些数学软件还存在系统庞大、使用不便、输出的结果不直观等弱点。 上述问题使我们不能不考虑研究有理矩阵在有理数域上相似对角化问题的算法，设计出能够精确解决该问题的专用程序，本文阐述我们为此所做的研究工作。 1 相关概念及理论依据 定义1[1] 定义2[1] 设A是数域F上的n阶矩阵，称fA（x）的根λ为矩阵A的特征根。 定义3 设A是数域F上的n阶矩阵，如果存在F上的可逆矩阵P与F上的对角形矩阵D，使得P-1AP=D，那么即称A在数域F上可以对角化。 定理1 定理2 定理3[1] 设A是数域F上一个n阶矩阵，A可以对角化的充要条件是： ①A的特征多项式的根都在F内。 ②对于A的每一特征根λ，秩（λ-A）=n-s， 这里s是λ的重数。 综合除法[1]：设f（x）=a0xn+a1xn-1+…+an-1x+an，用x-c除f（x）所得的商式为q（x）=b0xn-1+b1xn-2+…+bn-2x+bn-1。 且余式为r，那么a0=b0，b1=cb0+a1，b2=cb1+a2，…，bn-1=cbn-2+an-1，r=cbn-1+an。 为了方便，我们把有理数域上的矩阵称为有理矩阵。有理矩阵A在有理数域上的特征根叫作A的有理特征根。如果有理矩阵A在有理数域上可以相似对角化，那么即称A能有理相似对角化。 2 算法设计 2.1 算法设计思想 首先，之所以Matlab、Mathematica和Maple在计算有理矩阵有理相似对角化问题中产生了误差，是因为这些软件的四则运算是建立在实数基础上。因此要精确地解决有理相似对角化问题，就必须在分数运算的基础上设计算法及程序。但因为计算机程序语言系统的四则运算都定义在实数域上，所以算法必须设计分数运算功能，以此为基础进行相关处理。关于分数的存储问题，采用两个矩阵分别存储分子与分母较为方便。 其次，定理3给出了判定数域F上的n阶方阵在F上能否相似对角化的一般方法（如图1所示），因此，对于有理矩阵A，判定、求解其能否有理相似对角化应做好以下主要的工作： 2.1.1 求出A的特征多项式fA（x）=xI-A.此项工作所存在的困难是：因为xI-A是含有未知数x的行列式，如果按照行列式的计算方法直接计算，算法将非常复杂，程序设计也将非常困难.然而定理1为我们提供了不展开行列式的计算方法，据此方法求A的特征多项式，需做以下处理： ①计算Ai（i=2…n）及Ai的迹。 ②生成矩阵 ③求矩阵A1的逆矩阵A1-1。 ④计算A1-1A2即可得到fA（x）的各项系数。 如此处理即避免了直接计算行列式所带来的算法设计上的困难，为设计有理矩阵有理对角化算法开拓了最关键的一步。 2.1.2 判断A特征根是否都是有理数。对此工作，表面地看需要求出A的特征多项式的所有根，然而这是不可能的。不过根据定理2及综合除法原理，我们却能够求出A的所有的有理特征根及其重数，再判断其所有有理特征根的重数和是否等于n即可解决这一问题。要实现此设想，可分以下四步进行：①变