第005章 矩阵特征值及其计算.ppt

解：先由柔度矩阵求逆得到刚度矩阵（只是为了说明问题）。根据初始迭代向量的取法，经计算后所取的结果为 1 0 1 0 0 0 1 0 0 1 0 0 1 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 1 0 1 1 0 0 0 0 ?取6个初始迭代向量 例5.2 用子空间迭代法求例5.1系统的前三阶模态。 前6阶固有频率的迭代值 9.740774e+001 1.558175e+003 7.879525e+003 2.481355e+004 6.000000e+004 1.215997e+005 前6阶固有频率的准确值 9.740774e+001 1.558175e+003 7.879525e+003 2.481355e+004 6.000000e+004 1.215997e+005 很高的精度！ 迭代5次得： 迭代得到的前3阶振型 1阶 2阶 3阶 将梁分成20段时振型的迭代计算结果：振型已趋于光滑（迭代6次） 《振动力学》讲义 第5章 矩阵特征值及其计算 第5章 矩阵特征值及其计算 5.1 特征值问题 5.1.1 一般特征值问题 5.1.2 实对称矩阵的特征值问题 5.1.3 Hermitian矩阵的特征值问题 5.1.4 两个实对称矩阵形成的特征值问题 5.1.5 非对称实矩阵的特征值问题 5.2 特征解的性质 5.2.1 Rayleigh原理 5.2.2 约束系统的Rayleigh原理 5.2.3 特征值的极大极小性质 5.2.4 Hermitian矩阵的正定准则 5.2.5 Gerschgorin定理 5.3 特征解的计算方法 5.3.1 无阻尼系统特征值问题的Ritz减缩方法 5.3.2 无阻尼系统特征值问题的矩阵迭代法 5.3.3 无阻尼系统特征值问题的子空间迭代法 5.3.4 矩阵的变换和分解 5.3.5 实对称矩阵特征值问题的Jacobi法 5.3.6 一般矩阵特征值问题的QR法 解：我们将梁分成10段，共11个节点，将每段的质量均分到节点上，再将每个梁段视为无质量弹性梁。因此梁弯曲振动坐标为 y = [ y1, y2, …, y9]T 设梁的弯曲振型为 …… 0 1 9 10 l x y(x) 系统的质量、柔度矩阵分别为 例5.1 用Rayleigh法求均匀简支梁的基频估值。梁长l，质量m，抗弯刚度EI。 根据柔度计算方法和简支梁的挠度公式可得柔度系数为 这是柔度矩阵的上半三角矩阵的元素，下三角矩阵的元素可由柔度矩阵的对称性得到。 柔度矩阵为刚度矩阵的逆，系统的无阻尼振动方程，用柔度矩阵表示为： 上式两边同乘 M ，可得对称形式的特征值问题为 由上式得基频的估值为 所以，Rayleigh商的表达式为 基频的精确值为 即使只是将梁分成 2 段（很粗糙），算得基频的估值为 也有相当高的精度。