发明名称 一种确定多状态跟踪制导参数的方法
摘要 本发明公开了一种确定多状态跟踪制导参数的方法,该方法包括:建立飞行器再入飞行过程的多状态运动模型,并对所述多状态运动模型进行小偏差线性化处理,得到处理后的线性化方程;基于所述线性化方程,得到制导方程;基于线性二次调节器LQR,设计LQR跟踪控制器,得到相应的反馈控制律;根据多状态跟踪要求,选取控制器加权矩阵;根据所述控制器加权矩阵和制导方程,计算得到制导参数。通过使用本发明所提供的方法,可以根据所确定的制导参数实现对多个状态量的跟踪控制,降低多个状态量跟踪时的相互影响。
申请公布号 CN105867399B 申请公布日期 2017.05.03
申请号 CN201610244398.2 申请日期 2016.04.18
申请人 北京航天自动控制研究所 发明人 黄万伟;杨业;包为民;马卫华;祁振强;吴浩;刘毅;郭涛;梁禄扬;徐国强
分类号 G05D1/08(2006.01)I;G05D1/04(2006.01)I;G05D1/10(2006.01)I;G05D13/62(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 北京君恒知识产权代理事务所(普通合伙) 11466 代理人 张璐;黄启行
主权项 一种确定多状态跟踪制导参数的方法,其特征在于,该方法包括如下步骤:建立飞行器再入飞行过程的多状态运动模型,并对所述多状态运动模型进行小偏差线性化处理,得到处理后的线性化方程;基于所述线性化方程,得到制导方程;基于线性二次调节器LQR,设计LQR跟踪控制器,得到相应的反馈控制律;根据多状态跟踪要求,选取控制器加权矩阵;根据所述控制器加权矩阵和制导方程,计算得到制导参数;其中,所述建立飞行器再入飞行过程的多状态运动模型包括:根据纵向运动参数高度h、地心矢径r、速度v,当地弹道倾角θ,以及横向位置Z、航迹偏角ψ,建立飞行器再入飞行过程的有量纲的多状态运动模型;以及所述处理后的线性化方程为:<maths num="0001"><math><![CDATA[<mrow><mo>&lsqb;</mo><mi>&delta;</mi><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>&rsqb;</mo><mo>=</mo><mo>&lsqb;</mo><mi>A</mi><mo>&rsqb;</mo><mo>&lsqb;</mo><mi>&delta;</mi><mi>x</mi><mo>&rsqb;</mo><mo>+</mo><mo>&lsqb;</mo><mi>B</mi><mo>&rsqb;</mo><mo>&lsqb;</mo><mi>&delta;</mi><mi>u</mi><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0001226761100000011.GIF" wi="548" he="71" /></maths>其中,[δx]=[δr,δv,δθ,δψ,δZ]<sup>T</sup>=[r-r<sub>ref</sub>,v-v<sub>ref</sub>,θ-θ<sub>ref</sub>,ψ-ψ<sub>ref</sub>,Z-Z<sub>ref</sub>]<sup>T</sup>,δu=[δα,δγ<sub>V</sub>]<sup>T</sup>,A和B是由标准轨迹中特征点微分方程组的偏导数项组成的时变矩阵,r<sub>ref</sub>为参考地心矢径,v<sub>ref</sub>为参考速度,θ<sub>ref</sub>为参考弹道倾角,ψ<sub>ref</sub>为参考航迹偏角,Z<sub>ref</sub>为参考横向位置,α为攻角,γ<sub>V</sub>为倾侧角;以及所述制导方程为:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>&alpha;</mi><mrow><mi>c</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>&alpha;</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>+</mo><mi>&delta;</mi><mi>&alpha;</mi><mo>=</mo><msub><mi>&alpha;</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mi>&delta;</mi><mi>r</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mi>&delta;</mi><mi>v</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>3</mn></mrow></msub><mi>&delta;</mi><mi>&theta;</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>4</mn></mrow></msub><mi>&delta;</mi><mi>&psi;</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>5</mn></mrow></msub><mi>&delta;</mi><mi>Z</mi></mtd></mtr><mtr><mtd><msub><mi>&gamma;</mi><mrow><mi>V</mi><mi>c</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>&gamma;</mi><mrow><mi>V</mi><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>+</mo><mi>&delta;</mi><msub><mi>&gamma;</mi><mi>V</mi></msub><mo>=</mo><msub><mi>&gamma;</mi><mrow><mi>V</mi><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mi>&delta;</mi><mi>r</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mi>&delta;</mi><mi>v</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub><mi>&delta;</mi><mi>&theta;</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>4</mn></mrow></msub><mi>&delta;</mi><mi>&psi;</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>5</mn></mrow></msub><mi>&delta;</mi><mi>Z</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001226761100000012.GIF" wi="1414" he="159" /></maths>其中,α<sub>ref</sub>和γ<sub>V ref</sub>分别为预设的标准轨迹对应的参考攻角和参考倾侧角,α<sub>cx</sub>为当前攻角,γ<sub>V cx</sub>为当前倾侧角;K<sub>a1</sub>~K<sub>a5</sub>分别为攻角指令的第一至第五增益系数,K<sub>s1</sub>~K<sub>s5</sub>分别为倾侧角指令的第一至第五反馈系数,共同组成2×5的增益矩阵K<sub>d</sub>。
地址 100854 北京市海淀区永定路50号
您可能感兴趣的专利