Complexity
Volume 2020 (2020), Article ID 6645506, 12 pages
https://doi.org/10.1155/2020/6645506

Research Article

Chaotic Behaviors in a Nonlinear Game of Two-Level Green Supply Chain with Government Subsidies

Chang-Feng Zhu¹ and Qing-Rong Wang²

¹School of Traffic and Transportation, Lanzhou Jiaotong University, Lanzhou, Gansu 730070, China
²School of Electronic and Information Engineering, Lanzhou Jiaotong University, Lanzhou, Gansu 730070, China

Correspondence should be addressed to Qing-Rong Wang; wangqr@mail.lzjtu.cn

Received 12 October 2020; Revised 27 November 2020; Accepted 3 December 2020; Published 15 December 2020

Academic Editor: Y. T. Peter Bian

Copyright © 2020 Chang-Feng Zhu and Qing-Rong Wang. 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.

Abstract.

In this paper, a two-level green supply chain composed of a manufacturer and a retailer is taken as the background. Considering the consumer’s double consumption preference and the manufacturer’s green product R&D investment, a differential game model of the green supply chain under the government cost subsidy strategy is constructed. Firstly, the equilibrium points of the system are solved and their stability is discussed and analyzed. Secondly, the dynamic evolution process of Nash equilibrium under the parameters of green degree, green preference coefficient, retail channel preference coefficient, coefficient of the sensitivity of price, and adjustment speed are described by numerical simulation. The results show that the two ways of a system entering chaos are Flip bifurcation and N-S bifurcation, respectively, by 2D bifurcation graph, and it is also verified in 1D bifurcation diagram. When the bifurcation parameters are small, the system maintains Nash equilibrium stability. If the green degree of products is increased, the green preference coefficient will also increase; on the contrary, the retail preference coefficient will decrease. Research and development cost subsidy policy can effectively improve the green degree of products and increase the sales volume of products, so as to improve the profit of supply chain members.

1. Introduction

In recent years, more and more attention has been paid to environmental protection in the world, and corresponding environmental protection laws and regulations have been formulated. On the other hand, consumers’ awareness of environmental protection is also strengthening. These behaviors urge suppliers to implement green production, so as to reduce the environmental damage caused by people’s production activities [1–3]. The concepts of green economy and ecological protection have been gradually and deeply rooted in the hearts of the people, which have changed the production and business model of enterprises. Green production capacity has become an important indicator for enterprises to survive in the industry for a long time. Manufacturing enterprises improve the quality of products by using green and environment-friendly materials. Green products can not only attract the attention of consumers, so as to expand product demand and improve product competitiveness, but also reduce environmental pollution and create a good social image [4–6]. The State Council has also encouraged the green production of enterprises, guided the green consumption of the masses, and proposed to build a green manufacturing system. However, when manufacturing enterprises are faced with many obstacles in the process of green production, such as insufficient funds, technological innovation, and insufficient support for product research and development, the enthusiasm of green product production is affected. Therefore, the government has formulated the corresponding green subsidy policy in order to subsidize the enterprises producing green products and promote the green development of the supply chain.

Many scholars have done in-depth research on the green supply chain from many angles and achieved satisfactory research results. Ghosh and shah [7] studied the secondary supply chain composed of a manufacturer and a retailer. The demand of the market is determined by the price of the product and the green degree of the product. The relationship between the price of the product and the green degree of the product depends on the pricing between them. Li et al. [8] built a green supply chain game model under price strategy based on Stackelberg game. The paper mainly analyzed the advantages and disadvantages of competition between decentralized decision and centralized decision-making. The research results showed that when the green cost exceeds a certain critical value, manufacturers will open direct sales channels, and the price of green supply chain products under the centralized strategy is higher than that under decentralized decision-making. Literature [9–11] discussed and analyzed the impact of consumers’ green preference and environmental protection awareness on green supply chain decision-making. Huang et al. [2, 12] considered the pricing decision of closed-loop supply chain with the green degree of products under the same preference, and they found that green degree and fair preference will change the retail price, wholesale price, and waste recovery rate of products and also affect the total profit of the supply chain. Many scholars also studied the green supply chain system under government intervention [4]. For example, Sheng and other scholars [13–20] studied the impact of government subsidies on enterprise innovation. Madani et al. [21–24] considered the game model under the situation of green consumption subsidies and nongreen taxation and analyzed and discussed the influence of government fiscal strategy on the optimal strategy of supply chain members under centralized and decentralized decision-making. In this paper, we consider increasing innovation subsidies and construct effective measures that can make innovation subsidies become incentives for enterprises to produce green products.

