<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Numerical solution of the dynamic incentive problem in discrete time taking into account the learning curve effect</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Oleg Pavlov</string-name>
          <email>pavlov@ssau.ru</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Institute of economics and management Samara National Research University Samara</institution>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2020</year>
      </pub-date>
      <fpage>251</fpage>
      <lpage>254</lpage>
      <abstract>
        <p>-The paper considers the dynamic incentive problem in discrete time, taking into account the learning effect. The task is formulated as a dynamic game between the leader and the performers. To solve the problem, the principle of cost recovery is applied, which reduces the original task to the optimal control problem in discrete time. Numerical solutions of the problem for various models of learning curves are obtained using the Bellman dynamic programming method. Also, the study is conducted of the discount rate's impact on the solution of the incentive problem</p>
      </abstract>
      <kwd-group>
        <kwd>dynamic incentive problem</kwd>
        <kwd>inverse Stackelberg game</kwd>
        <kwd>learning curve effect</kwd>
        <kwd>Bellman dynamic programming</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Copyright © 2020 for this paper by its authors.</p>
      <p>Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0)
program. When making a decision, the center proceeds from
the principle of maximum guaranteed result. As a result, the
initial problem is transformed into the optimal control
problem.</p>
      <p>
        In this article basing on the approach [
        <xref ref-type="bibr" rid="ref1 ref12">1,12</xref>
        ], the dynamic
incentive problem of agents taking into account the learning
curve effect is formulated and numerically solved using the
Bellman dynamic programming method.
      </p>
      <p>II. STATEMENT AND ALGORITHM FOR SOLVING THE DYNAMIC</p>
      <p>INCENTIVE PROBLEM OF AGENTS</p>
      <p>
        A two-level dynamic manufacturing system consisting of
a center and n independent agents is considered. Agents
produce parts from which the finished product is then
assembled. Labor costs and financial incentives for agents
depend only on their own actions. This article applies the
principle of game decomposition [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], which allows to
consider the management of the i-th agent independently and
not to take into account the interaction of agents with each
other. The state of a dynamic production system depends on
the actions of agents, and the center affects the managed
system only through the payment of material remuneration to
agents.
      </p>
      <p>The dynamics of part production by the i-th agent is
described by a discrete equation:
 xt  xt 1  u t , t  1,T , 
where xt is the cumulative production volume of the part in
the time period t, t is the number of the time period, ut is the
production volume of the part in the period t, T is the
quantity of time periods considered.</p>
      <p>Before the start of mass production, we know the number
of manufactured parts, it is as follows:
 x0  X 0 . </p>
      <p>In the final time period, the cumulative volume of parts
must be equal to the specified as follows:</p>
      <p>xT  X 0  R , 
where R is the specified number of parts.</p>
      <p>Restrictions are imposed on the production volume of the
part:
 0  u t  X 0  R  xt 1 , t  1,T . </p>
      <p>The target function of the center is to maximize the
discounted total difference between the income from the
manufactured parts and the costs of the agent’s material
compensation:</p>
      <p>T 1
J p  
t 1 ( 1  r )t
[ p u t   ( xt )]  m a x , 
</p>
    </sec>
    <sec id="sec-2">
      <title>I. INTRODUCTION</title>
      <p>The article discusses the game dynamic task of the
executors performing the production task in the context of
new product development. The development of new
products at industrial enterprises is characterized by the
learning curve effect, which is that the time spent by
employees (laboriousness) on performing multiple repetitive
production operations is reduced.</p>
      <p>The task of executors stimulation is one of the most
important in the management theory. The management (the
center) should choose such an incentive system based on the
forecast of the agent’s actions in order to ensure the
fulfillment of their economic interests. The executor (the
agent) chooses an action (volume of work) based on his
economic interests.</p>
      <p>
        Dynamic problems of interaction of unequal players are
considered in the active systems theory [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], in the
information theory of hierarchical systems [
        <xref ref-type="bibr" rid="ref2 ref3">2–4</xref>
        ] and in the
dynamic games theory developed by international authors
[
        <xref ref-type="bibr" rid="ref10 ref11 ref5 ref6 ref7 ref8 ref9">5–11</xref>
        ]. It should be noted that the stimulation problem in
different theories has received various names. In the active
systems theory it is the incentive task, in publications of
foreign authors on game theory it is the inverse Stackelberg
game, in the information theory of hierarchical systems it is
the Germeyer game.
      </p>
      <p>
        The active systems theory [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ] offers the approach called
