这是为什么
mathematica吧
全部回复
仅看楼主
level 3
g = 9.8; ll = 1.5; ww = Sqrt[g/ll];
qt = {};
Do[w0 = v0/ll;
s = NDSolve[{q''[t] == -ww^2 Sin[q[t]], q[0] == 0, q'[0] == w0},
q, {t, 0, 10}]; a = (q /. s[[1]]);
qm1 = FindMaximum[a[t], {t, 0.5}];
qm2 = FindMaximum[a[t], {t, 3.3}];
tt = (t /. qm2[[2]]) - (t /. qm1[[2]]);
AppendTo[qt, {qm1[[1]] 180/Pi, tt}], {v0, 0.1, 5.4, 0.1}]
2022年06月08日 09点06分 1
吧务
level 10
你没贴出图的代码啊
2022年06月09日 14点06分 2
level 3
不好意思,我忘了
g1 = ListPlot[qt]
2022年06月10日 07点06分 3
吧务
level 10
我这里正常出图了
2022年06月11日 02点06分 4
1