发明名称 制导系统陀螺漂移的Q-SVR预测方法
摘要 本发明属于制导领域,具体涉及一种制导系统陀螺漂移的快速支持向量回归(Q-SVR)预测方法。主要应用于航空航天领域制导系统,目的是有效降低SVR法预测陀螺漂移的计算复杂度,依据历史漂移数据,快速预测陀螺漂移值。本预测方法,归一化陀螺漂移模型的拉格朗日乘子、惩罚因子、不敏感参数以及核参数为单一变量;将约束条件作为惩罚项重构入目标函数中;采用梯度下降方式单层迭代计算目标函数的更新方向和更新步长,给出最优值陀螺漂移预测值。简化了陀螺漂移模型,方便计算,且提升约2倍计算效率。
申请公布号 CN102538771B 申请公布日期 2014.07.23
申请号 CN201110332826.4 申请日期 2011.10.28
申请人 中国航天科工集团第三研究院第八三五七研究所 发明人 贾磊
分类号 G01C19/00(2013.01)I;G06F19/00(2011.01)I 主分类号 G01C19/00(2013.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 一种制导系统陀螺漂移的Q‑SVR预测方法,归一化陀螺漂移模型的拉格朗日乘子、惩罚因子、不敏感参数以及核参数为单一变量;将约束条件作为惩罚项重构入目标函数中;采用梯度下降方式单层迭代计算目标函数的更新方向和更新步长,给出最优值陀螺漂移预测值;具体包括如下步骤:步骤一:输入历史漂移数据(x<sub>i</sub>,y<sub>i</sub>),...,(x<sub>l</sub>,y<sub>l</sub>),预设收敛容忍度ω,障碍因子β,r<sub>0</sub>,变量初值X<sub>0</sub>;步骤二:计算初始目标函数值J(X<sub>0</sub>,r<sub>0</sub>),其中<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>l</mi></munderover><mrow><mo>(</mo><msub><mover><mi>&alpha;</mi><mo>^</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mover><mi>&alpha;</mi><mo>^</mo></mover><mi>j</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>)</mo></mrow><mover><mi>K</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>i</mi></mrow><mi>l</mi></munderover><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mover><mi>&alpha;</mi><mo>^</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>&epsiv;</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>l</mi></munderover><mrow><mo>(</mo><msub><mover><mi>&alpha;</mi><mo>^</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><msub><mi>r</mi><mi>k</mi></msub></mfrac><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>T</mi></msup><mover><mi>Y</mi><mo>~</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>r</mi><mi>k</mi></msub><mfrac><mn>1</mn><mi>X</mi></mfrac></mrow>]]></math><img file="FDA0000474627600000011.GIF" wi="1874" he="148" /></maths>其中:<img file="FDA0000474627600000012.GIF" wi="66" he="73" />α<sub>i</sub>为拉格朗日乘子,i=1,2,...,l其中,l为历史漂移数据个数;<img file="FDA0000474627600000013.GIF" wi="190" he="97" />为修正核函数,<img file="FDA0000474627600000014.GIF" wi="610" he="97" />K(x<sub>i</sub>,x<sub>j</sub>)为核函数,σ<sub>ij</sub>、C为惩罚因子,<img file="FDA0000474627600000015.GIF" wi="191" he="97" />的参数表示为Θ=(C,θ<sub>1</sub>,…,θ<sub>d</sub>)<sup>T</sup>,θ<sub>1</sub>,…,θ<sub>d</sub>为核参数值,d为核参数个数;目标变量<img file="FDA0000474627600000016.GIF" wi="450" he="96" />ε为不敏感度参数,α<sub>A</sub>=(α<sub>1</sub>,…,α<sub>l</sub>)<sup>T</sup>,<img file="FDA0000474627600000017.GIF" wi="380" he="82" />变量<img file="FDA0000474627600000018.GIF" wi="996" he="83" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>Y</mi><mo>~</mo></mover><mo>=</mo><msup><mrow><mo>(</mo><msup><msub><mi>Y</mi><mn>1</mn></msub><mi>T</mi></msup><mo>,</mo><msup><msub><mi>Y</mi><mn>2</mn></msub><mi>T</mi></msup><mo>,</mo><msup><msub><mi>T</mi><mn>3</mn></msub><mi>T</mi></msup><mo>)</mo></mrow><mi>T</mi></msup><mo>;</mo></mrow>]]></math><img file="FDA0000474627600000019.GIF" wi="375" he="96" /></maths>r<sub>k</sub>为障碍因子序列,满足{r<sub>k</sub>|r<sub>0</sub>>0,r<sub>k+1</sub>=β,0<β<1};步骤三:最小化目标函数J(X,r),求解最优值;步骤四:若有当前目标函数梯度值大于预设的收敛容忍度,则步骤终止,得到最优预测值。
地址 300141 天津市河北区黄纬路69号