In conclusion, the existing literature on green supply chain mainly focused on the green degree of products, government subsidies, and preferential policies of supply chain members. Little literature established a green supply chain model based on consumers’ double consumption preference, and few scholars studied the impact of product greenness and government subsidy policies on enterprise profits from the perspective of dynamics. Therefore, this paper introduces the green preference and channel preference into the green supply chain management under the government subsidy strategy and constructs the green supply chain differential game model under the Government R&D cost subsidy. This paper discusses the impact of double consumption and government subsidy policy on product green degree, product price, and supply chain member profit.

2. Establishment of Model

Consider a green supply chain system consisting of a manufacturer and a retailer. The manufacturer sells through two ways: one is to use traditional channels for wholesale sale and the other is to sell directly online. After comprehensive consideration of market factors, the manufacturer decides that the direct selling price and wholesale price of the product are and , respectively, and the retailer determines the retail price according to the manufacturer’s decision information and market demand. In addition, in order to encourage manufacturers to produce green products, the government decided to take the R&D cost subsidy strategy to intervene with manufacturers. In order to facilitate the solution and the establishment of the model, we make the following assumptions about the model:

 ① Manufacturers not only need to pay fixed production costs to produce green products but also need to purchase new equipment or technological innovation to increase additional investment in green products. We assume that the relationship between R&D costs of green products produced by manufacturers and greenness is , where is the R&D cost coefficient of green products. In addition, it is assumed that the unit cost of the product produced by the manufacturer is . ② Suppose that consumers have a preference for both channel and green. By the description of the demand function in [7, 25], we assumed that the demand of the retail channel and direct channel is a linear function and the cross price elasticity coefficient is symmetric. Thus, the demand functions of retail channel and direct channel are, respectively, expressed as

Among them, represents the potential market capacity, represents the retail channel preference coefficient, represents the green preference coefficient, is the price sensitivity coefficient, is the cross price sensitivity coefficient, and there are . In order to better analyze the characteristics of the green supply chain, referring to [26], is defined as the efficiency coefficient of product greening. The higher the value of is, the more the public likes green products, or it can be interpreted as that products with the same green quality need to pay less research and development costs.

In order to encourage manufacturers to improve the green quality of their products, the government considers subsidy strategies for green products [7, 27]. The government subsidizes the R&D cost of the manufacturer’s products. Assuming that the R&D cost subsidy coefficient is , the government subsidy expenditure under the R&D cost subsidy strategy is .

In order to ensure the theoretical and practical significance of the supply chain, the model needs to meet the following constraints:

Therefore, when the manufacturing enterprise considers the public’s green preference and channel preference and the government subsidizes the manufacturer according to the product R&D cost, the profit functions of the retailer and the manufacturer are, respectively,

By substituting equations (1) and (2) into (4) and (5), the profit functions of retailers and manufacturers under the dual preference and government subsidy strategy are, respectively,

By taking partial derivatives of equations (6) and (7) with respect to and , respectively, we can obtain the optimal first-order condition for the profits of retailers and manufacturers under dual preferences as follows:

If both the retailer and the manufacturer are bounded rational, that is, the two sides of the game do not fully understand each other’s relevant information. Then, they determine the price of products in period according to the local estimation of the marginal target . In other words, in period, if , then the manufacturer and retailer will increase the product price in period; on the contrary, if , then the manufacturer and retailer will decrease the product price in period. According to this idea, the dynamic adjustment mechanism [28, 29] is introduced as follows:

ｗhere is the price adjustment speed of the retailer and the manufacturer, respectively. Combining equations (8), (9), and (10), a two-dimensional nonlinear discrete difference equation is obtained as

For the convenience of calculation, let and ; then, equation (11) can be written in the following form:

