# Patent

## Title
Control system, oscillation control device, and control signal generation method

## Abstract

There is proposed a control system and the like that can realize a CPG (Central Pattern Generator) network having high controllability. A control system (CPG network) includes a plurality of oscillation control devices (CPG) and one target signal generating device (rhythm generator (RG)). Each CPG and the RG are described by the Van der Pol equation, and an amplitude and a period of an output waveform are substantially independently controlled by an external signal. In order to control a phase difference between the CPGs, the period of each CPG is temporarily controlled through connections that are only conjunctions between the each CPG and the RG.

## Background
TECHNICAL FIELD
The present invention relates to a control system, an oscillation control device, and a control signal generation method, and particularly to a control system and the like including a target signal generating device that generates a target signal and a plurality of oscillation control devices that generate control signals.
BACKGROUND ART
It is said that locomotion of a living organism such as an animal and a human being is generated and controlled by a pattern generator, that is, a neural oscillator CPG (Central Pattern Generator) network that is said to exist in a spinal cord. Hereinafter, a plurality of nonlinear oscillators included in the CPG network are referred to as the CPG in order that the CPG network and the CPG are distinguished from each other.
Many models such as a Matsuoka model have been proposed as a CPG network model, and the CPG network model is applied to real machines (robots) (see Non-Patent Documents 1, 2, and 3). Taga shows that the CPG network and a masculo-skeletal system generate mutual entrainment through sensory input, and Taga performs a simulation of two-legged locomotion. Kimura expands the Taga model to produce a quadruped locomotion robot and realizes some locomotion through irregular terrain by utilizing an ability of entrainment with the CPG network and musclo-skeletal system. As described above, the CPG network is a pattern generator that can explain abilities of irregular ground locomotion and the like which are possessed by a living organism.
The conventional CPG network model is one that simulates a structure and an activity of a nerve. In the conventional CPG network model, the CPGs are mutually connected through connection weights, and the nerve indicating a fatigue degree or the nerve indicating an internal state is designed in one CPG.
On the other hand, a Phaselock Techniques is known as a phase control method in which a phase can be controlled without use of the mutual connection between oscillators. A Phase-Locked Loop (PLL) is a specific example of the Phaselock Techniques and consists of a Phase Detector (PD), a Lowpass Filter (LF), and a Voltage Controlled Oscillator (VCO). A phase difference between two signals is detected in the PD, and a period of the oscillator is controlled in the VCO if needed. Assuming that a control signal is a signal outputted from the VCO while a target signal is a signal that becomes a target, in an operating principle of the PLL, the phase difference between the control signal and the target signal is detected by the PD, the output signal of the PD is inputted to the VCO through the LF, and the period of the control signal is controlled by the VCO until the value of the VCO becomes zero.
Volkovskii et al. propose a CPG network model in which the PLL is used (see Non-Patent Document 4). This CPG network model is not an oscillating system which has a limit cycle, because a sine function is used as the CPG model. Hoppensteadt et al. propose a PLL network model (see Non-Patent Document 5). In this model, a sinusoidal function is used as the VCO, and a model concerning pattern recognition is proposed as application of the PLL network model.
A Van der Pol (VDP) equation is a mathematical model that can describe the theory of vacuum-tube oscillators (see Non-Patent Document 6). There have been proposed some CPG networks in which the VDP equation is used (see Non-Patent Document 7).

    [Non-Patent Document 1] K. Matsuoka, “The dynamic model of binocular rivalry”, Biological Cybernetics, Vol. 49, pp. 201-208, 1984.
    [Non-Patent Document 2] G. Taga and two others, “Self-organized control of bipedal locomotion by neural oscillators in unpredictable environment”, Biological Cybernetics, Vol. 65, pp. 147-159, 1991.
    [Non-Patent Document 3] H. Kimura and two others, “Adaptive Dynamic Walking of a Quadruped Robot on Natural Ground Based on Biological Concepts”, International Journal of Robotics Research, Vol. 26, pp. 475-490, 2007.
    [Non-Patent Document 4] A Volkovskii and five others, “Analog electronic model of the lobster pyloric central pattern generator”, Journal of Physics, Conference Series, Vol. 23, pp. 47-57, 2005.
    [Non-Patent Document 5] Frank C. Hoppensteadt and one other, “Pattern Recognition Via Synchronization in Phase-Locked Loop Neural Networks”, IEEE Transactions on neural networks, Vol. 11, No. 3, 2000.
    [Non-Patent Document 6] Van der Pol, “On relaxation oscillations”, Phil. Mag., No. 2, pp. 987-993, 1926.
    [Non-Patent Document 7] Max S. Dutra and two others, “Modeling of a bipedal locomotor using coupled nonlinear oscillators of Van der Pol”, Biological Cybernetics, No. 88, pp. 286-292, 2003.


