Abstract:
A method automatically determines, controller parameters of a plurality of controllers for controlling a process controlled system that is composed of a plurality of coupled control systems. The controller parameters of the respective controllers are calculated from the measured responses or the individual control systems, being calculated by a computer for equivalent transfer functions that consider couplings between control systems in the, process controlled system. This occurs in an iterative method, in which a multi-dimensional, non-linear equation system is supplied step-by-step to an optimized solution A, a solution is individually calculated for each non-linear equation, that enters into the other equations. A different standard method for solving non-linear equation systems can also be utilized. A primary application is envisioned for: setting a furnace having a plurality of heating zones.

Description:
BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     The invention is directed to a method for setting a plurality of controllers for controlling a process controlled system composed of coupled control systems. 
     2. Description of the Related Art 
     In commonly used control systems without decoupling filters, as can often be encountered in practice, independent control systems are assumed for simplification for using manual initialization. Existing couplings between control systems can usually not be clearly recognized and are therefore frequently ignored. When, however, such couplings between control systems exist in application, then the initialization or readjustment of a control system affects other control systems in the overall control system. 
     In this situation the manual setting is accomplished in a plurality of steps, in which parameters of a controller are modified in order, to quickly compensate a rated value different when setting a new rated value, among other things moreover when, an overshooting can occurs in the control circuits, the controller parameters must be readjusted in the coupled control systems, so that the control circuits work satisfactorily. Each-therefor adjustment of the controller parameters in a control system clearly effects a modification of the control circuits in the control systems that have a coupling (of different extents) to the adjusted control system. When a process controlled system contains a plurality of control systems, the initialization is a lengthy and tiresome one must wait for the stationary condition of the overall process control system before every modification and control circuits that are already set are readjusted with progressive initialization due to the modifications of other controllers. 
     SUMMARY OF THE INVENTION 
     In a processed controlled system composed of a plurality of control systems, the object of the invention is to undertake a setting of the in which therefore, only measured values for the transient response of the systems are employed in order to be able to set the parameters of the controllers all at once. 
     This object is achieved a method for setting a plurality of controllers for controlling a process controlled system having coupled control systems, comprising the steps of determining a transfer function from a type of a respective controller, implementing, for each control circuit, measurements of setting signals and output signals of the process controlled system, determining, from these measurements, transfer functions for the coupled control systems contained in the process controlled system and their couplings, calculating, for each controller, an equivalent path that utilizes the couplings with other control systems in the process controlled system, and iteratively determining, for each respective coupled control system, controller parameters utilizing couplings with the other control systems for each respective controller until a predetermined quality has been reached. 
     With the inventive method, the respective controller parameters are calculated for each controller in a process controlled system composed of a plurality of control systems, in which the measured values for the transient response of the process controlled system are necessary in order to automatically determine all of the control parameters. The initialization of, for example, a furnace that is composed of a plurality of heating zones is thus accomplished in one step, since the manual iteration is eliminated because this initialization is automatically implemented on a computer. This implementation is possible because the couplings are also taken into consideration in the inventive method and, in contrast to the manual setting method, the controller parameters are modified based on the behavior of multiple control system. A considerable time saving is achieved because the required iterations for optimizing the controllers are implemented in a computer. Transient response events of temperature control systems can be simulated significantly faster on computer. 
     A further advantage of the invention is that controller settings can be reproduced and are thus less dependent on the experience level of the person undertaking the setting. 
     In general, automatic optimization is more economical than the described, manual iteration method, since resources are not wasted. 
     The automatically determined controller settings lie in the immediate proximity of a possible optimum setting, which enables an economical operation of the system. 
     Further, the automatic design method supplies a “reference” that can be employed for judging an optimization potential of a complex system. 
     It is an additional advantage of the invention that the noise behavior is also taken into consideration in the automatic controller setting. This contrast with the manual controller settings which are frequently only undertaken for the management behavior. 
     Further, it is also possible to reduce the plurality of required iterations on the computer by utilizing modern numerical methods. 
     A transfer function is determined for each controller A PID controller has the typical controller parameters “proportional part P”, “integrating part I” and “differentiating part D” available to it. In control technology, the corresponding transfer function for different types of controllers are adequately known to a person skilled in the art, usually in frequency-dependent presentation with Laplace operators s=jω. 
     As mentioned, measurements are implemented at the process controlled system, in which a transient behavior is registered for each control system to be taken into consideration. The transient behavior of each control system in the process controlled system is thus available and one knows the type of controller. An equivalent system is now calculated that takes couplings of the respective control system with other control systems into consideration. By knowing of this control system, it is possible to calculate the controller parameters for each controller of a control circuit in iterative steps, in which a separate presentation of the controller parameters dependent on other controller parameters exists for each control system and, as a result of the computer, the proper controller parameters for each control circuit are respectively identified with the controller parameters of the coupled control systems that have already been calculated. These parameters enter into the other equations that comprise a coupling to the control system. In a next step, the controller parameters of the next controller are optimized, etc. 
     One development of the invention is comprised in calculating the controller parameters for the respective controller of a control system in that a zero point search of a presentation present in implicit form that contains the transfer functions of all controllers, systems and couplings is implemented, whereby an iterative Newton method can also be applied. 
     Another development of the invention is comprised in implementing the setting of the plurality of controllers of the process controlled system with non-linear controller parameterization methods, preferably with a neural network (neural controller parameterization method). 
     One application of the inventive method is in an initialization of a furnace composed of a plurality of heating zones, e.g., for burning bricks, in which fixed demands are made of the temperature curve within the furnace. 
     Other applications of the invention are for combined pressure/flow controls in supply lines or temperature/humidity controls in air conditioning systems. 
     One development of the invention is comprised in the use of a processed control system that is composed of at least two control systems and is controlled with at least two controllers. 
     Developments of the invention are discussed below. 
    
    
     Exemplary embodiments of the invention are presented in greater detail on the basis of the following figures. 
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is a block diagram that shows a process controlled system having two manipulated quantities and two regulating variables; 
     FIG. 2 is a block diagram that shows a process controlled system having two controllers and without a decoupling filter; 
     FIG. 3 is a block diagram that illustrates the situation before initialization of the second control circuit; 
     FIG. 4 is a block diagram that shows a general dependency of a control system on a plurality of coupled control systems; 
     FIG. 5 is a block flow diagram that illustrates steps of a method for controller setting. 
    
    
     The invention is explained below on the basis of a process controlled system having 2 control systems, i.e., with 2 inputs and two outputs. The expansion for the case of more than 2 inputs and outputs is described afterwards. 
     DESCRIPTION OF THE PREFERRED EMBODIMENT 
     Description and Identification of a Multi-Quantity System 
     FIG. 1 shows a process controlled system RS having two manipulated variables u 1  and u 2  and two regulating variables y 1  and Y 2 . 
     For an exemplary water faucet having a hot water wheel and a cold water wheel, the quantity u 1  and u 2  respectively characterizes the position of the hot water or, respectively, of the cold water wheel. The regulating variables y 1  and y 2  respectively denote the water temperature and the volume stream of the water (water flow). 
     Possible applications of such a system include control circuits for pressure and temperature of a boiler, a steam turbine, a furnace with a plurality of heating zones, (e.g., for burning bricks), air temperature and air humidity for air conditioning controls. 
     The following is valid in the frequency domain given a linear, time-invariant multiple-quantity system with two manipulated quantities and two regulating variables:                  [             Y   1          (   s   )                   Y   2          (   s   )             ]     =       [             H   11          (   s   )               H   12          (   s   )                   H   21          (   s   )               H   22          (   s   )             ]     ·     [             U   1          (   s   )                   U   2          (   s   )             ]         ,           (   1   )                                
     where 
     s denotes the Laplace variable (s=jω) 
     Y 1 (s), u i (s) denotes the Laplace transform of y i  and u i  (i=1,2); 
     H 12 , H 22  sub-transfer functions of the diagonals (main path, main branches); 
     H 12 , H 21  denotes sub-transfer functions of the subsidiary diagonals (couplings). 
     The sub-transfer functions are obtained by evaluating equation (1) for U 1 (s) or, respectively, for U 2 (s), i.e.,            H   11          (   s   )       =           Y   1          (   s   )           U   1          (   s   )                         U   2          (   s   )       ≡   0            
              H   12          (   s   )       =           Y   1          (   s   )           U   2          (   s   )                         U   1          (   s   )       ≡   0            
              H   21          (   s   )       =           Y   2          (   s   )           U   1          (   s   )                         U   2          (   s   )       ≡   0                         and        
            H   22          (   s   )         =           Y   2          (   s   )           U   2          (   s   )                         U   1          (   s   )       ≡   0                                                               
     A measurement of the system frequency responses for s=jω in a given application occurs in that the manipulated variable u 1  is adjusted first and the manipulated variable u 2  is retained. Both regulating variables y 1  and y 2  are registered. Subsequently, the manipulated variable u 1  is retained and the manipulated variable u 2  is adjusted. Regulating variables y 1  and y 2  are again registered. Overall, accordingly, four reply signals to four input signals for the four sub-transfer functions H 11 (s), H 12 (s), H 21 (s) and H 22 (s) in the indicated sequence. 
     DECENTRALIZED DESIGN WITHOUT DECOUPLING FILTER 
     FIG. 2 shows the structure of a process controlled system RS without a decoupling filter, having a controller one REG 1 , a controller two REG 2  reference quantities w 1  and w 2 , and two control circuits RK 1  and RK 2 . This control structure is frequently employed in practice since it is often assumed in the initialization that two independent process controlled systems, i.e., the two main paths having the transfer functions H 11 (s) and H 22 (s) are present. Existing couplings with the transfer functions H 12 (s) and H 21 (s) can often not be clearly recognized or are often simply left out of consideration. In many instances, this control structure is selected because the design, the initialization, and the maintenance of preceding decoupling filters is complicated. This, for example, is the case when couplings change frequently due to remodelings, utilization of different operating means (catalysts) or fluctuating product quantities, or when an installation of decoupling filters is not technically possible in an existing automated process. This is the case when economical discrete controllers are utilized in a system for cost reasons or when the process quantities affected lie in system parts that are controlled with separate systems. 
     The problems with the initialization of a process controlled system without decoupling filters are discussed in greater detail below. It is assumed for this purpose that the control circuit RK 1  for the first main path has already been closed, and the initialization of the second control circuit is to take place (also see FIG.  3 ). The transfer function of the equivalent path H 20 (s) is important for the initialization of the second control circuit RK 2 , whereby            H   20          (   s   )       =         Y   2          (   s   )           U   2          (   s   )                                
     the transfer is calculated below. The argument “(s)” is omitted below for the sake of clarity. 
     The transfer functions of the two controllers REG 1  and REG 2  are referenced below as H R1  and H R2 . The reference quantity w 1  is to be set to zero for the calculation. Following from equation (1) is                  Y   2       U   2       =             H   21          U   1       +       H   22          U   2           U   2       =       H   22     +       H   21            U   1       U   2                     (1a)                                
     The quantity U 1  is to be eliminated from this equation (1 a ), i.e., to be expressed by U 2 . Initially following is: 
     
       
           U   1   =−H   R1   Y   1   =−H   R1   [H   11   U   1   +H   12 U 2 ]. 
       
     
     By solving for U 1 , 
     
       
         
           U 
           1 
           +H 
           R1 
           H 
           11 
           U 
           1 
           =−H 
           R1 
           H 
           12 
           U 
           2 
         
       
     
     
       
         
           
             
               U 
               1 
             
             = 
             
               
                 
                   - 
                   
                     H 
                     R1 
                   
                 
                  
                 
                   H 
                   12 
                 
                  
                 
                   U 
                   2 
                 
               
               
                 1 
                 + 
                 
                   
                     H 
                     R1 
                   
                    
                   
                     H 
                     11 
                   
                 
               
             
           
         
                 
         
             
         
      
     
     and insertion into equation (1 a ), the equation                H   20     =       H   22     -         H   12          H   21          H   R1         1   +       H   R1          H   11                     (   2   )                                
     follows. The controller two REG 2 , illustrated by its transfer function H R2 , thus does not “see” the main path H 22  in the initialization but the equivalent path H 20 , which differs from H 22  given couplings H 11  and H 21 , that do not disappear and that is dependent on the setting of the controller one REG 1 . 
     Analogous to equation (1), the equivalent path H 10 , which is dependent on the controller two REG 2 , follows for the controller one REG 1  instead of the main path H 11 :                H   10     =       H   11     -         H   12          H   21          H   R2         1   +       H   R2          H   22                     (   3   )                                
     The following difficulty thus derives in the initialization of the control structure: when, given an open, second control circuit (H R2 =0), the controller one REG 1  is parameterized with the transfer function H R1 , then this is set for the main path H 11 , since the back term in equation (3) is eliminated with H R2 =0. 
     The subsequent parameterization of the controller two REG 2  is based on equivalent path H 10  according to equation (2). When the second control circuit RK 2  is closed, then the behavior of the first control circuit RK 1  is detuned, since H R2 =0 is no longer valid. The first control circuit RK 1  is readjusted, which is turn leads to a detuning of the second control circuit RK 2  due to the change of H R1  in equation (2). 
     The initialization is thus time-consuming since, readjustment must be repeatedly carried out and, one must furthermore always wait for the stationary condition. 
     The invention enables a setting (parameterization) of the controller parameters without manual repeated adjustments and without waiting for the stationary condition. An initialization tool is required as basis wherein a controller is allocated to a control system, as expressed by the operator equation 
     
       
           H   R   =ℑ{H}   (4) 
       
     
     This allocation can take place with a neural network having a preceding data pre-precessing: in the data pre-processing, characteristic quantities are first formed from the measured transfer function of the control system; the neural network then determines the appertaining controller parameters from these characteristic quantities (neural parameterization). 
     After an adequate number of iterations, the initialization method then approximates—in the boundary case—the equations 
       H   R1   =ℑ{H   10 }  (5) 
     
       
           H   R1   =ℑ{H   20 }  (6) 
       
     
     i.e., the two controllers are respectively designed for equivalent paths. By insertion of equation (3) and equation (2) into equations (5) and (6), the following follows:                H   R1     =             {       H   11     -         H   12          H   21          H   R2         1   +       H   R2          H   22             }               (   7   )                   and                   H   R2       =             {       H   22     -         H   12          H   21          H   R1         1   +       H   R1          H   11             }         ,                             i   .   e   .              0     =               {       H   11     -         H   12          H   21          H   R2         1   +       H   R2          H   22             }       -     H   R1                                 and                 0     =               {       H   22     -         H   12          H   21          H   R1         1   +       H   R1          H   11             }       -       H   R2     .               (   8   )                                
     Since H 11 , H 12 , H 21  and H 22  are given and H R1  and H R2  are sought, it becomes clear that equation (7) and equation (8) cannot be solved independently of one another but represent a non-linear equation system for determining H R1  and H R2 . 
     When PID controllers are employed both as controller one REG 1  as well as controller two REG 2 , then                    H   R1          (   s   )       =       K   P1          (       1     s                   T   N1         +   1   +       s                   T   V1         1   +     s                     T   V1     /   V             )         ,           (   9   )                     H   R2          (   s   )       =       K   P2          (       1     s                   T   N2         +   1   +       s                   T   V2         1   +     s                     T   V2     /   V             )         ,           (   10   )                                
     whereby K P1 , K P2 , T N1 , T N2 , T V1 , T V2  represent parameters of the PID controllers REG 1  and REG 2  (index 1 or, respectively, 2 for the number of the controller). A parameter V is permanently prescribed, i.e., is not optimized. In the example cited here, equations (7) and (8) represent a non-linear equation system with 6 equations for 6 unknowns. 
     The following steps are required for calculating the controller H R1  given initially predetermined controller coefficients K P2 , T N2  and T V2 : 
     1. Calculation of the controller frequency response according to equation (10) dependent on the controller parameters K P2 , T N2  and T V2 ; 
     2. Calculation of the equivalent path frequency response H 10  by insertion of the frequency responses H R2  and of the path (H 11 , H 12 ,H 21  and H 22 ) for s=jω in equation (3); 
     3. Controller calculation based on the frequency response H 10 (jω) of the equivalent path (this step is formally covered with equation (5)). 
     These steps are combined to form a sub-program. The interface of the sub-program is formally described by          [           K   P1               T   N1               T   V1           ]     =       F        (       H   11     ,     H   12     ,     H   21     ,     H   22     ,     K   P2     ,     T   N2     ,     T   V2       )       .                            
     The sub-program explains equation (5) for the case of a parameterization of PID controllers. This sub-program can likewise be employed for the presentation of equation (6), i.e.,          [           K   P2               T   N2               T   V2           ]     =       F        (       H   22     ,     H   21     ,     H   12     ,     H   11     ,     K   P1     ,     T   N1     ,     T   V1       )       .                            
     A compressed presentation can likewise be recited, i.e.,                [           K   P1               T   N1               T   V1               K   P2               T   N2               T   V2           ]     =       [           F        (       H   11     ,     H   12     ,     H   21     ,     H   22     ,     K   P2     ,     T   N2     ,     T   V2       )                 F        (       H   22     ,     H   21     ,     H   12     ,     H   11     ,     K   P1     ,     T   N1     ,     T   V1       )             ]     .             (   11   )                                
     This non-linear equation system can be solved with an arbitrary standard method, for example with a Newton method. 
     The above comments can be generalized to the case of a process controlled system having n inputs and n outputs. The following equation system is formulated for this purpose (PID controllers are again employed):                [           K   P1               T   N1               T   V1               K   P2               T   N2               T   V2             ⋮             K   Pn               T   Nn               T   Vn           ]     =       [           F        (       H   11     ,     H   12     ,   …              ,     H   nn     ,     K   P2     ,     T   N2     ,     T   V2     ,     K   P3     ,     T   N3     ,     T   V3     ,   …              ,     K   Pn     ,     T   Nn     ,     T   Vn       )                 F        (       H   22     ,     H   21     ,   …              ,     H   nn     ,     K   P1     ,     T   N1     ,     T   V1     ,     K   P3     ,     T   N3     ,     T   V3     ,   …              ,     K   Pn     ,     T   Nn     ,     T   Vn       )               ⋮             F        (       H   nn     ,     H   n2     ,   …              ,     H   11     ,     K   P2     ,     T   N2     ,     T   V2     ,     K   P3     ,     T   N3     ,     T   V3     ,   …              ,     K   P1     ,     T   N1     ,     T   V1       )             ]     .             (   12   )                                
     As in the case of the process controlled system having two inputs and two outputs, each line on the right side of equation (12) stands for a sub-program that is called with different parameter values. The parameters of a k th  line (1≦k≦n) is obtained from the first line, in that “k” is written instead of “1” and “1” is written instead of “k” in the indices without modifying the sequence of the parameters. 
     Equivalent path H k0  is calculated below in the general case with reference to FIG.  4 . 
     A generalized control system is given: 
     
       
             Y = H · U     (13) 
       
     
     
       
         
           
             
               
                 
                   
                     
                       with 
                        
                       
                           
                       
                        
                       
                         Y 
                         _ 
                       
                     
                     = 
                     
                       
                         
                           [ 
                           
                             
                               
                                 
                                   Y 
                                   1 
                                 
                               
                             
                             
                               
                                 
                                   Y 
                                   2 
                                 
                               
                             
                             
                               
                                 ⋮ 
                               
                             
                             
                               
                                 
                                   Y 
                                   n 
                                 
                               
                             
                           
                           ] 
                         
                          
                         
                             
                         
                          
                         and 
                          
                         
                           
                               
                           
                            
                           
                               
                           
                         
                          
                         
                           U 
                           _ 
                         
                       
                       = 
                       
                         [ 
                         
                           
                             
                               
                                 u 
                                 1 
                               
                             
                           
                           
                             
                               
                                 U 
                                 2 
                               
                             
                           
                           
                             
                               ⋮ 
                             
                           
                           
                             
                               
                                 U 
                                 n 
                               
                             
                           
                         
                         ] 
                       
                     
                   
                   , 
                 
               
               
                 
                   (13a,b) 
                 
               
             
             
               
                 
                   
                     and 
                      
                     
                       
                           
                       
                        
                       
                           
                       
                     
                      
                     
                       H 
                       _ 
                     
                   
                   = 
                   
                     
                       [ 
                       
                         
                           
                             
                               H 
                               11 
                             
                           
                           
                             
                               H 
                               12 
                             
                           
                           
                             … 
                           
                           
                             
                               H 
                               
                                 1 
                                  
                                 n 
                               
                             
                           
                         
                         
                           
                             
                               H 
                               21 
                             
                           
                           
                             
                               H 
                               22 
                             
                           
                           
                             … 
                           
                           
                             
                               H 
                               
                                 2 
                                  
                                 n 
                               
                             
                           
                         
                         
                           
                             ⋮ 
                           
                           
                             ⋮ 
                           
                           
                             ⋰ 
                           
                           
                             ⋮ 
                           
                         
                         
                           
                             
                               H 
                               n1 
                             
                           
                           
                             
                               H 
                               n2 
                             
                           
                           
                             … 
                           
                           
                             
                               H 
                               nn 
                             
                           
                         
                       
                       ] 
                     
                     . 
                   
                 
               
               
                 
                   (13c) 
                 
               
             
           
         
                 
         
             
         
      
     
     The transfer function of the equivalent path H k0  is sought:                H   k0     =       Y   k       u   k               (   15   )                                
     which is the transfer function that derives given an open k th  control circuit and other control circuits that are otherwise closed. 
     The following vectors are introduced for this purpose (the raised ‘T’ indicates that it relates to a transposed vector): 
         Y     k   =[Y   1    . . . Y   k−1    Y   k+1    . . . Y   n ] T   (16) 
     
       
             U     k   =[U   1    . . . U   k−1    U   k+1    . . . U   n ] T   (17) 
       
     
     
       
             h     sk   =[H   1k    H   2k    . . . H   k−1,k    H   k+1,k    . . . H   nk ] T   (18) 
       
     
     
       
             h     kz   =[H   k1    H   k2    . . . H   K,k−1    H   k,k+1    . . . H   kn ]  (19), 
       
     
     where the following apply 
     Y k  output signals without Y k , 
     U k  setting signals without U k , 
     h sk  the k th  column of H without H kk ; and 
     h kz  the k th  row of H without H kk . 
     Furthermore, a remainder matrix without a k th  row and without a k th  column is calculated:                               H   _     k     =       [           H     1   ,   1             H     1   ,   2           …         H     1   ,     k   -   1               H     1   ,     k   +   1             …         H     1   ,   n               ⋮       ⋮       ⋰       ⋮       ⋮       ⋰       ⋮             H       k   -   1     ,   1             H       k   -   1     ,   2           …         H       k   -   1     ,     k   -   1               H       k   -   1     ,     k   +   1             …         H       k   -   1     ,   n                 H       k   +   1     ,   1             H       k   +   1     ,   2           …         H       k   +   1     ,     k   -   1               H       k   +   1     ,     k   +   1             …         H       k   +   1     ,   n               ⋮       ⋮       ⋰       ⋮       ⋮       ⋰       ⋮             H     n   ,   1             H     n   ,   2           …         H     n   ,     k   -   1               H     n   ,     k   +   1             …         H     n   ,   n             ]     .               (   20   )                                
     What follows is: 
     
       
           Y   k   =H   k1   U   1   +H   k2   U   2   +. . . +H   kn   U   n   =H   kk   U   k   +h   kz   U   k   (21). 
       
     
     To eliminate the vector U k , U k  can be expressed via the controller Y k , i.e., with 
     
       
           H   rk =diag[ H   R,1   H   R,2    . . . H   R,k−1   H   R,k+1    . . . H   R,n ]  (22) 
       
     
     applies 
         U     k   =− H     Rk   · Y     k   (23) 
     and with 
     
       
             Y     k   = H     k     U     K   + h     sk   U   k   (24) 
       
     
     then it follows that 
     
       
         
           
             U 
           
           k 
           =− H   
           Rk 
           
             H 
           
           k 
           
             U 
           
           k 
           − 
         
       
     
     
       
           H   Rk     h     sk   U   k ⇄ 
       
     
     
       
         ⇄(   E + H     Rk     H     k )   U     k   =   
       
     
     
       
         − H   Rk     h     sk   U   k   ⇄   
       
     
     
       
         ⇄   U     k =−(   E +   
       
     
     
       
           H   Rk     H     k ) −1     H     Rk     h     sk   U   k   (25) 
       
     
     When equation (25) is inserted into equation (21) the following derives: 
     
       
           Y   k   =H   kk   U   k   −h   kz ( E+H   RK   H   K ) −1   
       
     
     
       
         H RK   h   sk   U   K   (26) 
       
     
     When, further, equation (15) is co-considered for the sought equivalent path, then it follows that 
     
       
           H   k0   =H   kk   −h   kz ( E+H   RK   H   K ) −1   
       
     
     
       
         H RK   h   sk   (27) 
       
     
     equation (3) is obtained for n=2 from equation (27) given k=1 and equation (2) is obtained given k=2. 
     The equations 
     
       
           k =1: 0 ℑ{h   11   − h     1z (   E + H     R1     H     1 ) −1   
       
     
     
       
         
           H 
           R1 
           
             h 
           
           s1 
           }−H 
           R1 
         
       
     
     
       
         k=2: 0 =ℑ{H   22   −h   2z (   E +   
       
     
     
       
           H   R2     H     2 ) −1     H     R2     h     s2   }−   
       
     
     
       
         H R2   (28) 
       
     
     . . . 
     
       
           k=n : 0 =ℑ{H   nn   − h     nz (   E +   
       
     
     
       
           H   Rn     H     n ) −1     H     Rn     h     sn   }−H   Rn   
       
     
     arise for the case of n controllers. 
     FIG. 5 shows the steps of the invention. In step 5 a , the transfer functions of the controllers are determined in general dependent on the controller parameters (see equations (9) and (10)). Subsequently, measurements of the reply signals for each control system including all couplings are implemented in step 5 b  and a path transfer matrix H(s) is calculated, (see equation (1).) The couplings between the control systems contained in the process controlled system are taken into consideration in the form of an equivalent system for each control circuit in step 5 c  (see equations (2) and (3)). Finally, the controller parameters for each controller are identified in step 5 d  from the equivalent functions, in that the non-linear equation system that has arisen is iteratively solved (see equation (11)). To that end, the individual equations are brought into a form for zero place calculation (see equation (11) after subtracting the right-hand side), so that this zero place calculation for varying controller parameters is implemented in the form of a sub-program on the computer specifically provided for this purpose. The result is returned from the sub-program and represents all of the sought controller parameters.