Upper Limb Spasticity Evaluation Based on Stretch Reflex Threshold: Method, Device, Validity and Reliability
-
摘要: 针对目前临床痉挛评定主观性大,信度与效度有待提高的问题,提出了一种新的基于牵张反射阈值的痉挛评定方法,利用关节角加速度变化判定牵张反射阈值,通过设计相应装置检验了牵张反射阈值在上肢痉挛评定中的信度与效度,并探讨了加速度变化在上肢痉挛评定中的信度与效度.4位检查者利用改良Ashworth量表(Modified Ashworth scale,MAS)及该痉挛检测装置对招募的22例伴随有上肢痉挛症状的受试者进行了痉挛评定.将评定过程中采集的牵张反射阈值以及加速度平均变化值(Acceleration mean variance,AMV)与MAS评分进行相关性分析,显示牵张反射阈值数据与MAS评分显著相关,相关性满足(r=-0.831~-0.953,P < 0.05),AMV与MAS评分相关性满足(r=0.665~0.900,P < 0.05).它们重测信度分别满足(r=0.890~0.962,P < 0.05)和(r=0.632~0.928,P < 0.05).实验结果表明该方法及装置可为痉挛评定提供一种实用的定量分析手段.
-
关键词:
- 痉挛 /
- 牵张反射阈值 /
- 加速度 /
- 改良Ashworth量表
Abstract: Clinical methods for spasticity evaluation are criticized as subjective, and their validity and reliability remain to improve. This study is to present a new evaluation method for accurate and reliable evaluation of spasticity based on a stretch reflex threshold which is obtained from the acceleration and angle data. A corresponding device is designed to test the validity and reliability of the stretch reflex threshold in spasticity assessment, and further study is conducted to discuss the validity and reliability of acceleration mean variance (AMV). Spasticity is evaluated in 22 subjects with upper limb spasticity using the constructed portable device and the modified Ashworth scale (MAS) by four evaluators. The results of correlation analysis show that there is a strong negative correlation (r=-0.831~-0.953, P < 0.05) between MAS and the stretch reflex threshold. Also, the correlation coefficient between the MAS and AMV is (r=0.665~0.900, P < 0.05). Their test-retest reliabilities are (r=0.890~0.962, P < 0.05) and (r=0.632~0.928, P < 0.05). These results demonstrate that the method and device may provide a quantitative and practical way for spasticity measurement.-
Key words:
- Spasticity /
- stretch reflex threshold /
- acceleration /
- modified Ashworth scale (MAS)
-
众所周知, 聚类, 支持向量机及神经元网络等常见分类方法都属于浅层分类方法, 在处理蕴藏隐含信息的样本分类问题方面还存在不足.传统的聚类方法对于高维数据来说, 数据样本较低维数据聚类时分布更为稀疏, 且每个数据间的距离都可能相当, 因此难以找到聚类中心, 从而不容易进行分类[1]; SVM属于有监督学习算法, 在处理小样本分类时有一定优势, 然而该方法过于依赖样本数据尺度, 且算法复杂度较高.并且SVM中核函数的选择往往决定了分类的精度和收敛速度, 分类结果存在不确定性[2]; 此外, 由于神经网络缺乏预训练机制, 难以深度挖掘数据中的隐含信息[3].然而深度学习方法以"无监督训练-有监督调解全局网络参数"的框架, 从理论上避免了传统神经元网络易陷入局部极值的缺点[4], 且在挖掘数据隐含信息方面具有独特优势, 尤其是在面临大规模样本数据的时候, 有更加突出的表现.常见的深度学习分为自动编码器, 卷积神经网络和深度置信网络等模型.通过查阅文献, 深度置信网络建模方法在图像处理[5-6]、软测量技术[7]、计算智能[8]等诸多领域得到成功应用, 建模精确度普遍有所提升, 上述理论的成功应用, 为构建基于高维数据非线性特征提取的深度置信网络分类器提供了重要的理论和方法支撑.值得一提的是, 面向含非线性特征的高维数据的深度置信网络建模问题, 信息量上的冗余往往给网络带来不必要的负荷.因此预先对样本数据做特征提取十分必要.
过程变量的特征提取的目的是找到数据之间的线性以及非线性关系表达, 而后利用提取的低维特征数据表征原有的高维数据.故数据之间的线性以及非线性关系的提取是提取特征的关键.常见的过程数据特征提取方法有主成分分析方法(PCA), 独立主元分析(ICA), 偏最小二乘法(PLS)等.其中, PCA利用高斯分布数据的特征, 将数据映射到正交的低维子空间上, 保留数据的特征[9]; ICA根据已经存在的统计值, 进行独立主成分正交变换[10]; PLS利用线性拟合对多变量建模, 减少变量个数[11].以上方法在数据满足高斯分布和有线性关系的情况下适用, 且效果很好, 但是, 在一类多变量数据且变量分布不定, 且存在非线性关系时, 以上方法并不奏效.所以, 本文应用一类基于输入训练神经网络表征非线性主元分析的方法, 旨在解决在多变量过程中的非线性特征提取问题, 并且实现数据降维, 为后续构建一类新的深度置信网络提供数据预处理的方法支撑.
空气固体细微污染物PM2.5的形成, 受众多复杂因素影响(已知影响因素超过20种)[12-13].就产生过程而言, PM2.5可以由污染源直接排出(称为一次粒子), 也可以是各污染源排出的气态污染物经过冷凝或在大气中发生复杂的化学反应而生成(称为二次粒子).特别地, 在已知的众多理化因素中, 有别于湿度、风速、降雨等, O$_3$属于驱散因子, 其浓度与PM2.5浓度之间呈指数衰减规律, 此外, 其他因素(光照等)与PM2.5浓度的关系还有待探索[14].因此, PM2.5浓度预测是一类典型的数据维度高, 且数据含非线性特征的建模问题, 传统的基于浅层学习的数据驱动建模方法[15-17]在预测精度上还有待提升, 且不具备对PM2.5浓度影响因素进行诊断的功能.
受上述讨论启发, 针对过程变量数据维数高, 且含复杂非线性特征, 数据间隐含信息难以利用等特点, 本文提出一类基于非线性特征提取的深度置信网模型, 旨在解决高维数据非线性特征提取以及数据特征中隐含信息挖掘的问题, 并对影响模型输出的关键变量进行诊断.最后, 以一类具体的多变量建模和诊断问题讨论所提方法的应用.本文结构安排如下:第1节展示了基于非线性特征提取的深度置信网络的建模过程; 第2节基于信息熵理论, 对改进后的深度置信网络的建模复杂度优势进行论证; 第3节以河北省某市的PM2.5监测数据为对象, 验证本文所提方法的有效性; 第4节给出结论与工作展望.
1. 基于非线性特征提取的深度置信网络
本节提出一类基于非线特征提取的深度置信网网络模型.基于非线性PCA提取原始数据特征, 实现数据预处理.同时计算各变量的统计量, 作为影响因素诊断依据.同时, 将预处理后的数据作为深度置信网的输入以构建预测模型.改进的深度置信网络结构在下文中具体介绍.
1.1 数据预处理
高维多变量过程数据(维度为$N$)之间存在的线性关系可以利用主成分分析的方法, 进行数据特征提取, 实现降维目的.然而, 数据之间存在复杂的非线性关系时, 理论上同样可以利用$A$个主元($A<N$)就可以反映出过程的主要信息.非线性PCA就是一种对${X}$的估计量${\hat{X}}$的非线性表示即:
$ \begin{equation} {X}={\hat{X}}+{E}=F({T_N})+{E} \end{equation} $
(1) 其中, $\hat{X}$是$X$的估计矩阵, ${E}$是残差矩阵, $F(\cdot)$是一个非线性函数, ${T_N}$我们称之为非线性主元得分矩阵.基于Tan等提出的输入训练(Input-training, IT)神经网络的方法[18].本文将IT网络的输入作为非线性主元得分矩阵, IT网络的输出作为原始样本的估计值, 网络调节权值的时候, 不仅调节网络内部的参数, 输入也随之变换.当网络训练完成的时候, 便可以得到${T_N}$, 同时也得到了非线性函数$F(\cdot)$.本文采用三层的输入训练网络, 如图 1所示.
整体网络采用快速下降法调节网络间的连接权值.网络的目标函数为$J$:
$ \begin{eqnarray*} J=\frac{1}{2}\sum\limits_{j=1}^{N}\sum\limits_{i=1}^{M}(X_{ij}-{\hat{X}}_{ij})^2=\frac{1}{2}\sum\limits_{j=1}^{N}{e}^2_j \end{eqnarray*} $
(2) $i$表示变量的维度, $j$表示输出数据的组数, $e_j$则表示每组训练数据的绝对误差.
1.2 深度置信网模型
深度学习是Hinton等在2006年提出的一种基于概率表达网络模型[19].深度学习的技术可以分为两部分:第一部分是利用无监督的学习来预训练每一层, 第二部分是全网络自上而下的微调权值.由于无监督的方式, 使用所有无标签数据, 所以过程变量包含监督学习所不能表达的隐含信息.本文所提出基于深度置信网的预测模型中, 网络输入是上一级降维后的非线性主元得分矩阵, 输出是预测分类结果.其网络结构如下图所示.本文采用三层的输入训练网络, 如图 1所示.
自上而下是多层的限制性玻尔兹曼机, 隐含层中每一层的输出, 作为下一层的输入.在这个训练阶段, 在可视层会产生一个向量${\pmb v} $, 通过它将值传递到隐层.反过来, 可视层的输入会去重构原始的输入信号[20].我们定义联合概率分布:
$ \begin{align} P({\pmb v}\mid {\pmb h})=\, &\dfrac{{\rm exp}(-E({\pmb v}, {\pmb h}))}{{\pmb z}}= \nonumber\\ &\dfrac{1}{{\pmb z}}\prod\limits_{ij}{\rm e}^{W_{ij}v_ih_j}\prod\limits_{i}{\rm e}^{b_iv_j}\prod\limits_{j}{\rm e}^{a_jh_j} \end{align} $
(3) 其中${\pmb z}$为:
$ \begin{equation} {\pmb z}=\sum\limits_{{\pmb h}, {\pmb v}}{\rm exp}(-E({\pmb v}, {\pmb h})) \end{equation} $
(4) $v_i$表示可视层第$i$个节点的输出, $h_j$表示隐含层第$j$个节点的输出.整体网络的参数$\theta=\{W, a, b\}$, $W$是权值参数$a$和$b$分别表示可视层和隐含层的偏置变量.给定可视层的前提下, 隐含层的概率为:
$ \begin{align} P({\pmb h}\mid{\pmb v})=\, &\prod\limits_{j}p(h_j\mid {\pmb v})p(h_j=1\mid {\pmb v})= \nonumber\\ & \dfrac{1}{1+{\rm exp}(\sum\limits_{i}W_{ij}v_i-a_j)} \end{align} $
(5) 这样我们就建立了可视层与隐含层之间的概率表示.同样的, 隐含层之间的概率可表示为:
$ \begin{align} P({\pmb v}, h_1, &h_2, h_3)=\nonumber \\ &P({\pmb v}\mid h_1)P(h_1\mid h_2)P(h_2 \mid h_3)\end{align} $
(6) 对于RBM的学习算法我们采用梯度衰减法.可视层表达的是输入数据的特征, 所以学习算法的目标函数是将可视层的概率最大化.所以有如下最大似然的概率表示:
$ \begin{align} \frac{\partial {\rm log_e}P({\pmb v})}{\partial \theta }=\, &\frac{\partial {\rm log_e}\sum\limits_{\pmb h}P({\pmb v}, {\pmb h})}{\partial \theta}=\nonumber\\ & \frac{\sum\limits_{\pmb h}{\rm e}^{-energy({\pmb v}, {\pmb h}) \frac{\partial(-energy({\pmb v}, {\pmb h}) }{\partial \theta}}}{\sum\limits_{\pmb h}-energy({\pmb v}, {\pmb h}))} - \nonumber\\ & \frac{\sum\limits_{\pmb v}\sum\limits_{\pmb h}{\rm e}^{(-energy({\pmb v}, {\pmb h})\frac{\partial energy({\pmb v}, {\pmb h}) }{\partial \theta})} }{\sum\limits_{\pmb v}\sum\limits_{\pmb h}(-energy({\pmb v}, {\pmb h}))} \end{align} $
(7) 对于标准化的高斯RBM,
$ \begin{align} energy({\pmb v}, {\pmb h})= \frac{1}{2}{\pmb v}^{\rm T}{\pmb v}-{ a}^{\rm T}{\pmb v}-{ b}^{\rm T}{\pmb h}-{\pmb h}^{\rm T}W{\pmb v}\nonumber\end{align} $
得到:
$ \begin{align} \frac{\partial P({\pmb v})}{\partial \theta }=\, &\sum\limits_{\pmb h}P({\pmb h}\mid {\pmb v}) \frac{\partial (-energy({\pmb v}, {\pmb h}) )}{\partial \theta}-\nonumber\\ & \sum\limits_{\pmb v}\sum\limits_{\pmb h}P({\pmb v}, {\pmb h})\frac{\partial (-energy({\pmb v}, {\pmb h})) }{\partial \theta} \end{align} $
(8) 由以上的表达式中, 我们可以将前一部分定义为激励部分, 表示为$v$节点的期望值表示; 后一部分作为抑制部分, 表示在联合概率下的期望表示.
1.3 影响因素诊断
在完成上述数据降维与非线性特征提取之后, 对影响因素进行诊断.其中, 本文采用偏导数表示该变量对${T_N}$变化的贡献率大小, 因此对于某一个数据$X_0$, 它对的贡献率[21] $K$为:
$ \begin{equation} K=\frac{\partial {T_N}}{\partial {X}}\mid _{{X=X_0}} \end{equation} $
(9) 1.4 改进的深度置信网络建模流程
基于非线性特征提取的深度置信网络的建立步骤为:
1) 通过机理确定变量.
2) 对输入数据进行移除异常值以及零均值归一化.
3) 设计IT网络以及深度置信网的网络结构.
4) 选择数据训练IT网络, 得到非线性PCA降维模型, 并计算各变量的统计量, 作为影响因素诊断依据.
5) 将降维数据输入深度置信网训练网络.
6) 用检测数据对整个模型进行检验.如果效果不满意, 则返回3).
整体分类器模型结构如图 3所示.
2. 深度置信网络复杂度分析
为深刻揭示本文所提改进型DBN在网络结构和算法复杂度方面的优势, 本节从如下两个方面进行分析:
1) 网络结构复杂度
信息熵的概念是1958年香农借鉴热力学上分子混乱程度来描述信息源含信息量的不确定度.从信息学的角度出发, 可以论证所提方法在优化网络结构上的优势, 采用隐含层的信息熵来体现网络的结构性和组织性[22].武妍等在论述提高网络泛化能力优化网络结构中提出通过正则化(惩罚函数)的方法, 来控制网络的"有效复杂度"[23]. Deco等通过构建基于互信息熵的正则函数, 来等效网络的"有效复杂度", 并进行网络结构优化.其中输入层和隐含层之间的互信息熵[24]定义为:
$ \begin{eqnarray*} H=-\sum\limits_{j=1}^Qc_j{\rm log}c_j+\frac{l}{P}\sum\limits_{l=1}^Q\sum\limits_{j=1}^Qc_{jl}{\rm log}c_{jl} \end{eqnarray*} $
(10) 其中, $P$为输入样本数, $Q$为隐含层节点数, $c_{jl}$为第$l$个样本对第$j$个隐含单元的归一化输出, $c_j$为平均值.熵的单位取决于定义用到对数的底, 当底数为2, 熵的单位是bit; 当底数为e, 熵的单位是nat; 而当底数为10, 熵的单位是Hart.
定理1. 面向具有相同特征的样本数据设计的两个训练深度网络net1和net2, 若网络"有效复杂度"相同($H_{\rm net1}$=$H_{\rm net2}$), 当网络的输入层节点$P_{\rm net1}$ $<$ $P_{\rm net2}$时, 则有, 网络的隐含层节点总和$Q_{\rm net1}$ $<$ $Q_{\rm net2}$.
证明. 假设原DBN网络(net1)的互信息熵函数已是最小化, 其中第一层RBM完全反映了输入层和隐含层的互信息.根据信息熵原理, 则有[25]:
$ \begin{eqnarray*} -\sum\limits_{j=1}^Qc_j{\rm log}c_j=-\sum\limits_{j=1}^{Q_l}\frac{l}{Q_l}{\rm log}\frac{l}{Q_l}={\rm log}Q_l \end{eqnarray*} $
(11) $Q_l$代表隐含层第一层的节点数.将式(11)代入式(10)中可得:
$ \begin{equation} H={\rm log}Q_l+\frac{l}{P}\sum\limits_{l=1}^{P}\sum\limits_j^{Q_l}c_{jl}{\rm log}c_{jl} \end{equation} $
(12) 基于DBN原理, 本文提出的改进型DBN网络(net2)应使每一个RBM都能完全重构输入变量, 因此, 也应使所有互信息熵最小化, 则有改进方法后的互信息熵为$H'$:
$ \begin{equation} H'={\rm log}Q'_l+\frac{l}{P'}\sum\limits_{l=1}^{P'}\sum\limits_{j'}^{Q'_l}c_{j'l}{\rm log}c_{j'l} \end{equation} $
(13) 又因为, 如完全重构原始输入变量, (由于假设NPCA完全提取了原来样本数据中的特征信息, 因此, net2中第一层RBM依然为求解隐含层节点到原始样本信息的映射关系), 则必有:
$ \begin{equation} \sum\limits_{j'}^{Q'_l}c_{j'l}{\rm log}c_{j'l}=\sum\limits_j^{Q_l}c_{jl}{\rm log}c_{jl} \end{equation} $
(14) 此外, 因为同样满足互信息熵最小化, (对于同一样本数据, 我们采用同种DBN网络结构进行信息映射时, "有效复杂度"应该相等.也就是正则函数相等), 即$H=H'$, 因此当$P'\leq P$时, 则必有$Q'_l\leq Q_l$.同理, 后续隐含层之间的RBM节点个数同样具有此规律.因此可得, $Q_{\rm net2}=Q'_1+Q'_2+\cdots+Q'_n\leq Q_{\rm net1}=Q_1+Q_2+\cdots+Q_n$ ($n$为网络的隐含层总层数).综上可以得到改进后的网络总节点存在$S_{\rm net1}<S_{\rm net2}$.
2) 算法复杂度分析
算法的复杂度就是对算法计算所需要的时间和空间的一种度量[25].一般将算法的复杂度分为时间复杂度和空间复杂度.时间复杂度是以算法结构主体执行循环次数为依据, 空间复杂度以程序主体占据空间为依据[26].一个算法中的语句执行次数称为语句频度或时间频度, 记为$T(n)$, 若有某个辅助函数$f(n)$, 使得当$n$趋近于无穷大时, $T(n)/f(n)$的极限值为不等于零的常数, 则称$f(n)$是$T(n)$的同数量级函数, 记作$T(n)={\rm O}(f(n))$, 称${\rm O}(f(n))$为算法的渐进时间复杂度, 简称时间复杂度.用O代表一个算法的计算复杂度, 算法中的循环语句是算法的主体, 若算法中含有并列的算法, 则将并列的算法复杂度相加.例如:
for $i=1:n$
$x=x+1;$
end
for $i=1:n$
for $j=1:n$
$x=x+1;$
end
end
第一个for循环的复杂度为O$(n)$, 第二个循环的复杂度为O$(n2)$, 则整个算法的复杂度为O$(n+n2)={\rm O}(n2)$.
定理2. 假设存在一个DBN网络, 其结构为含有$n$层隐含层, 隐含层节点数为$[h_1, h_2, \cdots, h_n]$.则存在一类基于ITNN神经网络的非线性特征提取机制的DBN (假设该机制可以有效保证数据的互信息熵不变), 当DBN主结构的节点数可以减小到$[h_1', h_2', \cdots, h_n']$, 预处理IT网络部分隐含层节点数为$h_0$, 且满足ITNN的隐含层节点数$h_0<\beta$, 则有:改进后DBN的时间复杂度降低小于原DBN.其中,
$ \begin{equation} \beta =\prod\limits_{1}^{n}h_n-\prod\limits_{1}^{n}h'_n \end{equation} $
(15) 证明.
1) 传统DBN算法的伪代码如下:
for (1:DBN的隐含层第一层节点数$h_1$)
for (1:DBN的隐含层第二层节点数$h_2$)
$\cdots$
for (1:DBN的隐含层第$n$层节点数$h_n$)
计算网络的目标函数是否符合要求
end
$\cdots$
end
end
所以DBN的复杂度为O(DBN)= $\prod_{1}^{n}h_n$
2) NPCA-DBN算法的伪代码如下:
for (1:NPCA的隐含层节点数$h_0$)
计算输入训练网络的目标函数是否符合要求
end
for (1:DBN的隐含层第一层节点数$h_1'$)
for (1:DBN的隐含层第二层节点数$h_2'$)
$\cdots$
for (1:DBN的隐含层第$n$层节点数$h_n'$)
计算网络的目标函数是否符合要求
end
$\cdots$
end
end
故所提算法的复杂度为O(NPCA-DBN)$=h0+\prod_{1}^{n}h'_n$若$h0<\prod_{1}^{n}h_n-\prod_{1}^{n}h'_n$, 则有: O(NPCA-DBN)-O(DBN)\, $<$\, 0, 即改进型DBN的时间复杂度降低.此外, 由于改进型DBN降低了原DBN结构中的隐含层节点数, 则有效降低算法的空间复杂度. \hfill$\square$ \section{实例研究} PM2.5预测和影响因素诊断涉及的变量众多, 而且影响变量之间多存在关联, 本节给出了PM2.5浓度预测与超标影响因素诊断方法并进行了数值验证.
3. 实例研究
PM2.5预测和影响因素诊断涉及的变量众多,而且影响变量之间多存在关联,本节给出了PM2.5浓度预测与超标影响因素诊断方法并进行了数值验证.
3.1 PM2.5的预测和影响因素诊断方法
基于第二节所提混合分类器模型, 选用相关污染物和气象因素作为输入, 提取主元非线性特征之后, 输入深度置信网, 来进行预测, 并根据统计量信息诊断PM2.5浓度超标原因.算法建立的步骤如下:
1) 选择历史数据, 并建立非线性PCA和深度置信网的模型.
2) 训练模型.
3) 检测模型效果.
4) 得出预测及诊断结果.
如图 4所示.
3.2 数据来源
为验证本文混合模型的有效性, 采取河北省某市地表水厂, 华电二区和胶片厂三个检测点于2014年11月至2015年4月间的监测数据作为实验数据.其中, 为分析检测数据, 依据文献[27]选取相关污染物如: PM10, SO$_2$, NO, CO, O$_3$, 气象数据如:风速, 风向, 温度, 湿度, 相关空气指数数:空气指数AQI.实验采用500个训练样本, 100个检测样本, 模型训练次数设置为50 000次.
3.3 PM2.5浓度预测
1) 网络结构
基于本文所提出的改进DBN模型, 利用历史数据, 进行PM2.5的浓度预测, 本文采用实验的方式获得模型的结构, 并与传统的预测DBN模型进行对比.在参考其他文献以及经验规则的基础上, 通过实验获得改进DBN的网络结构, 如图 5所示.
图 5中: xx-xx-xx为隐含层的结构, 代表DBN三层主结构中的隐含层和内部节点分配.可见试验后得到DBN主结构隐含层的节点数结构为10-6-6为本次使用的网络结构, 其中数据预处理阶段采用的浅层学习网络采用试验方法得到有一层隐含层节点, 非线性节点数为10.对比传统DBN网络结构, 两者间的对比关系如表 1所示.
表 1 网络结构对比Table 1 The comparison of the network structure模型 结构 隐含层节点数 总节点数 算法总空间复杂度 NPCA-DBN (6-10-10) + (6-10-6-6-1) 32 55 $6\times 10\times 10+6\times 10\times 6\times 6\times 1$ DBN 10-12-10-10-1 32 43 $10\times 12\times 10\times 10$ 其中(6-10-10) + (6-10-6-6-1)代表网络整体结构, 对于预处理阶段的浅层网络有6-10-10的网络结构, DBN主结构的输入层为6个节点, 隐含层为三层, 第一层是10个节点, 第二层和第三层为6个节点, 一个输出的结构, 由于改进的DBN的两部分的节点不在同一个网络嵌套中, 故为两个部分的复杂度相加.由上表我们可以清楚地看出改进的DBN模型在主结构中的深层网络中, 大大减少了非线性节点的个数, 从而在算法复杂度上实现数量级上的减小.
2) 建模精度对比实验
预测阶段采用检测输出的平均相对误差MRE (Mean relative error)来表示预测的精度.
$ \begin{equation} {\rm MRE}=\frac{\sum\limits_{j=1}^{m}\left(\frac{\left |X_{obs, j}-X_{exp, j} \right |} {X_{exp, j}}\right)}{m} \end{equation} $
(16) 其中, $m$是检测数据的样本数. $X_{obs, j}$表示检测数据的输出值, $X_{exp, j}$表示检测数据的真值.平均相对误差反映出了在预测上偏离真值的平均水平.为清晰展现本文所提DBN的优势, 以华电二区监测点为例, 图 6和图 7分别给出了改进DBN与传统的DBN、SVM和PLS在预测效果上的对比结果和建模误差趋势.
在图 6中, 横坐标为监测数据的100个采样点, 纵坐标为PM2.5的浓度.其中*代表模型输出的预测值, o代表实际值.我们可以直观地看出, 改进DBN的模型预测效果更佳, 同时比传统的分类方法在精度上有所提升.
图 7中横坐标代表 100个监测时间点, 纵坐标代表各个预测值的相对误差.由图 6和图 7可以清晰地展现出, 改进DBN模型的预测精度并没有因为降低输入的维度而降低.通过对比分析, 我们可以得到:首先, 传统DBN没有经过降维预处理, 预测精度不高; 另外, 由于DBN在处理海量数据建模时有显著的优势, 而本实验训练样本为500个, 因此, 所提改进的DBN在预测精度上与SVM提升不明显; 其次, 与ANN(BP)方法对比, 是因为BP只有一个隐含层, 属于浅层学习, 训练网络深度不足; 最后, PLS适用于处理线性模型的预测问题, 对非线性关系的建模精度欠佳.此外, 基于多种非线性特征提取机制下的复合分类预测方法, 表 2给出了该市地表水厂, 华电二区和胶片厂周边PM2.5浓度的建模精度和收敛速度对比.
表 2 建模精度与收敛速度对比Table 2 The comparison of the network structure监测点 指标 NPCA-DBN NPCA-ANN NPCA-SVM NPCA-PLS DBN ANN SVM PLS 地表 MRE ($\times10^{-2}$) 13.32 22.21 13.14 26.82 17.92 23.40 12.19 24.54 水厂 训练时间(s) 44 16 180 46 89 33 349 94 华电 MRE ($\times10^{-2}$) 14.57 25.15 13.04 29.48 17.01 24.16 10.22 27.16 二区 训练时间(s) 37 12 211 49 90 38 401 103 胶片 MRE ($\times10^{-2}$) 10.51 26.49 11.09 33.16 12.77 23.32 12.73 30.06 厂 训练时间(s) 42 16 198 57 108 42 399 108 由表 2我们可以得到, 本文所提改进DBN在建模精度和收敛速度上都有较大提升, NPCA数据预处理算法通过提取数据之间的非线性特征, 得到原有数据的非线性表达, 对于网络化表达的机器学习算法可以提高建模精度, 并降低训练时间.特别的, 同样采用NPCA数据预处理机制的复合分类方法, 对于基于线性化拟合的浅层学习算法随着训练时间的降低, 建模精度有所下降.
3.4 PM2.5浓度影响因素诊断
基于所提出改进DBN的影响因素诊断方法, 在实验中, 我们将空气质量指数, PM10浓度, SO$_2$、CO、NO$_2$、O$_3$气体浓度, 风向, 风速, 相对湿度, 温度等10个过程变量作为诊断部分的输入变量.由于该地区的特殊性, 在4 000多次的采样数据中, 有二分之一采样点数据的PM2.5浓度都高于$100\, {\rm \mu g/m^3}$, 所以, 为了展现模型对PM2.5影响因素诊断结果, 我们实验设置的PM2.5预测限为$200{\rm \mu g/m^3}$, 即处于重度污染的情况下, 计算输入变量中对于结果的贡献率[28].并用贡献图的方式表达影响因素诊断结果.我国对PM2.5浓度级别划分如表 3所示[29].
表 3 PM2.5浓度级别Table 3 The PM2.5 concentration level浓度范围(${\rm \mu g/m^3}$) 级别 优良级别 0$\, \sim\, $50 1级 优 50$\, \sim\, $100 2级 良 101$\, \sim\, $150 3级 轻度污染 151$\, \sim\, $200 4级 中度污染 201$\, \sim\, $ 5级 重度污染 在历史数据中选定所有未超限数据对应的输入, 求平均水平代表未超限数据组合作为参考输入变量集: reference (152, 168, 63, 1.55, 50, 42, 163, 2.26, 0.63, 3.88)针对华电二区监测区域, 在图 8中, 我们以超限组第20组贡献图为例说明诊断过程.
由贡献图可以看出, 第七个变量对结果的贡献最大.我们观察验证第83组数据输入变量集为: $X_{20}= $ (179, 207, 79, 1.78, 73, 177, 266, 3.07, 0.44, 7.34)第6个变量的相对偏差为最大, 因此诊断结果为:造成此次污染物浓度过高的首要原因是风速的原因.
3.5 小结
从上述实验可以看出, 本文提出的改进的DBN模型在预测效果上并没有使精度降低, 同时, 加快了模型的收敛速度.并且在超标诊断中, 平均超标检测率达到$85\, \%$, 能够有效地诊断出PM2.5浓度超标的主要因素.
4. 结论
本文提出的基于非线性特征提取的DBN模型能够有效完成含复杂非线性特征关系高维数据的预测建模诊断任务.基于信息熵理论, 证明了本文所提DBN模型相比传统DBN, 能够在不降低建模精度的同时, 达到降低网络和算法复杂度的优势, 对于深度学习理论在海量数据挖掘中的应用具有重要理论意义.将所提建模方法应用到一类PM2.5浓度预测与诊断问题中, 并与传统DBN、SVM、ANN、PLS等分类方法和含NPCA数据预处理机制的复合分类方法做了详细对比, 验证了所提方法的优势与正确性.需要说明的是, 本文采取基于数据驱动的方法对PM2.5进行浓度预测和影响因素诊断, 在PM2.5的形成机理上还未做过多的分析, 在未来的研究中将深入探讨PM2.5浓度变化机理.此外, 由于本文数据来源于特定城市的采样点, 因此在方法的适用性方面还要做深入的研究.下一步的工作将分为以下两部分进行: 1)理论方面, 面向深度置信网络结构本身的优化方法的研究, 研究自适应样本数据特征的网络模型结构. 2)应用方面, 尝试将所提方法应用到复杂流程工业的建模和诊断问题中.
-
表 1 受试者个人相关信息
Table 1 Demographic characteristics of the participants
病例 年龄 病程 性别 患侧 患病部位 痉挛肌群 ROM ($^{\circ})$ MAS Test MAS Retest Test Retest E1 E2 E1 E2 A S1 64 6月 男 右 左侧脑梗塞 屈肘肌群 96 104 0 0 0 0 S2 60 7月 女 右 左侧脑梗死 屈肘肌群 120 126 0 0 0 0 S3 66 3月 男 左 右侧脑梗死 伸肘肌群 100 104 2 2 2 2 S4 68 3月 男 左 右侧脑栓塞 屈肘肌群 132 124 1 1 1 1 S5 66 4月 女 左 右侧脑梗塞 屈肘肌群 113 114 1.5 1.5 1.5 1.5 S6 37 2月 男 左 右侧脑出血 屈肘肌群 114 117 1 1 1 1 S7 73 3月 男 右 左侧脑梗塞 屈肘肌群 105 113 0 0 0 0 S8 51 2月 男 左 右侧脑梗塞 屈肘肌群 94 99 1 1 1 1 S9 45 2月 女 左 右侧脑出血 屈肘肌群 100 101 2 2 2 2 S10 30 2月 女 左 右侧脑梗塞 屈肘肌群 106 106 1.5 1.5 1.5 1.5 S11 62 6月 男 左 右侧脑梗塞 伸肘肌群 130 135 3 3 3 3 S12 68 3月 男 左 右侧脑栓塞 伸肘肌群 123 125 3 3 3 3 S13 37 2月 男 左 右侧脑出血 伸肘肌群 114 117 2 2 2 2 S14 51 2月 男 左 右侧脑梗塞 伸肘肌群 98 100 1.5 1.5 1.5 1.5 Test Retest E3 E4 E3 E4 B S15 32 14月 男 左 右侧脑出血 屈肘肌群 128 129 1.5 1.5 1.5 1.5 S16 37 15月 男 右 左侧脑梗塞 屈肘肌群 118 120 2 2 2 2 S17 39 2年 男 右 左侧脑出血 屈肘肌群 103 106 0 0 0 0 S18 21 2年 男 右 左脑挫裂伤 屈肘肌群 110 110 3 3 3 3 S19 41 15天 男 左 右侧脑出血 屈肘肌群 101 100 1.5 1.5 1.5 1.5 S20 30 9月 男 左 右侧脑出血 屈肘肌群 110 115 1 1 1 1 S21 46 1年 男 右 左侧脑出血 屈肘肌群 143 145 1.5 1.5 1.5 1.5 S22 54 1月 男 右 左侧脑梗死 屈肘肌群 132 136 1 1 1 1 A:安徽省立医院, B:安徽中医药大学第一附属医院; ROM: range of motion; E1 $\sim$ E4: Evaluator1 $\sim$ 4, 检查者1 $\sim$ 4; S1 $\sim$ S22: Subject1 $\sim$ 22, 受试者1 $\sim$ 22; MAS Test:第一次MAS评分, MAS Retest:三天后MAS评分. 表 2 第一次以及三天后测试结果
Table 2 Experimental result at test and at retest
Test Retest 受试者 $\theta_{ p}$ $\theta_{ p}$/$\theta_{\rm rom}$ AMV $\theta_{ p}$ $\theta_{ p}$/$\theta_{\rm rom}$ AMV $\theta_{ p}$ $\theta_{ p}$/$\theta_{\rm rom}$ AMV $\theta_{ p}$ $\theta_{ p}$/$\theta_{\rm rom}$ AMV A E1 E2 E1 E2 S1 96 1.00 0 96 1.00 0 104 1.00 0 104 1.00 0 S2 120 1.00 0 120 1.00 0 126 1.00 0 126 1.00 0 S3 38 0.38 0.698 34 0.34 0.423 50 0.48 0.596 46 0.44 0.624 S4 64 0.48 0.412 93 0.70 0.392 70 0.56 0.321 80.63 0.65 0.469 S5 64 0.57 0.910 50 0.44 0.862 58 0.51 0.867 68 0.60 0.453 S6 75 0.66 0.192 64 0.56 0.182 62 0.53 0.260 75 0.64 0.175 S7 105 1.00 0 105 1.00 0 113 1.00 0 113 1.00 0 S8 71 0.76 0.526 56 0.60 0.421 54 0.55 1.501 47 0.47 0.479 S9 36 0.36 0.636 30 0.30 0.584 27 0.27 0.673 23 0.23 0.702 S10 50 0.47 0.498 44 0.42 0.455 53 0.50 0.426 55 0.52 0.452 S11 27 0.21 0.565 33 0.25 0.627 22 0.16 1.421 40 0.30 0.526 S12 23 0.19 1.268 25 0.20 0.746 46 0.37 0.837 29 0.23 1.105 S13 22 0.19 0.412 30 0.26 0.536 33 0.28 0.459 40 0.34 0.381 S14 42 0.43 0.526 46 0.47 0.469 49 0.49 0.450 55 0.55 0.502 B E3 E4 E3 E4 S15 54 0.42 0.380 80 0.63 0.485 64 0.50 0.325 65 0.50 0.416 S16 57 0.38 0.440 48 0.41 0.648 68 0.57 0.574 51 0.43 0.464 S17 103 1.00 0 103 1.00 0 106 1.00 0 106 1.00 0 S18 38 0.35 1.172 28 0.25 0.947 29 0.26 0.804 37 0.34 2.828 S19 42 0.42 0.206 50 0.50 0.258 55 0.55 0.249 37 0.37 0.225 S20 63 0.57 0.364 60 0.55 0.426 74 0.64 0.378 68 0.59 0.301 S21 60 0.42 0.628 70 0.49 0.450 53 0.37 0.459 50 0.34 0.416 S22 74 0.56 0.488 87 0.66 0.598 66 0.49 0.450 77 0.57 0.459 注: $\theta_{p}$($^{\circ}$); AMV(m/s$^{2}$). 表 3 $\theta_{ p}$/$\theta_{\rm rom}$及Amv与MAS的相关性
Table 3 The correlation between the MAS scores and $\theta_{ p}$/$\theta _{\rm rom}$ and between the MAS scores and AMV
相关性分析结果 装置数据 MAS评分 相关系数 $P$ Test Retest Test Retest E1 $\theta_{ p}$/$\theta_{\rm rom}$ $-$0.944 $-$0.918 $P<$ 0.05 $P<$ 0.05 AMV 0.821 0.665 $P<$ 0.05 $P<$ 0.05 E2 $\theta_{ p}$/$\theta_{\rm rom}$ $-$0.953 $-$0.931 $P<$ 0.05 $P<$ 0.05 AMV 0.841 0.864 $P<$ 0.05 $P<$ 0.05 E3 $\theta_{ p}$/$\theta_{\rm rom}$ $-$0.855 $-$0.846 $P<$ 0.05 $P<$ 0.05 AMV 0.857 0.900 $P<$ 0.05 $P<$ 0.05 E4 $\theta_{ p}$/$\theta_{\rm rom}$ $-$0.940 $-$0.831 $P<$ 0.05 $P<$ 0.05 AMV 0.873 0.813 $P<$ 0.05 $P<$ 0.05 表 4 $\theta_{ p}$/$\theta_{\rm rom}$及加速度重测相关系数
Table 4 The test-retest reliability of $\theta_{ p}$/$\theta_{\rm rom}$ and AMV
相关性分析结果 重测相关参数 相关系数 $P$值 E1 $\theta_{ p}$/$\theta_{\rm rom}$ 0.938 $P<$ 0.05 AMV 0.632 $P<$ 0.05 E2 $\theta_{ p}$/$\theta_{\rm rom}$ 0.962 $P<$ 0.05 AMV 0.824 $P<$ 0.05 E3 $\theta_{ p}$/$\theta_{\rm rom}$ 0.890 $P<$ 0.05 AMV 0.928 $P<$ 0.05 E4 $\theta_{ p}$/$\theta_{\rm rom}$ 0.912 $P<$ 0.05 AMV 0.796 $P<$ 0.05 表 5 装置结果验证
Table 5 Result validation of device
受试者 年龄
(岁)患病部位 患侧 装置结果 装置预测MAS MAS评分 $\theta_{ p}$/$\theta_{\rm rom}$ AMV $\sqrt{(\Delta \theta)^{2}}$ $\sqrt{\left({\Delta a} \right)^{2}}$ $\sqrt{(\Delta \theta)^{2}+\left({\Delta a} \right)^{2}}$ S1 57 右侧脑出血 左侧 0.65 0.320 1 1 1 1 S2 46 左侧脑干梗塞 右侧 0.49 0.435 1.5 2 1.5 1.5 -
[1] Bar-On L, Molenaers G, Aertbeliën E, Monari D, Feys H, Desloovere K. The relation between spasticity and muscle behavior during the swing phase of gait in children with cerebral palsy. Research in Developmental Disabilities, 2014, 35(12):3354-3364 doi: 10.1016/j.ridd.2014.07.053 [2] 王广志.肌肉痉挛定量评估的研究进展.现代康复, 2000, 4(5):650-652 http://news.medlive.cn/ns/info-progress/show-131231_181.htmlWang Guang-Zhi. The advancing of quantitative assessment of spasticity. Modern Rehabilitation, 2000, 4(5):650-652 http://news.medlive.cn/ns/info-progress/show-131231_181.html [3] Feldman R G, Young R R, Koella W P. Spasticity: Disordered Motor Control. Chicago: Year Book Medical Publishers, 1980. 485-495 [4] Ishikawa S, Okamoto S, Isogai K, Akiyama Y, Yanagihara N, Yamada Y. Wearable dummy to simulate joint impairment: severity-based assessment of simulated spasticity of knee joint. In: Proceedings of the 2013 IEEE/SICE International Symposium on System Integration. Kobe, Japan: IEEE, 2013. 300-305 [5] Gracies J M, Brashear A, Jech R, McAllister P, Banach M, Valkovic P, Walker H, Marciniak C, Deltombe T, Skoromets A, Khatkova S, Edgley S, Gul F, Catus F, De Fer Beatrice B, Vilain C, Picaut P. Safety and efficacy of abobotulinumtoxinA for hemiparesis in adults with upper limb spasticity after stroke or traumatic brain injury:a double-blind randomised controlled trial. The Lancet Neurology, 2015, 14(10):992-1001 doi: 10.1016/S1474-4422(15)00216-1 [6] Kim J-H, Park H-S, Damiano D L. Accuracy and reliability of haptic spasticity assessment using HESS (Haptic Elbow Spasticity Simulator). In: Proceedings of the 33rd Annual International Conference of the IEEE EMBS. Boston, USA: IEEE, 2001. 8527-8530 [7] 李雪萍, 程凯, 周俊, 杨婷, 于俊龙, 陈安亮, 张红飞.表面肌电联合等速测试评定肌痉挛的临床研究.中国现代医学杂志, 2010, 20(4):605-608, 611 http://www.docin.com/p-1694078342.htmlLi Xue-Ping, Cheng Kai, Zhou Jun, Yang Ting, Yu Jun-Long, Chen An-Liang, Zhang Hong-Fei. Clinical assessment of muscle spasticity in stroke patients with surface electromyography and isokinetic assessment. China Journal of Modern Medicine, 2010, 20(4):605-608, 611 http://www.docin.com/p-1694078342.html [8] Klingels K, De Cock P, Molenaers G, Desloovere K, Huenaerts C, Jaspers E, Feys H. Upper limb motor and sensory impairments in children with hemiplegic cerebral palsy. Can they be measured reliably? Disability and Rehabilitation, 2010, 32(5):409-416 doi: 10.3109/09638280903171469?needAccess=true [9] Fleuren J F M, Voerman G E, Erren-Wolters C V, Snoek G J, Rietman J S, Hermens H J, Nene A V. Stop using the Ashworth Scale for the assessment of spasticity. Journal of Neurology, Neurosurgery & Psychiatry, 2010, 81(1):46-52 http://jnnp.bmj.com/content/jnnp/81/1/46.full.pdf [10] Mehrholz J, Wagner K, Meissner D, Grundmann K, Zange C, Koch R, Pohl M. Reliability of the modified Tardieu scale and the modified Ashworth scale in adult patients with severe brain injury:a comparison study. Clinical Rehabilitation, 2015, 19(7):751-759 https://www.researchgate.net/publication/7516824_Reliability_of_the_Modified_Tardieu_Scale_and_the_Modified_Ashworth_Scale_in_adult_patients_with_severe_brain_injury_A_comparison_study [11] Kaya T, Goksel K A, Gunaydin R, Koc A, Altundal E U. Inter-rater reliability of the modified Ashworth scale and modified modified Ashworth scale in assessing poststroke elbow flexor spasticity. International Journal of Rehabilitation Research, 2011, 34(1):59-64 doi: 10.1097/MRR.0b013e32833d6cdf [12] Kim K S, Seo J H, Song C G. Portable measurement system for the objective evaluation of the spasticity of hemiplegic patients based on the tonic stretch reflex threshold. Medical Engineering & Physics, 2011, 33(1):62-69 https://www.researchgate.net/publication/47382015_Portable_measurement_system_for_the_objective_evaluation_of_the_spasticity_of_hemiplegic_patients_based_on_the_tonic_stretch_reflex_threshold [13] Abdollahi I, Taghizadeh A, Shakeri H, Eivazi M, Jaberzadeh S. The relationship between isokinetic muscle strength and spasticity in the lower limbs of stroke patients. Journal of Bodywork and Movement Therapies, 2014, 19(2):284-290 https://www.sciencedirect.com/science/article/pii/S1360859214001296 [14] Wu Y N, Park H-S, Ren Y P, Deborah G-S, Chen J J, Zhang L Q. Measurement of elbow spasticity in stroke patients using a manual spasticity evaluator. In: Proceedings of the 28th IEEE EMBS Annual International Conference. New York, USA: IEEE, 2006. 3974-3977 [15] 宋凡, 张峰, 朱玉连, 沈建, 白玉龙, 李放, 胡永善, 吴毅.等速测试指标与改良Ashworth法用于评定肌痉挛的相关性研究.中国康复医学杂志, 2008, 23(7):615-617 http://www.doc88.com/p-9425621312124.htmlSong Fan, Zhang Feng, Zhu Yu-Lian, Shen Jian, Bai Yu-Long, Li Fang, Hu Yong-Shan, Wu Yi. Research of the correlation between the isokinetic parameters and the modified Ashworth scores in spasticity measurement. Chinese Journal of Rehabilitation Medicine, 2008, 23(7):615-617 http://www.doc88.com/p-9425621312124.html [16] Fujisawa T, Takagi M, Takahashi Y, Inoue K, Terada T, Kawakami Y, Komeda T. Basic research on the upper limb patient simulator. In: Proceedings of the 10th IEEE International Conference on Rehabilitation Robotics. Noordwijk, The Netherlands: IEEE, 2007. 48-51 [17] Park H-S, Kim J-H, Damiano D L. Development of a haptic elbow spasticity simulator (HESS) for improving accuracy and reliability of clinical assessment of spasticity. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2012, 20(3):361-370 doi: 10.1109/TNSRE.2012.2195330 [18] Takhashi Y, Komeda T, Koyama H, Yamamoto S, Arimatsu T, Kawakami Y, Inoue K, Ito Y. Development of an upper limb patient simulator for physical therapy exercise. In: Proceedings of the 2011 IEEE International Conference on Rehabilitation Robotics. Zurich, Switzerland: IEEE, 2011. 1117-1120 [19] 李林, 纪仲秋, 李艳霞, 刘威彤.关节角度重置法、运动最小阈值测量法和力量重现法的相关性研究.天津体育学院学报, 2016, 31(1):36-40 http://c.wanfangdata.com.cn/periodical/tjtyxyxb/2016-1.aspxLi Lin, Ji Zhong-Qiu, Li Yan-Xia, Liu Wei-Tong. Correlation between joint position sense reproduce test, the threshold of detection of movement test and force sense reproduce test. Journal of TUS, 2016, 31(1):36-40 http://c.wanfangdata.com.cn/periodical/tjtyxyxb/2016-1.aspx [20] Calota A, Feldman A G, Levin M F. Spasticity measurement based on tonic stretch reflex threshold in stroke using a portable device. Clinical Neurophysiology, 2008, 119(10):2329-2337 doi: 10.1016/j.clinph.2008.07.215 [21] Pandyan A D, Price C I M, Rodgers H, Barnes M P, Johnson G R. Biomechanical examination of a commonly used measure of spasticity. Clinical Biomechanics, 2001, 16(10):859-865 doi: 10.1016/S0268-0033(01)00084-5 [22] Pandyan A D, Price C I M, Barnes M P, Johnson G R. A biomechanical investigation into the validity of the modified Ashworth scale as a measure of elbow spasticity. Clinical Rehabilitation, 2003, 17(3):290-294 doi: 10.1191/0269215503cr610oa [23] McCrea P H, Eng J J, Hodgson A J. Linear spring-damper model of the hypertonic elbow:reliability and validity. Journal of Neuroscience Methods, 2003, 128(1-2):121-128 doi: 10.1016/S0165-0270(03)00169-9 [24] Grow D I, Wu M N, Locastro M J, Arora S K, Bastian A J, Okamura A M. Haptic simulation of elbow joint spasticity. In:Proceedings of the 2008 Symposium on Haptic Interfaces for Virtual Environment and Teleoperator Systems. Nevada, USA:IEEE, 2008. 475-476 [25] 韩京清.从PID技术到"自抗扰控制"技术.控制工程, 2002, 9(3):13-18 http://mall.cnki.net/magazine/Article/JZDF200203004.htmHan Jing-Qing. From PID technique to active disturbances rejection control technique. Control Engineering of China, 2002, 9(3):13-18 http://mall.cnki.net/magazine/Article/JZDF200203004.htm [26] 韩建达, 谈大龙, 蒋新松.直接驱动机器人关节加速度反馈解耦控制.自动化学报, 2000, 26(3):289-295 http://www.aas.net.cn/CN/abstract/abstract16059.shtmlHan Jian-Da, Tan Da-Long, Jiang Xin-Song. Joint acceleration feedback control for direct-drive robot decoupling. Acta Automatica Sinica, 2000, 26(3):289-295 http://www.aas.net.cn/CN/abstract/abstract16059.shtml 期刊类型引用(6)
1. 王创,陈文杰,陈伟海,孙先涛,林岩. 适应人体重心起伏的悬吊减重康复系统设计. 北京航空航天大学学报. 2024(07): 2322-2330 . 百度学术
2. 徐丹,黎华茂,杨慧,卫星,刘曦明,丁然,黄礼群,许萌. 悬吊运动疗法在青年军人膝关节前交叉韧带重建同期半月板缝合术后康复中的作用. 创伤外科杂志. 2023(08): 579-585 . 百度学术
3. 高海波,牛福亮,刘振,于海涛,李楠. 悬吊式微低重力环境模拟技术研究现状与展望. 航空学报. 2021(01): 80-99 . 百度学术
4. 于宁波,邹武林. 有限频域约束下串联弹性驱动器的刚度控制. 控制理论与应用. 2019(05): 711-719 . 百度学术
5. 张从鹏,李小龙,解毅. 动态减重康复机器人控制系统开发. 机械工程师. 2019(09): 14-16 . 百度学术
6. 申琳. 国际平衡能力研究前沿热点分析与演进. 中国学校体育(高等教育). 2018(01): 56-60 . 百度学术
其他类型引用(7)
-