以悬臂梁为例，用有限元法编程生成质量矩阵和刚度矩阵（采用平面梁单元，单元数在10到30之间）。
使用matlab的用于广义特征值问题分析的函数（eig(A,B,'chol')，该命令使用乔莱斯基分解的方法，计算对称矩阵A和对称正定矩阵B的广义特征值问题。A对应刚度矩阵，B对应质量矩阵）计算梁的前三阶固有频率和振型，计算得到梁的前三阶固有频率和振型，并画出振型图。(Taking cantilever beam as an example, the mass matrix and stiffness matrix are programmed by finite element method (using planar beam elements, with units ranging from 10 to 30).　　　　　　Using Matlab's function (eig (A, B,'chol') for the generalized eigenvalue problem analysis, the command uses the Jo Lessky decomposition method to calculate the generalized eigenvalue problem of the symmetric matrix A and the symmetric positive definite matrix B. A corresponds to the stiffness matrix and B corresponding to the mass matrix. The first three natural frequencies and modes of the beam are calculated, and the first three natural frequencies and modes of the beams are calculated, and the vibration pattern is drawn.)