求解时变高阶微分方程组

xjzuo · 发表于 2007-6-28 20:31

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

octopussheng · 发表于 2007-6-28 20:33

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

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

vib · 发表于 2007-9-6 16:59

原帖由 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 · 发表于 2007-9-7 08:40

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

gcy1234567 · 发表于 2009-9-3 11:53

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

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 求解时变高阶微分方程组

回复 #19 vib 的帖子

评分