CN103886215B - Walking ability analyzing method and device based on muscle collaboration - Google Patents
Walking ability analyzing method and device based on muscle collaboration Download PDFInfo
- 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
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
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:
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:
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:
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:
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
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:
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:
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:
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:
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:
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.
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)
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)
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 |
-
2014
- 2014-04-04 CN CN201410136557.8A patent/CN103886215B/en not_active Expired - Fee Related
Patent Citations (4)
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)
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 |