3. Stability Analysis of the Model

Let ; then, we could directly get four equilibrium solutions of mapping (12) as follows: , , , and .

Notice that is the trial fixed point, and are boundary fixed points, and is the internal fixed point. In order to guarantee that all the equilibrium points in system (12) have practical economic significance, the four equilibrium points are nonnegative. From the definition of parameters, we know , , , and . Therefore, we can get , , , and . Thus, it is concluded that all four equilibrium points are nonnegative, and all equilibrium points are of economic significance. For the sake of obtaining the global stable conditions of fixed points of mapping (12), then, the Jacobian matrix of mapping (12) at any point is given as follows:

Proposition 1. is an unstable source.

Proof. The Jacobian matrix of is

Obviously, is a diagonal matrix and eigenvalues of are elements located on diagonal; that is, and . Because of , , , so ; thus, is an unstable node.

Proposition 2. When , the equilibrium point is a saddle point, and when , the equilibrium point is an unstable node.

Proof. By substituting equilibrium point into equation (13), the Jacobian matrix of can be expressed as

The eigenvalues of matrix (15) are and . From the definition of the previous parameters and proof, we can know that , , and . So, there is . When , , equilibrium point is a saddle point; when , , then is an unstable node.

Proposition 3. When , the equilibrium point is a saddle point; when , the equilibrium point is an unstable node.
In order to avoid redundancy, the proof process is no longer described.

Proposition 4. If and only if condition (19) holds, the Nash equilibrium is stable.

Proof. The Jacobian matrix of is

The specific forms of trace and determinant of are given as

According to the Jury condition, the local stability conditions of the Nash equilibrium point are as follows:

And by the definition of the parameters, we know , , and , so . Obviously, the second condition in the inequality system (18) holds. Then, inequality group (18) can be rewritten as

Through the above analysis, the Nash equilibrium is stable if and only if the inequality set (19) is satisfied. Fix the initial value as and select parameters , , , , , , , and . Under this set of parameters, the stable region of equilibrium point with respect to price adjustment speed is shown in the blue area in Figure 1. In this stable region, the prices and of manufacturer and retailer converge to the Nash equilibrium after many games. If the value of the adjustment speed or is increased, the point will jump out of the stable region, and complex dynamic phenomena will occur in the later process.

Figure 1: Local stable region of the equilibrium point .

4. Numerical Bifurcation Analysis

As we all know, people mainly study the complex dynamic behavior of the nonlinear dynamic system by numerical simulation, and the tools used in the process of numerical simulation mainly include 1D bifurcation diagram, 2D bifurcation diagram, Lyapunov exponent diagram [30, 31], and phase diagram. In this section, the influence of the parameter changing on the stability of the Nash equilibrium of system (12) is discussed by numerical analysis.

4.1. The Influence of Price Sensitive Coefficient on the System

In order to study the effect of price sensitive coefficient on the system, we select a group of parameters such as , , , , , , and . When the value of the parameter increases gradually, 2D bifurcation diagrams of system (12) regarding the adjusted speeds and are obtained, as shown in Figure 2. The different colors in the 2D bifurcation diagram represent different periods (i.e., the color displayed by the color bar on the right side of the bifurcation diagram). The brown region represents the Nash equilibrium stable region, namely, 1-period, light green represents 2-period, dark green represents 3-period, orange represents 4-period, and so on, black represents quasiperiodic state or chaos, and white area represents escape state. These three states cannot be calculated at present. The stable region of the Nash equilibrium point with respect to adjustment speeds and is shown in the brown area in Figure 2. In the parameter plane , the manufacturer’s direct selling price and retailer’s selling price will converge to the Nash equilibrium point after many games. If we increase the value of or , it will cause or to jump out of the stable region. From Figures 2(a), 2(b), 2(c), and 2(d), we observe that the system has two different paths into chaos. One is that, with the increase of parameter values and , system (12) successively passes through the brown region (1-cycle), green region (2-cycle), orange region (4-cycle), light green region (8-cycle), and finally enters the black region; in other words, the system enters into the chaotic state by Flip bifurcation. When manufacturers and retailers determine their respective adjustment speed according to the above path, the system will enter into chaos through a period-doubling, which also means that a supply chain system composed of manufacturers and retailers will end up in chaos. With the increasing speed of adjustment, manufacturer or retailer may not survive and the supply chain will collapse. The other is that when the parameters enter the black region from the brown region (1-cycle) through the green region (2 cycles), the system enters into quasiperiod and finally enters chaos through the Neimark–Sacker bifurcation. From these four pictures, we find that the path of the system into chaos is the same. The different value of price sensitivity coefficient makes the unique graph of the score fork, and we also feel the complex dynamic phenomenon. Combined with the 2D bifurcation diagram and the above analysis, it is found that the price sensitive coefficient will not affect the bifurcation way of the system, but will change the dynamic phenomenon of the fork. At the same time, we find that, in order to make the manufacturer and retailer obtain larger profits, the value of adjustment speed should not be too large. Once the threshold value is exceeded, the system may collapse. Therefore, choosing a smaller adjustment speed can promote the stable development of the green supply chain system composed of manufacturer and retailer.

