CN103886215B - Walking ability analyzing method and device based on muscle collaboration - Google Patents

Walking ability analyzing method and device based on muscle collaboration Download PDF

Info

Publication number
CN103886215B
CN103886215B CN201410136557.8A CN201410136557A CN103886215B CN 103886215 B CN103886215 B CN 103886215B CN 201410136557 A CN201410136557 A CN 201410136557A CN 103886215 B CN103886215 B CN 103886215B
Authority
CN
China
Prior art keywords
gait
muscle
walking ability
cooperative mode
data
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.)
Expired - Fee Related
Application number
CN201410136557.8A
Other languages
Chinese (zh)
Other versions
CN103886215A (en
Inventor
刘玲
李飞
陈香
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201410136557.8A priority Critical patent/CN103886215B/en
Publication of CN103886215A publication Critical patent/CN103886215A/en
Application granted granted Critical
Publication of CN103886215B publication Critical patent/CN103886215B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

The invention provides a walking ability analyzing method and device based on muscle collaboration. The walking ability analyzing device based on muscle collaboration comprises a multi-channel surface myoelectricity acquisition module, a data analysis and storage module and a gait data analyzing and processing module. According to the walking ability calculating method based on muscle collaboration, the non-negative matrix factorization algorithm is adopted to extract gait muscle collaboration from a gait cycle surface myoelectricity signal, an adult gait muscle collaboration mode serves as a standard template, and the muscle collaboration mode of a subject is graded by means of a Synergy Comprehensive Calculation (SCC) template based on gait muscle collaboration mode similarity, so that characteristic parameters reflecting the gait of the subject are obtained. According to the walking ability analyzing method based on muscle collaboration, the difference of the walking abilities of different people is disclosed from the aspect of brain central nervous system control, and an objective walking ability parameter can be obtained. The objective walking ability parameter reflects the effect of control of the central nervous system over the gait.

Description

A kind of walking ability collaborative based on muscle analyzes method and apparatus
Technical field
The present invention is gait analysis technology based on sensor technology and bio signal treatment technology, refers in particular to a kind of based on flesh The walking ability analysis method and apparatus that meat is collaborative, belongs to sensor technology and bio signal processing technology field.
Background technology
Gait analysis has become one of analysis indispensable means of human motion system.Three-dimensional gait analysis technology at present There is two large problems: first gait analysis system uses dedicated video equipment and computer equipment, expensive, operates professional By force, it is difficult to universal.
Surface electromyogram signal (Surface Electromyography, SEMG) is obtained from body surface by surface electrode Electromyographic signal, gait data based on surface electromyogram signal obtains and need not special experimental site, can overcome gait analysis The shortcoming that video capture device is easily subject to environmental disturbances;And surface electromyogram signal can not only provide the essential information of gait, Can also extract muscle from surface electromyogram signal to work in coordination with, calculating for gait provides objective effective parameter.Therefore it provides it is a kind of Analyze method and apparatus based on the walking ability that muscle is collaborative to have great importance.
Summary of the invention
The deficiency existed for current gait analysis technology, it is contemplated that propose a kind of walking energy collaborative based on muscle Power analyzes method and apparatus.
A kind of walking ability collaborative based on muscle is analyzed method and is included:
(1) utilize Non-negative Matrix Factorization (Non-negative Matrix Factorization, NMF) algorithm from gait Period surface electromyographic signal (Surface Electromyography, SEMG) carries the gait muscle Cooperative Mode of step both legs. Refer specifically to NMF that (gait cycle Matrix of envelope is pretreatment and gait cycle to the gait cycle Matrix of envelope of experimenter's both legs Gait cycle surface electromyogram signal after segmentation) carry out Non-negative Matrix Factorization to extract left and right leg muscle Harmonious Matrix L, R.Will L, R are expressed as:
L=Lm×l={ l1,l2…ll}
R=Rm×r={ r1,r2…rr}
Wherein column vector li,rj(i=1,2 ... l;J=1,2 ... r) representing the muscle Cooperative Mode of left and right lower limb respectively, l, r divide The gait muscle of left and right lower limb Biao Shi not work in coordination with quantity, m represents the quantity of one group of muscle selected in gait analysis;Use adult The muscle Cooperative Mode Criterion template set T that left and right lower limb is total:
T=Tm×s={ t1,t2…ts}
Wherein column vector tk(k=1,2 ... s) represent the flesh that standard gait muscle Cooperative Mode, i.e. adult's both legs are total Meat Cooperative Mode;S represents the muscle Cooperative Mode number that the quantity of standard gait muscle Cooperative Mode, i.e. adult left and right lower limb are total Amount;
(2) coordination with the synthesis based on gait muscle Cooperative Mode similarity is utilized to calculate (Synergy Comprehensive Calculation, SCC) calculating of model realization experimenter's walking ability index;SCC model considers experimenter's gait flesh Three factors of the symmetry that the quantity of meat Cooperative Mode, structure and left and right lower limb show in walking;And assist based on muscle Same walking ability analyzes method, and the angle controlled from nervous system discloses the brain coordination control strategy to gait motion, The result of calculation of this model can inherently reflect the brain coordination control ability to gait motion.
Described SCC model calculates walking ability and refers to calibration method, including step calculated below:
(1) by Pearson correlation coefficients computing formula R, (x y) calculates experimenter left and right each gait muscle Cooperative Mode of lower limb li,rj(i=1,2 ... l;J=1,2 ... r) with each standard gait muscle Cooperative Mode tk(k=1,2 ... coefficient R between s) (tk,li),R(tk,rj).The similarity degree Sim (T, L), Sim of definition left and right leg muscle Harmonious Matrix L, R and standard form collection T (T, R), its computing formula is as follows:
S i m ( T , L ) = Σ i = 1 l F ( T , l i )
S i m ( T , R ) = Σ j = 1 r F ( T , r j )
F(T,li)=max (R (tk,li)) k=1,2 ... s
F(T,rj)=max (R (tk,rj)) k=1,2 ... s
F(T,li), F (T, rj) represent left and right leg muscle Cooperative Mode l respectivelyi,rjSimilarity degree with standard form collection T. To each F (T, li), F (T, rj) it is multiplied by score coefficient B=100/s, then left and right lower limb similarity score Sl,SrIt is respectively as follows:
Sl=B*Sim (T, L)
Sr=B*Sim (T, R)
Sl,SrValue between 0-100.
(2) can reflect that left and right lower limb is in gait due to the correlation coefficient between the lower limb gait muscle Cooperative Mode of experimenter left and right The symmetry shown in motion, therefore in symmetry score Sym defined in SCC model.Use Pearson's phase in the model (x y) calculates experimenter left and right lower limb each gait muscle Cooperative Mode l to close coefficient formulas Ri,rj(i=1,2 ... l;J=1, 2 ... the coefficient R (r between r)j,li) and define left and right leg muscle Harmonious Matrix L, the similarity degree Sim (R, L) between R, Its computing formula is:
S i m ( R , L ) = &Sigma; i = 1 l F ( R , l i ) r &GreaterEqual; l &Sigma; j = 1 r F ( L , r j ) r < l
F(R,li)=max (R (rj,li)) j=1,2 ... r
F(L,rj)=max (R (rj,li)) i=1,2 ... l
F(R,li) represent left leg muscle Cooperative Mode l in this SCC modeliSimilarity degree with right leg muscle Harmonious Matrix R; F(L,rj) then represent right leg muscle Cooperative Mode rjSimilarity degree with left leg muscle Harmonious Matrix L.The meter of symmetry score Sym Calculation formula is as follows:
Sym=A*Sim (R, L)
Wherein (l, r), the value of Sym is between 0-100 for symmetry score coefficient A=100/min.
Finally, S is soughtl,Sr, the meansigma methods of Sym, obtain overall scores S:
S=(Sl+Sr+Sym)/3
S is the result of calculation of SCC model, and value is between 0-100.
A kind of device realizing the above-mentioned walking ability analysis method collaborative based on muscle, employing techniques below scheme: This device is made up of multi-channel surface myoelectric acquisition module, data parsing memory module, gait data analysis and processing module.Wherein Multi-channel surface myoelectric acquisition module is for gathering the SEMG of both legs in gait motion, and will believe by the way of wired or wireless Number it is sent to data parsing memory module;Data parsing memory module resolves gait data in real time, and stores;Gait number It is analyzed processing according to analysis and processing module gait data of load store from data parsing memory module, calculates walking ability Index.
Described gait analysis processing module include gait data pretreatment unit 1., gait cycle cutting unit 2., gait Muscle work in coordination with extraction unit 3., walking ability indicator calculating unit 4., display control unit 5.;The most 4. refer to as walking ability The core cell that mark calculates;Display control unit is 5. for arranging the running parameter of system, duty, and result Display etc..
The most 4. described walking ability index calculates is four serial computing unit, and wherein gait data pretreatment unit is 1. Process for the SEMG of every passage is filtered, goes average and rectification, to eliminate baseline noise in gait data acquisition etc. Impact, finally exports the envelope signal of SEMG.2. gait cycle cutting unit uses gait cycle based on acceleration signal to divide Cut algorithm and the envelope of SEMG after pretreatment is carried out the segmentation of gait cycle;Passage SEMG every to all gait cycles after segmentation Amplitude is normalized, poor to eliminate the signal amplitude between the different experimenters caused by factors such as electrode placement positions Different, finally by down-sampled for each passage of all gait cycle SEMG signals being partitioned into for uniform length, to be walked State cycle Matrix of envelope.Gait muscle is worked in coordination with extraction unit and 3. the gait cycle Matrix of envelope being partitioned into is carried out nonnegative matrix and divide Solve, calculate left and right lower limb gait muscle Cooperative Mode and the quantity of muscle Cooperative Mode thereof.Walking ability indicator calculating unit is the most then Use the calculating of SCC model realization walking ability index.
Described data parsing memory module and gait data analysis and processing module can be at common computer or band WIFI Smart mobile phone, the mobile terminal such as panel computer realizes, and therefore this walking ability analytical equipment is not limited by use occasion, can In hospital, rehabilitation institution, the environment such as family uses.
Accompanying drawing explanation
Fig. 1 is the structural representation of apparatus of the present invention;
Fig. 2 is surface myoelectric electrode setting method schematic diagram in the embodiment of the present invention;
Fig. 3 is the general frame calculating walking ability index in the present invention;
Fig. 4 is the schematic diagram of gait cycle signal segmentation in the embodiment of the present invention;
Fig. 5 is SCC model calculation flow chart in the embodiment of the present invention.
Detailed description of the invention
Illustrating below in conjunction with the accompanying drawings, providing one embodiment of the present of invention, concrete implementation process is as follows:
1., according to step 1 shown in Fig. 3, in this case study on implementation, we are to choose left and right lower limb tibialis anterior (TA), flatfish (LG), vastus lateralis (VL) outside flesh (SO), gastrocnemius, rectus femoris (RF), semitendinosus m. (SE), caput longum musculi bicipitis femoris (BF) and wealthy Tensor fascia femoris (TF) is as a example by totally 16 pieces of muscle are as the acquisition position of gait SEMG signal.Therefore the multi-channel surface shown in Fig. 1 Myoelectricity acquisition module (1) should include at least 16 channel surface electromyographic electrodes (each 8 passages of left and right lower limb).As in figure 2 it is shown, this embodiment Electrode used therein uses two strip electrode rods with skin contact, and during placement, two electrode connecting lines are consistent with muscle fiber direction of travel, Shave off fine hair at electrode placement positions in advance and use alcohol swab wiping.Use in this embodiment based on acceleration signal simultaneously The gait cycle segmentation of ACC, therefore multi-channel surface myoelectric acquisition module is also integrated with 2 passage 3-axis acceleration sensors, real It is individually positioned in below the lower limb knee of experimenter left and right at tibia during executing.Experimenter holds according to the speed of oneself the most comfortable Row walking task, in implementation process, experimenter need to move ahead in level ground N number of gait cycle, during acquisition module synchronous acquisition walking SEMG and ACC signal and send by the way of wired to the data parsing memory module (2) shown in Fig. 1.
2. the data parsing memory module (2) of device realizes on a common computer, first according to the concrete frame of data Form carries out real-time parsing, then parsed data is stored in file.
3. gait data analysis and processing module (3) realizes on a common computer.First its display control unit starts literary composition Part is loaded into function, loads pending gait data, has loaded when file has loaded interface prompt, and display control unit opens Dynamic data analysis, detailed process is as follows:
(1) original gait data initially enters pretreatment unit, carries out gait signal according to step 2 as shown in Figure 3 and locates in advance Reason: in this case study on implementation, first use cut-off frequency be 40Hz, the FIR high pass filter on the 50 rank gait data to experimenter Be filtered, then go average, rectification, finally to use cut-off frequency be 10Hz, the FIR low pass filter on 50 rank again filter with Obtain the gait SEMG envelope as shown in B in Fig. 4.
(2) gait SEMG envelope data enters gait cycle cutting unit, utilizes ACC signal to split it: such as C in Fig. 4 Shown in, in walking, acceleration signal ACC is quasi-periodic time varying signal, utilizes based on window function the most in this embodiment The peak point of peak-value detection method detection left and right lower limb acceleration transducer Z axis signal, is respectively mapped to left and right lower limb by peak point Gait SEMG envelope, is partitioned into N number of gait cycle envelope signal.B, C, D in Fig. 4 illustrates and utilizes gait ACC signal to split The detailed process of gait cycle.Finally by unified for every passage of N number of gait cycle envelope signal down-sampled for length n=500, and This passage is normalized by the maximum utilizing the every passage of gait cycle, obtains for extracting the gait week that muscle is collaborative Phase Matrix of envelope, as shown in the E in Fig. 4, the corresponding muscle passage of every a line of matrix.N number of step by lower limb single after pretreatment State cycle Matrix of envelope is designated as EMGo={ Oi, i=1,2 ... N}, OiBeing the matrix of a m × n, m represents selected in gait analysis The quantity of the one group of muscle taken, m=8 in this embodiment.
(3) gait muscle works in coordination with extraction unit to gait cycle Matrix of envelope EMGoProcess, the step 3 of Fig. 3 know, Gait muscle works in coordination with extraction unit can extract the muscle Cooperative Mode of left and right lower limb respectively, and calculates the quantity of Cooperative Mode, and it is detailed Thin process is as follows:
1) to each gait cycle matrix Oi, according to Non-negative Matrix Factorization NMF algorithm, calculating muscle and working in coordination with quantity is j (j=1,2 ... m) corresponding muscle Harmonious Matrix Wi m×j,
2) gait muscle Cooperative Mode quantity s is sought: utilize Wi m×j,Trying to achieve muscle Cooperative Mode quantity is reconstruction during j Gait pattern matrix stackWherein rebuild gait pattern matrixThen institute is calculated There is EMGjWith EMG before decompositionoSquare error matrix stack V={Vj, j=1,2 ... m}, wherein VjComputing formula is:
V j = 1 N &Sigma; i = 1 N ( 1 - ( O i - R i j ) 2 O i 2 ) 0 < V j < 1
VjRepresent that muscle works in coordination with reconstruction gait pattern matrix stack EMG when quantity is jjWith original gait cycle matrix stack EMGo Mean Square Error.According to priori, gait muscle Cooperative Mode quantity s is defined as V under t inspection (p < 0.05)j's Value works in coordination with quantity j more than the minimum muscle corresponding to 0.95, and the muscle Harmonious Matrix of experimenter is:
W m &times; s = { w 1 , w 2 ... w s } = 1 N &Sigma; i = 1 N W i m &times; s
3) experimenter left and right leg muscle Harmonious Matrix is expressed as:
L=Lm×l={ l1,l2…ll}
R=Rm×r={ r1,r2…rr}
Wherein column vector li,rj(i=1,2 ... l;J=1,2 ... r) representing the muscle Cooperative Mode of left and right lower limb respectively, l, r divide The gait muscle of left and right lower limb Biao Shi not work in coordination with quantity.Extract in advance about several adult according to method 1-3 described in this embodiment The muscle Cooperative Mode Criterion template set T that lower limb is total:
T=Tm×s={ t1,t2…ts}
Wherein column vector tk(k=1,2 ... s) represent standard gait muscle Cooperative Mode, s represents that standard gait muscle is worked in coordination with The quantity of pattern, and standard form is stored in Data Analysis Services module.
(4) walking ability indicator calculating unit is according to SCC model, calculates walking ability according to flow chart as shown in Figure 5: First gait muscle is worked in coordination with result of calculation L of extraction unit, R, l, r and standard form collection T and standard gait muscle works in coordination with mould Formula quantity s inputs this unit;Calculate left and right lower limb respective similarity score S the most respectivelyl,Sr;According to the comparative result of l, r, meter Calculate left and right lower limb symmetry score Sym in walking;Last comprehensive Sl,Sr, Sym obtains walking ability index.Its concrete calculating Process is as follows:
1) left and right lower limb similarity score S is calculatedl,SrShi Xianyong Pearson correlation coefficients computing formula R (x, y) calculate tested Person left and right lower limb each gait muscle Cooperative Mode li,rj(i=1,2 ... l;J=1,2 ... r) work in coordination with mould with each standard gait muscle Formula tk(k=1,2 ... coefficient R (t between s)k,li),R(tk,rj), try to achieve l by equation belowi,rjWith standard form collection T Similarity degree F (T, li), F (T, rj):
F(T,li)=max (R (tk,li)) k=1,2 ... s
F(T,rj)=max (R (tk,rj)) k=1,2 ... s
R ( x , y ) = &Sigma; i = 1 s i z e ( x ) ( x i - x &OverBar; ) ( y i - y &OverBar; ) &Sigma; i = 1 s i z e ( x ) ( x i - x &OverBar; ) 2 &Sigma; i = 1 s i z e ( x ) ( y i - y &OverBar; ) 2
Then to each F (T, li), F (T, rj) being multiplied by score coefficient B=100/s, then left and right lower limb is similar to standard form Property score Sl,SrIt is respectively as follows:
S l = B &Sigma; i = 1 l F ( T , l i )
S r = B &Sigma; j = 1 r F ( T , r j )
WhereinRepresent the similarity degree Sim (T, L), Sim of L, R and T in SCC model respectively (T,R);Sl,SrValue between 0-100.
2), when calculating left and right lower limb symmetry score, first by Pearson correlation coefficients computing formula R, (x y) calculates experimenter left Right lower limb each gait muscle Cooperative Mode li,rj(i=1,2 ... l;J=1,2 ... the coefficient R (r between r)j,li), when r >= L, left and right leg muscle Harmonious Matrix L, the similarity degree Sim (R, L) between R be:
S i m ( R , L ) = &Sigma; i = 1 l F ( R , l i )
F(R,li)=max (R (rj,li)) j=1,2 ... r
When r < during l, left and right leg muscle Harmonious Matrix L, the similarity degree Sim (R, L) between R be:
S i m ( R , L ) = &Sigma; j = 1 l F ( L , r j )
F(L,rj)=max (R (rj,li)) i=1,2 ... j
Sim (R, L) and symmetry score coefficient A=100/min (l, r) is multiplied, then symmetry score Sym is:
Sym=A*Sim (R, L)
Last walking ability indicator calculating unit calculates Sl,Sr, the meansigma methods of Sym, obtain overall scores S:
S=(Sl+Sr+Sym)/3
The walking ability index S calculated by the method can reflect the difference of walking ability between Different Individual.

