作业帮 > 数学 > 作业

计算单摆运动一周的过程 MATLAB

来源:学生作业帮 编辑:拍题作业网作业帮 分类:数学作业 时间:2024/05/16 14:07:33
计算单摆运动一周的过程 MATLAB
(dθ)^2 /(dt)^2 =-0.98 sinθ
θ(0)=π/4
dθ(0)/dt=0
求大神给一个代码,急用啊..555
%单摆角度很小时,sin(θ)≈θy=dsolve('D2y = -0.98*y','y(0)=pi/4,Dy(0)=0');syms t;figure;ezplot(t,y,[0 10]);
%θ较大时,只能求数值解tspan=[0 10];y0=[pi/4 0];[t,x]=ode45(@odefun,tspan,y0);figure;plot(t,x(:,1));
function y = odefun(t,x);y = [-0.98*sin(x(2)); x(1)];