# Patent

## Title
Control system for walking assist device

## Abstract

A system capable of controlling an autonomous walking assist rhythm, while following changes in a walking rhythm of a walker at the same time. A control system 100 in accordance with the present invention generates a first oscillator x that tugs and is tugged by a walking oscillator (hip joint angular velocity) φH′ such that an intrinsic angular velocity is reflected. Then, a new intrinsic angular velocity ωm is determined on the basis of the difference between a first phase difference θHM between the first oscillator x and the walking oscillator φH′ and a desired phase difference θd. Furthermore, a second oscillator y is generated that tugs and is tugged by the walking oscillator φH′ such that the intrinsic angular velocity ωm is reflected, and has a second phase difference θhm, which is closer to the desired phase difference θd than the first phase difference θHM is, with respect to the walking oscillator φH′. Then, based on the second oscillator y and the walking oscillator (hip joint angle) φH, a walking assist oscillator (torque imparted to the walker by a walking assist unit 200) T is generated.

## Background
TECHNICAL FIELD
The present invention relates to a control system for a walking assist device.
BACKGROUND ART
There has been known a device adapted to assist a walker with his/her walking by imparting torques about leg joints (hip joints, knee joints, and foot joints) to the walker.
In the walking assist device, if a walking assist rhythm (a rhythm for imparting torque to a walker) is programmed beforehand, then the walker is forced to follow the walking rhythm based on the programmed walking assist rhythm. Therefore, if the walking assist rhythm is slow, then the walker has to walk slowly even if the walker wishes to walk faster. Conversely, if the walking assist rhythm is fast, then the walker has to walk fast even if the walker wishes to walk slowly. Hence, the walker is very likely to develop uncomfortableness, feeling that the walker is forced to walk in a manner against his/her will.
A walking assist rhythm arranged to fully follow a walking rhythm of a walker enables the walker to walk at a walking rhythm based on his/her own will.
In other words, if the walking assist rhythm is controlled to be faster in response to a change in the movement of legs when the walker moves his/her legs faster, speeding up the walking rhythm, then the walker will be able to walk faster as he/she intends to. Similarly, if the walking assist rhythm is controlled to be slower in response to a change in the movement of the legs when the walker moves his/her legs slowly, slowing down the walking rhythm, then the walker will be able to walk slowly as he/she intends to.
On the other hand, it is important for a walker wearing the walking assist device to feel that his/her walking is assisted since he/she is using the walking assist device. Hence, in walking assistance fully led by a walker, the walker may feel uncomfortable with the walk using the walking assist device due to such feeling of not being assisted, even if the walker can walk at an intended rhythm.
If a walking rhythm generated by a walker himself/herself suddenly changes or if the walking rhythm is unnatural from normal walk, then the walking based on the walking rhythm is furthered, possibly leading to an increase in burden to the mind and body of the walker.
Accordingly, the present invention is intended to achieve a solution by providing a system capable of controlling an autonomous walking assist rhythm, while following changes in a walking rhythm of a walker at the same time.
DISCLOSURE OF INVENTION
A control system for a walking assist device in accordance with the present invention to solve the aforesaid problem is characterized by comprising a walking oscillator measuring means for measuring a walking oscillator of a walker whose walk is assisted by a walking assist device, a first oscillator generating means for generating a first oscillator that tugs and is tugged by the walking oscillator measured by the walking oscillator measuring means such that an intrinsic angular velocity is reflected, a determining means for determining a new intrinsic angular velocity on the basis of a difference between a first phase difference between the first oscillator and the walking oscillator and a desired phase difference, a second oscillator generating means for generating a second oscillator that tugs and is tugged by the walking oscillator measured by the walking oscillator measuring means such that the intrinsic angular velocity determined by the determining means is reflected, and has a second phase difference, which is closer to the desired phase difference than the first phase difference is, with respect to the walking oscillator, and a walking assist oscillator generating means for generating a walking assist oscillator of the walking assist device on the basis of the second oscillator and the walking oscillator measured by the walking oscillator measuring means.
According to the present invention, it is possible to generate a oscillator having a proper phase difference from the viewpoint of walking assistance with respect to a walking oscillator in which a walking rhythm is reflected, and then a walking assist oscillator is generated on the basis of the oscillator, thereby allowing a walking assist rhythm to be optimized.
The “oscillator” used in the present invention means a realistic or virtual parameter that vibrates (temporally changes) at a certain rhythm (angular velocity). For example, the “walking oscillator” means a leg joint angle or a leg joint angular velocity or the like of a walker that vibrates at a rhythm in which a walking rhythm is reflected. The “walking assist oscillator” means an imparted torque or the like around leg joints applied to the walker that vibrates with a rhythm in which a walking assist rhythm is reflected.
More specifically, first, {circle around (1)} “the first oscillator” that tugs and is tugged by a walking oscillator such that “the intrinsic angular velocity” is reflected is generated. The first oscillator vibrates with an autonomous rhythm in which “the intrinsic angular velocity” is reflected while harmonizing with the rhythm of the walking oscillator due to the “mutual tugging” effect.
However, although the first oscillator harmonizes with the rhythm of the walking oscillator, the first phase difference relative to the walking oscillator is not necessarily a phase difference appropriate for assisting walk. Hence, if a walking assist oscillator is determined on the basis of the first oscillator, then a phase difference ( ̃a first phase difference) inappropriate for assisting walk may be generated between a walking rhythm and a walking assisting rhythm, leading to inappropriate walking with the walking assistance. Accordingly, it is necessary to generate a new oscillator having a phase difference appropriate for walking assistance relative to the walking oscillator on the basis of the first oscillator.
Hence, {circle around (2)} a new “intrinsic angular velocity” is determined, and then “a second oscillator” is generated that tugs and is tugged by the walking oscillator such that the intrinsic angular velocity is reflected, and has “a second phase difference,” which is closer to “a desired phase difference” than the first phase difference is, with respect to the walking oscillator.
The rhythm of the second oscillator has a phase difference ( ̃a second phase difference) that is “a desired phase difference” or one close thereto relative to a walking rhythm ( ̃rhythm of the walking oscillator). Thus, if the walking assist oscillator is determined on the basis of the second oscillator, then the phase difference between the walking rhythm and the walking assist rhythm approaches to the desired phase difference.
Then, {circle around (3)} a walking assist oscillator is generated on the basis of the second oscillator and the walking oscillator. This generates the walking assist oscillator so as to realize a walking rhythm and the desired phase difference. Furthermore, the walk of the walker can be assisted at a rhythm based on a walking condition reflected in the walking oscillator and the walking rhythm.
Accordingly, the present invention makes it possible to achieve a walking assist rhythm that {circle around (1)} follows changes in the walking rhythm ( ̃walking oscillator rhythm) of a walker, while {circle around (2)} maintaining autonomy that has the walking rhythm and “a desired phase difference,” and {circle around (3)} achieving smooth walking assistance according to a walking condition of the walker.
As described above, a walking assist rhythm harmonizes with a change in a walking rhythm when the walking rhythm changes, and the walking rhythm harmonizes with the walking assist rhythm, thus allowing the walker (human body) and a walking assist device (machine) to harmonize with each other (mutual concession). This arrangement makes it possible to achieve appropriate walking assistance, while allowing a walker to reasonably feel that his/her walk is being assisted by the walking assist device. Moreover, even if a walking rhythm suddenly changes, the walking assist rhythm does not fully follow the sudden change, making it possible to avoid furthering walking assistance and walking that cause burden to the mind and body of the walker.
The present invention is characterized by being further equipped with a walking condition determining means for determining the walking condition of the walker on the basis of a walking oscillator measured by the walking oscillator measuring means, and a desired phase difference determining means for determining a desired phase difference on the basis of the walking condition of the walker determined by the walking condition determining means.
The present invention is characterized by being further equipped with a physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a walking condition determining means for determining the walking condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and a desired phase difference determining means for determining a desired phase difference on the basis of a walking condition of the walker determined by the walking condition determining means.
According to the present invention, it is possible to achieve walk with proper walking assistance responsive to changes in a walking condition by bringing a phase difference between a walking rhythm and a walking assist rhythm close to a desired phase difference that changes with the walking condition of the walker.
The present invention is characterized in that the walking condition determining means determines whether a walking condition of a walker is a flat walk condition, an ascending walk condition, or a descending walk condition, and the desired phase difference determining means basically determines a desired phase difference to be 0 for the flat walk condition, − for the ascending walk condition, and + for the descending walk condition.
According to the present invention, when the desired phase difference is determined to be +, the walking rhythm and the walking assist rhythm will have a phase difference of the determined value (>0) or a phase difference close thereto, enabling the walker to walk while leading the walking assist device. Conversely, when the desired phase difference is determined to be −, the walking rhythm and the walking assist rhythm will have a phase difference of the determined value (<0) or a phase difference close thereto, enabling the walker to walk while being led by the walking assist device.
Thus, the walker can walk by leading the walking assist device in the descending walk condition wherein burden is relatively small, while by being led by the walking assist device in the ascending walk condition wherein burden is relatively large. The desired phase difference can be changed according to walker's intention or other factors rather than being fixed to 0, + or − according to a walking condition.
Furthermore, the present invention is characterized in that the walking condition determining means determines the speed of walk as a walking condition, and the desired phase difference determining means basically determines a desired phase difference to be − if the walk is fast, or + if the walk is slow.
The present invention enables a walker to walk by leading the walking assist device for “slow walk” in which burden is relatively small, while by being led by the walking assist device for “fast walk” in which burden is relatively large.
The present invention is characterized by comprising a physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a physiological condition determining means for determining the physiological condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and a desired phase difference determining means for determining a desired phase difference on the basis of a physiological condition of the walker determined by the physiological condition determining means.
According to the present invention, it is possible to achieve walk with proper walking assistance responsive to changes in a physiological condition by bringing a phase difference between a walking rhythm and a walking assist rhythm close to a desired phase difference that changes with the physiological condition of the walker.
The present invention is characterized in that the physiological condition determining means determines fatigue of a walker as a physiological condition of the walker, and basically determines a desired phase difference to be + if the fatigue of the walker is small, or − if the fatigue of the walker is large.
The present invention enables a walker to walk by leading the walking assist device if his/her fatigue is small, or by being led by the walking assist device if his/her fatigue is large.
Furthermore, according to the present invention, the first oscillator generating means generates a plurality of first oscillators corresponding to the motions of a plurality of virtual objects such that correlation among the plurality of the objects is reflected, and the second oscillator generating means generates a plurality of second oscillators corresponding to the motions of the plurality of virtual objects such that correlation among the plurality of the objects is reflected.
The present invention makes it possible to generate the first and second oscillators so that walking assist oscillators can be generated to further fit actual walk of a walker by adjusting an intensity level of correlation among virtual plural objects.
More specifically, if, for example, right and left legs or a plurality of joints of the same leg are assumed as virtual plural objects, then “the first oscillator” and “the second oscillator” and furthermore “a walking assist oscillator” are generated such that a qualitative correlation between the right and left legs, such as alternate, longitudinal motion, and qualitative correlation or the like among joints of the same leg, such as periods, phase differences, etc. of leg motions about hip joints and leg motions about knee joints, are reflected. With this arrangement, a generated walking assist rhythm based on a generated walking assist oscillator can be adjusted to a rhythm appropriate for a walker on the basis of the qualitative correlation.
The present invention is characterized by comprising a walking condition determining means for determining a walking condition of a walker on the basis of a walking oscillator measured by a walking oscillator measuring means and a first correlation adjusting means for adjusting a correlation among a plurality of virtual objects related to generation of first and second oscillators on the basis of a walking condition of the walker determined by the walking condition determining means.
Furthermore, the present invention is characterized by comprising a physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a walking condition determining means for determining the walking condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and a first correlation adjusting means for adjusting a correlation among a plurality of virtual objects related to the generation of the first and second oscillators on the basis of the walking condition of the walker determined by the walking condition determining means.
According to the present invention, a walking assist rhythm can be properly adjusted on the basis of a qualitative correlation of right and left legs or the like assumed as the plurality of virtual objects, as described above, and walking assistance at an optimum rhythm based on a walking condition can be achieved in real time such that a circumstance in which the qualitative correlation changes as “the walking condition” of the walker changes is reflected.
The present invention is characterized by comprising a physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a physiological condition determining means for determining the physiological condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and a first correlation adjusting means for adjusting a correlation among a plurality of virtual objects related to the generation of the first and second oscillators on the basis of the physiological condition of the walker determined by the physiological condition determining means.
According to the present invention, a walking assist rhythm can be properly adjusted on the basis of a qualitative correlation of right and left legs or the like assumed as the plurality of virtual objects, as described above, and walking assistance at an optimum rhythm based on a physiological condition can be achieved in real time such that a circumstance in which the qualitative correlation changes as “the physiological condition” of the walker changes is reflected.
The present invention is characterized in that a determining means determines an intrinsic angular velocity for bringing a second phase difference between a second oscillator and a walking oscillator closer to a desired phase difference than a first phase difference between a first oscillator and a walking oscillator is, according to an internal model representing a phase relationship between a virtual walking oscillator and a virtual walking assist oscillator.
The present invention is characterized in that a determining means is provided with a correlation factor determining means for determining correlation factors of a virtual walking oscillator and a virtual walking assist oscillator such that a second phase difference of the virtual walking oscillator and the virtual walking assist oscillator approaches to a first phase difference according to an internal model, a first angular velocity determining means for determining an angular velocity of the virtual walking oscillator on the basis of the correlation factors, and a second angular velocity determining means for determining an angular velocity of the virtual walking assist oscillator as an intrinsic angular velocity related to the generation of the second oscillator such that the second phase difference approaches to a desired phase difference on the basis of the angular velocity of the virtual walking oscillator.
According to the present invention, a phase difference ( ̃a second phase difference) between a virtual walking rhythm and a virtual walking assist rhythm is brought closer to “a desired phase difference” than a first phase difference is. This makes it possible to bring a phase difference between an actual walking assist rhythm and an actual walking rhythm close to “the desired phase difference” suited to the walking provided with walking assistance.
The present invention is characterized in that the walking assist oscillator generating means generates a walking assist oscillator corresponding to behaviors of a plurality of virtual neural elements such that a correlation among the plurality of virtual neural elements is reflected.
According to the present invention, it is possible to achieve walking assistance that is more suitable to actual walking of a walker by adjusting the intensity of correlation among a plurality of virtual neural elements.
Furthermore, the present invention is characterized by comprising a walking condition determining means for determining a walking condition of a walker on the basis of a walking oscillator measured by a walking oscillator measuring means, and a second correlation adjusting means for adjusting a correlation among a plurality of virtual neural elements related to the generation of a walking assist oscillator on the basis of the walking condition of the walker determined by the walking condition determining means.
Furthermore, the present invention is characterized by comprising a physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a walking condition determining means for determining the walking condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and a second correlation adjusting means for adjusting a correlation among a plurality of virtual neural elements related to the generation of the a walking assist oscillator on the basis of the walking condition of the walker determined by the walking condition determining means.
According to the present invention, a walking assist rhythm can be properly adjusted by adjusting the intensity of correlation among the plurality of virtual neural elements, as described above, and walking assistance at an optimum rhythm based on a walking condition can be achieved in real time such that a circumstance in which the intensity of correlation changes as “the walking condition” of the walker changes is reflected.
The present invention is characterized by comprising a physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a physiological condition determining means for determining the physiological condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and a second correlation adjusting means for adjusting a correlation among a plurality of virtual neural elements related to the generation of the walking assist oscillator on the basis of the walking condition of the walker determined by the walking condition determining means.
According to the present invention, a walking assist rhythm can be properly adjusted by adjusting the intensity of correlation among the plurality of virtual neural elements, as described above, and walking assistance at an optimum rhythm based on a physiological condition can be achieved in real time such that a circumstance in which the intensity of correlation changes as “the physiological condition” of the walker changes is reflected.
The present invention is characterized by comprising a storing means for storing a correspondence relationship between a walking condition and a tracing pattern drawn in an n-dimensional space by an n number of walking oscillators, wherein the walking condition determining means determines a walking condition on the basis of the correspondence relationship stored by the storing means and the tracing pattern drawn in the n-dimensional space by the n number of walking oscillators measured by the walking oscillator measuring means.
According to the present invention, a walking condition can be determined with high accuracy on the basis of a certain relationship (found by the inventors of the application concerned) between a walking condition and a trace drawn by the n number of walking oscillators in the n-dimensional space.
Furthermore, the present invention is characterized in that the first and second oscillator generating means generate the first and second oscillators according to a nonlinear oscillator model represented by the Van der Pol equation that includes a walking oscillator measured by the walking oscillator measuring means.
According to the present invention, mutual tugging between walking oscillators and the first and second oscillators can be adjusted by adjusting a term included in the Van der Pol equation.
Furthermore, the present invention is characterized in that the walking oscillator measuring means measures the first oscillators, the second oscillators, and at least one of various oscillators including a joint angle and angular velocity of a walker that periodically changes at a rhythm corresponding to a walking rhythm, as a walking oscillator for generating a walking assist oscillator.
According to the present invention, the first oscillators, the second oscillators, and the walking assist oscillators can be generated on the basis of a joint angular velocity and angles (walking oscillators) that vibrate such that an actual walking rhythm is accurately reflected.
The present invention is characterized in that the walking assist oscillator generating means generates, as a walking assist oscillator, a torque about a leg joint imparted to a walker, or a oscillator that periodically changes in a form perceivable by a walker in response to a torque and a change in torque.
According to the present invention, a torque or the like about a leg joint can be imparted to a walker at a rhythm harmonized with a walking rhythm of the walker with a desired phase difference.
The present invention is characterized in that the walking oscillator measuring means measures, as a walking oscillator for determining a walking condition, at least one of a leg joint angle, a leg joint angular velocity, a leg joint angular acceleration, and the position of a part of a leg of a walker.
According to the present invention, a walking condition is determined on the basis of a leg joint angular velocity or the like closely related to the walking condition, permitting improved accuracy of the determination to be achieved.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a constructional explanatory diagram of an embodiment of a control system in accordance with the present invention.
FIG. 2 is a functional explanatory diagram of the embodiment of the control system in accordance with the present invention.
FIG. 3 is an explanatory diagram related to a correlation between two oscillators.
FIG. 4 is an explanatory diagram (part 1) of a walking assist function according to the embodiment of the control system in accordance with the present invention.
FIG. 5 is an explanatory diagram (part 2) of a walking assist function according to the embodiment of the control system in accordance with the present invention.
FIG. 6 is an explanatory diagram (part 3) of a walking assist function according to the embodiment of the control system in accordance with the present invention.
FIG. 7 is an explanatory diagram (part 4) of a walking assist function according to the embodiment of the control system in accordance with the present invention.
FIG. 8 is an explanatory diagram (part 5) of a walking assist function according to the embodiment of the control system in accordance with the present invention.
FIG. 9 is an explanatory diagram (part 6) of a walking assist function according to the embodiment of the control system in accordance with the present invention.
FIG. 10 is an explanatory diagram of a neural oscillator model of another embodiment of the control system in accordance with the present invention.



