地震波的计算程序出错请教

跨越冬天

这是我用newmark-beta方法写的动力方程求解位移加速的的程序，可是有错误，请教高手帮忙看看，谢谢

[ 本帖最后由 eight 于 2007-9-28 20:13 编辑 ]

eight

原帖由 跨越冬天 于 2007-9-28 19:45 发表

                               
登录/注册后可看大图

这是我用newmark-beta方法写的动力方程求解位移加速的的程序，可是有错误，请教高手帮忙看看，谢谢

永远不要指望别人通过运行你的程序来猜测你遇到的问题

花如月

2点建议：

（1）把数据存到一个文件里，然后读取。这样会使程序简练些
（2）程序写得有些乱，定义了一堆参数有些很烦琐。建议在程序上再多下功夫

里边的K没有定义时，会有下边的错误提示：

??? Reference to a cleared variable K.
Error in ==> wavefx at 22
Keff = c0*M + c1*C+K;% 形成等效刚度矩阵  

具体的意义我也搞不清楚。当加上K=0.5（我随意指定的），后程序没有报错，调试通过。

[ 本帖最后由 花如月 于 2007-9-28 20:55 编辑 ]