主权项 |
一种确定多状态跟踪制导参数的方法,其特征在于,该方法包括如下步骤:建立飞行器再入飞行过程的多状态运动模型,并对所述多状态运动模型进行小偏差线性化处理,得到处理后的线性化方程;基于所述线性化方程,得到制导方程;基于线性二次调节器LQR,设计LQR跟踪控制器,得到相应的反馈控制律;根据多状态跟踪要求,选取控制器加权矩阵;根据所述控制器加权矩阵和制导方程,计算得到制导参数;其中,所述建立飞行器再入飞行过程的多状态运动模型包括:根据纵向运动参数高度h、地心矢径r、速度v,当地弹道倾角θ,以及横向位置Z、航迹偏角ψ,建立飞行器再入飞行过程的有量纲的多状态运动模型;以及所述处理后的线性化方程为:<maths num="0001"><math><![CDATA[<mrow><mo>[</mo><mi>δ</mi><mover><mi>x</mi><mo>·</mo></mover><mo>]</mo><mo>=</mo><mo>[</mo><mi>A</mi><mo>]</mo><mo>[</mo><mi>δ</mi><mi>x</mi><mo>]</mo><mo>+</mo><mo>[</mo><mi>B</mi><mo>]</mo><mo>[</mo><mi>δ</mi><mi>u</mi><mo>]</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>α</mi><mrow><mi>c</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>α</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>+</mo><mi>δ</mi><mi>α</mi><mo>=</mo><msub><mi>α</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>δ</mi><mi>r</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mi>δ</mi><mi>v</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>3</mn></mrow></msub><mi>δ</mi><mi>θ</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>4</mn></mrow></msub><mi>δ</mi><mi>ψ</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>a</mi><mn>5</mn></mrow></msub><mi>δ</mi><mi>Z</mi></mtd></mtr><mtr><mtd><msub><mi>γ</mi><mrow><mi>V</mi><mi>c</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>γ</mi><mrow><mi>V</mi><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>+</mo><mi>δ</mi><msub><mi>γ</mi><mi>V</mi></msub><mo>=</mo><msub><mi>γ</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>δ</mi><mi>r</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mi>δ</mi><mi>v</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub><mi>δ</mi><mi>θ</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>4</mn></mrow></msub><mi>δ</mi><mi>ψ</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>s</mi><mn>5</mn></mrow></msub><mi>δ</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>。 |