BEST MODE FOR CARRYING OUT THE INVENTION
An embodiment of the control system for a walking assist device in accordance with the present invention will be explained in conjunction with the accompanying drawings.
In the following explanation, apostrophes (′) and double-apostrophes (′′) attached to variables (except for circled numerals) represent time-dependent one time differentiation and two time differentiation, respectively, of the variables. Subscripts L and R will be attached to parameters in order to distinguish left and right legs or the like of a walker. However, the subscripts L and R will be omitted, as appropriate, when there is no need to distinguish right and left.
A walking assist device 200 to be controlled by a control system 100 shown in FIG. 1 is equipped with supporters 202 attached to the abdomen and right and left thighs of a walker, rigid members 204 made integral with the supporters 202, actuators 210 located on sides of the waist of the walker to impart torques about hip joints to the walker through the supporters 202 and the rigid members 204, angle sensors 220 that output signals based on the angles of the right and left hip joints, and a battery 240 built in a backpack 230 carried on the back of the walker to serve as a power source of the actuators 210, etc.
The control system 100 is constructed of a CPU (arithmetic processing unit), a ROM, a RAM or other types of memories (not shown) built in the backpack 230. The control unit 100 is further provided with a walking oscillator measuring unit 110, a first oscillator generating unit 120, a first phase difference measuring unit 130, a walking condition determining unit 140, a desired phase difference determining unit 150, a determining unit 160, a second oscillator generating unit 170, a neural oscillator generating unit 180, and a walking assist oscillator generating unit 190.
The walking oscillator measuring unit 110 includes a hip joint angular velocity measuring unit 112, and a hip joint angle measuring unit 114. The hip joint angular velocity measuring unit 112 measures a hip joint angular velocity φH′ on the basis of an output of the angle sensor 220. The hip joint angle measuring unit 114 measures a hip joint angle φH on the basis of an output of the angle sensor 220.
The first oscillator generating unit 120 generates first oscillators (to be discussed later) x that tug with each other such that the hip joint angular velocity φH′ that changes with an angular velocity ωH measured by the hip joint angular velocity measuring unit 112 and an intrinsic angular velocity ωM are reflected. The first oscillators x are generated according to a nonlinear oscillator model expressed by a Van der Pol equation, which will be discussed later.
The first phase difference measuring unit 130 measures a first phase difference θHM between the hip joint angular velocity φH′ that changes with the angular velocity ωH and the first oscillator x in which the intrinsic angular velocity ωM included in the Van der Pol equation is reflected.
The walking condition determining unit 140 is equipped with a memory (not shown) storing a correspondence relationship between “a walking condition” and “a tracing pattern drawn in an n-dimensional space (plane)” by an n number of walking oscillators including the hip joint angular velocity φH′. The walking condition determining unit 140 also determines a “walking condition” on the basis of the “correspondence relationship” stored in the memory and the “tracing pattern drawn in an n-dimensional space by the n number of walking oscillators” measured by the walking oscillator measuring unit 110.
The desired phase difference determining unit 150 determines a desired phase difference θd on the basis of a “walking condition” determined by the walking condition determining unit 140.
The determining unit 160 is equipped with a second phase difference determining unit 162, a correlation factor determining unit 164, a first angular velocity determining unit 166, and a second angular velocity determining unit 168.
The second phase difference measuring unit 162 determines a second phase difference θhm(=θh−θm) between a virtual walking oscillator θh and a virtual walking assist oscillator θm according to an “internal model” representing a phase relationship between the virtual walking oscillator (phase) θh and the virtual walking assist oscillator (phase) θm.
The correlation factor determining unit 164 determines a correlation factor ε between the virtual walking oscillator θh and the virtual walking assist oscillator θm such that the second phase difference θhm (=θh−θm) approaches to a first phase difference θHM.
The first angular velocity determining unit 166 determines an angular velocity ωh of the virtual walking oscillator θh on the basis of the correlation factor ε.
Based on the angular velocity ωh of the virtual walking oscillator θh, the second angular velocity determining unit 168 determines an angular velocity ωm of the virtual walking assist oscillator θm such that the second phase difference θhm approaches to the desired phase difference θd determined by the desired phase difference determining unit 150.
The second oscillator. generating unit 170 generates second oscillators (to be discussed later) y that tug with each other such that the hip joint angular velocity (the walking oscillator) φH′ measured by the hip joint angular velocity measuring unit 112 and an angular velocity ωm of the virtual walking assist oscillator θm, determined by the second angular velocity determining unit 166 are reflected as an intrinsic angular velocity (ωM). The second oscillator y is generated according to a nonlinear oscillator model expressed by the Van der Pol equation, as in the case of the first oscillator x.
The neural oscillator generating unit 180 generates an output (neural oscillator) z of each neural element according to a neural oscillator model (to be discussed later) representing behaviors of a plurality of neural elements (neurons) related to walking on the basis of the second oscillator y and the hip joint angle (walking oscillator) φH measured by the hip joint angle measuring unit 114.
The walking assist oscillator generating unit 190 has a torque generating unit 192 that generates a torque (walking assist oscillator) T imparted to a walker by the walking assist device 200 on the basis of the output z of each neural element.
The function of the control system 100 of the walking assist device 200 will be explained in conjunction with FIG. 2 and FIG. 3.
First, based on an output of an angle sensor 220 (refer to arrow {circle around (1)} in FIG. 2), angular velocities φHL′ and φHR′ about left and right hip joints of the walker are measured by the hip joint angular velocity measuring unit 112 (s1 in FIG. 2).
Next, first oscillators XL and XR are determined (s2 in FIG. 2) by the first oscillator generating unit 120 on the basis of the hip joint angular velocities φHL′ and φHR′ (refer to arrow {circle around (1)}-1 in FIG. 2) according to the Van der Pol equation represented by expressions (1a) and (1b) given below.




  
    
      
        
          x
          L
          ′′
        
        =
        
          
            
              ξ
              ⁡
              
                (
                
                  1
                  -
                  
                    x
                    L
                    2
                  
                
                )
              
            
            ⁢
            
              
            
            ⁢
            
              x
              L
              ′
            
          
          -
          
            
              ω
              ML
              2
            
            ⁢
            
              
            
            ⁢
            
              x
              L
            
          
          +
          
            g
            ⁡
            
              (
              
                
                  x
                  L
                
                -
                
                  x
                  R
                
              
              )
            
          
          +
          
            
              c
              1
            
            ⁢
            
              
            
            ⁢
            
              φ
              HL
              ′
            
          
        
      
    
    
      
        (
        
          1
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          x
          R
          ′′
        
        =
        
          
            
              ξ
              ⁡
              
                (
                
                  1
                  -
                  
                    x
                    R
                    2
                  
                
                )
              
            
            ⁢
            
              
            
            ⁢
            
              x
              R
              ′
            
          
          -
          
            
              ω
              MR
              2
            
            ⁢
            
              
            
            ⁢
            
              x
              R
            
          
          +
          
            g
            ⁡
            
              (
              
                
                  x
                  R
                
                -
                
                  x
                  L
                
              
              )
            
          
          +
          
            
              c
              1
            
            ⁢
            
              
            
            ⁢
            
              φ
              HR
              ′
            
          
        
      
    
    
      
        (
        
          1
          ⁢
          b
        
        )
      
    
  




where ξ denotes a coefficient (>0) determined such that the first oscillator x and its time-dependent one differentiation x′ draw a stable limit period on an x-x′ plane. ωM denotes an intrinsic angular velocity. g denotes a correlation factor of right/left leg. c1 denotes a feedback coefficient. The intrinsic angular velocity ωM may be arbitrarily set within a range in which it does not deviate significantly from an actual walking assist rhythm by the walking assist device 200.
The first oscillators xL and xR are determined by the Runge-Kutta method. The first oscillators xL and xR represent virtual walking assist rhythms related to the left and right legs, respectively. Because of the “mutual tugging”, which is one of the properties of the Van der Pol equation, the first oscillators xL and xR have properties in which they vibrate with an autonomous rhythm in which the “intrinsic angular velocities” ωML and ωMR are reflected, while harmonizing with the hip joint angular velocities (walking oscillators) φHL′ and ωHR′ (refer to FIG. 3). The hip joint angular velocities (walking oscillators) φHL′ and φHR′ vibrate with a rhythm that is substantially equivalent to an actual walking rhythm (refer to FIG. 3.).
Subsequently, the phase difference measuring unit 130 determines (s3 in FIG. 2) first phase differences θHML (=θHL−θML), θHMR (=θHR−θMR) between the hip joint angular velocities (walking oscillators) φHL′ and φHR′ (phases θHL, θHR/refer to arrow {circle around (1)}-2 in FIG. 2) and the first oscillators xL and xR (phases θML, θMR/refer to arrow {circle around (2)} in FIG. 2). Specifically, as shown in FIG. 3, a phase difference θHM(=θH−θM) between a phase θH of the hip joint angular velocity φH′ and a phase θM of the first oscillator x is determined on the basis of a time difference between a point at which φH′=0 and φH′′>0 ( . . . tid, tid+1, . . . ) and a point at which x=0 and x′>0.
Furthermore, the walking condition determining unit 140 determines a walking condition (s4) on the basis of an n number (n =1, 2, . . . ) of walking oscillators (refer to arrow {circle around (1)}-3 in FIG. 2) that include hip joint angular velocities φH′ measured by the walking condition determining unit 110.
Specifically, first, the walking condition determining unit 140 reads, from a memory (not shown), a correspondence relationship between “walking condition” and “tracing patterns drawn in an n-dimensional space by an n number of walking oscillators.” Then, the walking condition determining unit 140 determines a “walking condition” on the basis of the correspondence relationship and the tracing patterns drawn in the n-dimensional space by the n number of walking oscillators measured by the walking oscillator measuring unit 110.
As the walking oscillators for determining walking conditions, various parameters, such as the hip joint angle φH, or the angles or angular velocities of knee joins, foot joints, shoulder joints or elbow joints, angular accelerations, the position of a part of a leg of a walker, or a landing noise, breath sound, or intended utterance of the walker, which vary with a rhythm associated with a walking rhythm may be measured by the walking oscillator measuring unit 110. Walking conditions of the walker determined by the walking condition determining unit 140 include “flat walk condition,” “ascending walk condition,” and “descending walk condition.”
The desired phase difference determining unit 150 determines a desired phase difference θd (s5 in FIG. 2) on the basis of a “walking condition” (arrow {circle around (4)} in FIG. 2) determined by the walking condition determining unit 140. More specifically, the desired phase difference determining unit 150 basically determines the desired phase difference θd to be 0 for the flat walk condition, for the ascending walk condition (e.g., −0.5 [rad] or less), and + for the descending walk condition (e.g., +0.3 [rad] or more).
If the first phase difference θHM measured by the phase difference measuring unit 130 remains constant for the past three walk periods, then a second phase difference θhm (=θh−θm) between a virtual walking oscillator θh and a virtual walking assist oscillator θm is determined according to “an internal model” represented by the following expressions (2a) and (2b) by the second phase difference determining unit 162.

θh′=ωh+ε·sin(θm−θh)  (2a)

θm′=ωm+ε·sin(θh−θm)  (2b)

where ε denotes correlation factors of the virtual walking oscillator θh and the walking assist oscillator θm in the internal model. ωh denotes an intrinsic angular velocity of the virtual walking oscillator θh, and ωm denotes an intrinsic angular velocity of the virtual walking assist oscillator θm.
Specifically, the second phase difference determining unit 162 determines the second phase difference θhm (=θh−θm) between the virtual walking oscillator θh and the virtual walking assist oscillator θm according to the following expression (3).

θhm=sin−1[(ωh−ωm)/2ε]  (3)

Subsequently, the correlation factor determining unit 164 determines the correlation factor ε (s6.1 in FIG. 2) so that a difference θHM−θhm between the first phase difference θHM (refer to arrow {circle around (3)} in FIG. 2) and the second phase difference θhm is minimized.
Specifically, the correlation factor ε is determined one after another at discrete time ( . . . , tid−1, tid, tid+1, . . . /refer to FIG. 3) at which the hip joint angular velocity (walking oscillator) φH′ becomes zero, according to the following expressions (4a) and (4b).




  
    
      
        
          ɛ
          ⁡
          
            (
            
              
                t
                id
              
              +
              1
            
            )
          
        
        =
        
          
            ɛ
            ⁡
            
              (
              
                t
                id
              
              )
            
          
          -
          
            η
            ⁢
            
              
                {
                
                  
                    
                      V
                      1
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          id
                        
                        +
                        1
                      
                      )
                    
                  
                  -
                  
                    
                      V
                      1
                    
                    ⁡
                    
                      (
                      
                        t
                        id
                      
                      )
                    
                  
                
                }
              
              /
              
                {
                
                  
                    ɛ
                    ⁡
                    
                      (
                      
                        t
                        id
                      
                      )
                    
                  
                  -
                  
                    ɛ
                    ⁡
                    
                      (
                      
                        
                          t
                          id
                        
                        -
                        1
                      
                      )
                    
                  
                
                }
              
            
          
        
      
    
    
      
        (
        
          4
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            V
            1
          
          ⁡
          
            (
            
              
                t
                id
              
              +
              1
            
            )
          
        
        ≡
        
          
            (
            
              1
              /
              2
            
            )
          
          ⁢
          
            
          
          ⁢
          
            
              {
              
                
                  
                    θ
                    HM
                  
                  ⁡
                  
                    (
                    
                      
                        t
                        id
                      
                      +
                      1
                    
                    )
                  
                
                -
                
                  
                    θ
                    hm
                  
                  ⁡
                  
                    (
                    
                      t
                      id
                    
                    )
                  
                
              
              }
            
            2
          
        
      
    
    
      
        (
        
          4
          ⁢
          b
        
        )
      
    
  




where η denotes a coefficient representing stability of a potential V1 for bringing the first phase difference θHM and the second phase difference θhm close to each other.
Then, based on the correlation factor ε (refer to arrow {circle around (4)}-1 in FIG. 2), the first angular velocity determining unit 166 determines, according to the following expressions (5a) and (5b), an intrinsic angular velocity ωh of the virtual walking oscillator θh such that the difference of the first and second phase differences θHM−θhm is minimized under a condition wherein the intrinsic angular velocity ωm of the virtual walking assist oscillator θm remains constant without changing with time (s6.2 in FIG. 2).




  
    
      
        
          
            ω
            hL
          
          ⁡
          
            (
            
              t
              idL
            
            )
          
        
        =
        
          
            -
            α
          
          ×
          
            ∫
            
              
                d
                t
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  
                    
                      [
                      
                        
                          4
                          ⁢
                          
                            
                          
                          ⁢
                          ɛ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            
                              (
                              
                                t
                                idL
                              
                              )
                            
                            2
                          
                        
                        -
                        
                          
                            {
                            
                              
                                
                                  ω
                                  hL
                                
                                ⁡
                                
                                  (
                                  t
                                  )
                                
                              
                              -
                              
                                
                                  ω
                                  mL
                                
                                ⁡
                                
                                  (
                                  
                                    t
                                    idL
                                  
                                  )
                                
                              
                            
                            }
                          
                          2
                        
                      
                      ]
                    
                    
                      1
                      /
                      2
                    
                  
                  ×
                  
                    sin
                    ⁡
                    
                      [
                      
                        
                          
                            sin
                            
                              -
                              1
                            
                          
                          ⁢
                          
                            {
                            
                              
                                
                                  (
                                  
                                    
                                      
                                        ω
                                        hL
                                      
                                      ⁡
                                      
                                        (
                                        t
                                        )
                                      
                                    
                                    -
                                    
                                      
                                        ω
                                        mL
                                      
                                      ⁡
                                      
                                        (
                                        
                                          
                                            t
                                            idL
                                          
                                          -
                                          1
                                        
                                        )
                                      
                                    
                                  
                                  )
                                
                                /
                                2
                              
                              ⁢
                              
                                ɛ
                                ⁡
                                
                                  (
                                  
                                    t
                                    idL
                                  
                                  )
                                
                              
                            
                            }
                          
                        
                        -
                        
                          
                            θ
                            HM
                          
                          ⁡
                          
                            (
                            
                              t
                              idL
                            
                            )
                          
                        
                      
                      ]
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          5
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            ω
            hR
          
          ⁡
          
            (
            
              t
              idR
            
            )
          
        
        =
        
          
            -
            α
          
          ×
          
            ∫
            
              
                d
                t
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  
                    
                      [
                      
                        
                          4
                          ⁢
                          
                            
                          
                          ⁢
                          ɛ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            
                              (
                              
                                t
                                idR
                              
                              )
                            
                            2
                          
                        
                        -
                        
                          
                            {
                            
                              
                                
                                  ω
                                  hR
                                
                                ⁡
                                
                                  (
                                  t
                                  )
                                
                              
                              -
                              
                                
                                  ω
                                  mR
                                
                                ⁡
                                
                                  (
                                  
                                    t
                                    idR
                                  
                                  )
                                
                              
                            
                            }
                          
                          2
                        
                      
                      ]
                    
                    
                      1
                      /
                      2
                    
                  
                  ×
                  
                    sin
                    ⁡
                    
                      [
                      
                        
                          
                            sin
                            
                              -
                              1
                            
                          
                          ⁢
                          
                            {
                            
                              
                                
                                  (
                                  
                                    
                                      
                                        ω
                                        hR
                                      
                                      ⁡
                                      
                                        (
                                        t
                                        )
                                      
                                    
                                    -
                                    
                                      
                                        ω
                                        mR
                                      
                                      ⁡
                                      
                                        (
                                        
                                          
                                            t
                                            idR
                                          
                                          -
                                          1
                                        
                                        )
                                      
                                    
                                  
                                  )
                                
                                /
                                2
                              
                              ⁢
                              
                                ɛ
                                ⁡
                                
                                  (
                                  
                                    t
                                    idR
                                  
                                  )
                                
                              
                            
                            }
                          
                        
                        -
                        
                          
                            θ
                            HM
                          
                          ⁡
                          
                            (
                            
                              t
                              idR
                            
                            )
                          
                        
                      
                      ]
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          5
          ⁢
          b
        
        )
      
    
  




