Abstract
Starting from the leading role of public goods suppliers’ supply strategies, the government’s incentives and support to suppliers, and consumers’ influence on public opinion and other factors, the matrix game models between public goods suppliers and between the suppliers and the government are established, respectively. Considering the delay of government and supplier strategies, establish two delay differential equations between the two suppliers and between the supplier and the government, obtain the evolutionary game model of public goods supply, and solve the equilibrium point of the evolutionary game model with or without delay and its stability. Then, through model parameter analysis, verify the rationality of public goods supplier’s supply and government incentive strategy. Finally, through numerical simulation, it is verified that the delay does not affect the stability of the equilibrium point of the game party but changes the rate at which the game party’s strategy reaches a steady state, conducive to rapid, scientific, and reasonable decision-making between the suppliers of public goods and between the government and suppliers.
1. Introduction
In 1954, American economist Paul Samuelson pointed out in the Pure Theory of Public Expenditure that pure public goods are nonexclusive and noncompetitive. Public goods cannot be simply balanced through market allocation. The supply of public goods is the result of different stakeholders’ choice over a part of the social resource allocation strategy. Namely, public goods are a kind of product, where “public” means those products and services that can be shared by many people at the same time [1], such as science, education, culture, health, public facilities, environmental protection, diplomacy, national defense, etc. Public goods can be provided by the public sector, the private sector, or a combination of the two. The two have different orientations: public goods provided by the public sector focus on social benefits, while the private sector expects to obtain economic benefits.
It is precisely because of the different orientations of the two sectors that provide supply strategies, the social state will fall into a game dilemma. For example, when the government gives insufficient incentives to suppliers, the prisoner’s dilemma of insufficient public goods supply will arise. If the government encourages excessively, suppliers’ supply problems will arise again, like pig’s payoffs or chicken game dilemma. Literature [1–3] makes a comprehensive analysis. Literature [4, 5] uses population dynamics to study the two-predator-prey population equation, and it analyzes the application of delay differential equations in the public goods bidding game mechanism. Literature [6, 7] establishes a pair of evolutionary game between peasant entrepreneurs and village public goods providers and between the entrepreneurs and the government. Literature [8] analyzes the public goods game theory model of collective risk and adds growth factors to analyze the evolution of replication dynamics. Literature [9] uses the evolutionary game theory, introduces compensation coefficients based on the evolution of the government and social capital, builds a tripartite evolutionary game model for the government, enterprises, and the public, and analyzes the evolutionary equilibrium of the tripartite under different strategies. Literature [10–15] discusses the evolutionary game strategy equilibrium and other issues regarding local government itself and the coordinated control of air pollution, industrial buildings, campus online loan supervision, public transportation, information sharing, and enterprise knowledge sharing from the perspective of evolutionary game. Literature [16] aims at the evolutionary trend of the time-lagged evolutionary game of switching topology enterprise innovation, constructs enterprise cooperative innovation game network through enterprise nodeization, establishes the time-lagged evolutionary game model of switching topology enterprise innovation, and verifies the feasibility of equilibrium with simulation examples. Literature [17, 18] applies dynamics-related knowledge, the steady state of the players’ strategies in the game, and analyzes the impact of delay on decision-making through numerical simulation.
It can be seen from the existing literature that most scholars do research on the supply of public goods based on the evolutionary game model, however, there are few articles that really consider supply and encourage delay. This article starts with the factors that affect the game party’s income and establishes a matrix game model between public goods suppliers and between suppliers and the government. Starting from the delay of the game party’s strategy, this article establishes two delay differential equations between the suppliers and between the suppliers and the government and attempts to solve the equilibrium point of evolutionary game model with or without a delay and its stability. Then, through model parameter analysis, verify the rationality of public goods supplier’s supply and government incentive strategy. Finally, through numerical simulation, it is verified that the delay does not affect the stability of the equilibrium point of the game party but changes the rate at which the game party’s strategy reaches a steady state, conducive to rapid, scientific, and reasonable decision-making between the suppliers of public goods, such as coal mine, water, electricity, and gas, and between the government and suppliers.
2. Analysis of Game Model between Two Suppliers of Public Goods
2.1. Establish a Matrix Game Model between Public Goods Suppliers
According to the introduction analysis and literature [6], the assumptions made to establish a matrix game model between public goods suppliers, as shown in Table 1, are as follows:
Assumption 1. when supplier A and B choose to supply public goods, the cost they need to pay is >0, >0. The tangible benefits that public goods suppliers choose to bring to each other are as follows: >0, >0. At the same time, the hidden market benefits of the supplier who chooses the supply strategy are >0, >0.
Assumption 2. the government’s incentives for suppliers to choose supply strategies, such as rewards, subsidies, and preferential policies, are as follows: >0. When two suppliers chooses the supply strategy at the same time, the government’s incentives are evenly distributed, that is > .
Assumption 3. when only one supplier chooses the supply strategy, the costs and benefits remain the same, however, the overall revenue is less than 0, and the nonsuppliers’ product costs and revenues are 0. However, the intangible benefits of them losing market or public support because of nonsupply are far greater than the tangible benefits that they can obtain from other suppliers.
Assumption 4. when all suppliers do not supply, the supplier does not incur supply costs and benefits, i.e., the costs and benefits are 0.
Assumption 5. the supply probabilities of public goods provided by supplier A and B are , respectively, where 0≤ ≤1, and 0≤ ≤1.
Based on the matrix game model, the expected return of supplier A is as follows:The replication dynamic equation of supplier A 18,19is obtained as follows:Based on the matrix game model, the expected return of supplier B is as follows:The replication dynamic equation of supplier B is obtained as follows:The evolutionary game model between suppliers is obtained from (1) and (2).
2.2. Establish Delay Differential Equations
Since the supply strategy of the public goods supplier directly affects its own revenue and the expected revenue obtained in the current period is related to the game strategy of the previous competitors. The supplier will predict the current strategy based on the competitor’s previous game strategy, i.e., the delay of the competitor’s game strategy needs to be considered.

