Solve equation of motion using ode45!
조회 수: 2 (최근 30일)
이전 댓글 표시
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
채택된 답변
추가 답변 (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!