where α denotes a coefficient representing stability of a system.
Subsequently, based on the intrinsic angular velocity ωh of the virtual walking oscillator θh (refer to arrow {circle around (4)}-2 in FIG. 2), the intrinsic angular velocity ωm of the virtual walking assist oscillator θm is determined by the second angular velocity determining unit 168 according to the following expressions (6a) and (6b) such that the second phase difference θHM approaches to the desired phase difference θd (refer to arrow {circle around (5)} in FIG. 2) determined by the desired phase difference determining unit 150 (s6.3 in FIG. 2).




  
    
      
        
          
            ω
            mL
          
          ⁡
          
            (
            
              t
              idL
            
            )
          
        
        =
        
          β
          ⁢
          
            
          
          ⁢
          
            ∫
            
              
                d
                t
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  
                    [
                    
                      
                        4
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            ɛ
                            ⁡
                            
                              (
                              
                                t
                                idL
                              
                              )
                            
                          
                          2
                        
                      
                      -
                      
                        
                          {
                          
                            
                              
                                ω
                                hL
                              
                              ⁡
                              
                                (
                                
                                  t
                                  idL
                                
                                )
                              
                            
                            -
                            
                              
                                ω
                                mL
                              
                              ⁡
                              
                                (
                                t
                                )
                              
                            
                          
                          }
                        
                        2
                      
                    
                    ]
                  
                  ×
                  
                    sin
                    ⁡
                    
                      [
                      
                        
                          
                            sin
                            
                              -
                              1
                            
                          
                          ⁢
                          
                            {
                            
                              
                                
                                  (
                                  
                                    
                                      
                                        ω
                                        hL
                                      
                                      ⁡
                                      
                                        (
                                        
                                          t
                                          idL
                                        
                                        )
                                      
                                    
                                    -
                                    
                                      
                                        ω
                                        mL
                                      
                                      ⁡
                                      
                                        (
                                        t
                                        )
                                      
                                    
                                  
                                  )
                                
                                /
                                2
                              
                              ⁢
                              
                                ɛ
                                ⁡
                                
                                  (
                                  
                                    t
                                    idL
                                  
                                  )
                                
                              
                            
                            }
                          
                        
                        -
                        
                          θ
                          d
                        
                      
                      ]
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          6
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            ω
            mR
          
          ⁡
          
            (
            
              t
              idR
            
            )
          
        
        =
        
          β
          ⁢
          
            
          
          ⁢
          
            ∫
            
              
                d
                t
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  
                    [
                    
                      
                        4
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            ɛ
                            ⁡
                            
                              (
                              
                                t
                                idR
                              
                              )
                            
                          
                          2
                        
                      
                      -
                      
                        
                          {
                          
                            
                              
                                ω
                                hR
                              
                              ⁡
                              
                                (
                                
                                  t
                                  idR
                                
                                )
                              
                            
                            -
                            
                              
                                ω
                                mR
                              
                              ⁡
                              
                                (
                                t
                                )
                              
                            
                          
                          }
                        
                        2
                      
                    
                    ]
                  
                  ×
                  
                    sin
                    ⁡
                    
                      [
                      
                        
                          
                            sin
                            
                              -
                              1
                            
                          
                          ⁢
                          
                            {
                            
                              
                                
                                  (
                                  
                                    
                                      
                                        ω
                                        hR
                                      
                                      ⁡
                                      
                                        (
                                        
                                          t
                                          idR
                                        
                                        )
                                      
                                    
                                    -
                                    
                                      
                                        ω
                                        mR
                                      
                                      ⁡
                                      
                                        (
                                        t
                                        )
                                      
                                    
                                  
                                  )
                                
                                /
                                2
                              
                              ⁢
                              
                                ɛ
                                ⁡
                                
                                  (
                                  
                                    t
                                    idR
                                  
                                  )
                                
                              
                            
                            }
                          
                        
                        -
                        
                          θ
                          d
                        
                      
                      ]
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          6
          ⁢
          b
        
        )
      
    
  




