mathematica求解微分方程的数值解,做出图像,R = 0.008;Y = 47*2*3.14;z = 0.18;m = 0.14;k = 20;g = 10;f = 20;i = 0.143;e = 28;NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) + 0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/02 04:16:21
![mathematica求解微分方程的数值解,做出图像,R = 0.008;Y = 47*2*3.14;z = 0.18;m = 0.14;k = 20;g = 10;f = 20;i = 0.143;e = 28;NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) + 0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0](/uploads/image/z/9722789-53-9.jpg?t=mathematica%E6%B1%82%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%9A%84%E6%95%B0%E5%80%BC%E8%A7%A3%2C%E5%81%9A%E5%87%BA%E5%9B%BE%E5%83%8F%2CR+%3D+0.008%3BY+%3D+47%2A2%2A3.14%3Bz+%3D+0.18%3Bm+%3D+0.14%3Bk+%3D+20%3Bg+%3D+10%3Bf+%3D+20%3Bi+%3D+0.143%3Be+%3D+28%3BNDSolve%5B%7B0.5%2Am%2A%28R%5E4+%2Ay%5Bt%5D%5E2+%2Ay%27%5Bt%5D%5E2%29%2F%28R%2AY%5E2+%2B+z%5E2+-+R%2Ay%5Bt%5D%5E2%29+%2B+0.5%2Ai%2Ay%27%5Bt%5D%5E2+%2B+0.5%2Ak%2A%28R%2AY%5E2+%2B+z%5E2+-+R%2Ay%5Bt%5D%5E2%29+%2B+0)
mathematica求解微分方程的数值解,做出图像,R = 0.008;Y = 47*2*3.14;z = 0.18;m = 0.14;k = 20;g = 10;f = 20;i = 0.143;e = 28;NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) + 0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0
mathematica求解微分方程的数值解,做出图像,
R = 0.008;
Y = 47*2*3.14;
z = 0.18;
m = 0.14;
k = 20;
g = 10;
f = 20;
i = 0.143;
e = 28;
NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) +
0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0.5*k*z^2 -
k*z*Sqrt[R*Y^2 + z^2 - R*y[t]^2] + 0.5*f*(y[t] - Y)^2 -
m*g*(Sqrt[R*Y^2 + z^2 - R*y[t]^2] - z) == e,y[0] == 0,
y'[0] == 1},y[t],{t,10}];
Plot[Evaluate[y[t] /.%],{t,0,10}]
这是代码,我用的mathematica 8 看不出错误来.
mathematica求解微分方程的数值解,做出图像,R = 0.008;Y = 47*2*3.14;z = 0.18;m = 0.14;k = 20;g = 10;f = 20;i = 0.143;e = 28;NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) + 0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0
NDSolve里的语法错了,是{t,0,10}而不是{t,10};另外对于常微分方程,所需的边界条件数目,为它的最高阶导数的阶数,并且边界的导数阶数要小于方程的阶数.也就是说你这里的条件多了一个,把那个y'[0]==1去掉,并且,你这里的边界给的不合理,所以在计算开始不久就耗尽了全部的积分步长.试了一下也没找到合理的边界条件,你自己试试吧.