(a)

(b)

(c)

(d)

Figure 2: 2D bifurcation diagram of the system on the parameter plane , where (a) , (b) , (c) , and (d) .

4.2. The Influence of Green Degree on Retail Channel Preference Coefficient and Green Preference Coefficient

In order to better analyze the influence of green degree on retail channel preference coefficient and green preference coefficient on the system, we choose a set of parameter values as , , , , , , , and . Taking as the bifurcation variable, the 1D bifurcation diagrams and their corresponding maximum Lyapunov exponent are obtained as shown in Figures 3(a), 3(b), 3(c), and 3(d). When , the 1D bifurcation diagram as shown in Figure 3(a) is obtained. When , the system is in Nash equilibrium and steady state, which indicates that the green supply chain system develops steadily. When , the system begins bifurcation. When , the system is in a periodic state. When , the system is in a chaotic state, which indicates that the green supply chain system composed of a manufacturer and a retailer is in chaotic competition. When , the system experiences three-time bifurcations from period to chaos, which indicates that the manufacturer and retailer experience three price games. Similarly, when , 1D bifurcation diagram is shown in Figure 3(c). When , the system is in Nash equilibrium and stable state. When , the system bifurcates. When , the system is in a period state. When , the system alternates between chaos and periodic bifurcation. Figures 3(b) and 3(d) are the largest Lyapunov exponent graphs corresponding to Figures 3(a) and 3(c). Lyapunov exponent diagram is an important tool to study the nonlinear dynamic behavior. In the nonlinear dynamical system, the Lyapunov exponent graph describes the convergence or divergence of adjacent trajectories. It distinguishes the regular attractor from singular attractor. When the Lyapunov exponent is less than zero, we can judge that the attractor is a regular attractor. If at least one of the Lyapunov exponents is positive, then we can infer that there are singular attractors or chaotic attractors in the system. When and (Figures 3(e) and 3(g) are obtained, resp.), it can be seen that the inverse N-S bifurcation occurs first in the system, followed by positive N-S bifurcation. It is obvious from these two diagrams that the stability interval of Figure 3(g) is smaller than that of Figures 3(e), 3(f), and 3(h) which are the maximum Lyapunov exponent diagrams corresponding to Figures 3(e) and 3(g), respectively. The system experiences three-time bifurcation of entering chaos from the period, which indicates that the manufacturer and retailer experience three times of price game. When , it means that the system is in a stable state; when , it indicates that the system bifurcates at this point; when , the system enters into a chaotic state. It can be seen from the diagrams that when the green preference coefficient is small and the value of the retail channel preference coefficient is in the middle value, the system will be in a stable state. At the same time, we also find that if the green degree of the product is increased, the green preference coefficient will also increase. On the contrary, when the green degree of the product increases, the value of the retail preference coefficient will decrease.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 3: 1D bifurcation diagram with respect to green preference coefficient and retail preference channel coefficient and their corresponding maximum Lyapunov exponent diagrams. Select a set of parameter values as , , , , , , and , where (a) and , (c) and . Figures 3(b) and 3(d) are the maximum Lyapunov exponent diagrams corresponding to Figures 3(a) and 3(c), where (e) and , (g) and . Figures 3(f) and 3(h) are the maximum Lyapunov exponent graphs corresponding to Figures 3(e) and 3(g).

