发明名称 六旋翼无人机分立输入解耦装置及包含该装置的控制系统
摘要 一种六旋翼无人机分立输入解耦装置及包含该装置的控制系统。本发明涉及一种六旋翼无人机分立输入解耦装置,该装置包括系数矩阵Mco存储模块;结构矩阵Mrel存储模块;将系数矩阵与结构矩阵相乘并存储为映射矩阵Map的模块;映射矩阵降维模块;降维齐次线性方程求解模块;基础解向量归一化模块;差分模块;放大增益模块;求和模块;非线性解算模块;本发明使全驱动六旋翼无人机这个多输入多输出系统解耦为六个单输入单输出系统,从而实现了无人机六个自由度独立控制,其最大优点是解耦之后的系统形式简单易于设计控制算法,由虚拟附加控制量可以快速地求取各个电机实际转速,算法对电机转速发生波动具有一定的鲁棒性。
申请公布号 CN102323758A 申请公布日期 2012.01.18
申请号 CN201110216453.4 申请日期 2011.07.29
申请人 中国科学院长春光学精密机械与物理研究所 发明人 白越;侯至丞;宫勋;高庆嘉;孙强;续志军
分类号 G05B17/00(2006.01)I 主分类号 G05B17/00(2006.01)I
代理机构 长春菁华专利商标代理事务所 22210 代理人 王淑秋
主权项 1.一种六旋翼无人机分立输入解耦装置,其特征在于包括:系数矩阵M<sub>co</sub>存储模块(101);所述的系数矩阵M<sub>co</sub>=diag(k<sub>11</sub>,k<sub>22</sub>,....,k<sub>66</sub>);(1)结构矩阵M<sub>rel</sub>存储模块(102);所述的结构矩阵M<sub>rel</sub>={-1<sup>k</sup>AIV<sub>ij</sub>}<sub>6×6,k=1,2</sub>   (2)式(2)中,i代表旋翼编号,j代表自由度编号;i=1、2、3、……6,j=1、2、3、4、5、6;AIV<sub>ij</sub>为第i个旋翼动作时对第j个自由度的绝对影响值;设六旋翼无人机机体坐标系的x、y轴位于各连接杆所在的平面内并相互垂直,z轴通过六旋翼无人机的质心且垂直于各连接杆所在的平面;沿x、y、z轴正向移动自由度为正,沿x、y、z轴反向移动自由度为负;俯仰角θ为无人机机体绕y轴转动的角度,逆时针转动自由度为正,顺指针转动自由度为负;滚转角φ为无人机机体绕x轴转动的角度,逆时针转动自由度为正,顺时针转动自由度为负;偏航角ψ为无人机机体绕z轴转动的角度,逆时针转动自由度为正,顺时针转动自由度为负;则第i个旋翼动作时对x、y、z轴方向平动自由度的绝对影响值:<img file="FDA0000079786310000011.GIF" wi="1716" he="114" />第i个旋翼动作时对俯仰角θ、滚转角φ、偏航角ψ转动自由度的绝对影响值:<img file="FDA0000079786310000012.GIF" wi="1928" he="138" />当自由度向负方向变化时,k=1,当自由度向正方向变化时k=2;将系数矩阵与结构矩阵相乘并存储为映射矩阵Map的模块(103);映射矩阵降维模块(104):依次将映射矩阵Map中的某一行去掉,得到六个五行六列的行满秩矩阵Map′<sub>j</sub>(j=,1,2,3,4,5,6);降维齐次线性方程求解模块(105):求出六个降维齐次线性方程的基础解向量<img file="FDA0000079786310000013.GIF" wi="97" he="82" />所述降维齐次线性方程为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>Map</mi><mi>j</mi><mo>&prime;</mo></msubsup><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>1</mn><mi>j</mi></mrow><mn>2</mn></msubsup></mtd><mtd><msub><mi>&Delta;&Omega;</mi><mrow><mn>2</mn><mi>j</mi></mrow></msub></mtd><mtd><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>3</mn><mi>j</mi></mrow><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>4</mn><mi>j</mi></mrow><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>5</mn><mi>j</mi></mrow><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>6</mn><mi>j</mi></mrow><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>1</mn><mi>j</mi></mrow><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>2</mn><mi>j</mi></mrow><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>3</mn><mi>j</mi></mrow><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>4</mn><mi>j</mi></mrow><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>5</mn><mi>j</mi></mrow><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&Delta;&Omega;</mi><mrow><mn>6</mn><mi>j</mi></mrow><mn>2</mn></msubsup></mrow></math>]]></maths>分别为自由度虚拟控制量增量ΔU<sub>j</sub>引起的各旋翼转速平方增量;基础解向量归一化模块(106):将基础解向量<img file="FDA0000079786310000022.GIF" wi="63" he="82" />进行归一化运算得到<img file="FDA0000079786310000023.GIF" wi="94" he="172" />然后将<img file="FDA0000079786310000024.GIF" wi="72" he="172" />与自由度虚拟控制量增量ΔU<sub>j</sub>对应的映射矩阵的行a<sub>j</sub>相乘得到<img file="FDA0000079786310000025.GIF" wi="185" he="184" />最后将<img file="FDA0000079786310000026.GIF" wi="72" he="172" />送给求和模块(109),将<img file="FDA0000079786310000027.GIF" wi="160" he="184" />送给放大增益模块(108);差分模块(107):根据两个采样周期分别计算得到的自由度虚拟控制量<img file="FDA0000079786310000028.GIF" wi="78" he="58" />和<img file="FDA0000079786310000029.GIF" wi="74" he="58" />计算自由度虚拟控制量增量ΔU<sub>j</sub>,其中,t<sub>2</sub>=t<sub>1</sub>+Δt;放大增益模块(108):利用式(12)计算得到放大增益K<sub>j</sub>,然后将K<sub>j</sub>与<img file="FDA00000797863100000210.GIF" wi="72" he="172" />相乘,并将相乘的结果送到求和模块(109);<maths num="0003"><![CDATA[<math><mrow><msub><mi>&Delta;U</mi><mi>j</mi></msub><mo>=</mo><msub><mi>K</mi><mi>j</mi></msub><msub><mrow><mo>|</mo><mo>|</mo><msub><mover><mi>a</mi><mi>r</mi></mover><mi>j</mi></msub><mfrac><msub><mover><mi>&xi;</mi><mi>r</mi></mover><mi>j</mi></msub><mrow><mo>|</mo><msub><mover><mi>&xi;</mi><mi>r</mi></mover><mi>j</mi></msub><mo>|</mo></mrow></mfrac><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>求和模块(109):利用式(13)计算出各旋翼转速总平方增量;<img file="FDA00000797863100000212.GIF" wi="946" he="172" />非线性解算模块(110):利用式(14)得到各旋翼转速<maths num="0004"><![CDATA[<math><mrow><msub><mi>&Omega;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>&Delta;t</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msubsup><mi>&Omega;</mi><mi>i</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><msub><mi>&Delta;&Omega;</mi><mi>i</mi></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FDA00000797863100000214.GIF" wi="119" he="63" />为当前时刻第i个旋翼转速的平方,Ω<sub>i</sub>(t+Δt)为下一时刻第i个旋翼的转速。
地址 130033 吉林省长春市东南湖大路3888号
您可能感兴趣的专利