weixin_55982578的博客

Matlab 中的钟摆问题如图编辑器function[dx]=pendulum(t,x)% 公式变形%假设 Tc=0%x(1)=theta x(2)=d(theta)=dx(1)dx(1)=x(2);dx(2)=-g/l*sin(x(1))-b/(m*l*l)*x(2);dx=dx';% 让数据变成纵列%利用ode45公式计算 在 [1；0]角速度是1；初速度为0；在0到30秒里面[0 30]%ode45 公式：[t,y] = ode45(odefun,tspan,y0)%[t