DISCLOSURE OF THE INVENTION
Problems to be Solved by the Invention
However, in order to generate locomotion, it is necessary that the CPG be the oscillator (oscillating system which has a limit cycle: expressed by a nonlinear differential equation of second order or more) having the entrainment ability, and it is necessary that the CPG be able to control a period and an amplitude. The CPG network includes the plurality of CPGs, and it is necessary to control a phase difference between the CPGs. Further, it is necessary that the output values of the finally controlled CPGs be identical to one another in amplitudes and periods.
On the other hand, in the conventional CPG network models described in Non-Patent Documents 1, 2, 3, and the like, the CPGs expressed by a high-order nonlinear differential equation have a network structure in which the CPGs are mutually connected through connection weights. Therefore, an output waveform of each CPG is hardly controlled, which causes a problem in that the connection weights and the network structure are hardly designed. For example, the conventional connection weights are designed using trial and error or a neural network such as a genetic algorithm. However, in any case, an amount of calculation and a processing time become huge. When the CPG network model is mounted on a real machine (robot), it is desirable that the amplitude and period of the output waveform of each CPG and the phase difference between the CPGs are independently controlled.
The VCO corresponds to the CPG when the PLL is applied to the CPG network as described in Non-Patent Document 4. Here, it is assumed that the target signals having the desired phase differences can be designed, and it is assumed that the necessary PLLs are prepared. Then, it seemingly appears that the CPG network can be designed because the signals having arbitrary phase differences can be generated and controlled. However, there is a problem with respect to a PD designing method. That is, the PD designing method is classified into subtraction and multiplication. And the subtraction cannot be applied to the CPG because the sinusoidal wave such as a sine function is used as the VCO in the subtraction. Although periodic waves other than sinusoidal waves can also be used as the VCO in the multiplication, the phase is often controlled while the amplitude of each oscillator deviated. Therefore, the multiplication cannot be applied to the CPG network.
Further, the VDP equation is the oscillating system which has a limit cycle, and the VDP equation is the model in which the amplitude and the period can substantially independently be controlled. Accordingly, it is considered that the VDP equation is suitable to the CPG model. However, the models described in Non-Patent Document 7 and the like are the CPG network that does not retain controllability of the period and the amplitude, which is the feature of the VDP equation.
An object of the present invention is to propose a control system and the like that can realize a CPG network having high controllability.
Means for Solving the Problem
An invention as claimed in Claim 1 relates to a control system comprising a target signal generating device that generates a target signal and a plurality of oscillation control devices that generate control signals, wherein each of the oscillation control devices includes a signal generating means for generating a control signal based on a parameter concerning an amplitude and a parameter concerning a period and a regulating means for regulating the parameter concerning the amplitude and/or the parameter concerning the period based on the target signal and the control signal. Here, “a regulating means for regulating the parameter concerning the amplitude and/or the parameter concerning the period” means that the regulating means regulate either the parameter concerning the amplitude or the parameter concerning the period or both of them. The “and/or” is used in the following in the same way.
The invention according to Claim 2 is the control system according to Claim 1, wherein preferably the regulating means of each of the oscillation control devices is independent of a control signal generated by another oscillation control device.
The invention according to Claim 3 is the control system according to Claim 1 or 2, wherein preferably in an ith (i is a natural number of the number of oscillation control devices or less) oscillation control device, the signal generating means generates a signal xi based on equation (eq1), and the regulating means regulates Ai and/or Bi.
The invention according to Claim 4 is the control system according to any one of Claims 1 to 3, wherein preferably the target signal generating device generates a target signal Xi for the ith (i is a natural number of the number of oscillation control devices or less) oscillation control device based on equation (eq4).
The invention according to Claim 5 is the control system according to any one of Claims 1 to 4, wherein preferably the regulating means regulates the parameter concerning the amplitude and/or the parameter concerning the period based on a phase difference between the target signal and the control signal.
The invention according to Claim 6 is the control system according to any one of Claims 1 to 5, wherein preferably the regulating means includes a difference operation means for calculating a difference between the target signal and the control signal.
An invention as claimed in Claim 7 relates to an oscillation control device comprising a regulating means for independently regulating an amplitude and/or a period of a newly generated control signal based on a predetermined target signal and the generated control signal.
An invention as claimed in Claim 8 relates to a control signal generation method for generating a control signal, the method comprising a regulating step of regulating a parameter concerning an amplitude and/or a parameter concerning a period based on a predetermined target signal and a generated control signal, and a signal generating step of generating a new control signal based on the regulated parameter concerning the amplitude and the regulated parameter concerning the period.
The target signal may be generated so as to satisfy the condition defined by the designer. Particularly, in the invention set forth in claim 4, the combination of constants Ci1 and Ci2 of equation (eq4) may automatically be generated so as to satisfy the predetermined condition.
The regulating means of a certain oscillation control device may regulate the amplitude and/or period of the control signal that is to be newly generated based on, for example, the external signal or a control signal of another oscillation control device in addition to the provided target signal.
The present invention may be recognized as a program that realizes the invention set forth in each of the claims in a computer and a recording medium having the program recorded therein.




  
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          1
        
        ]
      
    
    
      
        
      
    
  
  
    
      
        
          
            
              
                
                  d
                  2
                
                ⁢
                
                  x
                  i
                
              
              
                d
                
                  t
                  2
                
              
            
            -
            
              2
              ⁢
              
                ɛ
                ⁡
                
                  (
                  
                    
                      A
                      i
                      2
                    
                    -
                    
                      x
                      i
                      2
                    
                  
                  )
                
              
              ⁢
              
                
                  d
                  
                    x
                    i
                  
                
                
                  d
                  t
                
              
            
            +
            
              
                B
                i
                2
              
              ⁢
              
                x
                i
              
            
          
          =
          0
        
        ,
      
    
    
      
        (
        
          eq
          ⁢
          
            
          
          ⁢
          1
        
        )
      
    
  
  
    
      
        
          
            
              
                
                  d
                  2
                
                ⁢
                
                  x
                  R
                
              
              
                d
                
                  t
                  2
                
              
            
            -
            
              2
              ⁢
              
                ɛ
                ⁡
                
                  (
                  
                    
                      A
                      2
                    
                    -
                    
                      x
                      R
                      2
                    
                  
                  )
                
              
              ⁢
              
                
                  d
                  
                    x
                    R
                  
                
                
                  d
                  t
                
              
            
            +
            
              
                B
                2
              
              ⁢
              
                x
                R
              
            
          
          =
          0
        
        ,
      
    
    
      
        (
        
          eq
          ⁢
          
            
          
          ⁢
          2
        
        )
      
    
  
  
    
      
        
          τ
          =
          
            
              
                max
                ⁡
                
                  (
                  
                    x
                    R
                  
                  )
                
              
              
                max
                ⁡
                
                  (
                  
                    
                      d
                      
                        x
                        R
                      
                    
                    
                      d
                      t
                    
                  
                  )
                
              
            
            =
            
              
                2
                ⁢
                A
              
              
                max
                ⁡
                
                  (
                  
                    
                      d
                      
                        x
                        R
                      
                    
                    
                      d
                      t
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        
          eq
          ⁢
          
            
          
          ⁢
          3
        
        )
      
    
  
  
    
      
        
          X
          i
        
        =
        
          
            
              c
              
                i
                ⁢
                
                  
                
                ⁢
                1
              
            
            ⁢
            
              x
              R
            
          
          +
          
            
              c
              
                i
                ⁢
                
                  
                
                ⁢
                2
              
            
            ⁢
            τ
            ⁢
            
              
                
                  d
                  
                    x
                    R
                  
                
                
                  d
                  t
                
              
              .
            
          
        
      
    
    
      
        (
        
          eq
          ⁢
          
            
          
          ⁢
          4
        
        )
      
    
  




where, Ai, Bi: parameter

    
    
        ε: nonlinear ratio
        A, B: parameter, and
        ci1, ci2: constant.
    
    


