Research Article | Open Access

# A Fractional Supervision Game Model of Multiple Stakeholders and Numerical Simulation

**Academic Editor:**Vladimir Turetsky

#### Abstract

Considering the popular use of a certain kind of supervision management problem in many fields, we firstly build an ordinary supervision game model of multiple stakeholders. Secondly, a fractional supervision game model is set up and solved based on the theory of fractional calculus and a predictor-corrector numerical approach. Thirdly, the methods of phase diagram and time series graph were applied to simulate and analyse the dynamic process of the fractional order game model. Results of numerical solutions are given to illustrate our conclusions and referred to the practice.

#### 1. Introduction

There is a typical kind of supervision game model used to describe the harmony, balance, and conflicts among the participating groups from different levels, such as coal mining workers, inspection team, and the government safety supervision department in colliery safety management system, or pollution companies, local environmental protection department and the supervisory department of environmental protection in environmental supervision system, and so forth. The characteristic of this kind of regulatory and supervision issues is that both the inspection and supervision behaviors are done to the same object (the producers). That is, not only does the internal quality inspection department of the enterprise check the standardization of the production process, but also, due to the social externality of the product, the superior administrative department or the community residents supervise the production process. In this way, workers, quality inspection departments and government departments, and other stakeholders will take a long-term dynamic game process. Many researchers have done a lot of study on this kind of problems [1–5] and have obtained plenty of theoretical achievements.

Actually, the supervision game is an extremely complex nonlinear system involving many factors. There are too many characteristics that cannot be described precisely by integral order calculus. Attributed to the fast development of computer technology and numerical simulation, the fractional order calculus has entered a brand new developing stage and has been already successfully applied in many fields of humanism social science researches except physics and engineering [6–9]. Fractional calculus can often make up for the deficiency of integer order calculus. The ever-broadening applications of fractional calculus make us believe that the research progress would be greatly promoted as fractional calculus is introduced into this subject. So in this paper, we try to set up a fractional supervision game model to study the dynamical system of multiple stakeholders.

The rest of the paper is organized as follows. In Section 2, we start with some concepts and assumptions and formulate an ordinary supervision game model of multiple stakeholders, which are necessary for future discussions. Section 3 contains some preliminaries of fractional calculus theory and introduces the Adams-Bashforth-Moulton (ABM) predictor-corrector approach, which are used in the next section. In Section 4, we set up a fractional game model to characterize the gambling process of those multiple stakeholders by the fractional order calculus; then we show an example and carry out numerical simulations. Phase diagram and time series graph are used to numerically simulate and analyse the dynamical process of the gambling system. We finally conclude our paper in Section 5.

#### 2. Assumption and Formulation of the Ordinary Supervision Game Model of Multiple Stakeholders

As we have introduced in the first section, there are many practical examples that can be abstracted into our proposed regulatory problem model. In this paper, we take the case of supervision game in colliery safety management system for example. There are actually two reasons for it: one is for the convenience and pertinence of narration; the other is that safety management is always the difficult and important factor of colliery management in China. As great dangers lie in coal industry, it is very urgent and necessary to decrease accident rate and improve our level of security management.

The gambling system of colliery safety management involves three interest groups from different levels, that is, coal mining workers, inspection team of the coal enterprise, and supervision department of the government, labeled as player I, player II, and player III, respectively. In order to build the game model, some assumptions and notations are given as follows (see Figure 1).

(1) Player I (coal mining workers) has two antagonistic strategies, doing jobs in the* right* way (i.e., operation behaviors complying with referenced regulatory requirements, standards, and codes, denoted as “”) or vice versa (some operation behaviors are not proper, even wrong, denoted as “”). We assume that the probability of choosing the strategy “” is and for choosing the strategy , where The job cost or input is when strategy “” is chosen, or else the job cost is . The parameters and (and some other similar parameters mentioned below) are percentages of unit industrial output.

(2) Player II (inspection team) also has two antagonistic strategies,* inspection and no inspection*, denoted as “” and “,” respectively. The cost of inspection is , and player I (mining workers) would be fined for choosing “.” The probability of choosing the strategy “” is .

