mathematica 解微分方程组的问题,谢谢了R = 0.008;r = 0.0125;Y = 47.0*2.0*Pi;Z = 0.09;m = 0.08;k = 9.0;g = 9.8;f = 0.000586;e = 28;NDsolve[{R^2*y[t]*y'[t] + z[t]*z'[t] == 0, R^2*Y^2 + Z^2 == R^2*(Y[t])^2 + (z[t])^2, 0.5*m*(z'[t])^2 + 0.
来源：学生作业帮助网 编辑：六六作业网 时间：2024/11/26 00:37:44
mathematica 解微分方程组的问题,谢谢了R = 0.008;r = 0.0125;Y = 47.0*2.0*Pi;Z = 0.09;m = 0.08;k = 9.0;g = 9.8;f = 0.000586;e = 28;NDsolve[{R^2*y[t]*y'[t] + z[t]*z'[t] == 0, R^2*Y^2 + Z^2 == R^2*(Y[t])^2 + (z[t])^2, 0.5*m*(z'[t])^2 + 0.
mathematica 解微分方程组的问题,谢谢了
R = 0.008;
r = 0.0125;
Y = 47.0*2.0*Pi;
Z = 0.09;
m = 0.08;
k = 9.0;
g = 9.8;
f = 0.000586;
e = 28;
NDsolve[{R^2*y[t]*y'[t] + z[t]*z'[t] == 0,
R^2*Y^2 + Z^2 == R^2*(Y[t])^2 + (z[t])^2,
0.5*m*(z'[t])^2 + 0.2*m*r^2*(y'[t])^2 + 0.5*k*(z[t] - Z)^2 +
0.5*f*(y[t] - Y)^2 - m*g*(z[t] - Z) == e, y[0] == Y,
z[0] == Z}, {z[t], y[t]}, {t, 0, 1}, MaxSteps -> 20000];
Plot[Evaluate[y[t] /. %], {t, 0, 1}]
之前问过的,我不太懂是哪里的问题,谢谢了
mathematica 解微分方程组的问题,谢谢了R = 0.008;r = 0.0125;Y = 47.0*2.0*Pi;Z = 0.09;m = 0.08;k = 9.0;g = 9.8;f = 0.000586;e = 28;NDsolve[{R^2*y[t]*y'[t] + z[t]*z'[t] == 0, R^2*Y^2 + Z^2 == R^2*(Y[t])^2 + (z[t])^2, 0.5*m*(z'[t])^2 + 0.
首先,是NDSolve而不是NDsolve,Mathematica是区分大小写的.
其次……其实你看看警告信息就该明白的——方程过定了,也就是,你给的方程数目过多了（两个未知数你给了三个方程）,试着去掉了那个代数方程（也就是第二个方程,不含导数的那个）之后,就可解了.不过积分步长又是早早耗光了,但是好歹算到了0.3多（图上只到了0.15,那是因为后面是复根）,比之前好很多嘛.
————
试着强制提升了一下计算精度,发现只能算到0.15左右.看来默认情况下能算到0.3是因为误差……
————
仔细一看,你的那个代数方程还写错了……是y[t]而不是Y[t].然后,去掉你的第一个方程,保留第二三个方程,这个问题就可解了!你试试!：