Solve equation of motion using ode45!

1 次查看（过去 30 天）
Alaa Hameed
Alaa Hameed 2017-12-11
编辑： Star Strider 2017-12-12
Hi, I am trying to solve an equation of motion using ode45. I saw others work, and how they implement the method, however, I did not get a reasonable answer. It seems so simple, but I do not know why I got such weird answers! Could anyone help me figure it out? My code is as follows:
main body:
t=0:.02:5;
X=[.002 ; 0];
[t,y]=ode45(@YP,t, X);
figure(1)
plot(t,y(1), '+')
figure(2)
plot(t,y(2), 'o')
Function:
function ydot= YP(t,x)
m=2459;k=107.178;zeta=.005;
wn=sqrt(k/m);
c=zeta*(2*m*wn);
ydot=[x(2);(-1/m)*(c*x(2)+k*x(1))];
  2 个评论
James Tursa
James Tursa 2017-12-11
What is the ode you are trying to solve?
Alaa Hameed
Alaa Hameed 2017-12-11
the free vibration : m*A+ c*V+ k*x=0 where, A:acceleration V:velocity x:displacement m:mass; c:damping coefficient; k:stiffness

请先登录，再进行评论。

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

采纳的回答

James Tursa
James Tursa 2017-12-12
figure;plot(t,y(:,1));grid on
figure;plot(t,y(:,2));grid on

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by