发明名称 一种无需加速度传感器组合补偿的无人机升降率测量方法
摘要 本发明公开了一种无需加速度传感器组合补偿的无人机升降率测量方法。所述的升降率测量方法主要是针对当前基于组合补偿原理的无人机升降率测量方法存在对加速度传感器信息过分依赖的弊端,一旦加速度传感器测量精度下降或者失效,此时将难以提供具有相当精度的升降率信息。基于此,本发明构造了一种仅依赖于无人机高度信息且具有级联形式的非线性跟踪微分器,同时采用“先预报,后微分”的原理来弥补升降率的相位损失,以此来解决缺少加速度传感器配置或失效情况下的无人机升降率的测量问题。本发明可用于无人机的自动控制飞行,特别是在加速度传感器失效等场合,将有利于提高飞行安全系数和高度控制精度。
申请公布号 CN103759761B 申请公布日期 2016.04.06
申请号 CN201410020172.5 申请日期 2014.01.16
申请人 北京航空航天大学 发明人 王宏伦;邵星灵;苏子康;向锦武
分类号 G01D21/00(2006.01)I;G05D1/00(2006.01)I 主分类号 G01D21/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文颖
主权项 一种无需加速度传感器组合补偿的无人机升降率测量方法,包括以下几个步骤:步骤一:由无线电高度表获取当前无人机的高度信息H(k);步骤二:将当前无人机的高度信息H(k)送入第一级非线性跟踪微分器,获得升降率信息z<sub>2</sub>(k);第一级非线性跟踪微分器为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>f</mi><mi>h</mi><mo>=</mo><mi>f</mi><mi>h</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><mi>H</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><msub><mi>z</mi><mn>2</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mi>r</mi><mo>,</mo><msub><mi>h</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mi>f</mi><mi>h</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000871698610000011.GIF" wi="908" he="230" /></maths>其中:fhan(x<sub>1</sub>,x<sub>2</sub>,r,h<sub>0</sub>)为最速控制综合函数,其公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>d</mi><mo>=</mo><msub><mi>rh</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mn>0</mn></msub><mo>=</mo><msub><mi>h</mi><mn>0</mn></msub><mi>d</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>h</mi><mn>0</mn></msub><msub><mi>x</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mn>0</mn></msub><mo>=</mo><msqrt><mrow><msup><mi>d</mi><mn>2</mn></msup><mo>+</mo><mn>8</mn><mi>r</mi><mo>|</mo><mi>y</mi><mo>|</mo></mrow></msqrt></mrow></mtd></mtr><mtr><mtd><mrow><mi>a</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><mo>(</mo><mrow><msub><mi>a</mi><mn>0</mn></msub><mo>-</mo><mi>d</mi></mrow><mo>)</mo></mrow><mn>2</mn></mfrac><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>,</mo><mo>|</mo><mi>y</mi><mo>|</mo><mo>&gt;</mo><msub><mi>d</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><mfrac><mi>y</mi><mi>h</mi></mfrac><mo>,</mo><mo>|</mo><mi>y</mi><mo>|</mo><mo>&le;</mo><msub><mi>d</mi><mn>0</mn></msub></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mi>h</mi><mi>a</mi><mi>n</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><mi>r</mi><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>,</mo><mo>|</mo><mi>a</mi><mo>|</mo><mo>&gt;</mo><mi>d</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>r</mi><mfrac><mi>a</mi><mi>d</mi></mfrac><mo>,</mo><mo>|</mo><mi>a</mi><mo>|</mo><mo>&le;</mo><mi>d</mi></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000871698610000012.GIF" wi="725" he="735" /></maths>式中:z<sub>1</sub>(k)为跟踪微分器对高度信息H(k)的最速跟踪,z<sub>2</sub>(k)为高度信息的近似微分,即升降率信息,记作<img file="FDA0000871698610000013.GIF" wi="143" he="78" />h为系统采样步长,r为加速度因子,h<sub>0</sub>为滤波因子;步骤三:将步骤二获得的升降率信息z<sub>2</sub>(k)送入第二级非线性跟踪微分器中,并采用先预报后微分的方法,获得升降率信息z<sub>3</sub>(k);第二级非线性跟踪微分器为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>H</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>h</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mi>h</mi><mo>=</mo><mi>f</mi><mi>h</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>3</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>H</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><msub><mi>z</mi><mn>4</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>h</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mi>f</mi><mi>h</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000871698610000014.GIF" wi="942" he="319" /></maths>式中,z<sub>3</sub>(k)为第二级跟踪微分器对第一级跟踪微分器输出的升降率信息z<sub>2</sub>(k)的最速跟踪,即最终获取的的升降率信息,z<sub>4</sub>(k)为升降率信息z<sub>3</sub>(k)的近似微分,即加速度信息;h<sub>1</sub>为预报步长,<img file="FDA0000871698610000021.GIF" wi="142" he="78" />为z<sub>2</sub>(k)的预测修正值,r<sub>1</sub>为第二级跟踪微分器的加速度因子;步骤四:飞行过程中针对无人机的高度信息H(k),重复上述步骤得到升降率信息z<sub>3</sub>(k)并送入飞控系统中。
地址 100191 北京市海淀区学院路37号