求助！！matlab解微分方程组！（急）

bruce0601

二自由度振动的弹簧摆  设弹簧原长为l,进度系数为k,弹簧变形为x1,摆角为x2,
已知动力学方程 求解。。。。。。。。求解 得出图形
万分感谢！！！！！！！谢谢！！！！！！

matlab 数学软件求解微分方程组  望好心的大虾指点！！！ 再次感谢！！！！

yuxin

[这个贴子最后由yuxin在 2007/05/11 08:15am 第 3 次编辑]

"只保留x1,x2 ， 及其导数的二次项，化简"后的方程求解
用Mathematica写出：
Solve[{D[x1, {t, 2}] + omg10^2 x1 == l D[x2, t]^2 - g/2 x2^2,
    D[x2, {t, 2}] + omg20^2 x2 == -2/l D[x1, t] D[x2, t] +
        omg20^2/2 x1 x2}, {x1, x2}]
结果：
\!\({{x1 -> 0, x2 -> 0}, {x1 -> 2,
      x2 -> \(-\(\(2\ \[ImaginaryI]\ omg10\)\/\@g\)\)}, {x1 -> 2,
      x2 -> \(2\ \[ImaginaryI]\ omg10\)\/\@g}}\)
没有了变量t,好象不太对,
尝试在方程中去掉导数项：
Solve[{omg10^2 x1 == -g/2 x2^2, omg20^2 x2 == omg20^2/2 x1 x2}, {x1, x2}]
结果相同：
\!\({{x1 -> 0, x2 -> 0}, {x1 -> 2,
      x2 -> \(-\(\(2\ \[ImaginaryI]\ omg10\)\/\@g\)\)}, {x1 -> 2,
      x2 -> \(2\ \[ImaginaryI]\ omg10\)\/\@g}}\)

bruce0601

我对Mathematica 不了解 真的谢谢你
可以用matlab编程吗？[br][br]-=-=-=-=- 以下内容由 bruce0601 在 时添加 -=-=-=-=-
请教 大虾 能做弹簧摆振动的相图吗？
谢谢！谢谢！

bruce0601

神啊！
救救我吧！
我现在在做  弹簧摆振动的相图可是就是做不出来
急用啊！！！
要不毕不了业了！！！
望高手指点啊！！！！