以悬臂梁为例，用有限元法编程生成质量矩阵和刚度矩阵（采用平面梁单元，单元数在10到30之间）。 使用matlab的用于广义特征值问题分析的函数（eig(A,B,'chol')，该命令使用乔莱斯基分解的方法，计算对称矩阵A和对称正定矩阵B的广义特征值问题。A对应刚度矩阵，B对应质量矩阵）计算梁的前三阶固有频率和振型，计算得到梁的前三阶固有频率和振型，并画出振型图。(Taking cantilever beam as an example,