Question

matlab二阶微分方程求解及画图

取m=10600kg k=31986497 N/m c=8.9×10^4kg/mx（0）=0； x'（0）=3.66
改区间为0~1
求全部代码，以及matlab版本

Answer 1

参考代码：

m=10600; k=31986497; c=8.9E4;
dy = @(t,y)[y(2); (k*y(1)-c*y(2)/abs(y(2))*y(2)^2)/m];
ode45(dy,[0 5],[0 3.66])

其中，蓝线为x，绿线为x'。

追问

"identifier" expected, "(" found. 我直接复制代码去运行报错了，那个@是不是还要先定义函数？

追答

你的MATLAB是什么版本？是不是6.x？
@用于定义匿名函数，需要7.0以上版本支持。