CAS (Computer Algebra System) 소개 비선형진동연구실 기계진동 2007-1 CAS 소개 CAS (Computer Algebra System) 소개 학부연구생 4학년 류대균

비선형진동연구실 기계진동 2007-1 CAS 종류 Matlab (Matrix Laboratory) 수치 해석, 행렬 연산, 신호 처리, 간편한 그래픽 기능 등을 통합하여 고성능의 수치 계산 및 결과의 가시화 기능을 제공하는 프로그래밍 언어. 행렬을 기초로 하여 수학과 공학 계산을 수행하는 언어로, 원래 Cleve Moler가 포트란으로 작성하였으나, 현재는 미국의 MathWorks사가 C++로 작성하였다. MATLAB 코딩 체계는 수학적인 기호와 간단한 C 문법으로 되어 있는데, 수학과 관련된 계산, 알고리즘 개발, 상황 모델링과 분석, 여러 가지 과학과 공학적인 그래픽 표현 등에서 이용된다. => 도서관 참고서적 이용 2. Mathematica => 인터넷 강의지원 3. Maple => 도서관 참고서적 이용

비선형진동연구실 기계진동 2007-1 Matlab의 기본연산(행렬)

비선형진동연구실 기계진동 2007-1 Matlab의 기본연산(행렬)

비선형진동연구실 기계진동 2007-1 Matlab의 기본연산(복소수, 함수의 계산과 그래프)

비선형진동연구실 기계진동 2007-1 Matlab의 기본연산(방정식의 해) ans = 1.5708 ans = -0.5000 + 0.5000i -0.5000 - 0.5000i ans = 1.0000 -2.0000 3.0000

비선형진동연구실 기계진동 2007-1 Matlab의 기본연산(고유치 문제)

비선형진동연구실 기계진동 2007-1 Matlab의 기본연산(미분방정식의 해)

비선형진동연구실 기계진동 2007-1 (예제 2.19번) 쿨롱 감쇠계의 자유진동 응답 2차 선형인 상미분 방정식(ODE)  2개의 1차 상미분 방정식(ODE)  연립하여 계산  행렬을 이용하면 편리함  Matlab 사용 (선형/비선형 계산가능)

비선형진동연구실 기계진동 2007-1 (예제 2.19번) 쿨롱 감쇠계의 자유진동 응답 (비선형인 2차 ODE) % Ex2_19.m % This program will use dfunc1.m tspan = [0: 0.05: 8]; => 시간과 간격 x0 = [5.0; 0.0]; => t=0, 초기값 [t, x] = ode23('dfunc1', tspan, x0); % ode23, ode45 => 부록F. 13 참고 plot(t, x(:, 1)); % x(:, 1) 1열의 모든 값 % xlabel('t'); ylabel('x(1)'); title('Example 2.19'); % dfunc1.m function f = dfunc1(t, x) f = zeros(2,1); => 모두 0 인 2 X 1 인 행렬 f(1;1) = x(2); => 1행1열의 함수지정 f(2;1) = -0.5 * 9.81 * sign(x(2)) - 200 * x(1) / 10; => 2행1열의 함수지정

비선형진동연구실 기계진동 2007-1 (예제 2.18번) 스프링-질량계의 자유진동 응답 % Ex2_18.m for i = 1: 101 t(i) = 6 * (i-1)/100; x(i) = 3.1048 * sin(5 * t(i) + 1.3102); x1(i) = 15.524 * cos(5 * t(i) + 1.3102); x2(i) = -77.62 * sin(5 * t(i) + 1.3102); end subplot(311); plot(t,x); ylabel('x(t)'); title('Example 2.18'); subplot(312); plot(t,x1); ylabel('x^.(t)'); subplot(313); plot(t,x2); xlabel('t'); ylabel('x^.^.(t)');

비선형진동연구실 기계진동 2007-1 (예제 2.17번) 정적변형에 따른 고유 진동수와 고유 진동주기 변화 % Ex2_17.m g = 9.81; for i = 1: 101 t(i) = 0.01 + (0.5-0.01) * (i-1)/100; w(i) = (g/t(i))^0.5; tao(i) = 2 * pi * (t(i)/g)^0.5; end plot(t,w); gtext('w_n'); hold on; plot(t,tao); gtext('T_n'); xlabel('Delta_s_t'); title('Example 2.17');

공업수학(1,2) - 2학년 동역학 - 2학년 2학기 기계진동 - 3학년 1학기 자동제어 - 3학년 1학기 비선형진동연구실 기계진동 2007-1 Matlab을 주로 사용하는 기계관련교과목 공업수학(1,2) - 2학년 동역학 - 2학년 2학기 기계진동 - 3학년 1학기 자동제어 - 3학년 1학기 동 및 제어실험 - 3학년 2학기 로봇공학 - 3학년 2학기

비선형진동연구실 기계진동 2007-1 끝맺음 CAS는 이론적인 증명과 확인 및 활용이 가능 - 경청해 주셔서 감사합니다.-