[go: up one dir, main page]

CN109330604A - A method and device for decomposing a human standing shaking signal - Google Patents

A method and device for decomposing a human standing shaking signal Download PDF

Info

Publication number
CN109330604A
CN109330604A CN201810989143.8A CN201810989143A CN109330604A CN 109330604 A CN109330604 A CN 109330604A CN 201810989143 A CN201810989143 A CN 201810989143A CN 109330604 A CN109330604 A CN 109330604A
Authority
CN
China
Prior art keywords
cop
acc
interp
trajectory
xacc
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
Application number
CN201810989143.8A
Other languages
Chinese (zh)
Other versions
CN109330604B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201810989143.8A priority Critical patent/CN109330604B/en
Publication of CN109330604A publication Critical patent/CN109330604A/en
Application granted granted Critical
Publication of CN109330604B publication Critical patent/CN109330604B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/1036Measuring load distribution, e.g. podologic studies
    • A61B5/1038Measuring plantar pressure during gait
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4005Detecting, measuring or recording for evaluating the nervous system for evaluating the sensory system
    • A61B5/4023Evaluating sense of balance

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physiology (AREA)
  • Dentistry (AREA)
  • Neurology (AREA)
  • Neurosurgery (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开了一种人体站立晃动信号分解方法及装置,本发明基于足下压力分布中心轨迹包含被动游移和主动调整两种运动机制的实际情况,通过相应的插值、海伦公式等数学方法,对轨迹本身进行了分解,得到分别代表两种运动机制的信号,可以有效地分离人体平衡控制过程中的被动游移和主动调整,为全面研究和评价人体平衡能力提供了可靠的工具,更为明确的展现了人体站立时姿态控制的运动过程。诸如视觉、前庭平衡能力的感官调节能力多与人体的主动调整能力相关,其作用更多反映在人体主动调整信号上;而人体的老化程度、肌肉整体稳定性则更多与人体的整体平衡能力相关,更多反映在人体被动游移信号上。

The invention discloses a method and a device for decomposing a human standing shaking signal. Based on the actual situation that the center trajectory of the pressure distribution under the foot includes two motion mechanisms, passive wandering and active adjustment, the trajectory can be analyzed by mathematical methods such as corresponding interpolation and Heron's formula. It decomposes itself, and obtains signals representing two movement mechanisms, which can effectively separate the passive wandering and active adjustment in the process of human balance control, and provide a reliable tool for comprehensive research and evaluation of human balance ability. The motion process of posture control when the human body is standing. The sensory adjustment ability such as vision and vestibular balance ability is mostly related to the active adjustment ability of the human body, and its role is more reflected in the active adjustment signal of the human body; while the aging degree of the human body and the overall stability of the muscles are more related to the overall balance ability of the human body. Correlation, more reflected in the passive wandering signal of the human body.

Description

A kind of human body, which is stood, shakes signal decomposition method and device
Technical field
The invention belongs to technical field of medical detection, it is related to a kind of human body and stands to shake signal decomposition method and device.
Background technique
The weight distribution for the support construction and the upper body lower part of the body that the biped of human body itself is stood determines that human body itself is being stood Unstability immediately, but also its adjustment process is extremely complex while improving balanced ability of human body.Various daily In activity, skeletal musculature, the proprioception of nervous system and human body and visual performance etc. are played in balanced adjustment Different effect.Meanwhile during human senility, with the decline of every function, this human body comprehensive of balanced capacity Function also can be relatively early appearance degeneration sign, therefore the balanced capacity for assessing human body can be used as a routine of old monitoring Index for examination, to realize the classification of physical condition and then implement different nurses.In addition, assessment balanced capacity also contributes to evading Fall risk reduces the injury being likely to occur.
The assessment of balanced capacity generally clinically passes through observation using means, doctors such as clinical observation method, scale Evaluation Methods Performance of the subject under different test items is scored, for example, transform to the performance of standing from sitting posture, stand still when Length, the smooth degree for turning round 360 degree etc..These method subjectivities are strong, and the assessment result of different doctors may have larger difference, difficult To quantify, objectively evaluate, and the relevant information for the main source of subject's equilibrium problem of more having no idea to provide, therefore it is answered Guiding value used in its clinical diagnosis while having compared with big limitation is also very limited.
Centre of body weight can generate involuntary swing when human body stationary vertical is stood, and underfoot then show as pressure distribution And the change at its center, the small sway of underfooting pressure distribution center reflect eyesight, vestibular organ and proprioception system And the influence to balance control such as kinematic system.It can be measured using balance ergograph, record underfooting pressure point when human body is stood The swinging track at cloth center, the measurement of correlation parameter of the track have become the quantitative criterion of assessment balanced ability of human body, can use In research balancing and Correlative Influence Factors.
Meanwhile the correlative study of current stage shows the track of the pressure distribution center of underfooting when human body stationary vertical is stood It is main comprising passively slowly vacillating relatively and the quickly adjustment of more active both locomotory mechanisms, the two are relatively independent Locomotory mechanism respectively with participating in each system of balance system, organ has different degrees of association.From the pressure of underfooting point Extracted respectively in the track at cloth center the relevant information with the two locomotory mechanisms facilitate deeper into research balancing, into And facilitate the assessment of balanced ability of human body.
Summary of the invention
Human body stationary vertical can not independently shake when standing, and underfooting pressure distribution therewith can change.Currently, related side Method is also limited to the direct parameter assessment to foot pressure distribution center track of standing, by calculate the area that track is covered with And movement velocity etc. obtains the evaluation to total balance of the body.But the complexity of balanced ability of human body itself cause so directly from Parameter that track mass motion obtains itself does not characterize total balance of the body regulating power strictly, and it is even more impossible to reflect that participating in balance adjusts The relevant information of each system of section.Importantly, containing human body relatively passively in foot pressure distribution center track It slowly vacillates and initiatively two kinds of locomotory mechanisms of quick adjustment, both locomotory mechanism reciprocal effects, and directly to track itself Parameterizing obtained information is that this two parts moves the description mixed, can not distinguish the same of two kinds of locomotory mechanisms very well When, it is easier to limit the accuracy of the balancing evaluation based on parameter.The present invention is based on such actual conditions, i.e. foot pressure Distribution center track includes passively to vacillate and two kinds of locomotory mechanisms of active accommodation, passes through the mathematics such as corresponding interpolation, Heron's formula Method decomposes track itself, obtains the signal for respectively representing two kinds of locomotory mechanism.The present invention is based on this to propose one kind Human body, which is stood, shakes signal decomposition method and device.
The purpose of the present invention is achieved through the following technical solutions: a kind of human body standing shaking signal decomposition method, Method includes the following steps:
1) foot Center of Pressure track COP (x, y) and time flag sequence T are recorded using test equipment, wherein x is a left side Right direction coordinate, y are front-rear direction coordinate;
2) according to known foot Center of Pressure track COP (x, y), the vector acceleration sequence of track movement is obtained Acc (xacc, yacc) and corresponding time flag sequence T, wherein xacc and yacc is respectively on left and right and front-rear direction Component of acceleration;
3) according to vector acceleration sequence Acc (xacc, yacc) and its corresponding time flag sequence T, two two-phases are calculated Minimum acceleration and its corresponding time label between adjacent sampled point, obtain corresponding sequence Acc_interval and T_interval;
4) acceleration rate threshold (general value 0.2-0.4) is set, screens numerical values recited in sequence Acc_interval and is less than The acceleration of threshold value and its corresponding time tag, obtain corresponding time tag sequence T_interp;
5) position for inscribing original foot Center of Pressure track COP (x, y) when each in sequence label T_interp It sets interpolation reconstruction and obtains representing passively slowly wandering track COP_interp (x_interp, y_interp), and track COP The difference of (x, y) and COP_interp (x_interp, y_interp) is the track COP_ quickly adjusted for representing active Resi (x_resi, y_resi), wherein x_interp, x_resi are the position of its track in the lateral direction, y_interp, y_ Resi is the position of its track in the longitudinal direction.
Further, the step 1) specifically includes following sub-step:
1.1) test equipment includes rectangular platform, four corners of rectangular platform according to after left front, left back, right, it is right before it is suitable Sequence is sequentially arranged four pressure sensors, and the collected signal of institute is followed successively by w1、w2、w3、w4, former and later two pressure sensors The distance of two pressure sensors of distance and left and right is respectively wd and lg;
In sampling instant tnThe coordinate of human body standing underfooting Center of Pressure when on platform:
X_instant=2lg (w1+w2-w3-w4)/(w1+w2+w3+w4)
Y_instant=2wd (w1-w2-w3+w4)/(w1+w2+w3+w4)
1.2) the coordinate COP_instant (x_instant, y_intant) for the instant underfooting Center of Pressure being calculated is passed through After crossing exponential smoothing, sampling instant t is obtainednCoordinate record to track COP (x, y), calculation formula is as follows:
Wherein, α is smoothing factor, takes 0.2-0.4, (xn-1,yn-1) and (xn,yn) it is respectively tn-1And tnThe track at moment Coordinate on COP (x, y).
Further, the step 2) specifically includes following sub-step:
2.1) left and right of the coordinate of track COP (x, y), front-rear direction are carried out at Savitzky-Golay filtering respectively Reason;
2.2) track after filtering processing is calculated separately respectively according to left and right, front-rear direction referring to label T of corresponding time The numerical gradient of a time obtains the velocity vector sequence Vel (xvel, yvel) of track;
2.3) to velocity vector sequence Vel (xvel, yvel) repeat step 2.1) and 2.2) in track COP (x, Y) operation executed, obtains the vector acceleration sequence Acc (xacc, yacc) of track.
Further, the step 3) specifically includes following sub-step:
3.1) t in Acc (xacc, yacc) is calculatedn-1And tnThe adjacent two vector accelerations Acc at momentn-1(xaccn-1, yaccn-1) and Accn(xaccn,yaccn) terminal be wired to the common starting point of two vectors vertical line H length h, calculation formula is such as Under:
λ=(| | Accn-1||+||Accn||+||Accn-1-Accn||)/2
Wherein, λ is auxiliary parameter;
3.2) vertical o'clock of vertical line H is calculated separately to two vector Accn-1(xaccn-1,yaccn-1) and Accn(xaccn,yaccn) Terminal distance, calculation formula is as follows:
Wherein, α1And α2It is the vertical point of vertical line H respectively to Accn-1(xaccn-1,yaccn-1) and Accn(xaccn,yaccn) The distance of terminal, it is positive and negative to represent point position of hanging down, it is all that timing indicates at vertical o'clock on two terminal lines, if there is a negative value, On the reverse extending line of its corresponding two terminals line;
If 3.3) α1·α2> 0, then acceleration representated by the vertical line is the minimum acceleration in the section, and numerical value is accmin=h calculates corresponding time point, and calculation formula is as follows:
If α1·α2< 0, then Accn-1(xaccn-1,yaccn-1) or Accn(xaccn,yaccn) be the section minimum Acceleration determines the numerical value acc of minimum acceleration more afterwardsminAnd corresponding time tmin
3.4) to vector acceleration sequence Acc (xacc, yacc), adjacent vector acceleration repeats above-mentioned calculating two-by-two, obtains Complete acceleration value sequence Acc_interval and corresponding time T_interval.
Further, the step 5) specifically includes following sub-step:
5.1) the upper track COP (x, y) of each moment in time tag sequence T_interp is calculated by cubic spline interpolation Position coordinates, form discrete loci, be denoted as COP_rough;
5.2) connection COP_rough is calculated by cubic spline interpolation, obtains it included in the time tag sequence T When the position coordinates that engrave, form continuous track, i.e. COP_interp (x_interp, y_interp);
5.3) COP_resi (x_ is calculated according to track COP (x, y) and COP_interp (x_interp, y_interp) Resi, y_resi), calculation formula is as follows:
Wherein, COP_interp (x_interp, y_interp) is passively slow wandering track, COP_resi (x_ Resi, y_resi) be active the track quickly adjusted.
The beneficial effects of the present invention are: the present invention devises a kind of human body standing shaking signal point for balanced ability of human body Method and device is solved, passively vacillating in total balance of the body control process and active accommodation can be efficiently separated, for research comprehensively Reliable tool is provided with evaluation balanced ability of human body.By foot Center of Pressure decomposing trajectories at passively vacillating and active accommodation The track of two kinds of locomotory mechanism more specifically presents the motion process of gesture stability when human body is stood.Such as vision, vestibular The sense organ regulating power of balanced capacity is mostly related to the active accommodation ability of human body, and effect is more reflected in human body active accommodation On signal;And the degree of aging of human body, muscle overall stability are then more related to the whole machine balancing ability of human body, it is more to reflect It passively vacillates on signal in human body.The present invention can individually extract the signal of both locomotory mechanisms.
Detailed description of the invention
Fig. 1 is that signal decomposition method flow chart is shaken in human body standing of the present invention;
Fig. 2 is the floor map of rectangular platform and measurement track of the invention, wherein subject is as represented by footprint Towards standing on platform, the circle of dotted line that quadrangle identifies 1,2,3,4 indicates four pressure sensors;
Fig. 3 is the schematic diagram of decomposing trajectories result of the invention, wherein COP, Interp, Resi are foot pressure respectively Center initial trace, passively slowly vacillate track and active quick adjusting track, dotted line is COP in (a), (c), solid line is Interp;
Fig. 4 is the structural block diagram for the balanced capacity detection device that the present invention is suitable for decomposing trajectories technology.
Specific embodiment
The following further describes the present invention in detail with reference to the accompanying drawings and specific embodiments.
A kind of human body proposed by the present invention, which is stood, shakes signal decomposition method, comprising the following steps:
1) foot Center of Pressure track COP (x, y) is recorded using test equipment and time flag sequence T (guarantees duration 10s or more, sample rate are higher than 50Hz), wherein x is left and right directions coordinate, and y is front-rear direction coordinate, and unit is millimeter (mm), when Between the unit second (s);
2) according to known foot Center of Pressure track COP (x, y), the vector acceleration sequence of track movement is obtained Acc (xacc, yacc) and corresponding time flag sequence T, wherein xacc and yacc is respectively on left and right and front-rear direction Component of acceleration;
3) according to vector acceleration sequence Acc (xacc, yacc) and its corresponding time flag sequence T, two two-phases are calculated Minimum acceleration and its corresponding time label between adjacent sampled point, obtain corresponding sequence Acc_interval and T_interval;
4) acceleration rate threshold (general value 10-50mm/s appropriate is set2, value 25mm/s herein2), screen sequence Numerical values recited is less than the acceleration and its corresponding time tag of threshold value in Acc_interval, obtains corresponding time tag sequence T_ interp;
5) position for inscribing original foot Center of Pressure track COP (x, y) when each in sequence label T_interp It sets interpolation reconstruction and obtains representing passively slowly wandering track COP_interp (x_interp, y_interp), and track COP The difference of (x, y) and COP_interp (x_interp, y_interp) is the track COP_ quickly adjusted for representing active Resi (x_resi, y_resi), wherein x_interp, x_resi are the position of its track in the lateral direction, y_interp, y_ Resi is the position of its track in the longitudinal direction.
Further, the step 1) specifically includes following sub-step:
1.1) test equipment includes rectangular platform, four corners of rectangular platform according to after left front, left back, right, it is right before it is suitable Sequence is sequentially arranged four pressure sensors, and the collected signal of institute is followed successively by w1、w2、w3、w4, former and later two pressure sensors The distance of two pressure sensors of distance and left and right is respectively wd and lg;In sampling instant tnHuman body standing foot when on platform The coordinate of lower Center of Pressure:
X_instant=2lg (w1+w2-w3-w4)/(w1+w2+w3+w4)
Y_instant=2wd (w1-w2-w3+w4)/(w1+w2+w3+w4)
1.2) the coordinate COP_instant (x_instant, y_intant) for the instant underfooting Center of Pressure being calculated is passed through After crossing exponential smoothing, sampling instant t is obtainednCoordinate record to track COP (x, y), calculation formula is as follows:
Wherein, α is smoothing factor (general value 0.2-0.4, herein value 0.25), (xn-1,yn-1) and (xn,yn) respectively For tn-1And tnCoordinate on the track COP (x, y) at moment.
Further, the step 2) specifically includes following sub-step:
2.1) left and right of the coordinate of track COP (x, y), front-rear direction are subjected to Savitzky-Golay filtering processing respectively (recommending window width is 7, recommendation order is 3);
2.2) track after filtering processing is calculated separately respectively according to left and right, front-rear direction referring to label T of corresponding time The numerical gradient of a time obtains the velocity vector sequence Vel (xvel, yvel) of track;
2.3) to velocity vector sequence Vel (xvel, yvel) repeat step 2.1) and 2.2) in track COP (x, Y) operation executed, obtains the vector acceleration sequence Acc (xacc, yacc) of track.
Further, the step 3) specifically includes following sub-step:
3.1) t in Acc (xacc, yacc) is calculatedn-1And tnThe adjacent two vector accelerations Acc at momentn-1(xaccn-1, yaccn-1) and Accn(xaccn,yaccn) terminal be wired to the common starting point of two vectors vertical line H length h, calculation formula is such as Under:
λ=(| | Accn-1||+||Accn||+||Accn-1-Accn||)/2
Wherein, λ is auxiliary parameter;
3.2) vertical o'clock of vertical line H is calculated separately to two vector Accn-1(xaccn-1,yaccn-1) and Accn(xaccn,yaccn) Terminal distance, calculation formula is as follows:
Wherein, α1And α2It is the vertical point of vertical line H respectively to Accn-1(xaccn-1,yaccn-1) and Accn(xaccn,yaccn) The distance of terminal, it is positive and negative to represent point position of hanging down, it is all that timing indicates at vertical o'clock on two terminal lines, if there is a negative value, On the reverse extending line of its corresponding two terminals line;
If 3.3) α1·α2> 0, then acceleration representated by the vertical line is the minimum acceleration in the section, and numerical value is accmin=h calculates corresponding time point, and calculation formula is as follows:
If α1·α2< 0, then Accn-1(xaccn-1,yaccn-1) or Accn(xaccn,yaccn) be the section minimum Acceleration determines the numerical value acc of minimum acceleration more afterwardsminAnd corresponding time tmin
3.4) to vector acceleration sequence Acc (xacc, yacc), adjacent vector acceleration repeats above-mentioned calculating two-by-two, obtains Complete acceleration value sequence Acc_interval and corresponding time T_interval.
Further, the step 5) specifically includes following sub-step:
5.1) the upper track COP (x, y) of each moment in time tag sequence T_interp is calculated by cubic spline interpolation Position coordinates, form discrete loci, be denoted as COP_rough;
5.2) connection COP_rough is calculated by cubic spline interpolation, obtains it included in the time tag sequence T When the position coordinates that engrave, form continuous track, i.e. COP_interp (x_interp, y_interp);
5.3) COP_resi (x_ is calculated according to track COP (x, y) and COP_interp (x_interp, y_interp) Resi, y_resi), calculation formula is as follows:
Wherein, COP_interp (x_interp, y_interp) is passively slow wandering track, COP_resi (x_ Resi, y_resi) be active the track quickly adjusted.
A kind of balanced capacity detection device suitable for decomposing trajectories technology proposed by the present invention, the detection device include: Medical D.C. regulated power supply, DC-DC isolation module, Voltage stabilizing module, rectangular platform, pressure sensor, signal imitation processing module, Multiple signals synchronize A/D module, microprocessor module and magnetic coupling isolation module;
The medical D.C. regulated power supply, DC-DC isolation module, Voltage stabilizing module are sequentially connected;
Four corners of the rectangular platform arrange that pressure sensor, four pressure sensors connect at signal imitation respectively After reason mould is filtered and amplifies, it is commonly connected to the synchronous A/D module of multiple signals;
The synchronous A/D module of the multiple signals, microprocessor module, magnetic coupling isolation module are sequentially connected;
The DC-DC isolation module is connect with four pressure sensors;The Voltage stabilizing module be signal imitation processing module, Multiple signals synchronize A/D module, microprocessor module and magnetic coupling isolation and provide burning voltage;
The synchronous A/D module of the multiple signals by input clamping protection, frequency overlapped-resistable filter, sample/hold amplifier with And high-speed ADC realizes the sampling of signal, guarantees the synchronization of four road signals;
The record for the data that the microprocessor module is used to complete to sample AD, storage, upload;
The magnetic coupling isolation module guarantees medical number for realizing the electrical isolation of data acquisition circuit and upstream port According to the electrical safety of acquisition equipment.
Signal processing method in the present invention can help to obtain more relevant slow with central nervous system control ability Ingredient and with the more relevant fast component track of peripheral neverous system control ability, by the discrete journey for calculating this two parts track The indexs of correlation such as degree and movement velocity can assess stability and its adjusting that the body of tested crowd itself is stood respectively The power of function, to provide more perfect balanced capacity evaluation.
EO EC
COP Interp Resi COP Interp Resi
RMSD ML 3.11 2.89 0.75 3.27 2.98 1.11a
(mm) (0.96) (0.92) (0.45) (1.08) (1.03) (0.73)
mVel ML 5.20 3.41 2.75 6.19a 3.53 3.94a
(mm/s) (1.63) (0.82) (1.55) (2.34) (0.84) (2.39)
RMSD AP 4.24 4.03 0.94 4.50 4.12 1.47a
(mm) (1.26) (1.27) (0.54) (1.41) (1.36) (0.93)
mVel AP 6.59 4.20 3.71 8.52a 4.65a 5.64a
(mm/s) (1.83) (0.97) (1.83) (2.81) (0.98) (3.04)
Upper table is the result of 124 normal personnel of actual measurement, in which: ML, AP respectively represent left and right and front-rear direction;EO, EC, which is represented, to open eyes and closes one's eyes;RMSD and mVel is respectively trajectory-based calculating parameter, RMSD full name root-mean-square Deviation, that is, root-mean-square-deviation, characterize the dispersion degree of data, and mVel is the average speed of track.In table, COP, Interp, Resi respectively indicate the original track COP, slow component track and fast component track in decomposition result, and subscript a is indicated The significance of difference (P < 0.05) between each each parameter eye opening of signal is closed one's eyes.We by closing one's eyes it can be found that separated by probation After the visual information for guaranteeing balance, increase, other than the increased variation of speed in addition to COP itself shakes range, slow component variation It is unobvious, and fast component variation is obvious, shaking amplitude and speed significantly increase, that is to say, that shaking degree after human body is closed one's eyes Raising the overwhelming majority be derived from fast component variations.This can illustrate to close one's eyes have no effect under human body midstance it is intrinsic Stability, but because it has separated visual information, human body transfers to keep standing balance by frequent attitude regulation.
Technical solution of the present invention and beneficial effect is described in detail in above-described specific embodiment, Ying Li Solution is not intended to restrict the invention the foregoing is merely presently most preferred embodiment of the invention, all in principle model of the invention Interior done any modification, supplementary, and equivalent replacement etc. are enclosed, should all be included in the protection scope of the present invention.