Here, suppose that the delay of public goods supplier B to A is and the delay of A to B is . Suppliers A and B evolve and copy the dynamic (6) to obtain the delay differential equation.
2.3. Analysis on the Stability of Equilibrium Points between Public Goods Providers
2.3.1. Stability of Equation (6) Equilibrium Points
(1)Equation (6) equilibrium points and Jacobian matrix Equation (6) equilibrium points are as follows: O(0,0), A (0,1), B (1,0), C (1,1), D (). A Jacobian matrix is obtained through (6).(2)Stability of Equation (6) From the Jacobian matrix trace and determinant symbol of equation (6), the stability results of equation (6) equilibrium points are shown in Table 2. According to evolutionary game theory and as can be seen from Table 2, among the five equilibrium points, only the Nash equilibrium point O (0,0) and B (1,1) are the ESS points, and C (0,1) and A (1,0) are the unstable points, while D D is a saddle point. According to Table 3, the replication dynamic relationship phase of the evolutionary game party can be obtained as shown in Figure 1.(3)Parameter Analysis Figure 1 shows that the equilibrium result depends on the initial state of the game. When the initial state is in the region ABCDA, the system converges to point B, i.e., evolutionary stability point (ESS). When the initial state is in the region OADCO, the system converges to point O, i.e., evolutionary stability point (ESS). Currently, there is an overall lack of public goods supply. Hence, it is clear that the (supply, supply) strategy combination is the best result. Consequently, to continuously expand the ABCDA area and shrink the OADCO area, the saddle point D must be moved downward. Namely, becomes smaller. Hence, it is necessary for public goods supplier to increase the probability of the supply strategy. Suppliers need to reduce the supply cost of public goods, increase the supply of public goods, or reduce the loss when public goods are not supplied.2.3.2. Stability of Equation (7) Equilibrium Points
(1)Equation (9) equilibrium points From replication dynamic (9), it is easy to find out that its equilibrium points are the same with (6), i.e., O(0,0), A(0,1), B(1,0), C(1,1), D . Since (9) is a replication dynamic equation of the mutual influence with two time delays, the sum of the two time delays, and will be used as a parameter in the following paragraphs to analyze the stability of the equilibrium points.(2)Analysis on the stability of equation (9) For the convenience of calculation, the matrix game model between the two suppliers is replaced, where A = + - + + , B = + - + + , C = -, D = - + + , E = - + + , F = -, M = B-- = , and N = = . Then, equation (9) is transformed into the following:
The equilibrium points of equation 10 are as follows: O(0,0), A(0, 1), B(1,0), C(1,1), and D().
Now, discuss the stability of equation (10) at the equilibrium point C (1,1).
At C (1,1), equation (10) is transformed into the following: R = −1, S = −.
A Jacobian matrix is obtained as follows:
The stability of the equilibrium point B (1,1) is studied through the stability of equation (11) at the origin (zero solution), i.e., R () = 0, S () = 0.
Hence,
The linear approximation of (11) is obtained as follows:
As A = + - + + , B = + - + + , C = -, D = - + + , E = - + + , F = -, M = B-- = , and N = = , substitute equation (14) to get the approximate equation of the original game model.
Assume that = = .
Then, equation (14) is transformed into the following:
Hence, the characteristic equation of (16) is as follows:where =
When = = 0, equation (17) is transformed into the following:
Therefore, the stability of equation (10) at B (1,1) is obtained by discussing the roots of characteristic equation (18). Assuming Δ ≠ 0, regarding the roots of equation (17), we have the following:
Theorem 1. When = []( = 0, 1, 2, ⋯⋯), equation (10) has a pair of pure imaginary root = i.
Prove: let = η+iξ, and substituting into equation (10), get:
. At the same time, using Euler’s formula, there are,
Separate the real part and the imaginary part of (19), and let both parts be 0.
So, = 2+ 2 = 2, i.e.,
Hence, equation (22) has a positive real solution, i.e., = i.
From (21), get the following:
= -
Hence,
i.e., ( = 0,1,2, ⋯), equation (18) has a pair of pure imaginary roots = i.
The conclusion is established.
Theorem 2. The real part of the solution of equation (18) satisfies 0, where = []( = 0,1,2, ⋯).
Prove: derive on both sides of equation - = 0 and get the following:
Let 2.
so, Where, R = ), and S = .
Therefore, the conclusion is established. (27) at B(1,1) is the evolutionary stable point of ESS.
Now, discuss the stability of (27) at equilibrium point D , where M = , N = , R = −, and S = −.
Then, (27) is transformed at ().
A Jacobian matrix is obtained as follows:
Similarly, study the stability of equilibrium point D by discussing the stability of equation (2)–(5) at origin (zero solution), i.e., R = 0 and S = 0.
So,
The linear approximation of equation (2)–(5) is obtained as follows:
As A = + - + + , B = + - + + , C = -, D = - + + , E = - + + , F = -, M = B-- = , and N = = , substitute equation (2)–(7) to get the approximate equation of the original game model.
Let Then, (16) is transformed into the following:
Hence, the characteristic equation of (32) is as follows:where >0, and = .
When = = 0, equation (33) is transformed into the following:
Therefore, the stability of (7) is obtained by discussing the roots of characteristic (33). Let Δ ≠ 0, regarding the roots of (33), we conclude that Theorem 1 and Theorem 2 are established.
Consequently, according to literature [19–23], there is an increasing sequence of real numbers ( = 1,2, ⋯), and all the solutions of (33) on τ∊ (0, ) have negative real parts. When = , there is a pair of pure imaginary solutions, and the others are solutions with strictly negative real parts. When > there is at least one solution with strictly positive real parts, i.e., when = , the origin of equation (18) is at the Hopf branch, which is asymptotically stable on ∊ (0, ) and unstable on > .
The above is the analysis of equilibrium points B(1,1) and D(). Other equilibrium points can be studied similarly, and it can be concluded that the stability of the equilibrium points of model (2-3) and model (2-4) are the same, except that there is a delay in the model (2-4). Moreover, the rate of equilibrium points evolving to a steady state is different from that of model (2-3), which provides a reference for public goods suppliers to provide rapid and reasonable supply strategies.
3. Game Analysis between Public Goods Suppliers and Government
3.1. Establish a Matrix Game Model between Public Goods Suppliers
To increase the supply of public goods, the government encourages suppliers to choose supply strategies through incentive policies. According to the introduction analysis and literature [6], the following assumptions are made to establish an evolutionary game model between public goods suppliers and the government, as shown in Table 4.
Assumption 1. when the government adopts incentive strategies, the social benefits, such as political achievements and public support, are as follows: >0. The cost paid by the government to encourage publicity and mobilization is as follows: >0.
Assumption 2. when suppliers adopt the supply strategy, the income obtained is divided into tangible income, >0, and invisible income, >0. At the same time, the government incentive is >0, however, suppliers have to bear the supply cost, i.e., >0.
Assumption 3. when the government adopts incentives and suppliers provide supplies, the total government revenue is as follows: −− <0. The total revenue of suppliers is as follows: + − + >0.
Assumption 4. when the government adopts a nonincentive strategy and suppliers provide supplies, the government's revenue and cost are 0. However, losing public the support and other social benefits is >0, i.e., the total government revenue is −, and the total revenue of suppliers is divided into − + <0.
Assumption 5. when suppliers choose not to supply and the government chooses an incentive strategy, suppliers will lose the opportunity revenue: >0, and the government's revenue is − <0.
Assumptions 6. when suppliers choose not to supply and the government chooses not to incentivize, suppliers will lose the opportunity revenue, i.e., >0. The government loses public support and other social benefits, i.e., .
Assumption 7. the probability of the public goods supplier choosing the supply strategy is , and the probability of the government providing incentives is . Among them, (0≤ ≤1) and ().
According to the assumptions, the model has pure strategy Nash equilibrium points (1, 1) and (0, 0), as well as mixed strategy Nash equilibrium points.Based on the payment matrix, suppliers’ expected return is as follows:Suppliers’ replication dynamic equation is obtained as follows:Based on the payment matrix, the government’s expected return is as follows:Government’ replication dynamic equation is obtained as follows:The evolutionary game model is derived from the replication dynamic equation of suppliers and the government:
3.2. Establish Two Delay Differential Equations
When the government provides incentive strategies, such as tax incentives and financial subsidies, to public goods suppliers, to make the incentive strategy more reasonable and achieve the expected benefits of social and economic development, the government needs to fully consider suppliers’ supply delay. At the same time, suppliers must consider their own development, and the delay of incentives given by the government should also be fully taken into account when selecting public goods supply strategies.
Here, suppose that the delay of public goods suppliers to the government’ incentives is , and the government’s supply delay to suppliers is . Based on the evolutionary game model (3-1) between suppliers and the government, a two delay differential equation is obtained as follows:
3.3. Analysis on the Stability of Equilibrium Points between Public Goods Suppliers and Government
3.3.1. Analysis on the Stability of Equations (42) and (43)
(1)Equilibrium points of (40) and (41) Equation (40) has the following equilibrium points: O(0,0), B(1,0), D(1,1), F(0,1), and P(). Meanwhile, it is known that equations (2)-(3) have the same equilibrium points, i.e., O(0,0), B(1,0), D(1,1), F(0,1), and P . Similarly, (41) is a replication dynamic equation of the mutual influence with two time delays. The sum of the two time delays, and will be used as parameters in the following paragraphs to analyze the stability of the equilibrium points.(2)Firstly, analyze the stability of the model (42) at P(). For the convenience of calculation, the game model (42) is replaced as shown, where , .
Let ф = − = -, = - = -. Then, model (44) is transformed into the following:
A Jacobian matrix is obtained as follows:
The stability of the equilibrium point P() is studied through the stability of equation (3)-(4) at the origin (zero solution).
The linear approximation of (45) is obtained as follows:
Substitute M = ,N = equation (49) to get the approximate equation of the original game model.
Let = <0, = >0.
Then, (47) is transformed into the following:
Hence, the characteristic equation of (48) is as follows: = <0, so the zero solution of equation (43) is unstable, i.e., the zero solution of equation (41) is unstable. At the same time, from the Jacobian matrix trace and determinant symbol of equation (41), the stability results of equation (41) equilibrium points are shown in Table 5.
For (41), it is known that its equilibrium points are the same with those of (40). Furthermore, (40 and 41) solve the same Jacobian matrix. Hence, the stability of (41) equilibrium points is the same, as shown in Table 5.
The equilibrium points of (40) are as follows: O(0,0), B(1,0), D(1,1), F(0,1), P(, ), and the replication dynamic phase of the evolutionary game is illustrated by Figure 2.

