新人求助，矩阵形式的微分方程组

这是我的代码，U = 34; b = 0.5; k = 0.098; \[Theta] = 11.84 + 9.87 Sin[k t]; a = -0.1;
NN = 6; Di = Table[Which[i == j + 1, 1/2/i, i == j - 1, -1/2/i, True, 0], {i, 6}, {j, 6}]
bn = Table[((-1)^(i - 1)*((NN + i - 1)!))/((NN - i - 1)!*(i!)^2), {i, NN - 1}];
bn = AppendTo[bn, -1]
d = Prepend[Table[0, {i, NN - 1}], 1/2]
c = Table[2/i, {i, NN}]
A = Di + Transpose[{d}].{bn} + Transpose[{c}].{d} + 1/2 Transpose[{c}].{bn} // MatrixForm
NDSolve[{A.\[Lambda]'[t] + U/b \[Lambda][t] == Transpose[{c}].(U D[\[Theta], t] + b (1/2 - a) D[\[Theta], t, 2]), [Lambda][0] == 0}, \[Lambda][ t], {t, 0, 2 Pi/k}]，最后一步求不出来，和matlab不太一样，不知道该怎么写，请大神指导下，我网上搜也搜不出个所以然来。。附上最后的方程。。感谢感谢

h那一项直接忽略就行

没有人吗。。。

搞定了。。。
\[Lambda] = {\[Lambda]1[t], \[Lambda]2[t], \[Lambda]3[t], \[Lambda]4[
t], \[Lambda]5[t], \[Lambda]6[t]};
equ111 = A.D[\[Lambda], t] + U/b Transpose[{\[Lambda]}] ==
Transpose[{c}] (34 D[\[Theta], t] +
0.5 (1/2 + 0.1) D[\[Theta], {t, 2}])
ans = NDSolve[{equ111, {\[Lambda]1[0], \[Lambda]2[0], \[Lambda]3[
0], \[Lambda]4[0], \[Lambda]5[0], \[Lambda]6[0]} == {0, 0, 0, 0,
0, 0}}, {\[Lambda]1[t], \[Lambda]2[t], \[Lambda]3[
t], \[Lambda]4[t], \[Lambda]5[t], \[Lambda]6[t]}, {t, 0.001,
2 Pi/k}]
还是语法不熟悉的锅

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频！

3回复贴，共1页

<<返回mathematica吧

分享到:

日	一	二	三	四	五	六