ardupilot开发 --- 无人机数学模型与控制律分解 篇
坐标系
- 机体系
FRD,右手系,原点取在几何中心或质心。
由于定义欧拉角;
用于建立动力学模型、运动学模型、动力分配模型;
是所有方程物理量需要归一化的基准坐标系。 - 导航坐标系
地固坐标系,NED。
是欧拉角的参考原点、基准系。
用于导航与制导。 - 全球系
定位传感器如GPS所在的坐标系,需要把GPS位置信息转换到导航坐标系NED下。 - 传感器坐标系
方向与机体系同,原点取传感器几何中心。
FRD。
用于卡尔曼滤波、状态估计。
数学模型
运动学模型
基于欧拉角
[px˙py˙pz˙]=[vxvyvz]
\begin{bmatrix}
\dot{p_x} \\
\dot{p_y} \\
\dot{p_z}
\end{bmatrix} =
\begin{bmatrix}
v_x \\
v_y \\
v_z
\end{bmatrix}
px˙py˙pz˙=vxvyvz
[ϕ˙θ˙ψ˙]=[1sinϕtanθcosϕtanθ0cosϕ−sinϕ0sinϕsecθcosϕsecθ][wxwywz]
\begin{bmatrix}
\dot{\phi} \\
\dot{\theta}\\
\dot{\psi}
\end{bmatrix}=
\begin{bmatrix}
1 & \sin\phi \tan\theta & \cos\phi \tan\theta \\
0 & \cos\phi & -\sin\phi \\
0 & \sin\phi \sec\theta & \cos\phi \sec\theta
\end{bmatrix}
\begin{bmatrix}
w_x \\
w_y\\
w_z
\end{bmatrix}
ϕ˙θ˙ψ˙=100sinϕtanθcosϕsinϕsecθcosϕtanθ−sinϕcosϕsecθwxwywz
其中:
- [px,py,pz],[vx,vy,vz][p_x,p_y,p_z],[v_x,v_y,v_z][px,py,pz],[vx,vy,vz]分别是导航坐标系NED下的位移矢量和速度矢量。
- ϕ,θ,ψ\phi ,\theta, \psiϕ,θ,ψ是欧拉角。注意,欧拉角的基准是导航系。
- wx,wy,wzw_x,w_y,w_zwx,wy,wz是机体系下的角速度,是陀螺仪的测量值。
基于四元素
[px˙py˙pz˙]=[vxvyvz]
\begin{bmatrix}
\dot{p_x} \\
\dot{p_y} \\
\dot{p_z}
\end{bmatrix} =
\begin{bmatrix}
v_x \\
v_y \\
v_z
\end{bmatrix}
px˙py˙pz˙=vxvyvz
[q˙0q˙1q˙2q˙3]=12[q0q1q2q3]⊗[0wxwywz]=12[qw−qx−qy−qzqxqw−qzqyqyqzqw−qxqz−qyqxqw][0ωxωyωz]
\begin{bmatrix}
\dot q_0 \\
\dot q_1 \\
\dot q_2\\
\dot q_3
\end{bmatrix} = \frac{1}{2}
\begin{bmatrix}
q_0 \\
q_1 \\
q_2\\
q_3
\end{bmatrix}
\otimes
\begin{bmatrix}
0 \\
w_x \\
w_y\\
w_z
\end{bmatrix}
= \frac{1}{2} \begin{bmatrix}
q_w & -q_x & -q_y & -q_z \\
q_x & q_w & -q_z & q_y \\
q_y & q_z & q_w & -q_x \\
q_z & -q_y & q_x & q_w
\end{bmatrix}
\begin{bmatrix}
0 \\
\omega_x \\
\omega_y \\
\omega_z
\end{bmatrix}
q˙0q˙1q˙2q˙3=21q0q1q2q3⊗0wxwywz=21qwqxqyqz−qxqwqz−qy−qy−qzqwqx−qzqy−qxqw0ωxωyωz
或者把上式拆成bw=[wx,wy,wz]^bw=[w_x,w_y,w_z]bw=[wx,wy,wz]作为输入的形式:
其中:
- q0q_0q0:标量部分。
- qv=[q1,q2,q3]Tq_v=[q_1,q_2,q_3]^Tqv=[q1,q2,q3]T:矢量部分。
- q02+q12+q22+q32=1q_0^2 + q_1^2 + q_2^2 + q_3^2 = 1q02+q12+q22+q32=1。