arXiv:cond-mat/0212203v1 [cond-mat.stat-mech] 9 Dec 2002 stepro ftegi yl hti liaeymea- the ultimately not is and that experiments, it interval cycle that stride gait stress these We the in of cycle. sured period gait con- the the they of is together locomotion and the system, trol motocontrol cou- is the system to nervous central pled the that assumes walk. (SCPG), we how determine is to it CPG single limb, a each than have for to rather oscillator to possible Thus, nonlinear enslaved separate activity. become a muscular having to rhythmic assembly single an a within property allow neurons systems, the for dynamical using nonlinear models, of gallop- synchronization dynamical loco- of and recent of cantering More trotting, modes ing. walking, various legs is, between that four distinguish motion, the to among de- order quadraped relations to phase in a oscillators correct nonlinear Therefore the four termine of [9]. coupling participating the process limb single required each locomotion a for the that used assumed in be gait oscillator for dy- CPGs nonlinear neurons nonlinear of early of models The namical network of output. capable syncopated a [8], a (CPG) through producing Generator Pattern regulated Central also a called is it but multi- weakly are 7]. [6, monofractal, fractal being than time fractal. rather these is that series, determined series was time it recently gait supported more human However, 4] the [3, that Griffin conclusion activity. and the West fractal by self-similar, studies a Subsequent is the walking time that of stride-interval suggesting phenomenon the correlation, that long-time Haus- demonstrate did exhibits [1] decade series last al. the et humans in of only dorff interval that and stride century [5], the a walking in over during 3-4% for of study known variation is to been there 4] has of 3, It 2, application data. [1, interval the techniques processing through data understand nonlinear to begun cently h oe htw rsn ee h ue CPG super the here, present we that model The process, voluntary a be to known is locomotion Human re- only have we which process complex a is Walking ASnmes 54.a8.8S,0.5D,87.17.Nn 05.45.Df, way. 05.45.-a,87.18.Sn, particular numbers: a PACS in walking locom of act the oscillat conscious control nonlinear ce the forced they cont a The with by is together simulated is output. locomotion and gait metronomic human syncopated system, The correlated motocontrol a metronome. the producing a of on m keyed capable long-range is the Moreover, gait gait. m the average become the fluctuations in the decrease of and nature fractal The fluctuations. eie.Tesrd nevltm eisi omlhmngai human normal in series time interval stride The regimes. epeetanniersohsi oe ftehmngi con gait human the of model stochastic nonlinear a present We 1 .INTRODUCTION I. rt colo EDprmn,Dk nvriy ..Bx902 Box P.O. University, Duke Department, EE of School Pratt 3 ahmtc iiin ryRsac ffie eerhTriang Research Office, Research Army Division, Mathematics 2 hsc eatet ueUiest,Dra,N 70.and 27701. NC Durham, University, Duke Department, Physics olna yaia oe fhmngait human of model dynamical nonlinear A rc .West J. Bruce Dtd etme 7 2018) 27, September (Dated: 1 , 2 , 3 n ioaScafetta. Nicola and lob ncntan htfre h siltrt follow to oscillator the frequency. fixed forces controlled particular that be a constraint may an that by behavior also complex oscillator, behav- Pol more chaotic der Van a the the allows fact, like oscillators, In nonlinear activity. of induced ior firing frequency neural potential frequency the gait the by the from way of slightly this differ period In may actual oscillator. stride the Pol der The with Van coincide cycle. the to gait assumed the der is mimics Van interval that forced [10] nonlinear un- oscillator a the Pol of only frequency regulate inner centers perturbed of neural gait impulses firing the complex correlated of the amplitude the the model that We the assuming by following system metronome. by a example, may of for gate frequency forced, human voluntarily of be dynamics also The activity. firing neural u,i eea,i ahcs hr sardcini the in reduction a is there case each in sig- general, change fluctuations, in series antipersistent but, to time interval persistent stochastic from stride nificantly, the the metronome, a of Moreover, by properties constrained paces. is normal pace the the of slow if that both above for pace, the increase fast to of and multifractality tends series, the partic- time as In interval well stride [7]. as regimes persistence, mode the gait ular, different the variability in fractal se- its time changes a and The non-stationary with 1/f-noise. sometimes the is multifractal of ries the that weakly that to is close shown fractality series main was time it interval [13], stride transform wavelet using tra set. data a second by a generate determined individuals paces to pace same fast minutes and a The 30 normal slow, for at average the walk hour. at fast set to one metronome and requested of normal then period slow, are a 10 a for for at series walking pace, time men interval walking. young stride of healthy rate the average contain the of data of dimension changes These fractal on the series time of [12] the dependence al. data the et These Hausdorff determine by to used [11]. and collected Physionet originally archives were domain public the yetmtn h ¨le xoet n hi spec- their and H¨older exponents the estimating By in available set data a on tested is model SCPG The rwt eidcetra oc associated force external periodic a with or scaatrzdb lgtymultifractal slightly by characterized is t mr nteeflcutosi otwhen lost is fluctuations these in emory r rnucdudrbt nincrease an both under pronounced ore ta evu ytmi ope to coupled is system nervous ntral to ftegi yl tef The itself. cycle gait the of otion rlsse navreyo gait of variety a in system trol oldb ewr fneurons of network a by rolled 1 uhm C27701. NC Durham, 91, 1 , ePr,N 27709. NC Park, le 2 2 long-term memory and an increase in randomness. tiple gait patterns, is shown to not depend on the detailed In Sec. 2 we give a short introduction to the phe- dynamics of the constituent nonlinear oscillators, nor on nomenon of locomotion, the traditional methods for mod- their inter-oscillator coupling strengths [8]. eling using the CPG, and review the data processing used As we mentioned in the Introduction, it has been to establish the fractal behavior of stride interval time known for over a century that there is a variation in series. Sec. 3 reviews the stochastic properties of the the stride interval of humans during walking of approxi- normal and metronomic gait under different various pace mately 3-4% [5]. This random variability has been shown velocities, slow, normal and fast. In Sec. 4 we present [1, 3, 4, 6, 12] to exhibit long-time correlations, and sug- the mathematical details of the SCPG model. In Sec. 5 gested that the phenomenon of walking is a self-similar, we compare the results of computation using the SCPG fractal, activity. The existence of fractal time series model with the phenomenological data. Finally, in Sec. better suggests that the nonlinear oscillators needed to 6 we draw some conclusions. model locomotion operates in the unstable, that is, in the chaotic regime. A stochastic version of a CPG was developed by Haus- II. CENTRAL PATTERN GENERATOR AND dorff et al. [6, 12] to capture the fractal properties of the LOCOMOTION inter-stride interval time series. This stochastic model was later extended by Ashkenazy et al. [17] to describe Walking consists of a sequence of steps. These steps the changing of gait dynamics as we develop from child- may be partitioned into two phases: a stance phase and hood to adulthood. The model is essentially a random a swing phase. The stance phase is initiated when a foot walk on a correlated chain, where each node of the chain strikes the ground and ends when it is lifted. The swing is a neural center of the kind discussed above, and with a phase is initiated when the foot is lifted and ends when different frequency. This random walk is found to gener- it strikes the ground again. The time to complete each ate a fractal process, with a multifractal width that de- phase varies with the stepping speed. A stride interval is pends parametrically on the range of the random walker’s the length of time from the start of one stance phase to step size. Ashkenazy et al. [17] focused on explaining the the start of the next stance phase. changes in the gait time series during maturation, using Traditionally the legged locomotion of animals is un- the stochastic CPG model. derstood through the use of a CPG, an intraspinal net- Herein we extend the previous models by assuming work of neurons capable of producing a syncopated out- that gait dynamics are regulated by a stochastic cor- put [8, 14]. The implicit assumption in such an interpre- related CPG similar to that of Ashkenazy et al. [17], tation is that a given limb moves in direct proportion to coupled to the nonlinear oscillators needed to model lo- the voltage generated in a specific part of the CPG. Ex- comotion in the unstable, forced and chaotic regimes. periments establishing the existence of a CPG have been We show that two parameters, the average frequency done on animals with spinal cord transections. Walking, f0 and the intensity A of the forcing component of the for example, in a mesencephalic cat, a cat with its brain nonlinear oscillator, are sufficient to determine both the stem sectioned rostral to the superior colliculus, is very fractal and multifractal variability of human gait under close to normal, on a flat, horizontal surface, when a sec- normal, stressed and metronomic conditions, using the tion of the midbrain is electrically stimulated. Stepping SCPG model. continues as long as a train of electrical pulses is used to drive the stepping. This is not a simple linear response process because changing the frequency of the driver has III. HUMAN GAIT ANALYSIS little effect of the walking cycle [15]. However, since the frequency of the stepping increases in proportion to the In this section we summarize the main fractal and amplitude of the stimulation, we can conclude that the multifractal characteristics of the stride interval of the variation in the stride interval of humans is related to the human gait data discussed in detail elsewhere [7]. We fluctuation of the amplitude of the impulses of the firing have analyzed the gait time series of 10 persons in the neural centers. three different conditions of slow, normal and fast walk- As Collins and Richmond [8] point out, in spite of the ing. Each time series is approximately one hour long studies establishing the existence of a CPG in the cen- for unconstrained walking, see, for example Fig. 1, for tral nervous system of quadrupeds, such direct evidence slow, fast and normal walking. Similarly, each time se- does not exist for a vertebrate CPG for biped locomotion. ries is approximately 30 minutes long for metronomically Consequently, these and other authors have turned to the constrained walking, see, for example, Fig. 2, for slow, construction of models, based on the coupling of non- fast and normal walking. Participants in the study had linear oscillators, the hard-wired CPG, to establish that no history of any neuromuscular, respiratory or cardio- the mathematical models are sufficiently robust to mimic vascular disorders. They were not taking any medica- the locomotion characteristics observed in the movements tions and had a mean age of 21.7 years (range: 18-29 of segmented bipeds [16], as well as in quadrupeds [9]. years); mean height 1.77 ± 0.08 meters and mean weight These characteristics, such as the switching among mul- 71.8 ± 10.7 kg. All subjects provided informed written 3 consent. Subjects walked continuously on level ground weakly multifractal [7]. However, the multifractal struc- around an obstacle-free, long (either 225 or 400 meters), ture is slightly more prominent for fast and slow gait than approximately oval path and the stride interval was mea- for normal gait. sured using ultra-thin, force sensitive switches taped in- Fig. 5 shows typical distributions of the H¨older expo- side one shoe. For the metronomic constrained walk- nents for metronome-constrained walking, which is little ing, the individuals were told only once, at the begin- different from the histograms in Fig. 3. Fig. 6 shows ning of their walk, to synchronize their steps with the the average distributions of the H¨older exponents for all metronome. More details regarding the collection of data 10 walkers. The figures clearly indicates that under the can be found at Physionet [11] from where the data were constraint of a metronome, the stride interval of human downloaded and in Ref. [7, 12]. gait become more random and the strong long-time per- The fractal and multifractal analysis of the data is done sistence of the 1/f-noise is lost for some individuals. The by studying the estimated distribution of the local H¨older data present a large variability of the H¨older exponents exponents using wavelet transforms [13, 18, 19, 20]. To from persistent to antipersistent fluctuations, that is, the better understand the meaning of the H¨older exponent h exponent spans the entire range of −1 signal for all three gait velocities. multifractal time series can only be distinguished from a In summary, the stride interval of human gait presents monofractal time series of the same length, if the width a complex behavior that depends on many factors. Walk- of its H¨older-exponent distribution is significantly larger ing is a strongly correlated neuronal and biomechanical than that of a correspondent monofractal time series [7]. phenomenon which may be strongly influenced by two Typical distributions of H¨older exponents, for uncon- different stress mechanisms; (a) a natural stress that in- strained walking of a single individual, are depicted in creases the correlation of the nervous system that regu- Fig. 3. Fig. 4 shows the average distributions of the lates the motion at the changing of the gait regime from a H¨older exponents for the cohort of 10 walkers. Figs. 3 normal relaxed condition, to a consciously forced slower and 4 show that stride interval time series for human or faster gait regime; (b) a psychological stress due to gait are characterized by strong persistent fractal prop- the constraint of following a fixed external cadence such erties very close to that of 1/f-noise, h ≈ 0. However, as a metronome. The metronome breaks the long-time normal gait is usually slightly less persistent than both correlation of the natural pace and generates a large frac- slow and fast gait. The slow gait has the most persis- tal variability of the gait regime. In the next section we tent fluctuations and may present non-stationary prop- present a multifractal CPG model able to reproduce these erties, h> 0. The slow gait fluctuation may deviate most properties. strongly from person to person. The higher values of the H¨older exponents for both slow and fast gait, relative to IV. THE SCPG MODEL FOR HUMAN GAIT normal gait, may be explained as due to a stress con- dition that increases the persistency and, therefore, the long-time correlation of the fluctuations. Moreover, the In this section we introduce a model of locomotion that regular curves of Fig. 4 show that unconstrained walking governs the stride interval time series for human gait. As is characterized by fractal properties that do not change anticipated in the previous sections the model has to sim- substantially from one individual to another. Finally, a ulate a CPG [8] capable of producing a syncopated corre- careful comparison of the widths of the distributions of lated output associated with a motocontrol process of the H¨older exponents for the different gaits with the widths gait cycle. Moreover, the model incorporates two sepa- for a corresponding monofractal noise of the same length rate and distinct stress mechanisms. One stress mecha- has proven that the stride interval of human gait is only nism, that has an internal origin, increases the correla- 4 tion of the time series due to the change in the velocity of we can mimic the biped’s gait with only one nonlinear os- the gait from normal to the slower or faster regimes. The cillator. In our model we use a well-known neuronal os- second stress mechanism, has an external origin, and de- cillator model, that is, the forced Van der Pol oscillator creases the long-time correlation of the time series under [8, 10] that is defined by the following equation the frequency constraint of a metronome. We model this 2 2 2 complex phenomenon by assuming that the intensity of x¨ + µ x − p x˙ + (2πfj) x = A sin (2πf0 t) . (1) the impulses of the firing neural centers regulate only the
inner virtual frequency of a forced Van der Pol oscillator The parameter p controls the amplitude of the oscilla- [10]. The observed stride interval is assumed to coincide tions, µ controls the degree of non linearity of the oscil- with the actual period of each cycle of the Van der Pol os- lator, fj is the inner virtual frequency of the oscillator cillator; a period that depends on the unperturbed inner during the j − th cycle that is related to the intensity frequency of the oscillator, the amplitude of the forcing of the the j − th neural fired impulse, and A and f0 are function and the frequency of the forcing function. respectively the strength and the frequency of the exter- nal driver. The frequency of the oscillator would be f if Since the frequency of the stepping increases in propor- A = 0. tion to the amplitude of the electric stimulation [15], we We notice that the non linear term, as well as the can assume that the time series of the intensity of the im- driver, induce the oscillator to move around a limit cycle. pulses fired by the neural centers is associated to a time The actual frequency of each cycle may differ from the series of a virtual frequencies{f }. So, in the spirit of j inner frequency f. We assume that at the conclusion of the model suggested by Ashkenazy et al. [17], we assume each cycle, a new cycle is initiated with a new inner fre- that the long-time correlated frequency of the SCPG is quency f produced by the stochastic CPG model while described by a random walk on a finite-size correlated j all other parameters are kept constant. However, the chain, where each node of the chain is a neural center simulated stride interval is given by the actual period of of the kind discussed above, that fires an impulse with each cycle of the Van der Pol oscillator. a particular intensity that is associated to a paticular We assume that the neural centers of the SCPG may virtual frequency. Ashkenazy et al. [17] focused on ex- fire impulses with different amplitudes that induce vir- plaining the multifractal changes in the gait time series tual frequencies {f } with finite-size correlations. Here, during maturation from childhood to adulthood, assum- i therefore, we model directly the time series of virtual fre- ing that neural maturation is parametrically associated quencies. The frequencies {f } are centered around the with the range C of the Brownian process that activates i driver frequency f0 according to the relation the nodes of the finite-size correlated chain of frequencies. Here, we adopt a different approach because we are in- fi = f0 + γ Xi (2) terested in modeling the gait for human adults operating under different conditions. We assume neural maturation where γ is a constant and Xi is a finite-size correlated and, therefore, the standard deviation C of the random variable, that is, walk process remains constant, whereas the strength of the correlation among the neural centers increases with r 2 = exp − . (3) the change of the velocity of the gait from the normal  r0  to the slower or faster regimes. The change of velocity is interpreted as a biological stress. Moreover, contrary The parameter r0 measures the spatial range of the cor- to Ashkenazy et al. [17] we do not add any noise to relations of the neural network. The chain of frequencies the output of each node to mimic biological noise. The fi = f0 +γ Xi is generated by a first-order autoregressive final output given by the actual frequencies of the gait process, also known as a linear Markov process [22], that cycle fluctuates due to the chaotic solutions of the non- is generated by the recursion equation linear oscillators in the SCPG, here that being the forced Van der Pol oscillator. The advantage of using chaos in Xi = a Xi−1 + εi , (4) the model, rather than noise, is that chaos is an intrinsic where 0 r the change of phase from the walk, trot, canter and gallop ρ(r)= 2 = a . (5) of the quadrupeds, but also the variability of the stride intervals observed in humans. In bipeds it is possible to A direct comparison between Eqs. (3) and (5) gives mimic the movements of the two legs with two nonlin- a = exp[−1/r0], so, we can easily generate a data se- ear coupled oscillators. However, because the geometry quence with the desired finite-size correlation value r0. of the bipeds’ gait, contrary to that for quadrupeds, is Following Ref. [17], we assume that a frequency is acti- unique and the two legs must be shifted by π rad in phase, vated by the position of a random walker given by the 5 function g(j) with j =1, 2,... , whose jump sizes follow a V. SIMULATED STRIDE INTERVAL GAIT Gaussian distribution of width C. The width of this dis- tribution, according to the interpretation of Ashkenazy In this section we present and comment on our com- et al. [17], is associated with the human neural age matu- puter simulations of the stride interval of human gait un- ration. This random walk mechanism allows us to obtain der a variety of conditions. For simplicity, we make use from the finite-time, correlated frequency series {fi}, a of the following values of the parameters. The frequency new time series of frequencies {fj} with i = g(j), charac- of normal gait is fixed at the experimentally determined terized by long-time correlations. Finally, the sequence value of f0,n = 1/1.1 Hz, so that the average period of frequencies {fj} is used in Eq. (1) recursively. of the normal gait is 1.1 second; the frequency of the Normal gait, characterized by the frequency f0,n, is slow and fast gait are respectively f0,s =1/1.45 Hz and assumed to occur when the individual is relaxed, conse- f0,f = 1/0.95 Hz, with an average period of 1.45 and quently the correlations between the neuronal centers are 0.95 seconds, respectively, that is similar to experimen- minimum. By implication, whether the gait increases or tally realized slow and fast human gaits shown in Fig. decreases in velocity, the correlations between the neu- 1. ronal centers increases. This increase in the stress is mod- Also the hopping-range parameter is chosen equal to eled by using the short-time correlation parameter r0 of that for adults [17], that is, C = 25 and kept constant. the stochastic CPG by assuming Moreover, we chose r0,n = 25 such that for f0 = f0,n we have r0 = 25, that coincides with the correspondend- 2 ing value found in Ref. [17]. To generate an artificial r0 = r0,n 1+ B (f0 − f0,n) , (6) sequence with a variability compatible to that of the ex- h i perimental sequence, we chose B = 50 in Eq. (6) and, in Eq. (2), γ =0.02 , that is a value compatible to the aver- where r0,n is the short-range correlation among the firing neural centers at the normal frequency gait and B is a age of the standard deviation of all the data analyzed by positive constant that measures the increasing of short- us [7], however, the value of γ may be smaller and may range correlation at the anomalous frequency gait. decrease with an increase in the frequency f0 and/or an increase in the intensity of the forcing amplitude A of Eq. We observe that in our experimental data set we mea- (1). So, we choose a frequency f0, calculate r0 via Eq. sure the stride intervals of the gait. It is true that the (6) and the Markovian parameter a, then we generate a frequency of walking may be associated with a long-time chain of frequencies {fi} via Eqs. (2) and (4) and, finally, correlated neural firing activity that induces virtual pace by using the random walk process to activate a particular frequency, nevertheless the walking is also constrained by frequency of the short-time correlated frequency neural the biomechanical motocontrol cycle that directly con- chain, we obtain the time series of the frequencies {fj } trols movement and produces the pace itself. The SCPG to use in the time evolution of the Van der Pol oscillator. introduced here incorporates both the neural firing ac- For simplicity, we keep constant the two parameters of tivity and the motocontrol constraints. the nonlinear component of the oscillator (1), µ = 1 and In summery, the SCPG model is based on the cou- p = 1. The only parameters allowed to change in the pling of a stochastic with a hard-wired CPG model and model are the mean frequency f0 that changes also the depends on many factors. The most important param- value of r0 via Eq. (6), and the intensity A of the driver eters are the short-correlation size r0 of Eq. (3), that of the Van der Pol oscillator (1). measures the correlation between the neuron centers of Fig. 7 shows the stride interval time series for slow, the stochastic CPG, the intensity A of the forcing driving normal and fast computer-simulated gaits using SCPG. component of the nonlinear oscillator of Eq. (1) and, of For the simulation of normal gait we use A = 1 and for course, the average frequency f0 of the actual pace that both slower and faster gait, we use A = 2. We suppose distinguish the slow, normal and fast gait regimes. The that the amplitude A of the driver of the Van der Pol other parameters, γ, C, µ and p may be, to a first-order oscillator (1) should be smaller for normal gait than that approximation, kept fixed. for either the slower or faster gaits, because in our inter- While the numerical simulations are left to the follow- pretation A measures the magnitude of the constraint to ing section, we can anticipate an interpretation of the two walk at a particular velocity. The amplitude A is smaller main parameters r0 and A. In fact, the short-correlation for normal gait because normal gait is the most relaxed, size r0 may be interpreted as a parameter that mea- spontaneous and, and consequently the most automatic sures the natural correlation between the neural centers of the three gaits. The figure shows that the SCPG model and such short-time correlation increases under particu- is able to reproduce a realistic persistence and volatility lar stress, for example, when the velocity of the gait is for the three gaits by simply changing the frequency of slower or faster than the normal relaxed situation. The the gait itself. In particular note the high volatility of intensity of the forcing driving component A may be as- the slow gait that is remarkably similar to that seen in sociated with the voluntary action of trying to follow a Fig.1. particular cadence and is expected to increase under a Fig. 8 shows the stride interval time series for metronomic constraint. slow, normal and fast metronome-triggered computer- 6 simulated gaits. We use the same frequency series gener- the metronome and, therefore, we should expect that the ated by the SCPG used to produce the sequences of Fig. normal gait is the most automatic and the least con- 7. We only change the intensity A of the driver of the strained by human consciousness. This is the reason we Van der Pol oscillator (1). We use for the normal gait have chosen A = 4 for the normal metronome-triggered A = 4 and for both slower and faster gait A = 8. Again gait. For both slower and faster metronome-triggered we suppose that the intensity A of the driver of the Van gaits we have chosen A = 8 to indicate a higher con- der Pol oscillator (1) should be smaller than that for both scious stress, that constrains gait at anomalous speeds. slower and faster gaits, because normal gait is the most Moreover, by comparing Fig. 9 and Fig. 10 and con- relaxed and spontaneous. By comparing Figs. 7 and 8 sidering that in both simulations we have used the same we note the increase in randomness, the loss of persis- value of the forcing parameter A for both slower and tency and the reduction in volatility; all effects that are faster gaits, we notice that the largest fractal shift oc- induced in the latter time series by increasing the value curs for the slower gait. This increased shift implies that of A and are found in the phenomenologic data shown in the slower gait is the more sensitive to a voluntary con- Figs. 1 and 2. straint and, so, and so the slower mode has the larger Fig. 9 shows histograms of distributions of the H¨older variability. In fact, our human experience and the super- exponents for the three computer-simulated gaits shown position of the distributions of H¨older exponents for the in Fig. 7. The calculation are done in the same way of 10 cohorts in Fig. 6, show a large fractal variability of the those used to produce the histograms in Fig. 3 for the slower gait. Finally, Fig. 6 reveals that few persons are experimental data, for details see Ref. [7]. The figure characterized by a strong antipersistent pace when asked shows that the SCPG model is able to generate artifi- to follow a metronome. According to the SCPG model , cial stride interval time series with statistical properties some people are not able to find a natural synchroniza- similar to the fractal and multifractal behavior of real tion and need to continuously adjust and readjust the data. By changing the gait mode from slow to normal, speed of their pace to match the beat of the metronome. the mean H¨older exponent h decreases. In the same way This changing of pace implies a very strong conscious act by changing the gait mode from normal to fast, the mean and, therefore, a very high value of the parameter A that H¨older exponent again increases, just as it does for the would produce a strong antipersistent signal. real data. According to the SCPG model, this increase in the scaling parameter is due to the increase of the in- ner short-time correlation among the neuronal centers, VI. CONCLUSION modeled by Eq. (6). Furthermore, this behavior is due to the biological stress of consciously walking at a speed We have introduced a new kind of CPG model. One that is different from the normal spontaneous speed. In able to mimic the complexity of the stride interval se- addition the multifractality of the gait time series slightly quences of human gait under the several conditions of increases for a walking rate different from normal. Here slow, normal and fast regimes for both walking freely again this effect is observed in the real stride interval data and keeping the beat of a metronome. The SCPG model and it is proven by a slight increase in the width of the is based on the assumption that human locomotion is histograms for fast and in particular slow gait. regulated by both the central nervous system and by a Fig. 10 shows the histograms of probability den- motocontrol system. A network of neurons produces a sity estimations of the H¨older exponents for the three correlated syncopated output that is correlated accord- metronome-triggered computer-simulated gaits shown in ing to the level of physiological stress and this network Fig. 8. The calculated points show that the SCPG is is coupled to the motocontrol process. The combina- able to generate artificial stride interval time series that tion of systems controls locomotion and the variability of present similar fractal and multifractal behaviors to those the gait cycle. It is the period of the gait cycle that is of real stride interval data taken under the constraint of measured in the data sets considered herein. Moreover, a metronome. By increasing the intensity A of the driver walking may be conditioned by a voluntary act as well, of the Van der Pol oscillator (1) the randomness of the for example, walking may be consciously forced following time series increases and it is possible to obtain a large the frequency of a metronome. We model the complex variety of time series, from those having antipersistent system generating the data by assuming that the corre- to those with persistent fractal properties. In the SCPG, lated firing activity of the neural centers regulate only the the parameter A measures the constraint of consciousness inner frequency of a forced Van der Pol oscillator. How- on the gait, and therefore the value of A has to increase ever, it is the forced Van der Pol oscillator that mimics if the walker is asked to synchronize his/her pace with the motocontrol mechanism of the gait cycle. The stride the frequency of a metronome. The figure suggests that interval is the actual period of each cycle of the forced the SCPG model is able to explain a number of other Van der Pol oscillator. In this way the gait frequency is properties of the metronome-triggered walking. Fig. 6 slightly different from the inner frequency induced by the shows that the usually normal metronome-triggered gait neural firing activity whose inpulse intensity are able to is that with the highest persistent fractal properties. Nor- generate only a potential, but not an actual frequency. mal gait is also the most natural under the constraint of The chaotic behavior of such a nonlinear oscillator, like 7 the Van der Pol oscillator, and the possibility to force the constraint. Finally, Ref. [24] reports that the stride in- frequency of the cycle with an external fixed frequency al- terval time series for elderly subjects and for subjects lows the SCPG model to generate time series that present with Huntington’s diseases are more random than for similar fractal and multifractal properties to that of the young healthy subjects. According to the SCPG model, human physiological stride interval data in all situations this may be explained by a decrease of the normal short- here analyzed. Moreover, by implementing the SCPG range correlation among the neural centers that may be with four coupled forced Van der Pol oscillators as in associated with a nervous degeneration caused by injury, Ref. [8], it should be possible to simulate the change of disease, or aging. This decrease in correlation may be phase between various modes of quadrupeds’ locomotion, modeled through r0,n of Eq. (6). However, the decrease that is, walking, trotting, cantering and galloping. of correlation in the gait of those subjects may also be The variety of complex behaviors is regulated by two associated with an increase of the amplitude A of the parameters, the average frequency f0 and the amplitude driving force of the Van der Pol oscillator, Eq. (1). In A of the driver of the Van der Pol oscillator. The fre- fact, those subjects may also consciously choose to walk quency f0 regulates the speed and may be associated with more carefully. a neuronal stress that increases the correlation among the ——– neural centers. The amplitude A may be associated with Acknowledgment: the voluntary action of trying to track a particular fre- N.S. thanks the Army Research Office for support under quency and it is expected to increase under a metronome grant DAAG5598D0002.

[1] J.M. Hausdorff, C.-K.Peng, Z. Ladin, J.Y. Wei and A.L. Control of Posture and Movements, Spring-Verlag, New Goldberger, “Is Walking a Random-walk - Evidence for York 2000. Long-range Correlations in Stride Interval of Human [15] M.D. Mann, The Nervous System and Behavior, Harper Gait,” J. Appl. Physiol. 78 (1), 349-358 (1995). & Row, Philadelphia (1981). [2] J.M. Hausdorff, L. Zemany, C.-K. Peng and A.L. Gold- [16] A.H. Cohen, S. rossignol, and S. Grillner, Editors, Neural berger, “Maturation of Gait Dynamics: Stride-to-stride control of rythmic movements in vertebrates, Wiley, New Variability and its Temporal Organization in Children,”J. York (1988). Appl. Physiol. 86(3), 1040 (1999). [17] Y. Ashkenazy, J.M. Hausdorff, P. Ivanov, A.L. Gold- [3] B.J. West and L. Griffin, “Allometric Control, Inverse berger and H.E. Stanley, “A Stochastic Model of Human Power Laws and Human Gait,” Chaos, Solitions & Frac- Gait Dynamics,” cond-mat/0103119 v1. tals 10 (9), 1519-1527 (1999). [18] S. G. Mallat, A Wavelet Tour of Signal Processing (2nd [4] B.J. West and L. Griffin, “Allometric Control of Human edition), Academic Press, Cambridge (1999). Gait,” Fractals 6 (2), 101-108 (1998). [19] J. F. Muzy, E. Bacry, A. Arneodo, “The Multifractal For- [5] Vierordt, Ueber das Gehen des Meschen in Gesunden und malism Revisited with Wavelets,” Int. J. Bifurc. Chaos 4, Kranken Zustaenden nach Selb-stregistrireden Methoden, No. 2, 245-302 (1994). Tuebigen, Germany (1881). [20] P. C. Ivanov, M. G. Rosenblum, L. A. Nunes Amaral, Z. [6] L. Griffin, D.J. West and B.J. West, “Random stride in- R. Struzik, S. Havlin, A. L. Goldberger and H. E. Stanley, tervals with memory” J. Biol. Phys. 26, 185-202 (2000). “Multifractality in Human Heartbeat Dynamics,” Nature [7] N. Scafetta, L. Griffin and B. J. West, “Holder expo- 399, 461-465 (1999). nent spectra for human gait,” unpublished (2002), cond- [21] B.B. Mandelbrot, The Fractal Geometry of Nature, Free- mat/0208028. man, New York, (1983). [8] J.J. Collins and S.A. Richmond, “Hard-wired Central [22] M.B. Priestley, Spectral analysis and time series, Aca- Pattern Generators for Quadrupedal Locomotion,” Biol. demic Press, UK 2001. Cyb. 71 (5), 375-385 (1994). [23] J.B. Bassigthwaighte, L.S. Liebowitch and B.J. West, [9] J.J. Collins and I.N. Stewart, “Coupled Nonlinear Oscil- Fractal Physiology, Oxford University Press, Oxford lators and the Symmetries of Animal Gaits,” J. Nonlinear (1994). Sci. 3, 349-392 (1993). [24] J.M. Hausdorff, Y. Ashkenazy, C-K. Peng, P.C. Ivanov, [10] J. Guckenheimer and P. Holmes, Nonlinear Oscillations, H.E. Stanley and A.L. Goldberger, “When Human Walk- Dynamical Systems, and Bifurcations of Vector Fields, ing Becomes Random Walking: Fractal Analysis and Springer-Verlag, New York USA 2002. Modeling of Gait Rhythm Fluctuations,” Physica A 302 [11] http://www.physionet.org/ . 138-147 (2001). [12] J.M. Hausdorff, P.L. Purdon, C.-K. Peng, Z. Ladin, J.Y. Wei, A.L. Goldberger, “Fractal Dynamics of Hu- man Gait: Stability of Long-range Correlations in Stride Interval Fluctuations,” J. Appl. Physiol. 80, 1448-1457, (1996). [13] Z. R. Struzik, “Determining Local Singularity Strengths and their Spectra with the Wavelet Transform,” Fractals, Vol. 8, No. 2, 163-179 (2000). [14] J. M. Winters and P. E. Crago, Biomechanics and Neural 8

1.7 slow gait 1.6 1.5 1.4 1.3 1.2 normal gait 1.1

stride interval (sec) 1 0.9 fast gait 0.8 0 10 20 30 40 50 60 time (min)

FIG. 1: Stride interval for slow, normal and fast gait. The period of time over which measurements were done is approximately one hour.

1.6 slow gait 1.5

1.4

1.3

1.2 normal gait 1.1 stride interval (sec) 1 fast gait

0.9 0 5 10 15 20 25 30 time (min)

FIG. 2: Stride intervals for slow, normal and fast gait for metronomic triggered walking. The total period of time is approximately 30 minutes. 9

8 fast gait 7 normal gait 6 slow gait 5 4

- p(h) 3 2 1 0 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 Holder exponent - h -

FIG. 3: Histogram and probability density estimation of the H¨older exponents: slow-star (h0 = 0.105, σ = 0.060), normal- triangle (h0 = −0.125, σ = 0.063) and fast-circle (h0 = −0.012, σ = 0.056) gait for a single individual. The fitting curves are Gaussian functions with average h0 and standard deviation σ. 10

6 Slow Gait

5 Normal Gait Fast Gait 4

3 - p(h) 2

1

0 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 Holder exponents - h -

FIG. 4: Histogram and probability density estimation of the H¨older exponents for the three walking groups are shown: slow- star, normal-triangle and fast-circle gait. Each curve is an average over the 10 cohorts in the experiment. By changing the gate mode from slow to normal the Holder exponents h decrease but from normal to fast they increase. There is also an increasing of the width of the distribution σ by moving from the normal to the slow or fast gait mode. The fitting curves are Gaussian functions: slow-star (h0 = 0.046, σ = 0.102), normal-triangle (h0 = −0.092, σ = 0.069) and fast-circle (h0 = −0.035, σ = 0.081) gait 11

8 slow fast norm 7 6 5 4 - p(h) 3 2 1 0 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 Holder exponent - h -

FIG. 5: Metronomic walking for a single individual. Histogram and probability density estimation of the H¨older exponents: slow-star (h0 = −0.765, σ = 0.064), normal-triangle (h0 = −0.204, σ = 0.064) and fast-circle (h0 = −0.436, σ = 0.066).

3.5 Fast gait 3 Normal gait Slow gait 2.5

2

1.5 - p(h)

1

0.5

0 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 Holder exponent - h -

FIG. 6: Metronomic walking. Histogram estimation of the H¨older exponents for the three walking groups: slow-star, normal- triangle and fast-circle gait. Each curve is an average over the 10 cohorts in the experiment. 12

1.7 slow gait 1.6 1.5 1.4 1.3 1.2 norm gait 1.1

stride interval (sec) 1 0.9 fast gait 0.8 0 10 20 30 40 50 60 time (min)

FIG. 7: Stride interval time series for slow, normal and fast computer-simulated gaits.

1.6 slow gait 1.5

1.4

1.3

1.2 norm gait

1.1 stride interval (sec)

1 fast gait

0.9 0 10 20 30 40 50 60 time (min)

FIG. 8: Stride interval time series for slow, normal and fast gait for metronome-triggered computer-simulated gait. 13

8 slow gait 7 norm gait 6 fast gait 5 4

- p(h) 3 2 1 0 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 Holder exponent - h -

FIG. 9: Histogram of probability density estimation of the H¨older exponents for computer-simulated gait: slow-star (h0 = 0.058, σ = 0.068), normal-triangle (h0 = −0.093, σ = 0.058) and fast-circle (h0 = −0.015, σ = 0.063) gait for a single individual.

8 slow gait 7 norm gait 6 fast gait 5 4

- p(h) 3 2 1 0 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 Holder exponent - h -

FIG. 10: Histogram and probability density estimation of the H¨older exponents for metronome-triggered computer-simulated gait: slow-star (h0 = −0.516, σ = 0.067), normal-triangle (h0 = −0.276, σ = 0.059) and fast-circle (h0 = −0.373, σ = 0.063) gait for a single individual.