Measurement on the Complexity Entropy of Dynamic Game Models for Innovative Enterprises under Two Kinds of Government Subsidies

Abstract: In this paper, setting the high-tech industry as the background, we build a dynamic duopoly game model in two cases with different government subsidies based on the innovation inputs and outputs, respectively. We analyze the equilibrium solution and stability conditions of the system, and study the dynamic evolution of the system under the conditions of different system parameters by the numerical simulation method. The simulation results show that both innovation subsidy policies have positive effects on firms’ innovation activities. Besides, improving the level of innovation can encourage firms to innovate. It also shows that an exaggerated adjusting speed of innovation outputs may cause complicated dynamic phenomena such as bifurcation and chaos, which means that the system has relatively higher entropy than that in a stable state. The degree of the government innovation subsidies is also shown to impact the stability and entropy of the system.


Introduction
Nowadays, with the rapid improvement of the global scientific and technological level, there is a growing need for innovation to be considered in many aspects of enterprises, such as production, operation and management.The use of innovation changes the complexity characteristics as well as the entropy of the economic systems, which makes the study of innovation activities, government policies as well as dynamic evolution and entropy of systems in this new context intriguing.Schumpeter was the first to integrate the concept of innovation into economic research.He proposed that innovation builds a new production function and introduces a new combination of production elements and production conditions into the production system.After him, many scholars have carried out in-depth studies in the field of innovation.González and Pazó [1] took the price and product innovation investment as decision variables in a multiplayer game model, and analyzed the optimal innovative input decisions of enterprises in different situations.Hasnas et al. [2] considered the effect of innovation spillover; they built a game model in the context of open-innovation and found that whoever gets more benefits from the innovation spillover can obtain higher profits than its competitors.Toivanen et al. [3] found that innovation plays a positive role in improving enterprises' market value.Li and Ma [4] considered both competition and cooperation in their model and studied the impact of research and development in technology on the production costs.Fontana and Nesta [5] studied the effect of product innovation on enterprise's survival and found that successful product innovation can improve the survival rates.Lambertini and Mantovani [6] investigated the timing of adoption of product and process innovation using a differential game where firms can invest in both activities.The studies above provide some references for establishing enterprise games and determining the optimal innovation inputs or outputs.These researchers also proved that innovation activities would have a positive influence for enterprises, so studying the innovation activities of enterprises is of great significance in practice.
In the long run, an innovation strategy will improve the competitiveness of enterprises, promote economic development and enhance countries' strength.However, innovation has the characteristics of high investment, high risk and long cycle, so there is no doubt it represents an enormous challenge for enterprises, especially for the small or medium enterprises, to be innovative, thus their enthusiasm for applying innovation strategies is greatly reduced, so governments have implemented a series of policies in order to improve enterprises' enthusiasm for innovation, such as direct government funding, tax incentives, the protection of intellectual property rights, government procurement policies and so on.These can all reduce the costs of innovation, improve enterprises' expected profits and promote innovation activities.In recent years, scholars did a series of studies on the relationships between government policies and innovation mechanisms.Park [7] analyzed the effectiveness of government subsidies and their impact on productivity in R&D activities through an empirical study of 6900 government-sponsored projects.Catozzella and Vivarelli [8] discussed the impact of innovation subsidies on both enterprises' input and output, and the results show that innovation subsidies have negative effects on innovation output.Huang et al. [9] studied the effectiveness of government subsidies which cultivate innovation ability, by applying a stochastic frontier analysis against the background of Chinese manufacturing.Fölster [10] discussed the incentive effect of subsidies on R&D in the presence of enterprise cooperation.Un and Montoro-Sanchez [11] did empirical research on the service industry and found that public funds provide necessary innovation resources for enterprises, thus they can effectively improve the enthusiasm for innovation.Kang and Park [12] studied the South Korean biotech industry, where they found that government R&D subsidies and the cooperative relations between enterprises and research institutions both have positive impacts on enterprises' innovation activities.Guo et al. [13] studied the effects of innovation funds on small and medium technology-based firms.Hinloopen [14,15] compared two kinds of policies in promoting R&D activity.One policy allows firms to cooperate in R&D and the other provides R&D subsidies, and the study revealed that providing R&D subsidies was more effective in promoting R&D activity.Kleer [16] held the view that government subsidy policies for R&D could provide an effective signal for private investors, and then they could make the right investment decisions.The study showed that the subsidy accompanied by a quality signal could lead to better selected private investments.Lerner [17] did an empirical study on the Small Business Innovation Research program in the United States.His study showed that the firms who got awards from the government grew significantly faster than matched firms, which means the awards had a positive influence.Meuleman and De Maeseneire [18] examined the impact of R&D subsidies on the access to external equity and debt financing of small firms by using a unique Belgian dataset.They found the R&D subsidies brought a positive effect both on SME quality and access to long-term debt.Takalo and Tanayama [19] considered the presence of financial constraints, they established a theoretical model and analyzed the effect of a government R&D subsidy program.The results show that the R&D subsidy policies may be welfare-improving under certain conditions.The studies above mainly discuss the impact of certain kinds of government subsidies on innovation activities.However, it's rare for such researches to cast light on the effect of different kinds of subsidy forms on enterprise innovation.
In recent years, many scholars have adopted dynamic oligopoly games in research on economic and management, and analyzed the dynamic evolution of systems under the conditions of repeated games.At the same time, some scholars have taken entropy theory and chaos theory into the field of economic and management, which can effectively analyze the stability of systems.Zapart [20] used entropy theory to predict financial time series and found that there is a weak trading advantage in financial forecasts of foreign exchange currency futures initiated in low entropy regions while it's extremely difficult to predict time series in high entropy regions.Han et al. [21] established a duopoly game model with double delays in the hydropower market, and they analyzed the influence of time delay parameter on the entropy and stability of the system.In this paper, we propose two different forms of government subsidies for innovation input and innovation output, respectively, and we establish a dynamic game model between two enterprises with bounded rationality.In this process, we analyze the complexity of the model and discuss the effects of two different forms of government subsidies on stimulating the innovation activities.Besides, we analyze how the innovation decisions' adjusting speed, the degree of government subsidies and the beneficial coefficient of innovation influence the equilibrium, entropy and stability of the system.

