2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种鲁棒的基于对抗结构的生物特征ROI提取方法

刘凤 刘浩哲 张文天 陈嘉树 沈琳琳 王磊

岳程斐, 霍涛, 陈雪芹, 沈强, 曹喜滨. 航天器姿态受限的协同势函数族设计方法. 自动化学报, 2024, 50(1): 54−65 doi: 10.16383/j.aas.c230409
引用本文: 刘凤, 刘浩哲, 张文天, 陈嘉树, 沈琳琳, 王磊. 一种鲁棒的基于对抗结构的生物特征ROI提取方法. 自动化学报, 2023, 49(6): 1339−1353 doi: 10.16383/j.aas.c200156
Yue Cheng-Fei, Huo Tao, Chen Xue-Qin, Shen Qiang, Cao Xi-Bin. Synergistic potential functions for constrained attitude control of rigid spacecraft. Acta Automatica Sinica, 2024, 50(1): 54−65 doi: 10.16383/j.aas.c230409
Citation: Liu Feng, Liu Hao-Zhe, Zhang Wen-Tian, Chen Jia-Shu, Shen Lin-lin, Wang Lei. A robust ROI extraction method for biometrics using adversarial structure. Acta Automatica Sinica, 2023, 49(6): 1339−1353 doi: 10.16383/j.aas.c200156

一种鲁棒的基于对抗结构的生物特征ROI提取方法

doi: 10.16383/j.aas.c200156
基金项目: 国家自然科学基金(62076163, 91959108, 61672357), 深圳市基础研究基金(JCYJ20190808163401646, JCYJ20180305125822769), 腾讯“犀牛鸟”深圳大学青年教师科学研究基金资助
详细信息
    作者简介:

    刘凤:深圳大学计算机与软件学院副教授. 2014年获得香港理工大学计算机系的计算机科学博士学位. 主要研究方向为模式识别和图像处理以及相关技术在指纹领域中的应用. 本文通信作者.E-mail: feng.liu@szu.edu.cn

    刘浩哲:深圳大学硕士研究生. 主要研究方向为计算机视觉和模式识别. E-mail: liuhaozhe2019@email.szu.edu.cn

    张文天:深圳大学硕士研究生. 主要研究方向为模式识别和生物特征识别. E-mail: zhangwentianml@gmail.com

    陈嘉树:英国纽卡斯尔大学硕士研究生. 主要研究方向为计算机视觉与生物识别. E-mail: gaasyu.chan@gmail.com

    沈琳琳:获得英国诺丁汉大学博士学位. 现为深圳市“鹏城学者” 特聘教授、英国诺丁汉大学计算机学院荣誉教授. 主要研究方向为深度学习理论及其在人脸识别/分析以及医学图像分析上的应用. E-mail: llshen@szu.edu.cn

    王磊:中国科学院深圳先进技术研究院副研究员. 主要研究方向为图像变换, 计算机视觉, 视觉语义理解, 视频分析, 深度学习. E-mail: lei.wang1@siat.ac.cn

A Robust ROI Extraction Method for Biometrics Using Adversarial Structure