(3) The game strategies of player III (supervision department of the government) are* supervision or no supervision*, denoted as “” and “.” If III found some mistakes of I but II takes no inspection, I and II should pay fines of and , respectively. The cost of supervision is , and the probability of choosing supervision strategy “” is .

Considering the actual situation, we assumed that the parameters mentioned above satisfy .

Based on the decision tree shown in Figure 1, we obtain the expection payoff (EP) of the three players. The EP of player I is

The of player II isThe of player III isBy the first-order optimality condition, we getThe mixed strategy Nash equilibrium (MSNE) is given by the solution of three equations above, which is the optimal game strategy choice as we assume that all the participants are fully rational [2, 10]. The MSNE point is not described in detail in this paper because the formula is so complicated that it is hard to know the real meaning of the results, and we prefer to believe that people are bounded rational in the reality.

#### 3. Preliminaries

##### 3.1. Fractional Order Calculus

As irregularities are common in nature and fractal theory has been highly realized in many fields of nonlinear science in nowadays, fractional order calculus is increasingly used as a kind of mathematical tool to describe exactly the essence of the problem. There are several commonly used definitions of fractional order derivative [9, 11–16]. We utilize Caputo fractional order derivative in this paper. In contrast to the Riemann-Liouville fractional derivative, when solving differential equations using Caputo’s definition, it is not necessary to define the fractional order initial conditions. The fractional derivatives of a continuous function are defined bywhere is the smallest integer number, larger than .

##### 3.2. Predictor-Corrector Approach for Fractional Order Calculus

Several numerical methods can be used to solve the fractional differential equations [14, 17–24]. Among them, Adams-Bashforth-Moulton (ABM) predictor-corrector approach has the merits of high precision and good stability, which was developed by Diethelm and other researchers [18, 25, 26] and successfully applied to many systems, such as [12, 14, 19, 27]. We take a fractional differential equation with initial conditions:System (6) is equivalent to following Volterra integral equation:Set the uniform grid () for some integer , and assume that there is a unique solution of the function on some interval ; system (6) can be approximated to the following difference schemes:where

The Adams-Moulton (AM) corrector method is formed by Formulae ((8) and (9)). Then, the predicted value is determined by the Adams-Bashforth (AB) method:whereSo, the Adams-Bashforth-Moulton (ABM) predictor-corrector approach for fractional order calculus has been achieved by formulae (8), (9), (10), and (11). The error of this numerical algorithm is expected as , whereThe two parts of ABM predictor-corrector approach, that is, the predictor method and the corrector method, are independent of each other. Almost any combination of them as we required can become a modifications of the algorithm [14, 18].

#### 4. Fractional Order Game Model and Numerical Simulation

By the traditional game theory, replicator dynamics is an important concept to describe the evolutionary advantage of population in evolutionary game theory [2, 10, 28–30]. The evolutionary advantage of a population is reflected in the proportion of participants; that is, the growth rate of each player ( for player I) equals to the difference between its adaptability and average adaptability, while the adaptability is represented as ER in our study. The most general form of replication equation iswhere, is the vector of the distribution of types in the population, represents the growth rate of type , and the right side of (13) is the difference between the fitness of type and the average population fitness, which reflects the evolutionary advantage of type .

So in the general way, the following replication equations are obtained:

##### 4.1. The Fractional Order Game Model of Colliery Safety Management

While Caputo fractional order derivative is substituted for ordinary integer derivatives, a fractional order game model is set up as follows.

In this part, the ABM predictor-corrector approach will be used to solve the fractional dynamical system (15) with given initial valuesIn order to make a more simple and intuitive qualitative analysis of the character of solutions, the phase diagram [12, 30, 31] and time series graph [12, 32, 33] will be used to depict and analyse the dynamic process of the fractional order game model.

##### 4.2. Numerical Simulation of the Fractional Game Model

