请大家帮忙看看Mathematica的一个方程

scleilei

请大家帮忙看看Mathematica的一个方程

NDSolve[{(0.0371 + 0.0183 b[t] + 0.0030 Cos[c[t]] + 0.061 Cos[b[
                t] + c[t]]) a''[t] + (0.0046 + 0.0091
                Cos[b[t]] + 0.0030 Cos[c[t]] +
                      0.0030 Cos[b[t] + c[t]]) b''[t] + (
          0.00076172 + 0.0015 Cos[c[t]] + 0.0030 Cos[b[t] + c[t]])
                c''[t] + 1.4625 Cos[a[t]] + 0.3656 Cos[a[
                t] + b[t]] + 0.1219 Cos[a[t] + b[
                      t] + c[t]] == 1.9550 - 0.0196 a'[
          t], (0.0046 + 0.0091 Cos[b[t]] + 0.0030 Cos[c[t]] + 0.0030 Cos[b[
                t] + c[t]]) a''[t] + (0.0042 + 0.0030
                      Cos[c[t]]) b''[t] + (0.00076172 + 0.0015 Cos[c[t]]) \
c''[t] + 0.3656 Cos[a[t] + b[t]] + 0.1219 Cos[b[t] + b[t] + c[t]] ==
          0.1625 - 0.0161 b[t], (
              0.00076172 + 0.0030 Cos[b[t] + c[t]] + 0.0015
        Cos[c[t]]) a''[t] + (0.00076172 + 0.0015 Cos[c[
      t]]) b''[t] + 0.0010 c''[t] + 0.1219 Cos[a[t] + b[t] + c[
                t]] == 0.2300 - 0.0026 c[t], a[0] == 0,
                 b[0] == 0, c[0] == Pi/2, a'[0] == 0, b'[0] == 0, c'[0] ==
               0}, {a[t], b[t], c[t]}, {t, 0, 1}];
Plot[a[t], {t, 0, 1}];

请问我要作出a-t,b-t,c-t和a'-t,b'-t,c'-t之间的曲线图,上面的程序有什么问题吗?

windfly

有两个问题。
一个是你的t取值范围好像只能限定在{0，0.12}

另外一个问题是,你的plot函数应该改为：

1. Plot[Evaluate[a[t]/.First[%]],{t,0,0.12}]

复制代码

scleilei

但是出来的曲线图结果应该是向上的才是正确的,就象y=x^2这条曲线右边的部分一样.
那我哪里错了呢?
请再帮我看一下

scleilei

还有我的其他几条曲线(a'-t,b'-t,c'-t)怎么出来呢 ?