主权项 |
一种减小航天动力学系统状态敏感度的方法,其特征在于:包括如下步骤,A步骤,建立系统动力学方程: <mrow> <mover> <mi>x</mi> <mo>·</mo> </mover> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>u</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>其中,x是维数为n的状态变量,u是控制变量,t是时间变量,n为自然数;B步骤,设任意t∈[t0,tf],其中t0是初始时刻,tf是终点时刻,系统状态记为x(t)=X(t|t0,x0),x0是初始时刻t0的状态;考虑系统状态敏感度的最优性能指标定义如下:J=J1+c0·J2; (2)其中,J1是最优性能指标,J2是系统状态敏感度指标;c0是惩罚因子,用于调节灵敏度指标所占的权重;C步骤,定义任意时刻t的状态x对初始时刻t0状态x0的敏感度矩阵S(t|t0)如下: <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>∂</mo> <mi>x</mi> </mrow> <mrow> <mo>∂</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>易验证敏感度矩阵S(t|t0)具有如下特性: <mrow> <mfrac> <mrow> <mo>∂</mo> <mi>S</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>∂</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>∂</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>∂</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>∂</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>∂</mo> <mi>u</mi> </mrow> </mfrac> <mfrac> <mrow> <mo>∂</mo> <mi>u</mi> </mrow> <mrow> <mo>∂</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>·</mo> <mi>S</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>S(t0|t0)=I; (5)S(t2|t1)=S(t2|t0)·S(t1|t0)-1; (6)其中,I是单位矩阵,(·)-1为求逆运算;D步骤,根据式(4)、(5)、(6)来求解式(3)所述系统状态敏感度矩阵S;E步骤,根据系统状态敏感度矩阵S定义系统状态敏感度指标J2如下:J2=∑g(Si,j(tf,t)) (7)其中g是以敏感度矩阵为自变量的函数,i,j均为自然数;F步骤,结合式(1)、式(4),采用优化算法求取式(2)所定义的考虑系统状态敏感度的最优性能指标J。 |