专利摘要

本发明涉及一种无拖曳卫星相对位移通道的精细抗干扰滤波方法，针对含有大气阻力、太阳光压、执行机构噪声、量测噪声以及模型不确定性项等多源干扰的无拖曳卫星的相对位移通道；首先，针对无拖曳卫星相对位移通道所受多源干扰的数学表征，将上述多源干扰分类并建模；其次，建立含有多源干扰的无拖曳卫星相对位移通道的数学模型；再次，构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器；最后，基于凸优化算法求解精细抗干扰滤波器的增益矩阵；本发明具有抗干扰能力强、滤波精度高等优点，可用于无拖曳卫星相对位移通道中。

权利要求

1.一种无拖曳卫星相对位移通道的精细抗干扰滤波方法，其特征在于包括以下步骤：首先，针对无拖曳卫星相对位移通道所受多源干扰的数学表征，将多源干扰分类并建模；其次，建立含有多源干扰的无拖曳卫星相对位移通道的数学模型；再次，构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器；最后，基于凸优化算法求解精细抗干扰滤波器的增益矩阵；具体步骤如下：

(1)无拖曳系统相对位移通道所受的多源干扰分类与建模：

大气阻力是无拖曳卫星相对位移通道面临的主要干扰，对大气阻力建模并将其描述成状态空间模型的形式：