According to Table 5 and Figure 2,(1)When the initial state of the game system is in the OAPGO area, saddle point O is stable along the OF direction, i.e., when the suppliers’ supply probability is < and continues to decrease, the government will provide more incentives to suppliers.(2)When the initial state of the game system is in the ABCPA area, saddle point B is stable along the BO direction, i.e., the probability that the government provides incentives is < and PG continues to decrease, the intensity of suppliers to provide supply will reduce.(3)When the initial state of the game system is in the PCDEP area, saddle point D is stable along the DB direction, i.e., when the supply probability of public goods suppliers is > and PE continues to increase, the intensity of government to offer incentives will reduce.(4)When the initial state of the game system is in the PEFGP area, the saddle point F is stable along the FD direction, i.e., the probability that the government provides incentives is > and PG continues to increase, the intensity of suppliers to provide supply will increase.
3.3.2. Parameter Analysis
The following parts will use the replication dynamic phase of evolutionary game to analyze the parameters, illustrating the practical significance of public goods suppliers’ selection of supply strategies and government’s incentive strategies.(1)From the government’s perspective, the most ideal evolutionary trend state is in the CDEP area, namely, public goods suppliers continue to increase the probability of supply, and the government continues to reduce support. At this time, suppliers are more proactive and enthusiastic, which may form an endogenous mechanism for suppliers to provide supply. Furthermore, it means the government expects that the central point P will gradually move downward, the value of will remain unchanged or decrease, and the value of will also decrease, thereby expanding the CDEP area. On the premise of < <1,0< <1, with all the parameters above 0, and decrease, and increase, making the value of bigger and smaller. Center point was moved downward. Suppliers reduce the supply cost , the government increases incentives, and suppliers increase supply revenue, which can fuel the enthusiasm of suppliers. The government can gradually reduce the incentives to suppliers based on the delay of revenue .(2)From suppliers’ perspective, the most ideal evolution trend state is in the EFGP area, i.e., while suppliers continue to increase the probability of supply, they can get incentives from the government. Therefore, public goods suppliers’ expectation P gradually moves to the lower right to continuously expand the EFGP area. On the premise of 0< <1 and 0< <1, with all the parameters above 0, and decrease, while and increase, making the value of bigger and smaller. Center point P moves to the lower right, i.e., reducing suppliers’ supply cost and the government cost , increasing suppliers’ revenue, or enhancing public support, for the government can mobilize the supply enthusiasm of public goods suppliers. At the same time, according to the delay of revenue , increase the probability of government incentives.(3)For OAPG and ABCP area, its space will become smaller, which is a state that neither the government nor public goods suppliers want. It will not be discussed in the article.
4. Numerical Simulation
4.1. Numerical Simulation of Equations (2)-(3)
= 2.5, = 1.5, = 4.0, = 4.5, = 1.5, = 4, = , and = 3. Let A = + - + + , B = + - + + , C = -, D = - + + , E = - + + , and F = -.
A = 6.5, B = 5, C = 1.5, D = 4, E = 5.5, F = 1.5, G = 0, and H = 0 are the initial values of the parameter, and the time length is 10. Systems with different supply probabilities are obtained (Figure 3(a), 3(b)).

