2.845

2023影响因子

(CJCR)

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

留言板

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

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

前交叉韧带断裂后足底压力特征的聚类分析

李晓理 黄红拾 王杰 于媛媛 敖英芳

王晓红, 赵于前, 廖苗, 邹北骥. 基于多尺度2D Gabor小波的视网膜血管自动分割. 自动化学报, 2015, 41(5): 970-980. doi: 10.16383/j.aas.2015.c140185
引用本文: 李晓理, 黄红拾, 王杰, 于媛媛, 敖英芳. 前交叉韧带断裂后足底压力特征的聚类分析. 自动化学报, 2017, 43(3): 418-429. doi: 10.16383/j.aas.2017.c160197
WANG Xiao-Hong, ZHAO Yu-Qian, LIAO Miao, ZOU Bei-Ji. Automatic Segmentation for Retinal Vessel Based on Multi-scale 2D Gabor Wavelet. ACTA AUTOMATICA SINICA, 2015, 41(5): 970-980. doi: 10.16383/j.aas.2015.c140185
Citation: LI Xiao-Li, HUANG Hong-Shi, WANG Jie, YU Yuan-Yuan, AO Ying-Fang. Cluster Analysis of Plantar Pressure Characteristics after Anterior Cruciate Ligament Deficiency. ACTA AUTOMATICA SINICA, 2017, 43(3): 418-429. doi: 10.16383/j.aas.2017.c160197

前交叉韧带断裂后足底压力特征的聚类分析

doi: 10.16383/j.aas.2017.c160197
基金项目: 

国家自然科学基金 61473034

低温重点实验室开放基金 CRYO 201316

国家自然科学基金 61673053

北京大学医-信交叉建设孵化基金 BMU2016-12

北京市科技新星交叉学科项目 Z161100004916041

高等学校博士学科点专项科研基金 2013000611008

低温重点实验室开放基金 TIPC

低温重点实验室开放基金 CAS

详细信息
    作者简介:

    李晓理 北京工业大学电子信息与控制工程学院教授.主要研究方向为多模型控制和自适应控制.E-mail:lixiaolibjut@bjut.edu.cn

    黄红拾 北京大学第三医院运动医学研究所副主任医师.主要研究方向为运动医学/康复, 运动能力的训练和提高, 运动损伤防护的中西医结合康复技术, 运动生物力学、康复技术和辅具功能研发.本文共同第一作者.E-mail:huanghs@bjmu.edu.cn

    于媛媛 北京大学第三医院运动医学研究所康复治疗师.主要研究方向为运动康复, 运动能力的训练和提高, 运动生物力学测试、康复技术和辅具功能研发.E-mail:yyy1091012@126.com

    通讯作者:

    王杰 北京科技大学自动化学院硕士研究生.主要研究方向为模糊控制与聚类算法.本文共同通信作者.E-mail:yizhe9288@163.com

    敖英芳 北京大学运动医学研究所教授、主任医师.主要研究方向为运动创伤学, 为膝关节运动损伤、交叉韧带与软骨损伤的修复与重建.本文通信作者.E-mail:yingfang.ao@vip.sina.com

Cluster Analysis of Plantar Pressure Characteristics after Anterior Cruciate Ligament Deficiency

Funds: 

National Natural Science Foundation of China 61473034

Opening Foundation of Key Laboratory of Cryogenics CRYO 201316

National Natural Science Foundation of China 61673053

the Incubation Grant for Medicine and Information Sciences of Peking University BMU2016-12

Beijing Nova Programme Interdisciplinary Cooperation Project Z161100004916041

Specialized Research Fund for the Doctoral Program of Higher Education 2013000611008

Opening Foundation of Key Laboratory of Cryogenics TIPC

Opening Foundation of Key Laboratory of Cryogenics CAS

