Equation of motion of Non linear pendulum

Question

KLETECH MOTORSPORTS 2020-11-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/636140-equation-of-motion-of-non-linear-pendulum

评论： Bjorn Gustavsson 2020-11-4

tspan=[0,2*pi];
u0=[pi/4;0];
[t,u]=ode45(@pendulum,tspan,u0,[]);
plot(t,u(:,1),'b-', 'LineWidth', 2)
xlabel('time')
ylabel('angle')
function zdot = pend (t,z)  
wsq = 1.56; % specify a value of w"2
t = time;
z = [z (l) ; z ( 2)] = [theta; thetadot] 
zdot = [z (2) ; -wsq* sin (z(l))]; 
end

This doesn't seem to be working. I'm trying to plot the discplacement versus time and velocity versus time. Any idea as to what the problem mightbe?

I got this from Rudra pratap. I've attached a link to the book, it's on page 159 (173 by docs)

https://lpuguidecom.files.wordpress.com/2017/05/p1-11385-scripts-and-matrix-operations.pdf

2 个评论
显示无隐藏无

KLETECH MOTORSPORTS 2020-11-4

( Second Order non linear ODE )

KSSV 2020-11-4

Refer this: https://in.mathworks.com/matlabcentral/fileexchange/33082-simplependulum

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bjorn Gustavsson 2020-11-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/636140-equation-of-motion-of-non-linear-pendulum#answer_533875

在 MATLAB Online 中打开

Note that the function you use for the input-argument to ode45 is pendulum while the function you implemented for the pendulum has the name pend. They should be the same. When that's fixed there seems to be no major problems with your code:

u0=[pi/4;0];
pend = @(t,z) [z(2);-1.56*sin(z(1))];
figure
tspan=[0,12*pi]; % just a slightly longer time-span
[t,u]=ode45(pend,tspan,u0,[]);
plot(t,u)

HTH

3 个评论
显示 1更早的评论隐藏 1更早的评论

KLETECH MOTORSPORTS 2020-11-4

i got this error

Bjorn Gustavsson 2020-11-4

Ok, that is because your script is named pend.m and then matlab doesn't like the naming-confusion with a variable with that name in the script. You should be fine if you rename the variable to someting like pendODE for example.

HTH

请先登录，再进行评论。