空间惯性传感器数据驱动自适应非对称约束控制.pdf

空间惯性传感器数据驱动自适应非对称约束控制孙笑云1,2，吴树范1,2，沈 强1,2（1.上海交通大学 航空航天学院，上海 200240；2.上海市引力波探测前沿科学研究基地，上海 200240）摘 要：在空间引力波探测航天器平台系统控制框架下，针对空间惯性传感器这一关键载荷具备非线性未建模动态且存在性能约束的高精度控制问题，提出一种基于数据驱动理论的数据驱动自适应控制方案，实现空间惯性传感器动力学系统作为一类非仿射非全局Lipstchiz连续系统时的精确稳定控制目标。基于模糊规则建立系统附加不确定性估计器，利用其通用逼近特性保证估计误差的有界性。基于控制障碍函数（Control Barrier Function，BLF）构建非对称性能约束，利用BLF设计的控制器实现闭环信号的非对称约束控制。根据收缩映射原理及离散时间系统Lyapunov理论分析各闭环信号及自适应估计的有界性，数值仿真验证了该数据驱动自适应非对称约束控制方案的可行性和有效性。关键词：控制障碍函数；自适应控制；数据驱动；空间引力波探测；非对称约束中图分类号：V448.2 文献标识码：A 文章编号：2096-9287（2023）03-0322-12DOI:10.15982/j.issn.2096-9287.2023.20220094引用格式：孙笑云，吴树范，沈强.空间惯性传感器数据驱动自适应非对称约束控制J.深空探测学报（中英文），2023，10（3）：322-333.Reference format:SUN X Y，WU S F，SHEN Q.Data driven-based asymmetric constrained control for spaceinertia sensorJ.Journal of Deep Space Exploration，2023，10（3）：322-333.引言1015m/s2/Hz1 nm/Hz空间惯性传感器是一类航天器系统关键载荷，广泛应用于空间探测领域。其工作原理为悬浮在航天器外壳内，为在轨航天器提供惯性基准1。低频段空间引力波探测是以航天器平台超静、超稳为前提开展的深空探测任务，利用空间惯性传感器提供高精度惯性基准，是探测工程实施的重要保障2-3。近年来，随着低频段空间引力波探测任务的发展及探测水平的提高，航天器平台中用于探测引力波信号的敏感轴在mHz频段的残余扰动加速度需优于，同时对于航天器平台的位移控制精度需优于4。针对低频段空间引力波探测这一特殊背景，研究空间惯性传感器在复杂干扰下的控制问题也是任务实施的一项重要研究内容。当前，针对空间惯性传感器控制方法的独立研究还较少，多数研究内容将该载荷的航天器系统平台作为一个研究整体，并认为空间惯性传感器的控制方法设计应包含在航天器系统平台的无拖曳控制系统（Drag-Free and Attitude Control System，DFACS）中，即将空间惯性传感器中仅用于内部稳定的自由度H作为非敏感轴，在控制方案的设计中予以考虑由航天器外部扰动及机械接口形变、传感器刚度形变及敏感轴交叉等作用等引起的耦合现象5。在既往研究中，静电悬浮控制技术的发展往往依托于性能指标而展开，在无拖曳控制系统的设计中进行统一讨论，如文献68提出的基于定量反馈理论及混合灵敏度的控制器设计方法，将性能指标直接转化成设计准则，用于满足不同的性能规范并实现科学探测模式下全频段能量消耗的降低。然而，随着高精度空间惯性传感器研制进度及建模理论的进一步推进，实际惯性传感器系统中的强烈非线性及未建模动态难以忽略，基于线性方法的上述控制方案均具有一定的局限性。尽管现有线性控制方法能够实现空间引力波探测任务中空间惯性传感器非敏感轴的指标要求，但随着空间惯性传感器系统允许的残余扰动性能指标要求逐渐提高，现有线性控制方法均以默认解耦的单输入单输出系统为设计基础，实际系统中存在的刚度耦合及电压驱动耦合关系在未得到解耦或抑制前惯性传感器系统难以满足新的性能约束，因此寻求一种可实现对复杂非线性、未知动态得以实现良好估计及抑制作用的非线性控制方案是当前 收稿日期：2022-10-25 修回日期：2023-05-16基金项目：科技部国家重点研发计划（2020YFC2200800，2021YFC2202600，2022YFC2204800）；青年科学基金项目（62103275）；上海市自然科学基金面上项目（20ZR1427000）第 10 卷 第 3 期深 空 探 测 学 报（中英文）Vol.10 No.32023 年 6 月Journal of Deep Space ExplorationJune 2023空间惯性传感器闭环控制器设计的前沿研究方向。在惯性传感器这一特殊的高精度载荷控制问题研究中，硬件内部的电回路数据链传输实际为固定采样频率的离散时间系统形式，因此在连续时间系统内实现的控制方法，如文献910中提出的自适应控制方案，即使其考虑了非线性动态的影响，也仍然难以在实际的惯性传感器系统中实现转化和运用。考虑最坏的情况，即由于电压驱动耦合作用及位移传感、交流执行等诸多现象的叠加，惯性传感器的动力学模型是非仿射且未知的，此时一类基于离散采样数据的数据驱动控制方案对该系统的有效性值得探讨和研究。数据驱动控制方案是考虑实际采样数字系统（如硬件在环）仅能获取离散量测输出信息时，利用量测数据的迭代实现闭环系统反馈控制的应用方案。目前，数据驱动控制方案已在随机衰减系统、卫星姿态动力学系统中获得相应的应用，因控制器设计仅需相邻已知的离散量测数据进行一步迭代，不涉及积分、微分器的引入，针对实际硬件系统拥有较好的可行性及简易性。文献1113提出一种基于自抗扰控制的数据驱动控制方案，运用扩张状态观测器实现外连续扰动的抑制，该方案针对惯性传感器控制系统的设计具有一定的启发作用，但仍需扩展至任意连续甚至非连续干扰存在的情况。性能指标是衡量空间惯性传感器作为惯性基准的重要尺度，在低频段空间引力波探测任务中，针对位移响应的性能指标要求也提升至较为精密的程度2,14-15。然而，在全过程控制中，仍然希望空间惯性传感器的静电悬浮控制系统能对被控对象的位移响应起到主动约束的作用，且该约束可以是非对称的、时变的。控制障碍函数（Control Barrier Function，BLF）是解决这一类约束控制需求的有效方法，文献1617提出的各类障碍函数控制方案解决了一类具有延迟、忽略初值的非对称时变约束控制问题，而上述方案在空间惯性传感器控制问题中也尚未得到考察和应用。为实现惯性传感器静电悬浮控制系统中仅能获得离散量测数据时的高精度控制问题，考虑性能约束要求及由静电驱动引起的非线性任意连续扰动抑制，利用数据驱动自适应控制方法，本文提出一种基于模糊规则及BLF的空间惯性传感器数据驱动约束控制方案，保证惯性传感器系统静电悬浮控制全过程满足性能约束。文章创新点如下：1）提出一种基于模糊规则的数据驱动自适应模型估计方法，用于解决空间惯性传感器任意连续未建模动力学的估计与抑制，相比现有成果将数据驱动估计方法自Lipstchiz连续目标估计推导到非Lipstchiz连续的情况；2）考虑空间惯性传感器实际任务需求下的性能约束问题，基于BLF设计一种非对称约束控制方案，基于BLF的输入补偿使得闭环系统控制全过程位移响应位于性能约束内；3）考虑惯性传感器实际控制对象的离散时间系统形式，基于数据驱动方案设计的离散时间系统自适应控制方案是在空间惯性传感器控制问题中的首次应用。文章章节安排如下：第1节对空间惯性传感器静电悬浮控制问题进行描述，给出静电力驱动关系和基于三体动力学的建模过程，并基于I/O量测数据对连续时间系统模型进行本地离散线性化；第2节给出控制系统设计结果，描述基于性能指标建立的非对称时变性能约束及BLF建立过程，设计基于模糊规则的自适应未知动态估计器及数据驱动自适应控制器；第3节数值仿真，对比验证该控制方案效果与控制系统性能，第4节给出结论。1 问题描述本节将对双质量块空间惯性传感器进行动力学建模，给出相应模型假设并实现本地动力学线性化。1.1 静电力驱动关系本文的任务背景参考LISA Pathfinder空间引力波探测试验验证任务14。在空间引力波探测航天器系统进行无拖曳控制的过程中，文献1820指出空间惯性传感器的电压驱动模式随DFACS的工作模式发生变化，本文进行的研究则是在DFACS执行科学模式时，空间惯性传感器随之执行的高分辨率模式（High Resolu-tion Mode，HR）下进行，此时的惯性传感器动力学系统已随达到临界稳定18，对HR模式下惯性传感器的控制系统设计则旨在提高各非敏感轴位移精度。x,作为空间惯性传感器动力学模型建立的初步，首先需描述HR模式下传感器各极板广义静电力及由控制输入决定的理想驱动力的转换关系。根据文献1819，在给出的极板广义静电力解算描述下，忽略电场边缘效应和驱动电子器件输出阻抗，给出由对置极板电压决定的一对平动及转动自由度上的实际静电力表达形式Fx=120Axd2x(2V21x2V22x+Ucouple1)+0Axd3x(2V21x+2V22x+2V21+2V22+Ucouple2)F=120Axd2x(2V212V22+Ucouple3)+0Axd3x(2V21x+2V22x+2V21+2V22+Ucouple4)(1)第 3 期孙笑云，等：空间惯性传感器数据驱动自适应非对称约束控制323FxF0AxdxV1,2x,UcoupleV1,2x,其中：、为各自由度对应的极板静电力；、分别为真空介电常数、极板面积及极板间距；为对置的各极板在平动及转动自由度的电势；为平动与转动自由度间电势贡献的耦合项，可利用表达为Ucouple1=2(V1x+V2x)(V1V2)Ucouple2=Ucouple4=2(V1xV2x)(V1V2)Ucouple3=2(V1xV2x)(V1+V2)V1,2x,x,Kxx=Kxx,req,K=K,reqfx、fV1x、V2x、V1、V2其次需描述惯性传感器内部电压驱动关系，即理想输入力向各极板电压的解算。仍以方向为例，在高分辨率模式下，规定极板电压驱动刚度为特定的需求刚度，即定义为。在传感器接受输入力时，给出各极板上平动及转动自由度反馈电压的表达形式V1x=12dx0Ax2fx+dx?Kxx,req?V2x=12dx0Ax2fx+dx?Kxx,req?V1=12Rx()dx0Axf+dxRx()2R22x()?K,req?V2=12Rx()dx0Axf+dxRx()2R22x()?K,req?(2)Rx()R22x()x其中：、是各极板电容与线位移、角位移的一、二阶偏导数的对应系数。上述反馈电压加载在各极板后，其产生的静电力将作用于惯性传感器的检验质量上，用于补偿检验质量因附加干扰产生的运动学偏差，从而实现检验质量的稳定惯性基准作用20。注意到理想驱动力与实际静电力之间，由于电压驱动算法的引入将具有如下的转换关系Fx=120Axd2x2dxfx0Ax+Ucouple1(V,Vx)+0Axd3x2d2x?Kxx,req?0Ax+Ucouple2(V,Vx)xF=120Axd2x2dxfx0Rx()Ax+Ucouple3(V,Vx)+0Axd3x2d2x?K,req?0R22x()Ax+Ucouple4(V,Vx)(3)fx,这也构成了传感器动力学模型建立过程中主要的非线性来源之一。在动力学模型建立中，由于控制输入与理想驱动力等价，上述过程也将直接引入动力学模型输入部分的建立中。1.2 空间惯性传感器动力学模型在研究空间惯性传感器动力学模型建立前，由于需考虑双质量块与航天器外壳的运动耦合作用，将首先以LISA Pathfinder空间引力波探测任务为样例1,21-22建立以三体动力学为基础的无拖曳控制系统动力学模型。执行探测任务的航天器构造如图1所示，由两个对置的惯性质量块TM1、TM2及航天器自身构成3个相互关联的物理实体。根据文献21中的分析，该多体动力学方程可以描述为下述近似二阶形式Mq q=KLTPq+Khq+fstray+factuation+Mu usc(4)factuationfstrayKh、KLTPqqMq、Muusc其中：为低频悬浮执行器作用力；为其它外部或内部来源构成的总噪声力；为惯性传感器耦合刚度；是检验质量块相对航天器的真实位移，需考虑每个检验质量外壳、航天器接口等构建的机械形变；则为相对位移的标称值；均为敏感度矩阵，用于衡量检验质量内部及检验质量对航天器的绝对运动；为航天器的绝对运动。TM1TM2r1xzIyIxIryr2z探测航天器rO1rO2 图 1 LISA Pathfinder航天器多体动力学构型18Fig.1 Diagram of multi-bodies dynamics for LISA Pathfinder spacecraft18 qusc将质量块相对位移 及航天器绝对运动描述为q=rT1,T1,rT2,T2T,usc=rT,TTMq、Mu同时给出敏感度矩阵的表达形式Mq=m1E33033033033033I1033033033033m2E33033033033033I2Mu=m1E33m1 ro1033I1m1E33m2 ro1033I2324深空探测学报（中英文）2023年m1m2I1I2 ro1、ro2ro1=ro1,x,ro1,y,ro1,zTro2=ro2,x,ro2,y,ro2,zT其中：、为惯性质量块质量及转动惯量；为斜对称交叉矩阵，由标称位置矢量及给出 ro1=0ro1,zro1,yro1,z0ro1,xro1,yro1,x0 ro2=0ro2,zro2,yro2,z0ro2,xro2,yro2,x0q=q忽略各接口及构件机械形变，认为，同时假设系统动力学仅存在小角度旋转，给出系统动力学的简化形式 r1 1 r2 2=033E33033033033033T1BT1B ro1E33033033033033T1B033E33033033T2BT2B ro2033033E33033033T2B033033033E33aa11a22T1BT2B其中：、为在标称位置下，航天器到检验质量块的位置变换矩阵。定义航天器和检验质量加速度=I1t,1=I11t1,2=I12t2a=1mf,a1=1m1f1,a2=1m2f2m、m1、m2、I、I1、I2t、t1、t2、f、f1、f2其中：分别为航天器及检验质量块的质量及转动惯量；为合外力矩和合外力。DDF、DSUS本文针对空间惯性传感器静电悬浮控制问题的研究建立在科学模式下，在LISA Pathfinder任务规则中，将这种模式称为科学探测模式1或测试模式M323。敏感轴及非敏感轴选取的原则为选取检验质量1的3个平动自由度、1个转动自由度及检验质量2的两个平动自由度施行无拖曳控制，其余6个自由度施行静电悬浮控制10。引入实现上述选取规则的坐标选择矩阵，将三体动力学系统重新表述为 qDF qSUS=BATT063DDFB1DDFB2DSUSB1DSUSB2(aSCaTM)(5)qDF、qSUSB1、B2、BATT其中：分别为无拖曳控制与静电悬浮控制坐标。为参数矩阵，定义为BATT=033E33，B1=T1BT1B ro10T1BT2BT2B ro20T2B，B2=E33,aSC、aTMaSC=(aT,T)T，aTM=(aT1,T1,aT2,T2)TuT、uS则用于描述航天器平台与检验质量所受合外力及力矩，。认为合外力及力矩由控制器输入，外界干扰dSC、dTM，及检验质量刚度变形构成，将式（5）中的动力学模型整理为 qDF qSUS=BATT063063BDFE66066BSUS066E66(uTuS1uS2)+dSCdTM1dTM2+0630630630662DF0660662C2SUS(qDFqSUS)(6)BDF=DDFB1，BSUS=DSUSB1uS1=DDFB2uS，uS2=DSUSB2uSdTM1、dTM22DF、2SUS2C其中：；悬浮控制输入；为系统输入噪声；对角矩阵与交叉耦合矩阵为刚度矩阵。仅针对静电悬浮控制通道进行分析，对式（6）中描述的运动学模型进行简化，表达为 qSUS=DSUSB1DSUSB2aSCaTMT(7)qSUS=1,1,x2,1,1,2TaSC、aTM其中：。根据式（6）将表达为aSC=uT+dSC,aTM=uS+dT+2DF002SUSqdSC、dTuT、uSx=qSUS,qSUST其中：为外界干扰；分别为微推力器控制力及悬浮控制自由度极板静电力。定义状态量，则计及干扰及静电力输入驱动关系式（3）的静电悬浮控制回路标准形式可表述为 x=Ax+B(f(uS,qSUS)+dTM)(8)dTMA、BA=0I662SUS0,B=0I66A其中：为系统总噪声；为系统状态参数矩阵，。在该简化的动力学模型中，状态矩阵中的交叉耦合刚度矩阵存在有界的慢时变余项，电压驱动关系及外源性扰动也将对原有的线性状态空间模型引入较为强烈的非线性动态。1.3 非仿射离散时间动力学模型及局部线性化在实际的空间惯性传感器静电悬浮控制实现中，精确的连续时间动力学模型由于电回路信号采集等因素不能准确获得，实际的控制系统设计过程也将基于离散的I/O量测输出数据。在空间关键载荷乃至航天器姿态系统的控制问题中，根据文献17，为便于控制系统的分析设计，也可利用量测数据将动力学模型式（8）转换为如下的非仿射离散时间系统表达形式x(k+1)=f(x(k),x(knx),uS(k),uS(knu)+h(d(k1),d(knd)(9)h(d(k1),d(knd)其中：为非线性系统干扰；第 3 期孙笑云，等：空间惯性传感器数据驱动自适应非对称约束控制325f(x(k),x(knx),uS(k),uS(knu)则为包含不确定性的非线性动力学标称模型。为便于问题研究，需对模型式（8）进行如下假设。假设假设1 输入的变化率存在如下约束?uSi(k)uSi(k1)?uuS(k)=uS1(k),uS2(k),;uSn(k)T，ui 0ii=1,2,n其中：为第个输入变化率绝对值的上界，。x(k+1)=f(x(k),x(knx)k1、k2假设假设2 在该未知的非线性系统中，标称状态方程满足全局Lipschitz条件，即对任意两个采样点，都有f(x1(k1),xN(k1)f(x1(k2),xN(k2)L1x1(k1)x1(k2)+LNxN(k1)xN(k2)(10)Lmm 1,NN=nx+nd+nu+3f()其中：为假设存在的Lipschitz常数，。同时，对控制输入的偏导数存在，不为零，符号不变。h(d(k1),d(knd)h()h1()h2()假设假设3 有界的系统干扰非Lipschitz连续，即任意连续。在实际问题中，存在因电压驱动耦合引起的与状态相关的非Lipschitz连续扰动及自相关的有界外界扰动，即h(d(k1),d(knd)=h1(x(k),x(knx)+h2(d0(k1),d0(knd0)(11)?h2(d0(k1),d0(knd0)?Eh同时应满足。根据上述假设，将非仿射离散时间空间惯性传感器动力学模型式（9）改写为下述差分形式x(k+1)=f(k),uS(k)f(k1),uS(k1)f(k),uS(k1)+f(k),uS(k1)+h(k)h(k1)=f()uS(k)uS(k)+(k)(12)(k)=x(k),x(knx),uS(k),uS(knu)(k)=d(k1),d(knd)(k)=f(k1),uS(k1)+f(k),uS(k1)+h(k)h(k1)x(k+1)x(k+1)x(k)T(k)=f()uS(k)T(k)其中：；为差分算子，满足。定义，认为有界，且存在上界，则系统转换为如下的动力学局部线性形式x(k+1)=T(k)uS(k)+(k)(13)T(k)(k)f()为实现上述模型的稳定控制，控制器的设计需同时实现未知增益的估计和任意连续未知动态的逼近，从而使各闭环信号有界。考虑到标称动力学模型分量的全局Lipstchiz连续性，标称控制输(k)入基于引入数据驱动理论的自适应控制方案，而任意连续未知动态则基于模糊规则进行逼近。同时，控制过程需根据性能指标要求设计非对称约束，运用BLF确保控制全过程位移响应位于约束内。2 基于数据驱动的模糊自适应约束控制方案设计(k)T(k)本节将首先给出任意连续未知动态的模糊估计器，其次设计的估计自适应律，给出闭环控制器，验证各闭环信号的有界性。2.1 基于模糊规则的自适应未知动态估计器(k)1(k)2(k)注意到连续未知动态由任意连续扰动分量及有界外界扰动分量构成(k)=1(k)+2(k)(14)2(k)2Eh1(k)1(k)根据式（12），。对于，由于其由Lipstchiz连续的标称模型差分项及代表耦合作用的有界任意连续扰动项构成，仍然有界。考虑模糊系统的通用逼近特性，设计如下IF-THEN模糊规则：If x1 Fl1,x2 Fl2,xn Fln,then y Gl.x0=x1,xnT Rn、y Rl=1,2,NNFln、GlhFli(xi)、hGl(y)Fln、Gl其中：分别为模糊系统输入及输出；，为模糊规则数；为模糊集。令分别为模糊集的对应隶属度函数，那么模糊系统输出可以表达为y(x)=Nl=1lni=1hFli(xi)Nl=1ni=1hFli(xi)(15)l=maxyRhGl(y)其中：中心值，由此定义模糊基函数l(x)=ni=1hFli(xi)Nl=1ni=1hFli(xi),l=1,2,N(16)T=1,N=1,N(x)=1,NT1(k)令，将任意连续扰动分量基于上述模糊规则改写为如下形式1(k)=T(k)(k)(17)0 0根据模糊基函数的定义，容易得到。假设输入变量处于内的紧集中，对于任意，存在模糊规则式（17）及最优权值，使得326深空探测学报（中英文）2023年sup(k)1(k)*T(k)(18)由此将本地动力学线性化模型式（13）改写为x(k+1)=T(k)uS(k)+T(k)(k)(19)T(k)、T(k)T(k)、T(k)T(k)其中：均未知。运用自适应律对进行估计。对于，首先给出如下目标函数J1(T(k)=?x(k)T(k)uS(k1)?2+1?T(k)T(k1)?(20)1 0?T(k)T(k1)?T(k)T(k)J1(T(k)其中：代表的权重系数；为的自适应估计。最小化并推导得出T(k)=T(k1)+1(x(k)T(k1)uS(k1)uTS(k1)1+uS(k1)2(21)1 0T(k)其中：为自适应增益。类似的，为估计得到，给出如下目标函数J2(T(k)=?x(k)T(k)(k1)?2+2?T(k)T(k1)?(22)2 0?T(k)T(k1)?T(K)T(k)J2(T(k)其中：代表的权重系数，为的自适应估计。最小化并推导得出T(k)=T(k1)+2(x(k)T(k1)(k1)T(k1)2+?(k1)?2(23)2 0其中：为自适应增益。为分析自适应估计式（20）与式（22）的有界性，首先给出如下引理：(k)引理引理1（收缩映射原理）（收缩映射原理）考虑如下迭代状态(k+1)=(k)(k)+(k)(24)(k)RMM(k)RM(k)1,limk(k)=0limk(k)=0其中：为映射矩阵：为外界输入。如果，那么。根据上述引理给出自适应估计式（21）与式（23）的稳定性分析过程。给出如下定理。0 1，2 2定理定理1 对于式（13）表达的非仿射非线性离散时间系统模型，当未知参数的自适应估计按式（21）及式（23）表达时，给定自适应增益，各自适应估计误差有界，并给定误差上界limk?Ti(k)?=b11a1,limk?T(k)?=b21a2其中a1=?I1uS(k1)uTS(k1)1+uS(k1)2?b1=1(?i?+2Eh)21+2a2=?I2(k1)T(k1)2+?(k1)?2?b2=2(u+2Eh)22Ti(k)T(k)iii代表第 行的估计误差；为最优权值的第 行向量。Ti(k)证明证明 首先证明的有界性。回顾自适应律式（21），将估计误差推导为如下表达形式Ti(k)=Ti(k1)+Ti(k)11+uS(k1)2(i(k1)Ti(k1)uS(k1)Ti(k1)uS(k1)uTS(k1)=Ti(k1)(I1uS(k1)uTS(k1)1+uS(k1)2)1i(k1)uTS(k1)1+uS(k1)2+Ti(k)(25)i(k1)ik1Ti(k)=Ti(k)Ti(k1)Ti(k)其中：为 第 行的第次迭代，。根据的有界性，?1i(k1)uTS(k1)1+uS(k1)2+Ti(k)?1(1i(k)+2i(k)21?+21(?i?+2Eh)21+2(26)1i,2i1,2i其中：代表的第 行。注意到?Ti(k1)(I1uS(k1)uTS(k1)1+uS(k1)2)?Ti(k1)?I1uS(k1)uTS(k1)1+uS(k1)2?(27)令a1=?I1uS(k1)uTS(k1)1+uS(k1)2?b1=1(?i?+2Eh)21+20 1 00 a1 1Ti(k)limk?Ti(k)?=b11a1考虑到，满足，从而根据引理1可以得到的有界性，其收敛上界可表达为。第 3 期孙笑云，等：空间惯性传感器数据驱动自适应非对称约束控制327T(k)类似的，也可得到模糊权值估计误差的有界性。根据自适应律式（23）得出T(k)=T(k1)I2(k1)T(k1)2+?(k1)?22(T(k1)uS(k1)+2(k)T(k1)2+?(k1)?2(28)(k)(k)=(k1)式（28）用到了的不变性，即。类比式（26），对式（28）中部分项进一步推导得出?2(T(k1)uS(k1)+2(k)T(k1)2+?(k1)?2?2(u+2Eh)22=b2(29)令a2=?I2(k1)T(k1)2+?(k1)?2?0 2 00 a2 0其中：。将式（33）代入式（32）并推导得出V1(k)M1xT(k)(T(k)uS(k)+1(k)+121M1xT(k)x(k)x(k)+2M1xT(k)x(k)x(k)+1(34)1=1/41+E2h/42其中：。给出下述控制器v(k)=T(k)(k)+x(k)H(k)+?T(k)?uSi(k)=uSi(k1)+usgn(vi(k),|vi(k)|uuSi(k1)+vi(k),|vi(k)|u(35)H(k)=121M1xT(k)x(k)x(k)+2M1xT(k)x(k)x(k)uSi(k)、vi(k)i 0?uSi(k)?|v(k)|uS(k)=G(k)v(k)0 G(k)1其中：；为第 个输入及输入变化率的取值；为反馈增益，满足。注意到，令，则。将控制器（35）代入（34）得到V1(k)M11xT(k)x(k)+M1(I1)xT(k)(H(k)+1(k)+M11xT(k)(1(k)T(k)(k)+1M1(1+12)xT(k)x(k)+M112(?b21a2(k)?+)2+M1I12(H(k)+?(k)?+)2+1M1(1+12)xT(k)x(k)+M1(b21a2)2+4M1(?H(k)?+)2+1(36)1=T(k)G(k)+T(k)1 1I1 20?(k)?1?H(k)?(x,xd)(x,xd)其中：，显然，从而。式（36）的推导还利用了定理1的结论及。考虑，注意到为一标量连续328深空探测学报（中英文）2023年:=(x,xd):V 0(x,xd)函数，考虑存在一紧集，在紧集上，存在的最大值，从而式（36）可转化为V1(k)M1(1+12)xT(k)x(k)+M1(b21a2)2+4M1(+?*?)2+1(37)1/21V pV1(k)在时，在上，满足V1(k)0p+0(38)其中：0=M1(1+12),0=M1(b21a2)2+4M1(+?*?)2+1.对式（38）两边积分，得到0 V(k)00+(V(0)00)e0k(39)V(k)LV p由此得到。由于为一不变集，系统各闭环信号半全局一致有界。由于已得到BLF的有界e(k)e(k)(H1(k),H2(k)H1(k),H2(k)性，表明系统跟踪误差满足，而代表可任意选取的时变非对称性能边界。根据（39）可进一步得到eT(k)e(k)2(H1(k)+e(k)T(H2(k)e(k)00+(V(0)00)e0k(40)e(k)(H1(k),H2(k)根据Youngs不等式，在紧集上e(k)12(H1(k)+H2(k)00+(V(0)00)e0k(41)k 随着采样步数，可推导得出limke(k)12(H1(k)+H2(k)00(42)0式（42）进一步表明了跟踪误差信号的一致有界性，且随着的合理选取，跟踪误差信号的稳态误差也将减小。根据本节所述闭环控制器的设计结果与描述，给出该数据驱动约束控制方案结构框图，如图2所示。性能约束e(k)(k)障碍函数H1(k),H2(k)模糊系统故障/不确定性估计器交流反馈控制器自适应反馈增益估计器(k)(k)T(k)(k)(k)us(k)x(k)零阶保持器惯性传感器动力学 图 2 控制结构Fig.2 Control framework 3 仿真结果与分析 3.1 仿真参数设置数值仿真验证基于样例空间引力波探测任务给出的空间惯性传感器非线性动力学模型7,20。由于需考察闭环系统各状态响应的性能指标，部分数值仿真结果展示为频域形式。数值仿真通过与QFT控制方案7的对比验证了该数据驱动自适应约束控制方案在2个分别代表平动和转动的非敏感轴静电悬浮自由度上针对位移噪声抑制的控制效果，设置仿真时间为3 000 s，步长为0.1 s。数值仿真结果如图35所示。在空间惯性传感器动力学系统中，存在以如下成型滤波器表达的随机过程噪声第 3 期孙笑云，等：空间惯性传感器数据驱动自适应非对称约束控制329L=K1s+z1s+p1+K2(s+p2)(s+p3)(s+p4)(s+p5)(43)x2,1p1，p2，p3，p4，p5，z1，K1，K2其中：作用于两敏感坐标轴上的噪声相关参数的取值由文献23给出。在动力学模型有关的参数中，选取带有摄动的刚度矩阵2DF=2SUS=107(1+5%)11.191.351.350.004 25001.359.551.350.004 25001.351.3524.120.004 250026.08726.08726.08730.640000009.551.3500001.3524.12.N=3在控制器参数的选取中，选取模糊集数量，并将模糊基函数按下式表述i(x)=HFli(x)/Nl=1HFli(x)(44)其中hFli(xi(k)=exp(xi+71012l32)2x2，1=0.011=2=11051=2=0.0011=10141=1010 x2=101=5.5在两通道上的其它控制器参数整定为，。两通道预设的非对称时变性能边界函数选取为H1x2(k)=(0.22exp(0.01k)+0.02cos(0.003k)+1.9)1012H2x2(k)=(0.52exp(0.003k)+0.02cos(0.003k)0.4)1012H11(k)=(5exp(0.001k)+0.8cos(0.001k)+1)1014H21(k)=(0.2exp(0.003k)+0.07cos(0.000 5k)+1.45)1013 3.2 仿真结果根据数值仿真结果，展开如下分析：1 nm/Hzxd=01）图35展示了该数据驱动模糊自适应约束控制方案的性能及与常规线性控制方法的控制效果对比。其中，图3展示了该方案下两自由度位移响应的频域幅值谱密度响应，在预设的干扰及动力学参数下，该方案的位移响应优于样例空间引力波探测任务给出的性能指标。图4则展示了两自由度内该方案与QFT方案的时域控制效果对比，仿真结果表明，本文提出的方案对于扰动及非线性不确定性具有明显的抑制作用，但针对的跟踪目标，两种方案均具有一定的稳态误差，这一结果也与理论分析相对应。图5表明了两通道的位移响应位于时变且非对阵的性能约束内，满足预设的性能指标要求。104103102101100101探测频段/Hz10201018101610141012幅值谱密度/（m（rad）Hz1/2）位移噪声x2 轴z2 轴 图 3 非敏感轴频域位移响应Fig.3 Response of frequency domain on non-sensitive axis 05001 00015002 0002 5003 0003210QFT 方案自适应方案05001 0001 5002 0002 5003 00010505QFT 方案自适应方案时间/s时间/s1 轴位移1013 radx2 轴位移1012 m 图 4 与QFT方案的位移响应效果对比Fig.4 Displacement response compared with QFT-based scheme 05001 0001 5002 0002 5003 0002.01.51.00.50 x2轴响应性能约束05001 0001 5002 0002 5003 0005051015201轴响应性能约束时间/s时间/sx2 轴位移1012 m1 轴位移1014 rad330深空探测学报（中英文）2023年05001 0001 5002 0002 5003 0002.01.51.00.50 x2轴响应性能约束05001 0001 5002 0002 5003 0005051015201轴响应性能约束时间/s时间/sx2 轴位移1012 m1 轴位移1014 rad 图 5 性能约束的验证Fig.5 Verification of performance constraints 2）图6给出了基于自适应方案得到的模糊权值及Lipschitz常数的估计结果。在两个非敏感轴通道，仿真曲线可表明各项自适应估计的有界性，即随着控制过程的递进，各估计有自适应趋稳的趋势。结合与线性方案对比的数值仿真结果，可以得出，在该静电悬浮控制问题中，系统中存在较为明显的未建模Lipschitz连续及任意连续分量，也表明针对非线性项及未建模动态的估计及反馈对控制效果的整体提升作用是显著的。05001 0001 5002 0002 5003 00010505x2 轴1 轴05001 0001 5002 0002 5003 0002101x2 轴1 轴时间/s时间/s/1016/1014 图 6 各自适应变量估计结果Fig.6 Result of each adaptive estimation 4 结论在考虑硬件数据传输、实际系统运行在离散时间系统内的前提下，本文提出了一种数据驱动离散时间模糊自适应约束控制方案，用于解决空间引力波探测任务中惯性传感器静电悬浮高精度控制问题，运用数据驱动方案及模糊系统实现系统内任意连续未知分量的估计与反馈，同时基于控制障碍函数未系统位移响应设计非对称时变性能约束，使闭环控制响应满足预设的控制性能。理论分析和数值仿真结果均验证了各闭环响应的有界性，数值仿真表明了该方案对预设性能指标的满足。参考文献FICHTER W，GATH P，VITALE S，et al.LISA Pathfinder drag-freecontrol and system implicationsJ.Classical and Quantum Gravity，2005，22（10）：S139.1MOBLEY F，FOUNTAIN G，SADILEK A，et al.Electromagneticsuspension for the tip-II satelliteJ.IEEE Transactions on Magnetics，1975，11（6）：1712-1716.2刘伟，高扬.空间引力波探测中无拖曳控制方法研究J.中国科学：物理学 力学 天文学，2020，50（7）：112-122.LIU W，GAO Y，Drag-free control methods for space-basedgravitational-wave detectionJ.SCIENTIA SINICA Physica，Mechanica