Claims (7)

1.一种人体站立晃动信号分解方法,其特征在于,该方法包括以下步骤:1. a human body standing shaking signal decomposition method, is characterized in that, the method comprises the following steps: 1)使用测试设备记录足下压力中心轨迹COP(x,y)以及时间标记序列T,其中x为左右方向坐标,y为前后方向坐标;1) Use the test equipment to record the foot pressure center track COP(x, y) and the time mark sequence T, where x is the left and right direction coordinates, and y is the front and rear direction coordinates; 2)根据已知的足下压力中心轨迹COP(x,y),得到该轨迹运动的加速度向量序列Acc(xacc,yacc)以及对应的时间标记序列T,其中xacc和yacc分别为左右和前后方向上的加速度分量;2) According to the known foot pressure center trajectory COP(x, y), the acceleration vector sequence Acc(xacc, yacc) of the trajectory motion and the corresponding time stamp sequence T are obtained, where xacc and yacc are the left and right and front and rear directions respectively. The acceleration component of ; 3)根据加速度向量序列Acc(xacc,yacc)及其对应的时间标记序列T,计算两两相邻采样点间的最小加速度及其对应时间标记,得到对应序列Acc_interval和T_interval;3) According to the acceleration vector sequence Acc(xacc, yacc) and its corresponding time stamp sequence T, calculate the minimum acceleration between two adjacent sampling points and its corresponding time stamp, and obtain the corresponding sequences Acc_interval and T_interval; 4)设定加速度阈值,筛选序列Acc_interval中数值大小小于阈值的加速度及其对应时间标签,得到对应时间标签序列T_interp;4) Set the acceleration threshold, filter the accelerations whose numerical value is less than the threshold in the sequence Acc_interval and their corresponding time labels, and obtain the corresponding time label sequence T_interp; 5)将原始的足下压力中心轨迹COP(x,y)在标签序列T_interp中各个时刻下的位置插值重构得到代表被动的缓慢游移的轨迹COP_interp(x_interp,y_interp),而轨迹COP(x,y)与COP_interp(x_interp,y_interp)的差值即为代表主动的快速调节的轨迹COP_resi(x_resi,y_resi),其中x_interp、x_resi为其轨迹在左右方向上的位置,y_interp、y_resi为其轨迹在前后方向上的位置。5) Interpolate and reconstruct the position of the original underfoot pressure center trajectory COP(x, y) at each moment in the label sequence T_interp to obtain the trajectory COP_interp(x_interp, y_interp) representing the passive slow movement, and the trajectory COP(x, y) ) and the difference between COP_interp(x_interp, y_interp) is the trajectory COP_resi(x_resi, y_resi) representing active rapid adjustment, where x_interp and x_resi are the positions of the trajectory in the left and right directions, and y_interp and y_resi are the trajectory in the front and rear directions. position on. 2.根据权利要求1所述的一种人体站立晃动信号分解方法,其特征在于,所述步骤1)具体包括以下子步骤:2. a kind of human body standing shaking signal decomposition method according to claim 1, is characterized in that, described step 1) specifically comprises following sub-step: 1.1)测试设备包括矩形平台,矩形平台的四个边角按照左前、左后、右后、右前顺序依次布置四个压力传感器,所采集到的信号依次为w1、w2、w3、w4,前后两个压力传感器的距离以及左右两个压力传感器的距离分别为wd和lg;在采样时刻tn人体站立在平台上时脚下压力中心的坐标:1.1) The test equipment includes a rectangular platform. Four pressure sensors are arranged on the four corners of the rectangular platform in the order of front left, rear left, rear right, and front right. The collected signals are w 1 , w 2 , w 3 , w 4. The distance between the front and rear pressure sensors and the distance between the left and right pressure sensors are wd and lg respectively; the coordinates of the pressure center under the foot when the human body stands on the platform at the sampling time t n : x_instant=2·lg·(w1+w2-w3-w4)/(w1+w2+w3+w4)x_instant=2·lg·(w 1 +w 2 -w 3 -w 4 )/(w 1 +w 2 +w 3 +w 4 ) y_instant=2·wd·(w1-w2-w3+w4)/(w1+w2+w3+w4)y_instant=2·wd·(w 1 -w 2 -w 3 +w 4 )/(w 1 +w 2 +w 3 +w 4 ) 1.2)计算得到的即时脚下压力中心的坐标COP_instant(x_instant,y_intant)经过指数平滑后,得到采样时刻tn的坐标记录到轨迹COP(x,y),计算公式如下:1.2) After exponential smoothing, the coordinates COP_instant(x_instant, y_intant) of the immediate pressure center under the foot are calculated to obtain the coordinates of the sampling time t n and record them in the trajectory COP(x, y). The calculation formula is as follows: 其中,α为平滑系数,取0.2-0.4,(xn-1,yn-1)和(xn,yn)分别为tn-1和tn时刻的轨迹COP(x,y)上的坐标。Among them, α is the smoothing coefficient, taking 0.2-0.4, (x n-1 , y n-1 ) and (x n , y n ) are the trajectory COP(x, y) at the time t n-1 and t n respectively coordinate of. 3.根据权利要求1所述的一种人体站立晃动信号分解方法,其特征在于,所述步骤2)具体包括以下子步骤:3. a kind of human body standing shaking signal decomposition method according to claim 1, is characterized in that, described step 2) specifically comprises following sub-step: 2.1)将轨迹COP(x,y)的坐标的左右、前后方向分别进行Savitzky-Golay滤波处理;2.1) Perform Savitzky-Golay filtering processing on the left and right and front and rear directions of the coordinates of the trajectory COP(x, y) respectively; 2.2)将滤波处理后的轨迹按照左右、前后方向参照对应的时间标记T分别计算各个时间的数值梯度,得到轨迹的速度向量序列Vel(xvel,yvel);2.2) Calculate the numerical gradient of each time with reference to the corresponding time mark T according to the left and right and front and rear directions of the filtered trajectory to obtain the velocity vector sequence Vel(xvel, yvel) of the trajectory; 2.3)对速度向量序列Vel(xvel,yvel)重复执行步骤2.1)和2.2)中对轨迹COP(x,y)执行的操作,得到轨迹的加速度向量序列Acc(xacc,yacc)。2.3) Repeat the operations performed on the trajectory COP(x, y) in steps 2.1) and 2.2) for the velocity vector sequence Vel(xvel, yvel) to obtain the trajectory acceleration vector sequence Acc(xacc, yacc). 4.根据权利要求1所述的一种人体站立晃动信号分解方法,其特征在于,所述步骤3)具体包括以下子步骤:4. a kind of human body standing shaking signal decomposition method according to claim 1, is characterized in that, described step 3) specifically comprises following sub-step: 3.1)计算Acc(xacc,yacc)中tn-1和tn时刻的相邻两加速度向量Accn-1(xaccn-1,yaccn-1)和Accn(xaccn,yaccn)的终点连线到两向量共同起点的垂线H的长度h,计算公式如下:3.1) Calculate the two adjacent acceleration vectors Acc n-1 (xacc n-1 , yacc n-1 ) and Acc n (xacc n , yacc n ) at the time t n-1 and t n in Acc (xacc, yacc) The length h of the vertical line H connecting the end point to the common starting point of the two vectors is calculated as follows: λ=(||Accn-1||+||Accn||+||Accn-1-Accn||)/2λ=(||Acc n-1 ||+||Acc n ||+||Acc n-1 -Acc n ||)/2 其中,λ为辅助参量;Among them, λ is an auxiliary parameter; 3.2)分别计算垂线H的垂点到两向量Accn-1(xaccn-1,yaccn-1)和Accn(xaccn,yaccn)的终点的距离,计算公式如下:3.2) Calculate the distance from the vertical point of the vertical line H to the end points of the two vectors Acc n-1 (xacc n-1 , yacc n-1 ) and Acc n (xacc n , yacc n ) respectively, and the calculation formula is as follows: 其中,α1和α2分别是垂线H的垂点到Accn-1(xaccn-1,yaccn-1)和Accn(xaccn,yaccn)的终点的距离,正负代表垂点位置,同为正时表示垂点在两终点连线上,若有一个负值,则在其对应的两终点连线的反向延长线上;Among them, α 1 and α 2 are the distances from the vertical point of the vertical line H to the end points of Acc n-1 (xacc n-1 , yacc n-1 ) and Acc n (xacc n , yacc n ), respectively, positive and negative represent the vertical Point position, when both are positive, it means that the vertical point is on the line connecting the two end points, if there is a negative value, it is on the reverse extension line of the corresponding line connecting the two end points; 3.3)如果α1·α2&gt;0,则该垂线所代表的加速度为该区间的最小加速度,数值为accmin=h,计算对应的时间点,计算公式如下:3.3) If α 1 ·α 2 &gt; 0, the acceleration represented by the vertical line is the minimum acceleration in the interval, the value is acc min =h, and the corresponding time point is calculated, and the calculation formula is as follows: 如果α1·α2&lt;0,则Accn-1(xaccn-1,yaccn-1)或者Accn(xaccn,yaccn)为该区间的最小加速度,比较后确定最小加速度的数值accmin以及对应的时间tminIf α 1 ·α 2 &lt;0, then Acc n-1 (xacc n-1 , yacc n-1 ) or Acc n (xacc n , yacc n ) is the minimum acceleration in the interval, and the value of the minimum acceleration is determined after comparison acc min and the corresponding time t min ; 3.4)对加速度向量序列Acc(xacc,yacc)两两相邻加速度向量重复上述计算,得到完整的加速度数值序列Acc_interval和对应的时间T_interval。3.4) Repeat the above calculation for the adjacent acceleration vectors of the acceleration vector sequence Acc(xacc, yacc) to obtain the complete acceleration numerical sequence Acc_interval and the corresponding time T_interval. 5.根据权利要求1所述的一种人体站立晃动信号分解方法,其特征在于,所述步骤4)中加速度阈值取10-50mm/s25 . The method for decomposing a standing shaking signal of a human body according to claim 1 , wherein the acceleration threshold in the step 4) is 10-50 mm/s 2 . 6 . 6.根据权利要求1所述的一种人体站立晃动信号分解方法,其特征在于,所述步骤5)具体包括以下子步骤:6. a kind of human body standing shaking signal decomposition method according to claim 1, is characterized in that, described step 5) specifically comprises following sub-step: 5.1)通过三次样条插值计算时间标签序列T_interp中各个时刻上轨迹COP(x,y)的位置坐标,组成离散轨迹,记为COP_rough;5.1) Calculate the position coordinates of the trajectory COP(x, y) at each moment in the time label sequence T_interp through cubic spline interpolation to form a discrete trajectory, denoted as COP_rough; 5.2)通过三次样条插值计算连接COP_rough,得到其在时间标签序列T中所包含的时刻上的位置坐标,形成连续的轨迹,即COP_interp(x_interp,y_interp);5.2) Calculate and connect COP_rough through cubic spline interpolation to obtain its position coordinates at the moment contained in the time label sequence T, forming a continuous trajectory, namely COP_interp(x_interp, y_interp); 5.3)根据轨迹COP(x,y)和COP_interp(x_interp,y_interp)计算COP_resi(x_resi,y_resi),计算公式如下:5.3) Calculate COP_resi(x_resi,y_resi) according to the trajectory COP(x,y) and COP_interp(x_interp,y_interp), and the calculation formula is as follows: 其中,COP_interp(x_interp,y_interp)为被动的缓慢游移的轨迹,COP_resi(x_resi,y_resi)为主动的快速调节的轨迹。Among them, COP_interp(x_interp, y_interp) is a passive slow-moving trajectory, and COP_resi(x_resi, y_resi) is an active fast-adjusting trajectory. 7.一种适用于轨迹分解技术的平衡能力检测装置,其特征在于,该检测装置包括:医用直流稳压电源、DC-DC隔离模块、稳压模块、矩形平台、压力传感器、信号模拟处理模块、多路信号同步AD模块、微处理器模块和磁耦合隔离模块;7. A balance ability detection device suitable for trajectory decomposition technology, characterized in that the detection device comprises: a medical DC stabilized power supply, a DC-DC isolation module, a voltage stabilizer module, a rectangular platform, a pressure sensor, and a signal analog processing module , Multi-channel signal synchronization AD module, microprocessor module and magnetic coupling isolation module; 所述医用直流稳压电源、DC-DC隔离模块、稳压模块依次连接;The medical DC stabilized power supply, the DC-DC isolation module, and the voltage stabilizer module are connected in sequence; 所述矩形平台的四个边角均布置压力传感器,四个压力传感器分别接信号模拟处理模进行滤波和放大后,共同连接到多路信号同步AD模块;The four corners of the rectangular platform are arranged with pressure sensors, and the four pressure sensors are respectively connected to the signal analog processing module for filtering and amplification, and then jointly connected to the multi-channel signal synchronization AD module; 所述多路信号同步AD模块、微处理器模块、磁耦合隔离模块依次连接;The multi-channel signal synchronization AD module, the microprocessor module and the magnetic coupling isolation module are connected in sequence; 所述DC-DC隔离模块与四个压力传感器连接;所述稳压模块为信号模拟处理模块、多路信号同步AD模块、微处理器模块和磁耦合隔离提供稳定电压;The DC-DC isolation module is connected with four pressure sensors; the voltage stabilization module provides stable voltage for the signal analog processing module, the multi-channel signal synchronous AD module, the microprocessor module and the magnetic coupling isolation; 所述多路信号同步AD模块通过输入箝位保护、抗混叠滤波器、采样保持放大器以及高速ADC实现信号的采样,保证四路信号的同步;The multi-channel signal synchronization AD module realizes signal sampling through input clamping protection, anti-aliasing filter, sample and hold amplifier and high-speed ADC, so as to ensure the synchronization of four channels of signals; 所述微处理器模块用于完成对AD采样得到的数据的记录、存储、上传;The microprocessor module is used to complete the recording, storage and uploading of the data obtained by AD sampling; 所述磁耦合隔离模块用于实现数据采集电路与上游端口的电气隔离,保证医用数据采集设备的电气安全。The magnetic coupling isolation module is used to realize the electrical isolation of the data acquisition circuit and the upstream port, so as to ensure the electrical safety of the medical data acquisition equipment.
CN201810989143.8A 2018-08-28 2018-08-28 Human body standing shaking signal decomposition method and device Active CN109330604B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810989143.8A CN109330604B (en) 2018-08-28 2018-08-28 Human body standing shaking signal decomposition method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810989143.8A CN109330604B (en) 2018-08-28 2018-08-28 Human body standing shaking signal decomposition method and device