4.3. The Influence of Other Parameters on the System

In order to further analyze the influence of other parameters on system (12), we select , , , , , and , respectively, as the bifurcation variables and obtain the 1D bifurcation diagrams and the maximum Lyapunov exponent diagrams, as shown in Figures 4. In Figures 4(a) and 4(b), we can see that the internal equilibrium point is stable when the value of parameters and is relatively small. However, as the value of parameter continues to increase, the internal equilibrium is no longer stable. Suddenly, irregular dynamic behaviors appear, including N-S bifurcation, periodic window, and chaos. In Figure 4(b), the system experiences 2-period, 4-period, 2-period, and finally enters chaos; that is, the system enters chaos through Flip bifurcation. It can be seen from Figures 4(c), 4(d), and 4(f) that the stable 2-periodic loop loses its stability through the N-S bifurcation, and the system finally enters chaos. In Figure 4(e), the bifurcation diagram with respect to the green degree is depicted. We can see that the period-doubling bifurcation and N-S bifurcation appear simultaneously in the process of bifurcation. In addition, we also find that periodic windows appear in the process of bifurcation in Figures 3(c), 3(d), 3(e), and 3(f). In fact, the periodic window corresponds to the part of the Lyapunov exponent closed to zero. When the Lyapunov exponent is equal to zero, the system bifurcates at this point; when the Lyapunov exponent is greater than zero, the system enters into a chaotic state. It is found that when the bifurcation parameters are small, the system is in a steady state. In this case, the supply chain system composed of a single manufacturer and a single retailer develops stably. At the same time, it is found that increasing the green degree of products can make the system more stable. However, due to the limitation of research and development costs, manufacturers will not increase the green degree of products to meet the green preference of consumers. In other words, increasing the green degree of products can improve the sales of products and the profits of supply chain members.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: 1D bifurcation diagrams and maximum Lyapunov exponent diagrams of system (12) with fixed parameters are , , , and . (a) Taking as a bifurcation parameter, the other parameters are , , , , and . (b) Taking as a bifurcation parameter, the other parameters are , , , , and . (c) Taking as bifurcation parameter, the other parameters are , , , and . (d) Taking as a bifurcation parameter, the other parameters are , , , , and . (e) Taking as a bifurcation parameter, the other parameters are , , ,, and . (f) Taking as a bifurcation parameter, the other parameters are , , , , and .

5. Global Dynamics Analysis of the System

People often pay attention to the final state of the system in nonlinear dynamics. The final state reflects the final state of the system. The final state can be expressed by the attractor. In the process of studying the nonlinear dynamic system, it is found that the change of parameters will cause variation in number and structure of attractors. Then, the influence of cross price sensitivity coefficient on the final behavior of system (12) is discussed through the evolution process of the attractor.

Parameters , , , , , , , , and are fixed and is taken as the variable to obtain the attractor evolution diagrams, as shown in Figure 5. When , a 2-periodic attractor is observed in Figure 5(a). With the increase in the cross price sensitivity coefficient , the 2-period attractor evolves into periodic cycles with rough edges as shown in Figure 5(b). When , the original 2-period ring with rough edges recedes and gradually evolved into a larger 2-period ring. When the value of parameter continues to increase, the attractor becomes larger and larger. Until , the 2-period cycles break to form many very small attractors, which are called phase-locking, as shown in Figure 5(e). The broken attractors recombine to form two attractors which look like a “ring” as shown in Figure 5(f). It is found that the number and shape of attractors change significantly when the cross price sensitivity coefficient increases. When , two strange attractors joint together, as shown in Figure 5(g). With the further increase of the parameter , two interlocking chaotic attractors evolve into a chaotic attractor as shown in Figure 5(h). Finally, the system is in a chaotic state, which means that the supply chain composed of a manufacturer and a retailer is in a chaotic competition scenario. Therefore, it shows that increasing the value of cross price sensitivity coefficient will make the market more and more chaotic, which is against the win-win intention of manufacturer and retailer. Therefore, choosing a smaller cross price sensitivity coefficient is conducive to the common development of manufacturers and retailers.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5: Diagrams of attractors evolution. (a) , (b) , (c) , (d) , (e) , (f) , (g) , and (h) .