the principle of agent’s cost compensation. The center pays
material remuneration to the agent, compensating his costs,
in the case of choosing the optimal planned trajectory of the
center and does not pay material compensation otherwise.
The initial problem is divided into two tasks: the choice of
the incentive system and the solution of the optimal control
problem. In [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ], results are presented that generalize the
theorems from the monograph [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ].
      </p>
      <p>
        The hierarchical systems theory [
        <xref ref-type="bibr" rid="ref2 ref3">2–4</xref>
        ] suggests the
approach that uses the center’s choice of the program of joint
actions with the agent and punishment for deviation from this



where p is the part price,  ( xt ) is the center incentive
function, r is the center discount rate.
      </p>
      <p>The incentive function of the center is a rule in
accordance with which a material remuneration is assigned
to the agent for the amount of work performed. The center
manages the production process through the mechanism of
material incentives  ( xt ) , economically encouraging agents
to fulfill the planned production volumes.</p>
      <p>The discount rate helps to take into account the time
preferences of the center (agent) for the cost of cash flows.
The more distant in time the cash flow, the cheaper it is for
the center (agent).</p>
      <p>The target function of the agent is to maximize the
discounted total difference between material remuneration
and labor costs, expressed in monetary form:</p>
      <p>[  ( xt )  C t ( u t , xt 1 )]  m a x ,  
where r is the agent discount rate, C t ( u t , xt 1 ) is the agent
labor costs.</p>
      <p>Agent labor costs are determined by the following
equation:
 C t ( u t , xt 1 )  s сt u t  
where s is the cost of one hour per agent, ct is the
laboriousness of manufacturing the part.</p>
      <p>
        The dependence of the part laboriousness on the
cumulative production volume is described by various
models of the learning curve given in [
        <xref ref-type="bibr" rid="ref13">13</xref>
        ]-[
        <xref ref-type="bibr" rid="ref15">15</xref>
        ].
      </p>
      <p>In accordance with his economic interests, the agent
selects parts production volumes that maximize his target
function (2). The center’s task is to choose the optimal
incentive system in which the agent will produce such parts
production volumes that maximize the center target function
(1).</p>
      <p>
        To solve the formulated control problem, the principle of
cost compensation is applied [
        <xref ref-type="bibr" rid="ref1 ref12">1, 12</xref>
        ]. The solution algorithm
consists in dividing the initial problem into two tasks:
choosing a compensatory incentive system and solving the
optimal control problem with the objective function equal to
the difference between the center’s income and the agent’s
labor costs.
      </p>
      <p>1. The choice of a compensatory incentive system.</p>
      <p>The center selects a compensatory incentive system,
which consists in compensating the agent costs in the case of
choosing the optimal planned production volume of the
center xtopt and the absence of material payments otherwise:
  ( xt )   C t ( u t , xt 1 ), е с л и xt  xtopt , д л я  t  1,T , 
 0 , е с л и xt  xtopt , д л я  t  1,T .</p>
      <p>2. The solution of the optimal control problem with the
target function equal to the difference between the center
income and the agent labor costs.</p>
      <p>To encourage the agent to choose the planned production
volume, the center pays a material remuneration equal to the
agent costs:




</p>
    </sec>
    <sec id="sec-3">
      <title>Power-based labour input model:</title>
      <p>ct  4 2 ,6 4 xt01,3 . 
Exponential labor input model:</p>
    </sec>
    <sec id="sec-4">
      <title>Logistic labor input model:</title>
      <p>сt  9 ,1 7  6 ,1 6 e 0 ,03 xt1 . 
 1 
 1  0 , 0 1 7 e 0 ,05 xt1  . 
сt  5 5 ,1 0  3 6 , 6 1 </p>
      <p>To solve the problem, the following data was used: the
number of time periods T=12 months, the production volume
of parts R=240 pcs., production experience before serial
 production x0  1 pcs. The discrete step of changing the
parts production volume when implementing the dynamic
programming method is 1 pcs.</p>
      <p>Numerical solutions of the optimal control problem for
power-based, exponential and logistic models of labor input
are presented in Fig. 1-3. The figures show the optimal
trajectories of cumulative production volumes for various
discount rates.
  ( xt )  C t ( u t , xt 1 )  </p>
      <p>We substitute the formula (4) into the target function of