where β denotes a coefficient representing stability of a system.
Next, second oscillators yL and yR are determined (s7 in FIG. 2) by the second oscillator generating unit 170 on the basis of the hip joint angular velocities φHL′ and φHR′ (refer to arrow {circle around (1)}-4 in FIG. 2) according to the Van der Pol equation represented by expressions (7a) and (7b) given below that include intrinsic angular velocities ωmL and ωmR determined by the second angular velocity determining unit 168 (refer to arrow {circle around (4)}-3 in FIG. 2).




  
    
      
        
          y
          L
          ′′
        
        =
        
          
            
              ξ
              ⁡
              
                (
                
                  1
                  -
                  
                    y
                    L
                    2
                  
                
                )
              
            
            ⁢
            
              
            
            ⁢
            
              y
              L
              ′
            
          
          -
          
            
              ω
              mL
            
            ⁢
            
              
            
            ⁢
            
              y
              L
            
          
          +
          
            g
            ⁡
            
              (
              
                
                  y
                  L
                
                -
                
                  y
                  R
                
              
              )
            
          
          +
          
            
              c
              1
            
            ⁢
            
              
            
            ⁢
            
              φ
              HL
              ′
            
          
        
      
    
    
      
        (
        
          7
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          y
          R
          ′′
        
        =
        
          
            
              ξ
              ⁡
              
                (
                
                  1
                  -
                  
                    y
                    R
                    2
                  
                
                )
              
            
            ⁢
            
              
            
            ⁢
            
              y
              R
              ′
            
          
          -
          
            
              ω
              mR
            
            ⁢
            
              
            
            ⁢
            
              y
              R
            
          
          +
          
            g
            ⁡
            
              (
              
                
                  y
                  R
                
                -
                
                  y
                  L
                
              
              )
            
          
          +
          
            
              c
              1
            
            ⁢
            
              
            
            ⁢
            
              φ
              HR
              ′
            
          
        
      
    
    
      
        (
        
          7
          ⁢
          b
        
        )
      
    
  




The second oscillators yL and yR are determined by the Runge-Kutta method, as in the case of the first oscillators xL and xR. The second oscillators yL and yR represent virtual walking assist rhythms related to the left and right legs, respectively. Because of the “mutual tugging”, which is one of the properties of the Van der Pol equation, the second oscillators yL and yR have properties in which they vibrate with an autonomous rhythm in which the “intrinsic angular velocities” ωmL and ωmR are reflected, while harmonizing with the rhythms of the hip joint angular velocities (walking oscillators) φHL′ and φHR′ (refer to FIG. 3).
Subsequently, a neural oscillator z is generated by the neural oscillator generating unit 180 (s8 in FIG. 2) according to a “neural oscillator model” on the basis of the second oscillators yL and yR (refer to arrow {circle around (7)} in FIG. 2) and the left and right hip joint angles φL and φR (refer to arrow {circle around (1)} ′ in FIG. 2) measured by the hip joint angle measuring unit 114. The left and right hip joint angles φL and φR of the walker are measured by the hip joint angle measuring unit 114 in parallel to the measurement of the hip joint angular velocity φH′ by the hip joint angular velocity measuring unit 112 (s1′ in FIG. 2).
Specifically, outputs zL+ and zL− of neural elements L+ and L− that govern the motions of the left thigh in a bending direction (to the front) and a stretching direction (to the rear) and outputs zR+ and zR− of neural elements R+ and R− that govern the motions of the right thigh in the bending direction and the stretching direction are determined according to the following expressions (8a) through (8f).




  
    
      
        
          
            τ
            
              L
              +
            
          
          ⁢
          
            u
            
              L
              +
            
            ′
          
        
        =
        
          
            -
            
              u
              
                L
                +
              
            
          
          +
          
            
              w
              
                L
                +
                
                  /
                  L
                
                -
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                L
                -
              
            
          
          +
          
            
              w
              
                L
                +
                
                  /
                  R
                
                +
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                R
                +
              
            
          
          -
          
            λ
            ⁢
            
              
            
            ⁢
            
              v
              
                L
                +
              
            
          
          +
          
            
              k
              +
            
            ⁢
            
              
            
            ⁢
            
              φ
              R
            
          
          -
          
            c
            ⁢
            
              
            
            ⁢
            
              y
              R
            
          
        
      
    
    
      
        (
        
          8
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            τ
            
              L
              -
            
          
          ⁢
          
            u
            
              L
              -
            
            ′
          
        
        =
        
          
            -
            
              u
              
                L
                -
              
            
          
          +
          
            
              w
              
                L
                -
                
                  /
                  L
                
                +
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                L
                +
              
            
          
          +
          
            
              w
              
                L
                -
                
                  /
                  R
                
                -
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                R
                -
              
            
          
          -
          
            λ
            ⁢
            
              
            
            ⁢
            
              v
              
                L
                -
              
            
          
          -
          
            
              k
              -
            
            ⁢
            
              
            
            ⁢
            
              φ
              R
            
          
          +
          
            c
            ⁢
            
              
            
            ⁢
            
              y
              L
            
          
        
      
    
    
      
        (
        
          8
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            τ
            
              R
              +
            
          
          ⁢
          
            u
            
              R
              +
            
            ′
          
        
        =
        
          
            -
            
              u
              
                R
                +
              
            
          
          +
          
            
              w
              
                R
                +
                
                  /
                  L
                
                +
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                L
                +
              
            
          
          +
          
            
              w
              
                R
                +
                
                  /
                  R
                
                -
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                R
                -
              
            
          
          -
          
            λ
            ⁢
            
              
            
            ⁢
            
              v
              
                R
                +
              
            
          
          +
          
            
              k
              +
            
            ⁢
            
              
            
            ⁢
            
              φ
              L
            
          
          +
          
            c
            ⁢
            
              
            
            ⁢
            
              y
              L
            
          
        
      
    
    
      
        (
        
          8
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            τ
            
              R
              -
            
          
          ⁢
          
            u
            
              R
              -
            
            ′
          
        
        =
        
          
            -
            
              u
              
                R
                -
              
            
          
          +
          
            
              w
              
                R
                -
                
                  /
                  L
                
                -
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                L
                -
              
            
          
          +
          
            
              w
              
                R
                -
                
                  /
                  R
                
                +
              
            
            ⁢
            
              
            
            ⁢
            
              z
              
                R
                +
              
            
          
          -
          
            λ
            ⁢
            
              
            
            ⁢
            
              v
              
                R
                -
              
            
          
          +
          
            
              k
              -
            
            ⁢
            
              
            
            ⁢
            
              φ
              R
            
          
          -
          
            c
            ⁢
            
              
            
            ⁢
            
              y
              R
            
          
        
      
    
    
      
        (
        
          8
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            
              
                
                  σ
                  i
                
                ⁢
                
                  v
                  i
                  ′
                
              
              =
              
                
                  -
                  
                    v
                    i
                  
                
                +
                
                  z
                  i
                
              
            
          
          
            
              (
              
                
                  i
                  =
                  
                    L
                    +
                  
                
                ,
                
                  L
                  -
                
                ,
                
                  R
                  +
                
                ,
                
                  R
                  -
                
              
              )
            
          
        
      
    
    
      
        (
        
          8
          ⁢
          e
        
        )
      
    
  
  
    
      
        
          z
          i
        
        =
        
          max
          ⁡
          
            (
            
              0
              ,
              
                u
                i
              
            
            )
          
        
      
    
    
      
        (
        
          8
          ⁢
          f
        
        )
      
    
  




where ui denotes a condition variable associated with a change in membrane potential of a neural element i, vi denotes a self inhibitor in which adaptation effect of the neural element i is reflected, τi and σi are time constants of the condition variable ui and the self inhibitor vi, wi/j (<0) denotes an association constant of the neural elements i and j, λ denotes a habituation coefficient, k+ and k− denote feedback coefficients, and c denotes a bias coefficient.
Next, torques about hip joints (walking assist oscillators) TL and TR imparted to the walker by the right and left actuators 220 are generated (s9 in FIG. 2) by the torque generating unit 192 on the basis of an output yi of the neural element i (refer to arrow {circle around (8)} in FIG. 2) according to expressions (9a) and (9b) given below.

TL=p+zL+−p−zL−  (9a)

TR=−p+zR++p−zR−  (9b)

where p+ and p− denote activation coefficients.
Electric power based on the generated torques TL and TR is supplied from the battery 230 to the right and left actuators 210, so that the torques TL and TR are imparted to the walker through the actuators 210 (refer to arrow {circle around (9)} in FIG. 2).
Thereafter, the above processing (s1 through s9 in FIG. 2) is repeated, causing the walker to walk with the walking assistance by the walking assist device 200.
According to the control system 100, it is possible to generate a second oscillator y having a proper phase difference ( ̃desired phase difference θd) from the viewpoint of walking assistance with respect to the hip joint angular velocity (walking oscillator) φH′ in which a walking rhythm is reflected, and then the torque (walking assist oscillator) T is generated on the basis of the second oscillator y, thereby allowing a walking assist rhythm to be optimized.
To be more specific, first, (1) a first oscillator x that tugs and is tugged by the hip joint angular velocity (walking oscillator) φH′ in which the intrinsic angular velocity ωM is reflected is generated (s2 in FIG. 2). The first oscillator x vibrates with an autonomous rhythm in which the intrinsic angular velocity ωM is reflected, while harmonizing with the rhythm of the hip joint angular velocity (walking oscillator) φH′ due to the “mutual tugging” effect.
However, although the first oscillator x harmonizes with the rhythm ( ̃walking rhythm) of the hip joint angular velocity (the walking oscillator) φH′, the first phase difference θHM relative to the hip joint angular velocity φH′ is not necessarily close to the desired phase difference θd appropriate for assisting walk. Hence, if the torque (walking assist oscillator) T is determined on the basis of the first oscillator x, then a phase difference ( ̃a first phase difference θHM) inappropriate for assisting walk may be generated between a walking rhythm and a walking assisting rhythm, leading to inappropriate walking with the walking assistance. Accordingly, it is necessary to generate a new oscillator having a phase difference appropriate for walking assistance relative to the hip joint angular velocity (walking oscillator) φH′ on the basis of the first oscillator x.
Hence, (2) an angular velocity (a new intrinsic angular velocity) ωm of a virtual walking assist oscillator θm is determined according to an “internal model” (s6.1 through s6.3 in FIG. 2). Furthermore, a second oscillator y is generated that tugs and is tugged by the hip joint angular velocity (the walking oscillator) φH′ such that the intrinsic angular velocity ωm is reflected, and has a second phase difference θhm which is closer to the desired phase difference θd than the first phase difference θHM is, with respect to the walking oscillator (s7 in FIG. 2).
Generating the second oscillator y is equivalent to correcting the intrinsic angular velocity ωM to a new intrinsic angular velocity ωM so as to bring the first phase difference θHM closer to the desired phase difference θd (refer to the above expressions (1a) and (1b)), and then regenerating the first oscillator x such that the corrected intrinsic angular velocity ωm is reflected (refer to the above expressions (7a) and (7b)).
The rhythm of the second oscillator y thus generated has the desired phase difference θd or a phase difference close thereto ( ̃a second phase difference θhm) with respect to a walking rhythm ( ̃rhythm of the walking oscillator φH). Thus, if the torque (the walking assist oscillator) T is determined on the basis of the second oscillator y, then the phase difference between the walking rhythm and the walking assist rhythm approaches to the desired phase difference θd.
Thus, (3) an output (neural oscillator) zi of the neural element i is generated according to the “neural element model” on the basis of the second oscillator y and the hip joint angle (the walking oscillator) φH (s8 in FIG. 2), and then the torque (walking assist oscillator) T is generated on the basis of each output zi (s9 in FIG. 2). This generates the torque (the walking assist oscillator) T so as to allow the walking rhythm ( ̃the rhythm of the walking oscillator φH′) and the desired phase difference θd to be achieved. Furthermore, the walk of the walker can be assisted at a rhythm based on a walking condition reflected in the walking rhythm.
Accordingly, the present invention makes it possible to achieve a walking assist rhythm that (1) follows changes in the walking rhythm ( ̃the rhythm of the walking oscillator φH′) of a walker, while (2) maintaining autonomy that has the walking rhythm and a desired phase difference θd, and (3) achieving smooth walking assistance according to a walking condition of the walker.
As described above, a walking assist rhythm harmonizes with a change in a walking rhythm when the walking rhythm changes, and the walking rhythm harmonizes with the walking assist rhythm, thus allowing the walker (human body) and the walking assist device (machine) 200 to harmonize with each other (mutual concession). This arrangement makes it possible to achieve appropriate walking assistance while allowing a walker to reasonably feel that his/her walk is being assisted by the walking assist device 200. Moreover, even if a walking rhythm suddenly changes, the walking assist rhythm does not fully follow the sudden change, making it possible to avoid furthering walking assistance and walking that cause burden to the mind and body of the walker.
Moreover, the phase difference between a walking rhythm and a walking assist rhythm is brought close to the desired phase difference θd that changes according the walking condition of a walker (refer to arrow {circle around (4)} in FIG. 2) (refer to s5, etc. in FIG. 2). This makes it possible to accomplish walking with proper walking assistance responsive to changes in the walking condition of a walker.
Furthermore, when the desired phase difference θd is determined to be + (refer to s5 in FIG. 2), the walking rhythm and the walking assist rhythm will have a phase difference of the determined value (>0) or a phase difference close thereto, enabling the walker to walk by leading the walking assist device 200 (walker-led walking). Conversely, when the desired phase difference θd is determined to be—(refer to s5 in FIG. 2), the walking rhythm and the walking assist rhythm will have a phase difference of the determined value (<0) or a phase difference close thereto, enabling the walker to walk by being led by the walking assist device 200 (walk led by the walking assist device).
Thus, the walker can walk by leading the walking assist device 200 in “a descending walk condition” wherein burden is relatively small, while by being led by the walking assist device 200 in “an ascending walk condition” wherein burden is relatively large. The desired phase difference can be changed according to walker's intention or other factors rather than being fixed to 0, + or − according to a walking condition.
Referring now to FIG. 4 through FIG. 9, the results of experiments of walk using the walking assist device 200 controlled by the control system 100 will be explained.
FIG. 4, FIG. 5, and FIG. 6 show how the phase difference between the hip joint angular velocity (walking oscillator) φH′ and the torque (walking assist oscillator) T changes when the desired phase difference θd is set to 0.0[rad], −0.2[rad], and −0.5[rad], respectively. The solid lines indicate phase differences on the right leg, while the dashed lines indicate phase differences on the left leg.
As shown in FIG. 4 through FIG. 6, the phase differences exhibit relatively large changes immediately after time t0 when walking assistance begins. Then, the phase differences gradually converge to the desired phase difference θd respectively set. As is obvious from the results, according to the control system 100, the phase difference between the walking rhythm ( ̃vibration rhythm of the hip joint angular velocity φH′) and the walking assist rhythm ( ̃vibration rhythm of the torque T) can be brought close to the desired phase difference θd determined from the viewpoint of walking assistance with a proper rhythm.
FIG. 7 and FIG. 8 show changes in walking period TH (solid lines) and walking assist period TM (dashed lines). As shown in FIG. 7, if a walker suddenly increases the walking speed at time t1, resulting in a shorter walking period TH, then the walking assist period TM shortens accordingly. If the walker suddenly decreases the walking speed at time t2, resulting in a longer walking period TH, then the walking assist period TM is prolonged accordingly. It is understood, however, that the walking assist period TM does not completely follow the walking period TH. Instead, the walking assist period TM is slightly changed, and the walking period TH changes to follow TM after the change.
As is obvious from the above, the control system 100 allows a walking assist rhythm to harmonize with a change in a walking rhythm when the walking rhythm changes, and the walking rhythm harmonizes with the walking assist rhythm, thus allowing the walker (human body) and a walking assist device (machine) 200 to harmonize with each other (mutual concession). Thus, it is possible to achieve appropriate walking assistance while allowing a walker to reasonably feel that his/her walk is being assisted by the walking assist device 200. Moreover, even if a walking rhythm suddenly changes, the walking assist rhythm does not fully follow the sudden change. This arrangement is expected to allow avoidance of furthering of walking assistance and walking that cause burden to the mind and body of the walker.
The same applies to a case where the walking period TH (solid lines) slowly changes, as shown in FIG. 8, rather than changing suddenly.
FIG. 9 shows how a walk symmetry parameter A indicating symmetric properties of motions of the right and left legs changes when one leg is loaded with a weight ( ̃11 kg), the control of the walking assist device 200 by the control system 100 is intermittently turned ON/OFF, and a walker makes a U turn at time t3. The walk symmetry parameter A is represented by a difference between the landing period of one leg and the landing period of the other leg. The walk symmetry parameter A approaches to zero when both legs are moving with equivalent landing periods, whereas it shifts away from zero if a difference is developed between the landing periods of the two legs, as in a case where one leg is dragged or the like.
As is obvious from FIG. 9, while the control of the walking assist device 200 by the control system 100 is ON, an absolute value of the walk symmetry parameter A is smaller, as a whole, than while the control is OFF. This indicates that the landing rhythms of the two legs are harmonized by the control system 100, whereas the one leg loaded with the weight would be dragged when the walker walks, resulting in an increased difference between the landing periods of the two legs.
The landing rhythms of the two legs are harmonized, as described above, because terms including correlation factors g of the right and left legs are included in the Van der Pol equations of the above expressions (1a), (1b), (7a), and (7b). More specifically, the first oscillators xL and xR associated with virtual left and right motions tug with each other, and similarly, the second oscillators yL and yR associated with virtual left and right motions tug with each other. Hence, even if only one leg is loaded (due to an injury or the like), adjusting the correlation factors g makes it possible to further reduce the absolute value of the walk symmetry parameter A and to achieve walking similar to normal walking with the same landing rhythms for both right and left legs.
In the present embodiment, the torque T has been imparted about hip joints after the hip joint angle φH and the angular velocity φH′ are measured. As another embodiment, however, a torque may be imparted about knee joints after a knee joint angle and angular velocity are measured, or a torque may be imparted about a foot joint after a foot joint angle and angular velocity are measured.
In the present embodiment, the torques about a pair of joints have been imparted to the walker. However, as another embodiment, torques about a plurality of pairs of joints (the right and left hips+knee joints, the right and left hips+foot joints, the right and left knees+foot joints, and the right and left hips+knees+foot joints) may be imparted to a walker. For instance, if torques are to be imparted about all joints (hip+knee+foot joints) of the right and left legs, then the torque to be imparted about each joint may be determined according to a neural oscillator model representing the behaviors of neural elements of a plurality of pairs related to stretching and bending for each joint, which interact with each other, as shown in FIG. 10.
As the number of walking oscillators to be measured increases, the number of correlation terms in the nonlinear differential equations (refer to the above expressions (1a), (1b), (7a), and (7b)) related to the generation of the first and second oscillators x and y of the Van der Pol equations or the like, and in the nonlinear differential equations (refer to the above expressions (8a) through (8f)) related to the generation of neural oscillators increases. Nevertheless, more accurate walking assistance based on the motions of diverse parts of the body of a walker can be achieved by adjusting the correlation factors.
In the present embodiment, the first oscillators x and the second oscillators y have been generated according to the Van der Pol equations represented by the equations (1a), (1b), (7a), and (7b). As another embodiment, however, the first oscillators x and the second oscillators y may be generated according Van der Pol equations of a form different from the above expressions, or according to any forms of nonlinear differential equations that permit mutual tugging among a plurality of oscillators.
In the present embodiment, the hip joint angular velocity φH′ has been measured as the “walking oscillator” related to the generation of the first and second oscillators x and y (refer to arrows {circle around (1)}-1 and {circle around (1)}-4 in FIG. 2). As another embodiment, however, various parameters, such as the hip joint angle φH, or the angles or angular velocities of knee joins, foot joints, shoulder joints or elbow joints, or a landing noise, breath sound, or intended utterance of the walker, which vary with a rhythm associated with a walking rhythm may be measured as the “walking oscillators” related to the generation of the first and second oscillators x and y.
In the present embodiment, the hip joint angle φH has been measured as the “walking oscillator” related to the generation of the torque (walking assist oscillator) T (refer to arrow {circle around (6)} in FIG. 2). As another embodiment, however, various parameters, such as the hip joint angular velocity φH′, or the angles or angular velocities of knee joins, foot joints, shoulder joints or elbow joints, or a landing noise, breath sound, or intended utterance of the walker, which vary with a rhythm associated with a walking rhythm may be measured as the “walking oscillators” related to the generation of the torque T.
In the present embodiment, the torques T about joints to be imparted to a walker through the actuators 210 have been generated as the “walking assist oscillators” (refer to s8 in FIG. 2). In another embodiment, the torques T as the “walking assist oscillators” may be replaced by, for example, periodic sounds that can be acoustically perceived by a walker through an audio device, such as a headphone (not shown), periodic light or a sign that can be perceived through a visual device (not shown), such as goggles, or periodic tapping (knocking) that can be perceived, via a device or the like, by a walker in terms of sensation at a part of his/her body, such as the back or shoulder.
In the present embodiment, walking conditions have been determined on the basis of the walking oscillators measured by the walking oscillator measuring unit 110. As another embodiment, physiological parameters representing a physiological condition of a walker (e.g., the potentials of muscles that move when the walker walks) may be measured by a physiological parameter measuring unit (not shown), and then a “walking condition” may be determined by the walking condition determining unit 140 on the basis of the pattern of the physiological parameters.
In the present embodiment, a “walking condition” has been determined as “the ascending walk condition,” “the flat walk condition,” or “the descending walk condition.” In another embodiment, as a “walking condition,” the speed of walking may be determined separately or in addition to the above. In this embodiment, the desired phase difference determining unit 150 may basically determine the desired phase difference θd to be − (e.g., −0.2 [rad] or less) if walking is fast, or + (e.g., +0.5 [rad] or more) if walking is slow.
The present embodiment enables the walker to walk by leading the walking assist device 200 in “slow walking” wherein burden is relatively small, while by being led by the walking assist device 200 in “fast walking” wherein burden is relatively large.
If “the descending walk condition” for which the desired phase difference θd is basically determined to be + is combined with “the fast walking” for which the desired phase difference θd is basically determined to be −, then a basic desired phase difference θd may be determined by fuzzy control that permit these two conditions to be matched.
Another embodiment of the control system 100 in accordance with the present invention may be provided with “a first correlation adjusting unit” for adjusting the correlation factors g (refer to the above expressions (1a), (1b), (7a), and (7b)) of the right and left legs (a plurality of virtual objects) related to the generation of the first and second oscillators x and y according to the “walking condition” determined by the walking condition determining unit 140.
The present embodiment makes it possible to properly adjust a walking assist rhythm on the basis of a qualitative correlation between the right and left legs, and also to achieve walking assistance in real time at an optimum rhythm for a walking condition, taking into account the fact that the qualitative correlation changes as the walking condition changes.
Still another embodiment of the control system 100 in accordance with the present invention may be provided with “a second correlation adjusting unit” for adjusting a correlation factor wi/j (refer to the above expressions (8a) through (8d)) among a plurality of virtual neural elements i and j (refer to FIG. 2) related to the generation of walking assist oscillators, such as the torques T, according to the “walking condition” determined by the walking condition determining unit 140.
The present embodiment makes it possible to properly adjust a walking assist rhythm by adjusting the correlation intensity of a plurality of virtual neural elements, and also to achieve walking assistance in real time at an optimum rhythm for a walking condition, taking into account the fact that the correlation intensity changes as the walking condition changes.
Still another embodiment of the control system 100 in accordance with the present invention may be provided with a “physiological condition determining unit” for determining the “physiological condition” of a walker on the basis of a walking assist oscillator measured by the walking assist oscillator measuring unit 110, or the amounts of other conditions (e.g., the heart rate, the respiration rate, the concentration of lactic acid, oxygen, etc. of blood, the amount of perspiration, and the number of blinks of the walker). Then, the desired phase difference θd may be determined by the desired phase difference determining unit 150 on the basis of the “physiological condition” determined by the physiological condition determining unit.
According to the present embodiment, the phase difference between a walking rhythm and a walking assist rhythm can be brought close to the desired phase difference θd that varies according to the physiological condition of the walker, thus making possible to accomplish walking with proper walking assistance responsive to changes in physiological condition.
Alternatively, “the physiological condition determining unit” may determine the “fatigue” of a walker, and then “the desired phase difference determining unit” may determine the desired phase difference θd so that it is basically determined to be 0 for a flat physiological condition, − (e.g., −0.5[rad] or less) for an ascending physiological condition, and + (e.g., +0.3[rad] or more) for a descending physiological condition.
The present embodiment enables a walker to walk by leading the walking assist device 200 when the level of his/her fatigue is low, and to walk by being led by the walking assist device 200 when the level of his/her fatigue is high. The desired phase difference can be changed according to walker's intention or other factors rather than being fixed to + or − according to the physiological condition.
Yet another embodiment of the control system 100 in accordance with the present invention may be provided with a “physiological condition determining unit” for determining the physiological condition of the walker on the basis of a “walking assist oscillator” measured by the walking assist oscillator measuring unit 110, and “a first correlation adjusting unit” for adjusting correlation factors g (refer to the above expressions (1a), (1b), (7a), and (7b)) of the right and left legs (a plurality of virtual objects) related to the generation of the first and second oscillators x and y according to the “physiological condition” determined by the physiological condition determining unit.
The present embodiment makes it possible to properly adjust a walking assist rhythm on the basis of a qualitative correlation between the right and left legs, and also to achieve walking assistance in real time at an optimum rhythm for a physiological condition, taking into account the circumstance that the qualitative correlation changes as the physiological condition changes.
Still another embodiment of the control system 100 in accordance with the present invention may be provided with “a physiological condition determining unit” for determining a “physiological condition” of a walker on the basis of “walking assist oscillators” measured by the walking assist oscillator measuring unit 110, and “a second correlation adjusting unit” for adjusting a correlation factor wi/j (refer to the above expressions (8a) through (8d)) among a plurality of virtual neural elements i and j (refer to FIG. 2) related to the generation of walking assist oscillators according to the “physiological condition” determined by the physiological condition determining unit.
The present embodiment makes it possible to properly adjust a walking assist rhythm by adjusting the correlation intensity of a plurality of virtual neural elements, and also to achieve walking assistance in real time at an optimum rhythm based on physiological conditions, taking into account the circumstance that the correlation intensity changes as the physiological condition changes.
Another embodiment of the control system 100 in accordance with the present invention may adopt the following construction. First, the functions of the determining unit 160 and the second oscillator generating unit 170 are temporarily suspended. Then, the first oscillator x is directly supplied to the neural oscillator generating unit 180 (directly proceeding from arrow {circle around (2)} to arrow {circle around (7)} in FIG. 2), and the neural oscillator z and the torque T are determined (s8 and s9 in FIG. 2). Subsequently, a third phase difference determining unit (not shown) measures the phase difference between the hip joint angular velocity (walking oscillator) φH′ and the torque (walking assist oscillator) T. Further, the desired phase difference determining unit 150 determines the desired phase difference θd that is identical or close to the aforesaid phase difference. Then, the functions of the determining unit 160 and the second oscillator generating unit 170 are restored, and the second oscillator y is generated (refer to s6.1 through s6,3 and s7 in FIG. 2) on the basis of the latest desired phase difference θd (refer to arrow {circle around (5)} in FIG. 2). Lastly, the neural oscillator z and the torque T are generated (refer to s8 and s9 in FIG. 2) on the basis of the second oscillator y (refer to arrow {circle around (7)} in FIG. 2).
According to the present embodiment, temporarily suspending the functions of the determining unit 160 and the second oscillator generating unit 170 generally causes the phase difference between the hip joint angular velocity (walking oscillator) φH′ and the torque (walking assist oscillator) T to move away from the desired phase difference θd up to that moment. In addition, the phase difference converges to an optimum phase difference based on the harmony between a walker and the walking assist device 200. Thereafter, a new desired phase difference θd is determined that is identical or close to the aforesaid phase difference, and then the functions of the determining unit 160 and the second oscillator generating unit 170 are restored to control the torque (walking assist oscillator) T so as to bring the walk and the walking assistance close to the new desired phase difference θd.
As described above, the new desired phase difference θd is an optimum phase difference (or a phase difference approximate thereto) from the viewpoint of harmony between a walker and the walking assist device 200. This enables the walker to walk with a rhythm optimally matched with a walking assist rhythm by the walking assistance provided by the walking assist device 200.

## Claims

1. A control system comprising:
walking oscillator measuring means for measuring a walking oscillator of a walker whose walk is assisted by a walking assist device;
first oscillator generating means for generating a first oscillator that tugs and is tugged by a walking oscillator measured by the walking oscillator measuring means such that an intrinsic angular velocity is reflected;
determining means for determining a new intrinsic angular velocity on the basis of a difference between a first phase difference between the first oscillator and the walking oscillator and a desired phase difference;
second oscillator generating means for generating a second oscillator that tugs and is tugged by a walking oscillator measured by the walking oscillator measuring means such that the intrinsic angular velocity determined by the determining means is reflected, and has a second phase difference, which is closer to the desired phase difference than the first phase difference is, with respect to the walking oscillator; and
walking assist oscillator generating means for generating a walking assist oscillator of the walking assist device on the basis of the second oscillator and the walking oscillator measured by the walking oscillator measuring means.

2. The control system according to claim 1, comprising: walking condition determining means for determining the walking condition of the walker on the basis of a walking oscillator measured by the walking oscillator measuring means; and desired phase difference determining means for determining a desired phase difference on the basis of the walking condition of the walker determined by the walking condition determining means.

3. The control system according to claim 1, comprising:
physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker; walking condition determining means for determining a walking condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means, and desired phase difference determining means for determining a desired phase difference on the basis of a walking condition of the walker determined by the walking condition determining means.

4. The control system according to claim 2, wherein
the walking condition determining means determines whether a walking condition of a walker is a flat walk condition, an ascending walk condition, or a descending walk condition, and
the desired phase difference determining means basically determines a desired phase difference to be 0 for the flat walk condition, − for the ascending walk condition, and + for the descending walk condition.

5. The control system according to claim 2, wherein
the walking condition determining means determines the speed of walk as a walking condition, and
the desired phase difference determining means basically determines a desired phase difference to be − if the walk is fast, or + if the walk is slow.

6. The control system according to claim 1, comprising:
physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker; a physiological condition determining means for determining a physiological condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means; and
desired phase difference determining means for determining a desired phase difference on the basis of a physiological condition of the walker determined by the physiological condition determining means.

7. The control system according to claim 6, wherein the physiological condition determining means determines fatigue of a walker as a physiological condition of the walker, and basically determines a desired phase difference to be + if the fatigue of the walker is small, or − if the fatigue of the walker is large.

8. The control system according to claim 1, wherein
the first oscillator generating means generates a plurality of first oscillators corresponding to the motions of a plurality of objects such that correlation among a plurality of virtual objects is reflected, and
the second oscillator generating means generates a plurality of second oscillators corresponding to the motions of the plurality of the objects such that correlation among a plurality of virtual objects is reflected.

9. The control system according to claim 8, comprising:
walking condition determining means for determining a walking condition of a walker on the basis of a walking oscillator measured by a walking oscillator measuring means; and
first correlation adjusting means for adjusting a correlation among a plurality of virtual objects related to generation of first and second oscillators on the basis of a walking condition of the walker determined by the walking condition determining means.

10. The control system according to claim 8, comprising:
physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker; walking condition determining means for determining a walking condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means; and
first correlation adjusting means for adjusting a correlation among a plurality of virtual objects related to the generation of the first and second oscillators on the basis of the walking condition of the walker determined by the walking condition determining means.

11. The control system according to claim 8, comprising:
physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker, a physiological condition determining means for determining the physiological condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means; and
first correlation adjusting means for adjusting a correlation among a plurality of virtual objects related to the generation of the first and second oscillators on the basis of the physiological condition of the walker determined by the physiological condition determining means.

12. The control system according to claim 1, wherein determining means determines an intrinsic angular velocity for bringing a second phase difference between a second oscillator and a walking oscillator closer to a desired phase difference than a first phase difference between a first oscillator and a walking oscillator according to an internal model representing a phase relationship between a virtual walking oscillator and a virtual walking assist oscillator.

13. The control system according to claim 12, wherein the determining means comprises:
correlation factor determining means for determining correlation factors of a virtual walking oscillator and a virtual walking assist oscillator such that the second phase difference of the virtual walking oscillator and the virtual walking assist oscillator approaches to the first phase difference according to the internal model;
first angular velocity determining means for determining an angular velocity of the virtual walking oscillator on the basis of the correlation factors; and
second angular velocity determining means for determining an angular velocity of the virtual walking assist oscillator as an intrinsic angular velocity related to the generation of the second oscillator such that the second phase difference approaches to a desired phase difference on the basis of the angular velocity of the virtual walking oscillator.

14. The control system according to claim 1, wherein the walking assist oscillator generating means generates walking assist oscillators corresponding to behaviors of a plurality of neural elements such that a correlation among the plurality of virtual neural elements is reflected.

15. The control system according to claim 14, comprising:
walking condition determining means for determining a walking condition of a walker on the basis of a walking oscillator measured by a walking oscillator measuring means; and
second correlation adjusting means for adjusting a correlation among a plurality of virtual neural elements related to the generation of a walking assist oscillator on the basis of the walking condition of the walker determined by the walking condition determining means.

16. The control system according to claim 14, comprising:
physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker; walking condition determining means for determining the walking condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means; and
second correlation adjusting means for adjusting a correlation among a plurality of virtual neural elements related to the generation of the a walking assist oscillator on the basis of the walking condition of the walker determined by the walking condition determining means.

17. The control system according to claim 14, comprising:
physiological parameter measuring means for measuring physiological parameters indicating a physiological condition of a walker; physiological condition determining means for determining the physiological condition of the walker on the basis of the physiological parameters measured by the physiological parameter measuring means; and
second correlation adjusting means for adjusting a correlation among a plurality of virtual neural elements related to the generation of the walking assist oscillator on the basis of the walking condition of the walker determined by the walking condition determining means.

18. The control system according to claim 2, comprising: storing means for storing a correspondence relationship between a walking condition and a tracing pattern drawn by an n number of walking oscillators in an n-dimensional space, wherein the walking condition determining means determines a walking condition on the basis of the correspondence relationship stored by the storing means and the tracing pattern drawn in the n-dimensional space by the n number of walking oscillators measured by the walking oscillator measuring means.

19. The control system according to claim 1, wherein the first and second oscillator generating means generate the first and second oscillators according to a nonlinear oscillator model represented by the Van der Pol equation that includes a walking oscillator measured by the walking oscillator measuring means.

20. The control system according to claim 1, wherein the walking oscillator measuring means measures the first oscillator, the second oscillator, and at least one of various oscillators including a joint angle and angular velocity of a walker that periodically changes at a rhythm corresponding to a walking rhythm, as a walking oscillator for generating a walking assist oscillator.

21. The control system according to claim 1, wherein the walking assist oscillator generating means generates, as a walking assist oscillator, a torque about a leg joint imparted to a walker, or a oscillator that periodically changes in a form perceivable by a walker in response to a torque and a change in torque.

22. The control system according to claim 2 wherein the walking oscillator measuring means measures, as a walking oscillator for determining a walking condition, at least one of a leg joint angle, a leg joint angular velocity, a leg joint angular acceleration, and the position of a part of a leg of a walker.

23. The control system according to claim 3, wherein
the walking condition determining means determines whether a walking condition of a walker is a flat walk condition, an ascending walk condition, or a descending walk condition, and
the desired phase difference determining means basically determines a desired phase difference to be 0 for the flat walk condition, − for the ascending walk condition, and + for the descending walk condition.

24. The control system according to claim 3, wherein
the walking condition determining means determines the speed of walk as a walking condition, and
the desired phase difference determining means basically determines a desired phase difference to be − if the walk is fast, or + if the walk is slow.

25. The control system according to claim 9, comprising: storing means for storing a correspondence relationship between a walking condition and a tracing pattern drawn by an n number of walking oscillators in an n-dimensional space, wherein the walking condition determining means determines a walking condition on the basis of the correspondence relationship stored by the storing means and the tracing pattern drawn in the n-dimensional space by the n number of walking oscillators measured by the walking oscillator measuring means.

26. The control system according to claim 15, comprising: storing means for storing a correspondence relationship between a walking condition and a tracing pattern drawn by an n number of walking oscillators in an n-dimensional space, wherein the walking condition determining means determines a walking condition on the basis of the correspondence relationship stored by the storing means and the tracing pattern drawn in the n-dimensional space by the n number of walking oscillators measured by the walking oscillator measuring means.

27. The control system according to claim 9, wherein the walking oscillator measuring means measures, as a walking oscillator for determining a walking condition, at least one of a leg joint angle, a leg joint angular velocity, a leg joint angular acceleration, and the position of a part of a leg of a walker.

28. The control system according to claim 15, wherein the walking oscillator measuring means measures, as a walking oscillator for determining a walking condition, at least one of a leg joint angle, a leg joint angular velocity, a leg joint angular acceleration, and the position of a part of a leg of a walker.

29. The control system according to claim 18, wherein the walking oscillator measuring means measures, as a walking oscillator for determining a walking condition, at least one of a leg joint angle, a leg joint angular velocity, a leg joint angular acceleration, and the position of a part of a leg of a walker.