Effects of the Invention
According to the invention set forth in each of the claims of the present application (hereinafter referred to as the “present invention”), each oscillation control device independently includes the oscillation system, and each oscillation control device can be controlled by the network structure, which allows the realization of the CPG network having the high controllability.
In addition, in the CPG network realized by the present invention, the CPGs are not mutually connected through the weight connection. Therefore, for example, the connection structure of the CPG network can uniquely be designed based on a phase transition diagram in each gait and the like, and therefore the design of the connection weights can be simplified.
In addition, according to claim 3, as the CPG network which is described by the VDP equations holds the feature of the VDP equation, the amplitude, period, and phase difference of the output signal from each CPG can be substantially independently controlled. Therefore, the CPG network realizes the high controllability CPG network in which the amplitude, the period, and the phase difference are easily designed and controlled.
In addition, according to claim 6, the amplitudes, periods, and phases of the target signal (signal that becomes a target) and control signal (signal to control) are aligned through subtraction processing of output value, which allows the realization of the CPG network.
For example, the present invention can be applied to a robot in which a high priority is set on the controllability.



BRIEF DESCRIPTION OF THE DRAWINGS
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
FIG. 1 shows a schematic block diagram illustrating a control system 1 according to an embodiment of the present invention.
FIG. 2 illustrates simulation results in which an amplitude and a period of an output waveform of a VDP equation are measured.
FIG. 3 shows a block diagram illustrating a CPGi 21 (corresponding to oscillation control devices 51, . . . , 5n of FIG. 1), which is a proposed CPG model.
FIG. 4 illustrates phase transitions in typical gaits of quadruped locomotion animal.
FIG. 5 shows a block diagram illustrating an RG 31 (corresponding to target signal generating device 3 of FIG. 1), which is a proposed rhythm generator.
FIG. 6 illustrates an example of a CPG network structure configured to generate four kinds of gaits.
FIG. 7 is a schematic block diagram illustrating a control system 41 according to another embodiment of the present invention.
FIG. 8 illustrates transitions of output signals xi and bi of each CPG in a walk mode.
FIG. 9 shows the transition of the output signal xi of each CPG. Fig. (a), (b) and (c) show trot mode, bound mode, and gallop mode, respectively.
FIG. 10 shows the transition of the output signal xi of each CPG with various parameter arrangements of A and B in the walk mode.
FIG. 11 illustrates phase transition diagrams in gaits of modes except for the walk, trot, bound, and gallop modes. These gates simultaneously satisfy a condition of the walk, trot, bound, and gallop modes in FIG. 4 and a condition obtained based on some constraints which are derived by biological gaits.
FIG. 12 shows the transition of (a) the output signal xi of each CPG and (b) bi, with a combination of automatically generated ci1 and ci2.
FIG. 13 shows a schematic block diagram illustrating an oscillation control device 51i (i is a natural number equal to nor less) according to still another embodiment of the present invention.
FIG. 14 shows the transition of (a) the output signal xi of each CPG and (b) bi, when phase difference control is performed with an RG up to 1500 steps and performed with an external signal after 1500 steps.
FIG. 15 shows the transition of (a) the output signal xi of each CPG and (b) bi, when the bi is controlled so that the phase difference control is performed with the RG up to 1500 steps and that the sum of the control signal xi except the control signal x2 becomes zero after 1500 steps.



DESCRIPTION OF REFERENCE NUMERALS


    
    
        1 control system
        3 target signal generating device
        51, . . . , 5n oscillation control device
        91, . . . , 9n regulating unit
        111, . . . , 11n signal generating unit
        21 CPGi 
        23 regulating unit
        25 VDPi 
        31 RG
        33 VDPR 
        43 target signal generating device
        471, . . . , 47n oscillation control device
        51i oscillation control device
        53i regulating unit
        55i signal generating unit
    
    


