求解时变高阶微分方程组

xjzuo

如果问题有进展,欢迎来论坛交流你的想法.
我们论坛相当欢迎解决问题后回贴的人,并会适当加分.

octopussheng

这是肯定的，有结果马上和大家分享的

我现在怎么说也算一个忠实的振动论坛的”泡友“了， 呵呵

vib

原帖由 octopussheng 于 2007-6-28 09:41 发表

                               
登录/注册后可看大图

请问xjzuo，按理说这个两自由度二阶微分方程要写成x'=Ax格式，微分方程应该写成
dy=[dx1,dx2,dx3,dx4]
dx1=x2;
dx2=-x2-x1-x1*cot(t)-x3*sin(t)+cos(t);
dx3=x4;
dx4=-x4-x3-x3*sin(t)-x1*sin(t)+si ...

初值一般都设为全0吗？

[ 本帖最后由 ChaChing 于 2010-4-10 23:49 编辑 ]

无水1324

这个不一定，可以根据实际系统调整，
0不行，全部换为1，或者更大，这个没有准则

gcy1234567

:@D 都是高手~学习了（正好在研究Matlab）