$\left\{\begin{array}{c}\stackrel{·}{\tau }\left(t\right)=M\tau \left(t\right)+\delta \left(t\right)\\ d_0\left(t\right)=N\tau \left(t\right)\end{array}\right.$

其中，d₀(t)为大气阻力干扰，τ(t)为大气阻力数学模型的状态变量，M和N为系数矩阵，δ(t)为模型不确定性项；

太阳光压和模型不确定性项δ(t)的L₂范数均有界，将其描述为能量有界干扰，即 其中d₁(t)表示太阳光压和模型不确定性项；

执行机构噪声w(t)及量测噪声v(t)均为高斯零均值不相关白噪声，其统计特性满足：

E{w(k)w^T(j)}＝Q_kδ_k-j

E{v(k)v^T(j)}＝R_kδ_k-j

E{w(k)v^T(j)}＝0

其中，E{*}表示*的数学期望，k、j为t的任意两个取值，t∈R；

其中，Q_k,R_k分别表示k时刻执行机构噪声和量测噪声的协方差矩阵；δ_k-j是Kronecker-δ函数，即如果k＝j则δ_k-j＝1，否则δ_k-j＝0；

(2)建立含有多源干扰的无拖曳卫星相对位移通道的数学模型：

$\left\{\begin{array}{c}\stackrel{·}{x}\left(t\right)=Ax\left(t\right)+Ff\left(x\right(t),t)+B\left[u\right(t)+w(t)+d_0(t)+d_1(t\left)\right]\\ y\left(t\right)=Cx\left(t\right)+v\left(t\right)\end{array}\right.$

其中，状态变量 x₁(t)＝r(t)， r(t)为无拖曳卫星的质心与检测质量块的质心的相对位移，系数矩阵 C＝[I_3×30_3×3]， m_sc为无拖曳卫星质量，m_tm为检测质量块质量，ω₀为轨道角速度，K_trans为无拖曳卫星与内部检测质量块之间耦合的水平弹性系数，D_trans为无拖曳卫星与内部检测质量块之间耦合的水平阻尼系数，u(t)为无拖曳卫星所受到的控制力，w(t)和v(t)分别为执行机构噪声和量测噪声， 为已知非线性项，其中，ω(t)为卫星的绝对角速度，ω^×(t)为相应的叉乘矩阵；假设非线性项满足Lipschitz条件，即对任意的两个系统状态x₁(t)，x₂(t)，存在已知的矩阵U使得下列不等式成立：

||f(x₁(t),t)-f(x₂(t),t)||≤||U(x₁(t)-x₂(t))||

(3)建立增广的含有多源干扰的无拖曳卫星相对位移通道的数学模型：

$\left\{\begin{array}{c}\stackrel{·}{\stackrel{‾}{x}}\left(t\right)=\bar{A}\bar{x}\left(t\right)+\bar{F}f\left(\bar{x}\right(t),t)+\bar{B}\left[u\right(t)+w(t)+d_1(t\left)\right]+\bar{G}\delta \left(t\right)\\ \bar{y}\left(t\right)=\bar{C}\bar{x}\left(t\right)+v\left(t\right)\end{array}\right.$

其中，状态变量 量测输出 系数矩阵

(4)构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器：

$\left\{\begin{array}{c}\stackrel{\stackrel{·}{^}}{\bar{x}}\left(t\right)=\bar{A}\hat{\bar{x}}\left(t\right)+\bar{F}f\left(\hat{\bar{x}}\right(t),t)+\bar{B}u\left(t\right)+L\left[\bar{y}\right(t)-\hat{\bar{y}}(t\left)\right]\\ \hat{\bar{y}}\left(t\right)=\bar{C}\hat{\bar{x}}\left(t\right)\end{array}\right.$

其中， 为 的估计值， 为 的估计值，L为针对无拖曳卫星相对位移通道的精细抗干扰滤波器的增益矩阵；

(5)基于凸优化算法求解精细抗干扰滤波器的增益矩阵L：

${\mathrm{min\gamma }}_3$$\begin{array}{c}\left[\begin{array}{cccccc}sym(P\bar{A}-P_L\bar{C})& P\bar{F}& P\bar{B}& P\bar{G}& U^T& T_1^T\\ *& -\frac{1}{{\lambda }^2}I& 0& 0& 0& 0\\ *& *& -{\gamma }_1^{2}I& 0& 0& 0\\ *& *& *& -{\gamma }_2^{2}I& 0& 0\\ *& *& *& *& -{\lambda }^{2}I& 0\\ *& *& *& *& *& -I\end{array}\right]<0\\ \left[\begin{array}{ccc}sym(P\bar{A}-P_L\bar{C})& \bar{B}& L\\ *& -I& 0\\ *& *& -I\end{array}\right]<0\\ \left[\begin{array}{cc}-Z& T_{2}P\\ *& -P\end{array}\right]<0\\ Trace\left(Z\right)<{\gamma }_3^2\end{array}$

其中，T₁和T₂为相应的加权矩阵，P和Z为正定对称阵， 为增广的无拖曳系统模型的系数矩阵，U为满足Lipschitz条件的已知矩阵，λ为任意正实数，γ₁为能量有界干扰d₁(t)到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数的H_∞范数的上限，γ₂为能量有界干扰δ(t)到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数的H_∞范数的上限，γ₃为高斯白噪声w(t)和v(t)到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数的H₂范数；求解矩阵不等式可得未知矩阵P_L，针对无拖曳卫星的精细抗干扰滤波器的增益阵L＝P^-1P_L。

说明书

技术领域

本发明涉及一种无拖曳卫星相对位移通道的精细抗干扰滤波方法，可用于无拖曳卫星相对位移通道的抗干扰滤波。

背景技术

无拖曳卫星是一种采用无拖曳技术的卫星，通过控制卫星本体和检测质量块达到高宁静的运动特性，能够达到航天领域对超静平台的要求；根据无拖曳卫星的原理，可将其数学模型分为卫星本体姿态通道、相对位移通道以及检查质量块姿态通道。其中，无拖曳卫星相对位移通道设计的基本思想是使无拖曳卫星跟随内部检测质量块运动，将无拖曳卫星相对位移通道受到的多源干扰尽可能抵消，从而使无拖曳卫星运行在接近纯引力作用下的轨道上；而无拖曳卫星运行在近地轨道中，面临大气阻力、太阳光压、执行机构噪声、量测噪声以及模型不确定性项等多源干扰，并且由于多源干扰的影响致使无拖曳卫星相对位移通道的无法达到其高控制精度的要求。

目前专门针对无拖曳卫星相对位移通道的滤波方法的研究较少，被采用过的滤波方法有传统的卡尔曼滤波、H_∞滤波以及基于干扰观测器的滤波。但是传统的滤波方法大多只是针对单一干扰类型，如卡尔曼滤波可以有效抑制高斯零均值不相关白噪声的干扰，H_∞滤波可以解决噪声统计特性未知的非线性系统滤波问题，基于干扰观测器的滤波方法可以有效估计并补偿已知干扰模型的外部干扰。而对于含有多源干扰的无拖曳卫星相对位移通道，仅使用针对单一干扰的类型的滤波方法，使其滤波精度不高，达不到无拖曳卫星相对位移通道的高精度要求。因此，目前的针对无拖曳卫星相对位移通道的滤波方法存在只能针对单一干扰进行滤波，滤波精度较低，保守性较大的问题。

发明内容

本发明的技术解决问题是：克服现有的针对无拖曳卫星相对位移通道的滤波方法只能针对单一干扰进行滤波，滤波精度较低，保守性较大的问题，设计一种无拖曳卫星相对位移通道的精细抗干扰滤波方法。

本发明的技术解决方案为：针对含有大气阻力、太阳光压、执行机构噪声、量测噪声以及模型不确定性项等多源干扰的无拖曳卫星相对位移通道。首先，针对无拖曳卫星相对位移通道所受多源干扰的数学表征并将多源干扰分类并建模；其次，建立含有多源干扰的无拖曳卫星相对位移通道的数学模型；再次，构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器；最后，基于凸优化算法求解精细抗干扰滤波器的增益矩阵。具体实施步骤如下：

(1)无拖曳系统相对位移通道所受的多源干扰分类与建模：

大气阻力是无拖曳卫星相对位移通道面临的主要干扰，可对大气阻力建模并将其描述成状态空间模型的形式：

τ·(t)=Mτ(t)+δ(t)d0(t)=Nτ(t)]]>

其中，d₀(t)为大气阻力干扰，τ(t)为大气阻力数学模型的状态变量，M和N为系数矩阵，δ(t)为模型不确定性项；

太阳光压和模型不确定性项δ(t)的L₂范数均有界，将其描述为能量有界干扰，即 其中d₁(t)表示太阳光压和模型不确定性项；

执行机构噪声w(t)及量测噪声v(t)均为高斯零均值不相关白噪声，其统计特性满足：

E{w(k)w^T(j)}＝Q_kδ_k-j

E{v(k)v^T(j)}＝R_kδ_k-j

E{w(k)v^T(j)}＝0

其中，Q_k,R_k分别表示k时刻执行机构噪声和量测噪声的协方差矩阵；δ_k-j是Kronecker-δ函数，即如果k＝j则δ_k-j＝1，否则δ_k-j＝0；

(2)建立含有多源干扰的无拖曳卫星相对位移通道的数学模型：

x·(t)=Ax(t)+Ff(x(t),t)+B[u(t)+w(t)+d0(t)+d1(t)]y(t)=Cx(t)+v(t)]]>

