发明名称 一种基于角增量和比力输出的划桨速度计算方法
摘要 本发明公开了一种基于角增量和比力输出的划桨速度计算方法,包括:(1)在速度算法单次更新周期(t<sub>m-1</sub>,t<sub>m</sub>)内,采用N个间隔,分别对角增量陀螺仪和比力加速度计的输出采样,得到角增量输出为Δθ<sub>i</sub>,i=1,2,...,N;比力输出为f<sub>j</sub>,j=1,2,...N+1;(2)计算单次算法更新周期(t<sub>m-1</sub>,t<sub>m</sub>)内机体的角增量Δθ<sub>m</sub>和速度增量△V<sub>m</sub>;(3)计算旋转效应补偿值ΔV<sub>rotm</sub>;(4)计算划桨效应补偿值ΔV<sub>sculm</sub>;(5)计算(t<sub>m-1</sub>,t<sub>m</sub>)周期内比力积分的补偿值ΔV<sub>sfm</sub>;(6)根据速度更新方程计算t<sub>m</sub>时刻的机体速度V<sub>m</sub>。本发明可根据角增量陀螺仪和比力角速度计的输出直接高精度的计算出机体的速度,具有方法简单、精度高的优势。
申请公布号 CN104677356A 申请公布日期 2015.06.03
申请号 CN201510118855.9 申请日期 2015.03.16
申请人 黄磊 发明人 黄磊
分类号 G01C21/16(2006.01)I;G01C19/64(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 代理人
主权项 一种基于角增量和比力输出的划桨速度计算方法,其特征在于,包括以下步骤:(1)在速度算法单次更新周期(t<sub>m‑1</sub>,t<sub>m</sub>)内,采用N个间隔,角增量陀螺仪输出N次,比力加速度计输出N+1次,所述的N为≥2的自然数,分别对角增量陀螺仪和比力加速度计输出采样,得到角增量输出为Δθ<sub>i</sub>,i=1,2,...,N;比力输出为f<sub>j</sub>,j=1,2,...N+1;(2)计算单次算法更新周期(t<sub>m‑1</sub>,t<sub>m</sub>)内机体的角增量Δθ<sub>m</sub>和速度增量ΔV<sub>m</sub>:Δθ<sub>m</sub>=Δθ<sub>1</sub>+Δθ<sub>2</sub>+...+Δθ<sub>N</sub>,所述的Δθ<sub>m</sub>是(t<sub>m‑1</sub>,t<sub>m</sub>)周期内角增量的输出Δθ<sub>1</sub>、Δθ<sub>1</sub>...Δθ<sub>N</sub>之和,ΔV<sub>m</sub>根据比力加速度计的输出采样f<sub>j</sub>经过数字积分公式计算得到;(3)计算旋转效应补偿值Δv<sub>rotm</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Delta;V</mi><mi>rotm</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&Delta;&theta;</mi><mi>m</mi></msub><mo>&times;</mo><msub><mi>&Delta;V</mi><mi>m</mi></msub><mo>;</mo></mrow>]]></math><img file="FSA0000114809130000011.GIF" wi="423" he="111" /></maths>(4)计算划桨效应补偿值ΔV<sub>sculm</sub>:<img file="FSA0000114809130000014.GIF" wi="725" he="129" />所述的l=2,...N,k<sub>l‑1</sub>为划桨优化算法系数,H为速度算法更新周期,H=t<sub>m</sub>‑t<sub>m‑1</sub>;(5)计算(t<sub>m‑1</sub>,t<sub>m</sub>)周期内比力积分的补偿值ΔV<sub>sfm</sub>:ΔV<sub>sfm</sub>=ΔV<sub>m</sub>+ΔV<sub>rotm</sub>+ΔV<sub>sculm</sub>;(6)根据速度更新方程计算t<sub>m</sub>时刻的机体速度V<sub>m</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>V</mi><mi>m</mi></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>&Delta;V</mi><mi>sfm</mi></msub><mo>+</mo><msub><mi>&Delta;V</mi><mrow><mi>g</mi><mo>/</mo><mi>corn</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>V</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>[</mo><msub><mi>&Delta;V</mi><mi>m</mi></msub><mo>+</mo><msub><mrow><mi>&Delta;</mi><mover><mi>V</mi><mo>^</mo></mover></mrow><mi>rotm</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mover><mi>V</mi><mo>^</mo></mover><mi>sculm</mi></msub><mo>]</mo><mo>+</mo><msub><mi>&Delta;V</mi><mrow><mi>g</mi><mo>/</mo><mi>corn</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>V</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>[</mo><msub><mi>&Delta;V</mi><mi>m</mi></msub><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&Delta;&theta;</mi><mi>m</mi></msub><mo>&times;</mo><msub><mi>&Delta;V</mi><mi>m</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>2</mn></mrow><mi>N</mi></munderover><msub><mi>k</mi><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>&Delta;&theta;</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>f</mi><mi>l</mi></msub><mo>-</mo><msub><mi>&Delta;&theta;</mi><mi>l</mi></msub><mo>&times;</mo><msub><mi>f</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>H</mi><mo>]</mo><mo>+</mo><mo>[</mo><mi>g</mi><mo>-</mo><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>&omega;</mi></mrow><mi>ie</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>en</mi></msub><mo>)</mo></mrow><mo>&times;</mo><msub><mi>V</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>]</mo><mi>H</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FSA0000114809130000013.GIF" wi="1700" he="287" /></maths>所述的,V<sub>m‑1</sub>为t<sub>m‑1</sub>时刻机体速度,C<sub>m‑1</sub>为t<sub>m‑1</sub>时刻方向余弦阵,g为重力加速度,ω<sub>ie</sub>为地理系对惯性系的角速率,ω<sub>en</sub>为导航系对地理系的角速率。
地址 210037 江苏省南京市龙蟠路159号南京林业大学机电学院自动化系