主权项 |
一种无需加速度传感器组合补偿的无人机升降率测量方法,包括以下几个步骤:步骤一:由无线电高度表获取当前无人机的高度信息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>·</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>·</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>></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>≤</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>></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>≤</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>·</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>·</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>·</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>·</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>·</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)并送入飞控系统中。 |