<div>안녕하세요 매틀랩 공부중인 공대 학생입니다,</div> <div><br></div> <div>진동학을 매틀랩으로 해석하는 코드를 책에 있는 그대로 복붙 했는데 실행이 안되네요 ㅠㅠ</div> <div><br></div> <div>파일 두개 만들어서 했는데 실행이 안됩니다..</div> <div><br></div> <div>혹시 아시는분 있으시면 알려주심 감사하겠습니다!!</div> <div> <p>////////////////////////////////////////</p> <p>% M5_8.m</p> <p>close all; clear </p> <p>x1=1; dotx1=1; x2=0; dotx2=0; </p> <p>x0=[x1; dotx1; x2; dotx2]; t=0:0.01:10; </p> <p>[t,x]=ode23('song5_8',t,x0); plot(t, x(:,1), 'linewidth', 2); hold on </p> <p>plot(t, x(:,3), 'r:'); hold on plot(t, 0*t, 'black')</p> <p>axis([0 10 -1.5 1.5]) xlabel('Time (s)'); xlabel('Displacement (m)'); grid legend('x_1(t)', 'x_2(t)');</p> <p>///////////////////////////////////////</p> <p>% song5_8.m</p> <p>funtion y=song5_8(t,x) y=zeros(4,1); m1=1; m2=1; c1=0; c2=0; c3=0; k1=4; k2=4; k3=4; f1=0; f2=0; w1=3; w2=5;</p> <p>y(1)=x(2); </p> <p>y(2)=(-(k1+k2)*x(1)-(c1+c2)*x(2)+k2*x(3)+c2*x(4)+f1*cos(w1*t))/m1; </p> <p>y(3)=y(4); </p> <p>y(4)=(k2*x(1)+c2*x(2)-(k2+k3)*x(3)-(c2+c3)*x(4)+f2*cos(w2*t))/m2;</p> <p>//////////////////////////////////////////////////////</p></div> <div><br></div>