the center, taking into account (3):
 J p  T 1 [ p  s сt ] u t  m a x . 
t 1 ( 1  r )t</p>
      <p>Since the part price p is constant, the center can increase
his profit only by minimizing the total cost of paying the
agent’s material remuneration. The target function of the
center will take the following form:





T 1
t 1 ( 1  r )t sсt u t  m in . </p>
      <p>J p  </p>
      <p>Thus, the initial dynamic incentive task is reduced to the
optimal control problem:</p>
      <p>T 1
t 1 ( 1  r )t sсt u t  m in . 
J p  

 xt  xt 1  u t , t  1,T ,  
 x0  X 0 ,  
 xT  X 0  R ,  
 0  u t  X 0  R  xt 1 , t  1,T .  </p>
      <p>The center’s task is to select the optimal production
volumes of parts u topt , taking into account restrictions (9),
under which the production process (6) will switch from the
initial state (7) to the final state (8) and the minimum of the
center’s target function (5) will be achieved.</p>
      <p>
        The formulated optimal control problem (5)-(9) was
solved using the Bellman dynamic programming method
[
        <xref ref-type="bibr" rid="ref16">16</xref>
        ], implemented in the pascal programming language.
      </p>
      <p>III. THE RESULTS OF THE NUMERICAL SOLUTION OF THE</p>
      <p>DYNAMIC INCENTIVE PROBLEM OF AGENTS</p>
      <p>The numerical solution of the optimal control problem is
carried out on the example of the production of parts of the
enterprise Salut JSC. According to the enterprise data,
regression models of the of laboriousness manufacturing
parts are constructed: power, exponential and logistic.</p>
      <p>250
