#진동학 #매트랩 

 

안녕하세요. 이전에 작성했던 글을 잘못 올린 것 같아서 다시 올립니다. 갓준표 사이트와 유튜브를 한번 훑어봐도 내용이 잘 안나와서 질문드립니다.

과도감쇠운동에서 매트랩으로 표현하실때 syms를 안쓰고 그냥 a1 , a2값을 집어넣어 계산을 하셨어 가지고 제가 한번 syms를 사용하여 타이핑을 쳐서 계산을 도출해봤습니다.

k=1;

m=1;

c=2.2;

x0=1;

v0=0;

t = linspace(0,10);

 

wn=sqrt(k/m);

xi=c/(2*m*wn);

wd=sqrt(1-xi^2)*wn;

%a1 = (-v0+wn*x0*(-xi+sqrt(xi^2-1)))/(2*wn*sqrt(xi^2-1));

%a2 = (v0+wn*x0*(xi+sqrt(xi^2-1)))/(2*wn*sqrt(xi^2-1));

 syms a1 a2

 f = a1+a2-x0;

 g = x0-a1-a2;

 [a1 a2 ] = solve(f,g);

x = exp(-xi*wn*t).*(a1*exp(-wn*sqrt(xi^2-1)*t)+a2*exp(wn*sqrt(xi^2-1)*t));

plot(x);

이렇게 과도감쇠운동에서 syms를 사용하여 사용을 했더니 영상에서 결과그래프가 다르게 나와서 질문드립니다. 혹시 어느쪽에서 다르게 틀리게 하여 이상하게 나왔는지 모르겠습니다. 전체적인 타이핑에서 다른건 linspace 범위만 (0,10)으로만 변경하고 다른건 없습니다. 영상에서 적은 a1, a2값은 주석 처리하여 결과에 반영이 안되게 했습니다. 다르게 많이 시도해봤는데 저 혼자서는 답이 안나와서 질문드립니다.