=Paper= {{Paper |id=Vol-1638/Paper67 |storemode=property |title=Transferable utility distribution algorithm for multicriteria control in strongly coupled system with priorities |pdfUrl=https://ceur-ws.org/Vol-1638/Paper67.pdf |volume=Vol-1638 |authors=Mikhail I. Geraskin }} ==Transferable utility distribution algorithm for multicriteria control in strongly coupled system with priorities == https://ceur-ws.org/Vol-1638/Paper67.pdf
Mathematical Modeling


   TRANSFERABLE UTILITY DISTRIBUTION
ALGORITHM FOR MULTICRITERIA CONTROL IN
STRONGLY COUPLED SYSTEM WITH PRIORITIES

                                       M.I. Geraskin

                   Samara National Research University, Samara, Russia



       Abstract. The problem of multicriteria control on the basis of aggregate utility
       maximizing is considered. The system utility distribution algorithm in case of
       unequal criteria set by priority vector, has been developed for strongly coupled
       technical and organized systems with transferable utility. The Pareto efficiency
       is proved for the distribution generated on the algorithm basis.

       Keywords: multicriteria control, aggregate utility, strongly coupled system,
       transferable utility, distribution algorithm, criteria priorities


       Citation: Geraskin MI. Transferable utility distribution algorithm for mul-
       ticriteria control in strongly coupled system with priorities. CEUR Workshop
       Proceedings, 2016; 1638: 542-551. DOI: 10.18287/1613-0073-2016-1638-542-
       551


Introduction

For technical and organizational systems there are now two ways to solve the problem
of multicriteria optimal control, efficiency of which is evaluated with the help of sev-
eral criteria. The first way is based on the Pareto efficiency criterion [1] with setting-
up of a corresponding set. Another way is Neumann-Morgenstern expected utility
maximization [2] in one form or another [3-6]. The latter approach ensures the
uniqueness of the solution of the problem. However, its use is justified in case the
strongly coupled system model is adequate [7] for the particular group of criteria. The
system is regarded as strongly coupled in case the criteria of its elements are compati-
ble with the axioms [8] of preference (reflexivity, connection, transitivity, etc.). It
allows for the criteria vector the reduction towards the scalar problem of aggregate
utility optimization. Solution of this problem is Pareto-efficient [10] in case it corre-
sponds to the minimax criterion [9]. Utility aggregation may be based on physical and
technical properties of systems. In particular, in the terminal control problems [11]
deviations from the phase coordinates target values are aggregatively set out in the
form of a metric (distance) in the criteria hyperspace. In the problems of multicriteria
statistical estimation [12] the moments of random variables are set out in the from of
the variation coefficient. In general, the criteria aggregation is based on the properties

Information Technology and Nanotechnology (ITNT-2016)                                     542
Mathematical Modeling                                                 Geraskin MI. Transferable…


of the system integrity and hierarchal pattern [13], for example, interrelated technolo-
gies implementing systems. Due to the production process uniformity the elements of
criteria are reduced to the aggregated criterion of the system.
The special case of a strongly coupled system is the transferable utility system with
criteria expressed in a single measuring element, for example, in case of homogene-
ous technologies [14] performance optimization. In general any system of criteria may
be formally reduced to the transferable system by means of setting out [9] criteria in a
non-dimensional form by means of corresponding normalization. Transferable utility
allows for the system elements criteria vector values redistribution (transfer), corre-
sponding to a certain Pareto-efficient control. Consequently, the problem of mul-
ticriteria control is reduced to the problem of optimal utility distribution between the
system elements on the basis of the aggregated criterion. Solution of the latter prob-
lem will also be Pareto-efficient.
Transferable utility distribution algorithms are intended for the systems with elements
having equisignificant criteria (of anonymous elements). In this case the aggregation
process doesn't take the criteria priorities into account. In particular, Pareto-efficiency
is justified [15] in case of an algorithm that establishes the minimum between the
element optimum and average non-distributed utility of the system. In case the system
elements criteria have different priorities, the distribution algorithms are reduced to
the median multicriteria choice [16]. However, in a general Pareto- inefficient case,
mechanisms of anonymous symmetric coalitions [17] have been developed. In particu-
lar cases these mechanisms are Pareto-efficient. We shall hereinafter consider the
problem of development of a Pareto-efficient algorithm of transferable utility distribu-
tion between elements with priorities on the basis of the results [14], obtained for
anonymous elements.


1       Problem definition

Let us consider the system with the elements efficiency determined by the criteria
                             
vector f  f k u , k  K . Let the scalar optimum u k0 for the criteria of k-element be
determined on the basis of the condition

u k0  Arg max f k u , f k0  f k uk0 , k  K ,                                              (1)
         u k U k


where u – control, U – feasible region, f k u , k  K – k-element efficiency criteri-
on.
Let us form the criteria minimums vector, required for further normalization, from the
values obtained in case of scalar optimums of other elements:

         jK \ k
                     
f kmin  min f k u 0j , k  K .                                                                  (2)

Subject to (1), (2) the range of criteria values variation equals to f k u k   f k0 , f kmin  .




Information Technology and Nanotechnology (ITNT-2016)                                           543
Mathematical Modeling                                         Geraskin MI. Transferable…


Let us introduce the hypothesis of benevolence [7], in accordance to which the ele-
ments of a strongly coupled system maximize the utility of the whole system. Control
is determined on the basis of the maximum value of the sum of criteria in case of
optimums of elements. In this context the unique solution of a multicriteria problem is
formed. This solution is Pareto-efficient according to the Herneyer's theorem [1,9]. Let
us denote this solution by u * and define it as follows:

                                             
u *  Arg max  max f j uk , f k*  f k u * , k  K ,                             (3)
              jK  k k       
          kK       u U



where f k* – k-criterion value in case of the strongly coupled system optimum.
Let us define the strongly coupled system maximum aggregate utility F u *  in case
of control (3) as exceeding of criteria optimums over their minimum values:

              
F u *   f k*  f kmin .                                                            (4)
         k K


Let us introduce the system utility distribution vector X  xk , k  K  , belonging to
the acceptable set

                                                  
                                         
FX   X  x k , k  K :  x k  F u * , X  Rk .
                          kK                     

Let us denote the criteria vector, formed as the result of aggregate utility distribution
 f k xk  , and, with account taken of (4), let us define it as follows:

f k xk   f kmin  xk , k  K .                                                     (5)

Criteria functions (5) constitute the efficiency criteria of utility distribution between
elements in case of Pareto-efficient control (3). Let us formulate the multicriteria
problem of optimal utility distribution as follows:

                    
X *  Arg max f k x k , k  K .                                                    (6)
         X FX


Let us introduce the criteria priorities vector  k , k  K , the components of which de-
termine their relative significance within the framework of the system aggregate utili-
ty. As hereinafter the problem (6) is formulated as minimization of criteria deviations
from their scalar maximums, it is expedient to pass on to inverse priorities
 k  1   k , k  K , that characterize the relevant significance of other criteria in
comparison with this particular criteria. Let us define the vector of inverse priorities
as follows

   k  R , k  K ,   k  1 .                                                  (7)
                         kK




Information Technology and Nanotechnology (ITNT-2016)                                544
Mathematical Modeling                                                       Geraskin MI. Transferable…


As it follows from (1), (3) that deviations of criteria (5) from scalar maximums (1) are
negative, let us define the normalized values of criteria f k xk ,  k  as squared relevant
                                                             ~

deviations of criteria (5) and (1) with account taken of inverse priorities:

                    f x   f k0 
                                        2

f k xk ,  k    k k mink
~
                                0 
                                     , k  K.                                                        (8)
                     k fk  fk 

Firstly, normalization (8) reduces criteria to the range f k xk ,  k   0,1 that makes it
                                                         ~

possible to operate on the criteria (8) as transferable ones; secondly, it prescribes the
utility transitivity f xk ,  k   f xi ,  k  k   i , k  K \ i in respect of which higher prior-
                     ~               ~

ity criteria have lesser utility loss.
Let us consider a strongly coupled system, in which control is defined according to
(3). It also ensures utility transferability with the help of (8). It allows for [7] the re-
duction of the multicriteria problem (6) towards the scalar optimization problem by
means of introduction of an aggregated utility function in the form of the relative
departures product (8):

X *  Arg min  f k  xk ,  k  .
                ~                                                                                     (9)
          X FX       kK


Criterion in the problem (9) constitutes an aggregated function of the system utility.
Moreover, as it is in a normalized form (8) that partial criteria components form part
of (9), the effect distribution that satisfies (9) provides [9] with the maximum efficient
solution of the multicriteria problem (6) in respect of minimax

     
f X  min max f k  xk ,  k  ,
~* *          ~
             X FX       kK


that corresponds, as it is shown in [10], to Pareto efficiency.


2        Utility distribution algorithm

Let us formulate the system (4) utility distribution that is optimal in respect of criteria
(9) in the form of the following assertion: distribution

                                                 f  f 0i   k f i*0i 
xk* 
            1
                    k 
                                
                        F u *   f i*   0 k
                                                         i              
                                                                          ,k K
                                                                                                    (10)
        1                      iK \ k  iK \ k
                                                                         
          iK \ k   i

is the solution of the problem (9).
Proof: differentiating Lagrange function, written down for the problem (9) with con-
sideration of (8)




Information Technology and Nanotechnology (ITNT-2016)                                                545
Mathematical Modeling                                                                                                                                                      Geraskin MI. Transferable…


                                               f x   f k0 
                                                                                                                                                           2
                                                                               
L   f k xk ,  k     F   xk     k k mink
        ~
                                                                    F   xk ,
                                                              0 
    kK                         kK            
                                        kK  k kf       f k             kK   

we arrive at the system of necessary conditions of optimality

                    k f k xk   f k0          f x   f     0, k  K                              0 2
L/xk  2
                    f      min
                                         f 
                                                 f  f 
                                               0 2
                                                   iK \ k
                                                                             i       i         i
                                                                                             min
                                                                                                           k
                                                                                                           0 2
                                                                                                                                                                                                           (11)
                       k    k                 k                                  i       i                 k



L/  F   xk  0 .                                                                                                                                                                                       (12)
                    k K


Eliminating Lagrange multiplier from (11), we get

 k f k xk   f k0       f x   f    f x   f    f x   f  , k , n  K , n  k ,
                                                                                             0 2                                           0                                                         0 2

                                                 i
                                                          i         i                   i                 n           n      n           n                    i       i           i             i

 f k   k
          min
                   f      f  f   f  f   f  f 
                           k
                            0 2
                                iK \ k
                                                      i        i
                                                                   min
                                                                                     i
                                                                                         0 2
                                                                                                               i       n
                                                                                                                        min
                                                                                                                                        n
                                                                                                                                         0 2
                                                                                                                                                 iK \ n
                                                                                                                                                                   i       i
                                                                                                                                                                               min
                                                                                                                                                                                             i
                                                                                                                                                                                                 0 2




whence we arrive at the system

 f x   f   f x   f    f x   f   f x   f  , k , n  K , n  k ,
     k   k     k           k
                            0
                                              i   i        i             i
                                                                             0 2
                                                                                                   n   n           n              n
                                                                                                                                   0
                                                                                                                                                  i   i    i                   i
                                                                                                                                                                                   0 2

                                iK \ k                                                                                                iK \ n


or

 k f k xk   f k0   n f n xn   f n0 , n  k  K (13)

Inserting system utility (4) and utility of elements (5) into (13) we arrive at the distri-
bution (10). Differentiating (11), we shall get the system of sufficient conditions of
optimality

L 
  //  f x   f   0, k  K ,
                           i    i         i           i
                                                          0 2

  xk
       f  f 
             iK \ k        k       i
                                        min
                                                  i
                                                      0 2



that are met in case of any fi xi , fi 0 , fi min .
As distribution (10) maximizes the aggregate criterion (9), consequently, in accord-
ance with the results of [10], it is Pareto-efficient.
On the basis of the solution (10) of the multicriteria problem (9) let us lay down the
aggregated utility distribution algorithm in the form of the series of the following
steps.

1. Input of the given data: functions of criteria f  f k u , k  K and restrictions de-                                                                                              
   fining the ranges of acceptable controls Uk.
2. Determining scalar optimums of criteria (1) and minimums of criteria (2).
3. Determining optimums of criteria in the system (3) and utility of the system 4 (4).



Information Technology and Nanotechnology (ITNT-2016)                                                                                                                                                      546
Mathematical Modeling                                                Geraskin MI. Transferable…


4. Defining the priorities vector (7).
5. Calculating the utility distribution (10) and the optimal criteria vector (5).
The algorithm realization requires for the criteria feasible region to be not empty
                                                                                
U i  U k  , k  K \ i and also for the functions f  f k u , k  K to have finite
maximums in corresponding feasible regions U k . These conditions are met in practi-
cal problems, because as a rule the regions U k , are compact and usually they coincide
for different criteria of one system.


3         Utility optimal distributions modelling

Modelling was done in the context of the three-element system of electronics market,
the elements of which are retailer, bank, insurer, with the functions f k u , k  K
having the form

         
f k u k  p k u kbk 1  c k u k , k  K ,

where f k , u k , c k – k-element utility (profit), output and marginal costs, p k , bk –price
trends parameters on corresponding markets. Scalar optimums of elements (1) were
defined according to the formulas:
                       1
        ck       bk
u 0
                  ,k K
                         .
    pk bk  1 
    k




Parameters of the model, criteria (1) scalar optimums, criteria (2) minimums, criteria
(3) system optimums are given in Table 1; system (4) utility equaled F=6989. Calcu-
lation of parameters of the model is given [14, 18] for the electronics retailer LLC
“Eldorado”, insurer OJSC “Insurance company “Renaissance”” and bank LLC
“Home Credit and Finance Bank”, that were assigned indices k=1,2,3 corresponding-
ly. In the considered system element 1 gets maximum utility (utility leader), parame-
ter u10 for which corresponds to the volume of sales of goods and is expressed in phys-
ical terms, and utility obtained by elements 2,3 (utility outsiders), are significantly
lower. However, their participation is required for the system to be integral. That is
why the system utility distribution problem becomes. Outputs u20 , u30 correspond to
the volume of insurance and lending and are expressed in monetary terms.
Let us suppose that elements may assign their priorities themselves in the course of
integration into the system, following two variants of strategic behaviour.
Firstly, lets consider the variant of non-cooperative game [19], elements wherein don't
form coalitions, they assign priorities simultaneously and independently. In this case
the Cournot game model [20], symbolized by the letter «K», has the form
                        
 k  K   arg max  k , k  K and leads to the priorities vector  k  K    *K  , k  K . The
               k 1




Information Technology and Nanotechnology (ITNT-2016)                                          547
Mathematical Modeling                                                              Geraskin MI. Transferable…


resulting distribution x k* K  , k  K (fig. 1) corresponds to the elements equisignifi-
cance. In this case the element 1 gets the predominant share of distributed utility, and
the share of utility distributed in favour of the elements 2,3 is reduced concurrently
with the growth of the equilibrium values of priorities  *K  . Model of game with the
institutional leader sensu Stackelberg [21], symbolized by the letter «S», takes the

                         k 1
                                  
form  k S   arg max  k ,  l S   arg max  l , k  K \ l , where l-element is the leader;
                                                   k  l 1

the model results in the vector of priorities  k S    *S  ,  l*S    *S  , k  K \ l . Distribu-
tion x k* S  , k  K (fig. 1) shows that the leader sensu Stackelberg (l=3) gets the ad-
vantage in spite of being the utility outsider.

                            Table 1. Parameters of the model and criteria vector

     Parameter of the                                            Element
        model                               1                         2                          3
   bk                                   -0,09                     -0,19                     -0,165
   ck                                  22000                       0,06                      0,025
    pk                                 50000                     0,780                       0,350
   u k0                                  3210                   240548                    2962123
    f k0                                6985,                          3,4                   14,6
                                            0
    f k*                                6985,                          3,0                   12,0
                                            0
    f kmin                                7,0                          1,2                     2,4
    f f
     k
      *
             k
              min
                                        6978,                          1,8                     9,6
                                           0
                    x
                         8000
                         7000
                         6000
                         5000
                         4000
                         3000
                         2000
                         1000
                            0
                        -1000
                        -2000
                                 0      0,05      0,1           0,15         0,2   0,25     0,3 γ

                                         x*1(K)                  x*2(K)            x*3(K)
                                         x*1(S)                  x*2(S)            x*3(S)


             Fig. 1. Utility distribution structure in case of non-cooperative behaviour
                        1K 2K 3K, S
                                             
                                             1  
                                                2 ,
                                                 
                                                 S S
                                                    
                                                    3  
                                                        
                                                        1
                                                        S 
                                                          , 
                                                            1



Information Technology and Nanotechnology (ITNT-2016)                                                    548
Mathematical Modeling                                                    Geraskin MI. Transferable…


Secondly, let us consider the variant of cooperative game [22] symbolized with the
letter “S”. In this case elements may form coalitions. The game model of the form
                                                      
 k C   arg max  k ,  l C   arg max  l , l  L , k  K \ L , where index l  L desig-
               k 1                    k  l 1

nates the elements of coalition, leads to the vector of priorities of the form
 l C    *C  ,  *C    k*C  , k  K \ l . Distribution x k*C  , k  K (fig. 2) shows that in
comparison with non-cooperative behaviour (fig.1), the coalition of elements (1,2)
significantly increases the utility leader's utility and increases the outsider's utility
insubstantially. Coalition of utility outsiders (2,3) leads to a drastic utility redistribu-
tion in their favour.
                 x
                        9000
                        8000
                        7000
                        6000
                        5000
                        4000
                        3000
                        2000
                        1000
                           0
                       -1000 0       0,1         0,2    0,3      0,4      0,5        0,6   γ

                                 x*1(С),(1,2)            x*2(С),(1,2)           x*3(С),(1,2)
                                 x*1(С),(2,3)            x*2(С),(2,3)           x*3(С),(2,3)


               Fig. 2. Utility distribution structure in case of cooperative behaviour
              
              1
              C 1
                ,
                2   
                         
                              
                              , 
                               2
                               C 1
                                 ,
                                 2   
                                      
                                          
                                           3
                                           C
                                            ,1
                                             ,
                                             2
                                              1        2
                                                       C 2
                                                         ,
                                                         3       3
                                                                 C 2
                                                                   ,
                                                                   3            1
                                                                                C 2
                                                                                  ,
                                                                                  3



In all the cases examined the system utility was fully distributed between the ele-
ments. This confirms Pareto efficiency of distributions that are formed on the
basis of the developed algorithm. Modelling showed that the developed algorithm
allows to analyze the impact of the priorities vector on the distribution structure
and to build the vector of the priorities that express different aims of the system
integration.


Conclusion

The algorithm of solving multicriteria problems of control in strongly coupled sys-
tems has been developed. This algorithm allows to determine the optimal system
aggregated utility distribution in case of inequisignificant efficiency criteria on the
basis of the minimum of deviations of criteria from scalar optimums. The algorithm
application sphere includes technical and organizational systems, efficiency criteria of
which are quantifiable and preference relations of which have been determined. In
case the system elements utilities are transferable, which is in general ensured by the



Information Technology and Nanotechnology (ITNT-2016)                                                549
Mathematical Modeling                                             Geraskin MI. Transferable…


criteria normalization. The algorithm is based on the normalized criteria aggregation
and in case of the utility function optimization with account taken of priorities, it
leads to the minimax distribution (of the guaranteed result). Consequently, the algo-
rithm determines the Pareto-efficient solution of the multicriteria problem.
The use of algorithm allows to reduce the multicriteria problem of optimization in the
control space to the set of scalar problems of optimal control and the problem of mul-
ticriteria distribution in the criteria space. As the result the unique solution is estab-
lished, and, moreover, solving of the multicriteria problem is significantly simplified.


References
 1. Nogin VD. Decision making in multicriteria area: numerical approach. Moscow: “Fizmat-
    lit” Publisher, 2005. [in Russian]
 2. Neumann G, Morgenstern O. Theory of games and economic behavior. Moscow: “Nauka”
    Publisher, 1970. [in Russian]
 3. Saaty TL. Decision making in functions and feedback. Analitical grids. Moscow: “LKI”
    Publisher, 2008. [in Russian]
 4. Fishhof VG. Goltein BG., Shapiro ZR. Subjective expected utility: decision making mod-
    el. Procedures for multicriteria evaluation facilities. Moscow: “VPIISI” Publisher, 1984; 9:
    24-46. [in Russian]
 5. Caverny JP, Bar-Hillel M, Barron FN, Jungerman H. Contributions to decision research.
    North-Holland, 1993.
 6. Smith GR, Speiser F. Logical decision: multimeasure decision analysis software. Golden.
    CO: PDQ Printing, 1991.
 7. Novikov D. Theory of Control in Organizations. New York, Nova Science Publishers,
    2013.
 8. Roy B. Multy criteria methodology for decision aiding. Dordrecht: Kluwer academic pub-
    lisher, 1996.
 9. Mashunin YK. Methods and models of vector optimization. Moscow: “Nauka” Publisher,
    1986. [in Russian]
10. Homenuk VV. Elements of the theory of multi-objective optimization. Moscow: “Nauka”
    Publisher, 1983 [in Russian]
11. Geraskin MI. Lazarev YN. Terminal Control descent of the spacecraft in the atmosphere
    and restrictions on movement modes. Proceedings of the Academy of Sciences. Theory
    and Control Systems. 2001; 5: 168-174.
12. Balakin VL, Geraskin MI, Lazarev YN. Research spacecraft maneuvering capabilities
    when driving on a suborbital trajectory. Aerospace Engineering and Technology, 1997; 2:
    46-48. [in Russian]
13. Geraskin MI. Optimization models of management of hierarchical corporate systems with
    inter-corporate interactions. Control sciences, 2010; 5: 28-38. [in Russian]
14. Geraskin MI, Manakhov VV. Optimization of interactions in the multi-agent tightly cou-
    pled system “retailer – bank – insurer”. Control Science, 2015; 4: 9-18. [in Russian]
15. Burkov VN, Gorgidze II, Novikov DA. Models and mechanisms for the distribution of
    costs and revenues in a market economy. Moscow: Institute of Control Sciences of Rus-
    sian Academy of Sciences Publisher, 1997. [in Russian]




Information Technology and Nanotechnology (ITNT-2016)                                       550
Mathematical Modeling                                            Geraskin MI. Transferable…


16. Burkov VN, Iskakov MB, Korgin NA. Application of generalized median schemes for the
    construction of the manipulated mechanism multicriterion active expertise. Control scienc-
    es, 2008; 4: 38-47. [in Russian]
17. Bondarik VN, Korgin NA. Mechanisms of resource allocation on the basis of manipulated
    symmetrical anonymous voting procedures with delegation. Control sciences, 2012; 5: 26-
    32. [in Russian]
18. Geraskin MI, Manakhov VV. Analysis of demand curves in the stock and financial mar-
    kets of monopolistic competition, Actual Problems of Economics and Law, 2016; 2(10):
    80-92. [in Russian]
19. Germeyer YuB, Shtilman MS. Non-cooperative repetitive games with arbitrary discount-
    ing. DAN Publisher, 1975. [in Russian]
20. Cournot AA. Researches into the Mathematical Principles of the Theory of Wealth. Lon-
    don, Hafner, 1960. Original 1838.
21. Stackelberg H. Market Structure and Equilibrium: 1st Edition. Translation into English,
    Bazin, Urch & Hill, Springer, 2011.
22. Aumann RJ. A survey of cooperative games without side payments. Essays in Mathemati-
    cal Economics, in Honor of Oskar Morgenstern, 2015: 3-27.




Information Technology and Nanotechnology (ITNT-2016)                                     551