e
iltvau ino .)scp125000</p>
      <p>(
um tcu e100
lc d m
a ro lu 50
itpm p vo 0
O</p>
      <p>From an analysis of Fig. 1-2, it follows that for a
powerbased and exponential model of labour input, a convex curve
is the optimal trajectory of the cumulative production
volume. The optimal strategy of the center is the
redistribution of large production volumes of parts for the
last time periods in which the production laboriousness of
parts is less than in the initial ones.</p>
      <p>With an increase in the discount rate, the center’s strategy
to redistribute large production volumes of parts for the last
time periods intensifies. This is due to the “cheaper” cost of
the money that the center pays to the agent as a material
reward in remote time periods. With large discount rates, the
effect of deferring the production of parts from the initial
time periods to later ones occurs. It is economically
advantageous for the center to postpone the production of
parts to late time periods, since in this case its total
discounted costs will be minimal.</p>
      <p>250
e
v
ilta n .)s 200
u io cp 150
um tcu (e
lc d m 100
a ro lu
itm p vo 50
p
O</p>
      <p>0
250
e
iltvau ino .)scp210500</p>
      <p>(
um tcu em100
lc d u
a ro l
itm p vo 50
p
O
0
r=0%
r=10%
r=20%
r=30%
r=40%
r=0%
r=5%
r=10%
r=15%
r=0%
r=10%
r=20%
r=30%
r=40%</p>
      <p>Analyzing Fig. 3, we conclude that for the logistic
laboriousness model in the absence of discounting (r=0%),
the optimal trajectory of cumulative production volume is the
logistic curve. The optimal trajectory of the cumulative
production volume consists of two sections: concave and
convex.</p>
      <p>Fig. 4 shows the optimal trajectories of production
volumes for various discount rates r. The optimal strategy of
the center in the absence of discounting (r = 0%) is:
reduction of production volumes for the concave section of
the optimal trajectory of the cumulative production volume
and increase in production volumes for the convex section of
the trajectory. The minimum of production volume
corresponds to the inflection point of the optimal trajectory
of the cumulative production volume.
r=0%
r=5%
r=10%
r=15%
100
n
it .) 80
o
c s
duo (cp 60
r
lap eum40
im lo 20
t v
p
O
0
1
3
5
7
9</p>
      <p>11</p>
      <p>Time period (months)</p>
      <p>When discounting for the logistic model of laboriousness
is taken into account, the effect of postponing the parts
production from the initial time periods to later ones is also
observed. Discounting leads to the appearance of the
cumulative production volume of an additional convex
section in the initial time periods on the optimal trajectory.</p>
      <p>The optimal trajectory of the cumulative production
volume is transformed into a curve of three sections: convex,
concave and convex. The optimal strategy of the center is: on
convex sections of the trajectory to increase production
volumes, on concave sections - to decrease. Inflection points
correspond to extreme values of production volumes.</p>
    </sec>
    <sec id="sec-5">
      <title>IV. CONCLUSION</title>
      <p>The paper considers the dynamic executors incentive task
in discrete time, taking into account the learning curve effect.
To solve the problem, the principle of cost compensation has
been applied, which consists in dividing the original problem
into two tasks: choosing a compensatory incentive system
and solving the optimal control problem with the objective
function equal to the difference between the income of the
center and the labor costs of the agent.</p>
      <p>Using the Bellman dynamic programming method,
numerical solutions of the optimal control problem are
obtained for various laboriousness models. The study of the
impact of the discount rate on the solution of the incentive
problem was conducted.</p>
      <p>Based on a numerical study, the following conclusions
are formulated:</p>
      <p>1. The optimal strategy of the center for the power
based and exponential learning curves models is to
redistribute large production volumes of parts to the last time
periods in which the production laboriousness of parts is less
than in the initial ones.</p>
      <p>2. The consideration of discounting for the power - based
and exponential learning curves models leads to an even
greater redistribution of the production volumes of parts over
the last time periods.</p>
      <p>3. Taking into account the discounting for all the
considered learning curves models leads to the effect of
postponing production from initial periods to later ones.</p>
      <p>4. The optimal trajectory of the cumulative production
volume in the case of the logistic learning curve model is a
curve consisting of several convex and concave sections. The
optimal strategy of the center is to increase production
volumes on convex sections of the trajectory, and to decrease
production volumes on concave sections. Inflection points
correspond to extreme values of production volumes.</p>
      <p>5. Taking into account the discounting for the logistic
learning curve model leads to a redistribution of production
volumes of parts in the middle and recent time periods.</p>
    </sec>
    <sec id="sec-6">
      <title>ACKNOWLEDGMENT</title>
      <p>The reported study was funded by RFBR and Samara
region according to the research project № 17-46-630606.</p>
      <p>Moscow: Izdatelstvo</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>D.A.</given-names>
            <surname>Novikov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.I.</given-names>
            <surname>Smirnov</surname>
          </string-name>
          and
          <string-name>
            <given-names>T.E.</given-names>
            <surname>Shokhina</surname>
          </string-name>
          , “
          <article-title>Mechanisms of Dynamic Active Systems Control</article-title>
          ,” Moscow: IPU RAN,
          <year>2002</year>
          , 124 p.
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>V.A.</given-names>
            <surname>Gorelik</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.A.</given-names>
            <surname>Gorelov</surname>
          </string-name>
          and
          <string-name>
            <given-names>A.F.</given-names>
            <surname>Kononenko</surname>
          </string-name>
          , “
          <source>Analysis of Conflict Situations in Control Systems,” Moscow: Radio i svyaz</source>
          ,
          <year>1991</year>
          , 228 p.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>V.A.</given-names>
            <surname>Gorelik</surname>
          </string-name>
          and
          <string-name>
            <given-names>A.F.</given-names>
            <surname>Kononenko</surname>
          </string-name>
          , “
          <article-title>Game-theoretic Models of Decision Making in Ecological-economic Systems</article-title>
          ,” Moscow: Radio i svyaz,
          <year>1982</year>
          , 144 p.
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          <string-name>
            <surname>Hierarchical</surname>
            <given-names>games</given-names>
          </string-name>
          ,
          <source>” Automation and Remote Control</source>
          , vol.
          <volume>2</volume>
          , pp.
          <fpage>89</fpage>
          -
          <lpage>106</lpage>
          ,
          <year>2015</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>T.</given-names>
            <surname>Basar</surname>
          </string-name>
          and
          <string-name>
            <given-names>G.J.</given-names>
            <surname>Olsder</surname>
          </string-name>
          , “Dynamic Noncooperative Game Theory,” Philadelphia: SIAM,
          <year>1999</year>
          , 519 p.
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>E.</given-names>
            <surname>Dockner</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Jorgensen</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.V.</given-names>
            <surname>Long</surname>
          </string-name>
          and G. Sorger, “Differential Games in Economics and Management Science,” Cambridge: Cambridge University Press,
          <year>2000</year>
          , 382 p.
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>Y.-C.</given-names>
            <surname>Ho</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Luh</surname>
          </string-name>
          and
          <string-name>
            <given-names>R.</given-names>
            <surname>Muralidharan</surname>
          </string-name>
          , “Information Structure, Stackelberg Games, and Incentive Controllability,”
          <source>IEEE Trans. Automat. Control.</source>
          , vol.
          <volume>26</volume>
          , no.
          <issue>2</issue>
          , pp.
          <fpage>454</fpage>
          -
          <lpage>460</lpage>
          ,
          <year>1981</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>G.J.</given-names>
            <surname>Olsder</surname>
          </string-name>
          , “
          <article-title>Phenomena in inverse Stackelberg games</article-title>
          .
          <source>Part</source>
          <volume>1</volume>
          : static problems,” J. Optim.
          <source>Theory Appl.</source>
          , vol.
          <volume>143</volume>
          , no.
          <issue>3</issue>
          , pp.
          <fpage>589</fpage>
          -
          <lpage>600</lpage>
          ,
          <year>2009</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>G.J.</given-names>
            <surname>Olsder</surname>
          </string-name>
          , “
          <article-title>Phenomena in inverse Stackelberg games</article-title>
          .
          <source>Part</source>
          <volume>2</volume>
          : dynamic problems,” J. Optim.
          <source>Theory Appl.</source>
          , vol.
          <volume>143</volume>
          , no.
          <issue>3</issue>
          , pp.
          <fpage>601</fpage>
          -
          <lpage>618</lpage>
          ,
          <year>2009</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>N.</given-names>
            <surname>Groot</surname>
          </string-name>
          , B. De Schutter and
          <string-name>
            <given-names>H.</given-names>
            <surname>Hellendoorn</surname>
          </string-name>
          , “
          <article-title>Reverse Stackelberg games. Part I: basic framework</article-title>
          ,
          <source>” Proc. of the IEEE Int. Conf. on Control Applications</source>
          , pp.
          <fpage>421</fpage>
          -
          <lpage>426</lpage>
          ,
          <year>2012</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>N.</given-names>
            <surname>Groot</surname>
          </string-name>
          , B. De Schutter and
          <string-name>
            <given-names>H.</given-names>
            <surname>Hellendoorn</surname>
          </string-name>
          , “
          <article-title>Reverse Stackelberg games. Part II: results</article-title>
          and open issues,
          <source>” Proc. of the IEEE Int. Conf. on Control Applications</source>
          , pp.
          <fpage>427</fpage>
          -
          <lpage>432</lpage>
          ,
          <year>2012</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>D.B.</given-names>
            <surname>Rokhlin</surname>
          </string-name>
          and
          <string-name>
            <given-names>G.A.</given-names>
            <surname>Ougolnitsky</surname>
          </string-name>
          , “
          <article-title>Stackelberg equilibrium in a dynamic stimulation model with complete information,” Automation and Remote Control</article-title>
          , vol.
          <volume>79</volume>
          , no.
          <issue>4</issue>
          , pp.
          <fpage>701</fpage>
          -
          <lpage>712</lpage>
          ,
          <year>2018</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>A.</given-names>
            <surname>Badiru</surname>
          </string-name>
          , “
          <article-title>Computational survey of univariate and multivariate learning curve models</article-title>
          ,
          <source>” IEEE Transactions on Engineering Management</source>
          , vol.
          <volume>39</volume>
          , no.
          <issue>2</issue>
          , pp.
          <fpage>176</fpage>
          -
          <lpage>188</lpage>
          ,
          <year>1992</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <given-names>L.E.</given-names>
            <surname>Yelle</surname>
          </string-name>
          , “
          <article-title>The learning curve: historical review and comprehensive survey,” Decision Sciences, vol</article-title>
          .
          <volume>10</volume>
          , no.
          <issue>2</issue>
          , pp.
          <fpage>302</fpage>
          -
          <lpage>328</lpage>
          ,
          <year>1979</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Jaber</surname>
          </string-name>
          , “Learning Curves: Theory, Models, and Applications,” Boca Raton: CRC Press,
          <year>2011</year>
          , 476 p.
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          [16]
          <string-name>
            <given-names>R.</given-names>
            <surname>Bellman</surname>
          </string-name>
          , “Dynamic inostrannoy literatury,
          <year>1960</year>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>