According to the model assumptions made in Section 2, the parameters of the game model (15) and initial values of variables are supposed asFormula (17) means if the labor cost of player I’s normal work is 1 unit, his inactive or improper actions only cast 0.2 unit and the costs of player II’s inspection and player III’s supervision are 0.2 and 0.25 unit, respectively, and both player I and player II should pay fine of 2 units when their improper behaviors were found. Initial values of variables (see formula (18)) are the chances of their strategy selection at initial time.

Set , and the fractional orders of differential equations are valued as ; the simulation results of the fractional game model are shown by the phase diagram (see Figure 2) and time series graphs (see Figures 3, 4, and 5).

As the fractional orders are changed to , we obtain the similar phase diagram and time series graph (set for more distinct pictures), as shown in Figures 6, and 7.

Comparing the two groups of graphs (see Figures 2, 3, 4, 5, and Figures 6, 7), we could find out that the fractional order game model possesses similar dynamic features just as the one-order game model (see formula (14)). The trajectory of the dynamical game system is a closed curve around the center. The production safety input of workers, the intensity of the inspection, and the intensity of the government’s supervision affect each other and form an interactive evolutionary process. In general, there is no evolutionary stable strategy (ESS) in this game model, and the dynamical system is not self-controlled. Therefore, it is necessary to restrain the behaviors of inspectors and supervisors.

Under certain conditions (see formulae (17), (18)), the proportion of player I choosing safe operation () maintains at a certain height with high stability (see Figures 3 and 7), while the proportions of player II choosing safety inspection () and of player III choosing the strategy of supervision () remain a certain cyclical volatility (see Figures 4, 5, and 7). These results conform to the facts to a certain extent. It is noticeable that mining workers must operate safely at every moment. However, the inspections or supervision could be done periodically instead of being done all the time. It should be pointed out that the inspections or supervision must maintain a certain degree of level. Otherwise, the strategy of player I maybe changes sharply.

##### 4.3. Some Special Cases of the Fractional Game System

Next, we consider three special cases; that is, there are rules and mandatory provisions on frequency of the inspections or the supervision, or both of them. Although these cases are special, they very likely occur in practice.

###### 4.3.1. Fixing

Firstly, let us fix , which means there are mandatory provisions in the coal enterprises on the frequency of the inspections. Figure 8 shows the phase diagrams of choices of player I and player III ( and ) as we fix to 0.1, 0.3, and 0.4, respectively. It can be concluded from their plots that along with the increase of , the value will more likely approach to 1. When exceeds a certain value, especially, will approach to 1 quickly and directly (see Figure 8(c)).

**(a) Fixing**

**(b) Fixing**

**(c) Fixing**

###### 4.3.2. Fixing

Secondly, let us fix , which means there are rules on the frequency of supervision. Figure 9 shows the phase diagrams of choices of player I and player II ( and ) as we fix to 0.1, 0.3, and 0.5, respectively. Some results of this case are very similar to the above case. It also can be concluded that will approach to 1 with the increase of , and as exceeds some value, will approach to 1 quickly (see Figure 9(c)).

**(a) Fixing**

**(b) Fixing**

**(c) Fixing**

###### 4.3.3. Fixing the Values of and

Thirdly, let us fix , which means there are rules or mandatory provisions on both the frequency of inspections and supervision. Figure 10 shows some time series diagrams of (choices of player I) as we fix to 0.2 and fix to 0.2, 0.25, and 0.3, respectively. It can be concluded that will approach to 1 as long as and exceed some values.

**(a) Fixing**

**(b) Fixing**

**(c) Fixing**

By comparing the three graphics (see Figures 8, 9, and 10), we come to the conclusion that the effect of regulations and supervision would be more significant if there are rules or mandatory provisions on both the frequency of inspections and supervision.

#### 5. Conclusion

In this paper, a fractional supervision game model of multiple stakeholders was built, which could be used to characterize a kind of typical problem in supervision management. The fractional model was solved by the ABM predictor-corrector calculation method. The methods of phase diagram and time series graph were applied to simulate the dynamic process of the fractional order game model. Some results of numerical solutions analysed the behaviors of those interest groups related.