BEST MODE FOR CARRYING OUT THE INVENTION
Embodiments of the present invention will be described below with reference to the drawings.
First Embodiment
FIG. 1 shows a schematic block diagram illustrating a control system 1 according to an embodiment of the present invention.
The control system 1 includes a target signal generating device 3 that generates a target signal, and n oscillation control devices 51, . . . , 5n that generate control signals.
In the embodiment, the target signal generating device 3 and the oscillation control devices 51, . . . , 5n are designed with a VDP equation. First, the VDP equation will be described.
The VDP equation, proposed by Van der Pol in 1926, describes the theory of vacuum-tube oscillators (see Non-Patent Document 6). The VDP equation is one of the simplest models of oscillating systems which have limit cycles, and the VDP equation is expressed by equation (1).
Here, A and B are parameters. ε indicates a nonlinear ratio, When the parameter ε is small (0<ε<<1), the solution of the VDP equation is a sustained sinusoidal oscillation. In the case of ε=0, the VDP equation becomes a harmonic oscillator. According to an approximation method (see Luis A. Pipes and one other, “Mathematics for engineers and physicists”, McGrawHill Education, Third Edition, 1970) proposed by Kryloff and Bogoliuboff, a solution of the harmonic oscillator can be expressed by equation (2) using a and φ. Equation (3) and equation (4) are obtained using a solution of Kryloff and Bogoliuboff.
As can be seen from equation (3), an amplitude of an output value of the VDP equation reaches a stable state in the case of da/dt=0, and the amplitude becomes a=2 A. In an actual simulation, however, the amplitude does not always exactly become 2 A. This is because the approximation is used in this solution.
As can be seen from equation (4), a phase of the output value of the VDP equation does not depend on time.
When a and φ are constant, it can be seen from equation (2) that a period of an output waveform of the VDP equation can be controlled by the parameter B.
Therefore, according to equation (2) and equation (3), for arbitrary constants a and φ, the amplitude and period of the output value of the VDP equation can be independently controlled by the parameters A and B, respectively, the period is inversely proportional to the parameter B, and the amplitude becomes 2 A.
FIG. 2 illustrates simulation results in which the amplitude and period of the output waveform of the VDP equation are measured. ε=0.2 is used, and x=0.5 and dx/dt=0.1 are used as initial values. The VDP equation operated as the damped oscillation in the case of A=0. Oscillatory solution are always obtained from the VDP equation in the range of A<B. As can be seen from FIG. 2, the amplitude and period of the output waveform of the VDP equation can independently be controlled by the parameters A and B. As described above, the VDP equation has a limit cycle, and has the independent controllability of the amplitude and period. Therefore, the VDP equation is useful for designing the CPG network having independent controllability.
Referring to FIG. 1, based on the parameters A and B obtained from the outside, the target signal generating device 3 generates a target signal Xi for the oscillation control device 5i (i is a natural number equal to n or less) by using the VDP equation. In order to align the amplitudes and periods of the oscillation control devices 5i, amplitude and period control terms of the VDP equations are uniformly controlled by the target signals Xi. By this processing, the CPG network in which the waveforms can uniformly be controlled according to an external change can be designed.
The parameters A and B are externally provided to the oscillation control device 5i, and the oscillation control device 5i generates a control signal xi from equation (5) based on parameters Ai and Bi obtained by regulating the parameters A and B. The signal control device 5i includes a regulating unit 9i and a signal generating unit 11i. The regulating unit 9i regulates at least one of the parameters Ai and Bi based on the target signal Xi and the generated control signal xi. The signal generating unit 11i newly generates a control signal xi based on the parameter regulated by the regulating unit 9i.
In FIG. 1, the oscillation control devices 5i, . . . , 5n are independent of one another and generate control signals using target signals X1, . . . , Xn respectively, generated by the target signal generating device 3. In each of the oscillation control devices 51, . . . , 5n, the amplitude and period of the control signal to be generated can substantially independently be regulated. Therefore, for example, in realizing a gait, each leg can independently be controlled by providing the oscillation control device corresponding to each leg, and the control system 1 can deal with the arbitrary number of legs such as two legs and three legs. Thus, the control system 1 can deal with the change of the number of oscillation control devices, and the present invention can realize a robust system. Such a system can be realized by focusing on the nature of motion such as gaits.




  
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          2
        
        ]
      
    
    
      
        
      
    
  
  
    
      
        
          
            
              
                
                  d
                  2
                
                ⁢
                x
              
              
                d
                
                  t
                  2
                
              
            
            -
            
              2
              ⁢
              
                ɛ
                ⁡
                
                  (
                  
                    
                      A
                      2
                    
                    -
                    
                      x
                      2
                    
                  
                  )
                
              
              ⁢
              
                
                  d
                  x
                
                
                  d
                  t
                
              
            
            +
            
              
                B
                2
              
              ⁢
              x
            
          
          =
          0
        
        ,
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          x
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                
                  Bt
                  +
                  φ
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        
          
            
              d
              a
            
            
              d
              t
            
          
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              ɛ
              ⁡
              
                (
                
                  
                    A
                    2
                  
                  -
                  
                    
                      a
                      2
                    
                    4
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        
          
            
              d
              φ
            
            
              d
              t
            
          
          =
          0
        
        ,
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
            
              
                d
                2
              
              ⁢
              
                x
                i
              
            
            
              d
              
                t
                2
              
            
          
          -
          
            2
            ⁢
            
              ɛ
              ⁡
              
                (
                
                  
                    A
                    i
                    2
                  
                  -
                  
                    x
                    i
                    2
                  
                
                )
              
            
            ⁢
            
              
                d
                
                  x
                  i
                
              
              
                d
                t
              
            
          
          +
          
            
              B
              i
              2
            
            ⁢
            
              x
              i
            
          
        
        =
        0.
      
    
    
      
        (
        5
        )
      
    
  




The design of the oscillation control device 5i in which the VDP equation is used will specifically be described below. More generally, the case where only the parameter Bi concerning the period is regulated in designing n CPGs (corresponding to oscillation control devices 5i of FIG. 1), in which the VDP equation is used, constituting the CPG network (corresponding to the control system 1 of FIG. 1) will specifically be described. In other words, in the following description, it is assumed that the parameter Ai concerning the amplitude in each CPG is not regulated and is equal to the externally provided parameter A. As described above, in the present invention, only the parameter concerning the period may be regulated and both of the parameters are not necessarily regulated.
FIG. 3 shows a block diagram illustrating a CPGi 21, which is a proposed CPG model. The CPGi 21 includes a regulating unit 23 corresponding to the regulating unit 9i of FIG. 1 and a VDPi 25 corresponding to the signal generating unit 11i of FIG. 1.
A VDP equation (VDPi) of an ith CPGi is expressed by equation (6).
In this embodiment, the following method is used in order that a phase difference between the CPGs is controlled keeping the controllability of the amplitude and period, which is the feature of the VDP equation. It is assumed that the signal having the desired phase difference from that of the output value xi of each CPG is the target signal Xi and Xi is known. A method will be described in which the period of the CPGi is tentatively controlled so that the phase difference between xi and Xi becomes zero, thereby controlling the phase difference between the CPGs.
The difference between xi and Xi is calculated, and only the period of each CPG is controlled so that the difference becomes zero. This operation can be expressed by equation (7). Here, bi is a parameter that tentatively controls the period of the output waveform of the ith CPGi and k is a coefficient that is used to determine how much the difference between xi and Xi is reflected in the period control. In the following description, it is assumed that the coefficient k is always set to 1.
Next, the parameter Bi (parameter concerning the period) of the CPGi expressed by equation (6) is controlled using bi in equation (8). Regarding the physical meaning of equation (8), the parameter B indicates a natural frequency of each CPG while the parameter bi indicates a phase deviation of the target signal Xi to the output waveform xi of the ith CPGi.
The processing expressed by equations (7) and (8) is repeated until the difference between xi and Xi is eliminated. Through this processing, when the target signal Xi is known, the output value xi can be generated as a signal having the phase identical to that of the target signal Xi, that is, the phase between the CPGs can be controlled.




  
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          3
        
        ]
      
    
    
      
        
      
    
  
  
    
      
        
          
            
              
                
                  d
                  2
                
                ⁢
                
                  x
                  i
                
              
              
                d
                
                  t
                  2
                
              
            
            -
            
              2
              ⁢
              
                ɛ
                ⁡
                
                  (
                  
                    
                      A
                      2
                    
                    -
                    
                      x
                      i
                      2
                    
                  
                  )
                
              
              ⁢
              
                
                  d
                  
                    x
                    i
                  
                
                
                  d
                  t
                
              
            
            +
            
              
                B
                i
                2
              
              ⁢
              
                x
                i
              
            
          
          =
          0
        
        ,
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          
            b
            i
          
          =
          
            k
            ⁡
            
              (
              
                
                  x
                  i
                
                -
                
                  X
                  i
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          B
          i
        
        =
        
          B
          +
          
            
              b
              i
            
            .
          
        
      
    
    
      
        (
        8
        )
      
    
  




The design of the target signal generating device 3 (oscillator that generates the target signal, hereinafter the oscillator is referred to as the rhythm generator (RG)) in which the VDP equation is used will specifically be described.
In order to set a unique target signal to each CPG, it is necessary to define the gait to be designed. In this embodiment, typical gait signals of the quadruped locomotion animal are generated as illustrated in FIG. 4. FIG. 4 illustrates phase transitions in typical gaits of the quadruped locomotion animal. FIG. 4(a) illustrates the phase transition diagram in a walk mode, FIG. 4(b) illustrates the phase transition diagram in a trot mode, FIG. 4(c) illustrates the phase transition diagram in a bound mode, and FIG. 4(d) illustrates the phase transition diagram in a gallop mode. In FIG. 4, the letters LF, RF, LH, and RH denote a left foreleg, a right foreleg, a left hind leg, and a right hind leg, respectively. The arrow and equal sign indicate the direction of the phase difference, and in phase, respectively.
As can be seen from FIG. 4, only four types of the phase differences of 0, π/2, π, and 3π/2 exist between the CPGs that should be controlled. A method for generating the target signal for expressing the phase difference will be described.
As described above in equation (7), xi−Xi=0 is required as a condition of the target signal after the phase difference between the output value xi and the target signal Xi of each CPG becomes zero. Besides, it is necessary that the target signals with the waveforms having the four types of phase differences of 0, π/2, π, and 3π/2 from the phase of one reference signal are generated.
A method for designing the rhythm generator with the VDP equation will be described.
FIG. 5 shows a block diagram illustrating an RG 31, which is the proposed rhythm generator. In FIG. 5, a VDPR 33 is the VDP equation constituting the RG, and an output signal of the VDPR 33 is expressed with a subscript R. As the RG is also designed using the VDP equation, the target signal Xi can be formed by using ±xR and ±τdxR/dt which are calculated in the VDP equation. Therefore, it is possible to express the four types of the phase differences described above from the phase of the reference signal xR. For example, letting X1=XR and X2=τdxR/dt leads to the phase difference of π/2 between x1 and x2 after the control. Here, it is assumed that the VDP equation is used within a range where the VDP equation has oscillatory solutions.
In FIG. 5, based on the phase transition diagram in each gait of FIG. 4, the target signal selecting unit 35 selects one output value Xi from the four inputted signals for each of the n CPGs. The values identical to those of the parameters A and B are inputted to each CPG, and one of ±τdxR/dt and ±xR selected by the target signal selecting unit 35 is used as the target signal Xi. This processing is expressed by a mathematical expression in equation (9).
Here, in this embodiment, ci1 and ci2 take 0, −1, or 1. Additionally one of ci1 and ci2 always takes zero, and ci1 and ci2 do not become zero at the same time. τ in equation (9) has two roles. One of the roles is a time constant to match necessarily the dimension of Xi expressed by one of xR and dxR/dt with the dimension of xR. The other role is to adjust the amplitudes of the output waveforms of dxR/dt to always keep the relationship of the amplitudes of xR and dxR/dt constant even if the parameter A and B change. This is to deal with the circumstances that the amplitudes of the output waveforms of the xR and dxR/dt are different, which can occur depending on the parameters A and B. Therefore, τ is calculated by equation (10). The max function indicates the amplitudes of the output waveform xR and τdxR/dt. The amplitude of the output waveform xR becomes 2 A in theory.




  
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          4
        
        ]
      
    
    
      
        
      
    
  
  
    
      
        
          
            X
            i
          
          =
          
            
              
                c
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              ⁢
              
                x
                R
              
            
            +
            
              
                c
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
              
              ⁢
              τ
              ⁢
              
                
                  d
                  
                    x
                    R
                  
                
                
                  d
                  t
                
              
            
          
        
        ,
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        τ
        =
        
          
            
              max
              ⁡
              
                (
                
                  x
                  R
                
                )
              
            
            
              max
              ⁡
              
                (
                
                  
                    d
                    
                      x
                      R
                    
                  
                  
                    d
                    t
                  
                
                )
              
            
          
          =
          
            
              
                2
                ⁢
                A
              
              
                max
                ⁡
                
                  (
                  
                    
                      d
                      
                        x
                        R
                      
                    
                    
                      d
                      t
                    
                  
                  )
                
              
            
            .
          
        
      
    
    
      
        (
        10
        )
      
    
  