Funds: Supported by National Natural Science Foundation of China (62076163, 91959108, 61672357), Shenzhen Fundamental Research Fund (JCYJ20190808163401646, JCYJ20180305125822769), and Tencent “Rhinoceros Birds” Scientific Research Foundation for Young Teachers of Shenzhen University
More Information
    Author Bio:

    LIU Feng Associate professor at the College of Computer Science and Software Engineering, Shenzhen University. She received her Ph.D. degree in computer science from the Department of Computing at Hong Kong Polytechnic University in 2014. Her research interest covers pattern recognition and image processing, and their applications to fingerprints. Corresponding author of this paper

    LIU Hao-Zhe Master student at Shenzhen University. His research interest covers computer vision and pattern recognition

    ZHANG Wen-Tian Master student at Shenzhen University. His research interest covers pattern recognition and biometrics

    CHEN Jia-Shu Master student at Newcastle University, UK. His research interest covers computer vision and biometrics

    SHEN Lin-Lin Received his Ph.D. degree from the University of Nottingham, Nottingham, UK. He is currently a Pengcheng Scholar Distinguished Professor at the College of Computer Science and Software Engineering, Shenzhen University. He is also a honorary professor at the School of Computer Science, University of Nottingham, UK. His research interest covers deep learning, facial recognition, analysis/synthesis and medical image processing

    WANG Lei Associate professor at Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences. His research interest covers image transforms, machine learning, computer vision, visual semantic understanding, video analysis, and deep learning

  • 摘要: 感兴趣区域(Region of interest, ROI) 提取在生物特征识别中, 常用于减少后续处理的计算消耗, 提高识别模型的准确性, 是生物识别系统中预处理的关键步骤. 针对生物识别数据, 提出了一种鲁棒的ROI提取方法. 方法使用语义分割模型作为基础, 通过增加全局感知模块, 与分割模型形成对抗结构, 为模型提供先验知识, 补充全局视觉模式信息, 解决了语义分割模型的末端收敛困难问题, 提高了模型的鲁棒性和泛化能力. 在传统二维(2D)指纹、人脸、三维(3D)指纹和指纹汗孔数据集中验证了方法的有效性. 实验结果表明, 相比于现有方法, 所提出的ROI提取方法更具鲁棒性和泛化能力, 精度最高.
  • 姿态控制是一个非线性控制问题, 在航空航天、机器人等领域受到了广泛关注和研究[1-2]. 用于描述刚体姿态的旋转矩阵构成了三维特殊正交群(Three-dimensional special orthogonal group, $ SO(3) $). 作为一个无边界的紧流形, $ SO(3) $不同胚于任何欧氏空间, 这导致在$ SO(3) $上, 不存在能够使姿态全局渐近收敛到目标姿态的连续时不变反馈控制律[3]. 受此拓扑性质约束, 使用光滑的反馈控制律能实现的最好结果是基于单一势函数获得几乎全局的收敛性[4-5]. 而“几乎全局”意味着该势函数在$ SO(3) $上存在一个测度为零的集合. 该集合中的点称为临界点, 当航天器初始姿态是临界点时, 势函数无法控制姿态收敛到目标姿态.

    为解决上述单一势函数存在的临界点问题, 混合控制被引入姿态控制问题中[6-7]. 混合控制是一种具有滞回特性的切换控制系统, 通过设计一族临界点互异的势函数, 系统状态可在不同的势函数之间进行切换, 从而改变系统状态和临界点的分布情况, 实现姿态控制目标. 具体而言, 当系统运动到当前势函数上不期望的临界点时, 会切换到另一个具有更低值且不处于临界点的势函数上进行控制. 这一类势函数被称为协同势函数族. 当单位矩阵是公共临界点时, 称此类协同势函数族为共单位元的协同势函数族. 文献[8]首次通过“角度扰动”方法构建了协同势函数, 但并未给出协同性的证明. 文献[9] 进一步给出了协同性对势函数参数的要求. 文献[10]给出了“角度扰动”方法中, 协同势函数族存在的充要条件, 并将协同势函数族应用到无角速度测量的航天器姿态控制系统中. 以上文献均需要对临界点进行复杂计算. 文献[11]放松了“共单位元”要求, 通过选择固定的参考向量, 设计出一族简单协同势函数族, 但缺点是其物理意义不明确. 此外, 上述文献都没有考虑航天器姿态机动中的姿态约束.

    在轨航天器姿态机动时, 存在多种约束, 其中一种典型约束是航天器的禁止或强制指向约束, 即要求航天器在机动中始终避开或始终指向某个方向. 例如, 航天器光敏感器的指向应当始终避开明亮天体[12-13]. 指向约束下实现姿态机动的方案大致可分为基于规划的跟踪控制方法[14-15] 和基于避障势函数方法[13, 16]两种. 基于规划的跟踪控制方法将原问题分解成2个子问题进行单独解决. 首先, 利用优化方法规划出合理的姿态机动路径; 然后, 设计跟踪控制器跟踪所给路径. 文献[15]整体性地考虑了姿态约束、机动过程的指向安全裕度和力矩饱和问题, 提出预设性能的“规划$+ $跟踪”控制方法, 系统性地实现了航天器姿态受限机动. 此类基于规划方法的姿态机动路径具备可预见性, 但缺点是结构复杂, 优化算法和控制器性能都会影响到航天器的实际指向. 相比之下, 基于避障势函数方法更加简单和有效. 该方法将目标姿态设计为极小值点, 并将势函数的梯度引入控制器设计, 当靠近姿态限制时, 势函数值和梯度急剧增加, 以此规避姿态限制. 文献[13]使用分数形式的避障势函数, 处理禁止指向约束问题; 文献[16]使用对数形式的避障势函数解决了同时具有禁止和强制指向约束的姿态机动问题. 然而, 势函数方法严重依赖势函数的凸性假设, 机动路径不可预见且单一势函数控制存在临界点问题, 航天器姿态可能会收敛到势函数的极小值, 而非目标姿态.

    为解决上述问题, 本文在考虑姿态约束情形下, 基于旋转矩阵设计用于混合控制的势函数族, 该势函数族能够对任意的初始姿态保证姿态收敛到目标姿态. 本文的主要贡献如下: 1)建立姿态禁止指向约束模型, 考虑机动安全裕度, 为每个指向约束设计软约束区域. 在该区域内, 将分数/对数形式排斥势函数引入协同势函数设计中, 设计保证全局收敛的避障势函数. 2)针对所设计的势函数, 给出软约束区域内临界点分布情况和势函数参数的关系, 指出当势函数参数满足给定不等式组时, 航天器不会在避障过程中陷入非期望的临界点. 3)将上述势函数引入控制器设计, 实现姿态受限下的全局收敛.

    本文内容结构如下: 第1节对使用到的符号和性质进行说明, 建立航天器的误差姿态运动模型和姿态约束模型; 第2节给出基于“角度扰动”构建协同势函数族的方法, 将该协同势函数与受限制姿态机动问题结合, 分别讨论分数形式和对数形式下, 势函数在软约束区域内临界点分布和势函数参数的关系; 第3节以比例−微分控制为例, 进行仿真验证.

    本文以旋转矩阵$ {\boldsymbol{R}} $作为姿态参数展开研究, 旋转矩阵形成的三维特殊正交群记作$SO(3) := \{ {\boldsymbol{R}} \in $ $ {{\bf{R}}^{3 \times 3}}|{{\boldsymbol{R}}^{\rm{T}}}{\boldsymbol{R}} = {\boldsymbol{R}}{{\boldsymbol{R}}^{\rm{T}}} = {{\boldsymbol{I}}},\det ({\boldsymbol{R}}) = 1 \}$. $ {{\boldsymbol{I}}}\in {{\bf{R}}^{3 \times 3}} $是单位矩阵. $ SO(3) $以矩阵乘法作为群乘法时, 是一个李群, 其对应的李代数用$ \mathfrak{so}\left( 3 \right) $表示, $ \mathfrak{so}\left( 3 \right) $的群元是三阶反对称矩阵, 即$\mathfrak{s}\mathfrak{o}\left( 3 \right) = \{ {\boldsymbol{S}} \in {{\bf{R}}^{3 \times 3}}: {{\boldsymbol{S}}^{\rm{T}}} = -{\boldsymbol{S}} \}$. 对于任意向量$ {\boldsymbol{x}} \in {{\bf{R}}^3} $, 定义线性映射${\left[ \cdot \right]^ \times }: {{\bf{R}}^3} \to \mathfrak{s}\mathfrak{o}\left( 3 \right)$和它的逆映射${\left[ \cdot \right]^ \vee }:\mathfrak{s}\mathfrak{o}\left( 3 \right) \to {{\bf{R}}^3}$为:

    $$ \begin{array}{l} {\left[ {\boldsymbol{x}} \right]^ \times } = \left[ {\begin{array}{*{20}{c}} 0&{ - {x_3}}&{{x_2}}\\ {{x_3}}&0&{ - {x_1}}\\ { - {x_2}}&{{x_1}}&0 \end{array}} \right] , {\left[ {{{\left[ {\boldsymbol{x}} \right]}^ \times }} \right]^ \vee } = {\boldsymbol{x}} \end{array} $$ (1)

    在欧氏空间中, $ n $维单位球面为${S^n}: = \{ {\boldsymbol{x}} \in {{\bf{R}}^{n + 1}} | {\boldsymbol{x}}^{\rm{T}}{\boldsymbol{x}} = 1 \}$. 为了强调矩阵$ {\boldsymbol{B}} $的某个元素时, 写${\boldsymbol{B}} = [ {{b_{ij}}} ]$. 对于任意向量$ {\boldsymbol{x}},{\boldsymbol{y}} \in {{\bf{R}}^n} $和矩阵${\boldsymbol{B}},{{\boldsymbol{C}}} \in {{\bf{R}}^{m \times n}}$, 内积定义为$\langle {{\boldsymbol{x}},{\boldsymbol{y}}} \rangle : = {\boldsymbol{x}}^{\rm{T}}{\boldsymbol{y}},\langle {{\boldsymbol{B}},{{\boldsymbol{C}}}} \rangle : = {\rm{tr}}( {{{\boldsymbol{B}}^{\rm{T}}}{{\boldsymbol{C}}}} )$. 向量的2范数是$ \| {\boldsymbol{x}} \| = \sqrt {\langle {{\boldsymbol{x}},{\boldsymbol{x}}} \rangle } $. 矩阵的${\rm{F }}$范数为${\| {\boldsymbol{B}} \|_{\rm{F}}} = \sqrt {\langle {{\boldsymbol{B}},{\boldsymbol{B}}} \rangle }$. 定义映射 $\psi :{{\bf{R}}^{3 \times 3}} \;\to\; {{\bf{R}}^3}$为$\psi ( {\boldsymbol{B}} ) \; $$= 0.5[ {\boldsymbol{B}} - {{\boldsymbol{B}}^{\rm{T}}} ]^ \vee$, $\psi ( { - {\boldsymbol{B}}} ) = \;\psi ( {{{\boldsymbol{B}}^{\rm{T}}}} ) = \;- \psi ( {\boldsymbol{B}} )$. 用$ \lambda_{\max}^{{\boldsymbol{A}}} $和$ \lambda _{\min}^{{\boldsymbol{A}}} $分别表示矩阵$ {{\boldsymbol{A}}} $的最大和最小特征值. 对于可微映射$ G:{\bf{R}} \to {{\bf{R}}^{m \times n}},H:{{\bf{R}}^{m \times n}} \to {\bf{R}} $, 定义$\dot G( t ) = { {\frac{{{\rm{d}}G( t )}}{{{\rm{d}}t}}} } = \left[\frac{{{\rm{d}}{G_{ij}}( t )}}{{{\rm{d}}t}}\right], {{\nabla H( {{\boldsymbol{X}}} )} } =\left [ \frac{{\partial H( {{\boldsymbol{X}}} )}}{{\partial {{{\boldsymbol{X}}}_{ij}}}}\right]$. 假设可微映射$ V $是一个复合映射, 即$ V = H \circ G $, 则根据链导法则, 有$ \dot V( t ) = \langle {\nabla H( {G( t )} ),\dot G( t )} \rangle $.

    对$ \forall {{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}},{{{\boldsymbol{R}}_3}} \in SO\left( 3 \right) $, 如果映射$\Phi :SO\left( 3 \right) \times SO\left( 3 \right) \to {{\bf{R}}_{ \ge 0}}$满足:

    1)$ \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}}} \right) = 0 \Leftrightarrow {{{\boldsymbol{R}}_1}} = {{{\boldsymbol{R}}_2}} $;

    2)$ \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}}} \right) = \Phi \left( {{{{\boldsymbol{R}}_2}},{{{\boldsymbol{R}}_1}}} \right) $;

    3)$ \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}}} \right) \le \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_3}}} \right) + \Phi \left( {{{{\boldsymbol{R}}_2}},{{{\boldsymbol{R}}_3}}} \right) $.

    则映射$ \Phi $被称为$ SO\left( 3 \right) $上的度量. 如果度量满足$\Phi \left( {{\boldsymbol{R}}{{{\boldsymbol{R}}_1}},{\boldsymbol{R}}{{{\boldsymbol{R}}_2}}} \right) = \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}}} \right), \forall {\boldsymbol{R}}\in SO\left( 3 \right)$, 则称之为左不变的; 如果满足$ \Phi \left( {{{{\boldsymbol{R}}_1}}{\boldsymbol{R}},{{{\boldsymbol{R}}_2}}{\boldsymbol{R}}} \right) = \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}}} \right), \forall {\boldsymbol{R}} \in SO\left( 3 \right) $, 则称之为右不变的. 本文用到的$ SO\left( 3 \right) $上的度量是角距离度量, 即:

    $$ \begin{split} \Phi \left( {{{{\boldsymbol{R}}_1}},{{{\boldsymbol{R}}_2}}} \right) = \;& \arccos \left( {\frac{1}{2}\left( {{\rm{tr}}\left( {{{{\boldsymbol{R}}_1}}{{\boldsymbol{R}}_2^{\rm{T}}}} \right) - 1} \right)} \right) = \\ & \mathop {\sup }\limits_{{\boldsymbol{u}} \in {S^2}} \arccos \left( {\langle {{{{\boldsymbol{R}}_1}}{\boldsymbol{u}},{{{\boldsymbol{R}}_2}}{\boldsymbol{u}}} \rangle } \right) \end{split} $$ (2)

    角距离度量是一种双平移不变度量[17]. 利用平移不变性, 角距离度量也可以写为:

    $$ \begin{array}{l} \Phi \left( {{{{\boldsymbol{R}}_1}}{{\boldsymbol{R}}_2^{\rm{T}}},{{\boldsymbol{I}}}} \right) = \mathop {\sup }\limits_{{\boldsymbol{u}} \in {S^2}} \arccos \left( {\langle {{{{\boldsymbol{R}}_1}}{{\boldsymbol{R}}_2^{\rm{T}}}{\boldsymbol{u}},{\boldsymbol{u}}} \rangle } \right) \end{array} $$ (3)

    任意旋转矩阵$ {\boldsymbol{R}} $相对于单位元$ {\boldsymbol{I}} $的角距离简记为${\Phi _{{\boldsymbol{R}}}} = \Phi ( {{\boldsymbol{R}},{\boldsymbol{I}}} )$, $ {\Phi _{{\boldsymbol{R}}}} \in [ {0,\pi } ] $是作用于向量时, 能使向量转动的最大角度. 给定一个转角$ \theta \in {\bf{R}} $和转轴$ {\boldsymbol{u}} \in {S^2} $, 则${{\rm{e}}^{\theta {{[ {\boldsymbol{u}} ]}^ \times }}} \in SO( 3 )$且${\boldsymbol{R}}( {\theta ,{\boldsymbol{u}}} ): = {{\rm{e}}^{\theta {{[ {\boldsymbol{u}} ]}^ \times }}} = {{\boldsymbol{I}}} + \sin \theta {[ {\boldsymbol{u}} ]^ \times } + ( {1 - \cos \theta } ){( {{{[ {\boldsymbol{u}} ]}^ \times }} )^2}$, 其中$ {{\rm{e}}^{\theta {{[ {\boldsymbol{u}} ]}^ \times }}} $是矩阵$ {\theta {{[ {\boldsymbol{u}} ]}^ \times }} $的矩阵指数. 任意给定的旋转矩阵都可以用$ {{\rm{e}}^{\theta {{[ {\boldsymbol{u}} ]}^ \times }}} $表示[18], 并且$\theta = \Phi_{{\boldsymbol{R}}} = \arccos ( 0.5({\rm{tr}}( {\boldsymbol{R}} ) - 1)).$

    令集合$ D \subset SO(3) $, 如果连续可微的函数$V: SO\left( 3 \right) \to {{\bf{R}}_{ \ge 0}}$满足:

    1)$ \forall {\boldsymbol{R}} \in SO\left( 3 \right)\backslash D,V\left( {\boldsymbol{R}} \right) > 0 $;

    2)$ V\left( D\right) = 0 $.

    则$ V $被称为$ SO(3) $上相对于集合$ D $的势函数.

    本节介绍对后续分析有帮助的性质和引理. 令$ {{\boldsymbol{A}}} $是一个实对称矩阵, 定义势函数:

    $$ \begin{array}{l} {P_{{{\boldsymbol A}}}}\left( {\boldsymbol{R}} \right) = {\rm{tr}}\left( {{{\boldsymbol{A}}}\left( {{{\boldsymbol{I}}} - {\boldsymbol{R}}} \right)} \right) \end{array} $$ (4)

    令$ {{\boldsymbol{A}}_1} = 0.5( {{\rm{tr}}( {{\boldsymbol{A}}} ){{\boldsymbol{I}}} - {{\boldsymbol{A}}}} ) $, ${{\boldsymbol{A}}_2} = {\rm{tr}}( {{{{\boldsymbol{A}}_1^2} }} ){{\boldsymbol{I}}} - 2{{{\boldsymbol{A}}_1^2} }$, ${{\boldsymbol{A}}_3} = 0.5( {{\rm{tr}}( {{{\boldsymbol{A}}_2} } ){{\boldsymbol{I}}} - {{\boldsymbol{A}}_2} } )$, 对$\forall {\boldsymbol{x}},{\boldsymbol{y}} \in {{\bf{R}}^3},{\boldsymbol{B}} \in {{\bf{R}}^{3 \times 3}}, {\boldsymbol{R}} \in SO( 3 ),$ 有下列性质成立:

    1)${\left[ [{{\boldsymbol{x}}]^\times {\boldsymbol{y}}} \right]^ \times } = {\boldsymbol{y}}{{\boldsymbol{x}}^{\rm{T}}} - {\boldsymbol{x}}{{\boldsymbol{y}}^{\rm{T}}}$;

    2)${\rm{tr}}( {{\boldsymbol{B}}{{[ {\boldsymbol{x}} ]}^ \times }} ) = - {{\boldsymbol{x}}^{\rm{T}}}{[ {{\boldsymbol{B}} - {{\boldsymbol{B}}^{\rm{T}}}} ]^ \vee } = - 2{{\boldsymbol{x}}^{\rm{T}}}\psi ( {\boldsymbol{B}} )$;

    3)$ 4\lambda _{\min }^{{{\boldsymbol{A}}_1}}h \le {P_{{\boldsymbol{A}}}}\left( {\boldsymbol{R}} \right) \le 4\lambda _{\max }^{{{\boldsymbol{A}}_1}}h $;

    4)${\left\| {\psi \left( {{{\boldsymbol{A}}}{{\boldsymbol{R}}}} \right)} \right\|^2} \le \left( {1 - {\xi ^2}h} \right)4\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }h$.

    式中, $\xi = {{\lambda _{\min }^{{{{{\boldsymbol{A}}}}_1} }} / {\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1} }}}, h = {\sin ^2}{\left( {\Phi \left( {{\boldsymbol{R}},{\boldsymbol{I}}} \right)/2} \right)}$. 性质1)和性质2)通过展开计算即可证明, 性质3)和性质4)的证明见文献[19].

    下面列出在本文中使用的引理.

    引理1. 令$ 0 < \alpha < \alpha + \Delta \alpha < {\pi }/{2}, x \in ( \alpha ,\alpha \;+ \Delta \alpha ] $, 则函数$ f( x ) = {2\sin x/( {\cos \alpha - \cos x} )} $的取值范围为$ [ {f( {\alpha + \Delta \alpha } ), + \infty } ) $.

    证明. 在区间$ ( {\alpha ,\alpha + \Delta \alpha } ] $上, 函数$ f( x ) $相对于$ x $的导数为$\frac{{\rm{d}}f( x )}{{\rm{d}}x }= {2( {\cos x\cos \alpha - 1} )/{{( {\cos \alpha - \cos x} )}^2}} < 0$. 因此, 函数$ f( x ) $严格递减, 即$f( {\alpha + \Delta \alpha } ) \le f( x ) < \infty$.

    引理2. 当常数$ a $满足$ 2{a^2} < 1 $时, 函数$f\left( x \right) = \sqrt {\left( {1 - {a^2}{{\sin }^2}\left( {x/2} \right)} \right){{\sin }^2}\left( {x/2} \right)} ,x \in \left( {0,\pi } \right]$的取值范围为$ \left( {f\left( 0 \right),f\left( \pi \right)} \right] $.

    证明. 在区间$ \left( {0,\pi }\right] $上, 函数$ f\left( x \right) $相对于$ x $的导数为$\frac{{\rm{d}}f\left( x \right)}{{\rm{d}}x }= \left( {1 - 2{a^2}{{\sin }^2}\left( {x/2} \right)} \right)\sin (x/4)$. 由于$2{a^2} < 1$, 函数$ f\left( x \right) $严格递增, 即$ f\left( 0 \right) < f\left( x \right) \le f\left( \pi \right) $.

    引理3. 令常数$ a,b $满足$ 0 \le a<\pi, 0<b<1 $, 则函数$ f\left( x \right) = x - 2b{\sin ^2}\left( {x/2} \right),x \in \left( {a,\pi } \right] $的取值范围为$ \left( {f(a),f(\pi) } \right] $.

    证明. 在区间$ \left[ {a,\pi }\right] $上, 函数$ f\left( x \right) $相对于$ x $的导数为$ \frac{{\rm{d}}f\left( x \right)}{{\rm{d}}x } = 1 - b\sin x $. 由于$ b < 1 $, 函数$ f\left( x \right) $严格递增, 即$ f\left( x \right) \ge f\left( a \right) $.

    假设航天器本体系相对于惯性系的姿态用$ {\boldsymbol{R}}_b \in SO(3) $表示, 姿态机动的目标姿态用$ {\boldsymbol{R}}_d \in SO(3) $表示. 令$ {\boldsymbol{\omega}} $表示航天器本体系下的角速度矢量, $ {\boldsymbol{\omega}}_d $表示期望角速度. 假设惯性矩阵$ {\boldsymbol{J}} $是对角阵, 则航天器的运动学和动力学方程为:

    $$ \begin{array}{l} {{\dot {{\boldsymbol{R}}}}_b} ={{\boldsymbol{R}}_b}{\left[ {\boldsymbol{\omega}} \right]^ \times } {} \end{array} $$ (5)
    $$ \begin{array}{l} {\boldsymbol{\boldsymbol{J}}} \dot{{\boldsymbol{\omega}}} + {\left[ {\boldsymbol{\omega}} \right]^ \times } {\boldsymbol{\boldsymbol{J}}} {\boldsymbol{\omega}} = {\boldsymbol{\tau}} {} \end{array} $$ (6)

    式中, $ {\boldsymbol{\tau}} $表示控制力矩. 定义的姿态误差为$ {\boldsymbol{R}}_e = {\boldsymbol{R}}_d^{\rm{T}}{\boldsymbol{R}}_b $, 定义角速度误差是$ {{\boldsymbol{\omega}} _e} = {\boldsymbol{\omega}} - {\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{\omega}} _d} $, 则误差形式的运动学方程和动力学方程为:

    $$ \begin{array}{l} {{\dot {{\boldsymbol{R}}}}_e} = {{\boldsymbol{R}}_e}{\left[ {\boldsymbol{\omega}}_e \right]^ \times } \end{array} $$ (7)
    $$ \begin{array}{l} {\boldsymbol{\boldsymbol{J}}}{{\dot {{\boldsymbol{\omega}}}}_e} = {\boldsymbol{H}}{{\boldsymbol{\omega}} _e} - {\boldsymbol{F}} + {\boldsymbol{\tau}} \end{array} $$ (8)

    式中, ${\boldsymbol{F}} = {\boldsymbol{\boldsymbol{J}}}{{\boldsymbol{R}}_e^{\rm{T}}}{{\dot {{\boldsymbol{\omega}}}}_d} +{[ {{{\boldsymbol{R}}_e^{\rm{T}}}{{\boldsymbol{\omega}} _d}} ]^ \times }{\boldsymbol{J}}{{\boldsymbol{R}}_e^{\rm{T}}}{{\boldsymbol{\omega}} _d}$, ${\boldsymbol{H}} = {[ {{\boldsymbol{J}}{{\boldsymbol{\omega}} _e}} ]^ \times } + {[ {{\boldsymbol{J}}{{\boldsymbol{R}}_e^{\rm{T}}}{{\boldsymbol{\omega}} _d}} ]^ \times } - {\boldsymbol{J}}{[ {{{\boldsymbol{R}}_e^{\rm{T}}}{{\boldsymbol{\omega}} _d}} ]^ \times } - {[ {{{\boldsymbol{R}}_e^{\rm{T}}}{{\boldsymbol{\omega}} _d}} ]^ \times }{\boldsymbol{J}}$是一个反对称矩阵. 姿态控制的目标是使航天器姿态$ {\boldsymbol{R}}_b $趋向于给定目标姿态$ {\boldsymbol{R}}_d $, 或等价地说, 是使航天器的姿态误差$ {\boldsymbol{R}}_e $趋向于$ {\boldsymbol{I}} $.

    航天器姿态机动过程中存在各种限制, 本文考虑一种几何视场约束, 即航天器上某个指向需要避开特定区域. 如图1 所示, 图1中下标$ B $表示本体系, 下标$ I $表示惯性系. 令单位向量$ {\boldsymbol{r}} \in {\bf{R}}^{3} $表示与航天器固连的敏感指向, 单位向量$ {\boldsymbol{v}} \in {\bf{R}}^{3} $表示惯性系下固定的不期望指向(例如明亮的星体), 航天器在姿态机动过程中, 单位向量$ {\boldsymbol{r}} $和$ {\boldsymbol{v}} $夹角应满足:

    图 1  姿态限制示意图
    Fig. 1  Illustration of attitude constraint
    $$ \langle {{{\boldsymbol{R}}_b}{\boldsymbol{r}},{\boldsymbol{v}}} \rangle \le \cos \left( {{\alpha} } \right), \alpha \in \left( {0,\frac{\pi }{2}} \right) $$ (9)

    因此, 本文中势函数设计的目标是在所设计的势函数指导下, 航天器的姿态误差$ {\boldsymbol{R}}_e $趋向于$ {\boldsymbol{I}} $, 同时始终满足式(9).

    利用内积的定义和性质2), 势函数$ V\left( {{\boldsymbol{R}}\left( t \right)} \right) $对时间的导数是:

    $$ \begin{split} \dot V\left( {{\boldsymbol{R}}\left( t \right)} \right) = \;& \langle {\nabla V\left( {{\boldsymbol{R}}\left( t \right)} \right),{\boldsymbol{R}}{{\left[ {\boldsymbol{\omega}} \right]}^ \times }} \rangle = \\ &{\rm{tr}}\left( {\nabla {V^{\rm{T}}}{\boldsymbol{R}}{{\left[ {\boldsymbol{\omega}} \right]}^ \times }} \right)=\\ & 2\langle {{\boldsymbol{\omega}},\psi \left( {{{\boldsymbol{R}}^{\rm{T}}}\nabla V} \right)} \rangle \end{split} $$ (10)

    由此可知, 当$ \psi \left( {{{\boldsymbol{R}}^{\rm{T}}}V} \right) = 0 $时, 无论角速度$ {\boldsymbol{\omega}} $为何值, $ V $不会发生任何改变, 即姿态停止向目标姿态集合$ D $收敛. 因此, 定义势函数$ V $的临界点如下:

    $$ \begin{array}{l} {\rm{Crit}}V = \left\{ {\left. {{\boldsymbol{R}} \in SO\left( 3 \right)} \right|\psi \left( {{{\boldsymbol{R}}^{\rm{T}}}\nabla V\left( {\boldsymbol{R}} \right)} \right) = 0} \right\} \end{array} $$ (11)

    例如, 对于势函数(4), 有下列结果:

    $$ \begin{array}{l} \begin{cases} \nabla {P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}} \right) = - {\boldsymbol{A}}\\ {\rm{Crit}}{P_{{{{\boldsymbol{A}}}}}} = \left\{ {\left. {{\boldsymbol{R}} \in SO\left( 3 \right)} \right|{{\boldsymbol{R}}^{\rm{T}}}{\boldsymbol{A}} = {\boldsymbol{A}}{\boldsymbol{R}}} \right\} \end{cases} \end{array} $$ (12)

    通过解上述方程可知, 势函数(4)共有4个临界点[8]:

    $$ \begin{array}{l} {\rm{Crit}}{P_{{{{\boldsymbol{A}}}}}} = \left\{ {\boldsymbol{I}} \right\} \cup {{\boldsymbol{R}}}\left( {\pi ,{\cal{E}}\left( {{\boldsymbol{A}}} \right)} \right) \end{array} $$ (13)

    式中, ${\cal{E}}\left( {{\boldsymbol{A}}} \right)$表示$ {\boldsymbol{A}} $的所有单位特征向量组成的集合.

    对任意的势函数$ V $, 都存在至少4个临界点[20], 基于势函数梯度的反馈控制律在靠近临界点时, 逐渐减小直至为零. 为了实现全局收敛的结果, 本文使用一族协同的共单位元势函数进行控制. 给定一个有限的指标集$Q \subset {{\bf N}}$, ${\left\{ {{V_q}} \right\}_{q \in {{Q}}}}$是一族势函数, 如果单位元$ {\boldsymbol{I}} $是所有势函数公共的临界点, 则称这样的一族势函数为共单位元的势函数. 对于势函数族${\left\{ {{V_q}} \right\}_{q \in {{Q}}}}$, 如果存在一个常数$ \delta > 0 $, 使得:

    $$ \begin{array}{l} \overline \delta : = \mathop {\min }\limits_{ \begin{subarray}{c} q \in {{Q}}\\ {\boldsymbol{R}} \in {\rm{Crit}}{V_q}\backslash \left\{ {\boldsymbol{I}} \right\} \end{subarray} } \left( {\mathop {\max }\limits_{\begin{subarray}{c} p \in {{Q}} \\ p \ne q\end{subarray} } {V_q}\left( {\boldsymbol{R}} \right) - {V_p}\left( {\boldsymbol{R}} \right)} \right) > \delta \end{array} $$ (14)

    则称势函数族${\left\{ {{V_q}} \right\}_{q \in {{Q}}}}$是协同的. 粗略地, “共单位元”意味着在这一族势函数控制下, 姿态有相同的收敛方向(单位元$ {\boldsymbol{I}} $); “协同”意味着在某个势函数临界时, 存在另一个非临界且取值更低的势函数, 如果切换为非临界的势函数进行控制, 姿态将继续向目标姿态收敛.

    为了说明协同势函数在姿态控制中的作用和切换逻辑, 假设航天器的一条姿态轨迹经过2个非期望临界点的邻域, 这2个临界点分别属于${P_{1,{{{\boldsymbol A}}}}}$和${P_{2,{{{\boldsymbol A}}}}}$, 则势函数随时间变化曲线和切换时各势函数与临界点关系如图2 所示, 其中状态$ N $和$ M $在第3节定义. 在单一势函数${P_{1,{{{\boldsymbol A}}}}}$(或${P_{2,{{{\boldsymbol A}}}}}$)的控制下, 轨迹在$ t_1 $(或$ t_2 $)时刻, 进入势函数所属的非期望临界点邻域, 此后基于势函数梯度的姿态控制项将减小, 甚至为零. 为避免该问题导致收敛过慢, 将在进入时刻进行切换. 如图2所示, 在$ t_3 $时刻, 姿态轨迹处于势函数${P_{1,{{{\boldsymbol A}}}}}$的非期望临界点, 势函数停止收敛, 此时存在势函数${P_{2,{{{{\boldsymbol{A}}}}}}}$在该点不临界且其取值更低. 因此, 在混合控制中, 根据势函数差值在$ t_2 $时刻提前进行控制切换, 由${P_{1,{{{\boldsymbol A}}}}}$跳转到${P_{2,{{{{\boldsymbol{A}}}}}}}$; 同样在$ t_1 $时刻, 由${P_{2,{{{{\boldsymbol{A}}}}}}}$跳转到${P_{1,{{{{\boldsymbol{A}}}}}}}$, 最终混合控制中势函数的变化曲线如图2 橙色粗线所示.

    图 2  使用协同势函数族的姿态控制示意图
    Fig. 2  Illustration of attitude control by means of synergistic potential functions

    为了构建协同势函数族, 定义映射$\mathcal{T}:SO\left( 3 \right) \to SO\left( 3 \right)$为:

    $$ \begin{array}{l} {\cal{T}}\left( {{\boldsymbol{R}},k,{P_{{\boldsymbol{A}}}},{\boldsymbol{u}}} \right) = {{\rm{e}}^{k{P_{{\boldsymbol{A}}}}\left( {\boldsymbol{R}} \right){{\left[ {\boldsymbol{u}} \right]}^ \times }}}{\boldsymbol{R}} \end{array} $$ (15)

    式中, $ k \in {\bf{R}},{\boldsymbol{u}} \in {S^2} $. 映射$ {\cal{T}} $的作用是对旋转矩阵$ {\boldsymbol{R}} $左乘“扰动”旋转矩阵, 该“扰动”矩阵由指数映射生成, 转轴$ {\boldsymbol{u}} $和转角$ k{P_{{\boldsymbol{A}}}}( {\boldsymbol{R}} ) $是待设计的轴角参数. 因此, 简记${\boldsymbol{R}}_c = {\cal{T}}( {{\boldsymbol{R}},k,{P_{{{\boldsymbol A}}}},{\boldsymbol{u}}} ) = {\boldsymbol{R}}_{{\cal{T}}}{\boldsymbol{R}}$, 其中$ {\boldsymbol{R}}_{{\cal{T}}} $是“扰动”旋转矩阵. 令$Q = \{ {1,2} \},q \in {{Q}}$, 当参数$ k $满足$0 < \sqrt 2 k{\| {\boldsymbol{A}} \|_{\rm{F}}} < 1$时, 协同势函数族$\{ {P_{q,{{{\boldsymbol A}}}}}({\boldsymbol{R}} ) \}_{q \in {{Q}}}$可以按如下形式构建[8]:

    $$ \begin{array}{l} {P_{q,{{{\boldsymbol A}}}}}\left({\boldsymbol{R}} \right) = {P_{{{\boldsymbol A}}}}\left( {{\cal{T}}\left( {{\boldsymbol{R}},{{\left( { - 1} \right)}^q}k,{\boldsymbol{u}},{P_{{{\boldsymbol A}}}}} \right)} \right) \end{array} $$ (16)

    并且, 有:

    $$ \begin{array}{l} \left\{ \begin{aligned} &\psi \left( {{{\boldsymbol{R}}^{\rm{T}}}\nabla {P_{q,{\boldsymbol{A}}}}\left( {\boldsymbol{R}} \right)} \right) = \Theta ^{\rm{T}}{\left( {\boldsymbol{R}} \right)}\psi \left( {{{{\boldsymbol{R}}}_c^{\rm{T}}}\nabla {P_{{{{\boldsymbol{A}}}}}}\left( {{\boldsymbol{R}}}_c \right)} \right)\\ &{\cal{T}}\left( {{\rm{Crit}}{P_{q,{{{{\boldsymbol{A}}}}}}}} \right) = {\rm{Crit}}{P_{{{{\boldsymbol{A}}}}}} \end{aligned} \right. \end{array} $$ (17)

    式中, $\Theta \left( {\boldsymbol{R}} \right) = {\boldsymbol{I}} + 2k{{\boldsymbol{R}}^{\rm{T}}}{\boldsymbol{u}}\psi {\left( {\nabla P_{{{{\boldsymbol{A}}}}}\left( {\boldsymbol{R}} \right){{\boldsymbol{R}}^{\rm{T}}}} \right)^{\rm{T}}}{\boldsymbol{R}}$.

    引理4. 当参数$ k $和$ {{\boldsymbol{A}}} $ 满足 $2k\lambda _{\max }^{{{\boldsymbol{A}}_1} } < 1$时, 有${\cal{T}}\left( {{\boldsymbol{R}},k,{P_{{{\boldsymbol A}}}},{\boldsymbol{u}}} \right) \ne {\boldsymbol{I}},\forall {\boldsymbol{R}} \in SO\left( 3 \right)\backslash {\boldsymbol{I}}$.

    证明(反证法). 假设$ \exists {\boldsymbol{R}} \in SO( 3 )\backslash {\boldsymbol{I}} $, 使得$ {\boldsymbol{R}}_c = {\boldsymbol{I}} $, 其中${\boldsymbol{R}}_c = {{\rm{e}}^{k{P_{{\boldsymbol{A}}}}( {\boldsymbol{R}} ){{[ {\boldsymbol{u}} ]}^ \times }}}{\boldsymbol{R}} = {\boldsymbol{R}}_{{\cal{T}}}{\boldsymbol{R}}$, 则有$ {\boldsymbol{R}}_{{\cal{T}}}={\boldsymbol{R}}^{\rm{T}} $. 考虑到$ \Phi_{{\boldsymbol{R}}^{\rm{T}}} = \Phi_{{\boldsymbol{R}}} $, 则假设成立时, 有$\Phi_{{\boldsymbol{R}}_{{\cal{T}}}}= k{P_{{{\boldsymbol A}}}}( {\boldsymbol{R}} ) = \Phi_{{\boldsymbol{R}}}$. 由于$ {\boldsymbol{R}} \ne {\boldsymbol{I}} $, $ \Phi_{{\boldsymbol{R}}} > 0 $, 再由第1.2节的性质3), 可得:

    $$ \begin{array}{l} 0 < 4k\lambda _{\min }^{{{{\boldsymbol A}}_1}}h \le k{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}} \right) \le 4k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1}}h \end{array} $$ (18)

    式中, $h = {\sin ^2} \left( {{\Phi _{{{\boldsymbol R}}}}/2} \right)$. 构造函数:

    $$ f\left( \Phi_{{{\boldsymbol R}}}\right) = \Phi_{{\boldsymbol{R}}} - 4k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1}}{\sin ^2}\left( \frac{{\Phi _{{{\boldsymbol R}}}}}{2} \right) $$ (19)

    根据引理3, 当$2k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1} } < 1$时, $f(\Phi_{{{\boldsymbol R}}}) > f(0)= 0$. 这与$k{P_{{{\boldsymbol A}}}}\left( {\boldsymbol{R}} \right) = \Phi_{{{\boldsymbol R}}}$矛盾, 因此假设不成立.

    引理4对映射$ {\cal{T}} $的参数要求如图3所示. $\overline {kP_{{{\boldsymbol A}}}\left( {\boldsymbol{R}} \right)}$和$\underline {kP_{{{\boldsymbol A}}}\left( {\boldsymbol{R}} \right)}$分别为根据性质3)计算得到的$kP_{{{\boldsymbol A}}}\left( {\boldsymbol{R}} \right)$的上限和下限随角距离$\Phi_{{{\boldsymbol R}}}$的变化曲线. 对$ \forall {\boldsymbol{R}} \in SO\left( 3 \right)\backslash {\boldsymbol{I}} $, 有$\Phi_{{{\boldsymbol R}}} > k{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}} \right)$始终成立, 即非单位元旋转矩阵经过映射$ {\cal{T}} $“扰动”后, 仍然是非单位元旋转矩阵, 不会到达临界点$ {\boldsymbol{I}} $.

    图 3  $k{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}} \right)$取值范围随$\Phi_{{{\boldsymbol R}}}$的变化曲线
    Fig. 3  The change curve of the range of$k{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}} \right)$with$\Phi_{{{\boldsymbol R}}}$

    假设航天器初始姿态位于姿态障碍以外, 航天器误差运动模型和势函数设计目标如第1.2节所述. 定义$ SO(3) $的一个子集为:

    $$ \begin{split} E = \;& \left.\{\left. {{{\boldsymbol{R}}_e} \in SO\left( 3 \right)} \right| \right.\\ &\left. \arccos \langle {{{\boldsymbol{R}}_d}{{\boldsymbol{R}}_e}{\boldsymbol{r}},{\boldsymbol{v}}} \rangle \in \left( {\alpha ,\alpha + \Delta \alpha } \right) \right.\} \end{split} $$ (20)

    该子集定义了姿态限制的软约束区域, 进入该区域时势函数将进行避障切换. 为在全局收敛情况下完成受限制姿态机动, 本文设计了如下形式的势函数用于姿态机动问题:

    $$ \begin{array}{l} V_q\left( {{{\boldsymbol{R}}_e}} \right) = \left\{ \begin{aligned} & {{P_{q,{{{\boldsymbol{A}}}}}}\left( {{{\boldsymbol{R}}_e}} \right){P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right)}, &{{{\boldsymbol{R}}_e} \in E}\\ &{{P_{q,{{{\boldsymbol{A}}}}}}\left( {{{\boldsymbol{R}}_e}} \right)}, &{{{\boldsymbol{R}}_e} \notin E} \end{aligned} \right. \end{array} $$ (21)

    式中, $ \alpha + \Delta \alpha < \pi /2, \Delta \alpha $是待设计参数, ${{P_{q,{{{\boldsymbol{A}}}}}}\left( {{{\boldsymbol{R}}_e}} \right)}$是针对目标姿态的吸引项势函数, ${P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right)$是待设计的用于避开姿态障碍的排斥项势函数, 以分数形式的${P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right)$为例:

    $$ {P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right) = \frac{b }{{{{\left( {\cos \alpha - {{\boldsymbol{r}}^{\rm{T}}}{\boldsymbol{R}}_b^{\rm{T}}{\boldsymbol{v}}} \right)}^a}}} = \frac{b }{{{{\left( {\cos \alpha - {{\boldsymbol{r}}^{\rm{T}}}{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}} \right)}^a}}} $$ (22)

    式中, $ {{\boldsymbol{v}}_d} = {\boldsymbol{R}}_d^{\rm{T}}{\boldsymbol{v}} $; $ a $和$b $是两个正的参数, 文献[13, 21-22]取$ a = 1 $或$ a = 2 $, 本文考虑更一般情况即只要求$ a > 0 $. 在实际的姿态机动中, 可能存在多个姿态限制, 本文方法适用于单敏感轴多约束问题. 具体来说, 为每个指向约束设计了不相交的软约束区域, 当航天器姿态进入软约束区域内, 即敏感轴在惯性系下指向$ {{\boldsymbol{R}}_b}{\boldsymbol{r}} $和禁止指向$ {\boldsymbol{v}} $夹角小于$ \alpha + \Delta \alpha $时, 切换为避障势函数, 因此只需讨论单个约束下如何实现姿态机动即可. 对于某个给定的姿态限制有下述引理成立.

    引理5. 对$\forall {{\boldsymbol{R}}_e} \in E$, $\Phi_{{\boldsymbol{R}}_e} \in [ {\varepsilon ,\pi } ], \varepsilon = \arccos \langle {{\boldsymbol{R}}_d}{\boldsymbol{r}}, {\boldsymbol{v}} \rangle - \left( {\alpha + \Delta \alpha } \right)$.

    证明. 旋转矩阵$ {\boldsymbol{R}} $作用于单位向量后, 会使单位向量发生旋转. 角距离的定义指出, 旋转前/后向量的夹角不超过旋转矩阵$ {\boldsymbol{R}} $的角距离. 对$ \forall {{\boldsymbol{R}}_e} \in E $, 敏感轴为了转动到目标指向, 转角至少是 $\arccos \langle {{\boldsymbol{R}}_d}{\boldsymbol{r}}, {\boldsymbol{v}} \rangle - \left( {\alpha + \Delta \alpha } \right)$. 因此引理5成立.

    集合$ E $在引理5 中的性质如图4 所示. 对$ \forall {{\boldsymbol{R}}_e} \notin E $, 姿态的全局收敛性已经由协同势函数族$\{ P_{{q}, {{{{\boldsymbol{A}}}}}}({\boldsymbol{R_e}} ) \}_{q \in Q}$保证[23]. 为保证全局收敛, 本文讨论避障区域内势函数族${\left\{ {V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right\}_{q \in {{Q}}}}$临界点分布情况.

    图 4  引理5的图示
    Fig. 4  Illustration of lemma 5

    定理1. 令函数$ g(x) = {\sin ^2}\left( {{x}/2} \right) $. 假设目标姿态$ {\boldsymbol{R}}_d $, 敏感轴指向终点及惯性下的敏感轴禁止指向$ {\boldsymbol{v}} $已知, ${{\boldsymbol{A}}_1}、 {{\boldsymbol{A}}_3}、 \xi$如第1.2 节所定义, $ \varepsilon $如引理5所定义. 令${h_1} = \sin( {\alpha + \Delta \alpha } )/( \cos \alpha - \cos ( \alpha \; + \Delta \alpha ) )$, ${h_2} = {\sin ^2}( { {\varepsilon - 4k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1}} g( \varepsilon ) })} /2$. 当非负参数$k、 \Delta \alpha$和三阶实对称矩阵${\boldsymbol{A}}$满足:

    $$ \left\{ \begin{aligned} &\sqrt 2 k{\left\| {\boldsymbol{A}} \right\|_{\rm{F}}} < 1 \\ &2k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1} } < 1 \\ &2{\xi ^2} < 1 \\ &\sqrt {\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}} \right)} < \frac{2a\lambda _{\min }^{{{{{\boldsymbol{A}}}}_1}}{h_1h_2}}{{1 + 4k\sqrt {\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}} \right)} }} \end{aligned} \right. $$ (23)

    势函数族${\left\{ {V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right\}_{q \in {{Q}}}}$在避障区域内无临界点. 在该势函数族控制下, 航天器可以避开姿态障碍从任意的初始姿态机动到目标姿态.

    证明. 为简化书写, 对证明中多次出现的符号简写如下:

    $$ \left\{ \begin{aligned} &\gamma = \arccos \langle {{\boldsymbol{R}}_b^{\rm{T}}{\boldsymbol{v}},{\boldsymbol{r}}} \rangle \in \left[ {\alpha ,\alpha + \Delta \alpha } \right]\\ &p = \cos \alpha - {{\boldsymbol{r}}^{\rm{T}}}{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d} = \cos \alpha - \cos \gamma \\ &{\boldsymbol{R}}_c = {{\rm{e}}^{{{\left( { - 1} \right)}^q}k{P_{{{\boldsymbol A}}}}\left( {\boldsymbol{R}}_e \right){{\left[ {\boldsymbol{u}} \right]}^ \times }}}{\boldsymbol{R}}_e = {\boldsymbol{R}}_{{\cal{T}}}{\boldsymbol{R}}_e \end{aligned} \right. $$ (24)

    对$ \forall {{\boldsymbol{R}}_e} \in E $, 根据势函数求导定义计算, 可得:

    $$ \begin{array}{l} \nabla V_q\left( {{{\boldsymbol{R}}_e}} \right) = {P_{{{O}}}}\nabla {P_{q,{{{\boldsymbol{A}}}}}} + {P_{q,{{{\boldsymbol{A}}}}}}\nabla {P_{{{O}}}} \end{array} $$ (25)

    于是

    $$ \begin{split} \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right) = \;& {} {P_{{{O}}}}\psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{q,{{{\boldsymbol{A}}}}}}} \right) + \\ & {} {P_{q,{{{\boldsymbol{A}}}}}}\psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{{{O}}}}} \right) \end{split} $$ (26)

    式中, $\psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{q,{{{\boldsymbol{A }}}}}}} \right)$由式(17)给出, $\nabla {P_{{{O}}}}$为:

    $$ \nabla {P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right) = - \frac{{ab }}{{{p^{a+1}}}}\left[ {\frac{{\partial p}}{{\partial {{\left( {{{\boldsymbol{R}}_e}} \right)}_{ij}}}}} \right] = \frac{{ab }}{{{p^{a+1}}}}{{\boldsymbol{v}}_d}{\boldsymbol{r}}^{\rm{T}} $$ (27)

    利用性质1)和映射$ \psi $的定义, 有:

    $$ \begin{split} \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right)} \right) =\;& {} \frac{{ab }}{{{p^{a+1}}}}\psi \left( {{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}{{\boldsymbol{r}}^{\rm{T}}}} \right)= \\ & \frac{{ab }}{{{p^{a+1}}}}{\left[ {{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}{{\boldsymbol{r}}^{\rm{T}}} - {\boldsymbol{r}}{\boldsymbol{v}}_d^{\rm{T}}{{\boldsymbol{R}}_e}} \right]^ \vee } = \\ & \frac{{ab }}{{{p^{a+1}}}}{\left[ {\boldsymbol{r}} \right]^ \times }\left( {{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}} \right)\\[-15pt] \end{split} $$ (28)

    将$ {\boldsymbol{R}}_c $和式(28)代入式(26), 可得:

    $$ \begin{split} \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla V_q^k\left( {{{\boldsymbol{R}}_e}} \right)} \right) = \;& \frac{b }{{{p^a}}} \biggl[\frac{{a{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}}_c \right) {\left[ {\boldsymbol{r}} \right]^ \times }\left( {{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}} \right) }}{{ {p} }}+ \\ & \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{q,{{{\boldsymbol{A}}}}}}} \right) ]\\[-15pt] \end{split} $$ (29)

    分别令${\boldsymbol{x}} = \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{q,{{\boldsymbol A}}}}} \right)$和 ${\boldsymbol{y}} = a{P_{{{\boldsymbol A}}}}\left( {\boldsymbol{R}}_c \right) {\left[ {\boldsymbol{r}} \right]^ \times }\;\times \left( {{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}} \right)/ p$.

    1)对$ \forall {{\boldsymbol{R}}_e} \in E $, 证明$ \left\| {\boldsymbol{y}} \right\| $存在下限.

    向量$ {\boldsymbol{y}} $的2范数为:

    $$ \left\| {\boldsymbol{y}} \right\| = \frac{{a{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}}_c \right)\sin \gamma }}{{{\cos \alpha - \cos \gamma } }} $$ (30)

    由于$ {\boldsymbol{R}}_{{\cal{T}}}^{\rm{T}}{\boldsymbol{R}}_c = {\boldsymbol{R}}_e $, $ \Phi_{{\boldsymbol{R}}_{{\cal{T}}}^{\rm{T}}} = \Phi_{{\boldsymbol{R}}_{{\cal{T}}}} $, 利用度量的三角不等式和平移不变性, 有:

    $$ \begin{array}{l} \Phi_{{\boldsymbol{R}}_e} \le \Phi \left( {{\boldsymbol{R}}_{{\cal{T}}}^{\rm{T}}{\boldsymbol{R}}_c,{\boldsymbol{R}}_c} \right) + \Phi_{{\boldsymbol{R}}_c} \le \Phi_{{\boldsymbol{R}}_{{\cal{T}}}} + \Phi_{{\boldsymbol{R}}_c} \end{array} $$ (31)

    从旋转矩阵的指数映射可知, $\Phi_{{\boldsymbol{R}}_{{\cal{T}}}} = k{P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}}_e \right)$. 当$2k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1} } < 1$时, 引理4成立, 于是将式(31)移项, 可得:

    $$ \begin{array}{l} \Phi_{{\boldsymbol{R}}_c} \ge \Phi_{{\boldsymbol{R}}_e} - k{P_{{{{\boldsymbol{A}}}}}}\left( {{{\boldsymbol{R}}_e}} \right) > 0 \end{array} $$ (32)

    利用性质3), 有:

    $$ \begin{array}{l} \Phi_{{\boldsymbol{R}}_c} \ge \Phi_{{\boldsymbol{R}}_e} - 4k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1}} g(\Phi_{{\boldsymbol{R}}_e}) \end{array} $$ (33)

    根据引理5, 对$ \forall {{\boldsymbol{R}}_e} \in E $, $\Phi_{{{\boldsymbol R}}_e} \in \left[ {\varepsilon ,\pi } \right]$. 利用引理3, 当参数$ k $和$ {\boldsymbol{A}} $满足$2k\lambda _{\max }^{{{\boldsymbol A}}_1 } < 1$时:

    $$ \begin{array}{l} \Phi_{{{\boldsymbol R}}_c} \ge \varepsilon - 4k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_1}}g\left( \varepsilon \right) = {\Phi _{\min }}\left( {{\boldsymbol{R}}_c,{\boldsymbol{I}}} \right) \end{array} $$ (34)

    因此, 对$ \forall {{\boldsymbol{R}}_e} \in E $, 势函数${P_{{{\boldsymbol A}}}}\left( {\boldsymbol{R}}_c \right)$存在最小值, 即:

    $$ {P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}}_c \right) {} \ge 4\lambda _{\min }^{{{{{\boldsymbol{A}}}}_1}}g(\Phi_{{\boldsymbol{R}}_c}) \ge 4\lambda _{\min }^{{{{{\boldsymbol{A}}}}_1}}h_2 > 0 $$ (35)

    利用引理1, 可知:

    $$ \begin{array}{l} \left\| {\boldsymbol{y}} \right\| \ge {\left\| {\boldsymbol{y}} \right\|_{\min }} = 4a\lambda _{\min }^{{{{\boldsymbol A}}_1}}{h_1}{h_2} \end{array} $$ (36)

    2)对$ \forall {{\boldsymbol{R}}_e} \in E $, 证明$ \left\| {\boldsymbol{x}} \right\| $存在上限.

    将式(17)代入$ {\boldsymbol{x}} $, 利用$ \psi \left( { - {\boldsymbol{B}}} \right) = \psi \left( {{{\boldsymbol{B}}^{\rm{T}}}} \right) $, 有:

    $$ \begin{split} {\boldsymbol{x}} =\;& {}{\left[ {{\boldsymbol{I}} + 2k{{\boldsymbol{R}}_e^{\rm{T}}}{\boldsymbol{u}} \psi{\left( - {{\boldsymbol{A}}}{{\boldsymbol{R}}_e^{\rm{T}}} \right)}^{\rm{T}} {{{\boldsymbol{R}}_e}}} \right]^{\rm{T}}}\psi{ \left(- {{\boldsymbol{R}}_c^{\rm{T}}}{{\boldsymbol{A}}}\right)} = \\ & \psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}}) + 2k{{\boldsymbol{R}}_e^{\rm{T}}}\psi{\left( - {{\boldsymbol{A}}}{{\boldsymbol{R}}_e^{\rm{T}}}\right)} {{\boldsymbol{u}}^{\rm{T}}}{{{\boldsymbol{R}}_e}}\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}})\\[-10pt] \end{split} $$ (37)

    利用范数三角不等式, 有:

    $$ \begin{split} \left\| {\boldsymbol{x}} \right\| \le\;& {} 2k\left\| {{{\boldsymbol{R}}_e^{\rm{T}}}\psi{ \left(- {{\boldsymbol{A}}}{{\boldsymbol{R}}_e^{\rm{T}}}\right)} {{\boldsymbol{u}}^{\rm{T}}}[{{{\boldsymbol{R}}_e}}\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}})]} \right\| \;+ \\ &\left\| {\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}}) } \right\| \end{split} $$ (38)

    式中, ${{\boldsymbol{u}}^{\rm{T}}}\left[ {{{{\boldsymbol{R}}_e}}\psi \left( {{{{}}}{{{\boldsymbol{AR}}_c}}} \right)} \right] \in {{\bf{R}}}$, 并且有:

    $$ \begin{split} {{\boldsymbol{u}}^{\rm{T}}}\left[ {{{{\boldsymbol{R}}_e}}\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}})} \right] =\;& \left\| {{{{\boldsymbol{R}}_e}}\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}})} \right\|\langle {{\boldsymbol{u}}, {{{{\boldsymbol{R}}_e}}\psi( {{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}})} }} \rangle \; \le\\ \;& \left\| {\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}})} \right\|\\[-10pt] \end{split} $$ (39)

    式(38)可写为:

    $$ \begin{array}{l} \left\| {\boldsymbol{x}} \right\| \le \left\| {\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}})} \right\| + 2k\left\| {\psi({{{\boldsymbol{A}}}{{\boldsymbol{R}}_e^{\rm{T}}}})} \right\|\left\| {\psi ( {{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}} )} \right\| \end{array} $$ (40)

    利用性质4), 有:

    $$ \begin{array}{l} \left\| {\psi ( {{{{{\boldsymbol{A}}}}}{{{\boldsymbol{R}}_c}}}) } \right\| \le \sqrt {4\lambda _{\max }^{{{{\boldsymbol A}}_3} }\left( {1 - {\xi ^2}g\left( {{\Phi _{{{{{\boldsymbol R}}_c}}}}} \right)} \right)g\left( {{\Phi _{{{{{\boldsymbol R}}_c}}}}} \right)} \end{array} $$ (41)
    $$ \begin{array}{l} \left\| {\psi ( {{{\boldsymbol{A}}}{{\boldsymbol{R}}_e^{\rm{T}}}} )} \right\| \le \sqrt {4\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}g\left( {{\Phi _{{{{{\boldsymbol R}}_e}}}}} \right)} \right)g\left( {{\Phi _{{{{{\boldsymbol R}}_e}}}}} \right)} \end{array} $$ (42)

    对$ \forall {{\boldsymbol{R}}_e} \in E $, $\Phi_{{{\boldsymbol R}}_e} \in \left[ {\varepsilon ,\pi } \right]$, $\Phi_{{{\boldsymbol R}}_c} \in \left[ {{\Phi _{\min }}\left( {{{\boldsymbol R}}_c,{\boldsymbol{I}}} \right) ,\pi } \right]$, 当$ 2{\xi ^2} < 1 $时, 利用引理2, 有:

    $$ \begin{array}{l} \left\| {\psi({{{\boldsymbol{A}}}{{{\boldsymbol{R}}_c}}}) } \right\| \le \sqrt {4\lambda _{\max }^{{{{\boldsymbol A}}_3} }\left( {1 - {\xi ^2}} \right)} \end{array} $$ (43)
    $$ \begin{array}{l} \left\| {\psi ({{{\boldsymbol{A}}}{{\boldsymbol{R}}_e^{\rm{T}}}}) } \right\| \le \sqrt {4\lambda _{\max }^{{{\boldsymbol A}}_3 }\left( {1 - {\xi ^2}} \right)} \end{array} $$ (44)

    因此, 可得$ \left\| {\boldsymbol{x}} \right\| $的上限为:

    $$ \begin{split} \left\| {\boldsymbol{x}} \right\|\le \;& \sqrt {4\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}} \right)}\; + \\ &8k\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}} \right) = {\left\| {\boldsymbol{x}} \right\|_{\max }} \end{split} $$ (45)

    3)讨论$ E $内势函数的临界点. 若:

    $$ \sqrt {\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}} \right)} < \frac{2a\lambda _{\min }^{{{{{\boldsymbol{A}}}}_1}}{h_1h_2}}{{1 + 4k\sqrt {\lambda _{\max }^{{{{{\boldsymbol{A}}}}_3} }\left( {1 - {\xi ^2}} \right)} }} $$ (46)

    则有:

    $$ \begin{array}{l} {\left\| {\boldsymbol{x}} \right\|_{\max }} < {\left\| {\boldsymbol{y}} \right\|_{\min }} \end{array} $$ (47)

    对于$ \forall {{\boldsymbol{R}}_e} \in E $, 有:

    $$ \begin{array}{l} \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right) \ne {\boldsymbol{0}} \end{array} $$ (48)

    首先, 势函数族${\left\{ {V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right\}_{q \in {{Q}}}}$在避障区域内无临界点; 然后, 在姿态机动过程中, 当惯性系下敏感轴指向$ {\boldsymbol{R}}_b{{\boldsymbol{r}}} $和禁止指向$ {\boldsymbol{v}} $的夹角趋于$ \alpha $时, ${P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right)$趋于无穷大, 因此姿态机动过程可以避开障碍; 最后, 任给初始姿态$ {{\boldsymbol{R}}_e}\left( 0 \right) $, 若该姿态是非临界点, 则在势函数的控制下向单位元$ {\boldsymbol{I}} $收敛. 若该姿态是临界点, 则$ {{\boldsymbol{R}}_e}\left( 0 \right) \notin E $, 在协同性要求下, 利用混合控制可实现向单位元$ {\boldsymbol{I}} $收敛[23]. 因此, 航天器可以避开姿态障碍, 从任意的初始姿态机动到目标姿态.

    在基于势函数的受限制姿态机动问题中, 除分数形式排斥势函数(21)外, 还有对数形式的排斥势函数[16, 24]:

    $$ {P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right) = {} a - b \ln {\frac{{\cos \alpha - {{\boldsymbol{r}}^{\rm{T}}}{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}}}{c}} $$ (49)

    式中, $a、b、c$是常数. 在第2.2节中, 本文得出在考虑姿态受限情形下姿态全局收敛的参数设计要求(23). 本节将证明对数形式排斥势函数也可以经过相似推导得出参数设计要求, 以此说明本文方法的可拓展性.

    对于式(49), 利用势函数导数相关定义计算, 可得:

    $$ \begin{split} \nabla {P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right) = \;& \left[ {\frac{\partial }{{\partial {{\left( {{{\boldsymbol{R}}_e}} \right)}_{ij}}}}\left( {a - b \ln \frac{p}{c}} \right)} \right] = \\ & {}\frac{{ - b }}{p}\left[ {\frac{\partial }{{\partial {{\left( {{{\boldsymbol{R}}_e}} \right)}_{ij}}}}\left( {\cos \alpha - {{\boldsymbol{r}}^{\rm{T}}}{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}} \right)} \right] = \\ &\frac{b }{p}{{\boldsymbol{v}}_d}{\boldsymbol{r}}^{\rm{T}} \\[-15pt]\end{split} $$ (50)

    因此, 有:

    $$ \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{{{O}}}}\left( {{{\boldsymbol{R}}_e}} \right)} \right) = {} \frac{{b }}{{{p}}}\psi \left( {{\boldsymbol{R}}_e^{\rm{T}}{{\boldsymbol{v}}_d}{{\boldsymbol{r}}^{\rm{T}}}} \right) = \frac{{b }}{{{p}}} \sin \gamma $$ (51)

    将式(51)代入式(25), 可得:

    $$ \begin{split} \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla V_q} \right) =\;& {P_{{{O}}}\left( {{{\boldsymbol{R}}_e}} \right)} \biggl[ \psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla {P_{q,{{\boldsymbol{A}}}}}} \right) +\\ & {P_{{\boldsymbol{A}}}}\left( {{{\boldsymbol{R}}_c}} \right)\frac{{b \sin \gamma}}{{p{P_{{{O}}}\left( {{{\boldsymbol{R}}_e}} \right)}}}] \end{split} $$ (52)

    令${\boldsymbol{y}} = b{P_{{{{\boldsymbol{A}}}}}}\left( {{{\boldsymbol{R}}_c}} \right)\sin \gamma/ \left( {p{P_{{{O}}}\left( {{{\boldsymbol{R}}_e}} \right)}}\right)$, 则:

    $$ \left\| {\boldsymbol{y}} \right\| = \frac{{b {P_{{{{\boldsymbol{A}}}}}}\left( {\boldsymbol{R}}_c \right)\sin \gamma }}{{\left( {\cos \alpha - \cos \gamma } \right)P_{{{O}}}\left( {{{\boldsymbol{R}}_e}} \right)}} $$ (53)

    取函数$ f(\gamma) $为:

    $$ f(\gamma) = \frac{{\sin \gamma }}{{\left( {\cos \alpha - \cos \gamma } \right)P_{{{O}}}\left( {{{\boldsymbol{R}}_e}} \right)}} $$ (54)

    由于

    $$ \frac{{{\rm{d}}f\left( \gamma \right)}}{{{\rm{d}}\gamma }} = \frac{{\left( {\cos \gamma \cos \alpha - 1} \right){P_{{{O}}}} + b {{\sin }^2}\gamma }}{{{{\left( {\left( {\cos \alpha - \cos \gamma } \right){P_{{{O}}}}} \right)}^2}}} $$ (55)

    注意到, 当$P_O \ge a \; > \;0, \;\cos ( {\alpha + \Delta \alpha } )\cos \alpha \; -$ $1 \le \cos \gamma \cos \alpha - 1 \le {\cos ^2}\alpha - 1 < 0, \;b {\sin ^2}\alpha \le b {\sin ^2}\gamma \;\le$ $ b {\sin ^2}( {\alpha + \Delta \alpha } ) $时, 存在合适的参数ab, 使得$\frac{ {\rm{d}}f( \gamma )}{ {\rm{d}}\gamma } \le 0$, 即$ f(\gamma) $在区间$ [ {\alpha ,\alpha + \Delta \alpha } ] $上递减, 则:

    $$ \begin{array}{l} f\left( \gamma \right) \ge {f_{\min }} = f\left( {\alpha + \Delta \alpha} \right) \end{array} $$ (56)

    因此, 式(54)有最小值, 利用和定理1类似的推理方式, 可得对数形式排斥势函数的参数设计要求.

    文献[10, 23]利用“角度扰动”方法构建协同势函数族, 将其应用在无约束的姿态机动控制中, 实现了姿态全局收敛. 本文将这种构造方法应用到受限姿态机动任务的势函数设计中, 实现了带有禁止指向约束的姿态全局收敛. 与其他处理受限姿态机动问题的势函数方法[1516, 24]相比, 本文方法的优点是不依赖势函数凸性假设, 可避开非期望的临界点实现全局收敛; 缺点是势函数切换会造成控制律突变, 不能处理强制指向约束. 第3节将通过仿真分析具体讨论本文方法的优缺点.

    本节将建立航天器姿态受限情形下的混合控制模型, 利用第2.2节提出的协同势函数族${\left\{ {V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right\}_{q \in {{Q}}}}$进行控制, 来验证方法有效性.

    基于势函数族${\left\{ {V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right\}_{q \in {{Q}}}}$, 定义下列两个函数:

    $$ \begin{array}{l} \rho \left( {{\boldsymbol{R}}_e} \right) = \mathop {\min }\limits_{q \in {{Q}}} {V_q}\left( {{{\boldsymbol{R}}_e}} \right) \end{array} $$ (57)
    $$ \begin{array}{l} \sigma \left( {{\boldsymbol{R}}_e} \right) = \left\{ {\left. {q \in {{Q}}} \right|{V_q}\left( {{{\boldsymbol{R}}_e}} \right) = \rho \left( {{{\boldsymbol{R}}_e}} \right)} \right\} \end{array} $$ (58)

    在混合控制中, 系统状态$ \left( {{\boldsymbol{R}}_e,{{\boldsymbol{\omega}} _e},q} \right) $被分成$M, N \in SO\left( 3 \right) \times {{\bf{R}}^{3}} \times {{Q}}$两个集合:

    $$ \begin{array}{l} M = \left\{ {\left. {\left( {{\boldsymbol{R}}_e,{{\boldsymbol{\omega}} _e},q} \right)} \right|{V_q}\left( {{{\boldsymbol{R}}_e}} \right) - \rho \left( {{{\boldsymbol{R}}_e}} \right) \le \delta } \right\} \end{array} $$ (59)
    $$ \begin{array}{l} N = \left\{ {\left. {\left( {{\boldsymbol{R}}_e,{{\boldsymbol{\omega}} _e},q} \right)} \right|{V_q}\left( {{{\boldsymbol{R}}_e}} \right) - \rho \left( {{{\boldsymbol{R}}_e}} \right) > \delta } \right\} \end{array} $$ (60)

    式中, $ \delta $为式(14)给定的协同间隙. 令$x = ( {\boldsymbol{R}}_e, {{\boldsymbol{\omega}} _e})$, 根据不同的系统状态, 当$ \left( {x,q} \right) \in M $时, 式(7)、式(8)航天器模型可写为连续模型:

    $$ \left\{ \begin{aligned} &\dot x = f\left( x \right) + {\boldsymbol{\tau}}\\ &\dot q = 0 \end{aligned} \right. $$ (61)

    当$\left( {x,q} \right) \in {{N}}$时, 可写为切换模型:

    $$ \left\{ \begin{aligned} &{x^ + } = x \\ &{q^ + } = \sigma \left( {{\boldsymbol{R}}_e} \right) \end{aligned} \right. $$ (62)

    即系统状态位于集合$ M $中时, 势函数指标不变, 航天器的状态遵循微分方程式(7)、式(8)连续变化. 文献[23] 指出, 协同势函数族的所有非单位元的临界点都落在集合$ N $中, 此时航天器状态不改变, 势函数指标发生切换, 切换到非临界且值最低的势函数进行控制. 设计比例−微分反馈控制律如下:

    $$ \begin{array}{l} {\boldsymbol{\tau}} = {\boldsymbol{F}} - {k_p}\psi \left( {{\boldsymbol{R}}_e^{\rm{T}}\nabla V_q\left( {{{\boldsymbol{R}}_e}} \right)} \right) - {k_d}{{\boldsymbol{\omega}} _e} \end{array} $$ (63)

    在该控制律下的稳定性分析见文献[23].

    为了验证避障和切换控制, 本节进行了3种仿真案例分析. 仿真中, 姿态限制假定为3个, 具体信息如表1 所示. 针对这3个姿态限制, 势函数中的参数$ \Delta \alpha $分别为10°、5°和8°. 假设航天器的转动惯量矩阵$ {\boldsymbol{J}} = {\rm{diag}}\left\{ 4,5,4.5 \right\}{\rm{kg}} \cdot {{\rm{m}}}^2 $. 根据式(23)要求, 势函数参数设计为${\boldsymbol{w}} =[ 0.3,0.4, 0.6]^{\rm{T}}$, ${\boldsymbol{u}} = {\boldsymbol{w}}/ \left\| {\boldsymbol{w}} \right\|,$ ${\boldsymbol{A}} = {\rm{diag}}\left\{ {0.3,0.4,0.6} \right\}, $ $k=0.25, a=0.7, $ $b = 1/20$. 在上述参数下, 对于每一个姿态约束, 不等式组(23)都能得到满足. 根据式(13)、函数(4), 在上述参数下的奇异点为 ${\rm{Crit}}{P_{{{\boldsymbol A}}}} = \left\{ {\boldsymbol{I}} \right\} \cup {{\boldsymbol{R}}}\left( {\pi ,{\boldsymbol{e}}_i} \right), i = 1,2,3$, 其中$ {\boldsymbol{e}}_i $是3维空间$ {{\bf{R}}^{3}} $的标准正交基. 协同势函数族$\left\{ {P_{q,{{{\boldsymbol A}}}}}\left({\boldsymbol{R}}_e \right) \right\}_{q \in {{Q}}}$的 6 个非单位元临界点可根据式(17)确定, 当姿态误差处于这些临界点时, 指向$ {\boldsymbol{r}} $在惯性系下的指向以三角符号表示, 见图5 ~ 7. 可以看出, 这些指向都落在集合$ E $之外.

    表 1  惯性系下姿态限制
    Table 1  Attitude constraints in the inertial frame
    序号 指向${\boldsymbol{v}}$ 角度$\alpha \ (^ \circ) $
    CZ1 $\left[ {0.5237,0.7208,0.4540} \right]^{\rm{T}}$ 20
    CZ2 $\left[ {-0.5530,0.7612,-0.3387} \right]^{\rm{T}}$ 15
    CZ3 $\left[ {-0.1488,-0.9393,0.3090} \right]^{\rm{T}}$ 25
    下载: 导出CSV 
    | 显示表格
    图 5  无切换下姿态机动仿真
    Fig. 5  Attitude maneuver simulation without switching
    图 6  避障切换下姿态机动仿真
    Fig. 6  Attitude maneuver simulation when avoiding attitude constraints
    图 7  临界点切换下姿态机动仿真
    Fig. 7  Attitude maneuver simulation when avoiding critical points

    控制器参数设置如下: 在案例1中, 设置$ \delta = 0.06 $, $ {k_p} = 1 $, $ {k_d} = 5 $, 仿真时间为$ 60\;{\rm{s}} $; 在案例2和案例3中, 设置$ \delta = 0.06 $, $ {k_p} = 1 $, $ {k_d} = 3 $, 仿真时间为40 s. 在3个案例中, 都假定航天器存在力矩饱和, 力矩饱和值为0.5 N·m.

    案例1. 无任何切换

    假设航天器的目标姿态、初始姿态误差和航天器在本体下固连的指向分别为:

    $$ \begin{split} &{\boldsymbol{R}}_d = \exp \left( {{{[ {{{[ {-1.4002,-2.0192,1.8282} ]}^{\rm{T}}}} ]}^ \times }} \right)\\ &{\boldsymbol{R}}_e(0) = \exp \left( {{{[ {{{[ {-2.8782,-0.5547,0.3698} ]}^{\rm{T}}}} ]}^ \times }} \right)\\ &{\boldsymbol{r}} = [ 0.5344,-0.6835,0.4973]^{\rm{T}} \end{split} $$

    案例1使用上述初始条件进行2组仿真实验, 航天器的初始指向位于姿态约束CZ2下方, 目标姿态位于姿态约束CZ2和 CZ3之间. 航天器初始姿态误差$ {\boldsymbol{R}}_e(0) $位于势函数${P_{1,{{{\boldsymbol A}}}}}\left({\boldsymbol{R}}_e \right)$的临界点邻域内, 2组仿真没有进行避障切换, 分别在势函数${P_{1,{{{\boldsymbol A}}}}}\left({\boldsymbol{R}}_e \right)$和${P_{2,{{{{\boldsymbol{A}}}}}}}\left({\boldsymbol{R}}_e \right)$控制下, 绕过CZ3运动至目标姿态. 在2组仿真实验中, 向量$ {\boldsymbol{r}} $的指向轨迹如图5(a)所示, 势函数值变化曲线如图5(b)所示. 第1组仿真实验从临界点邻域开始机动, 机动路径更长, 在$ t < 20\;{\rm{s}} $时, 势函数基本保持不变; 在$t = 50\;{\rm{s}}$时, 完成收敛. 相比之下, 第2组仿真实验从非临界点邻域开始姿态机动, 机动路径更短, 势函数收敛更快, 在$ t = 30\;{\rm{s}} $时, 完成收敛. 2组仿真实验的控制力矩和角速度误差变化曲线如图5(c)和5(d)所示. 由于第1组仿真从临界点邻域开始, 导致收敛缓慢, 因此图5(c)的峰值出现时间晚于图5(d). 案例1用于说明临界点问题的极端情况, 当姿态靠近势函数的临界点时, 收敛变慢. 根据式(10), 若恰好处于临界点, 则会完全停止收敛.

    案例2. 避障切换

    在案例2中, 本文进行了2组仿真实验. 在第1组仿真实验中, 假设航天器的目标姿态、初始姿态误差和航天器在本体下固连的指向分别为:

    $$ \begin{split} &{\boldsymbol{R}}_d = \exp \left( {{{[ {{{[ {-0.3012,-2.1271,-2.1086} ]}^{\rm{T}}}} ]}^ \times }} \right) \\ &{\boldsymbol{R}}_e(0) = \exp \left( {{{[ {{{[ {-0.4906,-1.9914,-1.0410} ]}^{\rm{T}}}} ]}^ \times }} \right)\\ &{\boldsymbol{r}} = [ 0.9753,-0.2156,-0.0472]^{\rm{T}} \end{split} $$

    在第2组仿真实验中, 假设航天器的目标姿态、初始姿态误差和航天器在本体下固连的指向分别为:

    $$ \begin{split} &{\boldsymbol{R}}_d = \exp \left( {{{[ {{{[ {-0.2456,-0.6270,0.5070} ]}^{\rm{T}}}} ]}^ \times }} \right) \\ &{\boldsymbol{R}}_e(0) = \exp \left( {{{[ {{[ {-1.3113,0.6093,-1.7902} ]}^{\rm{T}}}} ]}^ \times } \right)\\ &{\boldsymbol{r}} = \left[ -0.6849,0.3701,0.6276\right]^{\rm{T}} \end{split} $$

    在2组仿真实验中, 目标姿态下向量$ {\boldsymbol{r}} $的指向相同, 位于CZ2和CZ3之间; 而位于CZ1下方的起始指向不同. 在第1组仿真初值下, 向量$ {\boldsymbol{r}} $的指向轨迹如图6(a)中虚线所示, 轨迹由CZ2下方绕过CZ2, 到达目标点. 在绕过CZ2时, 存在1次转折, 这是由于势函数在$t = 6\; {\rm{s}}$时, 进行了避障切换, 规避姿态约束CZ2. 势函数随时间变化趋势如图6(b)中虚线所示, 在$ t = 20 $s后, 收敛到0, 姿态机动完成. 在第2组实验条件下, 向量$ {\boldsymbol{r}} $的指向轨迹如图6(a)中实线所示, 轨迹穿过CZ1和CZ2, 到达目标点, 分别在靠近CZ1和CZ2时有一次转折. 这是因为势函数在$ t = 1.6 $s时进行避障切换, 规避CZ1; 在$ t = 7 $s时进行避障切换, 规避 CZ2. 势函数变化趋势如图6(b)中实线所示, 在$ t = 20 $s后, 收敛到0, 姿态机动完成. 2组仿真的控制力矩和角速度误差变化曲线如图6(c)和6(d)所示, 每次势函数切换都导致了控制力矩和角速度突变, 控制力矩在切换时饱和. 案例2说明本文设计的势函数在规避姿态约束上的有效性, 以及切换控制对航天器状态的影响.

    案例3. 临界点切换

    假设航天器的目标姿态、初始姿态误差和航天器在本体下固连的指向分别为:

    $$ \begin{split} &{\boldsymbol{R}}_d = \exp \left( {{{[ {{{[ {1.6253,-1.3826,0.2402} ]}^{\rm{T}}}} ]}^ \times }} \right) \\ &{\boldsymbol{R}}_e(0) = \exp \left( {{{[ {{{[ {-1.2090,-1.8832,2.0448} ]}^{\rm{T}}}} ]}^ \times }} \right)\\ &{\boldsymbol{r}} = [ -0.2790,0.7829,0.5561]^{\rm{T}} \end{split} $$

    图7 所示, 在上述初始条件下, 向量$ {\boldsymbol{r}} $的初始指向位于CZ1左下方, 目标指向位于CZ2 和CZ3 之间. 该指向从 CZ2下方绕过CZ2, 到达目标指向. 航天器在12 s前, 基于势函数${P_{2,{{\boldsymbol A}}}}$进行姿态控制; 在$ t = 12 $s时, 姿态靠近临界点, 满足切换控制模型式(61)和式(62), 势函数发生切换, 切换到值更低的势函数${P_{1,{{\boldsymbol A}}}}$进行姿态控制. 如图7(c)所示, 控制力矩和角速度误差在势函数切换时突变. 在$ t = 25 $s后, 势函数收敛到0, 姿态机动完成. 案例3说明本文方法在规避非期望临界点时的一般情况, 如图2所示, 在进入非期望临界点某一给定的邻域时, 就通过切换势函数避开该邻域.

    综上所述, 1)由上述仿真实验可知, 为了规避非期望的临界点和姿态限制, 势函数的切换会导致控制力矩发生突变, 该突变值上限受势函数参数的影响, 在势函数参数选取时, 应在满足式(23)基础上, 兼顾力矩突变大小. 2)控制器的参数选取应兼顾力矩和机动时间与最终指向误差等收敛效果的需求. 3)势函数非期望临界点的分布是其固有性质, 仅与势函数参数有关, 与动力学或运动学的不确定性无关. 但不确定性的存在会影响航天器对本文算法所给出的期望控制力矩的执行, 进而影响航天器的实际机动路径, 此时无法保证航天器能够规避姿态约束和非期望临界点.

    在姿态受约束情况下, 基于旋转矩阵设计可规避临界点并实现姿态全局收敛的协同势函数. 为每一个姿态约束设计软约束区域, 在软约束区域内, 避障势函数由相对目标姿态的吸引项和相对姿态约束的排斥项构成, 当航天器姿态运动到该区域内, 切换为避障势函数. 针对常见的分数形式和对数形式排斥项, 讨论了本文的避障势函数的临界点分布和势函数参数的关系, 给出了避障区域内不存在临界点的参数要求. 仿真实验结果表明, 航天器能够在本文设计的势函数控制下规避姿态约束, 收敛到目标姿态. 由于本文进行了势函数切换, 在切换时控制律不连续, 未来将考虑如何改善控制律突变对系统的影响.

  • 图  1  基于PASCAL VOC 2011验证集的分割结果(第1行是以马作为提取目标的案例; 第2行是以飞行器作为提取目标的案例)

    Fig.  1  Sample segmentation results on the PASCAL VOC 2011 validation set (The first row shows the ROI extraction result for horse and the second row shows the result for aircraft extraction)

    图  2  拥有不同域信息的指纹图像

    Fig.  2  Samples of 2D fingerprint images in different domains

    图  3  基于语义分割的ROI提取模型, 模型分为两部分: 基础网络和分割网络

    Fig.  3  The flowchart of ROI extraction network based on semantic segmentation

    图  4  在生物特征识别中, 基于语义分割的ROI提取模型存在的问题(第1行是人脸提取的案例分析[3536]; 第2行是指纹ROI提取的案例分析[2124])

    Fig.  4  ROI extraction issues we observe on biometrics cases (The first row shows the ROI extraction result for face[3536] and the second row presents the result for fingerprint ROI extraction[2124])

    图  5  像素级损失函数的失效情况(分割结果II是分割结果I向左平移一个像素得到, 结果显示两个分割结果的交叉熵为264.80, L2为23.00)

    Fig.  5  Failure of pixel level loss functions (Result I translates one pixel to the left to get Result II. Cross entropy between I and II is 264.80 and L2 is 23.00)

    图  6  基于对抗结构的全局损失模块的结构图

    Fig.  6  Adversarial structure based global perceptual loss module

    图  7  基于全局损失函数的ROI提取模型

    Fig.  7  Overview of our proposed ROI extraction model

    图  8  一个3D指纹的带标注横截面样例图 ((a) 标注了生物组织结构的指纹横截面图像; (b)该横截面对应的ROI区域; (c) 指尖的生物结构[39])

    Fig.  8  An example of X-Z cross-section image labeled for 3D fingerprints ((a) The cross-section image marked with biological structure; (b) The labeled image mark with the ROI; (c) Physical structure of human skin[39])

    图  9  ROI提取模型的收敛折线图(第1列的评价指标为交并比; 第2列的评价指标为像素级准确率)

    Fig.  9  The convergent plots for ROI extraction model (The evaluation metric of the first column is mean IoU, that of the second column is pixel Acc.)

    图  10  不同训练次数下的2D指纹、人脸和3D指纹ROI的提取结果(从左至右依次是不同的迭代次数的模型分割结果. 上面一行是Baseline的分割结果, 下面一行是本文方法的分割结果)

    Fig.  10  The results for 2D fingerprint, face and 3D fingerprint ROI extraction with different iteration numbers (From left to right, there are the extraction results with different iteration numbers. The upper row corresponds to the extraction results of baseline, and the lower row shows the results of the proposed method)

    图  11  人脸ROI提取和2D指纹ROI提取结果

    Fig.  11  The results for face ROI extraction and 2D fingerprint ROI extraction

    图  12  基于全局感知模块的3D指纹ROI提取结果

    Fig.  12  A set of images which show the ROI extraction result of our proposed method for 3D fingerprint

    图  13  基于全局感知模块的汗孔提取结果

    Fig.  13  The ROI extraction result of our proposed method for pore extraction

    表  1  不同设置下的全局感知模块表现

    Table  1  Investigation of global perceptual loss module with different settings

    优化策略2D传统指纹(Pixel Acc.(%)/Mean IoU)人脸(Pixel Acc.(%)/Mean IoU)3D指纹(Pixel Acc.(%)/Mean IoU)
    本文方法Baseline本文方法Baseline本文方法Baseline
    损失函数IoU loss[2]92.07/0.863290.66/0.838092.05/0.857990.03/0.825496.97/0.885995.18/0.8640
    Lovasz loss[25]92.48/0.864893.14/0.882297.21/0.947596.71/0.938895.74/0.878895.69/0.8767
    L2 loss93.33/0.861389.33/0.821996.99/0.943496.90/0.942095.70/0.885094.14/0.8331
    CrossEntropy loss (base)92.58/0.860682.71/0.718097.06/0.942996.77/0.938996.13/0.897595.43/0.8719
    优化器AMSGrad[29]93.65/0.886392.39/0.867296.50/0.935396.17/0.928993.56/0.823090.45/0.7540
    Radam[30]92.72/0.869492.27/0.866596.72/0.939096.52/0.935095.77/0.880695.19/0.8676
    Adam (base)[27]92.58/0.860682.71/0.718097.06/0.942996.77/0.938996.13/0.897595.43/0.8719
    下载: 导出CSV

    表  2  2D指纹ROI提取实验结果

    Table  2  ROI extraction results of 2D fingerprints

    方法FVCs vs. NIST 29
    pixel Acc. (%)/Mean IoU
    NIST 29 vs. FVCs
    pixel Acc. (%)/Mean IoU
    平均值(Average)
    pixel Acc. (%)/Mean IoU
    Mean and variance based method[13]76.71/0.685277.23/0.755176.97/0.7202
    Orientation based method[12]75.37/0.753274.46/0.621374.92/0.6873
    Fourier based method[13]65.45/0.634965.45/0.634965.45/0.6349
    PSPNet[17]87.74/0.800079.41/0.720983.58/0.7605
    FCN[15]87.20/0.793275.77/0.673681.49/0.7334
    U-Net[16]85.83/0.783976.46/0.725181.15/0.7545
    Baseline82.71/0.718073.12/0.734177.92/0.7261
    Baseline+Dense-CRF[48]90.33/0.783578.30/0.734784.32/0.7591
    本文方法92.58/0.860680.29/0.746986.44/0.8038
    本文方法+Dense-CRF94.67/0.885282.73/0.785288.70/0.8352
    下载: 导出CSV

    表  3  人脸提取案例实验结果

    Table  3  ROI extraction results of face images

    方法 Pixel Acc. (%) Mean IoU
    PSPNet[17] 93.62 0.8803
    FCN[15] 95.90 0.9212
    U-Net[16] 95.55 0.9147
    Baseline 96.77 0.9389
    Baseline+Dense-CRF[48] 95.77 0.9712
    本文方法 97.06 0.9429
    本文方法+Dense-CRF 96.41 0.9734
    下载: 导出CSV

    表  4  3D指纹的ROI提取结果

    Table  4  ROI extraction results of 3D fingerprint

    方法 Pixel Acc. (%) Mean IoU
    PSPNet[17] 93.67 0.8296
    FCN[15] 94.62 0.8526
    U-Net[16] 94.82 0.8614
    Baseline 95.43 0.8719
    Baseline+Dense-CRF[48] 95.50 0.8718
    本文方法 96.13 0.8975
    本文方法+Dense-CRF 96.12 0.8898
    下载: 导出CSV

    表  5  指纹汗孔提取实验结果

    Table  5  Fingerprint pore extraction results

    $ R_T$(%)$ R_F$(%)
    Gabor Filter[44]75.90 (7.5)23.00 (8.2)
    Adapt. Dog[14]80.80 (6.5)22.20 (9.0)
    DAPM[14]84.80 (4.5)17.60 (6.3)
    Xu等[45]84.80 (4.5)17.60 (6.3)
    Labati等[46]84.69 (7.81)15.31 (6.2)
    DeepPore[47]93.09 (4.63)8.64 (4.15)
    DeepPore*96.33 (6.57)6.45 (17.22)
    Baseline97.48 (9.63)7.57 (5.85)
    本文方法98.30 (9.2927)7.83 (4.18)
    下载: 导出CSV
  • [1] Caelles S, Maninis K K, Pont-Tuset J, Leal-Taixé L, Cremers D, Van Gool L. One-shot video object segmentation. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, USA: IEEE, 2017. 5320−5329
    [2] Rahman A, Wang Y. Optimizing intersection-over-union in deep neural networks for image segmentation. In: Proceedings of the 12th International Symposium on Visual Computing (ISVC). Las Vegas, USA: Springer, 2016. 234−244
    [3] Wei Y C, Feng J S, Liang X D, Cheng M M, Zhao Y, Yan S C. Object region mining with adversarial erasing: A simple classification to semantic segmentation approach. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, USA: IEEE, 2017. 6488−6496
    [4] Xu X W, Lu Q, Yang L, Hu S, Chen D, Hu Y, et al. Quantization of fully convolutional networks for accurate biomedical image segmentation. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Salt Lake City, USA: IEEE, 2018. 8300−8308
    [5] 刘青山, 卢汉清, 马颂德. 综述人脸识别中的子空间方法. 自动化学报, 2003, 29(6): 900-911

    Liu Qing-Shan, Lu Han-Qing, Ma Song-De. A survey: Subspace analysis for face recognition. Acta Automatica Sinica, 2003, 29(6): 900-911
    [6] 高全学, 潘泉, 梁彦, 张洪才, 程咏梅. 基于描述特征的人脸识别研究. 自动化学报, 2006, 32(3): 386-392

    Gao Quan-Xue, Pan Quan, Liang Yan, Zhang Hong-Cai, Cheng Yong-Mei. Face recognition based on expressive features. Acta Automatica Sinica, 2006, 32(3): 386-392
    [7] Wang S, Zhang W W, Wang Y S. New features extraction and application in fingerprint segmentation. Acta Automatica Sinica, 2003, 29(4): 622-627
    [8] Yu C Q, Wang J B, Peng C, Gao C X, Yu G, Sang N. Learning a discriminative feature network for semantic segmentation. In: Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). Salt Lake City, USA: IEEE, 2018. 1857−1866
    [9] Hsu C Y, Yang C H, Wang H C. Multi-threshold level set model for image segmentation. EURASIP Journal on Advances in Signal Processing, 2010, 2010(1): Article No. 950438 doi: 10.1155/2010/950438
    [10] Taheri S, Ong S H, Chong V F H. Level-set segmentation of brain tumors using a threshold-based speed function. Image and Vision Computing, 2010, 28(1): 26-37 doi: 10.1016/j.imavis.2009.04.005
    [11] Xu A P, Wang L J, Feng S, Qu Y X. Threshold-based level set method of image segmentation. In: Proceedings of the 3rd International Conference on Intelligent Networks and Intelligent Systems. Shenyang, China: IEEE, 2010. 703−706
    [12] Feng J J, Zhou J, Jain A K. Orientation field estimation for latent fingerprint enhancement. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2013, 35(4): 925-940 doi: 10.1109/TPAMI.2012.155
    [13] Maltoni D, Maio D, Jain A K, Prabhakar S. Handbook of Fingerprint Recognition (2nd edition). London: Springer, 2009.
    [14] Zhao Q J, Zhang D, Zhang L, Luo N. Adaptive fingerprint pore modeling and extraction. Pattern Recognition, 2010, 43(8): 2833-2844 doi: 10.1016/j.patcog.2010.02.016
    [15] Long J, Shelhamer E, Darrell T. Fully convolutional networks for semantic segmentation. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Boston, USA: IEEE, 2015. 3431−3440
    [16] Ronneberger O, Fischer P, Brox T. U-Net: Convolutional networks for biomedical image segmentation. In: Proceedings of the 18th International Conference on Medical Image Computing and Computer-Assisted Intervention (MICCAI). Munich, Germany: Springer, 2015. 234−241
    [17] Zhao H S, Shi J P, Qi X J, Wang X G, Jia J Y. Pyramid scene parsing network. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, USA: IEEE, 2017. 6230−6239
    [18] Lazebnik S, Schmid C, Ponce J. Beyond bags of features: Spatial pyramid matching for recognizing natural scene categories. In: Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR). New York, USA: IEEE, 2006. 2169−2178
    [19] Chen L C, Zhu Y K, Papandreou G, Schroff F, Adam H. Encoder-decoder with atrous separable convolution for semantic image segmentation. In: Proceedings of the 15th European Conference on Computer Vision (ECCV). Munich, Germany: Springer, 2018. 833−851
    [20] Wang P Q, Chen P F, Yuan Y, Liu D, Huang Z H, Hou X D, et al. Understanding convolution for semantic segmentation. In: Proceedings of the IEEE Winter Conference on Applications of Computer Vision (WACV). Lake Tahoe, USA: IEEE, 2018. 1451−1460
    [21] Maio D, Maltoni D, Cappelli R, Wayman J L, Jain A K. FVC2000: Fingerprint verification competition. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2002, 24(3): 402-412 doi: 10.1109/34.990140
    [22] Maio D, Maltoni D, Cappelli R, Wayman J L, Jain A K. FVC2002: Second fingerprint verification competition. In: Proceedings of the International Conference on Pattern Recognition. Quebec City, Canada: IEEE, 2002. 811−814
    [23] Maio D, Maltoni D, Cappelli R, Wayman J L, Jain A K. FVC2004: Third fingerprint verification competition. In: Proceedings of the 1st International Conference on Biometric Authentication (ICBA). Hong Kong, China: Springer, 2004. 1−7
    [24] Watson C I. NIST Special Database 29: Plain and Rolled Images from Paired Fingerprint Cards, Internal Report 6801, National Institute of Standards and Technology, Gaithersburg, USA, 2001
    [25] Berman M, Triki A R, Blaschko M B. The lovasz-softmax loss: A tractable surrogate for the optimization of the intersection-over-union measure in neural networks. In: Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition. Salt Lake City, USA: IEEE, 2018. 4413−4421
    [26] Tieleman T, Hinton G. Lecture 6.5-rmsprop: Divide the gradient by a running average of its recent magnitude. COURSERA: Neural Networks for Machine Learning, 2012, 4(2): 26-31
    [27] Kingma D P, Ba J. Adam: A method for stochastic optimization. arXiv preprint arXiv: 1412.6980, 2014.
    [28] Wilson A C, Roelofs R, Stern M, Srebro N, Recht B. The marginal value of adaptive gradient methods in machine learning. In: Proceedings of the 31st International Conference on Neural Information Processing Systems. Long Beach, USA: Curran Associates Inc., 2017. 4151−4161
    [29] Luo L C, Xiong Y H, Liu Y, Sun X. Adaptive gradient methods with dynamic bound of learning rate. In: Proceedings of the 7th International Conference on Learning Representations (ICLR). New Orleans, USA: OpenReview.net, 2019.
    [30] Liu L Y, Jiang H M, He P C, Chen W Z, Liu X D, Gao J F, et al. On the variance of the adaptive learning rate and beyond. In: Proceedings of the 8th International Conference on Learning Representations (ICLR). Addis Ababa, Ethiopia: OpenReview.net, 2020.
    [31] Goodfellow I J, Pouget-Abadie J, Mirza M, Xu B, Warde-Farley D, Ozair S, et al. Generative adversarial nets. In: Proceedings of the 27th International Conference on Neural Information Processing Systems. Montreal, Canada: MIT Press, 2014. 2672−2680
    [32] Ganin Y, Lempitsky V. Unsupervised domain adaptation by backpropagation. In: Proceedings of the 32nd International Conference on International Conference on Machine Learning. Lille, France: JMLR.org, 2015. 1180−1189
    [33] Pan S J, Yang Q. A survey on transfer learning. IEEE Transactions on Knowledge and Data Engineering, 2010, 22(10): 1345-1359 doi: 10.1109/TKDE.2009.191
    [34] He K M, Zhang X Y, Ren S Q, Sun J. Deep residual learning for image recognition. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. Las Vegas, USA: IEEE, 2016. 770−778
    [35] Kaggle. AISegment.com-Matting human datasets [Online], available: https://www.kaggle.com/laurentmih/aisegmentcom-matting-human-datasets/, March 20, 2020
    [36] Shen X Y, Tao X, Gao H Y, Zhou C, Jia J Y. Deep automatic portrait matting. In: Proceedings of the 14th European Conference on Computer Vision (ECCV). Amsterdam, the Netherlands: Springer, 2016. 92−107
    [37] Arjovsky M, Chintala S, Bottou L. Wasserstein generative adversarial networks. In: Proceedings of the 34th International Conference on Machine Learning. Sydney, NSW, Australia: PMLR, 2017. 214−223
    [38] Maas A L, Hannun A Y, Ng A Y. Rectifier nonlinearities improve neural network acoustic models. In: Proceedings of the 30th International Conference on Machine Learning (ICML). Atlanta, USA: JMLR, 2013.
    [39] Madhero88. Layers of the skin [Online], available: https://en.wikipedia.org/wiki/File:Skin_layers.png, March 20, 2020
    [40] Liu F, Shen L L, Liu H Z, Shen C X, Liu G J, Liu Y H, et al. A-benchmark-database-using-optical-coherence-tomography-for-fingerprints [Online], available: https://github.com/CV-SZU/A-Benchmark-Database-using-Optical-Coherence-Tomography-for-Fingerprints, March 20, 2020
    [41] Liu F, Shen C X, Liu H Z, Liu G J, Liu Y H, Guo Z H, et al. A flexible touch-based fingerprint acquisition device and a benchmark database using optical coherence tomography. IEEE Transactions on Instrumentation and Measurement, 2020, 69(9): 6518-6529 doi: 10.1109/TIM.2020.2967513
    [42] Liu F, Liu G J, Wang X Z. High-accurate and robust fingerprint anti-spoofing system using optical coherence tomography. Expert Systems With Applications, 2019, 130: 31-44 doi: 10.1016/j.eswa.2019.03.053
    [43] Liu H Z, Zhang W T, Liu F, Qi Y. 3D fingerprint gender classification using deep learning. In: Proceedings of the 14th Chinese Conference on Biometric Recognition (CCBR). Zhuzhou, China: Springer, 2019. 37−45
    [44] Jain A, Chen Y, Demirkus M. Pores and ridges: Fingerprint matching using level 3 features. In: Proceedings of the 18th International Conference on Pattern Recognition (ICPR). Hong Kong, China: IEEE, 2006. 477−480
    [45] Xu Y R, Lu G M, Liu F, Li Y X. Fingerprint pore extraction based on multi-scale morphology. In: Proceedings of the 12th Chinese Conference on Biometric Recognition (CCBR). Shenzhen, China: Springer, 2017. 288−295
    [46] Labati R D, Genovese A, Muñoz E, Piuri V, Scotti F. A novel pore extraction method for heterogeneous fingerprint images using Convolutional Neural Networks. Pattern Recognition Letters, 2018, 113: 58-66 doi: 10.1016/j.patrec.2017.04.001
    [47] Jang H U, Kim D, Mun S M, Choi S, Lee H K. DeepPore: Fingerprint pore extraction using deep convolutional neural networks. IEEE Signal Processing Letters, 2017, 24(12): 1808-1812 doi: 10.1109/LSP.2017.2761454
    [48] Chen L C, Papandreou G, Kokkinos I, Murphy K, Yuille A L. Semantic image segmentation with deep convolutional nets and fully connected CRFs. arXiv preprint arXiv: 1412.7062, 2014.
    [49] Chen L C, Papandreou G, Kokkinos I, Murphy K, Yuille A L. DeepLab: Semantic image segmentation with deep convolutional nets, atrous convolution, and fully connected CRFs. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2018, 40(4): 834-848 doi: 10.1109/TPAMI.2017.2699184
  • 加载中
图(13) / 表(5)
计量
  • 文章访问数:  1594
  • HTML全文浏览量:  421
  • PDF下载量:  191
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-24
  • 录用日期:  2020-06-23
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2023-06-20

目录

/

返回文章
返回