Claims (9)

1. analyze method based on the walking ability that muscle is collaborative for one kind, it is characterised in that this walking ability analyze method first with Non-negative Matrix Factorization (Non-negative Matrix Factorization, NMF) algorithm is from gait cycle surface electromyogram signal (Surface Electromyography, SEMG) extracts the gait muscle Cooperative Mode of both legs;Then utilize based on gait The coordination with the synthesis of muscle Cooperative Mode similarity calculates (Synergy Comprehensive Calculation, SCC) model meter Calculate the walking ability of Different Individual;Wherein SCC model is using the stable gait muscle Cooperative Mode of adult as standard gait Muscle Cooperative Mode, the comprehensive similarity analyzing experimenter left and right leg muscle Cooperative Mode and standard gait muscle Cooperative Mode, And the symmetry that left and right lower limb shows in gait motion, calculate the characteristic parameter of reflection walking ability;
Described walking ability is analyzed method and is specifically included following steps:
(1) utilize Non-negative Matrix Factorization NMF algorithm from gait cycle surface electromyogram signal (Surface Electromyography, SEMG) extracting the gait muscle Cooperative Mode of both legs, the left and right leg muscle Harmonious Matrix L that will extract in, R is expressed as:
L=Lm×l={ l1,l2…ll}
R=Rm×r={ r1,r2…rr}
Wherein column vector li,rjThe muscle Cooperative Mode of expression left and right lower limb respectively, described i=1,2 ... l;J=1,2 ... r, l, r divide The gait muscle of left and right lower limb Biao Shi not work in coordination with quantity, m represents the quantity of one group of muscle selected in gait analysis;By adult The muscle Cooperative Mode that both legs have is as standard form collection T:
T=Tm×s={ t1,t2…ts}
T is the matrix of a m × s, column vector tkExpression standard gait muscle Cooperative Mode, described k=1,2 ... s, s represent standard The quantity of gait muscle Cooperative Mode;
(2) coordination with the synthesis based on gait muscle Cooperative Mode similarity is utilized to calculate (Synergy Comprehensive Calculation, SCC) model, calculating walking ability index S:
1) by Pearson correlation coefficients computing formula R, (x y) calculates experimenter left and right lower limb each gait muscle Cooperative Mode li,rj With each standard gait muscle Cooperative Mode tkBetween coefficient R (tk,li),R(tk,rj);Definition left and right leg muscle works in coordination with square Similarity degree Sim (T, L), the Sim (T, R) of battle array L, R and standard form collection T, its computing formula is as follows:
S i m ( T , L ) = &Sigma; i = 1 l F ( T , l i )
S i m ( T , R ) = &Sigma; j = 1 r F ( T , r j )
F(T,li)=max (R (tk,li)) k=1,2 ... s
F(T,rj)=max (R (tk,rj)) k=1,2 ... s
F(T,li), F (T, rj) represent left and right leg muscle Cooperative Mode l respectivelyi,rjWith the similarity degree of standard form collection T, to often Individual F (T, li), F (T, rj) it is multiplied by score coefficient B=100/s, then left and right lower limb similarity score Sl,SrIt is respectively as follows:
Sl=B*Sim (T, L)
Sr=B*Sim (T, R)
Sl,SrValue between 0-100;
2) in SCC model, by Pearson correlation coefficients computing formula R, (x y) calculates experimenter left and right each gait muscle of lower limb association With pattern li,rjBetween coefficient R (rj,li) and define left and right leg muscle Harmonious Matrix L, the similarity degree Sim between R (R, L), its computing formula is:
S i m ( R , L ) = &Sigma; i = 1 l F ( R , l i ) r &GreaterEqual; l &Sigma; j = 1 r F ( L , r j ) r < l
F(R,li)=max (R (rj,li)) j=1,2 ... r
F(L,rj)=max (R (rj,li)) i=1,2 ... l
F(R,li) represent left leg muscle Cooperative Mode l in this SCC modeliSimilarity degree with right leg muscle Harmonious Matrix R;F(L, rj) then represent right leg muscle Cooperative Mode rjWith the similarity degree of left leg muscle Harmonious Matrix L, the calculating of symmetry score Sym is public Formula is as follows:
Sym=A*Sim (R, L)
(l, r), the value of Sym is between 0-100 for the coefficient A=100/min of wherein said symmetry score;
3) last, seek Sl,Sr, the meansigma methods of Sym, obtain overall scores S:
S=(Sl+Sr+Sym)/3
S is the result of calculation of SCC model, and the value of S is between 0-100.
2., for performing a walking ability analytical equipment collaborative based on muscle for method described in claim 1, its feature exists In, described analytical equipment includes multi-channel surface myoelectric acquisition module (1), data parsing memory module (2), gait data analysis Processing module (3);Wherein multi-channel surface myoelectric acquisition module (1) communicates with data parsing memory module (2);Gait data divides The gait data of storage in data parsing memory module (2) is analyzed processing by analysis processing module (3), calculates walking ability and refers to Mark.
3. walking ability analytical equipment as claimed in claim 2, it is characterised in that described multi-channel surface myoelectric acquisition module (1) for gathering the surface electromyogram signal of both legs in gait, owing to gait motion relates to polylith lower limb muscles, module should at least be wrapped Surface myoelectric electrode containing 16 passages;The data of multi-channel surface myoelectric acquisition module are sent to by the way of wired or wireless Data parsing memory module (2);Data parsing memory module (2) completes the parsing storage of gait data in real time.
4. walking ability analytical equipment as claimed in claim 2, it is characterised in that described gait data analysis and processing module (3) include four serial computing unit: gait data pretreatment unit 1., gait cycle cutting unit 2., gait muscle works in coordination with Extraction unit 3., walking ability indicator calculating unit 4.;The most 4. the core cell calculated as walking ability index;Additionally Including display control unit 5., it is for arranging the running parameter of system, duty, and the display of result.
5. walking ability analytical equipment as claimed in claim 4, it is characterised in that gait data pretreatment unit is 1. from data Resolve and memory module (2) load original gait SEMG signal, and complete the smooth of every passage SEMG signal and filter preprocessing, To eliminate the impact of baseline noise in gait data acquisition, output is the envelope signal of SEMG.
6. walking ability analytical equipment as claimed in claim 4, it is characterised in that gait cycle cutting unit 2. use based on The gait cycle partitioning algorithm of acceleration signal carries out the segmentation of gait cycle to the envelope signal of SEMG after pretreatment;After segmentation The amplitude of passage SEMG every to all gait cycles is normalized, and causes not by electrode placement positions factor to eliminate With the signal amplitude difference between experimenter;Finally by down-sampled for all gait cycle signals being partitioned into for uniform length, defeated Go out gait cycle Matrix of envelope.
7. walking ability analytical equipment as claimed in claim 4, it is characterised in that gait muscle is worked in coordination with extraction unit and 3. used Both legs gait cycle Matrix of envelope is decomposed by Non-negative Matrix Factorization NMF algorithm, to extract both legs gait muscle Cooperative Mode, And export the muscle Harmonious Matrix L of left and right lower limb, Cooperative Mode quantity l of R and correspondence, r;It it is critical only that and utilizes NMF algorithm to carry The Harmonious Matrix taken rebuilds gait pattern, meets constraints when rebuilding the gait pattern error amount with original gait, then tries to achieve Quantity that gait muscle is collaborative and gait muscle Cooperative Mode.
8. walking ability analytical equipment as claimed in claim 4, it is characterised in that walking ability indicator calculating unit is first Load standard form collection T, then with SCC computation model, walking ability index is calculated, export result of calculation.
9. walking ability analytical equipment as claimed in claim 2, it is characterised in that described data parsing memory module (2) and Gait data analysis and processing module (3), at common computer or the smart mobile phone of band WIFI, panel computer realizes.
CN201410136557.8A 2014-04-04 2014-04-04 Walking ability analyzing method and device based on muscle collaboration Expired - Fee Related CN103886215B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410136557.8A CN103886215B (en) 2014-04-04 2014-04-04 Walking ability analyzing method and device based on muscle collaboration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410136557.8A CN103886215B (en) 2014-04-04 2014-04-04 Walking ability analyzing method and device based on muscle collaboration

