发明名称 一种基于扩展卡尔曼滤波的无人机角加速度估计方法及装置
摘要 本发明涉及基于扩展卡尔曼滤波的无人机角加速度估计方法及装置,方法包括:S1初始化,S2预测过程,S3校正过程。本发明提供的基于扩展卡尔曼滤波的角加速度估计方法,从含有噪声的低成本传感器中估计出角加速度状态量,并且针对实际工程应用中的传感器不完全测量的问题和如何提取方位角,提供了计算策略。此外,本发明提供的基于扩展卡尔曼滤波的角加速度估计方法不仅可以估计出角加速度,还可以将陀螺仪,磁罗盘和加速度计原始值进行滤波。
申请公布号 CN106597017A 申请公布日期 2017.04.26
申请号 CN201611169133.7 申请日期 2016.12.16
申请人 上海拓攻机器人有限公司 发明人 尹亮亮;龙诗科;李少斌;张羽
分类号 G01P15/14(2013.01)I;G01H17/00(2006.01)I 主分类号 G01P15/14(2013.01)I
代理机构 北京卓唐知识产权代理有限公司 11541 代理人 龚洁
主权项 一种基于扩展卡尔曼滤波的无人机角加速度估计方法,其特征在于,包括如下步骤:S1初始化,用以确定状态矢量和观测矢量;S2预测过程,根据上述状态矢量和观测矢量转化得到非线性动力学模型;基于非线性动力学模型进行初始估计,向前推算得到第一步的k+1先验状态估计:<maths num="0001"><math><![CDATA[<mrow><msubsup><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi></msub><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001183163170000011.GIF" wi="355" he="78" /></maths>向前推算出误差协方差:<img file="FDA0001183163170000012.GIF" wi="469" he="71" />其中,测量噪声矩阵Q<sub>k</sub>保持变量σ<sup>2</sup>作为状态量对角项,用以作为滤波器的变量,S3校正过程,建立观测方程,计算卡尔曼增益:<maths num="0002"><math><![CDATA[<mrow><msub><mi>K</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>H</mi><mi>k</mi></msub><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>R</mi><mi>k</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA0001183163170000013.GIF" wi="614" he="75" /></maths>由观测变量z<sub>k</sub>更新后验估计:<maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msubsup><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>K</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><msub><mi>H</mi><mi>k</mi></msub><msubsup><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001183163170000014.GIF" wi="549" he="71" /></maths>更新误差协方差:<maths num="0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msub><mi>K</mi><mi>k</mi></msub><msub><mi>H</mi><mi>k</mi></msub><mo>)</mo></mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>.</mo><mo>,</mo></mrow>]]></math><img file="FDA0001183163170000015.GIF" wi="421" he="71" /></maths>继续步骤S2,进行迭代。
地址 200000 上海市浦东新区泥城镇新城路2号24幢C1935室