6. Conclusion

In this paper, a two-level green supply chain composed of a single manufacturer and a single retailer is taken as the background. Firstly, considering the consumer’s double consumption preference and the manufacturer’s green product R&D investment, the game model of the green supply chain under the government cost subsidy strategy is constructed. Secondly, the dynamic evolution process of Nash equilibrium under the parameters of green degree, green preference coefficient, retail channel preference coefficient, coefficient of sensitivity of price, and adjustment speed is described by numerical simulation. The results show that the two ways of the system entering into chaos are Flip bifurcation and N-S bifurcation. The N-S bifurcation is mainly discussed in this paper. When the bifurcation parameters are small, the system maintains Nash equilibrium stability, and the stability of the system is destroyed if it exceeds a certain threshold value. Increasing the green degree of the product will increase the green preference coefficient, but the retail preference coefficient will decrease. The subsidy policy of R&D cost plays an irreplaceable role in encouraging the manufacturer to produce green products. Green sales can effectively improve the sales volume of products and realize the win-win benefits between manufacturers and retailers.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 71961016) and the Foundation of A Hundred Youth Talents Training Program of Lanzhou Jiaotong University.

References

Q. Zhu, J. Sarkis, and K.-H. Lai, “Institutional-based antecedents and performance outcomes of internal and external green supply chain management practices,” Journal of Purchasing and Supply Management, vol. 19, no. 2, pp. 106–117, 2013.

Z.-Y. Zhang, D.-X. Fu, and Q. Zhou, “Optimal decisions of a green supply chain under the joint action of fairness preference and subsidy to the manufacturer,” Discrete Dynamics in Nature and Society, vol. 2020, Article ID 9610503, 18 pages, 2020.

J. Li, Y. Yi, and H. Yang, “Information investment and sharing in a two-echelon supply chain under government subsidy and consumer preference for energy-saving products,” Mathematical Problems in Engineering, vol. 2020, Article ID 3935762, 18 pages, 2020.

Y. Tian, K. Govindan, and Q. Zhu, “A system dynamics model based on evolutionary game theory for green supply chain management diffusion among Chinese manufacturers,” Journal of Cleaner Production, vol. 80, pp. 96–105, 2014.

Q. H. Zhu and Y. J. Dou, “Evolutionary model between governments and core-enterprises in green supply chains,” Systems Engineering-Theory & Practice, vol. 27, no. 12, pp. 85–89, 2007.

X. Liu, K. Lin, L. Wang, and L. Ding, “Pricing decisions for a sustainable supply chain in the presence of potential strategic customers,” Sustainability, vol. 12, no. 4, p. 1655, 2020.

D. Ghosh and J. Shah, “A comparative analysis of greening policies across supply chain structures,” International Journal of Production Economics, vol. 135, no. 2, pp. 568–583, 2012.

B. Li, M. Zhu, Y. Jiang, and Z. Li, “Pricing policies of a competitive dual-channel green supply chain,” Journal of Cleaner Production, vol. 112, no. Part 3, pp. 2029–2042, 2016.

K. Chitra, “In search of the green consumers: a perceptual study,” Journal of Services Research, vol. 7, no. 1, pp. 173–291, 2007.

L. Zhang, J. Wang, and J. You, “Consumer environmental awareness and channel coordination with two substitutable products,” European Journal of Operational Research, vol. 241, no. 1, pp. 63–73, 2015.

D. Ghosh, J. Shah, and S. Swami, “Product greening and pricing strategies of firms under green sensitive consumer demand and environmental regulations,” Annals of Operations Research, vol. 290, no. 1-2, pp. 491–520, 2020.

J. Wen and Y. Li, “Profit distribution model of green building supply chain with fairness preferences and cap-and-trade policy,” IOP Conference Series: Earth and Environmental Science, vol. 237, Article ID 052043, 2019.

