四旋翼无人机建模与实现(三)

四旋翼无人机建模与实现(三)

本文是对GitHub上的一个开源的四旋翼无人机建模与编程实现的学习笔记,因为原作者的文档使用英文写的,花了很长的时间才看懂,因此,我觉得用自己的语言将这个开源项目分享出来是有意义的,爱因斯坦曾说过:“如果你不能简单地解释某一概念,说明你没有很好地掌握它”。因此便有了写这篇文档的想法。
原文链接: https://github.com/SKYnSPACE/AE450
本文的贡献在于将原作者文档中的公式做了更为详细的推导,并以更为通俗的语言来描述整个建模过程。


10、计算状态矩阵A和输入矩阵B

本节将详细介绍对控制模型进行线性化过程,也即计算状态矩阵A和输入矩阵B,首先,矩阵A和矩阵B是一个稀疏矩阵,将被控对象的整体模型转化为矩阵的形式如下,当然,这个只是矩阵的形式而已,对求解过程没啥用。
将上述模型写成矩阵的形式:(这里写成这样没意义,只是展示一下这是一个稀疏矩阵而已)
[ X ˙ Y ˙ Z ˙ X ¨ Y ¨ Z ¨ ϕ ˙ ψ ˙ θ ˙ p ˙ q ˙ r ˙ ] = [ 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 s i n ϕ t a n θ c o s ϕ t a n θ 0 0 0 0 0 0 0 0 0 0 c o s ϕ − s i n ϕ 0 0 0 0 0 0 0 0 0 0 s i n ϕ s e c ϕ c o s ϕ s e c θ 0 0 0 0 0 0 0 0 0 0 0 ( I x x − I z z ⋅ q ) / I x x 0 0 0 0 0 0 0 0 0 0 0 ( I z z − I x x ⋅ p ) / I y y 0 0 0 0 0 0 0 0 0 ( I x x − I y y ⋅ q ) / I z z 0 0 ] [ X Y Z X ˙ Y ˙ Z ˙ ϕ ψ θ p q r ] + [ 0 0 0 0 0 0 0 0 0 0 0 0 − 1 m ( s i n ψ s i n ϕ + c o s ψ s i n θ c o s ϕ ) 0 0 0 − 1 m ( − c o s ψ s i n ϕ + s i n ψ s i n θ c o s ϕ ) 0 0 0 − 1 m c o s θ c o s ϕ + g / T ∑ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 / I x x 0 0 0 0 1 / I y y 0 0 0 0 1 / I z z ] [ T ∑ M 1 M 2 M 3 ] \boldsymbol{\begin{bmatrix}\dot{X}\\ \dot{Y}\\\dot{Z}\\\ddot{X}\\\ddot{Y}\\\ddot{Z}\\ \dot{\phi}\\\dot{\psi}\\\dot{\theta}\\\dot{p}\\\dot{q}\\\dot{r}\end{bmatrix} =\begin{bmatrix} 0&0&0&1&0&0&0&0&0&0&0&0\\ 0&0&0&0&1&0&0&0&0&0&0&0\\ 0&0&0&0&0&1&0&0&0&0&0&0\\ 0&0&0&0&0&0&0&0&0&0&0&0\\ 0&0&0&0&0&0&0&0&0&0&0&0\\ 0&0&0&0&0&0&0&0&0&0&0&0\\ 0&0&0&0&0&0&0&0&0&1&sin\phi tan\theta&cos\phi tan\theta\\ 0&0&0&0&0&0&0&0&0&0&cos\phi&-sin\phi\\ 0&0&0&0&0&0&0&0&0&0&sin\phi sec\phi&cos\phi sec\theta\\ 0&0&0&0&0&0&0&0&0&0&0&(I_{xx}-I_{zz}\cdot q)/I_{xx}\\ 0&0&0&0&0&0&0&0&0&0&0&(I_{zz}-I_{xx}\cdot p)/I_{yy}\\ 0&0&0&0&0&0&0&0&0&(I_{xx}-I_{yy}\cdot q)/I_{zz}&0&0\\ \end{bmatrix} \begin{bmatrix}{X}\\ {Y}\\{Z}\\\dot{X}\\\dot{Y}\\\dot{Z}\\ {\phi}\\{\psi}\\{\theta}\\{p}\\{q}\\{r}\end{bmatrix}+ \begin{bmatrix} 0&0&0&0\\ 0&0&0&0&\\ 0&0&0&0\\ -\frac{1}{m}(sin\psi sin\phi+cos\psi sin\theta cos\phi)&0&0&0\\ -\frac{1}{m}(-cos\psi sin\phi+sin\psi sin\theta cos\phi)&0&0&0&\\ -\frac{1}{m}cos\theta cos\phi+g/T_{\sum}&0&0&0\\ 0&0&0&0\\ 0&0&0&0&\\ 0&0&0&0\\ 0&1/I_{xx}&0&0\\ 0&0&1/I_{yy}&0&\\ 0&0&0&1/I_{zz}\\ \end{bmatrix} \begin{bmatrix}T_{\sum}\\M_1\\M_2\\M_3\end{bmatrix} } X˙Y˙Z˙X¨Y¨Z¨ϕ˙ψ˙θ˙p˙q˙r˙=00000000000000000000000000000000000010000000000001000000000000100000000000000000000000000000000000000000000000000010000(IxxIyyq)/Izz000000sinϕtanθcosϕsinϕsecϕ000000000cosϕtanθsinϕcosϕsecθ(IxxIzzq)/Ixx(IzzIxxp)/Iyy0XYZX˙Y˙Z˙ϕψθpqr

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值