Publications (2)

Publication Number Publication Date
CN103886215A CN103886215A (en) 2014-06-25
CN103886215B true CN103886215B (en) 2017-01-11

Family

ID=50955105

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410136557.8A Expired - Fee Related CN103886215B (en) 2014-04-04 2014-04-04 Walking ability analyzing method and device based on muscle collaboration

Country Status (1)

Country Link
CN (1) CN103886215B (en)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105457164B (en) * 2016-01-06 2018-05-18 电子科技大学 Multi-channel function electrical stimulation method and system under a kind of muscle Cooperative Mode
CN105769186B (en) * 2016-03-23 2018-05-29 河北工业大学 Electromyographic signal collection position choosing method based on complex network
CN106067178B (en) * 2016-05-30 2018-11-23 杭州电子科技大学 A kind of continuous estimation method of hand joint movement based on muscle synergistic activation model
WO2020112986A1 (en) 2018-11-27 2020-06-04 Facebook Technologies, Inc. Methods and apparatus for autocalibration of a wearable electrode sensor system
US11635736B2 (en) * 2017-10-19 2023-04-25 Meta Platforms Technologies, Llc Systems and methods for identifying biological structures associated with neuromuscular source signals
CN107808694B (en) * 2016-09-08 2023-04-14 深圳先进技术研究院 System and method for judging walking ability of cerebral palsy child patient based on GMFCS
CN106334267A (en) * 2016-10-31 2017-01-18 徐钧 Neuromuscular electrical stimulation system with supplementary external skeleton
CN106474623A (en) * 2016-10-31 2017-03-08 何晓阔 Bionical multi-channel nerve muscle electric stimulation method and electric stimulation
CN106492345A (en) * 2016-10-31 2017-03-15 徐钧 Hemiplegia upper and lower extremities synergy movement electrical stimulation method and electric stimulation
CN106422172B (en) * 2016-11-22 2018-10-30 西安交通大学 Lower limb rehabilitation training system treadmill speed adaptive control method based on myoelectricity
JP6887274B2 (en) * 2017-03-21 2021-06-16 パナソニック株式会社 Walking support robot and walking support method
CN106991290B (en) * 2017-04-12 2020-03-31 中国科学技术大学 Gait ability assessment method and device
JP6590416B2 (en) 2017-07-10 2019-10-16 トヨタ自動車株式会社 Rehabilitation evaluation apparatus, operation method of rehabilitation evaluation apparatus, and rehabilitation evaluation program
US10772787B2 (en) * 2017-10-20 2020-09-15 Panasonic Intellectual Property Management Co., Ltd. Assistance apparatus, assistance method, and recording medium
US11961494B1 (en) 2019-03-29 2024-04-16 Meta Platforms Technologies, Llc Electromagnetic interference reduction in extended reality environments
US11907423B2 (en) 2019-11-25 2024-02-20 Meta Platforms Technologies, Llc Systems and methods for contextualized interactions with an environment
SG10201800971RA (en) * 2018-02-05 2019-09-27 Nec Asia Pacific Pte Ltd Method and system for motor function rehabilitation and monitoring a patient’s recovery
CN109620651B (en) * 2018-11-16 2020-03-31 中国科学技术大学 Intelligent auxiliary rehabilitation equipment based on synchronous brain and muscle electricity
CN109394230A (en) * 2018-12-07 2019-03-01 北京大学深圳研究生院 A kind of gait analysis system based on surface electromyogram signal
CN109262618B (en) * 2018-12-12 2022-04-01 武汉理工大学 Muscle cooperation-based upper limb multi-joint synchronous proportional myoelectric control method and system
CN110472595B (en) * 2019-08-20 2022-03-11 郑州大学 Electroencephalogram recognition model construction method and device and recognition method and device
CN110917577B (en) * 2019-11-27 2021-04-20 西安交通大学 Multi-stage lower limb training system and method utilizing muscle synergistic effect
CN111326252A (en) * 2020-02-18 2020-06-23 杭州电子科技大学 Method for analyzing motor function based on inter-muscle cooperation and inter-muscle coupling
CN111326232A (en) * 2020-03-18 2020-06-23 福建中医药大学 Method for identifying abnormal gait form based on three-dimensional gait analysis system
CN112025682B (en) * 2020-09-02 2022-04-19 中国科学技术大学 Flexible exoskeleton robot control method and system based on muscle cooperation theory

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6002960A (en) * 1996-08-27 1999-12-14 The Johns Hopkins University Passive, non-invasive method to quantify objectively the level and density of a neural blockade
CN201760063U (en) * 2010-04-08 2011-03-16 上海体育学院 Athletic training analysis system based on digital court and electromechanical signal
CN102210925A (en) * 2010-04-08 2011-10-12 上海体育学院 Sports training analysis system based on digital site and electromechanical signal
JP5158824B2 (en) * 2009-09-14 2013-03-06 国立大学法人大阪大学 Muscle synergy analysis method, muscle synergy analysis device, and muscle synergy interface

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6002960A (en) * 1996-08-27 1999-12-14 The Johns Hopkins University Passive, non-invasive method to quantify objectively the level and density of a neural blockade
JP5158824B2 (en) * 2009-09-14 2013-03-06 国立大学法人大阪大学 Muscle synergy analysis method, muscle synergy analysis device, and muscle synergy interface
CN201760063U (en) * 2010-04-08 2011-03-16 上海体育学院 Athletic training analysis system based on digital court and electromechanical signal
CN102210925A (en) * 2010-04-08 2011-10-12 上海体育学院 Sports training analysis system based on digital site and electromechanical signal

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Lower-limb muscle synergies in children with cerebral palsy;Fei Li et al.;《6th Annual International IEEE EMBS Conference on Neural Engineering》;20131108;第1226-1229页 *

