뛰는 놈 위에 나는 공대생

역학 관련 분야에서는 미분방정식이 자주 쓰이는데 이런 미분방정식을 푸는 방법은 미분방정식 수업에서 배웠을 것이다. 물론 쉬운 미분방정식은 푸는 방법이 바로 생각날 수도 있지만 그렇지 않을 때 간단하게 답을 체크하는 방법이 있다. 1. 울프람알파 사용 공대생들은 대부분 아는 울프람알파(wolfram alpha) 사이트에서 식을 검색해보는 것이다. https://www.wolframalpha.com/ Wolfram|Alpha: Making the world’s knowledge computable Wolfram|Alpha brings expert-level knowledge and capabilities to the broadest possible range of people—spanning all prof..

sym을 이용해서 symbolic 식을 쓰면 편리하게 기호가 들어간 식 계산을 수행할 수 있다. 그런데 그렇게 계산하던 중에 분수가 들어가는 경우 다음과 같이 분수가 제대로 계산되지 않고 길게 늘어진 것을 확인했다. 굳이 식을 확인하지 않으려면 무시하고 써도 상관은 없지만 보기 편하게 하기 위해서 다음과 같은 코드를 사용했다. pos라는 식이 저장된 변수를 각 소수에 대해 다섯자리만 나오게 하려면 vpa(pos,5) 다음과 같이 입력하면 된다. pos 대신 자신이 만든 식을 넣으면 된다. 그렇게 하면 분수가 계산되어서 위처럼 간단하게 식을 볼 수 있다.

1. Matrix에 부가적인 설명 넣기 matrix 아래에 글자를 넣고 싶을 수 있습니다. $\underset{\Sigma}{\underbrace{\begin{bmatrix}\sigma & \omega \\ -\omega & \sigma\end{bmatrix}}}$ 이렇게 하려면 행렬을 감싸는 brace가 필요하고 밑에 글자도 써야하므로 복잡한 구조를 가집니다. \underset{\Sigma}{\underbrace{\begin{bmatrix}\sigma & \omega \\ -\omega & \sigma\end{bmatrix}}} 2. matrix에 가로 세로로 선 긋기 & array 칸의 간격 늘리기 큰 행렬의 경우 block diagonal을 표현하고 싶을 때 유용한 기능입니다. 이렇게 실선으로 표현..

이 내용은 purdue 대학의 글쓰기 센터에서 알게된 comma를 쓰는 규칙에 대해 번역한 글이다. 영어권에서 살지 않는 이상 문장을 쓸 때 미묘하게 나는 comma 사용법을 잘 구별하지 못해서 이렇게 정리된 글을 찾아보았다. https://owl.purdue.edu/owl/general_writing/punctuation/commas/extended_rules_for_commas.html Extended Rules for Commas // Purdue Writing Lab The Purdue University Online Writing Lab serves writers from around the world and the Purdue University Writing Lab helps writers ..

선형대수학에서 여러 matrix category가 있는데 맨날 정의를 잊어버리고 각 개념의 포함관계를 잘 알 수 없었다. 그래서 책(아래 참고문헌 적어놓음)을 참고하여 이번에 정리를 하고자 한다. 첫 번째 분류 1. Square / Non-square matrix matrix의 행과 열의 길이가 같으면 square matrix라고 한다. 그렇지 않은 non-square matrix는 선형대수학을 공부하면서 거의 다루지 않지만 dynamics를 공부하다보면 흔하게 나온다. 일반적인 경우에 변수와 방정식의 갯수가 일치하기가 쉽지 않기 때문이다. 그래서 이런 non-square matrix도 분석할 수 있는 Singular Value Decomposition(SVD)를 수행할 수 있다. SVD는 square ..

Introduction 설명을 시작하기 앞서 우리가 구하고자 하는 Ordinary Differential Equation(ODE)는 다음과 같습니다. $y' = f(y,t)$ 저번 글에서 작성했던 RK4 식은 아래와 같습니다. $y^{n+1}=y^{n}+\frac{1}{6} k_{1}+\frac{1}{3}\left(k_{2}+k_{3}\right)+\frac{1}{6} k_{4}$ $k_{1}=\Delta t f\left(y^{n}, t^{n}\right)$ $k_{2}=\Delta t f\left(y^{n}+\frac{1}{2} k_{1}, t^{n}+\frac{\Delta t}{2}\right)$ $k_{3}=\Delta t f\left(y^{n}+\frac{1}{2} k_{2}, t^{n}+\fra..

앞으로는 코딩할 때 에러가 생기면 기록해놓기로 했다. 수많은 에러를 접하면서 한국어 뿐 아니라 영어로 된 해결법을 찾는게 어려워서 이렇게 글을 써놓으면 누군가에게 도움이 될 수 있지 않을까 한다. modulenotfounderror: no module named 'pytz.tzinfo' pandas 라이브러리를 import할 때 발생했던 문제이다. seaborn을 깔기 위해서 anaconda prompt에서 pip install seaborn을 하면서 pandas 버전도 바뀌었다. 그래서 pytz가 설치되었던 것 같은데도 pandas를 import할 때 문제가 발생했다. 해결방법은 anaconda prompt에서 conda install pytz를 하여 다시 설치했더니 해결되었다.

이전에 DATCOM에 관한 글을 써놓고 예시를 안 적은 것 같아서 내가 만든 코드의 예시를 보여주고자 한다. (도움이 되실지는 모르겠지만..) 그리고 그 .INP 파일로 DATCOM을 돌렸을 때 나오는 결과(.out 파일)도 보여주도록 한다. 다만 .out 파일은 엄청 길기 때문에 핵심적인 부분만 올린다. 1. INP 파일 예시 아래 코드는 .INP 파일의 예시이다. DAMP DERIV RAD DIM M CASEID Project-BODY $FLTCON NMACH = 1.0, MACH(1)=0.03, NALPHA = 20.0, ALSCHD(1) = -2.0, -1.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 12.0, 14.0, 16.0, 18..

ODE에 대한 Introduction에서 적분을 할 때 정해놓은 $\Delta t$가 아니라 그 안에서 여러 스텝을 거쳐서 적분을 하는 방법은 substep method이고, $t^{n+1}$라는 미래 시간을 계산하기 위해 현재와 그 이전 스텝을 쓰는 것을 multi step method라고 했습니다. 이번 글과 다음 글에서 소개할 방법은 multi step method의 범주에 드는 방법들입니다. 위의 그림처럼 $\Delta t$ 간격의 시간에서 $t^{n-1}$ 등 과거 스텝도 사용합니다. 이 방법에는 몇 가지 단점이 있습니다. 첫 번째, timestep size가 고정됩니다. 두 번째, self starting이 불가능합니다. 즉, 초기 조건을 알아도 그 초기 조건을 현재로 생각했을 때 이전 스텝도..

* 이 글의 그림에서 사용된 아이콘은 모두 flaticon의 Freepik에서 나온 것입니다. 아래 링크로 들어가면 찾아볼 수 있습니다. User icons created by Freepik - Flaticon Coordinated Turn 항공기의 coordinate turn은 slip, skid 없이 롤각으로 이루어지는 turn이다. 원운동하듯이 일정 요각 변화율($\dot{\psi}$)로 움직이고 있는 비행기가 있을 때 이 요각 변화율은 bank angle 덕분에 만들어진 것이다. 일정 요각 변화율로 turn을 하고 싶다면 다음 조건이 만족해야한다. $mU_{0}\dot{\psi} \cos\phi = mg\sin \phi$ $\Rightarrow \tan\phi =\frac{U_{0}\dot{\ps..