The results of the model show that the fractional order game model is scientific and reasonable and may provide new research methods and research directions for the original research of supervision management systems.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work is supported by National Nature Science Foundation of China (Grant no. 51574157): Research on Multilateral Game and Control Strategies in Coal Mine Safety System.

#### References

- J. Rasmussen, “Risk management in a dynamic society: a modelling problem,”
*Safety Science*, vol. 27, no. 2-3, pp. 183–213, 1997. View at: Publisher Site | Google Scholar - R. Lu, X. Wang, and D. Li, “Extensive form game model of multiple interest groups in colliery safety,”
*Journal of Qufu Normal University*, vol. 39, no. 4, pp. 22–24, 2013. View at: Google Scholar - Y. Fernando, H. H. Ng, and Y. Yusoff, “Activities, motives and external factors influencing food safety management system adoption in Malaysia,”
*Food Control*, vol. 41, no. 1, pp. 69–75, 2014. View at: Publisher Site | Google Scholar - R. Rocha, V. Mollo, and F. Daniellou, “Work debate spaces: a tool for developing a participatory safety management,”
*Applied Ergonomics*, vol. 46, pp. 107–114, 2015. View at: Publisher Site | Google Scholar - O. Akintola and G. Chikoko, “Factors influencing motivation and job satisfaction among supervisors of community health workers in marginalized communities in South Africa,”
*Human Resources for Health*, vol. 14, no. 1, article 54, 2016. View at: Publisher Site | Google Scholar - G. A. Anastassiou,
*Advances on Fractional Inequalities*, SpringerBriefs in Mathematics, Springer, New York, NY, USA, 2011. View at: Publisher Site | MathSciNet - J. A. Tenreiro Machado, “A probabilistic interpretation of the fractional-order differentiation,”
*Fractional Calculus & Applied Analysis*, vol. 6, no. 1, pp. 73–80, 2003. View at: Google Scholar | MathSciNet - A. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo,
*Theory and Applications of Fractional Differential Equations*, vol. 204, North-Holland Publishing, 2006. - I. Petráš, “Chaos in the fractional-order Volta's system: modeling and simulation,”
*Nonlinear Dynamics*, vol. 57, no. 1-2, pp. 157–170, 2009. View at: Publisher Site | Google Scholar - J. Maynard Smith, “The theory of games and the evolution of animal conflicts,”
*Journal of Theoretical Biology*, vol. 47, no. 1, pp. 209–221, 1974. View at: Publisher Site | Google Scholar | MathSciNet - Z. Bai, “On positive solutions of a nonlocal fractional boundary value problem,”
*Nonlinear Analysis. Theory, Methods & Applications*, vol. 72, no. 2, pp. 916–924, 2010. View at: Publisher Site | Google Scholar | MathSciNet - B. Xin, J. Ma, T. Chen, and Y. Liu, “A fractional model of labyrinth chaos and numerical analysis,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 34, no. 1, pp. 32–33, 2011. View at: Publisher Site | Google Scholar - Z. Bai and W. Sun, “Existence and multiplicity of positive solutions for singular fractional boundary value problems,”
*Computers & Mathematics with Applications*, vol. 63, no. 9, pp. 1369–1381, 2012. View at: Publisher Site | Google Scholar | MathSciNet - L. Zhao and W. Deng, “Jacobian-predictor-corrector approach for fractional differential equations,”
*Advances in Computational Mathematics*, vol. 40, no. 1, pp. 137–165, 2014. View at: Publisher Site | Google Scholar | Zentralblatt MATH - H. H. Dong, B. Y. Guo, and B. S. Yin, “Generalized fractional supertrace identity for Hamiltonian structure of NLS–MKdV hierarchy with self-consistent sources,”
*Analysis and Mathematical Physics*, vol. 6, no. 2, pp. 199–209, 2016. View at: Publisher Site | Google Scholar | Zentralblatt MATH - A. Atangana and D. Baleanu, “New fractional derivatives with non-local and non-singular kernel: theory and application to heat transfer model,”
*Thermal Science*, vol. 20, no. 2, pp. 763–769, 2016. View at: Publisher Site | Google Scholar - C. Yin, Y. Cheng, S.-M. Zhong, and Z. Bai, “Fractional-order switching type control law design for adaptive sliding mode technique of 3D fractional-order nonlinear systems,”
*Complexity*, vol. 21, no. 6, pp. 363–373, 2016. View at: Publisher Site | Google Scholar - K. Diethelm, N. J. Ford, and A. D. Freed, “A predictor-corrector approach for the numerical solution of fractional differential equations,”
*Nonlinear Dynamics*, vol. 29, no. 1–4, pp. 3–22, 2002. View at: Publisher Site | Google Scholar - B. Yu and X. Jiang, “A fractional anomalous diffusion model and numerical simulation for sodium ion transport in the intestinal wall,”
*Advances in Mathematical Physics*, vol. 2013, Article ID 479634, 8 pages, 2013. View at: Publisher Site | Google Scholar - Z. Wang, “A numerical method for delayed fractional-order differential equations,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 256071, 7 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - J. Chen, Z. Zeng, and P. Jiang, “Global Mittag-Leffler stability and synchronization of memristor-based fractional-order neural networks,”
*Neural Networks*, vol. 51, pp. 1–8, 2014. View at: Publisher Site | Google Scholar - Y. Cui, “Existence results for singular boundary value problem of nonlinear fractional differential equation,”
*Abstract and Applied Analysis*, vol. 2011, Article ID 605614, 9 pages, 2011. View at: Publisher Site | Google Scholar | MathSciNet - Y. Cui, “Uniqueness of solution for boundary value problems for fractional differential equations,”
*Applied Mathematics Letters*, vol. 51, pp. 48–54, 2016. View at: Publisher Site | Google Scholar | MathSciNet - Z. Bai, S. Zhang, S. Sun, and C. Yin, “Monotone iterative method for fractional differential equations,”
*Electronic Journal of Differential Equations*, vol. 2016, article 6, 2016. View at: Google Scholar - K. Diethelm, N. J. Ford, and A. D. Freed, “Detailed error analysis for a fractional Adams method,”
*Numerical Algorithms*, vol. 36, no. 1, pp. 31–52, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. Diethelm, “Increasing the efficiency of shooting methods for terminal value problems of fractional order,”
*Journal of Computational Physics*, vol. 293, pp. 135–141, 2015. View at: Publisher Site | Google Scholar | MathSciNet - B. Xin and Y. Li, “0-1 test for chaos in a fractional order financial system with investment incentive,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 876298, 10 pages, 2013. View at: Publisher Site | Google Scholar - J. W. Weibull, “Evolutionary game theory,”
*Current Biology*, vol. 9, no. 14, pp. 503–505, 1995. View at: Google Scholar - K. M. Page and M. A. Nowak, “Unifying evolutionary dynamics,”
*Journal of Theoretical Biology*, vol. 219, no. 1, pp. 93–98, 2002. View at: Publisher Site | Google Scholar - R. Lu, X. Wang, and D. Li, “Application of multi-population evolutionary game in colliery safety and supervision,”
*Mathematics in Practice and Theory*, vol. 46, no. 19, pp. 31–36, 2016. View at: Google Scholar - R. Shone,
*Economic Dynamics: Phase Diagrams and Their Economic Application*, Cambridge University Press, Cambridge, UK, 2nd edition, 2002. View at: Publisher Site | MathSciNet - W. Zhao, T. Zhang, X. Meng, and Y. Yang, “Dynamical analysis of a pest management model with saturated growth rate and state dependent impulsive effects,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 204642, 18 pages, 2013. View at: Publisher Site | Google Scholar - X. Meng and L. Zhang, “Evolutionary dynamics in a Lotka-Volterra competition model with impulsive periodic disturbance,”
*Mathematical Methods in the Applied Sciences*, vol. 39, no. 2, pp. 177–188, 2016. View at: Publisher Site | Google Scholar | Zentralblatt MATH

#### Copyright

Copyright © 2017 Rongwu Lu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.