矩阵法计算Normal+Form的关键问题.pdfVIP

第5誊增刊 非线性动力学学报 VoI5SJ 1999年6月 ofNonlinearinsci“cealld Jan．1998 生!些1 Dynamics Technology l 矩阵法计算NormaForm的关键问题+ 吴志强 胡海岩 (南京航空航天大学振动工程研究所南京．2Jfx)】6) 摘要 用矩阵涪求常微系统非半单分岔的NommlForm时．需要求解一系列的线性符 号方程组．其中．非菇振项系数的方程组的求解问题已经解决：共振项系数的方程组．还未 见一般的解法，因为这种方程组的系数矩阵是奇异的要求解这种方程组，会面临关键问题 是如何确定其系数矩阵的满秩子矩阵解决这一问题是本文的主要目的 关键词NormalForm，非半单分岔，最大秩 中囤分类号032．2-。- 1前言 在计算非线性常微分方程组NormalForm的方法中，矩阵法是一种最基本、晟直观 的方法．由于它只涉及到线性代数运算，而一般的计算机代数语言都提供了这种功能， 因而在计算机推导方面显得更具潜力． 如大家所知．矩阵法最大的缺点是，随着问题的复杂性增加．要求解的线性符号代 数方程组的维数会迅速增加，很快会超出人手算的能力。用现有的方法，即便是用计算 机也难以胜任．对非共振项来说，相应的代数方程组的系数矩阵是可逆的，它的求解不 成问题．但对共振项来说，相应的代数方程组的系数矩阵是不可逆的，它的求解会面临 困难． 要想对菸振项的方程组进行求解，首先要消元降维，为此必须先找出其系数矩阵的 最大的满秩子矩阵，然而就笔者所知，到目前为止还没有这样的方法． 文【l】提出的求解NormalForm的直接方法．由于对非线性函数向量的级数表达，采 用了适当的记法，并在此基础上对非线性项进行了适当分类，使得定位符号矩阵的满秩 子矩阵有了可能．文【l】中采用的非线性函数向量的级数表达为 删=∑g。u-， Imll2 其中耳．是系数向量， ：三‰．‰：一：．．”毫j ””=“，．“，．．．．．”，。 本文以文【l】工作为基础．探讨符号矩阵满秩子矩阵的定位问题．第二节简述文【l】中 的方法。第三节，一方面，巧妙地将耦合方程组系数矩阵满秩予矩阵的定位问题，转化 为一常数矩阵最大秩子矩阵的确定，另一方面．给出了常教矩阵最大秩子矩阵的定位方 法，从而解决了利用计算机简化该方程组遇到的关键问题第四节举例说明本方法的求 解NormalForm的步骤和特点，因受篇幅的限制而忽略．第五节是小结 。本文于1998年3月30日收到 210 非线性动力学学报 1998年 2方法简介 考虑一般的高维非线性常微系统 量=4】：+刷 (1) 其中X∈R“，／扛J：R”——÷JRl．则存在变换 x=Tou+H(11) ∞ Form形式： 将(1)化成Normal n=Jou+C俐 (3) 其中J。是由A的所有‰个零实部特征值组成的Jordan标准形矩阵，瓦是由零实部特 征值对应的特征向量(当特征值为非半单时，包括广义特征向量)组成的，，×no子矩阵．将 (2)、0)代入(1)．即得变换(2)的非线性项．阿倒及NormalForm(3)中非线性项C(u