The RG and each CPG are designed using the VDP equation. The structure of each CPG differs from the structure of the RG in that a feedback signal from a sensor mounted on the actual robot and the target signal Xi are inputted.
In this embodiment, although attention is paid to the Phaselock Techniques in order to solve the problem of mutual connection and connection weights of the conventional CPG network, subtraction processing of output values of the two oscillators is performed in order to align the amplitudes and phases of the two oscillators in this embodiment. The output value of the subtraction processing becomes zero when the amplitudes and phases of the target signal are equal to those of the control signal. However, the output value of the subtraction processing is not always proportional to the phase difference between the two oscillators. Therefore, the operation of the subtraction is not functional equivalent of the PD. That is, the Phaselock Techniques in which the subtraction processing is performed is not the PLL.
Next, the design of the CPG network for the quadruped locomotion animal will be described below with reference to the CPGi of FIG. 3 and the RG of FIG. 5.
FIG. 6 illustrates an example of a CPG network structure configured to generate four kinds of gaits. FIG. 6(a) illustrates the CPG network structure in the walk mode, FIG. 6(b) illustrates the CPG network structure in the trot mode, FIG. 6(c) illustrates the CPG network structure in the bound mode, and FIG. 6(d) illustrates the CPG network structure in the gallop mode.
First, one RG and the necessary number (n) of CPGs are prepared. In this embodiment, the four (n=4) CPGs are required because the locomotion signals for the quadruped locomotion are generated and controlled. The parameters A and B are inputted to the RG and each CPG in order to determine the amplitude and the period.
The output signal of each CPG is allocated as the control signal to each leg. In this embodiment, the CPG1, CPG2, CPG3, and CPG4 are made correspondent to the LF, RF, LH, and RH, respectively.
The target signal Xi of each CPG is determined based on the phase transition diagram in each gait of FIG. 4. A method for determining Xi (1≦i≦n) in the walk mode will be described below.
Because xi corresponds to LF which works as the reference signal for calculating the phase difference between CPG1 and the other CPGs, xR is selected as the target signal X1.
Because x2, x3, and x4 correspond to RF, LF and RH having the phase differences of π, 3π/2, and π/2, respectively, relative to LF, −xR, τdxR/dt, and −τdxR/dt are inputted as the target signals X2, X3, and X4, respectively.
When the CPG network in the walk mode is compared with equation (9), c11=1, c21=−1, c31=0, c41=0, c12=0, c22=0, c32=1, and c42=1 are obtained. Similarly, the network can uniquely be designed for the trot, bound, and gallop modes based on the phase transition diagram in each gait. Table 1 collectively shows the combinations of ci1 and ci2 in each gait.
















TABLE 1









walk

trot

bound

gallop

















i
xi
ci1
ci2
ci1
ci2
ci1
ci2
ci1
ci2



















1
x1
1
0
1
0
1
0
1
0


2
x2
−1
0
−1
0
1
0
0
−1