(a)

(b)
4.2. Numerical Simulation of Equations (2)–(4)
A = 6.5, B = 5, C = 1.5, D = 4, E = 5.5, F = 1.5, G = 0, and H = 0 are still the initial values of the parameter, and the time length is 10. Take = 0.2 and = 0.2, respectively (as shown in Figures 4(c), 4(d)), and when = 0.2 and = 0.5 (as shown in figures 4(e), 4(f)), the two time-delay systems are as follows:

(a)

(b)
4.3. Numerical Simulation of Equation (50)
and = d = 1 are the initial values of the parameter, and the time length is 10 (Table 6). When = = 0.2 (Figure 5(a)), = = 0.4 (Figure 5(b)), when = 0.5, = 0.2 (Figure 5(c)), and when = 0.2 and = 0.8 (Figure 5(d)), it is a deterministic system.

(a)

(b)

(c)

(d)
4.4. Numerical Simulation of Equation (50)
and = d = 1 are the initial values of the parameter, and the time length is 10. When = 0.2 and = 0.2, it is a system with different time delays (Figure 6(a), 6(b)).

(a)

(b)
Based on numerical simulation figures,(1)From the game process, regardless of whether there is a delay between suppliers, it does not affect the final stable state of the strategy [6] (see Figures 3 and 4). Finally, and are asymptotically stable at 1, which is the same as the previous conclusion. However, compared with deterministic systems, time delay systems have different rates from decision-making to a stable state of suppliers because of delays. Both are slower, and the main reason is that delays lead to a decrease in the number of times players learn to adjust strategies during the game. However, when the suppliers are in a delay, the longer the delay, the faster the evolution rate. See Figure 3(a), Figure 4(a), 4(b), Figures 3(b), and Figure 4(e), 4(f). It shows that in the game, because of excessive consideration of the competitor’s previous strategy, although the frequency of the continuous adjustment of the strategy is less, the time consumed to reach a stable strategy will be shortened because the strategy selection is more conservative and prudent.(2)In the game between public goods suppliers and the government, the cyclical evolutionary relationship between and is shown in Figure 5, however, it cannot reach a stable state. However, because of delay, the supply and incentive strategies between suppliers and the government have an uncertain probability of growth, as shown in Figure 6. It shows that as for the supply of public goods, if there is a delay between suppliers and the government, both sides will act cautiously, which can largely avoid loses caused by improper supply and incentive strategies.
5. Research Conclusions
5.1. Conclusion
The equilibrium results of the evolutionary game between public goods suppliers depend on the initial state of the game. In the absence of overall supply, the (supply, supply) strategy combination is the game result expected by suppliers. Therefore, suppliers should strive to reduce the supply cost of public goods, increase the supply revenue of public goods, and/or increase the loss of nonsupply, which can increase the probability of supplying public goods.
5.2. Conclusion
(1)Considering the evolutionary trend of the equation from the perspective of the government, as public goods suppliers continue to increase the probability of supply, and the government improves publicity, incentives, subsidies, and the supply revenue of public goods, the enthusiasm of public goods suppliers will be mobilized. Consequently, the government can gradually reduce incentives and other measures for public goods suppliers.(2)Considering the evolutionary game from the perspective of public goods suppliers, suppliers can obtain government incentives and others while continuously increasing the probability of public goods supply. Increasing the revenue of public goods suppliers, enhancing public support for the government, and raising the probability of government incentives can mobilize the enthusiasm of public goods suppliers to supply.5.3. Conclusion
The government and suppliers should view delay gains more rationally. Only when the government handles the supply incentive mechanism of public goods suppliers under the principle of satisfying public needs, can it handle the structural equilibrium of public goods supply. Otherwise, because of improper government incentives and other measures, there will be insufficient supply of public goods, resulting in a prisoner’s dilemma in the supply of public goods. At the same time, it is necessary to prevent the failure of incentive policies because of excessive government incentives and guard against the problem of pig’s payoffs or chicken game dilemma. In the government’s incentive regulation, how to obtain more time-sensitive data is of great significance to the state of public goods supply.
It can be drawn from the above conclusions that in the game of public goods supply, the government focuses on social benefits, while suppliers focus on obtaining economic benefits. It is of great significance for the government to encourage suppliers to provide public goods to improve the government’s social service level, enhance the city’s brand and humanistic literacy, and foster a beautiful living environment.
6. Further Studies
Based on the existing research on the evolutionary game delay of public goods supply, the next step will be to consider the game delay between public goods multisuppliers and the government and the multigroup game model for research.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research study was sponsored by the Project Funding for Science and Technology Program of Guizhou Province, the project number is [2019]1067; the Research Project Funding for Talent Introduction of Guizhou University, the project number is [2017]59; and the Teaching Reform Project of Guizhou Province, the project number is 2019008.