Publications (2)

Publication Number Publication Date
CN109330604A true CN109330604A (en) 2019-02-15
CN109330604B CN109330604B (en) 2020-06-09

Family

ID=65291939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810989143.8A Active CN109330604B (en) 2018-08-28 2018-08-28 Human body standing shaking signal decomposition method and device

Country Status (1)

Country Link
CN (1) CN109330604B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111528804A (en) * 2020-05-21 2020-08-14 上海体育学院 A detection device system for human balance ability

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5388591A (en) * 1992-12-09 1995-02-14 Trustees Of Boston University Method and apparatus for analyzing the human postural control system
CN1771886A (en) * 2005-08-31 2006-05-17 南京大学 Human body balance function detection method and training system
US20060293613A1 (en) * 2005-06-27 2006-12-28 Concept Development Group Method and Apparatus for Automated Monitoring and Tracking of the Trajectory of Patients' Center of Gravity Movements
EP2422698A1 (en) * 2010-08-26 2012-02-29 Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO Method and system for determining the walking or running speed of a person
CN103417216A (en) * 2013-04-08 2013-12-04 杭州电子科技大学 Pressure center nonlinear feature extraction method based on complexity
US20140330171A1 (en) * 2013-05-03 2014-11-06 National Taiwan University Method and device for monitoring postural and movement balance for fall prevention
CN104636599A (en) * 2014-12-31 2015-05-20 北京大学 Plantar pressure center track automatic comparison method
US9072463B2 (en) * 2009-01-27 2015-07-07 University Of Washington Prosthetic limb monitoring system
KR20150128067A (en) * 2014-05-08 2015-11-18 창명제어기술 (주) Method for correcting data of Center Of Pressure using image processing technique
CN106618498A (en) * 2016-12-14 2017-05-10 广州市章和电气设备有限公司 System and method for evaluating body balance
CN107169406A (en) * 2017-03-28 2017-09-15 中山大学 It is a kind of to extract the method that balance vibrates initial time by power platform Center of Pressure data
CN108185994A (en) * 2018-02-02 2018-06-22 杭州拜伦医疗科技有限公司 A kind of balanced capacity detection method and device based on nonlinear technology
CN108309236A (en) * 2018-01-15 2018-07-24 新绎健康科技有限公司 Total balance of the body appraisal procedure and system

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5388591A (en) * 1992-12-09 1995-02-14 Trustees Of Boston University Method and apparatus for analyzing the human postural control system
US20060293613A1 (en) * 2005-06-27 2006-12-28 Concept Development Group Method and Apparatus for Automated Monitoring and Tracking of the Trajectory of Patients' Center of Gravity Movements
CN1771886A (en) * 2005-08-31 2006-05-17 南京大学 Human body balance function detection method and training system
US9072463B2 (en) * 2009-01-27 2015-07-07 University Of Washington Prosthetic limb monitoring system
EP2422698A1 (en) * 2010-08-26 2012-02-29 Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO Method and system for determining the walking or running speed of a person
CN103417216A (en) * 2013-04-08 2013-12-04 杭州电子科技大学 Pressure center nonlinear feature extraction method based on complexity
US20140330171A1 (en) * 2013-05-03 2014-11-06 National Taiwan University Method and device for monitoring postural and movement balance for fall prevention
KR20150128067A (en) * 2014-05-08 2015-11-18 창명제어기술 (주) Method for correcting data of Center Of Pressure using image processing technique
CN104636599A (en) * 2014-12-31 2015-05-20 北京大学 Plantar pressure center track automatic comparison method
CN106618498A (en) * 2016-12-14 2017-05-10 广州市章和电气设备有限公司 System and method for evaluating body balance
CN107169406A (en) * 2017-03-28 2017-09-15 中山大学 It is a kind of to extract the method that balance vibrates initial time by power platform Center of Pressure data
CN108309236A (en) * 2018-01-15 2018-07-24 新绎健康科技有限公司 Total balance of the body appraisal procedure and system
CN108185994A (en) * 2018-02-02 2018-06-22 杭州拜伦医疗科技有限公司 A kind of balanced capacity detection method and device based on nonlinear technology

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MANCINI, MARTINA等: "Trunk accelerometry reveals postural instability in untreated Parkinson’s disease", 《PARKINSONISM & RELATED DISORDERS》 *
MOE-NILSSEN, R; HELBOSTAD, JL: "Trunk accelerometry as a measure of balance control during quiet standing", 《GAIT & POSTURE 》 *
SLOMKA, KAJETAN等: "The reliability of a rambling-trembling analysis of center of pressure measures", 《GAIT & POSTURE 》 *
ZATSIORSKY, VM, DUARTE, M: "Instant Equilibrium point and its migration in standing tasks rambling and trembling components of stabilogram", 《MOTOR CONTROL 》 *
贾桂锋: "人体平衡能力定量评估方法研究", 《中国博士学位论文全文数据库 医药卫生科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111528804A (en) * 2020-05-21 2020-08-14 上海体育学院 A detection device system for human balance ability
CN111528804B (en) * 2020-05-21 2023-02-17 上海体育学院 Detection device system for human body balance ability

