一种无拖曳卫星相对位移通道的精细抗干扰滤波方法
【技术领域】
[0001] 本发明设及一种无拖曳卫星相对位移通道的精细抗干扰滤波方法，可用于无拖曳 卫星相对位移通道的抗干扰滤波。
【背景技术】
[0002] 无拖曳卫星是一种采用无拖曳技术的卫星，通过控制卫星本体和检测质量块达到 高宁静的运动特性，能够达到航天领域对超静平台的要求;根据无拖曳卫星的原理，可将其 数学模型分为卫星本体姿态通道、相对位移通道W及检查质量块姿态通道。其中，无拖曳卫 星相对位移通道设计的基本思想是使无拖曳卫星跟随内部检测质量块运动，将无拖曳卫星 相对位移通道受到的多源干扰尽可能抵消，从而使无拖曳卫星运行在接近纯引力作用下的 轨道上;而无拖曳卫星运行在近地轨道中，面临大气阻力、太阳光压、执行机构噪声、量测噪 声W及模型不确定性项等多源干扰，并且由于多源干扰的影响致使无拖曳卫星相对位移通 道的无法达到其高控制精度的要求。
[0003] 目前专口针对无拖曳卫星相对位移通道的滤波方法的研究较少，被采用过的滤波 方法有传统的卡尔曼滤波、Hoo滤波W及基于干扰观测器的滤波。但是传统的滤波方法大多 只是针对单一干扰类型，如卡尔曼滤波可W有效抑制高斯零均值不相关白噪声的干扰，Hoo 滤波可W解决噪声统计特性未知的非线性系统滤波问题，基于干扰观测器的滤波方法可W 有效估计并补偿已知干扰模型的外部干扰。而对于含有多源干扰的无拖曳卫星相对位移通 道，仅使用针对单一干扰的类型的滤波方法，使其滤波精度不高，达不到无拖曳卫星相对位 移通道的高精度要求。因此，目前的针对无拖曳卫星相对位移通道的滤波方法存在只能针 对单一干扰进行滤波，滤波精度较低，保守性较大的问题。
【发明内容】
[0004] 本发明的技术解决问题是:克服现有的针对无拖曳卫星相对位移通道的滤波方法 只能针对单一干扰进行滤波，滤波精度较低，保守性较大的问题，设计一种无拖曳卫星相对 位移通道的精细抗干扰滤波方法。
[0005] 本发明的技术解决方案为:针对含有大气阻力、太阳光压、执行机构噪声、量测噪 声W及模型不确定性项等多源干扰的无拖曳卫星相对位移通道。首先，针对无拖曳卫星相 对位移通道所受多源干扰的数学表征并将多源干扰分类并建模;其次，建立含有多源干扰 的无拖曳卫星相对位移通道的数学模型；再次，构造针对无拖曳卫星相对位移通道的精细 抗干扰滤波器;最后，基于凸优化算法求解精细抗干扰滤波器的增益矩阵。具体实施步骤如 下：
[0006] (1)无拖曳系统相对位移通道所受的多源干扰分类与建模：
[0007] 大气阻力是无拖曳卫星相对位移通道面临的主要干扰，可对大气阻力建模并将其 描述成状态空间模型的形式： 命）=A/z"(/) +抓）
[000引 ' \d,(t) = Nrit)
[0009] 其中，do(t)为大气阻力干扰，T(t)为大气阻力数学模型的状态变量，M和N为系数 矩阵，S(t)为模型不确定性项；
[0010] 太阳光压和模型不确定性项s(t)的L2范数均有界，将其描述为能量有界干扰，即 '>,(〇|2曲< m，其中山(t)表示太阳光压和模型不确定性项； V.0
[0011] 执行机构噪声w(t)及量测噪声v(t)均为高斯零均值不相关白噪声，其统计特性满 足：
[0012] E{w(k)w^(j) I=QkSk-J
[0013] E{v(k)yT( j)}=化 Sk-j
[0014] E{w(k)v^(j)}=0
[0015] 其中，Qk，Rk分别表示k时刻执行机构噪声和量测噪声的协方差矩阵；Sk-J是 Kronecker-S 函数，即如果k = j 则 Sk-j = 1，否则 Sk-j = 0;
[0016] (2)建立含有多源干扰的无拖曳卫星相对位移通道的数学模型： 「00171 二 + 巧.('叫')，"斗公["(f)+ + 00。 -[.哦)三。(巧
[001引其中，状态变量邱)=[却(0 .培的了，Xi(t)=;r(t)，馬的.=.唯），;r(t)为无拖曳卫星 的质屯、与检测质量块的质屯、的相对位移，系数矩阵^
，msc为无拖曳卫星质量，Hltm为检测质量块质量， ? 0为轨道角速度，Ktrans为无拖曳卫星与内部检测质量块之间禪合的水平弹性系数，Dtrans为无 拖曳卫星与内部检测质量块之间禪合的水平阻尼系数，U(t)为无拖曳卫星所受到的控制力，W (tWPWt)分别为执行机构噪声和量测噪声，/W如)二2蚁龄-诚]带)+ 11> (of-(崎)3}/脚 为已知非线性项，其中，《(t)为卫星的绝对角速度，CO X (t)为相应的叉乘矩阵;假设非线 性项满足LipscMtZ条件，即对任意的两个系统状态Xi(t)，X2(t)，存在已知的矩阵U使得下 列不等式成立：
[0019] If(xi(t) ,t)-f(X2(t),t)I I < I IU(xi(t)-X2(t))
[0020] (3)建立增广的含有多源干扰的无拖曳卫星相对位移通道的数学模型： t(/')二方:(〇 + 巧(巧/)、〇 + 孩[w(0 + H.(0 + 攻(0] + (7()-(0
[00別]^ ' 、巧0 二巧0 + ^0
[002^ 其中，状态变量郭) = 量测输出W) = W。，系数矩阵
[0023] (4)构造针对无拖曳卫星相对位移通道的精细抗干扰滤波器：
[0024] 椅。二風')+ 巧巧(/).!') + 哀+ 化巧-巧/)] |.訂,'):二巧(0
[002引其中，託)为碟)的估计值，KO为京掛的估计值，L为针对无拖曳卫星相对位移通 道的精细抗干扰滤波器的增益矩阵；
[0026] (5)基于凸优化算法求解精细抗干扰滤波器的增益矩阵^
[0032] 其中，min表示最小值，Tl和T2为相应的加权矩阵，P和Z为正定对称阵，A为任意正实 数，丫 1为能量有界干扰di(〇到精细抗干扰滤波器的观测误差的参考输出的闭环传递函数 的Hoo范数的上限，丫 2为能量有界干扰S(t)到精细抗干扰滤波器的观测误差的参考输出的 闭环传递函数的Hoo范数的上限，丫 3为高斯白噪声W(t)和V(t)到精细抗干扰滤波器的观测 误差的参考输出的闭环传递函数的出范数;求解矩阵不等式可得未知矩阵Pl,针对无拖曳卫 星的精细抗干扰滤波器的增益阵L = P斗L。
[0033] 本发明与现有技术相比的优点在于：
[0034] (1)本发明充分考虑了无拖曳卫星相对位移通道含有的大气阻力、太阳光压、执行 机构噪声、量测噪声W及模型不确定性项等多源干扰，并进行分析建模，建立了含有多源干 扰的无拖曳卫星相对位移通道的数学模型；克服了现有无拖曳卫星相对位移通道的数学模 型未充分考虑多源干扰的问题；
[0035] (2)本发明构造了针对无拖曳卫星相对位移通道的精细抗干扰滤波器，结合干扰 观测器的思想，通过凸优化算法求解精细抗干扰滤波器的增益，利用出范数和Hoo范数，保证 滤波器的鲁棒性的同时可W有效补偿和抑制多源干扰;克服了卡尔曼滤波、Hoo滤波W及基 于干扰观测器的滤波等方法只能针对单一干扰进行滤波，滤波精度较低，保守性较大的问 题。
【附图说明】
[0036] 图1为针对无拖曳卫星相对位移通道的精细抗干扰滤波方法的流程框图。
【具体实施方式】
[0037] 下面结合附图及实施例对本发明进一步详细说明。
[0038] 如图1所示，本发明具体实现步骤如下：
[0039] (1)将无拖曳系统相对位移通道所受的多源干扰分类并建模：
[0040] 由于无拖曳卫星运行在近地轨道中，大气阻力是无拖曳卫星相对位移通道面临的 主要干扰，通过正弦曲线模拟大气阻力，可建立大气阻力的数学模型： 。、-+ 年 sin(巧/ +疼）+ Ar [0041 ]或片）二 < 幻,' + 4, sin(份/ + A') + A,' 岛：+ A sin(似/ +'疼）+ A:
[0042] 其中，d日(t)为大气阻力，ax，ay，az，Ax，Ay，Az为未知的幅值，(6x，(6y，(6z为未知的相 位，COd为大气阻力干扰的频率，Ax, Ay, Az为模型误差，再将其描述成状态空间模型的形 式： >(〇 = +抓）
[0043] VoW =轉)
[0044] 其中，T(t)为大气阻力模型的状态变量，M和N为系数矩阵，S(t)为模型不确定性项；
，N=[1 1 1 0 0 0]，S(t)为模型不确定性项
[0045] 太阳光压和模型不确定性项S(t)的L2范数均有界，将其描述为能量有界干扰 |4(〇||2成<00巧中di(t)表示太阳光压和模型不确定性项； a
[0046] 执行机构噪声及量测噪声均为高斯零均值不相关白噪声，其统计特性满足：
[0047] E{w(k)w^(j)} =QkSk-J [004引 E{v(k)vT(j)}=化Sk-j
[0049] E{w(k)v^(j)}=0
[0050] 其中，Qk，Rk分别表示k时刻执行机构噪声和量测噪声的协方差矩阵；Sk-J是 Kronecker-S 函数，即如果k = j 则 Sk-j = 1，否则 Sk-j = 0;
[0051] (2)建立含有多源干扰的无拖曳卫星相对位移通道的数学模型：
[0052] 无拖曳卫星相对位移通道设计的基本思想是使无拖曳卫星跟随内部检测质量块 运动，从而使无拖曳卫星运行在接近纯引力作用下的轨道上，假设无拖曳卫星的质屯、与内 部检测质量块空腔的质屯、重合，建立无拖曳卫星相对位移通道的动力学方程如下：
[0054]其中，r(t)为无拖曳卫星的质屯、与内部检测质量块质屯、的相对位移，《h为内部检 测质量块空腔的绝对角速度，Hltm和Hlsc分别为内部检测质量块和无拖曳卫星的质量，抗tm内部 检测质量块所受的引力，抗sc、Fcsc和Fdsc分别表示无拖曳卫星所受的引