3
x3
0
1
−1
0
−1
0
0
1


4
x4
0
−1
1
0
−1
0
−1
0









FIG. 7 is a schematic block diagram illustrating a control system 41 according to another embodiment of the present invention. As illustrated in FIG. 10, for example, a target signal generating device 43 may be the device which generates basic signals xR and τdxR/dt based on the VDR equation. A target signal selecting device 45i may be the device which perform the calculation of equation (9) for the signal xR and τdxR/dt using the constants ci1 and ci2. And n oscillation control devices 47i may be the device which generate the output signal xi based on the signals selected by the target signal selecting device 45i.
As still another embodiment, referring to FIG. 7, the target signal generating device 43 may be the device which generates the signals (for example, four signals xR, −xR, τdxR/dt, and −τdxR/dt) necessary for each of the oscillation control devices 47i, and the target signal selecting device 45i may be the device which select one signal from the signals if necessary.
Thus, the target signal generating device (rhythm generator) may be recognized as one that generates the basic signal as illustrated in FIG. 7 (in the embodiment described above, two-signal type of xR and τdxR/dt), the target signal generating device may be recognized as one that includes a portion that performs operation of multiplication by −1 to the basic signal to generate the necessary signal (in the embodiment described above, four-signal type), and further the target signal generating device may be recognized as one that includes a portion that specifies the signal to be provided to each oscillation control device (CPG) as illustrated in FIG. 5 (in the embodiment described above, n-signal type).
Second Embodiment
Next, referring to FIGS. 8, 9, and 10, effectiveness of the proposed model will be described by simulations in which Matlab is used. A 4th order Runge-Kutta method is used as the solution of the nonlinear differential equation.
FIG. 8 shows the simulation results of the transitions of the output signals xi and bi of each CPG in the walk mode. FIG. 8(a) illustrates the transition of the output signal xi of each CPG, and FIG. 8(b) illustrates the transition of bi. ε=0.2 and k=1 are used, and xR=0.1, x1=0.1, x2=−0.5, x3=0.3, x4=0.7, dxR/dt=0.1, dx1/dt=0.1, dx2/dt=0.3, dx3/dt=0.3, and dx4/dt=0.2 are used as initial values. A=0.5 and B=1 are used as the parameters to determine the amplitude and the frequency. bi indicates the phase deviation of xi from the phase of the target signal Xi. In these simulations, the operation in which the value of bi of 0.3 or less is omitted was added. The reason why this processing is performed is that bi in equation (8) does not always completely fall to zero, depending on the parameters A and B, even after the phases of xi and Xi are synchronous with each other. As can be seen from the output result, the period of each CPG is controlled by bi so that the phase difference between the output waveforms of the CPGs becomes π/2. As can also be seen from the output result, bi becomes zero after the output value of each CPG is controlled in each target signal.
FIG. 9 shows the transition of the output signal xi of each CPG. Fig. (a), (b) and (c) show trot, bound, and gallop modes, respectively. FIG. 9(a) illustrates the simulation result in the trot mode, FIG. 9(b) illustrates the simulation result in the bound mode, and FIG. 9(c) illustrates the simulation result in the gallop mode. ε, k, A, B, and the initial values of xR, dxR/dt, xi, and dxi/dt are the same with those of the simulation conditions in the walk mode. As can be seen from FIG. 9, the signals having the phase differences identical to those of the phase transition diagram in each gait of FIG. 4 are generated and controlled. The identical parameters A and B are used in the output results of FIG. 8(a), FIG. 9(a), FIG. 9(b), and FIG. 9(c), so that the amplitudes and the periods are identical in each gait. Regarding the output values, it can also be seen that the waveforms (x2, x3, and x4) having the desired phase differences relative to x1 that is the output value of the CPG1 are generated and controlled.
FIG. 10 illustrates experimental result in the walk mode under several values of A and B. FIG. 10(a) illustrates the experimental result in the case of A=0.5 and B=0.6, FIG. 10(b) illustrates the experimental result in the case of A=0.5 and B=1, and FIG. 10(c) illustrates the experimental result in the case of A=0.8 and B=1. The above-described values are used as e, k, and the initial values of xR, dxR/dt, xi, and dxi/dt. As can be seen from FIG. 10, the amplitude and the period of the output waveform of each CPG can be independently controlled even in the network structure. Besides, as described above, it can be seen that the amplitude of the output waveform of the VDP equation reaches double of the parameter Awhile the period is inversely proportional to the parameter B. Thus, it can be seen that the amplitude and period of the output waveform can be independently controlled by the parameters A and B. The measured result of a period T of each signal is (a) 105 (steps), (b) 62 (steps), and (c) 63 (steps).
As described above, in this embodiment, each CPG constituting the CPG network has a limit cycle, and only the phase difference between the CPGs is controlled by the RG. Further, each CPG is expressed using the VDP equation, which is the nonlinear differential equation where the amplitude and frequency of the output value can substantially independently be controlled. Therefore, the CPG network having the high controllability can be designed.
Third Embodiment
Another target signal selecting device (see the target signal selecting unit 35 of FIG. 5 and the target signal selecting devices 451, . . . , 45n of FIG. 7) of another embodiment will be described below.
The gait satisfying a certain condition can freely be switched by freely selecting the combination of coefficients ci1 and ci2 of equation (9) under the condition. Therefore, some rules are found from the combinations of ci1 and ci2 of the typical gaits (for example, see FIG. 4 and Table 1) and these rules are set as the conditions.
An example of the rules will be described with reference to Table 1. x1 is fixed as the reference signal (c11=1 and c12=0). From Table 1, the summation of ci2 and the summation of ci2 become zero (Σci1=Σci2=0). The sum of ci1 and ci2 becomes 1 or −1 (ci1+ci2=1 or −1). Therefore, as described later, three conditions are produced based on constraints (the following condition 1 and condition 2) in the CPG network and a constraint (condition 3) in the gait of the living organism, and the target signal selecting device determines the combination of ci1 and ci2 under the conditions.
The first condition (condition 1) is that the sum of all the generated waveforms is zero (Σci1=Σci2=0). The second condition (condition 2) is that xi and dxi/dt are not connected at the same time (ci1+ci2=1 or −1). The third condition (condition 3) is that the foreleg and the hind leg on the identical side do not become in phase (c21≠c41 and c31≠1).
Four gaits simultaneously satisfying these three conditions exist except the walk, trot, bound, and gallop modes (see FIG. 4 and Table 1). Table 2 illustrates the combinations of ci1 and ci2 of the four gaits except the walk, trot, bound, and gallop modes. FIG. 11 illustrates transition diagrams in the gaits of modes except for the walk, trot, bound, and gallop modes illustrated in FIG. 4.
Therefore, there are eight types of parameter sets satisfying the three conditions, and the gaits satisfying the conditions defined by the designer can be obtained.
















