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

第 5卷 增刊 非 线 性 动 力 学 学报 1998年 6月 JoamalofNonlinearDynamicshiScienceandTechnologv VJaonl159S981 矩阵法计算NormalForm的关键I61题* 吴志强 胡海岩 (南京航空航天大学振动工程研究所 南京，z1Ix1J6) 摘要 用矩阵法求常微系统非半单分岔的NormalForm时，需要求解一系列的线性符 号方程组，其中，非共振项系致的方程组的求解问题己经解决:共振项系数的方程组，还未 见一般的解法，因为这种方程组的系数矩阵是奇异的.要求解这种方程组，会面临关键问硬 是如何确定其系数矩阵的满秩子矩阵.解决这一问题是本文的主要 目的 关键词 NormalF。朋，非半单分岔，最大秩 中图分类号 夕22- I 前言 在计算非线性常微分方程组NormalForm的方法中，矩阵法是一种最基本、最直观 的方法 由于它只涉及到线性代数运算，而一般的计算机代数语言都提供了这种功能， 因而在计算机推导方面显得更具潜力. 如大家所知，矩阵法最大的缺点是，随着问题的复杂性增加，要求解的线性符号代 数方程组的维数会迅速增加，很快会超出人手算的能力.用现有的方法，即便是用计算 机也难以胜任对非共振项来说，相应的代数方程组的系数矩阵是可逆的，它的求解不 成问题但对共振项来说，相应的代数方程组的系数矩阵是不可逆的，它的求解会面临 困难. 要想对共振项的方程组进行求解，首先要消元降维，为此必须先找出其系数矩阵的 最大的满秩子矩阵，然而就笔者所知，到目前为止还没有这样的方法. 文[1]提出的求解NormalF。二的直接方法，由于对非线性函数向量的级数表达，采 用了适当的记法，并在此基础上对非线性项进行了适当分类，使得定位符号矩阵的满秩 子矩阵有了可能.文[1]中采用的非线性函数向量的级数表达为 Su()=Zg，um， 1.122 其中9.是系数向量， “=令1us…u}, 切 = IM m2 ’“ m. u脚=u 卜.u笋u乒。 本文以文1【]工作为基础。探讨符号矩阵满秩子矩阵的定位问题.第二节简述文[1]中 的方法，第三节，一方面，巧妙地将祸合方程组系数矩阵满秩子矩阵的定位问题，转化 为一常数矩阵最大秩子矩阵的确定，另一方面，给出了常数矩阵最大秩子矩阵的定位方 法，从而解决了利用计算机简化该方程组遇到的关键问题.第四节举例说明本方法的求 解NormalForm的步骤和特点，因受篇幅的限制而忽略第五节是小结. .本文于 1998年3月30日收到 .‘‘‘苗 ‘. ‘ 非 线性 动 力 学 学 报 1999年 2 方法简介 考虑一般的高维非线性常微系统 s=Ax+刀刁 (1) 其中二。R,fx():It--)R，则存在变换 _ x=Tou+H(u) (2) 将(1)化成NormalForm形式: u=J,。十C(u) (3) 其中J。是由A的所有no个零实部特征值组成的Jordan标准形矩阵，兀是由零实部特 征值对应的特征向量(当特征值为非半单时，包括广义特征向且)组成的nx凡子矩阵，将 (2)(3)代入(