J.-Y. Chen, S. Dimitrov, and H. Pun, “The impact of government subsidy on supply chains’ sustainability innovation,” Omega, vol. 86, pp. 42–58, 2019.

Z. Huang, G. Liao, and Z. Li, “Loaning scale and government subsidy for promoting green innovation,” Technological Forecasting and Social Change, vol. 144, pp. 148–156, 2019.

M. Franco, H. Haase, and A. Pereira, “Empirical study about the role of social networks in SME performance,” Journal of Systems and Information Technology, vol. 18, no. 4, pp. 383–403, 2016.

F.-M. Hsu, D.-J. Horng, and C.-C. Hsueh, “The effect of government-sponsored R & D programmes on additionality in recipient firms in Taiwan,” Technovation, vol. 29, no. 3, pp. 204–217, 2009.

W. Wang, Y. Zhang, Y. Li, X. Zhao, and M. Cheng, “Closed-loop supply chains under reward-penalty mechanism: retailer collection and asymmetric information,” Journal of Cleaner Production, vol. 142, no. Part 4, pp. 3938–3955, 2017.

W. Wang, Y. Zhang, K. Zhang, T. Bai, and J. Shang, “Reward-penalty mechanism for closed-loop supply chains under responsibility-sharing and different power structures,” International Journal of Production Economics, vol. 170, pp. 178–190, 2015.

S. Zhang, Y. Yu, Q. Zhu, C. M. Qiu, and A. Tian, “Green innovation mode under carbon tax and innovation subsidy: an evolutionary game analysis for portfolio policies,” Sustainability, vol. 12, no. 4, p. 1385, 2020.

X. Deng, “The incentive effect of modes of government subsidies on technological innovation of LED lighting enterprises in China,” Light & Engineering, vol. 24, no. 3, pp. 145–148, 2016.

S. R. Madani and M. Rasti-Barzoki, “Sustainable supply chain management with pricing, greening and governmental tariffs determining strategies: a game-theoretic approach,” Computers & Industrial Engineering, vol. 105, pp. 287–298, 2017.

D. Yang and T. Xiao, “Pricing and green level decisions of a green supply chain with governmental interventions under fuzzy uncertainties,” Journal of Cleaner Production, vol. 149, pp. 1174–1187, 2017.

S. Zhan, Z. Shu, and H. Jiang, “Research on two-echelon green supply chain decision under government subsidy,” American Journal of Industrial and Business Management, vol. 8, no. 3, pp. 487–495, 2018.

A. Xu and Z. Zhou, “A pricing model for governments’ subsidy in the green supply chain,” International Journal of Networking and Virtual Organisations, vol. 14, no. 1/2, p. 40, 2014.

Z. Liu, T. D. Anderson, and J. M. Cruz, “Consumer environmental awareness and competition in two-stage supply chains,” European Journal of Operational Research, vol. 218, no. 3, pp. 602–613, 2012.

W. Zhu and Y. He, “Green product design in supply chains under competition,” European Journal of Operational Research, vol. 258, no. 1, pp. 165–180, 2017.

S. Swami and J. Shah, “Channel coordination in green supply chain management,” Journal of the Operational Research Society, vol. 64, no. 3, pp. 336–351, 2013.

W. Zhou, J. Zhou, T. Chu, and H. Li, “A dynamic duopoly cournot model with R & D efforts and its dynamic behavior analysis,” Complexity, vol. 2020, Article ID 4580415, 19 pages, 2020.

L. Zhou, J. Zhou, T. Chu, and H. Li, “Dynamic positioning in ice,” Encyclopedia of Ocean Engineering, vol. 2020, Article ID 9634878, 8 pages, 2020.

L. Gori, M. Sodini, and L. Fanti, “A nonlinear cournot duopoly with advertising,” Chaos, Solitons & Fractals, vol. 79, pp. 178–190, 2015.

A. Agliari, N. Pecora, and A. Szuz, “Dynamical analysis of cournot oligopoly models: Neimark-Sacker Bifurcation and related mechanisms,” in Qualitative Theory of Dynamical Systems, Tools and Applications for Economic Modelling, Springer Proceedings in Complexity, Urbino, Italy, 2016.