计算单摆运动一周的过程 MATLAB(dθ)^2 /(dt)^2 =-0.98 sinθθ(0)=π/4dθ(0)/dt=0求大神给一个代码,急用啊..555

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 20:34:21

计算单摆运动一周的过程 MATLAB(dθ)^2 /(dt)^2 =-0.98 sinθθ(0)=π/4dθ(0)/dt=0求大神给一个代码,急用啊..555
计算单摆运动一周的过程 MATLAB
(dθ)^2 /(dt)^2 =-0.98 sinθ
θ(0)=π/4
dθ(0)/dt=0
求大神给一个代码,急用啊..555

计算单摆运动一周的过程 MATLAB(dθ)^2 /(dt)^2 =-0.98 sinθθ(0)=π/4dθ(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)];