其中，状态变量 x₁(t)＝r(t)， r(t)为无拖曳卫星的质心与检测质量块的质心的相对位移，系数矩阵 m_sc为无拖曳卫星质量，m_tm为检测质量块质量，ω₀为轨道角速度，K_trans为无拖曳卫星与内部检测质量块之间耦合的水平弹性系数，D_trans为无拖曳卫星与内部检测质量块之间耦合的水平阻尼系数，u(t)为无拖曳卫星所受到的控制力，w(t)和v(t)分别为执行机构噪声和量测噪声， 为已知非线性项，其中，ω(t)为卫星的绝对角速度，ω×(t)为相应的叉乘矩阵；假设非线性项满足Lipschitz条件，即对任意的两个系统状态x₁(t)，x₂(t)，存在已知的矩阵U使得下列不等式成立：

||f(x₁(t),t)-f(x₂(t),t)||≤||U(x₁(t)-x₂(t))||

(3)建立增广的含有多源干扰的无拖曳卫星相对位移通道的数学模型：

x‾·(t)=A‾x‾(t)+F‾f(x‾(t),t)+B‾[u(t)+w(t)+d1(t)]+G‾δ(t)y‾(t)=C‾x‾(t)+v(t)]]>

其中，状态变量 量测输出 系数矩阵

(4)构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器：

x‾^·(t)=A‾x‾^(t)+F‾f(x‾^(t),t)+B‾u(t)+L[y‾(t)-y‾^(t)]y‾^(t)=C‾x‾^(t)]]>

其中， 为 的估计值， 为 的估计值，L为针对无拖曳卫星相对位移通道的精细抗干扰滤波器的增益矩阵；

(5)基于凸优化算法求解精细抗干扰滤波器的增益矩阵L：

minγ₃

sym(PA‾-PLC‾)PF‾PB‾PG‾UTT1T*-1λ2I0000**-γ12I000***-γ22I00****-λ2I0*****-I<0]]>