Also Published As

Publication number Publication date
CN109330604B (en) 2020-06-09

Similar Documents

Publication Publication Date Title
Boto et al. Wearable neuroimaging: Combining and contrasting magnetoencephalography and electroencephalography
Noreika et al. 14 challenges and their solutions for conducting social neuroscience and longitudinal EEG research with infants
KR101247748B1 (en) Apparatus for analysing focus and nonfocus states and method therof
Murayama et al. Relationship between neural and hemodynamic signals during spontaneous activity studied with temporal kernel CCA
KR101535432B1 (en) Contents valuation system and contents valuating method using the system
MX2013007965A (en) Posturographic system using a balance board.
CN101720200A (en) Be used for participating in the method and system that reaction is measured and graded
JP2010520554A5 (en)
Sforza et al. Facial movement before and after masseteric-facial nerves anastomosis: a three-dimensional optoelectronic pilot study
Cudejko et al. Validity and reliability of accelerations and orientations measured using wearable sensors during functional activities
CN211355518U (en) A depression diagnosis device
Foy et al. Dynamic changes in EEG power spectral densities during NIH-toolbox flanker, dimensional change card sort test and episodic memory tests in young adults
Kaur et al. Using virtual reality to examine the neural and physiological anxiety-related responses to balance-demanding target-reaching leaning tasks
CN109330604A (en) A method and device for decomposing a human standing shaking signal
Harrison et al. EEG and fMRI agree: Mental arithmetic is the easiest form of imagery to detect
CN107997760B (en) Relevance evaluation method for human body physiological and psychological measurement data
Guo et al. Weakening of recognition memory caused by task-irrelevant emotional encoding context can be modulated by individuals' inhibitory control
CN107595295B (en) A kind of lower extremity movement resistive state recognition methods based on brain hemoglobin information
Mutha et al. Medical Image Fusion Solid Works Supporting Multiple Techniques with Feature-Level Transforms
Jervis et al. Residual ocular artefact subsequent to ocular artefact removal from the electroencephalogram
Solis et al. Frontoparietal network and neuropsychological measures in typically developing children
Hussain et al. Novel data fusion approach for drowsiness detection
Fisher Head-body differentiations in body image and skin resistance level.
Bhoomika Patel et al. Eeg data analysis for iq test using machine learning approaches: a survey
Akano et al. Investigating social cognition in infants and adults using dense array electroencephalography (dEEG)

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