Abstract:
A self-organizing control system suitable for nonlinear control of a physical object is described. The control system calculates the entropy production difference between a time differentiation (dS u /dt) of the entropy inside the controlled object and a time differentiation (dS c /dt) of the entropy given to the controlled object from a PID controller that controls the object. The entropy production difference is used to generate an evolving control rule by using the entropy production difference as a performance function for a genetic optimizer.

Description:
BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     The disclosed invention relates to generally control systems, and more particularly to control systems using artificial intelligence to control a nonlinear plant. 
     2. Description of the Related Art 
     Feedback control systems are widely used to maintain the output of a dynamic system at a desired value in spite of external disturbance forces that would move the output away from the desired value. For example, a household furnace controlled by a thermostat is an example of a feedback control system. The thermostat continuously measures the air temperature of the house, and when the temperature falls below a desired minimum temperature, the thermostat turns the furnace on. When the furnace has warmed the air above the desired minimum temperature, then the thermostat turns the furnace off. The thermostat-furnace system maintains the household temperature at a constant value in spite of external disturbances such as a drop in the outside air temperature. Similar types of feedback control are used in many applications. 
     A central component in a feedback control system is a controlled object, otherwise known as a process “plant,” whose output variable is to be controlled. In the above example, the plant is the house, the output variable is the air temperature of the house, and the disturbance is the flow of heat through the walls of the house. The plant is controlled by a control system. In the above example, the control system is the thermostat in combination with the furnace. The thermostat-furnace system uses simple on-off feedback control to maintain the temperature of the house. In many control environments, such as motor shaft position or motor speed control systems, simple on-off feedback control is insufficient. More advanced control systems rely on combinations of proportional feedback control, integral feedback control, and derivative feedback control. Feedback that is the sum of proportional plus integral plus derivative feedback is often referred to as PID control. 
     The PID control system is a linear control system that is based on a dynamic model of the plant. In classical control systems, a linear dynamic model is obtained in the form of dynamic equations, usually ordinary differential equations. The plant is assumed to be relatively linear, time invariant, and stable. However, many real-world plants are time varying, highly nonlinear, and unstable. For example, the dynamic model may contain parameters (e.g., masses, inductances, aerodynamic coefficients, etc.) which are either poorly known or depend on a changing environment. If the parameter variation is small and the dynamic model is stable, then the PID controller may be sufficient. However, if the parameter variation is large, or if the dynamic model is unstable, then it is common to add adaptation or intelligent (AI) control to the PID control system. 
     AI control systems use an optimizer, typically a nonlinear optimizer, to program the operation of the PID controller and thereby improve the overall operation of the control system. The optimizers used in many AI control systems rely on a genetic algorithm. Using a set of inputs, and a fitness function, the genetic algorithm works in a manner similar to process of evolution to arrive at a solution which is, hopefully, optimal. The genetic algorithm generates sets of chromosomes (corresponding to possible solutions) and then sorts the chromosomes by evaluating each solution using the fitness function. The fitness function determines where each solution ranks on a fitness scale. Chromosomes which are more fit, are those chromosomes which correspond to solutions that rate high on the fitness scale. Chromosomes which are less fit, are those chromosomes which correspond to solutions that rate low on the fitness scale. Chromosomes that are more fit are kept (survive) and chromosomes that are less fit are discarded (die). New chromosomes are created to replace the discarded chromosomes. The new chromosomes are created by crossing pieces of existing chromosomes and by introducing mutations. 
     The PID controller has a linear transfer function and thus is based upon a linearized equation of motion for the plant. Prior art genetic algorithms used to program PID controllers typically use simple fitness functions and thus do not solve the problem of poor controllability typically seen in linearization models. As is the case with most optimizers, the success or failure of the optimization often ultimately depends on the selection of the performance (fitness) function. 
     Evaluating the motion characteristics of a nonlinear plant is often difficult, in part due to the lack of a general analysis method. Conventionally, when controlling a plant with nonlinear motion characteristics, it is common to find certain equilibrium points of the plant and the motion characteristics of the plant are linearized in a vicinity near an equilibrium point. Control is then based on evaluating the pseudo (linearized) motion characteristics near the equilibrium point. This technique works poorly, if at all, for plants described by models that are unstable or dissipative. 
     SUMMARY OF THE INVENTION 
     The present invention solves these and other problems by providing a new AI control system. Unlike prior AI control systems, the new AI control system is self-organizing and uses a new fitness (performance) function ƒ that is based on the physical law of minimum entropy. The self-organizing control system may be used to control complex plants described by nonlinear, unstable, dissipative models. The self-organizing control system further provides a physical measure of control quality based on physical law of minimum production entropy in an intelligent control system and in complex dynamic behavior control of a plant. 
     In one embodiment, the invention includes a method for controlling a nonlinear object (a plant) by obtaining an entropy production difference between a time differentiation (dS u /dt) of the entropy of the plant and a time differentiation (dS c /dt) of the entropy provided to the plant from a controller. A genetic algorithm that uses the entropy production difference as a fitness (performance) function evolves a control rule for a low-level controller, such as a PID controller. The nonlinear stability characteristics of the plant are evaluated using a Lyapunov function. The evolved control rule may be corrected using further evolutions. 
     In some embodiments, the control method may also include evolving a control rule relative to a variable of the controller by means of a genetic algorithm. The genetic algorithm uses a fitness function based on a difference between a time differentiation of the entropy of the plant (dS u /dt) and a time differentiation (dS c /dt) of the entropy provided to the plant. The variable may be corrected by using the evolved control rule. 
     In another embodiment, the invention comprises an AI control system adapted to control a nonlinear plant. The AI control system includes a simulator configured to use a thermodynamic model of a nonlinear equation of motion for the plant. The thermodynamic model is based on a Lyapunov function (V), and the simulator uses the function V to analyze control for a state stability of the plant. The AI control system calculates an entropy production difference between a time differentiation of the entropy of said plant (dS u /dt) and a time differentiation (dS c /dt) of the entropy provided to the plant by a low-level controller that controls the plant. The entropy production difference is used by a genetic algorithm to obtain an adaptation function in which the entropy production difference is minimized. The genetic algorithm provides a teaching signal to a fuzzy logic classifier that determines a fuzzy rule by using a learning process. The fuzzy logic controller is also configured to form a control rule that sets a control variable of the low-level controller. 
     In one embodiment, the low-level controller is a linear controller such as a PID controller. The learning processes may be implemented by a fuzzy neural network configured to form a look-up table for the fuzzy rule. 
     In yet another embodiment, the invention comprises a new physical measure of control quality based on minimum production entropy and using this measure for a fitness function of genetic algorithm in optimal control system design. This method provides a local entropy feedback loop in the control system. The entropy feedback loop provides for optimal control structure design by relating stability of the plant (using a Lyapunov function) and controllability of the plant (based on production entropy of the control system). The control system is applicable to all control systems, including, for example, control systems for mechanical systems, biomechanical systems, robotics, electro-mechanical systems, etc. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The advantages and features of the disclosed invention will readily be appreciated by persons skilled in the art from the following detailed description when read in conjunction with the drawings listed below. 
     FIG. 1 is a block diagram showing an example of an AI control method in the prior art. 
     FIG. 2 is a block diagram showing an embodiment of an AI control method in accordance with one aspect of the present invention. 
     FIG. 3 is a front view showing a robot unicycle used in a simulation of the present invention. 
     FIG. 4 is a side view showing the robot unicycle shown in FIG.  3 . 
     FIG. 5 is a graph showing the time changes in angle for each part of the robot unicycle as a result of a simulation of an AI posture control system for the unicycle. 
     FIG. 6 is a graph showing the time changes in angular velocity for each part of the robot unicycle during the posture control simulation. 
     FIG. 7 is a graph showing a phase surface for the roll angle and angular velocity during the posture control simulation. 
     FIG. 8 is a graph showing a phase surface for the pitch angle and angular velocity during the posture control simulation. 
     FIG. 9 is a graph showing the time changes in angle for each part of the robot unicycle during a posture control simulation using an AI system of the prior art. 
     FIG. 10 is a graph showing the time changes in angular velocity during the prior art posture control simulation. 
     FIG. 11 is a graph showing a phase surface for the roll angle and angular velocity during the prior art posture control simulation. 
     FIG. 12 is a graph showing a phase surface for the pitch angle and angular velocity during the prior art posture control simulation. 
     FIG. 13 shows another example of the simulation for the posture control indicating the effect of the control method in the present invention. 
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
     A feedback control system is commonly used to control an output variable of a process or plant in the face of some disturbance. Linear feedback control systems typically use combinations of proportional feedback control, integral feedback control, and derivative feedback control. Feedback that is the sum of proportional plus integral plus derivative feedback is often referred to as PID control. The Laplace transform of an output u(s) of a PID controller is given by:                u        (   s   )       =         G        (   s   )            e        (   s   )         =       [       k   1     +       k   2     s     +       k   3        s       ]          e        (   s   )                   (   1   )                                
     In the above equation, G(s) is the transfer function of the PID controller, e(s) is the controller input, u(s) is the controller output, k 1  is the coefficient for proportional feedback, k 2  is the coefficient for integral feedback, and k 3  is the coefficient for derivative feedback. The coefficients k i  may be represented by a coefficient vector K, where K=[k 1 , k 2 , k 3 ]. The vector K is commonly called a Coefficient Gain Schedule (CGS). The values of the coefficients K used in the linear PID control system are based on a dynamic model of the plant. When the plant is unstable, nonlinear, and/or time-variant, then the coefficients in K are often controlled by an AI control system. 
     FIG. 1 shows a typical prior art AI control system  100 . An input y(t) of the control system  100  is provided to a plus input of an adder  104  and an output x(t) of a plant  110  is provided to a minus input of the adder  104 . An output of the adder  104  is provided as an error signal e(t) to an error signal input of a PID controller  106 . An output u(t) of the PID controller  106  is provided to a first input of an adder  108  and to a first input of a Genetic Algorithm (GA). A disturbance m(t) is provided to a second input of the adder  108 . An output u*(t) of the adder  108  is provided to an input of the plant  110 . The plant  110  has a transfer function H(s) and an output x(t), where x(t)⇄X(s) (where the symbol⇄denotes the Laplace transform) and X(s)=G(s)u*(s). An output of the genetic algorithm  116  is provided to an input of a Fuzzy logic Neural Network (FNN)  118  and an output of the fuzzy neural network  118  is provided to a Fuzzy Controller (FC)  120 . An output of the fuzzy controller  120  is a set of coefficients K, which are provided to a coefficient input of the PID controller  106 . 
     The combination of the genetic algorithm  216  and the entropy calculator  214  comprises a Simulation System of Control Quality  215 . The combination of the fuzzy neural network  218  and the fuzzy controller  220  comprises a Fuzzy Logic Classifier System FLCS  219 . The combination of the plant  210  and the adder  208  comprises a disturbed plant model  213 . The disturbed plant signal u*(t)=u(t)+m(t), and the disturbance m(t) are typically unobservable. 
     The error signal e(t) provided to the PID controller  106  is the difference between the desired plant output value y(t) and the actual plant output value x(t). The PID controller  106  is designed to minimize the error represented by e(t) (the error being the difference between the desired and actual output signal signals). The PID controller  106  minimizes the error e(t) by generating an output signal u(t) which will move the output signal x(t) from the plant  110  closer to the desired value. The genetic algorithm  116 , fuzzy neural network  118 , and fuzzy controller  120  monitor the error signal e(t) and modify the gain schedule K of the PID controller  106  in order to improve the operation of the PID controller  106 . 
     The PID controller  106  constitutes a reverse model relative to the plant  110 . The genetic algorithm  116  evolves an output signal α based on a performance function ƒ. Plural candidates for α are produced and these candidates are paired according to which plural chromosomes (parents) are produced. The chromosomes are evaluated and sorted from best to worst by using the performance function ƒ. After the evaluation for all parent chromosomes, good offspring chromosomes are selected from among the plural parent chromosomes, and some offspring chromosomes are randomly selected. The selected chromosomes are crossed so as to produce the parent chromosomes for the next generation. Mutation may also be provided. The second-generation parent chromosomes are also evaluated (sorted) and go through the same evolutionary process to produce the next-generation (i.e., third-generation) chromosomes. This evolutionary process is continued until it reaches a predetermined generation or the evaluation function ƒ finds a chromosome with a certain value. The outputs of the genetic algorithm are the chromosomes of the last generation. These chromosomes become input information a provided to the fuzzy neural network  118 . 
     In the fuzzy neural network  118 , a fuzzy rule to be used in the fuzzy controller  120  is selected from a set of rules. The selected rule is determined based on the input information α from the genetic algorithm  116 . Using the selected rule, the fuzzy controller  120  generates a gain schedule K for the PID controller  106 . The vector coefficient gain schedule K is provided to the PID controller  106  and thus adjusts the operation of the PID controller  106  so that the PID controller  106  is better able to minimize the error signal e(t). 
     Although the AI controller  100  is advantageous for accurate control in regions near linearized equilibrium points, the accuracy deteriorates in regions away from the linearized equilibrium points. Moreover, the AI controller  100  is typically slow or even unable to catch up with changes in the environment surrounding the plant  110 . The PID controller  106  has a linear transfer function G(s) and thus is based upon a linearized equation of motion for the plant  110 . Since the evaluation function ƒ used in the genetic algorithm  116  is only based on the information related to the input e(t) of the linear PID controller  106 , the controller  100  does not solve the problem of poor controllability typically seen in linearization models. Furthermore, the output results, both in the gain schedule K and the output x(t) often fluctuate greatly, depending on the nature of the performance function ƒ used in the genetic algorithm  116 . The genetic algorithm  116  is a nonlinear optimizer that optimizes the performance function ƒ. As is the case with most optimizers, the success or failure of the optimization often ultimately depends on the selection of the performance function ƒ. 
     The present invention solves these and other problems by providing a new AI control system  200  shown in FIG.  2 . Unlike prior AI control systems, the control system  200  is self-organizing and uses a new performance function ƒ which is based on the physical law of minimum entropy. An input y(t) of the control system  200  is provided to a plus input of an adder  204  and an output x(t) of a plant  210  is provided to a minus input of the adder  204 . An output of the adder  204  is provided as an error signal e(t) to an error signal input of a PID controller  206  and to an input of a fuzzy controller  220 . An output u(t) of the PID controller  206  is provided to a first input of an adder  208  and to a first input of an entropy calculator (EC)  214 . A disturbance m(t) is provided to a second input of the adder  208 . An output u*(t) of the adder  208  is provided to an input of the plant  210 . The plant  210  has a transfer function H(s) and an output x(t), such that X(s)=G(s)u*(s), where x(t)⇄X(s). The output x(t) is provided to a second input of the entropy calculator  214  and to the minus input of the adder  204 . An output of the entropy calculator  214  is provided to an input of a genetic algorithm  216  and an output of the genetic algorithm  216  is provided to an input of a Fuzzy logic Neural Network (FNN)  218 . An output of the fuzzy neural network  218  is provided to a rules selector input  222  of the fuzzy controller  220 . A Coefficient Gain Schedule (CGS) output  212  of the fuzzy controller  222  is provided to a gain schedule input of the PID  206 . 
     The error signal e(t) provided to the PID controller  206  is the difference between the desired plant output value y(t) and the actual plant output value x(t). The PID controller  206  is designed to minimize the error represented by e(t). The PID controller  206  minimizes the error e(t) by generating an output signal u(t) which will move the output signal x(t) from the plant  210  closer to the desired value. The fuzzy controller  220  monitors the error signal e(t) and modifies the gain schedule K of the PID controller  206  according to a fuzzy control rule selected by the fuzzy neural network  218 . 
     The genetic algorithm  216  provides a teaching signal K T  to the fuzzy neural network  218 . The teaching signal K T  is a global optimum solution of a coefficient gain schedule K generated by the genetic algorithm  216 . 
     The PID controller  206  constitutes a reverse model relative to the plant  210 . The genetic algorithm  216  evolves an output signal a based on a performance function ƒ. Plural candidates for α are produced and these candidates are paired by which plural chromosomes (parents) are produced. The chromosomes are evaluated and sorted from best to worst by using the performance function ƒ. After the evaluation for all parent chromosomes, good offspring chromosomes are selected from among the plural parent chromosomes, and some offspring chromosomes are randomly selected. The selected chromosomes are crossed so as to produce the parent chromosomes for the next generation. Mutation is also employed. The second-generation parent chromosomes are also evaluated (sorted) and go through the same evolutionary process to produce the next-generation (i.e., third-generation) chromosomes. This evolutionary process is continued until it reaches a predetermined generation or the evaluation function ƒ finds a chromosome with a certain value. Then, a component from a chromosome of the last generation becomes a last output, i.e., input information α provided to the fuzzy neural network  218 . 
     In the fuzzy neural network  218 , a fuzzy rule to be used in the fuzzy controller  220  is selected from a set of rules. The selected rule is determined based on the input information α from the genetic algorithm  216 . Using the selected rule, the fury controller  220  generates a gain schedule K for the PID controller  206 . This is provided to the PID controller  206  and thus adjusts the operation of the PID controller  206  so that the PID controller  206  is better able to minimize the error signal e(t). 
     The fitness function ƒ for the genetic algorithm is given by              f   =     min                        S          t                 (   2   )                                
     were                     S          t       =     (              S   c            t       -            S   u            t         )             (   3   )                                
     The quantity dS u /dt represents the rate of entropy production in the output x(t) of the plant  210 . The quantity dS c /dt represents the rate of entropy production in the output u(t) of the PID controller  206 . 
     Dynamic stability properties of the plant  210  near an equilibrium point can be determined by use of Lyapunov functions. Let V(x) be a continuously differentiable scalar function defined in a domain D⊂R n  that contains the origin. The function V(x) is said to be positive definite if V( 0 )=0 and V(x)&gt;0 for x≠0. The function V(x) is said to be positive semidefinite if V(x)≧0 for all x. A function V(x) is said to be negative definite or negative semidefinite if −V(x) is positive definite or positive semidefinite, respectively. The derivative of V along the trajectories {dot over (x)}=ƒ(x) is given by:                V        (   x   )       =         ∑     i   =   1     n                         ∂   V       ∂     x   i                x   .     i         =         ∂   V       ∂   x            f        (   x   )                   (   4   )                                
     where ∂V/∂x is a row vector whose ith component is ∂V/∂x i  and the components of the n-dimensional vector ƒ(x) are locally Lipschitz functions of x, defined for all x in the domain D. The Lyapunov stability theorem states that the origin is stable if there is a continuously differentiable positive definite function V(x) so that V(x) is negative definite. A function V(x) satisfying the conditions for stability is called a Lyapunov function. 
     Calculation of the Lyapunov dynamic stability and entropy production for a closed nonlinear mechanical system is demonstrated by using the Holmes-Rand (Duffing-Van der Pol) nonlinear oscillator as an example. The Holmes-Rand oscillator is described by the equation: 
     
       
           {umlaut over (x)}+ (α+β x   2 ) {dot over (x)}−γx+x   3 =0  (5) 
       
     
     where α, β, and γ are constant parameters. A Lyapunov function for the Holmes-Rand oscillator is given by:                V   =         1   2            x   .     2       +     U        (   x   )           ,                  where                                U     =         1   4          x   4       -       1   2        γ                   x   2                   (   6   )                                
     Entropy production d i S/dt for the Holmes-Rand oscillator is given by the equation:                       i        S          t       =       (     α   +     β                   x   2         )            x   .     2               (   7   )                                
     Equation 5 can be rewritten as:                  x   ¨     +       (     α   +     β                   x   2         )          x   .       +       ∂   U       ∂   x         =   0           (   8   )                                
     After multiplying both sides of the above equation by {dot over (x)}, then dV/dt can be calculated as:                     V          t       =           x   ¨          x   .       +         ∂   U       ∂   x            x   .         =       -     1   T                   i        S          t                   (   9   )                                
     where T is a normalizing factor. 
     An interrelation between a Lyapunov function and the entropy production in an open dynamic system can be established by assuming a Lyapunov function of the form              V   =       1   2            ∑     i   =   1     6                     (       q   i   2     +     S   2       )                 (   10   )                                
     where S=S u −S c  and q i =(α, γ, β, {dot over (α)},{dot over (γ)},{dot over (β)}). It is possible to introduce the entropy function S in the Lyapunov function V because entropy S is also a scalar function of time. Differentiation of V with respect to time gives:                     V          t       =         ∑     i   =   1     6                       q   i            q   .     i         +     S        S   .                 (   11   )                                
     In this case, q i =ψ i (q i ,τ, t), S=S u −S c , {dot over (S)}={dot over (S)} u −{dot over (S)} c  and thus:                     V          t       =         ∑     i   =   1     6                       q   i            ϕ   i          (       q   i     ,   τ   ,   t     )           +       (       S   u     -     S   c       )          (              S   u            t       -            S   c            t         )                 (   12   )                                
     A special case occurs when β=0 and the Holmes-Rand oscillator reduces to a force-free Duffing oscillator, wherein:                       i        S          t       =       -   α                       x   .     2                     (     Duffing                 oscillator     )               (   13   )                                
     A Van der Pol oscillator is described by the equation: 
     
       
           {umlaut over (x)}+ ( x   2 −1) {dot over (x)}+ =0  (14) 
       
     
     and the entropy production is given by:                       i        S          t       =       1   T          (       x   2     -   1     )            x   .     2                     (     Van                 der                 Pol                 oscillator     )               (   15   )                                
     For a micro-mobile robot in fluid, a mechanical model is given by:                    m   1            x   ¨     1       +       C   d          ρ   2          A   1                 x   .     1                 x   .     1       +       K   1          (       x   1     -     x   0     -       l   1          θ   0         )       -       K   2          (       x   2     -     x   1     -       l   2          θ   1         )         =   0           (   16   )                     m   2            x   ¨     2       +       C   d          ρ   2          A   2                 x   .     2                 x   .     2       +       K   2          (       x   2     -     x   1     -       l   2          θ   1         )       -       K   3          (       x   3     -     x   2     -       l   3          θ   2         )         =   0           (   17   )                     m   3            x   ¨     3       +       C   d          ρ   2          A   3                 x   .     3                 x   .     3       +       K   3          (       x   3     -     x   2     -       l   3          θ   2         )         =   0           (   18   )                                
     where:                θ     n   +   1       =         -     1   2            θ   n       +       3   2          1     l     n   +   1              (       x     n   +   1       -     x   n       )                 (   19   )                                
     Values for a particular micro-mobile robot are given in Table 1 below. 
     
       
         
               
               
               
             
           
               
                 TABLE 1 
               
               
                   
               
               
                 Item 
                 Value 
                 Units 
               
               
                   
               
             
             
               
                 m 1   
                 1.6 × 10 −7   
                 kg 
               
               
                 m 2   
                 1.4 × 10 −6   
                 kg 
               
               
                 m 3   
                 2.4 × 10 −6   
                 kg 
               
               
                 l 1   
                 2.0 × 10 −3   
                 m 
               
               
                 l 1   
                 4.0 × 10 −3   
                 m 
               
               
                 l 3   
                 4.0 × 10 −3   
                 m 
               
               
                 K 1   
                 61.1 
                 N/m 
               
               
                 K 2   
                 13.7 
                 N/m 
               
               
                 K 3   
                 23.5 
                 N/m 
               
               
                 A 1   
                 4.0 × 10 −6   
                 m 2   
               
               
                 A 2   
                 2.4 × 10 −6   
                 m 2   
               
               
                 A 3   
                 4.0 × 10 −6   
                 m 2   
               
               
                 C d   
                 1.12 
                 — 
               
               
                 ρ 
                 1000 
                 Kg/m 3   
               
               
                   
               
             
          
         
       
     
     Entropy production for the micro-mobile robot is given by the equation:                       S   i            t       =       ∑     n   =   1     3                       C     d                       ρ   2          A   n                 x   .     n                 x   .     n   2                 (   20   )                                
     and the Lyapunov Function is given by:              V   =         ∑     i   =   1     3                       m   i            x   .     i   2                     ρ   2         +       ∑     i   =   1     3                           K   i          (       x   i     -     x     i   -   1       -       l   i          θ     i   -   1           )       2     2       +       S   2     2               (   21   )                                
     where S=S i −S c  and S c  is the entropy of a controller with torque τ. 
     The necessary and sufficient conditions for Lyapunov stability of a plant is given by the relationship:                      ∑     i                                       q   i            ϕ   i          (       q   i     ,   τ   ,   t     )         &lt;       (       S   u     -     S   c       )          (              S   c            t       -            S   i            t         )         ,              S   c            t       &gt;              S   i            t          u               (   22   )                                
     According to the above equation, stability of a plant can be achieved with “negentropy” S c  (by Brillouin&#39;s terminology) where a change of negentropy dS c /dt in the control system  206  is subtracted from a change of entropy dSi/dt in the motion of the plant  210 . 
     The robust AI control system  200  provides improved control of mechanical systems in stochastic environments (e.g., active vibration control), intelligent robotics and electro-mechanical systems (e.g., mobile robot navigation, manipulators, collective mobile robot control), bio-mechanical systems (e.g., power assist systems, control of artificial replaced organs in medical systems as artificial lung ventilation), micro electromechanical systems (e.g., micro robots in fluids), etc. 
     The genetic algorithm realizes the search of optimal controllers with a simple structure using the principle of minimum entropy production. The fuzzy neural network controller offers a more flexible structure of controllers with a smaller torque, and the learning process produces less entropy. The fuzzy neural network controller gives a more flexible structure to controllers with smaller torque and the learning process produces less entropy than a genetic analyzer alone. Thus, an instinct mechanism produces less entropy than an intuition mechanism. However, necessary time for achieving an optimal control with learning process on fuzzy neural network (instinct) is larger than with the global search on genetic algorithm (intuition). 
     Realization of coordinated action between the look-up tables of the fuzzy controller  220  is accomplished by the genetic algorithm and the fuzzy neural network. In particular, the structure  200  provides a multimode fuzzy controller coupled with a linear or nonlinear neural network  218 . The control system  200  is a realization of a self-organizing AI control system with intuition and instinct. In the adaptive controller  200 , the feedback gains of the PID controller  210  are changed according to the quantum fuzzy logic, and approximate reasoning is provided by the use of nonlinear dynamic motion equations. 
     The fuzzy tuning rules for the gains k i  are shaped by the learning system in the fuzzy neural network  218  with acceleration of fuzzy rules on the basis of global inputs provided by the genetic algorithm  216 . The control system  200  is thus a hierarchical, two-level control system that is intelligent “in small.” The lower (execution) level is provided by a traditional PID controller  206 , and the upper (coordination) level is provided by a KB (with fuzzy inference module in the form of production rules with different model of fuzzy implication) and fuzzification and de-fuzzification components, respectively. 
     Parallel soft computing for fuzzy control is preferably based on quantum computing. The genetic algorithm  216  simulates an intuition mechanism of choosing the optimal structure of the PID controller  206  by using the fitness function, which is the measure of the entropy production, and the evolution function, which in this case is entropy. 
     Unicycle Robot with Self-Organizing Control 
     One embodiment of the self-organizing control system  200  is advantageously used in the simulation of a robot unicycle  300  shown in FIGS. 3 and 4. The unicycle  300  comprises a motor  304  driving a weighted yoke  302 . The motor  304  and a rate gyro sensor  306  are attached to a frame  318 . A right closing link mechanism  308  is attached between the frame  318  and an axle of a wheel  312 . A left closing link mechanism  316  is attached between the frame  318  and the axle of the wheel  312 . A motor  310  drives the right link mechanism  308  and a motor  314  drives the left link mechanism  316 . 
     A unicycle is an inherently unstable nonlinear system and simultaneous control of both longitudinal and lateral stability is needed to maintain the unicycle&#39;s postural stability. It is an unstable problem in three dimensions. However, a rider can achieve a postural stability on a unicycle, keep the wheel speed constant and change the unicycle&#39;s posture in the yaw direction by using a flexible body, good sensory systems, skill and intelligence. Thus, the unicycle is a reasonable testbed for the control system  200 . 
     Analysis of the control system  200  is based on a thermodynamic model representation of the dynamic equation of motion for the plant, described as an open dynamic system. Accordingly, a mechanical model with dissipative processes for the plant in  210  for a unicycle is given by the equation:                [           q   ¨             λ         ]     =         [           M        (   q   )               -     ∂   c         ∂   q                 E        (   q   )           0         ]       -   1                 [     τ   -       B        (   q   )            [       q   .     ,     q   .       ]       -       C        (   q   )            [       q   .     2     ]       -       D        (   q   )            [     q   .     ]       -     G        (   q   )       -     F        (     q   ,     q   .       )         ]                 (   23   )                                
     A thermodynamic description of the production of entropy in the control system  200  used in a unicycle can be expressed as:                [                  S   u            t                      Sc          t                        ]     =         [           M        (   q   )           0           1       0         ]       -   1                   [       τ   d     -       B        (   q   )            [       q   .     ,     q   .       ]       -       C        (   q   )            [       q   .     2     ]       -       D        (   q   )            [     q   .     ]       -     G        (   q   )       -     f        (     q   ,     q   .       )         ]          [           q   .             0         ]            
                 (   24   )                                
     where τ=(τ 105 ,0,0,0,0, τ θ2 ,0,τ θ4 ,τ n ) and dS c /dt is the rate of entropy production in the PID controller  206  and dSu/dt is the rate of entropy production of the plant  210  (the unicycle). E(q) is a 4×4 coefficient matrix for acceleration. F(q,{dot over (q)}) is a 4 dimensional vector containing coriolis and centrifugal terms. Both E and F are determined from the constraint equations on a closed-link loop acceleration. The matrix ∂c/∂q is a 4×4 matrix determined by constraint equations of the two closed link loops. The quantity λ is a 4-dimensional vector of Lagrangian multipliers; M(q) is a 9×9 matrix of masses; B(q) is a 9×36 matrix of coriolis coefficients; [{dot over (q)},{dot over (q)}] is a 36×1 vector of velocity products given by ({dot over (q)}, {dot over (q)})=[{dot over (Ψ)},{dot over (α)}, {dot over (Ψ)}{dot over (γ)}, . . . {dot over (θ)} 4 {dot over (η)}] T ; C(q) is a 9×9 matrix of centrifugal coefficients; [{dot over (q)} 2 ] is a 9×1 vector given by ({dot over (Ψ)} 2 ,{dot over (α)} 2 ,{dot over (γ)} 2 , {dot over (β)} 2 ,{dot over (θ)} 1   2 ,{dot over (θ)} 2   2 ,{dot over (θ)} 3   2 ,{dot over (θ)} 4   2 , {dot over (η)} 2 ) T ; D(q) is a 9×9 matrix of friction coefficients; [{dot over (q)}]=({dot over (Ψ)},{dot over (α)},{dot over (γ)}, {dot over (β)},{dot over (θ)} 1 ,{dot over (θ)} 2 ,{dot over (θ)} 3 ,{dot over (θ)} 4 , {dot over (η)}) T ; and G(q) is a 9×1 vector of gravity terms. The state vector q is defined as the vector [q]=(Ψ,α,γ,β,θ 1 ,θ 2 , θ 3 ,θ 4 ,η) T . The vector τ d  describes the dissipative parts of the control torque on the unicycle plant. 
     Stability analysis of the unicycle as a nonlinear plant  210  is accomplished, as before, by using asymptotic methods with Lyapunov functions and methods of qualitative physics based on correlation between Lyapunov and production entropy functions. The new approach for definition of the Lyapunov function described above is used. A Lyapunov function for the system is given by:              V   =       1   2          (         ∑     i   =   1     n                     q   i   2       +     S   2       )               (   25   )                                
     where S=S u −S c  and q i  is a generalized coordinate of the plant model. 
     The input u* is a random signal produced by the sum of the control signal u and the disturbance (noise) m(t). Calculation of entropy production performed in the entropy calculator  214 . The output of the entropy calculator is the full entropy production of the dynamic system, dS/dt, and the minimum of this function is used as a fitness function in the genetic algorithm  216 . The output of the genetic algorithm  216  is an optimal gain structure K T  for the PID controller  206 . The structure K T  is used as the teaching signal for the fuzzy neural network  218  in the FLCS  219 . 
     The FLCS  219  comprises the fuzzy neural network  218  and the fuzzy controller  220 . The fuzzy neural network  218  generates a Knowledge Base (KB) comprising a look-up table of production rules by using a neural network learning process. The KB output of the fuzzy neural network  218  comprises a Look-up table of Fuzzy Rules (LPTR) as input for the fuzzy controller  220 . The KB provides for self-organization of the fuzzy controller  220  in the FLCS  219 . The FLCS  219  provides the adaptation process of an intelligent fuzzy controller  220  structure to the random disturbance m(t) in the control signal u*(t). The Output Coefficient Gain Schedule (CGS) of the fuzzy controller is provided as an input to the PID-controller to program the coefficient gains k i  in accordance. This approach differs from procedural design of fuzzy hybrid PID-type controllers. This design method provides new concepts of a robust structure controller  210  with optimal learning, with minimum entropy production in the control process, and self-organization of the KB structure of the fuzzy controller  220 . 
     The genetic algorithm  216  is preferably a massy genetic algorithm with elite strategy. Genetic operators have the following parameters. For the selection, a roulette method (reduce rate=5) is used. A reproduction operator replaces some chromosome generations. Crossover (with a probability of 0.64) and mutation (with a probability of 0.08) determines, by random numbers, every chromosome&#39;s bit (i.e., whether crossover and/or mutation occur). Genes encode each parameter from the interval [ 0 ,10000] and the number of used chromosomes is 100. 
     The fuzzy neural network  218  is preferably a feedforward structure network with four layers. The structure of the fuzzy neural network  218  is preferably similar to the structure of the fuzzy controller  220 . The first layer in the fuzzy neural network  218  is comparable to fuzzification in the fuzzy controller  220 . The second layer in the fuzzy neural network  218  is comparable to an “If . . . ” part in the inference mechanism of the fuzzy controller  220 . The third layer in the fuzzy neural network is comparable to a “Then . . . ” part in the inference mechanism of the fuzzy controller  220 . Finally, the fourth layer in the fuzzy neural network is comparable to a de-fuzzification operator structure in the fuzzy controller  220 . The structure of the fuzzy neural network  218  may use a max-min method for the fuzzy inference and the fuzzy production rules. The membership functions of the fuzzy neural network  218  are preferably optimized by the back-propagation method with minimum information entropy as the criteria of accuracy. In one embodiment used with a unicycle simulation, the fuzzy neural network  218  has 2 inputs, 2 hidden layers (14 units in the first layer and 49 units in the second layer), 1 output, 7 linguistic terms for error and change of error, and 49 rules in the Knowledge Base. 
     One aspect of the simulation is a comparison of results between the genetic algorithm  216  and the fuzzy neural network  218  for the torque and PID controller gains with minimum entropy production. 
     FIG. 5 is a graph  500  of angle (in radians) as a function of time (in seconds) for a simulated unicycle with an AI controller  200 . The graph  500  shows a curve  502  showing the angle of the rotor, a curve  504  showing the angle of the wheel, a curve  506  showing yaw of the unicycle, a curve  508  showing roll of the unicycle, and a curve  510  showing pitch of the unicycle. 
     FIG. 6 is a graph  600  of angular velocity (in radians/second) as a function of time (in seconds). The graph  600  shows a curve  602  showing the angle of the rotor, a curve  604  showing the angle of the wheel, a curve  606  showing yaw of the unicycle, a curve  608  showing roll of the unicycle, and a curve  610  showing pitch of the unicycle. 
     FIG. 7 is a graph showing a phase surface for the roll angle (in radians) and angular velocity (in radians/second) as a result of the simulation for the posture control performed by the robot unicycle in FIG.  3 . 
     FIG. 8 is a graph showing a phase surface for the pitch angle (in radians) and angle speed (in radians/second) as a result of the simulation for the posture control performed by the robot unicycle in FIG.  3 . 
     FIG. 9 is a graph  900  showing the time changes in angle (in radians) as a function of time (in seconds) for each part in the unicycle as a result of the simulation for the posture control using the methods in the prior art. The graph  900  includes a curve  902  showing the angle of the rotor, a curve  904  showing the angle of the wheel, a curve  906  showing yaw of the unicycle, a curve  908  showing roll of the unicycle, and a curve  910  showing pitch of the unicycle. 
     FIG. 10 is a graph showing the time changes in angular velocity (in radians/second) versus time (in seconds) for each part of the unicycle as a result of the simulation for the posture control in the prior art. The graph  1000  includes a curve  1002  showing the angle of the rotor, a curve  1004  showing the angle of the wheel, a curve  1006  showing yaw of the unicycle, a curve  1008  showing roll of the unicycle, and a curve  1010  showing pitch of the unicycle. 
     FIG. 11 is a graph showing a phase surface for the roll angle and angle speed as a result of the simulation for the posture control in the prior art which is performed by the robot unicycle in FIG.  3 . 
     FIG. 12 is a graph  1200  illustrating a phase surface for the pitch angle and angle speed as a result of a simulation of the robot unicycle in FIG. 3, using the posture control in the prior art. 
     FIG. 13 is a graph  1300  illustrating another example of the simulation for the posture control indicating the effect of the control method in the present invention. The graph  1300  is a plot of angle (in radians) versus time (in seconds). The graph  1300  includes a curve  1302  showing roll without control, a curve  1304  showing roll with control, a curve  1306  showing pitch without control, and a curve  1308  showing pitch with control. 
     Although the foregoing has been a description and illustration of specific embodiments of the invention, various modifications and changes can be made thereto by persons skilled in the art, without departing from the scope and spirit of the invention as defined by the following claims.