sym(PA‾-PLC‾)B‾L*-I0**-I<0]]>

-ZT2P*-P<0]]>

Trace(Z)<γ32]]>

其中，min表示最小值，T₁和T₂为相应的加权矩阵，P和Z为正定对称阵，λ为任意正实数，γ₁为能量有界干扰d₁(t)到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数的H_∞范数的上限，γ₂为能量有界干扰δ(t)到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数的H_∞范数的上限，γ₃为高斯白噪声w(t)和v(t)到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数的H₂范数；求解矩阵不等式可得未知矩阵P_L，针对无拖曳卫星的精细抗干扰滤波器的增益阵L＝P^-1P_L。

本发明与现有技术相比的优点在于：

(1)本发明充分考虑了无拖曳卫星相对位移通道含有的大气阻力、太阳光压、执行机构噪声、量测噪声以及模型不确定性项等多源干扰，并进行分析建模，建立了含有多源干扰的无拖曳卫星相对位移通道的数学模型；克服了现有无拖曳卫星相对位移通道的数学模型未充分考虑多源干扰的问题；

(2)本发明构造了针对无拖曳卫星相对位移通道的精细抗干扰滤波器，结合干扰观测器的思想，通过凸优化算法求解精细抗干扰滤波器的增益，利用H₂范数和H_∞范数，保证滤波器的鲁棒性的同时可以有效补偿和抑制多源干扰；克服了卡尔曼滤波、H_∞滤波以及基于干扰观测器的滤波等方法只能针对单一干扰进行滤波，滤波精度较低，保守性较大的问题。

附图说明

图1为针对无拖曳卫星相对位移通道的精细抗干扰滤波方法的流程框图。

具体实施方式

下面结合附图及实施例对本发明进一步详细说明。

如图1所示，本发明具体实现步骤如下：

(1)将无拖曳系统相对位移通道所受的多源干扰分类并建模：

由于无拖曳卫星运行在近地轨道中，大气阻力是无拖曳卫星相对位移通道面临的主要干扰，通过正弦曲线模拟大气阻力，可建立大气阻力的数学模型：

d0(t)=ax+Axsin(ωdt+φx)+Δxay+Aysin(ωdt+φy)+Δyaz+Azsin(ωdt+φz)+Δz]]>

其中，d₀(t)为大气阻力，a_x,a_y,a_z，A_x,A_y,A_z为未知的幅值，φ_x,φ_y,φ_z为未知的相位，ω_d为大气阻力干扰的频率，Δ_x,Δ_y,Δ_z为模型误差，再将其描述成状态空间模型的形式：

τ·(t)=Mτ(t)+δ(t)d0(t)=Nτ(t)]]>

其中，τ(t)为大气阻力模型的状态变量，M和N为系数矩阵，δ(t)为模型不确定性项； N＝[111000]，δ(t)为模型不确定性项

太阳光压和模型不确定性项δ(t)的L₂范数均有界，将其描述为能量有界干扰 其中d₁(t)表示太阳光压和模型不确定性项；

执行机构噪声及量测噪声均为高斯零均值不相关白噪声，其统计特性满足：

E{w(k)w^T(j)}＝Q_kδ_k-j

E{v(k)v^T(j)}＝R_kδ_k-j

E{w(k)v^T(j)}＝0

其中，Q_k,R_k分别表示k时刻执行机构噪声和量测噪声的协方差矩阵；δ_k-j是Kronecker-δ函数，即如果k＝j则δ_k-j＝1，否则δ_k-j＝0；

(2)建立含有多源干扰的无拖曳卫星相对位移通道的数学模型：

无拖曳卫星相对位移通道设计的基本思想是使无拖曳卫星跟随内部检测质量块运动，从而使无拖曳卫星运行在接近纯引力作用下的轨道上，假设无拖曳卫星的质心与内部检测质量块空腔的质心重合，建立无拖曳卫星相对位移通道的动力学方程如下：

r··(t)+2ωh×r·(t)+ω·h×r(t)+ωh×[ωh×r(t)]=1mtm(FGtm+FSCtm)-1msc(FGsc+FCsc+FDsc+FTMsc)]]>