Also Published As

Publication number Publication date
CN103886215A (en) 2014-06-25

Similar Documents

Publication Publication Date Title
CN103886215B (en) Walking ability analyzing method and device based on muscle collaboration
CN104382595A (en) Upper limb rehabilitation system and method based on myoelectric signal and virtual reality interaction technology
CN104706349B (en) Electrocardiosignal construction method based on pulse wave signals
CN105326500B (en) Action identification method and equipment based on surface electromyogram signal
CN102319067B (en) Nerve feedback training instrument used for brain memory function improvement on basis of electroencephalogram
CN105997064B (en) A kind of discrimination method for human body lower limbs surface electromyogram signal
CN104978035B (en) Brain machine interface system and its implementation based on body-sensing electric stimulus inducing P300
CN104000586A (en) Stroke patient rehabilitation training system and method based on brain myoelectricity and virtual scene
CN110179643A (en) A kind of neck rehabilitation training system and training method based on annulus sensor
CN105708587B (en) A kind of the lower limb exoskeleton training method and system of the triggering of Mental imagery pattern brain-computer interface
CN106569604A (en) Audiovisual dual-mode semantic matching and semantic mismatch co-stimulus brain-computer interface paradigm
CN105249964B (en) Appraisal procedure is rebuild based on magneticencephalogram and the multi-modal brain function of diffusion tensor
CN103610447A (en) Mental workload online detection method based on forehead electroencephalogram signals
CN103735262A (en) Dual-tree complex wavelet and common spatial pattern combined electroencephalogram characteristic extraction method
CN102306303A (en) Electroencephalography signal characteristic extraction method based on small training samples
CN110339024A (en) Lower limb exoskeleton robot and its real-time gait switching method and storage device
CN108268844A (en) Movement recognition method and device based on surface electromyogram signal
CN107440716A (en) Human body lower limbs athletic performance classification discrimination method based on single channel electromyographic signal
CN102178524A (en) Synchronization likehood-based electroencephalograph and electromyography synergistic analyzing method
CN105534517A (en) Method for removing vehicle motion noise in three-lead electrocardiosignal
CN104887198A (en) Pain quantitative analysis system and method based on human body physiological signal multi-parameter fusion
CN104571504A (en) Online brain-machine interface method based on imaginary movement
CN114298089A (en) Multi-mode strength training assisting method and system
Andrade et al. On the relationship between features extracted from EMG and force for constant and dynamic protocols
CN110464348A (en) The continuous amount of exercise recognition methods of joint of lower extremity and system based on electromyography signal

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170111