More Information
    Author Bio:

    Professor at the College of Electronic Information and Control Engineering, Beijing University of Technology. His research interest covers multiple model control and adaptive control.

    Associate chief physician at Peking University Institute of Sports Medicine. His research interest covers sports medicine/rehabilitation, elite athletes sports training and performance enhancement methods, East-meets-West in rehabilitation techniques in sports injuries and prevention programs, biomechanics of sports and lower extremity injury, isokinetic measurement and training, brace and shoe functions. Co-first author of this paper.

    Therapist at Peking University Institute of Sports Medicine. Her research interest covers sports medicine/reh- abilitation, elite athletes sports training and performance enhancement methods, biomechanics of sports and lower extremity injury, isokinetic measurement and training.

    Corresponding author: WANG Jie Master student at the School of Automation and Electrical Engineering, University of Science and Technology Beijing. His research interest covers fuzzy control and clustering algorithm. Co-corresponding author of this paper.; AO Ying-Fang Professor, chief physician at Peking University Institute of Sports Medicine. His research interest covers sports injury, knee injury and cartilage repair, clinical studies and basic research on sports injuries, clinical studies on motor function following knee cruciate ligament injuries, repair and reconstruction of knee cruciate ligament. Corresponding author of this paper.
  • 摘要: 运动过程中,人体的步态特征可以在足底压力图像上有准确的记录,而这也就可以成为判断步态正常与否的一条有效依据.通过一组压力传感器阵列获取人体运动过程的足底压力分布数据,提取步态的运动学和动力学特性.在此基础上,采用极限学习机(Extreme learning machines,ELM)神经网络聚类算法对足底压力数据进行分析,完成正常与异常步态的分类辨识工作.本文从实际临床数据出发,对前交叉韧带断裂患者进行步态分析,并据医生的临床诊断结果进行校验.该方法在步态分析上取得了较为良好的效果,仿真结果表明了其有效性.
  • 工业控制系统(Industrial control system)广泛应用于国家基础设施, 工控系统一旦遭到破坏可能会造成难以估量的经济损失甚至人员伤亡.一个典型的工业控制网路可分为企业网络层、控制网络层和现场网络层, 如图 1所示.

    图 1  典型的工业控制系统网络架构
    Fig. 1  Typical network architecture for industrial control system

    近年来, 由于IT系统的软件和硬件技术不断集成到工控领域, IT领域的一些漏洞和后门也出现在工控系统上.以智能制造为核心的工业4.0, 在推动工业转型的同时也使得工控系统面临更多来自互联网的威胁[1]. 2010年, 世界首个网络超级破坏武器"震网"病毒[2]被检测出来, 证明了工控系统是可被攻击并被利用的.此后又相继爆发了"毒区"、"火焰"等病毒和一连串的工控系统入侵事件.工控信息安全受到越来越多的重视, 针对工控系统入侵检测的研究也成为目前的一个热点.

    由于工控保护机制的存在, 大的破坏性攻击容易被检测出来, 因此对入侵检测的研究越来越集中于隐蔽攻击[3-5].当前, 工控系统入侵检测的研究方法主要包含三类: 1)基于概率统计的方法[6-7]; 2)基于机器学习的方法[8-12]; 3)基于状态的方法[13-16].

    基于状态的方法因为更好地考虑了工控系统的物理特性, 检测准确度高, 目前大部分研究都是基于此类方法[13].基于状态的方法通过历史数据和专家经验定义系统的临界状态, 实时监控系统当前状态与临界状态的距离以判断系统是否处于危险之中.Fovino等[15]通过监控系统状态的变化来检测复杂攻击, 创立了系统的虚拟镜像作为系统的内部表示, 并用规则语言描述系统的临界状态, 但是仅考虑了离散输入输出数据, 没有考虑连续数据. Carcano等[17]进一步考虑了连续输入和输出, 但在输入输出数据量庞大的情况下, 很难定义系统的临界状态.针对这一问题, Khalili等[13]提出了一个系统化的解决方案SysDetect, 该方案基于Apriori算法, 通过专家经验的判定可显著减少算法下一次迭代产生的候选临界状态数.但是该方法依然不能很好地处理高维数据, 且在一定程度上依赖人工判定.

    针对当前工控入侵检测系统存在的问题, 本文提出一种新的基于状态迁移图的异常检测方法.利用工控系统与物理世界交互的特性, 采用数据驱动的方法进行建模, 即利用系统运行时系统变量的数据来建立检测模型.但考虑到数据维度可能过高, 不直接以系统变量数据来描述系统运行状态, 而是以相邻数据向量间的余弦相似度和欧氏距离来表征, 因而可以处理高维数据.状态迁移图刻画系统运行过程中的正常模型, 根据正常历史数据样本训练得出, 所以不需要依赖专家预先定义系统的临界状态.

    一个工业控制过程可简单地用图 2来表示.其中传感器和执行器易成为攻击者的攻击目标, 因为它们直接与物理过程相连, 一旦遭到攻击将会产生不可估量的后果, 例如"震网"病毒恶意篡改了伊朗核工厂离心机的转速, 使核工厂被迫停工.

    图 2  工业控制过程:传感器和执行器易成为攻击目标
    Fig. 2  Industrial control process: sensors and actuators are vulnerable targets

    由于工控系统的物理特性, 恶意数据注入成为一种简单且收益高的方式.攻击者通过攻击执行器和传感器, 注入恶意数据, 修改变量的值, 对系统造成破坏.

    系统变量分为传感变量和操纵变量, 传感变量的值由传感器测量取得, 而操纵变量则由控制器通过系统参数和传感变量值计算得出.记系统含$l$个测量变量, 用集合${\mathit{\boldsymbol{y}}} = \{ {{y_k}|k = 1, 2, \cdots , l} \}$表示, 含$p$个操纵变量, 用集合${\mathit{\boldsymbol{u}}} = \{ {{u_k}|k = 1, 2, \cdots , p} \}$表示.则系统变量集合可用${\mathit{\boldsymbol{s}}} = \{ {{s_j}|j = 1, 2, \cdots , l + p} \}$表示.根据攻击者试图破坏系统速度的快慢, 将恶意数据注入攻击分为快速注入和慢速注入.

    快速注入使得系统变量在短时间内发生较大改变, 以最快速度造成破坏, 表现为最大/最小值攻击.

    $ {\tilde s_i}\left( t \right) = \begin{cases} {s_i}\left( t \right)\text{, }{\kern 20pt} t \notin {T_a}\\ {s_i}^{\max }\text{, }{\kern 20pt} t \in {T_a} \end{cases} $

    (1)

    $ {\tilde s_i}\left( t \right) = \begin{cases} {s_i}\left( t \right)\text{, }{\kern 20pt} t \notin {T_a}\\ {s_i}^{\min }\text{, }{\kern 22pt} t \in {T_a} \end{cases} $

    (2)

    其中, ${T_a}$为攻击时段.

    慢速注入缓慢改变系统变量的值, 以达到潜伏的目的, 可以表现为偏置注入和几何注入.偏置注入连续注入多个较小的常量.

    $ {\tilde s_i}\left( t \right) = \begin{cases} {s_i}\left( t \right) + {c_i}\text{, }{\kern 20pt} t \in {T_a}\\ {s_i}\left( t \right)\text{, } t \notin {T_a} \end{cases} $

    (3)

    其中, 表示常量, 通常取值较小.几何注入逐渐增大系统变量的改变幅度, 具有一定的潜伏性, 同时不失破坏性, 可表现为指数形式的注入.

    $ {\tilde s_i}\left( t \right) = \begin{cases} {s_i}\left( t \right) + a{b^{t - {t_0}}}\text{, }{\kern 10pt} t \in {T_a}\\ {s_i}\left( t \right)\text{, }{\kern 49pt} t \notin {T_a} \end{cases} $

    (4)

    其中, $a$, $b$可以是常数, 也可以是变量, ${t_0}$为常数.

    以上只是列举了三种可能的攻击形式, 实际上恶意数据注入攻击可以表现为更多的形式.

    由于工控系统的诸多限制, 例如有限的内存、有限的计算能力、较高的实时性要求等, 工控入侵检测系统必须是轻量级的, 检测规则或检测模型应设置得相对简单.本文采用一个状态迁移图的检测模型, 状态迁移图的状态数与数据向量的维度无关, 因而可以处理高维的数据向量.

    实际的工控系统中, 状态变量可能有很多个, 因此数据向量的维度可能很高, 而现有的基于状态的入侵检测算法无法很好地处理数据维度较高的情况.考虑控制系统的运行模式, 控制器根据前一时刻传感变量的值, 计算出操纵变量的值, 操纵变量的值再作用于物理系统使之发生变化, 传感变量随着系统的运动而变化, 在下一个采样周期将传感变量的值传给控制器, 如此循环往复.通过分析这个过程可以发现, 系统变量在当前时刻的取值很大程度上取决于上一时刻的取值, 因此正常情况下, 相邻两个数据向量间的变化不会太大.

    相邻数据向量之间的变化反应了系统变量运行的时间特性, 本文以余弦相似度和欧氏距离来衡量这种变化.余弦相似度反应数据向量方向的变化, 常用于文档聚类和信息提取.欧氏距离衡量两个数据向量绝对距离的远近.

    选择作为系统状态表征的特征应该是相互关联的, 即相邻数据向量间的余弦相似度和欧氏距离应该是相互关联的.这样当余弦相似度在某个区间取值时, 欧氏距离只位于某些特定区间, 而不会在所有区间取值.如图 3所示, 假设余弦相似度取值区间为$A$, $B$, $C$, 欧氏距离取值区间为$D$, $E$, $F$, 二者在二维空间所有可能的区域为1 $\sim$ 9, 实际取值区域为1, 4, 5, 9.当余弦相似度取值区间为$B$时, 欧氏距离只在区间$D, E$取值, 而不会在$F$区间取值, 若欧氏距离在$F$区间取值就可以认为发生了异常.这样可以更加细致地刻画数据的轮廓, 而不是将所有可能的区域为1 $\sim$ 9都视为可接受的区域.余弦相似度和欧氏距离的关联性如图 4所示, 其中余弦相似度和欧氏距离采用Normal数据集计算得出(经过归一化处理), 图 4展现了前30组数据.从图 4可以看出, 余弦相似度和欧氏距离的变化基本相反, 当余弦相似度在高位取值时, 欧氏距离总在低位取值, 证明了系统变量数据向量间的余弦相似度和欧氏距离是相互关联的.

    图 3  余弦相似度和欧氏距离取值示意
    Fig. 3  Possible taking-value for cosine similarity and Euclidian distance
    图 4  余弦相似度和欧氏距离关联性示意
    Fig. 4  Relevance between cosine similarity and Euclidian distance

    本文以数据向量间的余弦相似度和欧氏距离组成的二元组来表征系统运行状态, 避免了直接对数据进行建模, 将高维的数据向量转化为二维数据, 有效降低了计算复杂度.同时, 考虑数据的时间特性, 更加符合工控系统的动态特性.

    为了求解出状态迁移图, 需要知道图的顶点和各顶点之间的转换关系.状态迁移图的顶点表示了系统可能的状态, 本文用一个状态矩阵求解状态迁移图, 需要知道图的顶点和各顶点之间的转换关系.状态迁移图的顶点表示了系统可能的状态, 本文以状态矩阵$STATE$表示可能的状态, $STATE$的每一行称为一个状态点, 一个状态点表示一个状态, 对应状态迁移图中的一个顶点.

    2.2.1   余弦相似度和欧氏距离的求解

    在某一时刻所有系统变量取值的集合称为一样本点, 例如${{\mathit{\boldsymbol{m}}}_i} = \{ {{s_{1, i}}, {s_{2, i}}, \cdots , {s_{l + p, i}}} \}$表示系统在$i$时刻的样本点, ${s_{k, t}} = {s_k}( t )$为第$k$个系统变量在$t$时刻的取值.为了计算方便, 用行向量的形式来表示样本点, 即${{\mathit{\boldsymbol{m}}}_i} = [ {{s_{1, i}}, {s_{2, i}}, \cdots , {s_{l + p, i}}} ]$.假设训练样本点数为$n$, 相邻样本点(经过z-score标准化的样本点)间的余弦相似度由式(5)计算得出.

    $ cossi{m_i} = \frac{{{{\mathit{\boldsymbol{m}}}_i} \times {{\mathit{\boldsymbol{m}}}}^{\rm T}_{i + 1}}}{{\left\| {{{\mathit{\boldsymbol{m}}}_i}} \right\| \times \left\| {{{\mathit{\boldsymbol{m}}}_{i + 1}}} \right\|}}, \quad i = 1, 2, \cdots , n - 1 $

    (5)

    ${\mathit{\boldsymbol{cossim}}}$是所有训练样本点任意相邻两点间余弦相似度组成的列向量.相邻样本点间的欧氏距离由式(6)计算得出.

    $ e{d_i} = \sqrt {{{{{\mathit{\boldsymbol{m}}}}}_i} \times {{\mathit{\boldsymbol{m}}}}^{\rm T}_{i + 1}} , \quad i = 1, 2, \cdots , n{\rm{ - }}1 $

    (6)

    令${\mathit{\boldsymbol{ed}}} = {[ {e{d_1}, e{d_2}, \cdots , e{d_{n - 1}}} ]^{\rm{T}}}$, ${\mathit{\boldsymbol{ed}}}$是所有训练样本点任意相邻两点间欧氏距离组成的列向量.

    2.2.2   cossim和ed量化

    由于${\mathit{\boldsymbol{cossim}}}$和${\mathit{\boldsymbol{ed}}}$中的元素都是实数, 这样通过二者确定的状态图中的状态数目将是无限个, 所以需要对它们进行量化处理, 使得${\mathit{\boldsymbol{cossim}}}$和${\mathit{\boldsymbol{ed}}}$中的元素只能取一些固定值.可以用一些量化阈值将${\mathit{\boldsymbol{cossim}}}$和${\mathit{\boldsymbol{ed}}}$的取值分成数个量化区间, 同一量化区间的点用相同的值来替换.本文用$intervals$表示${\mathit{\boldsymbol{cossim}}}$和${\mathit{\boldsymbol{ed}}}$量化区间的数量.

    为便于处理, 首先对${\mathit{\boldsymbol{cossim}}}$和${\mathit{\boldsymbol{ed}}}$进行归一化处理.将${\mathit{\boldsymbol{cossim}}}$归一化到$[-1, 1]$区间, 得到$\pmb{comsim}'$; ${\mathit{\boldsymbol{ed}}}$归一化到$[0, 1]$区间, 得到${\mathit{\boldsymbol{ed}}}'$.令$\pmb{comsim}'$的量化阈值集为

    $ {\mathit{\boldsymbol{r}}} = \left\{ {{r_0}, {r_1}, {r_2}, \cdots , {r_{\rm intervals}}} \right\} $

    (7)

    ${\mathit{\boldsymbol{ed}}}'$的量化阈值集为

    $ {\mathit{\boldsymbol{d}}} = \left\{ {{d_0}, {d_1}, {d_2}, \cdots , {d_{\rm intervals}}} \right\} $

    (8)

    按照各量化区间点数相等的原则计算${\mathit{\boldsymbol{r}}}$和${\mathit{\boldsymbol{d}}}$.定义函数$k:{\boldsymbol{\rm{R}}} \to {\boldsymbol{\rm{R}}}$, 表示某实数集合${{P}}$中元素的个数.定义函数

    $ f(x) = k(\{ j|j \le x, j \in {\mathit{\boldsymbol{cossim}}}'\} ) $

    $ g(x) = k(\{ j|j \le x, j \in {\mathit{\boldsymbol{ed}}}'\} ) $

    ${\mathit{\boldsymbol{r}}}$和${\mathit{\boldsymbol{d}}}$的计算如下:

    $ {r_k} = \begin{cases} - 1\text{, }\quad\ k = 0\\ \arg\min \limits_{x \in \left[ { -1, 1} \right]} f\left( x \right) - f\left( {{r_{k - 1}}} \right) = \dfrac{n}{{intervals}}\text{, }\\{\kern 54pt}\quad\ \, k = 1, 2, \cdots , intervals - 1\\ 1\text{, }\quad\ \, \;\;\;\;\;\; k = intervals \end{cases} $

    (9)

    $ {d_k} = \begin{cases} 0\text{, }\quad\ k = 0\\ \arg\min\limits_{x \in \left[ { 0, 1} \right]} g\left( x \right) - g\left( {{d_{k-1}}} \right) = \dfrac{n}{{intervals}}\text{, }\\{\kern 54pt}\quad\ \, k = 1, 2, \cdots , intervals - 1\\ 1\text{, }\quad\ \;\;\;\;\;\; k = intervals \end{cases} $

    (10)
    2.2.3   状态矩阵生成

    确定量化阈值后, 根据量化阈值对$\pmb{comsim}'$和${\mathit{\boldsymbol{ed}}}'$进行量化处理.

    $ \begin{align} & qcossi{m_i} = \arg\min\limits_{r \in {{\mathit{\boldsymbol{r}}}}} {\kern 1pt} cossi{m_i} \le r, \nonumber\\&\qquad\qquad\qquad\qquad\qquad\quad i = 1, 2, \cdots , n - 1 \end{align} $

    (11)

    $\pmb{qcomsim}'\!\!=\!\!{[ {qcossi{m_1}|qcossi{m_2}| \cdots |qcossi{m_{n-1}}} ]^{\rm{T}}}$为量化后的$\pmb{comsim}'$.

    $ qe{d_i} = \arg\min\limits_{d \in {\mathit{\boldsymbol{d}}}} e{d_i} \le d, \quad i = 1, 2, \cdots , n - 1 $

    (12)

    ${\mathit{\boldsymbol{qed}}}' = [ {qe{d_1}|qe{d_2}| \cdots |qe{d_{n - 1}}} ]$为量化后的${\mathit{\boldsymbol{ed}}}'$.最终的状态矩阵为$STATE = [ {\pmb{qcossim}'|{\mathit{\boldsymbol{qed}}}'} ]$.

    状态迁移图是一个有向图, 作为异常检测模型, 其顶点由$STATE$确定, 边用邻接矩阵表示.

    2.3.1   顶点确定

    由于对相邻样本点间的余弦相似度和欧氏距离进行了量化, 状态矩阵$STATE$中含有相同的状态点.相同的状态点对应状态迁移图中的同一个顶点, 应将其合并, 具体做法为:

    步骤 1.  用自然数对$STATE$中的每个状态点进行标号, 称为状态号, 并用向量${\mathit{\boldsymbol{sign}}}$来表示, 例如${\mathit{\boldsymbol{sign}}}\left( k \right)$表示$STATE$中的第$k$个状态点对应的状态号.状态号唯一标明了状态迁移图中的状态(即顶点).

    步骤 2.  令$STATE$第1个状态点的状态号为1, 即${\mathit{\boldsymbol{sign}}}(1) = 1$, 令状态数$sn = 1$.

    步骤 3.  从第2个状态点开始, 对第$k$个状态点遍历$STATE$中该状态点之前的状态点, 将该状态点与之前的状态点进行比较, 若找到相同的状态点$j$则停止遍历, 将该点的状态号设为第$j$个状态点的状态号, 即${\mathit{\boldsymbol{sign}}}(k) = {\mathit{\boldsymbol{sign}}}(j)$; 若没有找到相同的状态点, 则令, 将该点的状态号设为$sn$, 即${\mathit{\boldsymbol{sign}}}(k) = sn$.

    经过步骤1 $\sim$ 3后, $STATE$中相同的状态点具有相同的状态号, 最终的状态数为$sn$.状态迁移图的顶点对应$STATE$中状态号为1 $\sim$ $sn$的状态点.

    2.3.2   邻接矩阵确定

    用$NM$表示邻接矩阵, 邻接矩阵反映了状态迁移图中的状态转移关系, 若${NM_{i, j}} = 1$, 则状态迁移图中的第$j$个顶点到第$j$个顶点可达; 若${NM_{i, j}}$ $=$ $0$, 则第$i$个顶点到第$i$个顶点不可达.由于样本点是随时间增长逐渐采样而来的, 因此前一样本点到当前样本点总是可达的.相应地, $STATE$中前一状态点到当前状态点也总是可达的.根据这一原则可以求解出邻接矩阵$NM$, 具体做法为:

    步骤 1.  将$NM$置为全零.

    步骤 2.  根据前一样本点(状态点)到当前样本点(状态点)可达的原则, 令$NM({\mathit{\boldsymbol{sign}}}( i )$, ${\mathit{\boldsymbol{sign}}}( i$ $+$ $1 ))$ $= 1$, 由此计算出$NM$.

    确定了顶点和邻接矩阵后, 就能生成状态迁移图.

    状态迁移图生成后, 可以据此进行在线异常检测.基于状态迁移图的在线检测流程如图所示, 具体步骤如下:

    步骤 1.  获取当前样本点.

    步骤 2.  计算当前样本点与上一样本的余弦相似度$t\_cossi{m_i}$和欧氏距离$t\_e{d_i}$.

    步骤 3.  按照训练样本的标准对$t\_cossi{m_i}$和$t\_e{d_i}$进行归一化处理, 得到$qt\_cossi{m_i}$和$qt\_e{d_i}$.

    步骤 4.  遍历$STATE$, 判断$STATE$中是否有和当前点相同的点.若存在, 根据${\mathit{\boldsymbol{sign}}}$确定当前状态点对应的状态号.

    步骤 5.  若当前状态点对应的状态不在状态迁移图内, 产生告警, 记录异常类别为"异常1";若当前点对应的状态位于状态图内, 即$STATE$中有和当前点相同的点, 则根据状态迁移图判断上一节点对应状态到当前节点对应状态是否可达.若不可达, 则产生告警, 记录异常类别为"异常2";若可达, 继续下一个样本点的检测.

    目前工控领域没有一个标准的测试集, 而实际工控系统在攻击下的数据样本很难获取.针对这一情况, 本文采用田纳西-伊斯曼(Tennessee-eastman, TE)过程[18]的MATLAB模型作为仿真平台.根据第1节建立的恶意数据攻击模型, 选定斜坡注入和偏置注入两种攻击形式, 在特定时刻对TE过程实施攻击, 生成系统正常运行数据和攻击情形下的数据, 并用第2.4节描述的检测方法进行检测.

    图 5  基于状态迁移图的在线检测流程
    Fig. 5  Online detection process based on state transition graph

    TE过程的原型由Downs和Vogel在1993年根据一个真实的化工过程创建. TE过程主要部件有冷凝器、反应器、循环压缩机、汽提塔和气液分离器等.该过程包含两个放热反应和两个副反应, 生成产品$G$, $H$和副产品$F.$四种化学反应如下所示:

    $ {{A}}\left( { {g}} \right) + {{C}}\left( { {g}} \right) + {{D}}\left( { {g}} \right) \to {{G}}\left( {{ {liq}}} \right) $

    (13)

    $ {{A}}\left( { {g}} \right) + {{C}}\left( { {g}} \right) + {{E}}\left( { {g}} \right) \to {{H}}\left( {{ {liq}}} \right) $

    (14)

    $ {{A}}\left( { {g}} \right) + {{E}}\left( { {g}} \right) \to {{F}}\left( {{ {liq}}} \right) $

    (15)

    $ {{D}}\left( { {g}} \right) \to 2{{F}}\left( {{ {liq}}} \right) $

    (16)

    TE过程包含41个测量变量和12个操纵变量, 是一个多变量, 多数据, 复杂的非线性控制系统, 常用于多变量控制、非线性控制和故障诊断等领域. TE过程的MATALB模型由Ricker[19]给出, 该模型采用文献[20]中的控制策略, 仿真时间为48 h, 采样时间为3 min, 每次仿真共产生960组数据, 每组数据中都包含高斯白噪声. TE过程的工艺流程如图 6所示.

    图 6  基于状态迁移图的在线检测流程
    Fig. 6  Online detection process based on state transition graph

    反应器是TE过程的核心部件, 是发生化学反应的地方, 对温度的要求很高.因此本文选定攻击对象为反应器的温度传感器, 并根据第1节所述攻击模型, 选定偏置注入和斜坡注入两种攻击方式.

    $ {\tilde s_i}\left( t \right) = \begin{cases} {s_i}\left( t \right) + c\text{, }& t \in {T_a}\\ {s_i}\left( t \right)\text{, }& t \notin {T_a} \end{cases} $

    (17)

    $ {\tilde s_i}\left( t \right) = \begin{cases} {s_i}\left( t \right) + k(t - {t_0})\text{, }& t \in {T_a}\\ {s_i}\left( t \right)\text{, }& t \notin {T_a} \end{cases} $

    (18)

    其中, 偏置注入中$c$为常数, 而斜坡注入在${t_0}$时注入斜率为$k$的斜坡信号.反应器温度控制回路如图 7所示, 两个输入分别代表反应器温度初始设定值(用输入点9表示)和反应器温度(用xmeas 9表示), 输出为反应器冷却水流量(用xmv10表示).该回路通过反应器温度测量值和反应器温度初始设定值计算冷却水流量, 以此来控制反应器温度.为达到攻击反应器温度的目的, 在反应器温度输入端口增加Add模块, 将其与攻击信号叠加, 攻击信号通过延时模块在特定时刻发挥作用.图 8为加入攻击信号后的反应器温度控制回路, 延时模块时间设置为20 h, 攻击信号为常数, 代表偏置注入信号.攻击信号可以替换为其他信号, 以此设计更多的攻击形式.

    图 7  反应器温度控制回路
    Fig. 7  Control loop for reactor temperature
    图 8  加入攻击信号的反应器温度控制回路
    Fig. 8  Control loop for reactor temperature added with attack signal

    采集系统正常运行时各变量的数据, 生成数据集合Normal.选定式(17)中$c$为0.01, 0.1和1, 在系统仿真20 h加入注入攻击, 生成偏置注入下的数据集Dataset1, Dataset2, Dataset3.选定式(18)中$k$为0.01, 0.1和1, 在仿真20 h实施斜坡攻击, 生成斜坡注入下的数据集Dataset4, Dataset5, Dataset6.各数据集及其对应参数见表 1.这些数据集将用来测试本文所提的方法.

    表 1  测试数据集及其参数
    Table 1  The test data set and the corresponding parameters
    实验数据$c$$k$
    Normal00
    Dataset10.010
    Dataset20.10
    Dataset310
    Dataset400.01
    Dataset500.1
    Dataset601
    下载: 导出CSV 
    | 显示表格

    为直观感受加入攻击信号后对反应器温度的影响, 用MATLAB采集生成了反应器温度在各种情况下的温度-时间变化曲线, 分别如图 9 $\sim$ 11所示.图 9为正常工况条件下的温度-时间变化曲线, 可以看出反应器温度稳定在123℃左右. 图 10为斜坡注入$k=0.01$情况下的温度变化情况, 可以看到在20 h后反应器温度呈线性下降趋势. 图 11为偏置注入$c=0.1$情况下的温度-时间变化曲线, 反应器温度在20 h突然下降, 之后则稳定下来.

    图 9  正常工况下反应器温度随时间变化情况
    Fig. 9  Reactor temperature varies with time under normal condition
    图 10  斜坡注入工况($k=0.01$)下, 反应器温度随时间变化情况
    Fig. 10  Reactor temperature varies with time under ramp signal injection with $k$ set at 0.01
    图 11  偏置注入工况($c=0.1$)下, 反应器温度随时间变化情况
    Fig. 11  Reactor temperature varies with time under bias signal injection with $c$ set at 0.1

    由于状态迁移图通过正常数据训练得出, 所以用表 1中的Normal数据集来进行训练.根据第2.2节的分析, 生成的状态迁移图与余弦相似度和欧氏距离的量化区间数$intervals$相关, 因此本文不断改变$intervals$的取值, 进行了多组状态图测试.首先按照第2.2节和第2.3节的方式对Normal数据集进行处理, 然后设定量化区间数$intervals$, 计算状态矩阵, 最后生成状态迁移图.状态迁移图的顶点数和边数随$intervals$的变化情况如表 2所示.

    表 2  状态迁移图的顶点数和边数随$intervals$的变化
    Table 2  The nodes and triangles number of state transition graph varies with $intervals$
    $intervals$顶点数边数
    3948
    520147
    844315
    1073677
    1297782
    15143871
    下载: 导出CSV 
    | 显示表格

    表 2可以看出, 随着$intervals$的增大, 状态迁移图的顶点数和边数也不断增大.而$intervals$越大, 意味着将余弦相似度和欧氏距离的量化区间划分得越小, 由于状态迁移图的检测由顶点和边来决定, 所以顶点数和边数越多, 意味着检测规则的粒度越细, 相应地对异常会更加敏感, 此外还会额外增加资源开销. 图 12图 13分别显示了$intervals$为5和8时用MATLAB生成的状态迁移图.

    图 12  $intervals$为5时的状态迁移图
    Fig. 12  The state transition graph when $intervals$ is equal to 5
    图 13  $intervals$为8时的状态迁移图
    Fig. 13  The state transition graph when $intervals$ is equal to 8

    按照第2.4节的在线检测步骤, 用表 1的测试数据对状态迁移图检测模型的检测效果进行测试, 选定$intervals$为5, 8, 10, 15, 检测结果分别见表 3 $\sim$ 6.本文选用从攻击到正确检测到异常的时间以及误报率作为评价的标准.

    表 3  $intervals=5$时的检测结果
    Table 3  Detection results when $intervals$ is equal to 5
    Dataset1Dataset2Dataset3Dataset4Dataset5ataset6
    正确检测到异常的样本点数409401401436407402
    异常类别异常2异常2异常2异常2异常2异常2
    误报率(%)0.630.210.420.210.830.42
    下载: 导出CSV 
    | 显示表格
    表 4  $intervals=8$时的检测结果
    Table 4  Detection results when $intervals$ is equal to 8
    Dataset1Dataset2Dataset3Dataset4Dataset5Dataset6
    正确检测到异常的样本点数401401401411401401
    异常类别异常2异常2异常2异常2异常2异常2
    误报率(%)5.624.385.214.175.835.42
    下载: 导出CSV 
    | 显示表格
    表 5  $intervals=10$时的检测结果
    Table 5  Detection results when $intervals$ is equal to 10
    Dataset1Dataset2Dataset3Dataset4Dataset5Dataset6
    正确检测到异常的样本点数401401401401401401
    异常类别异常2异常2异常2异常2异常2异常2
    误报率(%)7.507.088.128.216.677.71
    下载: 导出CSV 
    | 显示表格
    表 6  $intervals=15$时的检测结果
    Table 6  Detection results when $intervals$ is equal to 15
    Dataset1Dataset2Dataset3Dataset4Dataset5Dataset6
    正确检测到异常的样本点数401401401401401401
    异常类别异常2异常1异常1异常1异常1异常1
    误报率(%)24.3726.2522.9230.2128.9626.04
    下载: 导出CSV 
    | 显示表格

    从检测时间来看, 在各种攻击情况下检测模型检测到异常的时间随着$intervals$的变化而变化.从检测最坏的结果对应的数据集来看(对应的数据集为Dataset4, 此时$c = 0$, $k = 0.01$), 当$intervals$为5和8时候, 分别在第436和第411个样本点检测到异常(在第20小时进行注入攻击, 对应的第一个异常样本点为第401个).而当$intervals$为10和15时, 检测模型在第401个样本点即检测到异常, 说明当$intervals$增大时, 检测模型对异常更加敏感, 检测异常的速度有所提升.

    从误报率来看, 当$intervals$逐渐增大时, 相应地误报率也随之增大.当$intervals$为5时误报率均不超过1 %.当$intervals$增大到8时, 误报率迅速增加到5 %左右.而当$intervals$取15时, 误报率超过20 %, 这时可认为状态迁移图不能正确检测异常.

    由于工控领域缺乏标准的数据集, 研究人员都用各自的模型和数据进行工控系统异常的研究, 很难复现, 难以对比各种方法的好坏.为了进一步评估状态迁移图检测的性能, 将其与控制系统常用的异常检测方法---主元分析法(Principal component analysis, PCA)进行比较. PCA采用${{\rm{T}}^2}$统计量和$\rm{SPE}$统计量作为检测异常的指标, 当样本的${{\rm{T}}^2}$统计量和$\rm{SPE}$统计量超过各自的控制限时, 检测到异常.选定状态迁移图模型检测结果最差的数据集Dataset4作为PCA方法的测试集, 其测试结果如图 14图 15所示.从图 14图 15可以看出, PCA虽然最终能检测到异常, 但是在第500个样本点才检测出来.相比于状态迁移图检测模型, 检测到异常的速度较慢.

    图 14  PCA方法${{\rm{T}}^2}$统计量
    Fig. 14  ${{\rm{T}}^2}$ statistic of PCA method
    图 15  PCA方法$\rm{SPE}$统计量
    Fig. 15  $\rm{SPE}$ statistic of PCA method

    本节从时空资源消耗以及检测性能两方面对状态迁移图检测模型进行讨论, 由于训练过程是通过历史数据离线进行的, 因此时空消耗只考虑检测时的情况.

    1) 从时间消耗上来看, 对新样本点的检测需要判断新样本点是否位于状态迁移图内, 即需要遍历一次状态图, 若状态迁移图含$n$个顶点, 则时间复杂度为O$(n)$.实际检测过程中, 顶点数与设置的量化区间数$intervals$有关.通过第3.4节的实验结果不难发现, 当$intervals$取5或8时已有较好的检测结果, 继续增大$intervals$只会增大资源消耗, 同时增加误报率.且$intervals$越大, 需要越多的训练样本进行训练.以$intervals$取10为例, 顶点数$n$最多为100, 而实际会小于100, 因此这一步骤需要的计算量很小.第二个步骤判断上一状态到当前状态是否可达, 只需查询邻接矩阵中对应的数值即可, 时间复杂度为O$(1)$.

    2) 从空间消耗上来看, 本文的异常检测方法需要存储量化阈值集合${\mathit{\boldsymbol{r}}}$和${\mathit{\boldsymbol{d}}}$, 状态矩阵$STATE$, 邻接矩阵$NM$, 其中邻接矩阵$NM$占用的存储空间为$n$的二阶次, ${\mathit{\boldsymbol{r}}}$和${\mathit{\boldsymbol{d}}}$, 状态矩阵$STATE$存储消耗为$n$的一阶次.但由于$n$取值很小(不超过100), 且邻接矩阵存储单位为比特, 当$intervals$取10时, 估计整体消耗的存储空间不超过1 MB.

    3) 从检测性能来看, 即便对于微小的攻击(对应数据集Dataset1和Dataset4), 检测模型也能较快检测出异常, 相比于常规的异常检测算法, 例如第3.4节的PCA算法, 状态迁移图检测模型能够更好地刻画系统运行的动态变化, 更快地检测出异常.但是当$intervals$较大(取10以上)时, 检测误报率较高.而当$intervals$较小(例如, 取5时)时, 不能快速地对微小攻击做出反应.因此需要谨慎选择$intervals$以平衡误报率和检测速度之间的关系.

    本文提出了一种新的基于状态迁移图的异常检测方法, 详细描述了其状态表示、状态图生成和在线检测过程, 建立了工控系统的恶意数据攻击模型, 并基于TE仿真模型进行了各项实验, 最后从时空资源消耗和检测性能两方面对方法进行了讨论.该方法具有较小的时空消耗, 且考虑了系统运行的动态特性, 能够快速检测出异常, 因而适用于资源受限和对实时性要求较高的工业控制系统.

  • 图  1  FootScan 足底压力测量平板系统

    Fig.  1  FootScan plantar pressure measurement

    图  2  整个着地周期的5个特征时刻

    Fig.  2  Five key moments of the whole stance phase

    图  3  ELM神经网络结构图

    Fig.  3  The structure of ELM neural network

    图  4  足底压力数据图

    Fig.  4  Plantar pressure map

    图  5  正常测试者与ACL断裂测试者的左侧足底压力COP轨迹线

    Fig.  5  Left feet's COP lines of normal and ACLD participators

    图  6  正常测试者与ACL断裂测试者的右侧足底压力COP轨迹线

    Fig.  6  Right feet's COP lines of normal and ACLD participators

    图  7  COP轨迹线放缩效果图

    Fig.  7  Zoom of COP line

    表  1  正常测试者左侧足底压力数据特征序列表

    Table  1  COP feature vectors of left feet in the normal group

    序号 组别 $C_1x$ $C_1y$ $C_2x$ $C_2y$ $C_3x$ $C_3y$ $C_4x$ $C_4y$ $C_5x$ $C_5y$
    1 $A$ 2 12 5 12 24 14 26 15 26 14
    2 3 11 4 12 17 12 23 12 27 14
    3 2 13 3 13 10 13 23 12 27 14
    4 3 13 5 13 15 11 23 12 28 13
    5 3 13 5 13 14 11 22 11 28 11
    6 $B$ 3 12 5 13 13 11 22 11 27 11
    7 3 13 4 13 13 10 22 10 28 13
    8 3 13 6 13 15 8 21 7 27 7
    9 3 13 5 13 14 11 22 11 28 11
    10 3 12 5 13 13 11 22 11 27 11
    11 $C$ 3 12 5 13 10 11 21 9 28 12
    12 3 10 5 9 12 9 21 10 27 15
    13 2 13 9 10 13 10 11 10 26 13
    14 2 9 4 9 8 9 21 9 27 15
    15 3 9 5 10 6 10 23 13 28 14
    16 $D$ 3 8 4 9 8 9 22 12 27 15
    17 3 8 5 9 13 8 22 10 28 16
    18 3 12 8 11 16 9 21 10 27 13
    19 2 13 7 12 13 10 22 10 27 13
    20 3 13 8 11 19 8 21 9 27 13
    21 $E$ 3 12 10 10 17 9 21 9 27 12
    22 3 14 5 14 13 12 6 14 27 9
    23 2 15 5 15 7 15 23 10 27 8
    24 3 15 8 13 14 9 10 11 28 7
    25 3 14 5 15 16 10 22 9 26 10
    下载: 导出CSV

    表  2  左膝ACL断裂测试者左侧足底压力数据特征序列表

    Table  2  COP feature vectors of left feet in the group of left side ACLD

    序号 组别 $C_1x$ $C_1y$ $C_2x$ $C_2y$ $C_3x$ $C_3y$ $C_4x$ $C_4y$ $C_5x$ $C_5y$
    1 $A$ 2 12 4 12 11 12 20 11 27 14
    2 2 11 9 9 15 9 21 9 27 15
    3 2 10 6 9 11 9 21 10 28 15
    4 3 10 11 7 15 9 21 9 28 15
    5 3 11 6 11 14 9 20 9 27 13
    6 $B$ 3 10 11 7 15 9 21 9 28 15
    7 3 10 11 6 13 7 20 8 27 15
    8 2 10 7 10 12 9 20 8 27 14
    9 4 12 5 13 13 11 20 11 27 13
    10 2 11 5 12 13 12 21 12 27 14
    11 $C$ 3 13 4 14 10 14 23 12 28 13
    12 3 11 5 11 12 12 22 13 28 14
    13 2 12 4 13 10 12 21 12 27 15
    14 2 13 5 13 12 12 22 12 28 14
    15 2 13 7 11 11 10 21 10 26 11
    16 $D$ 2 13 8 10 12 9 20 7 23 9
    17 3 13 6 11 14 9 23 9 24 9
    18 3 14 8 12 14 10 21 9 24 9
    19 1 13 6 12 15 7 21 8 24 10
    20 2 9 5 9 15 10 22 12 27 13
    21 $E$ 3 9 14 6 19 7 21 10 27 12
    22 2 10 8 7 14 8 21 11 27 13
    23 2 9 8 7 18 8 21 10 27 13
    24 3 11 6 11 14 8 20 9 27 13
    25 3 10 5 10 15 7 21 10 28 14
    下载: 导出CSV

    表  3  左侧足底压力数据聚类分析辨识结果

    Table  3  The result of left plantar pressure analysis

    序号 训练组别 训练数据个数 测试组别 测试数据个数 辨识准确率 (%)
    1 $A, B, C, D$ 40 $E$ 10 90
    2 $A, B, C, E$ 40 $D$ 10 50
    3 $A, B, D, E$ 40 $C$ 10 60
    4 $ A, C, D, E$ 40 $B$ 10 100
    5 $B, C, D, E$ 40 $A$ 10 80
    下载: 导出CSV

    表  4  正常测试者右侧足底压力数据特征序列表

    Table  4  COP feature vectors of right feet in the normal group

    序号 组别 $C_1x$ $C_1y$ $C_2x$ $C_2y$ $C_3x$ $C_3y$ $C_4x$ $C_4y$ $C_5x$ $C_5y$
    1 $A$ 3 8 5 7 14 8 7 7 27 6
    2 3 10 4 8 10 8 22 5 27 5
    3 2 8 4 8 12 8 21 9 27 7
    4 2 5 4 4 11 7 22 9 28 9
    5 2 4 4 4 13 8 22 10 27 9
    6 $B$ 3 5 5 4 5 4 22 10 28 10
    7 3 4 4 4 16 9 23 9 28 9
    8 3 5 4 5 13 8 22 9 28 8
    9 2 4 4 4 13 8 22 10 27 9
    10 3 5 5 4 5 4 22 10 28 10
    11 $C$ 3 5 4 4 12 7 22 9 28 10
    12 4 7 5 7 12 8 21 8 23 9
    13 2 7 6 6 16 7 22 6 26 3
    14 3 9 4 9 14 8 21 6 27 4
    15 2 5 5 5 13 6 21 7 27 4
    16 $D$ 3 6 4 6 12 7 21 8 27 5
    17 3 7 8 8 16 8 21 9 27 6
    18 2 8 7 10 12 10 22 10 28 7
    19 2 7 4 8 15 9 9 9 27 6
    20 3 7 9 8 15 8 22 7 28 4
    21 $E$ 3 9 9 12 16 12 22 9 27 5
    22 3 4 5 4 14 6 5 4 27 9
    23 3 7 5 7 13 10 23 13 27 14
    24 2 5 5 5 10 6 6 5 27 8
    25 3 6 4 6 7 6 4 6 27 8
    下载: 导出CSV

    表  5  右膝ACL断裂测试者右侧足底压力数据特征序列表

    Table  5  COP feature vectors of right feet in the group of right side ACLD

    序号 组别 $C_1x$ $C_1y$ $C_2x$ $C_2y$ $C_3x$ $C_3y$ $C_4x$ $C_4y$ $C_5x$ $C_5y$
    1 $A$ 2 4 8 8 17 10 21 10 27 8
    2 3 7 4 7 15 9 20 9 27 7
    3 2 6 5 6 17 10 21 9 27 8
    4 3 6 7 8 18 9 22 9 27 8
    5 2 6 7 8 17 10 21 10 27 8
    6 $B$ 2 6 7 7 13 7 23 8 27 8
    7 3 4 5 4 14 7 22 9 27 11
    8 2 4 7 6 12 6 23 8 27 9
    9 3 5 8 7 14 8 23 8 27 9
    10 2 6 5 6 15 7 22 8 27 8
    11 $C$ 2 4 5 5 12 6 21 9 27 8
    12 3 5 5 6 12 8 22 10 27 10
    13 2 7 6 7 16 9 14 10 27 6
    14 2 8 11 11 16 10 21 9 27 6
    15 3 7 5 8 16 9 21 9 27 5
    16 $D$ 2 6 6 7 13 9 21 9 27 4
    17 2 6 7 8 16 9 21 8 27 4
    18 3 7 6 8 17 10 22 9 28 6
    19 2 6 4 6 13 8 22 9 28 8
    20 2 7 4 7 10 7 21 9 27 7
    21 $E$ 3 9 6 9 12 10 21 10 27 6
    22 3 7 6 7 15 9 22 8 27 7
    23 3 6 6 6 17 8 23 7 28 6
    24 3 6 5 6 16 9 21 10 27 5
    25 3 7 8 8 15 10 22 11 28 8
    下载: 导出CSV

    表  6  右侧足底压力数据聚类分析辨识结果

    Table  6  The result of right plantar pressure analysis

    序号 训练组别 训练数据个数 测试组别 测试数据个数 辨识准确率 (%)
    1 $ A, B, C, D$ 40 $E$ 10 80
    2 $A, B, C, E$ 40 $D$ 10 70
    3 $ A, B, D, E$ 40 $C$ 10 80
    4 $ A, C, D, E$ 40 $B$ 10 70
    5 $ B, C, D, E$ 40 $A$ 10 80
    下载: 导出CSV
  • [1] Taylor A J, Menz H B, Keenan A M. The influence of walking speed on plantar pressure measurements using the two-step gait initiation protocol. The Foot, 2004, 14 (1): 49-55
    [2] Stolwijk N M, Duysens J, Louwerens J W K, Keijsers N L. Plantar pressure changes after long-distance walking. Medicine and Science in Sports and Exercise, 2010, 42 (12): 2264-2272
    [3] Lai Y C, Lin H S, Pan H F, Chang W N, Hsu C J, Renn J H. Impact of foot progression angle on the distribution of plantar pressure in normal children. Clinical Biomechanics, 2014, 29 (2): 196-200 doi: 10.1016/j.clinbiomech.2013.11.012
    [4] Keijsers N L W, Stolwijk N M, Nienhuis B, Duysens J. A new method to normalize plantar pressure measurements for foot size and foot progression angle. Journal of Biomechanics, 2009, 42 (1): 87-90 http://www.academia.edu/14125878/A_new_method_to_normalize_plantar_pressure_measurements_for_foot_size_and_foot_progression_angle
    [5] Kohonen T. The self-organizing map. Proceedings of the IEEE, 1990, 78 (9): 1464-1480
    [6] Kohonen T. Self-Organizing Maps. New York: Springer-Verlag, 2001.
    [7] Carpenter G A, Grossberg S, Rosen D B. Fuzzy ART: fast stable learning and categorization of analog patterns by an adaptive resonance system. Neural Networks 1991, 4 (06): 759-771 doi: 10.1016/0893-6080(91)90056-B
    [8] Huang G B, Wang D H, Lan Y. Extreme learning machines: a survey. International Journal of Machine Learning and Cybernetics, 2011, 2 (2): 107-122 doi: 10.1007/s13042-011-0019-y
    [9] Stolwijk N M, Duysens J, Louwerens J W, van de Ven Y H, Keijsers N L. Flat feet, happy feet? Comparison of the dynamic plantar pressure distribution and static medial foot geometry between Malawian and Dutch adults. PLoS One, 2013, 8 (2): e57209 doi: 10.1371/journal.pone.0057209
    [10] Cavanagh P R, Rodgers M M. The arch index: a useful measure from footprints. Journal of Biomechanics, 1987, 20 (5): 547-551 doi: 10.1016/0021-9290(87)90255-7
    [11] Andriacchi T P, Dyrby C O. Interactions between kinematics and loading during walking for the normal and ACL deficient knee. Journal of Biomechanics, 2005, 38 (2): 293-298 doi: 10.1016/j.jbiomech.2004.02.010
    [12] 黄红拾, 于媛媛, 郭秦炜, 徐雁, 敖英芳.前交叉韧带断裂患者步行时足底压力时间特征.中国运动医学杂志, 2015, 34 (3): 271-274 http://www.cnki.com.cn/Article/CJFDTOTAL-YDYX201503009.htm

    Huang Hong-Shi, Yu Yuan-Yuan, Guo Qin-Wei, Xu Yan, Ao Ying-Fang. Temporal characteristics of plantar pressure in patients with anterior cruciate ligament rupture during walking. Chinese Journal of Sports Medicine, 2015, 34 (3): 271-274 http://www.cnki.com.cn/Article/CJFDTOTAL-YDYX201503009.htm
    [13] Huang G B, Chen Y Q, Babri H A. Classification ability of single hidden layer feedforward neural networks. IEEE Transactions on Neural Networks, 2000, 11 (3): 799-801 doi: 10.1109/72.846750
    [14] Feng G R, Huang G B, Lin Q P, Gay R. Error minimized extreme learning machine with growth of hidden nodes and incremental learning. IEEE Transactions on Neural Networks, 2009, 20 (8): 1352-1357 https://www.researchgate.net/publication/26665344_Error_Minimized_Extreme_Learning_Machine_With_Growth_of_Hidden_Nodes_and_Incremental_Learning
    [15] Huang G B, Zhou H M, Ding X J, Zhang R. Extreme learning machine for regression and multiclass classification. IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics), 2012, 42 (2): 513-529 http://www.doc88.com/p-9197171944428.html
    [16] DeFrate L E, Papannagari R, Gill T J, Moses J M, Pathare N P, Li G A. The 6 degrees of freedom kinematics of the knee after anterior cruciate ligament deficiency: an in vivo imaging analysis. American Journal of Sports Medicine, 2006, 34 (8): 1240-1246 https://www.researchgate.net/publication/7145174_The_6_Degrees_of_Freedom_Kinematics_of_the_Knee_After_Anterior_Cruciate_Ligament_Deficiency_An_In_Vivo_Imaging_Analysis
    [17] Ikeda H, Kurosawa H, Kim S G. Quadriceps torque curve pattern in patients with anterior cruciate ligament injury. International Orthopaedics, 2002, 26 (6): 374-376 doi: 10.1007/s00264-002-0402-0
    [18] Chiu M C, Wu H C, Chang L Y, Wu M H. Center of pressure progression characteristics under the plantar region for elderly adults. Gait and Posture, 2013, 37 (3): 408-412 doi: 10.1016/j.gaitpost.2012.08.010
  • 期刊类型引用(20)

    1. 巩彬,安爱民,石耀科,杜先君. 基于IMODA自适应深度信念网络的复杂模拟电路故障诊断方法. 电子科技大学学报. 2024(03): 327-344 . 百度学术
    2. 金亮,闫银刚,杨庆新,刘素贞,张闯. 小样本条件下永磁同步电机深度迁移学习性能预测方法. 电工技术学报. 2023(18): 4921-4931 . 百度学术
    3. 朱小勇,陈胜. 基于ResNet-ViT的海战多目标态势感知. 信息与控制. 2023(05): 638-647 . 百度学术
    4. 吴智怀,王志伟. 基于SDAE的液压泵故障诊断方法研究. 冶金设备. 2023(05): 15-22 . 百度学术
    5. 李锦键,王兴贵,杨维满,赵玲霞. 基于改进递归深度信念网络的CSP电站短期出力预测. 太阳能学报. 2022(07): 225-232 . 百度学术
    6. 沈浩宇,江先志,冯涛. 二轮平衡车的神经元控制算法及仿真. 机床与液压. 2022(19): 159-166 . 百度学术
    7. 王功明,乔俊飞,关丽娜,贾庆山. 深度信念网络研究现状与展望. 自动化学报. 2021(01): 35-49 . 本站查看
    8. 石媛媛,裴志利,姜明洋. 基于深度信念网络的文本分类研究综述. 内蒙古民族大学学报(自然科学版). 2021(01): 44-50 . 百度学术
    9. 李娜娜,胡坚剑,顾军华,张亚娟. 深度置信网络优化模型在人才评价中的应用. 计算机工程. 2020(02): 80-87+102 . 百度学术
    10. 唐魏,郑源,潘虹,徐晶珺. 引入动态调节学习率的SAE轴承故障诊断研究. 计算机工程与应用. 2020(20): 264-269 . 百度学术
    11. 李少波,李传江,胡建军,张安思,杨静. 基于深度置信网络的机械设备故障诊断研究综述. 现代制造工程. 2020(10): 156-162+142 . 百度学术
    12. 史科,陆阳,刘广亮,毕翔,王辉. 基于多隐层Gibbs采样的深度信念网络训练方法. 自动化学报. 2019(05): 975-984 . 本站查看
    13. 高磊,范冰冰,黄穗. 基于残差的改进卷积神经网络图像分类算法. 计算机系统应用. 2019(07): 139-144 . 百度学术
    14. 张德,李国璋,王怀光,张峻宁. 位姿估计自适应学习率的改进. 电子测量与仪器学报. 2019(06): 51-58 . 百度学术
    15. 邢海霞,程乐. 一种基于强化学习的深度信念网络设计方法. 控制工程. 2019(11): 2115-2120 . 百度学术
    16. 索明何,程乐. 基于PLSR的深度信念网输出权值确定方法. 控制工程. 2018(04): 668-676 . 百度学术
    17. 王功明,乔俊飞,王磊. 一种能量函数意义下的生成式对抗网络. 自动化学报. 2018(05): 793-803 . 本站查看
    18. 牟善仲,徐天赐,符奥,王萌,白茹. 基于自适应深度学习模型的变压器故障诊断方法. 南方电网技术. 2018(10): 14-19 . 百度学术
    19. 刘冰尧,庞腾,雷菊阳. 基于深度信念网络的变压器故障诊断. 化工自动化及仪表. 2018(11): 873-878 . 百度学术
    20. 王献锋,张传雷,张善文,朱义海. 基于自适应判别深度置信网络的棉花病虫害预测. 农业工程学报. 2018(14): 157-164 . 百度学术

    其他类型引用(31)

  • 加载中
图(7) / 表(6)
计量
  • 文章访问数:  2492
  • HTML全文浏览量:  230
  • PDF下载量:  715
  • 被引次数: 51
出版历程
  • 收稿日期:  2016-02-29
  • 录用日期:  2016-08-15
  • 刊出日期:  2017-03-20

目录

/

返回文章
返回