其中，r(t)为无拖曳卫星的质心与内部检测质量块质心的相对位移，ω_h为内部检测质量块空腔的绝对角速度，m_tm和m_sc分别为内部检测质量块和无拖曳卫星的质量，F_Gtm内部检测质量块所受的引力，F_Gsc、F_Csc和F_Dsc分别表示无拖曳卫星所受的引力、控制力以及干扰力，将无拖曳卫星与内部检测质量块之间的静电耦合当作一个弹簧阻尼系统，则无拖曳卫星与内部检测质量块之间的静电耦合F_SCtm和F_TMsc可以表示为 其中K_trans为无拖曳卫星与内部检测质量块之间耦合的水平弹性系数(量级为10^-6N/m)，D_trans为水平阻尼系数(量级为10^-11N/m)；

令x₁(t)＝r， 则状态变量 那么含有多源干扰的无拖曳卫星相对位移通道的数学模型为：

x·(t)=Ax(t)+Ff(x(t),t)+B[u(t)+w(t)+d0(t)+d1(t)]y(t)=Cx(t)+v(t)]]>

其中，状态变量 x₁(t)＝r(t)， r(t)为无拖曳卫星的质心与内部检测质量块质心的相对位移，系数矩阵分别表示为 ω₀为轨道角速度，K_trans为无拖曳卫星与内部检测质量块之间耦合的水平弹性系数，D_trans为水平阻尼系数，u(t)＝F_Csc为无拖曳卫星所受到的控制力，w(t)和v(t)分别为执行机构噪声和量测噪声， 为已知非线性项，其中，ω(t)为卫星的绝对角速度，ω×(t)为相应的叉乘矩阵；假设非线性项满足Lipschitz条件，即对任意的两个系统状态x₁(t)，x₂(t)，存在已知的矩阵U使得下列不等式成立：

||f(x₁(t),t)-f(x₂(t),t)||≤||U(x₁(t)-x₂(t))||

(3)建立增广的含有多源干扰的无拖曳卫星相对位移通道的数学模型：

即将大气阻力的数学模型与含有多源的干扰无拖曳相对位移通道的数学模型联立，得：

x‾·(t)=A‾x‾(t)+F‾f(x‾(t),t)+B‾[u(t)+w(t)+d1(t)]+G‾δ(t)y‾(t)=C‾x‾(t)+v(t)]]>

其中，状态变量 量测输出 系数矩阵

(4)构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器：

x‾^·(t)=A‾x‾^(t)+F‾f(x‾^(t),t)+B‾u(t)+L[y‾(t)-y‾^(t)]y‾^(t)=C‾x‾^(t)]]>

其中， 为 的估计值， 为 的估计值，L为精细抗干扰滤波器的增益矩阵；

定义针对无拖曳卫星相对位移通道的精细抗干扰滤波器的观测误差为 则其观测误差的动态方程为：

e·x‾(t)=(A‾-LC‾)ex‾(t)+F‾[f(x‾(t),t)-f(x‾^(t),t)]+B‾[w(t)+d1(t)]-Lv(t)+G‾δ(t)z1(t)=T1ex‾(t)z2(t)=T2ex‾(t)]]>

其中，z₁(t)和z₂(t)为针对无拖曳卫星相对位移通道的精细抗干扰滤波器的观测误差的参考输出，T₁和T₂为相应的加权矩阵， 为上文中增广的无拖曳系统模型的系数矩阵，U为满足Lipschitz条件的已知矩阵，L为精细抗干扰滤波器的增益矩阵， 为非线性项， 为非线性估计值，d₁(t)和δ(t)为能量有界干扰，w(t)和v(t)为高斯白噪声。针对无拖曳卫星的精细抗干扰滤波方法的思想为：通过设计一个观测器增益L，在利用干扰观测器补偿大气阻力干扰的同时，使得从能量有界干扰d₁(t)和δ(t)到z₁(t)的闭环传递函数的H_∞范数不超过给定的上界γ₁和γ₂，而从高斯白噪声w(t)和v(t)到z₂(t)的闭环传递函数的H₂范数γ₃尽可能的小。

(5)基于凸优化算法求解精细抗干扰滤波器的增益矩阵L，针对无拖曳卫星相对位移通道的精细抗干扰滤波方法可以表述为下列凸优化问题：

给定γ₁>0，γ₂>0，如果存在正定对称阵P＝P^T>0

一种无拖曳卫星相对位移通道的精细抗干扰滤波方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。