TABLE 2









gait 1

gait 2

gait 3

gait 4

















i
xi
ci1
ci2
ci1
ci2
ci1
ci2
ci1
ci2



















1
x1
1
0
1
0
1
0
1
0


2
x2
−1
0
0
−1
0
1
0
1


3
x3
0
−1
−1
0
−1
0
0
−1


4
x4
0
1
0
1
0
−1
−1
0









Referring to FIG. 12, the effectiveness of the proposed model is illustrated by a simulation. The experiment and simulation has been conducted under the conditions that the combination of ci1 and ci2 simultaneously satisfying the conditions 1, 2, and 3 is generated using a random number. Here, the random number is selected from −1, 0, and 1. The CPG network structure is changed by generating the random number every 1500 step.
FIG. 12 shows the transition of (a) the output signal xi of each CPG and (b) bi, with a combination of automatically generated ci1 and ci2. In FIG. 12, a horizontal axis indicates the number of steps. The automatically generated combination of ci1 and ci2 corresponds to gallop mode up to 1500 steps, the mode of gait 4 in the range of 1501 to 3000 steps, and the trot mode in the range of 3001 to 4500 steps. As can be seen from FIG. 12, the gait satisfying the specified condition can automatically be generated. By using the eight types of signals generated under the conditions, stable locomotions can be expected because one leg moves while the other three legs always stand on the ground. Therefore, these signals controlled by these phase relationships can be regarded as the gaits.
Fourth Embodiment
Referring to FIGS. 13, 14 and 15, a mechanism that generates various phase differences of another embodiment will be described below. As described later, the generation of such various phase differences enables to slightly adjust the phase difference between the basic signals generated by the rhythm generator RG and adapt to the change of the number of legs.
FIG. 13 shows a schematic block diagram illustrating an oscillation control device 51i (i is a natural number equal to the number n of oscillation control devices or less in the system) according to still another embodiment of the present invention. Similarly to the oscillation control device 5i of FIG. 1, the external rhythm generator inputs the target signal Xi to the oscillation control device 51i. In FIG. 13, a signal (external signal) ei is externally inputted to a regulating unit 53i. Control signals x1, . . . , xi−1, xi+1, . . . , xn are inputted to the regulating unit 53i from other oscillation control devices.
First, the case where the regulating unit 53i regulates the phase difference using the external signal ei will be described. In this embodiment, Bi=B+bi+ei is used instead of equation (8). FIG. 14 shows the transition of (a) the output signal xi of each CPG and (b) bi, when the phase difference control is operated with an RG up to 1500 steps and operated with the external signal after 1500 steps. The gait becomes the bound mode in the range of 800 to 1500 steps. The regulating unit 53i regulates Bi based on the external signal ei (the external signal e2=0.2 in the range of 1500 to 1550 steps, and the external signal ei=0 in other ranges). As can be seen from FIG. 14, the phase difference is modified. Thus, as can be seen from FIG. 14, the waveform having the arbitrary phase difference can be generated by utilizing the external signal.
Then, the case where the regulating unit 53i regulates the phase difference using the control signal of other oscillation control devices when, for example, the number of oscillation control devices is changed will be described. A relationship that the summation of xi is zero (Σxi=0) always holds in a waveform having a deviation of (one period×1/n). Therefore, the waveform having the phase difference of the deviation of (one period×1/n) can be generated and controlled. FIG. 15 shows the transition of (a) the output signal xi of each CPG and (b) bi, when bi is controlled so that the phase difference control is performed with the RG up to 1500 steps and that the sum of the control signal xi becomes zero after 1501 steps. The gait is the walk mode in the range of 1000 to 1500 steps, and the phase difference is set to 2π/4. Now, it is assumed that the signal output unit for the right foreleg (x2) is broken. In such cases, only the left foreleg (x1), the left hind leg (x3), and the right hind leg (x4) can generate the control signals. Therefore, after 1500 steps, x2 is set to zero without use of the RG, and bi is controlled so that the summation of xi becomes zero, which allows the phase difference to be set to 2π/3. As can be seen from FIG. 15, the oscillation control device 51i can adapt to the change of the number of legs.
Here, for example, Japanese Patent Publication Laid-Open No. 2006-289602 discloses a robot apparatus in which periodic motion of a moving unit is controlled with an oscillator having the feature of entrainment. Although the phase differences are aligned with each other, the robot apparatus differs from the present invention in that the robot apparatus does not have the independent controllability of the amplitude and period, that the entrainment of the oscillator is utilized (in the present invention, the period can be controlled directly), and that the Matsuoka model is used in the signal generating unit.
The VDP equation used in the CPG network of the embodiments is used only in the range where the VDP equation has the stable state. That is, the amplitude and amplitude of the sinusoidal oscillation that can be expressed by the proposed model are restricted. However, the moving range of each leg of the actual quadruped locomotion robot is mechanically determined, and it is not necessary to move the leg beyond the moving range. Therefore, the problem can be solved by considering the relationship between the output range of the output waveform that can be expressed by one VDP equation and the moving range of each leg.
In the quadruped locomotion animals, locomotion speed is increased as the gait makes the transition in the order of the walk, trot, bound, and gallop modes. In contrast, in the proposed technique, because the gait speed (period) is kept constant even after the gait is changed, it is necessary that the parameter B in the VDP equation is controlled in order to determine the gait speed when the gait transition is generated.
However, in the embodiments, the amplitude and the period can be substantially independently controlled, and the locomotion speed (period) can be controlled while the phase difference and the amplitude are fixed. Therefore, the embodiments can deal with the gait transitions.
Furthermore in the actual locomotion animal, when one leg stumbles over something, other legs are influenced by the stumbled leg to sophisticatedly change motion, and thereby the animal continues the stable gait. In the embodiments, although the CPGs do not have the mutual connection, sensor feedback, for example, can be effectively used. The CPG network can be more effective by feeding a signal of sensing the body balance back to each CPG, and by having the CPGs influence one another based on the body balance. For example, when the external device such as the sensor is added, the gait period can be regulated while the gait is kept constant (locomotion speed can be regulated). Besides, the designer can generate a preferred gait by, for example, changing the length of stride like extending it or by defining an arbitrary energy function. Further, the lengths of stride of only a part of legs can be changed, or the signal can be generated in order to generate the locomotion adapting to the change of the number of legs. These achievements have been hardly realized in the conventional CPG network.

