发明名称 基于飞机运动方程的飞行数据“野值”检测与剔除方法
摘要 本发明涉及一种基于飞机运动方程的飞行数据“野值”检测与剔除方法,属于飞行力学和信息处理技术领域,用于检查飞行数据是否存在与其它测量值相差很大“野值”。为了克服现有技术方法仅仅依赖于自身测量数据而导致的“野值”误判断,本发明提出了一种基于飞机运动方程的逐步预测“野值”阈值判定和剔除方法,该方法根据状态的变化快慢,根据有效简化的输出响应来逐步预测飞行状态和参数,每步按照设置好的阈值判定和剔除测量中“野值”;由于在实际测量中快变状态或参数存在“野值”可能性比慢变量大,而且输出响应具有工程上的精度,本发明提出的方法对剔除飞行数据“野值”更有效。<pb pnum="1" />
申请公布号 CN106508018B 申请公布日期 2013.07.24
申请号 CN201010047879.7 申请日期 2010.03.18
申请人 西北工业大学 发明人 史忠科
分类号 G06F17/00(2006.01)I;G05D1/00(2006.01)I 主分类号 G06F17/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于飞机运动方程的飞行数据“野值”检测与剔除方法,其特征在于包括以下步骤:(a)在进行速度“野值”剔除时,令两个采样点之间的速度增量预测值为:ΔV=gT{n<sub>x</sub>(k)cos[α(k)]+n<sub>z</sub>(k)sin[α(k)]}当ΔV≥0,V<sub>0</sub>(k)+0.8ΔV≤V<sub>0</sub>(k+1)≤V<sub>0</sub>(k)+1.2ΔV时无“野值”;当ΔV<0,V<sub>0</sub>(k)+0.8ΔV≥V<sub>0</sub>(k+1)≥V<sub>0</sub>(k)+1.2ΔV时无“野值”;如果速度增量的实测值超出预测值过多,则需要进行“野值”剔除,否则认为没有“野值”;(b)在进行气流迎角“野值”剔除时,令两个采样点之间的气流迎角增量预测值为:Δα=T{q(k)+g[1+n<sub>z</sub>(k)]/V<sub>0</sub>(k)]}当Δα≥0,α(k)+0.8Δα≤α(k+1)≤α(k)+1.2Δα时无“野值”;当Δα<0,α(k)+0.8Δα≥α(k+1)≥α(k)+1.2Δα时无“野值”;如果气流迎角增量的实测值超出预测值过多,则需要进行“野值”剔除,否则认为没有“野值”;(c)在进行侧滑角“野值”剔除时,令两个采样点之间的侧滑角增量预测值为:<maths id="cmaths0001" num="0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>&beta;</mi><mo>=</mo><mi>T</mi><mo>{</mo><mfrac><mn>1</mn><mrow><msub><mi>V</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>gn</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>r</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mo>&lsqb;</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mi>p</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>s</mi><mi>i</mi><mi>n</mi><mo>&lsqb;</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow>]]></math><img file="FWW0000000051980000011.GIF" wi="1010" he="113" /></maths>当Δβ≥0,β(k)+0.8Δβ≤β(k+1)≤β(k)+1.2Δβ时无“野值”;当Δβ<0,β(k)+0.8Δβ≥β(k+1)≥β(k)+1.2Δβ时无“野值”;如果侧滑角增量的实测值超出预测值过多,则需要进行“野值”剔除,否则认为没有“野值”;(d)在进行三个欧拉角“野值”剔除时,令两个采样点之间的俯仰角增量预测值为:<img file="FWW0000000051980000012.GIF" wi="790" he="54" />当Δθ≥0,θ(k)+0.8Δθ≤θ(k+1)≤θ(k)+1.2Δθ时无“野值”;当Δθ<0,θ(k)+0.8Δθ≥θ(k+1)≥θ(k)+1.2Δθ时无“野值”;令两个采样点之间的滚转角增量预测值为:<img file="FWW0000000051980000013.GIF" wi="1221" he="54" />当<img file="FWW0000000051980000014.GIF" wi="953" he="53" />时无“野值”;当<img file="FWW0000000051980000015.GIF" wi="954" he="54" />时无“野值”;令两个采样点之间的偏航角增量预测值为:<img file="FWW0000000051980000016.GIF" wi="1019" he="55" /><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>u</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><mn>7</mn><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mn>8</mn><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mn>12</mn><mi>T</mi></mrow></mfrac></mrow></mtd><mtd><mrow><mi>u</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>V</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mo>&lsqb;</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo>&lsqb;</mo><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><mn>7</mn><mi>v</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mn>8</mn><mi>v</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mn>12</mn><mi>T</mi></mrow></mfrac></mrow></mtd><mtd><mrow><mi>v</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>V</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>sin</mi><mo>&lsqb;</mo><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>w</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><mn>7</mn><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mn>8</mn><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mn>12</mn><mi>T</mi></mrow></mfrac></mrow></mtd><mtd><mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>V</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>sin</mi><mo>&lsqb;</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo>&lsqb;</mo><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FWW0000000051980000028.GIF" wi="1737" he="399" /></maths>采样点处三个加速度预测值为:<img file="FWW0000000051980000029.GIF" wi="1216" he="215" />分别令<img file="FWW00000000519800000210.GIF" wi="504" he="60" />及实测值<img file="FWW00000000519800000211.GIF" wi="561" he="60" />当<img file="FWW00000000519800000212.GIF" wi="155" he="47" />时,<img file="FWW00000000519800000213.GIF" wi="493" he="51" />时无“野值”;当<img file="FWW00000000519800000214.GIF" wi="155" he="47" />时,<img file="FWW00000000519800000215.GIF" wi="493" he="52" />时无“野值”;如果三个加速度实测值超出预测值过多,则需要进行“野值”剔除,否则认为没有“野值”;以上式中,T是采样周期;g是重力加速度。pb pnum="2" />
地址 710072 陕西省西安市友谊西路127号