参考文献:基于 Lagrange 运动方程的悬臂梁的固有频率分析,黄永玉,赵永刚,赵伟东,不过该文献中应变能公式应把r改为2r-1,已经M矩阵最后一个数错误,其他没有大问题