## Claims

1. A control System comprising:
a plurality of oscillation control devices that generate control signals; and
a target signal generating device that generates a target signal for each of the oscillation control devices, the target signal entraining the oscillation control devices to synchronized movement, wherein
each of the oscillation control devices is an oscillating system which has a limit cycle and which operates independently of the other oscillation control devices, including:
a signal generating unit that generates a control signal based on a parameter concerning an amplitude and a parameter concerning a period of the control signal; and
a regulating unit that regulates the parameter concerning the amplitude and/or the parameter concerning the period based on the target signal and the control signal and thereby
an amplitude and a period of a control signal generated by the oscillation control device can be independently controlled,
assuming that Xi is the target signal for the ith (i is a natural number equal to a number of the oscillation control devices or less) oscillation control device,
the target signal generating device generates a signal xR obtained by a following equation:







  
    
      
        
          
            d
            2
          
          ⁢
          
            x
            R
          
        
        
          d
          
            t
            2
          
        
      
      -
      
        2
        ⁢
        
          ɛ
          ⁡
          
            (
            
              
                A
                2
              
              -
              
                x
                R
                2
              
            
            )
          
        
        ⁢
        
          
            d
            
              x
              R
            
          
          
            d
            t
          
        
      
      +
      
        
          B
          2
        
        ⁢
        
          x
          R
        
      
    
    =
    0
  
  ,




where A is a parameter concerning an amplitude of the target signal,
B is a parameter concerning a period of the target signal, and
ε is a predetermined constant indicating a nonlinear ratio, and
the target signal generating device generates the target signal Xi for the ith oscillation control device based on the following equation:





  
    
      X
      i
    
    =
    
      
        
          c
          
            i
            ⁢
            
              
            
            ⁢
            1
          
        
        ⁢
        
          x
          R
        
      
      +
      
        
          c
          
            i
            ⁢
            
              
            
            ⁢
            2
          
        
        ⁢
        τ
        ⁢
        
          
            d
            
              x
              R
            
          
          
            d
            t
          
        
      
    
  
  ,




where ci1 and ci2 are constants indicating phase differences of the ith oscillation control device, and constant expressed by the following equation:




  
    τ
    =
    
      
        
          max
          ⁡
          
            (
            
              x
              R
            
            )
          
        
        
          max
          ⁡
          
            (
            
              
                d
                
                  x
                  R
                
              
              
                d
                t
              
            
            )
          
        
      
      =
      
        
          2
          ⁢
          
            
          
          ⁢
          A
        
        
          max
          ⁡
          
            (
            
              
                d
                
                  x
                  R
                
              
              
                d
                t
              
            
            )
          
        
      
    
  
  ,




and
constants ci1 and ci2 satisfy three conditions, that is
ΣCi1=ΣCi2 as a first condition,
Ci1+CiZ=1 or −1 as a second condition,
C21≠C41 and C11≠C31 as a third condition.

2. The control system according to claim 1, wherein the regulating unit of the ith oscillation control device sets the parameter Ai concerning the amplitude to an externally provided parameter A,
the regulating unit includes a difference operation unit that calculates a difference between the control signal xi and the target signal Xi,
the regulating unit calculates a parameter bi by multiplying a coefficient k and a difference xi−Xi calculated by the difference operation unit, and the regulating unit calculates the parameter Bi concerning the period by adding the externally provided parameter B and bi.

3. The control system according to claim 1, wherein a number of a plurality of oscillation control devices is a number of plurality of legs of a robot,
the plurality of oscillation control devices correspond to the plurality of legs of the robot, respectively,
the control signal generated by each of the oscillation control devices is used to control motion of each leg of the robot to realize a gait, and
the regulating unit of each of the oscillation control devices performs regulating processing independently of a control signal generated by another oscillation control device.

4. A control signal generation method in a control system including a plurality of oscillation control devices that generate control signals, each of the oscillation control devices being an oscillating system which has a limit cycle and which operates independently of the other oscillation control devices, and a target signal generating device that generates a target signal for each of the plurality of the oscillation control devices,
the method comprising:
a step in which the target signal generating device generates a signal xR obtained by a following equation:





  
    
      
        
          
            d
            2
          
          ⁢
          
            x
            R
          
        
        
          d
          
            t
            2
          
        
      
      -
      
        2
        ⁢
        
          ɛ
          ⁡
          
            (
            
              
                A
                2
              
              -
              
                x
                R
                2
              
            
            )
          
        
        ⁢
        
          
            d
            
              x
              R
            
          
          
            d
            t
          
        
      
      +
      
        
          B
          2
        
        ⁢
        
          x
          R
        
      
    
    =
    0
  
  ,




where A: a parameter concerning an amplitude of the target signal,
B: a parameter concerning a period of the target signal, and
ε: a predetermined constant indicating a nonlinear ratio,
in which the target signal generating device generates a target signal Xi for the ith (i is a natural number equal to a number of the oscillation control devices or less) oscillation control device based on a following equation:






  
    
      X
      i
    
    =
    
      
        
          c
          
            i
            ⁢
            
              
            
            ⁢
            1
          
        
        ⁢
        
          x
          R
        
      
      +
      
        
          c
          
            i
            ⁢
            
              
            
            ⁢
            2
          
        
        ⁢
        τ
        ⁢
        
          
            d
            
              x
              R
            
          
          
            d
            t
          
        
      
    
  
  ,




where ci1 and Ci2: constants indicating phase differences of the ith oscillation control device, and
τ: a constant expressed by a following equation:





  
    τ
    =
    
      
        
          max
          ⁡
          
            (
            
              x
              R
            
            )
          
        
        
          max
          ⁡
          
            (
            
              
                d
                
                  x
                  R
                
              
              
                d
                t
              
            
            )
          
        
      
      =
      
        
          2
          ⁢
          
            
          
          ⁢
          A
        
        
          max
          ⁡
          
            (
            
              
                d
                
                  x
                  R
                
              
              
                d
                t
              
            
            )
          
        
      
    
  
  ,




and
in which a signal generating unit of the ith oscillation control device generates a control signal xi based on a parameter concerning an amplitude and a parameter concerning a period of the control signal;
a regulating step in which a regulating unit of the ith oscillation generating device regulates the parameter concerning the amplitude and/or the parameter concerning the period based on the target signal Xi and the control signal xi; and
a signal generating step in which the signal generating unit of the ith oscillation control device generates a new control signal xi based on the regulated parameter concerning the amplitude and the regulated parameter concerning the period, wherein
an amplitude and a period of the control signal xi to be generated by the oscillation control device can independently be controlled.