关于多体系统动力学的问题（传递矩阵法）
最近在学习南理工芮筱亭教授的《多体系统传递矩阵法及其应用》，在学习了第8章内容后，按照书上的思路和算法，根据例8.1提供的例子（平面三摆）用matlab编制了程序。但是在运行的时候，求出总的传递矩阵后，再求两端的状态模量时，提示矩阵是奇异的，运算结果将不准确。调了好几个参数后还是这样。请问这是什么原因造成的？因为程序调试一直没有调通，心中很是苦闷焦急。所以我附上了自己编的例8.1（平面三摆）的程序，希望各位达人能够在百忙中帮我分析程序错误的原因，给予纠正。不胜感激，谢谢！！！
顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶顶
页:
[1]