CN109645995A - 基于肌电模型和无迹卡尔曼滤波的关节运动估计方法 - Google Patents
基于肌电模型和无迹卡尔曼滤波的关节运动估计方法 Download PDFInfo
- Publication number
- CN109645995A CN109645995A CN201910038177.3A CN201910038177A CN109645995A CN 109645995 A CN109645995 A CN 109645995A CN 201910038177 A CN201910038177 A CN 201910038177A CN 109645995 A CN109645995 A CN 109645995A
- Authority
- CN
- China
- Prior art keywords
- joint
- myoelectricity
- model
- parameter
- kalman filtering
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1071—Measuring physical dimensions, e.g. size of the entire body or parts thereof measuring angles, e.g. using goniometers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1118—Determining activity level
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4528—Joints
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Physiology (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Rheumatology (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及到一种基于肌电模型和无迹卡尔曼滤波的关节运动估计方法,首先采集膝关节在连续运动状态下股二头肌、股四头肌、股外侧肌、股内侧肌、半腱肌、股薄肌的肌电信号和实时角度,对其进行带通滤波处理,并提取小波系数和均方根特征,然后使用一种结合了肌肉动力学、关节动力学、骨骼动力学和相关肌电特征的状态空间肌电模型,通过无迹卡尔曼滤波算法,得出Sigma采样集χi和权重Wi,然后进行进一步的预测,计算出系统状态变量和协方差矩阵P(k+1|k),迭代循环后,实现对膝关节连续运动的估计。该方法与传统的角度估计方法相比,减小了系统误差、累积误差和外部干扰的影响,精度高,稳定性好,对目标机动反应快速,有了明显的改进。
Description
技术领域
本发明属于模式识别领域,涉及一种肌电信号模式识别方法,特别涉及一种基于状态空间肌电模型和无迹卡尔曼滤波的关节连续运动估计方法。
背景技术
表面肌电信号(surface Electromyography,sEMG)是现今比较热门的前沿科学技术人机交互的输入信号源,它是一种非平稳的微弱信号,是肌肉兴奋时与其相关的运动单元共同产生的一组动作电位序列,它的特征区分明显,蕴含的信息丰富,且采集简单无创伤,成为目前人机交互技术中热门的研究领域。对于表面肌电信号的研究主要集中在特征提取和模式识别这两个过程。相应的研究成果也比较成熟,可以识别多个离散的动作类别。但在康复医疗机器人等领域更多时候需要的是对患者连续运动变量的预测,以实现对康复机器人的平滑柔顺控制。
传统的关节连续运动估计方法有提取肌电特征,然后通过神经网络建立sEMG与连续运动的回归模型。生理肌肉模型也是估计连续关节运动的一种方法。Buchanan等人提出了一种基于肌电信号的正向动力学模型,该模型由希尔肌肉模型(Hill muscle model,HMM)、肌肉激活动力学和关节正向动力学组成。该模型涉及多个生理参数,计算困难,实际应用受到限制。HMM是估计连续关节运动最常用的肌肉模型,但存在两个问题:一是HMM涉及到许多难以识别的复杂的生理参数,计算负担也比较大;二是HMM可以直接从sEMG信号计算关节力矩,但如果需要连续的关节运动估计,则还需要从力矩中计算运动状态。这通常会带来累计误差,降低预测精度。
而将HMM与关节正向动力学结合起来,并对模型参数进行简化替代的方法能够有效解决上述问题,该方法不需要计算关节力矩,而是可以直接从sEMG信号中计算关节运动。同时使用肌电特征组成测量方程作为反馈,并使用闭环预测算法,可较为准确的估计关节连续运动。
发明内容
本发明涉及到一种用于关节角度估计的状态空间肌电模型和无迹卡尔曼滤波方法,首先采集膝关节在连续运动状态下股二头肌、股四头肌、股外侧肌、股内侧肌、半腱肌、股薄肌的肌电信号和实时角度,对其进行带通滤波处理,并提取小波系数和均方根特征,然后使用一种结合了肌肉动力学、关节动力学、骨骼动力学和相关肌电特征的状态空间肌电模型,通过无迹卡尔曼滤波算法,得出Sigma采样集χi和权重Wi,然后进行进一步的预测,计算出系统状态变量和协方差矩阵P(k+1|k),迭代循环后,实现对膝关节连续运动的估计。
为了实现以上目的,本发明方法主要包括以下步骤:
步骤一,采集关节连续运动时相关肌肉的肌电信号,即通过肌电信号采集仪采集关节运动时相关肌肉的肌电信号,然后采用带通滤波的方法对原始信号进行预处理。
步骤二,根据希尔肌肉模型和关节动力学求出状态空间肌电模型的非线性表达式;该状态空间肌电模型首先对希尔肌肉模型进行参数代替及简化处理,简化处理后的离散时间的预测模型如下:
Ts是采样时间,是k时刻的关节角速度,θk是k时刻的关节位置,si是代替参数,均为常数。
然后提取均方根Xrms和小波系数αj,k组成测量方程作为状态反馈。然后按下式将肌电特征与关节运动拟合。
u取值1和2,是离线识别的固定参数,和是时间k的均方根和小波系数。
得到最终表达式:
其中ak=a(k),ωk是过程噪声,υk是测量噪声,T是采样时间,是关节角加速度,是关节角速度,θk是关节位置,是离线识别的固定参数,和是时间k的均方根和小波系数,si是均为常数的代替参数;
步骤三,根据步骤二的状态空间肌电模型,使用无迹卡尔曼滤波算法对膝关节连续运动进行估计。Sigma采样集χi和权重Wi定义如下:
其中χi是Sigma采样集,Wi是相应权重,是特征状态的均值,n是特征状态维数,P(k)是误差协方差矩阵,是调节参数,并假设ωk和υk均为高斯白噪声。
步骤四,对步骤三的Sigma采样集进行进一步预测,并计算出系统状态变量和协方差矩阵如下:
其中是特征状态变量,χi(k+1|k)是Sigma采样集,P(k+1|k)是误差协方差矩阵。
然后再次运用无迹变换,得出系统残差和卡尔曼增益矩阵如下:
其中Sk+1是系统残差,Kk+1是卡尔曼增益矩阵,是进一步预测的特征状态变量,γi(k+1|k)是进一步观测的采样集,R是噪声协方差矩阵。
步骤五:令k=k+1,对步骤四进行迭代循环,最终实现对膝关节连续运动的估计。
本发明设计的基于状态空间肌电模型和无迹卡尔曼滤波的关节连续运动估计,具有如下特点:
本发明建立的用于膝关节角度估计的状态空间肌电模型和无迹卡尔曼滤波方法,将正向动力学与希尔肌肉模型结合了起来,并对其进行了参数简化,使其可直接用于膝关节运动的估计,减小了累积误差。同时提取均方根和小波系数等肌电特征,建立测量方程,减小了系统误差和外部干扰,提高了关节预测的精度。使用的闭环预测算法,无迹卡尔曼滤波算法,精度高,稳定性好,对目标机动反应快速。该方法与传统的角度估计方法相比,在预测精度方面有了明显的改进。
附图说明
图1为本发明的实施流程图;
图2(a)为本发明有负重状态采集的肌电信号特征;
图2(b)为本发明无负重状态采集的肌电信号特征;
图3为使用发明的预测模型在无负重状态下的估计效果图;
图4为使用发明的预测模型在有负重状态下的估计效果图。
具体实施方式
如图1所示,本实施例包括如下步骤:
步骤一,采集膝关节连续运动时的肌电信号,具体是:四名志愿者坐在椅子上分别在有负重和无负重情况下进行膝关节屈伸运动,动作时长为10秒,通过Trigno肌电信号采集仪采集膝关节运动时相关肌肉的肌电信号,分别是股二头肌,股四头肌,股外侧肌,股内侧肌,半腱肌,股薄肌,然后采用带通滤波方法进行预处理。
步骤二,根据希尔肌肉模型和关节动力学求出状态空间肌电模型的非线性表达式,该状态空间肌电模型首先对希尔肌肉模型进行参数代替及简化处理,再提取均方根和小波系数肌电特征,组成测量方程作为状态反馈,最终与关节运动进行拟合,得到状态空间肌电模型的非线性表达式。
参数替代和模型结合后的加速度计算如下:
si是代替参数,均为常数。
对滤波处理后的肌电信号进行特征提取,提取均方根Xrms和小波系数αj,k如图2所示,组成测量方程作为状态反馈。然后按下式将肌电特征与关节运动拟合。
u取值1和2,是离线识别的固定参数,和是时间k的均方根和小波系数。
参数识别如表1所示:
表1有无负重情况下的参数
最终得到状态空间肌电模型的非线性表达式:
其中ak=a(k),ωk是过程噪声,υk是测量噪声,T是采样时间,是关节角加速度,是关节角速度,θk是关节位置,是离线识别的固定参数,和是时间k的均方根和小波系数,si是均为常数的代替参数。
步骤三,根据步骤二的状态空间肌电模型,使用无迹卡尔曼滤波算法对膝关节连续运动进行估计。Sigma采样集χi和权重Wi定义如下:
其中χi是Sigma采样集,Wi是相应权重,是特征状态的均值,n是特征状态维数,P(k)是误差协方差矩阵,是调节参数,并假设ωk和υk均为高斯白噪声。
步骤四,对步骤三的Sigma采样集进行进一步预测,并计算出系统状态变量和协方差矩阵如下:
其中是特征状态变量,χi(k+1|k)是Sigma采样集,P(k+1|k)是误差协方差矩阵。
然后再次运用无迹变换,得出系统残差和卡尔曼增益矩阵如下:
其中Sk+1是系统残差,Kk+1是卡尔曼增益矩阵,是进一步预测的特征状态变量,γi(k+1|k)是进一步观测的采样集,R是噪声协方差矩阵。
步骤五:令k=k+1,对步骤四进行迭代循环,最终实现对膝关节连续运动的估计,结果如图3和图4所示。
Claims (1)
1.基于肌电模型和无迹卡尔曼滤波的关节运动估计方法,其特征在于,该方法包括如下步骤:
步骤一,采集关节连续运动时相关肌肉的肌电信号,即通过肌电信号采集仪采集关节运动时相关肌肉的肌电信号,然后采用带通滤波的方法对原始信号进行预处理;
步骤二,根据希尔肌肉模型和关节动力学求出状态空间肌电模型的非线性表达式;该状态空间肌电模型首先对希尔肌肉模型进行参数代替及简化处理,简化处理后的离散时间的预测模型如下:
Ts是采样时间,是k时刻的关节角速度,θk是k时刻的关节位置,si是代替参数,均为常数;
然后提取均方根Xrms和小波系数αj,k组成测量方程作为状态反馈;然后按下式将肌电特征与关节运动拟合;
u取值1和2,是离线识别的固定参数,和是时间k的均方根和小波系数;
得到最终表达式:
其中ak=a(k),ωk是过程噪声,υk是测量噪声,T是采样时间,是关节角加速度,是关节角速度,θk是关节位置,是离线识别的固定参数,和是时间k的均方根和小波系数,si是均为常数的代替参数;
步骤三,根据步骤二的状态空间肌电模型,使用无迹卡尔曼滤波算法对膝关节连续运动进行估计;Sigma采样集χi和权重Wi定义如下:
其中χi是Sigma采样集,Wi是相应权重,是特征状态的均值,n是特征状态维数,P(k)是误差协方差矩阵,是调节参数,并假设ωk和υk均为高斯白噪声;
步骤四,对Sigma采样集进行进一步预测,并计算出系统状态变量和协方差矩阵如下:
其中是特征状态变量,χi(k+1|k)是Sigma采样集,P(k+1|k)是误差协方差矩阵;
然后再次运用无迹变换,得出系统残差和卡尔曼增益矩阵如下:
其中Sk+1是系统残差,Kk+1是卡尔曼增益矩阵,是进一步预测的特征状态变量,γi(k+1|k)是进一步观测的采样集,R是噪声协方差矩阵;
步骤五:令k=k+1,对步骤四进行迭代循环,最终完成对膝关节连续运动的估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910038177.3A CN109645995B (zh) | 2019-01-16 | 2019-01-16 | 基于肌电模型和无迹卡尔曼滤波的关节运动估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910038177.3A CN109645995B (zh) | 2019-01-16 | 2019-01-16 | 基于肌电模型和无迹卡尔曼滤波的关节运动估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109645995A true CN109645995A (zh) | 2019-04-19 |
CN109645995B CN109645995B (zh) | 2021-09-07 |
Family
ID=66120114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910038177.3A Active CN109645995B (zh) | 2019-01-16 | 2019-01-16 | 基于肌电模型和无迹卡尔曼滤波的关节运动估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109645995B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111258426A (zh) * | 2020-01-17 | 2020-06-09 | 杭州电子科技大学 | 基于肌电肌音模型和无迹粒子滤波的关节运动估计方法 |
CN111820902A (zh) * | 2020-06-29 | 2020-10-27 | 北京科技大学 | 一种基于活动度特征的踝关节韧带损伤智能决策系统 |
CN114371232A (zh) * | 2021-12-22 | 2022-04-19 | 天津国科医工科技发展有限公司 | 基于卡尔曼滤波算法的色谱滤波方法、装置、介质、系统 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060004299A1 (en) * | 2004-06-30 | 2006-01-05 | Honda Motor Co., Ltd. | Motion measurement method, motion measurement system, and motion measurement program |
US20080009771A1 (en) * | 2006-03-29 | 2008-01-10 | Joel Perry | Exoskeleton |
CN102764167A (zh) * | 2012-06-12 | 2012-11-07 | 天津大学 | 基于相关系数的肌电假肢控制源导联优化方法 |
CN103054585A (zh) * | 2013-01-21 | 2013-04-24 | 杭州电子科技大学 | 基于生物运动信息的上肢肩肘腕关节运动功能评价方法 |
CN103976739A (zh) * | 2014-05-04 | 2014-08-13 | 宁波麦思电子科技有限公司 | 穿戴式摔倒动态实时检测方法和装置 |
US20140277271A1 (en) * | 2013-03-14 | 2014-09-18 | The Chinese University Of Hong Kong | Device and methods for preventing knee sprain injuries |
CN106456000A (zh) * | 2014-05-30 | 2017-02-22 | 微软技术许可有限责任公司 | 生物测定信号的基于运动的估计 |
CN107622260A (zh) * | 2017-10-26 | 2018-01-23 | 杭州电子科技大学 | 基于多源生物信号的下肢步态相识别方法 |
CN107703756A (zh) * | 2017-11-03 | 2018-02-16 | 广州视源电子科技股份有限公司 | 动力学模型参数辨识方法、装置、计算机设备及存储介质 |
CN109086247A (zh) * | 2018-09-19 | 2018-12-25 | 合肥工业大学 | 基于双时间尺度无迹卡尔曼滤波的系统故障参数估计方法 |
CN109084772A (zh) * | 2018-07-25 | 2018-12-25 | 北京航天长征飞行器研究所 | 一种基于无迹卡尔曼的视线转率提取方法及系统 |
-
2019
- 2019-01-16 CN CN201910038177.3A patent/CN109645995B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060004299A1 (en) * | 2004-06-30 | 2006-01-05 | Honda Motor Co., Ltd. | Motion measurement method, motion measurement system, and motion measurement program |
US20080009771A1 (en) * | 2006-03-29 | 2008-01-10 | Joel Perry | Exoskeleton |
CN102764167A (zh) * | 2012-06-12 | 2012-11-07 | 天津大学 | 基于相关系数的肌电假肢控制源导联优化方法 |
CN103054585A (zh) * | 2013-01-21 | 2013-04-24 | 杭州电子科技大学 | 基于生物运动信息的上肢肩肘腕关节运动功能评价方法 |
US20140277271A1 (en) * | 2013-03-14 | 2014-09-18 | The Chinese University Of Hong Kong | Device and methods for preventing knee sprain injuries |
CN103976739A (zh) * | 2014-05-04 | 2014-08-13 | 宁波麦思电子科技有限公司 | 穿戴式摔倒动态实时检测方法和装置 |
CN106456000A (zh) * | 2014-05-30 | 2017-02-22 | 微软技术许可有限责任公司 | 生物测定信号的基于运动的估计 |
CN107622260A (zh) * | 2017-10-26 | 2018-01-23 | 杭州电子科技大学 | 基于多源生物信号的下肢步态相识别方法 |
CN107703756A (zh) * | 2017-11-03 | 2018-02-16 | 广州视源电子科技股份有限公司 | 动力学模型参数辨识方法、装置、计算机设备及存储介质 |
CN109084772A (zh) * | 2018-07-25 | 2018-12-25 | 北京航天长征飞行器研究所 | 一种基于无迹卡尔曼的视线转率提取方法及系统 |
CN109086247A (zh) * | 2018-09-19 | 2018-12-25 | 合肥工业大学 | 基于双时间尺度无迹卡尔曼滤波的系统故障参数估计方法 |
Non-Patent Citations (1)
Title |
---|
QI SHAO ET AL: "An EMG-driven model to eatimate muscle forces and joint moments in stroke patients", 《COMPUTERS IN BIOLOGY AND MEDICINE》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111258426A (zh) * | 2020-01-17 | 2020-06-09 | 杭州电子科技大学 | 基于肌电肌音模型和无迹粒子滤波的关节运动估计方法 |
CN111258426B (zh) * | 2020-01-17 | 2023-03-21 | 杭州电子科技大学 | 基于肌电肌音模型和无迹粒子滤波的关节运动估计方法 |
CN111820902A (zh) * | 2020-06-29 | 2020-10-27 | 北京科技大学 | 一种基于活动度特征的踝关节韧带损伤智能决策系统 |
CN114371232A (zh) * | 2021-12-22 | 2022-04-19 | 天津国科医工科技发展有限公司 | 基于卡尔曼滤波算法的色谱滤波方法、装置、介质、系统 |
CN114371232B (zh) * | 2021-12-22 | 2024-03-22 | 天津国科医工科技发展有限公司 | 基于卡尔曼滤波算法的色谱滤波方法、装置、介质、系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109645995B (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | sEMG-based joint force control for an upper-limb power-assist exoskeleton robot | |
Phinyomark et al. | Application of wavelet analysis in EMG feature extraction for pattern classification | |
CN109262618B (zh) | 基于肌肉协同的上肢多关节同步比例肌电控制方法与系统 | |
CN108743223B (zh) | 一种基于健侧生物电控的智能外骨骼康复机械手及方法 | |
CN109645995A (zh) | 基于肌电模型和无迹卡尔曼滤波的关节运动估计方法 | |
Chen et al. | Cross-comparison of EMG-to-force methods for multi-DoF finger force prediction using one-DoF training | |
Phinyomark et al. | Wavelet-based denoising algorithm for robust EMG pattern recognition | |
CN109657651A (zh) | 一种基于肌电信号的下肢膝关节连续运动估计方法 | |
JP7477309B2 (ja) | 生体信号が表す情報を識別するためのシステム | |
JP2010125287A (ja) | 指関節角度推定装置 | |
Wang et al. | Prediction of lower limb joint angle using sEMG based on GA-GRNN | |
CN113990441A (zh) | 一种基于生物动力学的下肢膝关节主动肌肌电拟合方法 | |
CN109758336A (zh) | 一种利用患者患侧sEMG进行主动康复治疗的控制方法 | |
CN112617858A (zh) | 一种多尺度补偿传递熵的皮层肌肉功能耦合方法 | |
Khorasani et al. | Brain control of an external device by extracting the highest force-related contents of local field potentials in freely moving rats | |
Zhuang et al. | Motion control of a four-wheel-independent-drive electric vehicle by motor imagery EEG based BCI system | |
Zhang et al. | sEMG-based shoulder-elbow composite motion pattern recognition and control methods for upper limb rehabilitation robot | |
Wang et al. | Lower limb motion recognition based on surface electromyography signals and its experimental verification on a novel multi-posture lower limb rehabilitation robots☆ | |
Li et al. | Prediction of knee joint moment by surface electromyography of the antagonistic and agonistic muscle pairs | |
CN105769186B (zh) | 基于复杂网络的肌电信号采集位置选取方法 | |
CN116831874A (zh) | 一种基于肌电信号的下肢康复器控制方法 | |
Herle | Movement intention detection from SEMG signals using time-domain features and discriminant analysis classifiers | |
Yu et al. | The research of sEMG movement pattern classification based on multiple fused wavelet function | |
Zhou et al. | sEMG-driven functional electrical stimulation tuning via muscle force | |
CN110151176A (zh) | 一种基于肌电信号的上肢肘关节连续运动估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190419 Assignee: HANGZHOU DUKANG TECHNOLOGY CO.,LTD. Assignor: HANGZHOU DIANZI University Contract record no.: X2022330000025 Denomination of invention: Joint motion estimation method based on EMG model and unscented Kalman filter Granted publication date: 20210907 License type: Common License Record date: 20220128 |