Model Analysis
Assuming that there are two firms in a high-tech industry (such as electronic information and the Internet), producing two similar products, where one has some substitutability compared to the other.The firms both carry out innovation activities in their production process.Referring to A-J's classical two-stage study [22], we divide the game into two stages, one is the innovation investment decision making stage and the other is the product pricing stage.We solve our game model using the classic backward induction method.We assume that p 1 (t) , p 2 (t) are respectively the prices of the products produced by the two firms at discrete periods t (t = 0, 1, 2, . ..), accordingly, the quantity demanded are Q 1 (t) , Q 2 (t).Assume M 1 (t) , M 2 (t) are their investment funds for innovation activities.It is obvious that a firm's innovation activities will increase its own product demand and threaten its rival, so Q I 1 (t) , Q I 2 (t) represent the increased demand resulting from their own innovation.We call them the innovation output in this paper.Because of the substitutability of the two products, assume θQ I 1 (t) , θQ I 2 (t) are the decreased demand due to a rival's innovation activity, so in period t, the change in demand resulting from the innovation is: Considering the marginal diminishing effect of the innovation benefits, we assume that the investment funds in the innovation activity and the innovation output have a quadratic relationship, thus: where λ (λ ≥ 0) is the innovation input parameter; the greater λ is, the more it will cost to get a unit of the innovation output, θ(0 < θ ≤ 1) represents the degree of substitution of the products, the greater θ is, the higher degree of the substitution will be, the products are complete substitutes when θ = 1.Equation (2) means that , which reflects the marginal diminishing effect of the inputs that conforms to the actual situation.Equation (2) assumes that the product demand has a linear relation with the price, considering the increase of product demand resulted from the product innovation, we can get the demand function of the two firms as: where a(a > 0) represents the market capacity, b(b > 0) is the price elasticity.Equation (3) means the government supports firms' innovation by giving a certain proportion of financial subsidies according to their investment funds in the innovation activity.The government, as the policy maker, maintains the same subsidy rates in different periods, and we use ω(0 < ω < 1) to represent the subsidy rates, so the government subsidies are ωλ(Q I i (t)) 2 , and the innovation expenses to be paid by the firms themselves are (1 Assume c 1 , c 2 are the costs of production, we can get the profit function in period t: Substituting Equation ( 3) into (4), then the objective profit function is:

Model Solving
As stated earlier, the game is divided into two stages, one is the innovation investment decision stage and the other is the product pricing stage.Because the investment funds in the innovation activity and the innovation output have a quadratic relationship, a change in the innovation output can lead to a change in the innovation input in the same direction, so in this game, the firms determine their innovation outputs Q I i (t) in the first stage, then, in the price decision stage they determine their product prices p i (t).We solve this game model by using the classic backward induction method.Firstly, two firms determine their prices by maximizing their profits in the second stage.Then, in the first stage we assume that both players are bounded rational [23], when the enterprises are making decisions at period t, because of the limitation of obtaining information, neither party of the competitors can learn about the exact innovation output decisions that its rival makes, so the decisions they make are based on the incomplete information, and the decisions are not made by maximizing the profits but by maximizing the marginal profit effect instead.

The Stage of Price Decision
Using the profits function of two firms in (5): Let ∂p 2 (t) = 0, then we can get the equilibrium price in period t: Putting Equation ( 7) into (6), we get: where Entropy 2016, 18, 424 5 of 20

The Stage of Innovation Output Decision
Two firms are both bounded rational, so their decisions for the next period are based on the innovation output in this period and the marginal profit effect is: where α, β(α > 0, β > 0) are respectively the adjustment parameter of innovation decisions of firm 1 and firm 2.

Analysis of the Equilibrium Point and the Stability of the System
Let Q I i (t + 1) = Q I i (t), we have assumed that both firms carry out innovation activity, so exclude Q I i (t) = 0, then the unique equilibrium point of the system is: where: is the unique equilibrium point.In order to analyze its stability, we will calculate the Jacobi matrix of the system, which is: and the characteristic equation of J Q I is: Entropy 2016, 18, 424 6 of 20 According to the Jury stability criterion, we can get the stability conditions as below: As there are too many parameters in the system, it is too difficult to analyze the system directly.In order to analyze the stability of the system more clearly, we do a numerical simulation by considering actual competition, and we will show the parameter range according to the Jury stability criterion in the next section.

The Influence of the Decision Parameters on the Stability and the Entropy of the System
In order to study the impact of adjustment parameters on the system evolution, we set β = 0.8, and we can describe the dynamic behavior and entropy of the system with varying α.Figures 1-3 show two firms' innovation outputs, product prices and the profits as α changes.We can find that when α takes a small value, the equilibrium point is stable, but with α increasing, bifurcation occurs in the system and then the system even falls into chaos.
We know that entropy measures the chaotic degree of the system, so it is not difficult to find that the entropy of the system is increasing with α.The entropy of the system shows the probability of some particular information, and it will increase with the increasing of the uncertainty of the information.When the entropy is high, the information is so uncertain that we will need more information to make it clear.
We can find in Figure 1 that when α is less than 0.85, the two firms' innovation output Q I * is a certain value, so the entropy is low.With α increasing, the value of the innovation output Q I * goes from one certain value to two values.Finally when α is greater than 1.08, it has multiple values.With the increase of the uncertainty of the system, the entropy is increasing, so we can come to the conclusion that a too large decision parameter will lead to a large entropy to the system, and so that the companies will have to get more information to make an optimal decision.Figure 4 shows the maximal Lyapunov exponent with α varying from 0 to 1.2.The maximal Lyapunov exponent reflects the state of the system.When it is equal to zero, it corresponds to the critical point of bifurcation in the system, and when it is larger than zero, it means the system falls into chaos with increased entropy.We can find in Figure 4 that when α takes a small value, the equilibrium point is stable, but with α increasing to a certain degree, the system falls into chaos, so the conclusion is the same as above.β has a similar impact on the system as α has, as we have analyzed before.Figure 5 presents the parameter basin plots of the system with respect to parameters (α, β), we use different colors to describe different states of the system, and that is, the stable state(red), cycles of periods 2 (yellow), 4 (green), 8 (pink), chaos (white) and divergence (black).
Figure 5 presents the parameter basin plots of the system with respect to parameters ( , ), we use different colors to describe different states of the system, and that is, the stable state(red), cycles of periods 2 (yellow), 4 (green), 8 (pink), chaos (white) and divergence (black).Figure 6 shows the system's sensitivity to the initial value.We set = 1.2, = 1.0, in this condition the system falls into chaos.The firms' initial innovation inputs are respectively taken as = 2.83, = 2.75 and = 2.831, = 2.75, then we get variation of the innovation output decisions from period 0 to 100.We can find that the system shows a significant difference after about 10 cycles of iteration, so when the system is in the state of chaos with high entropy, even a minor change in the initial value will cause a huge difference in later decisions.Figure 6 shows the system's sensitivity to the initial value.We set α = 1.2, β = 1.0, in this condition the system falls into chaos.The firms' initial innovation inputs are respectively taken as 75, then we get variation of the innovation output decisions from period 0 to 100.We can find that the system shows a significant difference after about 10 cycles of iteration, so when the system is in the state of chaos with high entropy, even a minor change in the initial value will cause a huge difference in later decisions.
Entropy 2016, 18, 424 8 of 21 Figure 5 presents the parameter basin plots of the system with respect to parameters ( , ), we use different colors to describe different states of the system, and that is, the stable state(red), cycles of periods 2 (yellow), 4 (green), 8 (pink), chaos (white) and divergence (black).Figure 6 shows the system's sensitivity to the initial value.We set = 1.2, = 1.0, in this condition the system falls into chaos.The firms' initial innovation inputs are respectively taken as = 2.83, = 2.75 and = 2.831, = 2.75, then we get variation of the innovation output decisions from period 0 to 100.We can find that the system shows a significant difference after about 10 cycles of iteration, so when the system is in the state of chaos with high entropy, even a minor change in the initial value will cause a huge difference in later decisions.By analyzing the influence of the decision parameters on the stability of the system, we can conclude that the firms should not take an overlarge adjustment parameter when making innovation output decisions, otherwise the system will fall into an unstable state with high entropy and have an adverse impact on the firms' correctly making decisions and their profits.

The Influence of Government Innovation Subsidy Rate on the Stability Region
We further set a = 8, b = 2, c 1 = 0.5, c 2 = 0.6, θ = 0.6, λ = 1 and get the 3-D stability region of the system with change of α, β, ω in Figure 7 according to Jury stability criterion which is given in Equation (13).We can find that, as the government innovation subsidy rate ω increases, the stability region increases in the area of the α − β plane.So it can be concluded that an appropriate increase in the government subsidy rate of innovation can improve the stability of the system and decrease entropy.By drawing the stability region about α, β with different government subsidy rates (Figure 8), we can draw the conclusion more clearly.
Entropy 2016, 18, 424 9 of 21 By analyzing the influence of the decision parameters on the stability of the system, we can conclude that the firms should not take an overlarge adjustment parameter when making innovation output decisions, otherwise the system will fall into an unstable state with high entropy and have an adverse impact on the firms' correctly making decisions and their profits.

The Influence of Government Innovation Subsidy Rate on the Stability Region
We further set = 8, = 2, = 0.5, = 0.6, = 0.6, = 1 and get the 3-D stability region of the system with change of , , ω in Figure 7 according to Jury stability criterion which is given in Equation (13).We can find that, as the government innovation subsidy rate increases, the stability region increases in the area of the − plane.So it can be concluded that an appropriate increase in the government subsidy rate of innovation can improve the stability of the system and decrease entropy.By drawing the stability region about , with different government subsidy rates (Figure 8), we can draw the conclusion more clearly.By analyzing the influence of the decision parameters on the stability of the system, we can conclude that the firms should not take an overlarge adjustment parameter when making innovation output decisions, otherwise the system will fall into an unstable state with high entropy and have an adverse impact on the firms' correctly making decisions and their profits.

The Influence of Government Innovation Subsidy Rate on the Stability Region
We further set = 8, = 2, = 0.5, = 0.6, = 0.6, = 1 and get the 3-D stability region of the system with change of , , ω in Figure 7 according to Jury stability criterion which is given in Equation (13).We can find that, as the government innovation subsidy rate increases, the stability region increases in the area of the − plane.So it can be concluded that an appropriate increase in the government subsidy rate of innovation can improve the stability of the system and decrease entropy.By drawing the stability region about , with different government subsidy rates (Figure 8), we can draw the conclusion more clearly.

The Effect of the Government Innovation Subsidy Rate on Firms' Innovation Activities
We also set a = 8; b = 2; c 1 = 0.5; c 2 = 0.6; θ = 0.6; λ = 1, then we can get the equilibrium solution Entropy 2016, 18, 424 Figures 9 and 10 show how two firms' innovation outputs vary with the changes in government innovation subsidy.We can conclude that the innovation outputs Q I i increase along with the increase in the subsidy rate, and there's a trend of increasing marginal.It indicates that the policy of government subsidies based on the innovation inputs plays a positive role in encouraging firms' innovation activities, but we can see that when the subsidy rate is close to 1, the equilibrium solution is either too large or too small.In this condition the government pays almost all of the innovation activity funds, so firms may increase their innovation input infinitely, which will place too much fiscal burden to the government, so we should avoid the occurrence of this phenomenon in practice.Therefore, the government should synthetically determine a reasonable rate of subsidy to encourage firms' innovation activities.Figures 9 and 10 show how two firms' innovation outputs vary with the changes in government innovation subsidy.We can conclude that the innovation outputs increase along with the increase in the subsidy rate, and there's a trend of increasing marginal.It indicates that the policy of government subsidies based on the innovation inputs plays a positive role in encouraging firms' innovation activities, but we can see that when the subsidy rate is close to 1, the equilibrium solution is either too large or too small.In this condition the government pays almost all of the innovation activity funds, so firms may increase their innovation input infinitely, which will place too much fiscal burden to the government, so we should avoid the occurrence of this phenomenon in practice.Therefore, the government should synthetically determine a reasonable rate of subsidy to encourage firms' innovation activities.We have proved that the government's innovation subsidy policies based on the innovation input had a positive effect on firms' innovation activities, and the government should decide a reasonable subsidy rate .Then we will analyze the optimal subsidy decision that the government should make.
The government's goal is to maximaze the social welfare , and the social welfare is the sum of the consumer surplus and the producer surplus.According to the study of Lopez and Naylor [24], the consumer surplus is defined as: We also set = 8; = 2; = 0.5; = 0.6; = 0. Figures 9 and 10 show how two firms' innovation outputs vary with the changes in government innovation subsidy.We can conclude that the innovation outputs increase along with the increase in the subsidy rate, and there's a trend of increasing marginal.It indicates that the policy of government subsidies based on the innovation inputs plays a positive role in encouraging firms' innovation activities, but we can see that when the subsidy rate is close to 1, the equilibrium solution is either too large or too small.In this condition the government pays almost all of the innovation activity funds, so firms may increase their innovation input infinitely, which will place too much fiscal burden to the government, so we should avoid the occurrence of this phenomenon in practice.Therefore, the government should synthetically determine a reasonable rate of subsidy to encourage firms' innovation activities.We have proved that the government's innovation subsidy policies based on the innovation input had a positive effect on firms' innovation activities, and the government should decide a reasonable subsidy rate .Then we will analyze the optimal subsidy decision that the government should make.
The government's goal is to maximaze the social welfare , and the social welfare is the sum of the consumer surplus and the producer surplus.According to the study of Lopez and Naylor [24], the consumer surplus is defined as: We have proved that the government's innovation subsidy policies based on the innovation input had a positive effect on firms' innovation activities, and the government should decide a reasonable subsidy rate ω.Then we will analyze the optimal subsidy decision that the government should make.
The government's goal is to maximaze the social welfare W, and the social welfare is the sum of the consumer surplus CS and the producer surplus.According to the study of Lopez and Naylor [24], the consumer surplus CS is defined as: in which U is the consumer's utility function.
We assume that consumer's utility function U has a quadratic relationship with the production [25], which means: Entropy Then we can get that the social welfare W is: Putting Equations ( 15) and ( 16) into ( 17), then we can get: As we have analyzed above, we obtain p 1 Figures 11 and 12 show how two firms' innovation outputs vary with the innovation input parameter.We can find the innovation outputs Q I i decrease along with the increase of parameter λ.We know that λ measures the benefits brought by innovation acticities, and the smaller λ is, the more benefits firms will obtain from the same capital spent on innovation.Therefore, firms are more willing to engage in innovation when λ is smaller, so we advise that countries should support the cultivation of innovative talents and firms can improve their innovation ability by introducing talent, which will improve the earning of innovation and the innovation level of the whole society.
willing to engage in innovation when is smaller, so we advise that countries should support the cultivation of innovative talents and firms can improve their innovation ability by introducing talent, which will improve the earning of innovation and the innovation level of the whole society.

The Model
In this section, we assume that the innovation subsidies given by the government are not based on the innovation inputs but on the innovation outputs, which means that government financial subsidies will be greater when there is higher innovation outputs of firms.
The elements in this model are the same as the last section except that the government's innovation subsidy is ( ( )), in which m is the subsidy coefficient of innovation outputs, which represents the subsidies given for each increase of the market demand .So we can get the profit function of the two firms as:

Model Solving
Same as illustrated above, we get two firms' product prices in period t by maximizing their profits, then we get the equilibrium price as: willing to engage in innovation when is smaller, so we advise that countries should support the cultivation of innovative talents and firms can improve their innovation ability by introducing talent, which will improve the earning of innovation and the innovation level of the whole society.

The Model
In this section, we assume that the innovation subsidies given by the government are not based on the innovation inputs but on the innovation outputs, which means that government financial subsidies will be greater when there is higher innovation outputs of firms.
The elements in this model are the same as the last section except that the government's innovation subsidy is ( ( )), in which m is the subsidy coefficient of innovation outputs, which represents the subsidies given for each increase of the market demand .So we can get the profit function of the two firms as:

Model Solving
Same as illustrated above, we get two firms' product prices in period t by maximizing their profits, then we get the equilibrium price as:

The Model
In this section, we assume that the innovation subsidies given by the government are not based on the innovation inputs but on the innovation outputs, which means that government financial subsidies will be greater when there is higher innovation outputs of firms.
The elements in this model are the same as the last section except that the government's innovation subsidy is m(Q I i (t)), in which m is the subsidy coefficient of innovation outputs, which represents the subsidies given for each increase of the market demand Q I .So we can get the profit function of the two firms as:

Model Solving
Same as illustrated above, we get two firms' product prices in period t by maximizing their profits, then we get the equilibrium price as: Entropy 2016, 18, 424 13 of 20 Putting ( 25) into (24), we get: In the stage of innovation output decision, the firms make decisions based on the marginal profits, then:

Analysis of the Equilibrium Point and the Stability of the System
Let Q I i (t + 1) = Q I i (t), we have assumed that both firms carry out innovation activities, so exclude Q I i (t) = 0, then the unique equilibrium point of the system is: In which We analyze the equilibrium point and the stability of system by numerical simulation similar to that of Section 2.4, and we will show the parameter basin according to the Jury stability criterion in the next section.

The Influence of the Decision Parameters on the Stability and the Entropy of the System
Setting β = 0.6, we can describe how the dynamic behavior of the system varies with α in Figures 13-15.It is similar with Section 2.4, in that as α is increasing, bifurcation occurs in the system and then the system even falls into chaos.
We can also find that the entropy of the system is as α is increasing.We can find in Figure 13 that when α is less than 0.69, the two firms' innovation output Q I * has a certain value, so the entropy is low.Then with α increasing, the value of the innovation output Q I * goes from one certain value to two values.Finally when α is greater than 0.89, it has multiple values.
The increase of uncertainty of the innovation output Q I * will lead to the increase of the entropy of the system.Then we can come to the conclusion that a too large decision parameter will lead to a large entropy to the system, so companies will have to get more information to make an optimal decision.
Figure 16 shows the maximal Lyapunov exponent with α varying from 0 to 1.0.It can also reflect the transformation of the system from stability to bifurcation, then to chaos.The entropy of the system also increases in this process.
certain value to two values.Finally when is greater than 0.89, it has multiple values.The increase of uncertainty of the innovation output * will lead to the increase of the entropy of the system.Then we can come to the conclusion that a too large decision parameter will lead to a large entropy to the system, so companies will have to get more information to make an optimal decision.Figure 16 shows the maximal Lyapunov exponent with α varying from 0 to 1.0.It can also reflect the transformation of the system from stability to bifurcation, then to chaos.The entropy of the system also increases in this process.certain value to two values.Finally when is greater than 0.89, it has multiple values.The increase of uncertainty of the innovation output * will lead to the increase of the entropy of the system.Then we can come to the conclusion that a too large decision parameter will lead to a large entropy to the system, so companies will have to get more information to make an optimal decision.Figure 16 shows the maximal Lyapunov exponent with α varying from 0 to 1.0.It can also reflect the transformation of the system from stability to bifurcation, then to chaos.The entropy of the system also increases in this process.The increase of uncertainty of the innovation output * will lead to the increase of the entropy of the system.Then we can come to the conclusion that a too large decision parameter will lead to a large entropy to the system, so companies will have to get more information to make an optimal decision.Figure 16 shows the maximal Lyapunov exponent with α varying from 0 to 1.0.It can also reflect the transformation of the system from stability to bifurcation, then to chaos.The entropy of the system also increases in this process.Figure 17 shows the parameter basin plots of the system with respect to parameters ( , ), where we can find that the system is stable when the decision parameters ( , ) both take small values.We use different colors to describe different states of the system, and that is, the stable state (red), cycles of periods 2 (yellow), 4 (green), 8 (pink), chaos (white) and divergence (black).Figure 17 shows the parameter basin plots of the system with respect to parameters (α, β), where we can find that the system is stable when the decision parameters (α, β) both take small values.We use different colors to describe different states of the system, and that is, the stable state (red), cycles of periods 2 (yellow), 4 (green), 8 (pink), chaos (white) and divergence (black).Figure 17 shows the parameter basin plots of the system with respect to parameters ( , ), where we can find that the system is stable when the decision parameters ( , ) both take small values.We use different colors to describe different states of the system, and that is, the stable state (red), cycles of periods 2 (yellow), 4 (green), 8 (pink), chaos (white) and divergence (black).Figure 18 shows the system's sensitivity to the initial value.We set = 1.0, = 0.9, in this condition the system falls into chaos.The firms' initial innovation inputs are respectively taken as = 1.58, = 1.55 and = 1.581, = 1.55, then we get variation of the innovation output decisions from period 0 to 100.We can find that the system shows a significant difference after about 18 cycles of iteration, so when the system is in the state of chaos with high entropy, even a minor change in the initial value will make a huge difference to later decisions.Figure 18 shows the system's sensitivity to the initial value.We set α = 1.0, β = 0.9, in this condition the system falls into chaos.The firms' initial innovation inputs are respectively taken as 55, then we get variation of the innovation output decisions from period 0 to 100.We can find that the system shows a significant difference after about 18 cycles of iteration, so when the system is in the state of chaos with high entropy, even a minor change in the initial value will make a huge difference to later decisions.We further set = 8; = 2; = 0.5; = 0.6; = 0.6; = 1 and get the 3-D stability region of the system with change of ( , , ) in Figure 19 according to the Jury stability criterion.We can find that as the subsidy coefficient of innovation outputs increases, the stability region decreases

The Influence of the Subsidy Coefficient of Innovation Outputs to the Stability Region
We further set a = 8; b = 2; c 1 = 0.5; c 2 = 0.6; θ = 0.6; λ = 1 and get the 3-D stability region of the system with change of (α, β, m) in Figure 19 according to the Jury stability criterion.We can find that as the subsidy coefficient of innovation outputs m increases, the stability region decreases in the area of the α − β plane, so it can be concluded that the stability region can be affected by the value of m, but what is different from the last section is that the increase of the government's subsidy coefficient will weaken the stability of the system and increase the entropy.We further set = 8; = 2; = 0.5; = 0.6; = 0.6; = 1 and get the 3-D stability region of the system with change of ( , , ) in Figure 19 according to the Jury stability criterion.We can find that as the subsidy coefficient of innovation outputs increases, the stability region decreases in the area of the − plane, so it can be concluded that the stability region can be affected by the value of , but what is different from the last section is that the increase of the government's subsidy coefficient will weaken the stability of the system and increase the entropy.We can also draw the conclusion more clearly by drawing the stability region about α, β with different government subsidy coefficients (Figure 20).We can also draw the conclusion more clearly by drawing the stability region about , with different government subsidy coefficients (Figure 20).We can draw the conclusion from the expression of the equilibrium point and Figures 21 and 22 that firms' decisions on innovation outputs are proportional to the subsidy coefficient, which indicates that the policy of government subsidies based on the innovation outputs also plays a positive role in encouraging firms' innovation activities, so governments should determine an appropriate subsidy coefficient by comprehensively considering the cost and price of the product and the effect on social welfare.
We can draw the conclusion from the expression of the equilibrium point and Figures 21 and 22 that firms' decisions on innovation outputs are proportional to the subsidy coefficient, which indicates that the policy of government subsidies based on the innovation outputs also plays a positive role in encouraging firms' innovation activities, so governments should determine an appropriate subsidy coefficient m by comprehensively considering the cost and price of the product and the effect on social welfare.
We can draw the conclusion from the expression of the equilibrium point and Figures 21 and 22 that firms' decisions on innovation outputs are proportional to the subsidy coefficient, which indicates that the policy of government subsidies based on the innovation outputs also plays a positive role in encouraging firms' innovation activities, so governments should determine an appropriate subsidy coefficient by comprehensively considering the cost and price of the product and the effect on social welfare.We can draw the same conclusion from Figures 23 and 24 that the innovation output decreases with the increasing of parameter .
We have proved that the government's innovation subsidy policies based on the innovation output have a positive effect on firms' innovation activities, and the government should determine an appropriate subsidy coefficient .Just as the study in Section 2.4.3, we can get the optimal subsidy coefficient of innovation output .We can draw the same conclusion from Figures 23 and 24 that the innovation output decreases with the increasing of parameter .
We have proved that the government's innovation subsidy policies based on the innovation output have a positive effect on firms' innovation activities, and the government should determine an appropriate subsidy coefficient .Just as the study in Section 2.4.3, we can get the optimal subsidy coefficient of innovation output .We have proved that the government's innovation subsidy policies based on the innovation output have a positive effect on firms' innovation activities, and the government should determine an appropriate subsidy coefficient m.Just as the study in Section 2.4.3, we can get the optimal subsidy coefficient of innovation output m.
We also set a = 8; b = 2; c 1 = 0.5; c 2 = 0.6; θ = 0.6; λ = 1 as above, then we can get that: Let dW dm = 0, then we get that m = 6.71.When 0 ≤ m < 6.71, d 2 W dm 2 < 0. When m > 6.71, d 2 W dm 2 > 0, so when 0 ≤ m < 6.71, we will get the biggest social utility when m = 0, but in this condition, there is no promoting effect on enterprises' innovation activities.When m > 6.71, the increase of m will not only promote enterprises' innovation activities, but also improve the social welfare, so the government should set the subsidy coefficient m larger than 6.71.

Comparison of the Effects of the Two Kinds of Subsidy Policies on Innovation
In Sections 2 and 3, we have analyzed the stability and the equilibrium solution of the system under two kinds of government subsidies, one is based on the innovation inputs and the other is based on the innovation outputs.We can conclude that: (1) In the aspect of stimulating firms' innovation activities, the two kinds of subsidy policies both play a positive role, so the government can encourage firms' innovation activities by increasing the subsidy rate ω or increasing the subsidy coefficient of innovation outputs m.What is different is that under the first policy, the relationship of the innovation output decision and the subsidy rate has a marginal increasing character, while under the second policy, firms' innovation outputs decision is proportional to the subsidy coefficient.This means that the first policy, which is based on the innovation inputs, has a more significant incentive effect on innovation activities than the second policy, which is based on the innovation outputs.(2) In the aspect of stability of the system, the system will bifurcate and even fall into chaos with increasing entropy when taking large decision parameters (α, β) in both models.Besides, the subsidy rate ω and subsidy coefficient m can both affect the stability region of the system.What is different is that the increase of the subsidy rate ω can enhance the stability of the system while the increase of the subsidy coefficient m will weaken the stability of the system and increase entropy.

Conclusions
In this paper, we propose two different forms of government subsidies based on innovation input and innovation output, respectively, and we establish a dynamic game model between two bounded rational enterprises.In the process, we analyze the complexity of the model, and we discuss the effect of two different forms of government subsidies on stimulating the innovation activities.At the same time, we analyze how the innovation decision's adjusting speed, the degree of government subsidies and the beneficial coefficient of innovation influence the equilibrium, entropy and stability of the system.We draw the following conclusions: (1) Both kinds of innovation subsidy policies proposed in this paper have positive effects on firms' innovation activities.Therefore, the government can increase firms' enthusiasm for innovation by giving an appropriate government subsidy, which is beneficial to build a good innovative economic environment.(2) Assumed to be bounded rationally, the two firms make decisions of innovation outputs according to the marginal profit effect, but in the decision making process, the decision parameters α, β should not be too large, or the system will fall into the unstable state where decisions fluctuate disorderly and the entropy of the system will increase which means the companies will need more information to make an optimal decision.(3) The degree of the government innovation subsidies will impact the stability of the system.Under the subsidy policy based on the innovation inputs, the increase in subsidy rate ω can decrease the entropy and enhance the stability of the system, but under the subsidy policy based on the innovation outputs, the increase in the subsidy coefficient m will increase the entropy and weaken the stability of the system, so the government should synthetically consider the effect of the innovation subsidy on innovation incentives, the stability of the system, the budget of fiscal expenditure and the social benefits, and then decide a rational subsidy rate ω or subsidy coefficient m. (4) The innovation input parameter λ measures the benefits brought by innovation acticities, and the smaller is λ, the more benefits firms will obtain from the same investment in innovation, which means a higher level of inovation.The results derived from both models in this paper indicate that firms are more willing to engage in innovation when λ is small, so we advise that countries should support the cultivation of innovative talents and firms can improve their innovation ability by introducing talents, which will improve the earning of innovation and the innovation level of the whole society.
There are two limitations in this paper that should be noticed, and these point to possible future research directions.Firstly, we followed the classic linear demand function when establishing the game model, but because of the randomness of the market, there may be some distinction between the real demands and the ideal demands, so in future research we will improve the demand function, for example, by adding a probability distribution function in the demand function to adapt to stochastic demand in the real market.Secondly, we assume that the players in the game are both bounded rationally.Either of them is the rival for the other, so neither of the two players can learn about the exact innovation output decision that its rival makes.This assumption is a good fit for a tit-for-tat game, but there are some other types of game where there exists information sharing between the players, which may bring a win-win result.We leave this for future research.

Figure 1 .
Figure 1.The innovation output bifurcation with change of α.Figure 1.The innovation output bifurcation with change of α.

Figure 1 .
Figure 1.The innovation output bifurcation with change of α.Figure 1.The innovation output bifurcation with change of α.

Figure 1 .
Figure 1.The innovation output bifurcation with change of α.

Figure 2 .
Figure 2. The product price bifurcation with change of α.

Figure 3 .
Figure 3.The profit bifurcation with change of α.

Figure 4 .
Figure 4.The maximal Lyapunov exponent with change of α.

Figure 2 .
Figure 2. The product price bifurcation with change of α.

Figure 1 .
Figure 1.The innovation output bifurcation with change of α.

Figure 2 .
Figure 2. The product price bifurcation with change of α.

Figure 3 .
Figure 3.The profit bifurcation with change of α.

Figure 4 .
Figure 4.The maximal Lyapunov exponent with change of α.

Figure 3 .
Figure 3.The profit bifurcation with change of α.

Figure 1 .
Figure 1.The innovation output bifurcation with change of α.

Figure 2 .
Figure 2. The product price bifurcation with change of α.

Figure 3 .
Figure 3.The profit bifurcation with change of α.

Figure 4 .
Figure 4.The maximal Lyapunov exponent with change of α.Figure 4. The maximal Lyapunov exponent with change of α.

Figure 4 .
Figure 4.The maximal Lyapunov exponent with change of α.Figure 4. The maximal Lyapunov exponent with change of α.

Figure 5 .
Figure 5.The parameter basin of the system.

Figure 6 .
Figure 6.The sensitivity to initial values of the system.(a) Comparison of the evolutions of the system taking two different initial values; (b) ∆ in the period of 0-100 taking two different initial values.

Figure 5 .
Figure 5.The parameter basin of the system.

Figure 5 .
Figure 5.The parameter basin of the system.

Figure 6 .
Figure 6.The sensitivity to initial values of the system.(a) Comparison of the evolutions of the system taking two different initial values; (b) ∆ in the period of 0-100 taking two different initial values.

Figure 6 .
Figure 6.The sensitivity to initial values of the system.(a) Comparison of the evolutions of the system taking two different initial values; (b) ∆Q I 1 in the period of 0-100 taking two different initial values.

Figure 7 .
Figure 7.The 3-D stability region of the system with change of , , .

Figure 8 .
Figure 8.The stability regions of the system with changes in , (set different values of ).

Figure 7 .
Figure 7.The 3-D stability region of the system with change of α, β, ω.

Figure 7 .
Figure 7.The 3-D stability region of the system with change of , , .

Figure 8 .
Figure 8.The stability regions of the system with changes in , (set different values of ).

Figure 8 .
Figure 8.The stability regions of the system with changes in α, β (set different values of ω).

Figure 9 .
Figure 9.The relationship of and .

Figure 10 .
Figure 10.The relationship of and .

Figure 9 .
Figure 9.The relationship of Q I 1 and ω.

Figure 9 .
Figure 9.The relationship of and .

Figure 10 .
Figure 10.The relationship of and .

Figure 10 .
Figure 10.The relationship of Q I 2 and ω.

Figure 11 .
Figure 11.The relationship of and .

Figure 12 .
Figure 12.The relationship of and .

Figure 11 .
Figure 11.The relationship of Q I 1 and λ.

Figure 11 .
Figure 11.The relationship of and .

Figure 12 .
Figure 12.The relationship of and .

Figure 12 .
Figure 12.The relationship of Q I 2 and λ.

Figure 13 .
Figure 13.The innovation output bifurcation with change of α.

Figure 14 .
Figure 14.The product price bifurcation with change of α.

Figure 15 .
Figure 15.The profit bifurcation with change of α.

Figure 13 .
Figure 13.The innovation output bifurcation with change of α.

Figure 13 .
Figure 13.The innovation output bifurcation with change of α.

Figure 14 .
Figure 14.The product price bifurcation with change of α.

Figure 15 .
Figure 15.The profit bifurcation with change of α.

Figure 14 .
Figure 14.The product price bifurcation with change of α.

Figure 13 .
Figure 13.The innovation output bifurcation with change of α.

Figure 14 .
Figure 14.The product price bifurcation with change of α.

Figure 15 .
Figure 15.The profit bifurcation with change of α.

Figure 16 .
Figure 16.The maximal Lyapunov exponent with change of α.

Figure 16 .
Figure 16.The maximal Lyapunov exponent with change of α.

Figure 17 .
Figure 17.The parameter basin of the system.

Figure 17 .
Figure 17.The parameter basin of the system.

Figure 18 .
Figure 18.The sensitivity to initial values of the system.(a) Comparison of the evolutions of the system taking two different initial values; (b) ∆ in the period of 0-100 taking two different initial values 3.4.2.The Influence of the Subsidy Coefficient of Innovation Outputs to the Stability Region

Figure 18 .
Figure 18.The sensitivity to initial values of the system.(a) Comparison of the evolutions of the system taking two different initial values; (b) ∆Q I 1 in the period of 0-100 taking two different initial values.

Figure 18 .
Figure 18.The sensitivity to initial values of the system.(a) Comparison of the evolutions of the system taking two different initial values; (b) ∆ in the period of 0-100 taking two different initial values 3.4.2.The Influence of the Subsidy Coefficient of Innovation Outputs to the Stability Region

Figure 19 .
Figure 19.The 3-D stability region of the system with change of , , .

Figure 19 .
Figure 19.The 3-D stability region of the system with change of α, β, m.

Figure 20 .
Figure 20.The stability regions of the system with changes in , (set different values of ).

Figure 20 .
Figure 20.The stability regions of the system with changes in α, β (set different values of m).

Figures 23 .
Figures 23.The relationship of and .

Figure 22 . 21 Figure 22 .
Figure 22.The relationship of Q I 2 and m.
Figures 23.Therelationship of and .Figure 23.relationship of Q I 1 and λ.

Figures 23 .
Figures 23.The relationship of and .

Figure 24 .Figure 24 .
Figure 24.The relationship of and .Figure 24.The relationship of Q I 2 and λ.