如果我已经用描述质量弹簧阻尼器系统的一阶微分得到了二阶微分方程,那么当我不知道一阶微分时如何使用欧拉方法绘制该方程?我想在MatLab中做到这一点。这是一个作业问题,所以我没有发布任何代码。.我只是想简要概述一下您的意图。
二阶导数为d2(t + 1)=(-1 / m)*(c * d1 + k * y)其中c,m,k为常数,y最初为1,d1为一阶微分从0开始,t是时间。
有任何想法吗?
谢谢 :)。
回答:
可以将二阶eqn转换为一阶微分方程组。
function dy = ex(y) dy = zeros(2,1); dy(1) = y(2); dy(2) = -c/m*y(2) - k/m*y(1); 由此,您可以使用Matlab的内置求解器。 ode23s可以正常工作:
[t,y] = ode23s(@ex, y0, tspan)
更多&回答...