Reference : Principles of GNSS Inertial and Multisensor Integrated Navigation Systems
Attitude : orientation of axes of one coordinate w.r.t another coordinate frame.
Coordinate Transformation Matrix, Rotation Matrix : Cαβ
xβ=Cαβxα : 임의의 벡터 x의 resolving axes를 α에서 β로 변환
위와 같이 Cαβ는 vector의 resolving frame을 α에서 β로 변환시키는 coordinate transformation Matrix이면서, β에서 α로의 Rotation Matrix. 2가지 의미를 갖는다.
Kinematics
xβαγ
- α : object frame. x의 motion이 기술되는 프레임.
- β : reference frame : motion의 기준이 되는 프레임.
- γ : resolving frame : motion이 표현되는 프레임.
예를 들어 x가 position(위치)에 대한 물리량이라고 하면, xβαγ는 β frame 기준 α frame의 위치를 γ frame으로 표현한 값이다.
1. Angular rate
wβαγ : β frame 기준 α frame의 회전속도를 γ frame으로 표현한 벡터.
wβαγ=wβδγ+wδαγ, wβαδ=Cγδwβαγ
α frame은 회전하고, β frame은 회전하지 않는다고 가정하면
Cβα(t+δt)=Cα(t)α(t+δt)Cβα(t)=exp(−[wβαα]×δt)Cβα(t)≈(I−[wβαα]×δt)Cβα(t)
C˙βα=limδt−>0(δtCβα(t+δt)−Cβα(t))=−[wβαα]×Cβα를 얻는다.
반대로 α frame은 회전하지 않고, β frame이 회전한다고 가정하면
C˙βα=−Cβα[wβαβ]×를 얻는다.
Cβα[wαββ]×Cαβ=[wαβα]×이므로
C˙βα=−Cβα[wβαβ]×=−[wβαα]×Cβα
2. Cartesian Position
rβαγ=(xβαγ,yβαγ,zβαγ) : β frame 기준 α frame의 origin의 위치를 γ frame으로 표현한 값.
rβαγ=−rαβγ, rβαγ=rδαγ−rδβγ=rβδγ+rδαγ
rβαδ=Cγδrβαγ
3. velocity
vβαγ=Cβγr˙βαβ : b frame기준 a 프레임의 위치의 변화율을 γ frame으로 표현한 값.
position의 경우와 달리, ref frame과 resolving frame이 서로 회전하는 경우(Cβγ=0˙) resolving frame에서 표현되는 position의 변화율과 velocity는 서로 다르다.
r˙βαγ=Cβγrβαβ˙=C˙βγrβαβ+Cβγr˙βαβ=C˙βγrβαβ+vβαγ
또한 position과 달리 velocity는 마이너스를 곱하고 reference frame과 object frame을 바꾼 값이 원래의 값과 일치하지 않는다.
vαβγ=−vβαγ−CαγC˙βαrβαβ
또한 position과 달리 velocity의 vector sum을 사용할 수 없다.
vβαγ=vβδγ+vδαγ
velocity의 resolving frame을 변환하는 것은 가능하다.
rate of change of position(position의 변화율)과는 다름을 주의.
vβαδ=Cγδvβαγ
4. acceleration
aβαγ=Cβγr¨βαβ
reference frame과 resolving frame이 서로 회전하는 경우 v˙와 r¨은 아래와 같다.
v˙βαγ=Cβγr˙βαβ˙=C˙βγr˙βαβ+Cβγr¨βαβ=C˙βγr˙βαβ+aβαγ
r¨βαγ=Cβγrβαβ¨=C¨βγrβαβ+2C˙βγr˙βαβ+aβαγ
angular rate에서 구한 식을 이용하면
C¨βγrβαβ=([wβγγ]×2−[wβγγ]˙×)rβαγ
Cβγ˙r˙βαβ=−[wβγγ]×2rβαγ−[wβγγ]×r˙βαγ를 얻을 수 있고, 두 식을 사용해 정리하면
r¨βαγ=−[wβγγ]×2rβαγ−2[wβγγ]×r˙βαγ−[wβγγ]˙×rβαγ+aβαγ
1, 2, 3번째 term은 각각 centrifugal, coriolis, euler pseudo-force에 관련된 값이다.
velocity와 마찬가지로 vector sum이 가능하지 않고, accel의 좌표계 변환은 가능하다.
aβαγ=aβδγ+aδαγ
aβαδ=Cγδaβαγ