关于多体系统传递矩阵法的问题

dosking

程序8.1.rar (2.16 KB, 下载次数: 7)

2013-1-20 11:51 上传

点击文件名下载附件

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

木兵

提高计算精度，试试。