Next Article in Journal
Prioritizing Critical Management Strategies to Improving Construction Productivity: Empirical Research in Singapore
Next Article in Special Issue
Formulating Service-Oriented Strategies for Servitization of Manufacturing Companies
Previous Article in Journal
Emotion Regulation and Employability: The Mediational Role of Ambition and a Protean Career among Unemployed People
Previous Article in Special Issue
Development of a New Green Indicator and Its Implementation in a Cyber–Physical System for a Green Supply Chain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Dynamic Pricing and Service Decisions in the Reverse Supply Chain considering Consumers’ Service Sensitivity

College of Economics and Management, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(22), 9348; https://doi.org/10.3390/su12229348
Submission received: 12 October 2020 / Revised: 1 November 2020 / Accepted: 7 November 2020 / Published: 11 November 2020

Abstract

:
In this work, considering consumers’ service sensitivity for the third-party collector in the reverse supply chain, a game model (including a manufacturer and a third-party collector) with dual-channel recycling structure was established. Based on the bounded rationality, we used dynamics theory to analyze the dynamic service and pricing decisions of the manufacturer and the third-party collector in the process of collecting waste products. The results showed that adjustment speeds of pricing and service decisions have great effect on the stability of recycling market, and the recycling market is most sensitive to the change of third-party collectors’ service decision. When the recycling market is in chaos, it has a strong initial value sensitivity. In addition, we found that increase of consumers’ service sensitivity and service cost parameter had negative effect on stability of the recycling market; there is a reverse relationship between price adjustment parameters and the service adjustment parameter. With the increase of service adjustment parameter, the bifurcation point of the system will appear more quickly. Finally, chaotic market can be controlled by the feedback control method.

1. Introduction

With the rapid development of technology, the speed of upgrading products becomes faster. A large number of products are eliminated every year. These waste products are usually abandoned casually or disposed of in inappropriate ways, which causes great damage to the environment and people’s health. Along with the shortage of natural resources, enterprises need to change their traditional production mode and consider a sustainable development mode. Thus, the closed-loop supply chain (CLSC) appears; in contrast to traditional supply chains, CLSC contains the process of collecting waste products from consumers and remanufacturing [1], which is regarded as the reverse supply chain [2]. Collecting and remanufacturing waste products can effectively reduce environmental pollution and bring economic benefits to enterprises [3]. In recent years, remanufacturing has been successful in many areas such as computers, cameras, and car batteries. For example, Apple Inc. makes great profits through the business of recycling used mobile phones, and it is constantly expanding its recycling business [4,5].
The collection of waste products is an important part of the reverse supply chain and ensures the continuation of the remanufacturing activities, and it can be done by different members in the supply chain. Savaskan [6] introduced three different recycling models (manufacturer recycling, retailer recycling, and the third-party recycling) and studied the forward supply chain pricing decision under different recycling channels. These three ways of recycling are common and effective in the real world. However, with the rapid development of production, these single recycling modes are gradually unable to meet the demand for recycling. Then, dual-channel recycling model appeared. Compared with the single-channel recycling structure, the dual-channel recycling structures tend to be more effective [7]. In this paper, we considered the dual-channel recycling structure. Specifically, we focus on the mode of joint recycling by the manufacturer and the third-party collector.
In recent years, many independent third-party collectors have appeared in the recycling market and make great contributions to the collection of waste products. They make profit by collecting or handling waste products. Furthermore, different from the recycling way of manufacturers and retailers, third-party collectors always provide different service for consumers. For example, Aihuishou.com not only bears the postage cost of consumers but also provides pick-up service [8]. Service is an important factor in supply chain, which can improve consumer satisfaction and increase profit [9], and it is also helpful to promote the demand of products [10,11,12]. The study of supply chain considering the factor of service has been an issue of concern. Ren et al. [13] analyzed the return strategy of customers when service level competition was considered in the model. Dan [14] found that the profit of the manufacturer and the sales service integrator are positively correlated with the consumer’s sensitivity coefficient to the service level. Therefore, besides the pricing decision, we also consider the service decision of the third-party collector in the reverse supply chain.
When studying pricing and service decisions in the supply chain, much of the literature assumed that participants are completely rational and pursue individual interests. Hence, the optimal decision is often obtained during the single period. Without a doubt, the supply chain is a complex economic system, and it is difficult for decision makers to know all the market information. Especially in the recycling market, the decision maker can hardly know the complete market information such as potential total waste products quantity, consumers’ different preferences, and strategies of competitors. Hence, the participants tend to be bounded rationality. Along with incomplete information, the behavior of participants in the recycling market tends to be complex and the result may not converge to the Nash equilibrium, and complex phenomena like bifurcation and chaos may occur [15]. In this case, traditional static optimization methods are not suitable and dynamic game model should be established to analyze the supply chain.
By the above discussion, we found that the recycling mode, service, and bounded rationality have great effect on the optimal operational decision in the supply chain. The main purpose of this paper is to analyze the dynamic change of pricing and service decisions under bounded rationality in the reverse supply chain and discuss how consumers’ service sensitivity affect the decision. Then, a dynamic model which described the recycling market was established. By analyzing the local stability of equilibrium point and simulation, this study discussed the dynamic evolution of pricing and service decisions and the influence of different adjusting parameters and parameters related to service on the stability of the recycling market.
In this study, we will address the following questions.
(1) In the dual-channel recycling supply chain, how do adjustment parameters of recycling price and service affect the stability of the recycling market in the dynamic evolution? Does the market have initial value sensitivity when in chaotic state?
(2) How does consumer service sensitivity and service cost parameter influence the stable region of the system?
(3) In a long-term game, how does the adjustment speed of service influence profits of the third-party collector and the manufacturer?
(4) How can the chaotic phenomena be controlled if chaos occurs?
The rest of the paper is organized as follows. After a review of relevant literatures in Section 2, the assumptions and notations are described in Section 3. Then, we analyze the equilibrium point of the game model and local stability range of the equilibrium point in Section 4. In Section 5, a numerical study is given to analyze the complexity of the market. A method of chaos control is given in Section 6, and in Section 7 we draw the relevant conclusions.

2. Literature Review

This section mainly reviews literature relevant to this study. Three main streams of research are related to our work: research on supply chain with dual-channel recycling structure, pricing or service decisions in supply chain, and dynamic game with bounded rationality in supply chain.
There are a growing number of research papers on dual-channel recycling supply chain in recent years. Different structures of recycling modes can affect the operation of the supply chain. Gong [8] analyzed optimal decisions under different recycling modes and channel power structures and found that hybrid recycling is more effective. Ranjbar [16] studied optimal pricing and recycling strategies under the competitive recycling of the retailer and the third party, and found that the overall recycling rate of the supply chain is the highest under the leadership of retailers. He [17] established a closed-loop supply chain model with competitive recycling and proposed two coordination mechanisms in the case of centralized decision-making and decentralized decision-making. Zhao et al. [18] established two models considering the competition in collecting channels, and analyzed the optimal choice of the collecting channel as well as the optimal operational decisions. Liu [19] considered recycling innovation input for the collector and found that the manufacturer prefers the hybrid recycling mode. Moreover, the collection rate was the highest under the hybrid model compared with other cases.
Besides the research on dual-channel recycling supply chain, the service decision in supply chain has attracted extensive attention from academia. With the improvement of people’s living standards, people’s sensitivity to services has also become higher. Dumrongsiri [20] thought that service has become an important factor which affects the behavior of the customer. Zhou [21] analyzed the impact of free riding on supply chain members’ service strategies and profits in the case of differential pricing and non-differential pricing. Yan and Pei [22] found that improved retail services effectively alleviated the channel competition and conflict and improved the supply chain performance in a competitive market. Wu [23] obtained optimal service and pricing decisions in a dual-channel reverse supply chain by using the Stackelberg model, and found that the improvement of service level can maximize the profit of the recycling center. In this paper, we discussed the service decision of the third-party collector and customers’ service sensitivity was also considered in the model.
Most of the above literature assumes that members in supply chains are completely rational. When considering the bounded rational of the members in supply chains, the optimal decisions tend to be more complicated. Some studies focused on the dynamic duopoly games of Bertrand and Cournot [24,25,26]. Elsadany [27] established a bounded rational duopoly game with delay to analyze its complexity and the local stability of the equilibrium points. Xue [28] established a dynamic model to identify the effect of retailer-led recycling and how contracts affect different retail modes. Li [29] studied a recycling price game model in the long-term competition and analyzed the stability of the system. Furthermore, he considered the fair concern of a recycler in a competitive reverse supply chain and analyzed how fair concern influenced the stable area [30].
By reviewing literature, although some aforementioned literature have drawn valuable conclusions in above three streams, most of the literature [10,12,14,21] relevant to the service focused on the service in the forward supply chain. While Wu [23] and Kong [31] considered the level of service in the reverse supply chain, they ignored the bounded rationality of the supply chain member. Thus, the research concerned with service decision and pricing decision in reverse supply chain under bounded rationality is valuable. Based on the above discussion, from the perspective of bounded rationality, this study considered the dynamic change of the pricing and service decision under the mode of joint recycling by the manufacturer and third-party collector. The present paper contributes to the literature in three aspects. First, to the best of our knowledge, this is the first research to consider the dynamic change of the third-party collector’s service decision in the reverse supply chain. Second, we focus on the dynamic pricing and service decisions in the process of recycling and analyze the stability of the recycling market under different parameters. Third, this research enriches the theory of the recycling market with a dual-channel recycling structure.

3. Model Notation and Assumptions

In this section, a game model which describes the process of recycling was established. As shown in Figure 1, a third-party collector and a manufacturer were contained in the model. Collection is done by both the manufacturer and the third-party collector. In this model, the manufacturer set up his own collecting channel to collect waste products from consumers at price P 2 , and the third-party collector collected waste products from consumers at price P 1 . Different from the way of collection by remanufacturers and retailers, professional third-party collectors also pay attention to service provided to the customers in order to obtain more waste products. Therefore, we assumed the service level provided by the third-party collector is s , and the cost for the third-party collector to achieve service level s is σ s 2 / 2 , where σ is the parameter of the service cost [32,33,34].
The notations in this study are listed below:
  • P 1 : Recycling price of the third-party collector per unit.
  • P 2 : Recycling price of the manufacturer per unit.
  • P m : The subsidy given by the manufacturer to the third-party collector per unit.
  • s : The service level of the third-party collector.
  • k 1 : The recycling quantity when the third-party collector’s recycling price is 0.
  • k 2 : The recycling quantity when the manufacturer’s recycling price is 0.
  • C 1 /   C 2 : The collection cost of the third-party collector/manufacturer per unit.
  • M : The production cost saved by waste product per unit.
  • π 3 p : The profit of the third-party collector.
  • π m : The profit of the manufacturer.
  • G 1 : The recycling quantity function of the third-party collector.
  • G 2 : The recycling quantity function of the manufacturer.
The following are some assumptions that need to be illustrated specifically.
Hypothesis 1 (H1).
The recycling quantity function is a linear function of the recovery price and service level.
The remanufacturer and the third-party collector can be regarded as duopoly in the recycling market; they decide recycling price and level of service at the same time. According to the study of Huang and Wang [35], we assumed that the recycling quantity function in the market was a linear function of the recycling price and there had price competition between the third-party collector and the manufacturer. In addition, following the literature [31,36,37], the improvement of service level has a positive effect on the recycling quantity function. Therefore, recycling quantity functions of the third-party collector and manufacturer are:
G 1 = k 1 + α P 1 + γ s β P 2
G 2 = k 2 + α P 2 β P 1
Here, the parameter α represents the consumers’ sensitivity to recycling price, the parameter β represents the price cross parameter between response channels, and the parameter γ represents consumers’ sensitivity to the service level of the third-party collector. We assumed α > β , which means the channel’s own price has a greater impact on the recycling quantity function.
Hypothesis 2 (H2).
We assume that all waste products collected by the manufacturer and the third-party collector are valuable. To ensure the collection and remanufacturing can bring benefit to both manufacturer and third-party collection, similar to study of Savaskan [6], the condition M > M a x { P 2 +   C 2 , P m } is always true. Respectively, for simplification of the analysis, we assume C 1 =   C 2 = c = 0 [30].
Hypothesis 3 (H3).
Neither the manufacturer nor the third-party collector can totally know the market information, and they are bounded rational.

4. Dynamic Model Construction and Analyzation

According to the above assumptions, the profit functions of the third-party collector and the manufacturer are:
π 3 p = G 1 ( P m P 1   C 1 ) σ s 2 2 = ( k 1 + α P 1 + γ s β P 2 ) ( P m P 1 ) σ s 2 2
π m = G 1 ( M P m ) + + G 2 ( M P 2   C 2 ) = ( k 1 + α P 1 + γ s β P 2 ) ( M P m ) + ( k 2 + α P 2 β P 1 ) ( M P 2 )
In the profit functions, the decision variables of the third-party collector are the recycling price P 1 and service level s . This paper is aimed to study the dynamic changes in the process of collection,   P m is regarded as a constant [30]. Therefore, the decision variable of the manufacturer is the recycling price P 2 .

4.1. Methodology

In the recycling market, the third-party collector and the manufacturer can be regarded as the duopoly in the reverse supply chain, and the game model is similar to the Bertrand model. In the dynamic process, let P 1 ( t ) , P 2 ( t ) , and s ( t ) represent recycling prices and service level during period t = 0 , 1 , 2 Then, profits of the third-party collector and the manufacturer in a single period are Equations (3) and (4). According to the literature [38,39], the recycling price and service level for period t + 1 are decided by solving the optimization problem:
{ P 1 ( t + 1 ) = a r g m a x π 3 p ( P 1 ( t ) , P 2 e ( t + 1 ) , s ( t ) ) s ( t + 1 ) = a r g m a x π 3 p ( P 1 ( t ) , P 2 e ( t + 1 ) , s ( t ) ) P 2 ( t + 1 ) = a r g m a x π m ( P 1 e ( t + 1 ) , P 2 ( t ) , s e ( t + 1 ) )
where P 2 e ( t + 1 ) represents the expectation of the third-party collector about the recycling price of the manufacturer, P 1 e ( t + 1 ) represents the expectation of the manufacturer about the recycling price of the third-party collector, and s e ( t + 1 ) represents the expectation of the manufacturer about the service level of the third-party collector.
Due to the bounded rationality and uncertainty of the recycling market, the manufacturer and third-party collector cannot completely know the market information. Thus, they will make decisions according to the local estimation of marginal profit. Such decision-making based on the bounded rationality is common in practice [40,41]. If the marginal profit in period (t) is positive, the manufacturer or the third-party collector will correspondingly increase the recycling price or service level in period (t + 1). Otherwise, the manufacturer or the third-party collector will decrease the recycling price or service level.

4.2. Model Construction

According to the above profit functions and methodology, the model can be constructed as follows.
{ P 1 ( t + 1 ) = P 1 ( t ) + g 1 P 1 ( t ) π 3 p P 1 s ( t + 1 ) = s ( t ) + g 2 s ( t ) π 3 p s P 2 ( t + 1 ) = P 2 ( t ) + g 3 P 2 ( t ) π m P 2
System (6) gives the dynamic price and service decision, and the decision variables directly relate to the positive parameter g i (i =1, 2, 3), which represents the price and service adjustment speed. By solving the first-order condition π 3 p P 1 ,   π 3 p s ,   π m P 2 , the marginal profit functions are:
{ π 3 p P 1 = α ( P m P 1 ) ( k 1 + α P 1 + γ s β P 2 ) π 3 p s = γ ( P m P 1 ) σ s π m P 2 = α ( M P 2 ) β ( M P m ) ( k 2 + α P 2 β P 1 )
When substituting expression (7) into (6), we gain the three-dimensional discrete dynamical system (8):
{ P 1 ( t + 1 ) = P 1 ( t ) + g 1 P 1 ( t ) [ α ( P m P 1 ) ( k 1 + α P 1 + γ s β P 2 ) ] s ( t + 1 ) = s ( t ) + g 2 s ( t ) [ γ ( P m P 1 ) σ s ]   P 2 ( t + 1 ) = P 2 ( t ) + g 3 P 2 ( t ) [ α ( M P 2 ) β ( M P m ) ( k 2 + α P 2 β P 1 ) ]

4.3. Equilibrium Points and Local Stability

According to the definition of the equilibrium points, in order to guarantee the practical significance, we assume that all the equilibrium points are non-negative. By setting P 1 ( t + 1 ) = P 1 ( t ) , s ( t + 1 ) = s ( t ) , and P 2 ( t + 1 ) = P 2 ( t ) , we obtained the following eight equilibrium points:
R 1 ( 0 , 0 , 0 ) ,     R 2 ( 0 , γ P m σ , 0 ) ,   R 3 ( 0 , 0 , H 2 α ) ,   R 4 ( 0 ,   γ P m σ ,   H 2 α ) ,
R 5 ( α P m k 1 2 α , 0 , 0 ) ,   R 6 ( α σ P m γ 2 P m σ k 1 2 α σ γ 2 , γ α P m + γ k 1 2 α σ γ 2 , 0 ) ,   R 7 ( T 4 α 2 β 2 ,   0 , N 4 α 2 β 2 ) ,
R 8 ( σ T 2 α γ 2 P m F ,   γ ( 4 α 2 P m β 2 P m T ) F ,   σ N γ 2 H β γ 2 P m F )
where   { H = M ( α β ) + β P m k 2 T = M ( α β β 2 ) + ( 2 α 2 + β 2 ) P m 2 α k 1 β k 2 N = 2 M α ( α β ) + 3 α β P m 2 α k 2 β k 1 F = 4 α 2 σ 2 α γ 2 β 2 σ
Obviously, R i (i = 1,2,3,4,5,6,7) are boundary equilibrium points. When F > 0 ,   σ T 2 α γ 2 P m > 0 , 4 α 2 P m β 2 P m T > 0 , and σ N γ 2 H β γ 2 P m > 0 ,   R 8 is the Nash equilibrium point.
In order to analyze the stability of the equilibrium point, we obtain the Jacobian matrix of system (8):
J = [ J 11 γ g 1 P 1 β g 1 P 1 γ g 2 s J 22 0 β g 3 P 2 0 J 33 ]
where   { J 11 = 1 + g 1 [ α ( P m 3 P 1 ) ( k 1 + α P 1 + γ s β P 2 ) ] J 22 = 1 + g 2 [ γ ( P m P 1 ) 2 σ s ]   J 33 = 1 + g 3 [ α ( M 3 P 2 ) β ( M P m ) ( k 2 + α P 2 β P 1 ) ]
Theorem 1.
The boundary equilibrium points, R i (i = 1,2,3,4,5,6,7) of system (8) are unstable equilibrium points.
Proof. 
In order to prove this result, we find the eigenvalues of the Jacobian matrix at each boundary equilibrium points R i (i = 1,2,3,4,5,6,7). The Jacobian matrix at R 1 is:
J ( R 1 ) = [ 1 + g 1 [ α P m k 1 ] 0 0 0 1 + g 2 γ P m 0 0 0 1 + g 3 [ α M β ( M P m ) k 2 ] ]
One eigenvalue of J ( R 1 ) is: λ 2 = 1 + g 2 γ P m > 1 , therefore, R 1 is an unstable equilibrium point.
The Jacobian matrix at R 2 is:
J ( R 2 ) = [ 1 + g 1 [ α P m ( k 1 + γ 2 P m σ ) ] 0 0 0 1 g 2 γ P m 0 0 0 1 + g 3 H ]
One eigenvalue of J ( R 2 ) is: λ 3 = 1 + g 3 H , since H > 0 , we obtain | λ 3 | > 1 , R 2 is an unstable equilibrium point.
The Jacobian matrix at R 3 is:
J ( R 3 ) = [ 1 + g 1 [ α P m k 1 + β H 2 a ] 0 0 0 1 + g 2 γ P m 0 β g 3 H 2 a 0 1 + g 3 [ α M β ( M P m ) k 2 2 H ] ]
One eigenvalue of   J ( R 3 ) is: λ 1 = 1 + g 1 [ α P m k 1 + β H 2 a ] , since α P m k 1 > 0 ,   H > 0 , we obtain | λ 1 | > 1 ,   R 3 is an unstable equilibrium point.
Similarly, other boundary equilibrium points can be proved to be unstable equilibrium points. □
Theorem 2.
If the system parameters satisfy F > 0 ,   σ T 2 α γ 2 P m > 0 , 4 α 2 P m β 2 P m T > 0 , σ N γ 2 H β γ 2 P m > 0 and the Jury conditions, the Nash equilibrium point R 8 is locally asymptotically stable.
Proof. 
In contrast to the boundary equilibrium, the stability of the Nash equilibrium point is difficult to determine by the eigenvalue because of the complexity of the calculation, but the local stability of the Nash equilibrium point can be obtained by using the Jury condition [42]. By substituting R 8 into (9), we can obtain:
J ( R 8 ) = [ J 11 γ g 1 P 1 * β g 1 P 1 * γ g 2 s * J 22 0 β g 3 P 2 * 0 J 33 ]
where   { J 11 = 1 + g 1 [ α ( P m 3 P 1 * ) ( k 1 + α P 1 * + α s * β P 2 * ) ] J 22 = 1 + g 2 [ γ ( P m P 1 * ) 2 σ s * ] J 33 = 1 + g 3 [ α ( M 3 P 2 * ) β ( M P m ) ( k 2 + α P 2 * β P 1 * ) ]
We gain the characteristic equation from matrix (10); the characteristic equation is as follows:
λ 3 + a 2 λ 2 + a 1 λ + a 0 = 0
where   { a 2 = ( J 11 + J 22 + J 33 ) a 1 = J 22 J 33 + J 11 J 22 + + J 11 J 33 γ 2 g 1 P 1 * g 2 s * β 2 g 1 P 1 * g 3 P 2 * a 0 = J 11 J 22 J 33 + β 2 g 1 P 1 * g 3 P 2 * J 22 + γ 2 g 1 P 1 * g 2 s * J 33
By using the Jury condition, the sufficient and necessary condition for the Nash equilibrium point R 8 to be stable is:
{       1 + a 2 + a 1 + a 0 > 0 1 + a 2 a 1 + a 0 < 0       | a 0 | < 1 | a 0 2 1 | > | a 0 a 2 a 1 |
The stable region of the Nash equilibrium point R 8 is defined by (12). However, beyond the stability region more complex phenomena such as bifurcation and chaos will occur. Moreover, the local stability of the system in R 8 can be affected by consumers’ service sensitivity parameter γ and service cost parameter σ . Based on inequalities (12), the three-dimensional stability domains of the system (8) are simulated when γ and σ take different values (as shown in Section 5.3). □

5. Numerical Simulation

In this section, numerical simulations were carried out to show the dynamic change in the system (8). The dynamic characteristics of the system will be analyzed by bifurcation diagram, maximum Lyapunov exponent, and strange attractor. Both the manufacturer and third-party collector can observe the influence of the adjustment on the system (8). The main purpose of the numerical simulation is to analyze the local stability of the equilibrium point and simulate the dynamic evolution of the pricing and service decisions under different adjustment parameters. Similar to the numerical simulation section in literatures [29,30,39,43], the values of parameter just need to meet the model assumptions. Hence, we set the parameter values as follows: α = 0.6 , β = 0.3 , γ = 3 , M = 40 , k 1 = 5 , k 2 = 3 ,   σ = 50 , P m = 30 ; these data were close to the real-world condition. According to these parameters, we can obtain the Nash equilibrium point R 8 is (12.8,1.03,18.2).

5.1. The Impact of Adjustment Parameters on the Stability of the System

When the adjustment parameter of the third-party collector g 1 changed in the range of   [ 0.02 , 0.18 ] , Figure 2 showed the bifurcation diagrams of P 1 and P 2 in (a) and s in (b) with respect to g 2 = 0.01 and g 3 = 0.01 . When g 1 < 0.1255 , the system (8) was in the stable region, the Nash equilibrium point was locally stable. When g 1 = 0.1255 , the system (8) will deviate from the stable region and the phenomenon of bifurcation appears. Corresponding to the maximum Lyapunov exponent in Figure 3, we found that the maximum Lyapunov exponent of the system (8) reached the first peak when the bifurcation occurred. When g 1   increased to the value of 0.157, the second bifurcation occurred in Figure 2 (corresponding to the second peak in Figure 3). When g 1 > 0.1639 , the maximum Lyapunov in Figure 3 was larger than 0, and system (8) occurred chaos and appeared complex dynamic states.
Figure 4 showed the bifurcation diagrams of P 1 and P 2 in (a) and s in (b) with respect to adjustment parameter g 2 while g 1 = 0.1 , g 3 = 0.01 . When the parameter g 2 belonged to the stable region ( g 2 [ 0.01 , 0.0252 ] ) , the Nash equilibrium point was locally stable. Corresponding to the maximum Lyapunov exponent in Figure 5, we found the maximum Lyapunov exponent of the system (8) reached the first peak when g 2 = 0.0252 . Then, the system (8) will deviate from the stable region and the first bifurcation point appeared. It was followed by periods 2, 3, 4… and gradually descended into chaos when g 2 > 0.0375 (the maximum Lyapunov exponent in Figure 5 is larger than 0).
Figure 6 showed the bifurcation diagrams of P 1 and P 2 in (a) and s in (b) with respect to adjustment parameter g 3 while g 1 = 0.1 , g 2 = 0.01 . When the parameter g 3 < 0.07 , the Nash equilibrium point was locally stable. It was found that the bifurcation point appeared when g 3 = 0.07 . Then, the system (8) will deviate from the stable region and gradually descended into chaos when g 3 > 0.1057 (the maximum Lyapunov exponent in Figure 7 is larger than 0).
Both the maximum Lyapunov exponent and the strange attractor are important indexes of the system dynamics. For example, Figure 8 showed the strange attractor under different value of g 1 corresponding to Figure 2. From Figure 8a, it can be found that when the system (7) was locally stable, the strange attractor will evolve to the Nash equilibrium point R 8 ( 12.8 , 1.03 , 18.2 ) . When g 1 > 0.1255 , like the situation in Figure 8b, the Nash equilibrium point begins to change into a branch state for the first time (the system appears the first bifurcation point corresponding to Figure 2). Figure 8c shows when g 1 > 0.1639 , the system (8) is in completely chaotic state and the stranger attractor is an irregular shape.
Through the above numerical simulation analysis, it can be concluded that the adjustment speeds g 1 ,   g 2 , and g 3 had great effect on the stability of system (8). If the manufacturer and third-party collector adjusted their price or service strategies too fast, the recycling market will be in a chaotic state. Once chaos occurred in the economic system, it was difficult for the game player to make a correct prediction. Furthermore, by comparing bifurcation diagrams of the system under different adjustment parameters, it was found that even if the service adjustment parameter of the third-party collector changed in a very small range, chaos may occur in the recycling market. Therefore, the system (8) is most sensitive to the change of the third-party’s service decision. The third-party collector should not adjust his service decision too quickly.

5.2. Sensitivity Analysis of Initial State in the System

In this section, we explored the sensitivity of initial value in the system (8). First, we compared the slight difference of initial value when the system (8) is in a stable state. From Section 4.1, it can be seen that when g 1 = 0.1 , g 2 = 0.01 , and   g 3 = 0.01 , the system was locally stable. Here we took a group of initial values as baseline data, i.e., P 1 = 11 , s = 0.7 , and   P 2 = 15 . Then, we added value of 0.001 to the three initial values separately. The difference in recycling prices and service after 30 iterations is shown in Figure 9. It was found that when the initial value of the baseline data was changed by 0.001, the system tended to a slight difference at the beginning. As iterations increased, the difference gradually became 0 and the value was equal to the Nash equilibrium point. It showed that even if there was a little deviation in the initial value when system (8) in the stable state, with the increase of iterations, the influence brought by the deviation could be eliminated.
Then, we studied the sensitivity of the initial value when the system (8) was in chaotic state. Based on Section 4.1, when g 1 = 0.18 , g 2 = 0.01 , and g 3 = 0.01 , the system (8) appeared chaotic. We still took the same baseline data, i.e., P 1 = 11 , s = 0.7 , and   P 2 = 15 . Figure 10 shows the sensitivity of the system (8) to initial values when in chaotic state. It was obvious that the value of recycling prices and service level in the system (8) tended to be significantly different over time. We can conclude that the chaotic system was highly sensitive to the initial value. Even if the initial value changed slightly, the result would change a lot. Therefore, the third-party collector and manufacturer should pay more attention to choosing the initial price and service strategies when the recycling market is in chaotic state.

5.3. The Impact of Consumers’ Service Sensitivity and Service Cost Parameter on the Stable Region of the System

Different from former literature, this study focused on the dynamic evolution of the third-party collector’s service decision in reverse supply chain. From Section 4.1, we found that the adjustment parameter of service decision had a great effect on the stability of the recycling market. Actually, consumers’ service sensitivity is also an important factor which influences the service decision. From the recycling quantity function of the third-party collector in Assumption 1, it is clear that the more sensitive consumers are to the service, the more quantity of waste products the third-party collector can collect. Hence, the value of consumers’ service sensitivity parameter also plays a significant role in the recycling market.
First, we analyzed the impact of consumers’ service sensitivity on the stable region of the system (8). Inequalities (11) showed the stable region of the system (8) under adjustment speed parameters g 1 , g 2 , and g 3 . Here, we took three different values of the service sensitivity parameter, i.e., γ = 1 , γ = 3 , and γ = 5 . We divided consumer service sensitivity into three levels corresponding to the different value of parameter γ , i.e., weak, middle, and strong. When other parameters are fixed, the stable region for the Nash equilibrium point under different values of γ are shown in Figure 11. From the stable region in Figure 11a–c, with the increase of parameter γ , the stable region of the system (8) decreased. Therefore, from the perspective of the recycling market, the increase in consumers’ service sensitivity is not conducive to stability of the market. By comparing the stable range of adjustment speed parameters under different parameter γ , we found that the increase of consumers’ service sensitivity can expand the stable range of adjustment parameters g 1 and g 3 while narrowing the stable range of the adjustment parameter g 2 . Therefore, if consumers have strong sensitivity to a third-party collector’s service, the third-party collector should maintain his service decision because the change in service decision may cause great fluctuations in the recycling market.
The Nash equilibrium points corresponding to Figure 11 are as follows: R 8 ( 15.29 ,   0.29 ,   18.82 ) in Figure 11a, R 8 ( 12.8 ,   1.03 ,   18.2 ) in Figure 11b, and R 8 ( 4 ,   2.6 ,   16 ) in Figure 11c. By analyzing the changes in the value of three Nash equilibrium points, we found that when consumers’ service sensitivity increased, the recycling price of the manufacturer and the third-party collector decreased while the service level increased. The result indicates that when consumers have a strong sensitivity to the service provided by the third-party collector, in order to obtain more waste products from consumers, the third-party collector tends to improve his service level. Under this situation, the third-party collector does not need to set a high recycling price to attract consumers, and the third-party collector will reduce his recycling price. In the process of collection, there occurs price competition between the manufacturer and the third-party collector. When the recycling price of the third-party decreases, the price competition from the third-party becomes weak and the manufacturer also reduces his recycling price.
Besides the parameter of consumer service sensitivity, the value of the service cost parameter σ can affect the stable region of the system (8) at the same time. When other parameters are fixed, the stable region for the Nash equilibrium point under different values of σ are shown in Figure 12. The stable region of the system (8) decreased when the parameter   σ increased. By comparing the stable range of adjustment speed parameters under different value of parameter σ , we found that the increase of service cost parameter can narrow the stable range of adjustment parameters g 1 and g 3 while expanding the stable range of the adjustment parameter g 2 . Furthermore, the Nash equilibrium point corresponding to Figure 12 changed from R 8 ( 12.8 ,   1.03 ,   18.2 ) in Figure 12a to R 8 ( 14.3 ,   0.47 ,   18.57 ) in Figure 12b. The result indicates that the third-party collector will reduce the types of services provided to consumers if service cost is high. In this situation, both the manufacturer and the third-party collector will increase their recycling price to obtain more waste products. Since the value of the service cost parameter is dependent on the third-party collector, the third-party collector can take management, innovation, and other means to reduce the service cost and make the recycling market more stable.

5.4. The Impact of Service Adjustment Parameter on the Profit in Supply Chain

From the previous analysis, we found that the value of adjustment parameters, consumer service sensitivity, initial value, and service cost parameter are crucial to the stability of the market. In this section, we analyzed how the increase of the service adjustment speed influences game players’ profits. When parameter g 1 changes, Figure 13a,b showed the dynamic change of game players’ profit when g 2 = 0.01   and   g 3 = 0.01 , and Figure 13c,d showed the dynamic change of game players’ profit when g 2 = 0.02   and   g 3 = 0.01 . The result showed that no matter what the service adjustment parameter g 2 is, when the market is in a stable state, the profits of the manufacturer and third-party collector depend on the value of the Nash equilibrium point R 8 . With increase of the adjustment parameter g 1 , the profit of the manufacturer and third-party collector began to diverge and gradually become a chaotic state. Furthermore, when the adjustment parameter of service increased to 0.2, the bifurcation point of Figure 13a,b shifted to the left in Figure 13c,d. The result shows that there was a reverse relationship between parameters g 1 and g 2 . According to the recycling quantity function in Assumption 1, the increase change of service leads to the great changes in recycling quantity of the third-party collector. In order to maintain profit, the third-party collector tends to accelerate its recycling price adjustment. As a result, it may lead the market into the chaotic state more quickly. Hence, when the third-party collector adjusted service decision more quickly, the collector should control recycling price from a long-term perspective and take scientific methods to avoid chaos.
When parameter g 3 changes, Figure 14a,b shows the dynamic change of game players’ profit when g 1 = 0.1   and   g 2 = 0.01 , and Figure 13c,d shows the dynamic change of game players’ profit when g 1 = 0.1   and   g 2 = 0.02 . The result shows that when the system (8) is in stable state, the profits of the manufacturer and third-party collector depend on the value of the Nash equilibrium point R 8 . With the increase of the adjustment parameter g 3 , the profit of the manufacturer and third-party collector began to diverge gradually toward a chaotic state. Furthermore, when the adjustment parameter of service increased to 0.2, the bifurcation point of Figure 14a,b shifted to the left in Figure 14c,d. There was also a reverse relationship between the parameters g 2 and g 3 . Therefore, when the service adjustment parameter of the third-party collector increases, the stable range of adjustment parameters g 1 and g 3 will narrow. Once chaos appears in the market, the profit of the manufacturer and the third-party collector are difficult to predict.

5.5. Discussion

Under bounded rationality, the pricing and service decisions in the reverse supply chain tend to be complicated in the dynamic evolution. During the decision-making process, the participants adjust their decisions according to the marginal profit. According to the numerical simulation, the recycling market is always stable when the adjustment parameters meet the inequalities (12). Once the speed of adjustment is too fast, the recycling market will get out of the stable state and gradually become chaos. The chaotic system has strong initial sensitivity, which is similar to the conclusion of Zhang [39] and Xi [44]. Different from Li’s study [29], besides the pricing decision, the third-party collectors’ service decisions also affect the stability of the recycling market. Furthermore, the market is more sensitive to the change in service decision compared with the pricing decision. In addition, the stable area of the system (8) has strong connection with the consumers’ service sensitivity and the service cost parameter. The stable region of the system (8) decreases when these two parameters increase. By analyzing the change in the Nash equilibrium point, the increase of consumer service sensitivity promotes the improvement in third-party collectors’ service and weakens the recycling price competition between the manufacturer and the third-party collector. On the contrary, the increase of service cost parameter reduces the service level of the third-party and strengthens the recycling price competition between the manufacturer and the third-party collector, both of whom need to increase recycling price to obtain more waste products.
When discussing the profit situation in the reverse supply chain, we found an interesting result that there was a reverse relationship between the adjustment parameters of pricing decision and the adjustment parameter of service level, the increase of service adjustment will lead to the recycling market loss of control and the occurrence of chaos. In a chaotic state, the profit of the manufacturer and the third-party collector become disordered. The chaotic state is harmful to the development of the manufacturer and the third-party. In the next section, we will find measures to avoid chaos.

6. Chaos Control

Chaos is harmful to an economic system. When the market is in a chaotic state, members of the system are unable to accurately analyze the current market and make right price and service decisions. Therefore, certain measures should be taken to prevent and avoid the phenomenon of chaos. There are some methods to avoid chaos, for example, Zhang [43] used the delayed feedback control method to control the chaos in a nitrogen emissions system. In this paper, we used the variable feedback control method [45] to control the recycling market.
We expressed the discrete dynamic system (8) without control as:
{ P 1 ( t + 1 ) = g 1 ( P 1 ( t ) , s ( t ) , P 2 ( t ) ) s ( t + 1 ) = g 2 ( P 1 ( t ) , s ( t ) , P 2 ( t ) ) P 2 ( t + 1 ) = g 3 ( P 1 ( t ) , s ( t ) , P 2 ( t ) )
By adding control parameter L into system (13), the new system can be expressed as:
{ P 1 ( t + 1 ) = ( 1 L ) g 1 ( P 1 ( t ) , s ( t ) , P 2 ( t ) ) + L P 1 ( t ) s ( t + 1 ) = ( 1 L ) g 2 ( P 1 ( t ) , s ( t ) , P 2 ( t ) ) + L s ( t ) P 2 ( t + 1 ) = ( 1 L ) g 3 ( P 1 ( t ) , s ( t ) , P 2 ( t ) ) + L P 2 ( t )
From Section 5.1, it is clear that when g 1 = 0.18 , g 2 = 0.01 , and   g 3 = 0.01 , the system (8) becomes chaos. By using the variable feedback control method, Figure 15 shows the bifurcation diagram for control parameter L . As is shown in Figure 15, L 0 = 0.285 is corresponding to the bifurcation point of the system (14). When L < L 0 , the system (14) is still in the state of bifurcation or chaos. The smaller L is, the less stable the system (14) is. When L = 0 , the system (14) is completely in chaos. When L > L 0 , the system (14) tends to be stable.
In order to further illustrate the effectiveness of the control method, we drew the strange attractors of system (14) under the conditions of L = 0.1 and L = 0.4 . As is shown in Figure 16a, when L = 0.1 < L 0 , system (14) was still in a chaotic state. From Figure 16b, when L = 0.4 > L 0 , system (13) was effectively controlled, and the strange attractor will evolve to the Nash equilibrium point R 8 . The control parameter L can be understood as an internal or external factor. Internal factors can be understood as the learning ability of enterprises. In order to strengthen the control of chaos, the enterprise can predict the market information through big data analysis. External factors can be understood when enterprises lack learning ability, the government enters the recycling market by means of policy regulation to maintain the market order.

7. Conclusions

This study established a game model with dual-channel recycling structure. With the rapid development of the independent third-party collector, we considered the situation that the manufacturer and the third-party collector collect waste products at the same time. Different from manufacturers and retailers, third-party collectors always provide different service for consumers. Hence, we considered the service decision of the third-party collector. Based on bounded rationality, we analyzed the dynamic change of the service and price decisions in the recycling market. Then, we analyzed equilibrium points of the system and stable condition for the Nash equilibrium point. In the section of numerical simulation, we analyzed the dynamic characteristics of the system. Finally, we used the feedback control method to control the chaotic market.
The conclusions are as follows. (1) The adjustment speeds of the third-party and the manufacturer have great effect on the stability of the system. If the adjustment parameter is too large, the system will evolve to chaotic state. Furthermore, the market is most sensitive to the change of third-party collectors’ service decisions. (2) When the recycling market is in a stable state, the slight deviation of the initial value will be eliminated by increasing iterations. When the market is in a chaotic state, the small deviation of the initial value will be gradually amplified with the increase of iterations. The chaotic system has a strong sensitivity to the initial value. (3) The stable region of the dynamic system shrinks with the increase of consumer service sensitivity and the service cost parameter. Furthermore, the increase of consumer service sensitivity can expand the stable range of price adjustment parameters while narrowing the stable range of the service adjustment parameter. On the contrary, the increase of the service cost parameter can narrow the stable range of price adjustment parameters while expanding the stable range of the service adjustment parameter. (4) There was a reverse relationship between price adjustment parameters and the service adjustment parameter. With the increase of service adjustment parameter, the bifurcation point of the system will appear more quickly. (4) Chaotic systems can be controlled by the feedback control method.

Author Contributions

J.Z.: Supervision, methodology, funding acquisition, writing—original draft. C.Z.: Validation, conceptualization, writing—review and editing, software. J.Z. and C.Z. should be regarded as the co-first authors. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Fundamental Research Funds for the Central Universities (Grant No. NS2019045) and the National Natural Science Foundation of China (No. 71101071).

Acknowledgments

The authors would like to express our sincere appreciation for the valuable comments made by anonymous reviewers. The authors also thank editing service provider for editing this manuscript.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Guide, V.D.R.; Van Wassenhove, L.N. The Evolution of Closed-Loop Supply Chain Research. Oper. Res. 2009, 57, 10–18. [Google Scholar] [CrossRef]
  2. Wu, H.; Han, X.; Yang, Q.; Pu, X. Production and coordination decisions in a closed-loop supply chain with remanufacturing cost disruptions when retailers compete. J. Intell. Manuf. 2015, 29, 227–235. [Google Scholar] [CrossRef]
  3. Nnorom, I.C.; Osibanjo, O. Overview of prospects in adopting remanufacturing of end-of-life electronic products in the developing countries. Int. J. Innov. Manag. Technol. 2010, 1, 328–338. [Google Scholar]
  4. Chen, J.; Chen, B. Competing with customer returns policies. Int. J. Prod. Res. 2016, 54, 2093–2107. [Google Scholar] [CrossRef]
  5. Zhu, C.; Zhang, J. Study on pricing and coordination of closed-loop supply chain considering retailer’s dual-channel and differences in recycled products quality. Comput. Integr. Manuf. Syst. 2020, 2020, 1–19. [Google Scholar]
  6. Savaskan, R.C.; Bhattacharya, S.; Van Wassenhove, L.N. Closed-Loop Supply Chain Models with Product Remanufacturing. Manag. Sci. 2004, 50, 239–252. [Google Scholar] [CrossRef] [Green Version]
  7. Gong, Y.; Chen, M.; Zhuang, Y. Decision-Making and Performance Analysis of Closed-Loop Supply Chain under Different Recycling Modes and Channel Power Structures. Sustainability 2019, 11, 6413. [Google Scholar] [CrossRef] [Green Version]
  8. Wang, N.; Yang, S.; He, Q.; Jia, T. Competitive dual-collecting regarding consumer behavior and coordination in closed-loop supply chain. Comput. Ind. Eng. 2020, 144, 106481. [Google Scholar] [CrossRef]
  9. Lightfoot, H.; Baines, T.; Smart, P. The servitization of manufacturing: A systematic literature review of interdependent trends. Int. J. Oper. Prod. Manag. 2013, 33, 1408–1434. [Google Scholar] [CrossRef] [Green Version]
  10. Xia, Y.; Gilbert, S.M. Strategic interactions between channel structure and demand enhancing services. Eur. J. Oper. Res. 2007, 181, 252–265. [Google Scholar] [CrossRef]
  11. Chen, X.; Li, L.; Zhou, M. Manufacturer’s pricing strategy for supply chain with warranty period-dependent demand. Omega 2012, 40, 807–816. [Google Scholar] [CrossRef]
  12. Yang, D.; Xiao, T. Coordination of a supply chain with loss-averse consumers in service quality. Int. J. Prod. Res. 2017, 55, 3411–3430. [Google Scholar] [CrossRef]
  13. Ren, L.; He, Y.; Song, H. Price and Service Competition of Dual-Channel Supply Chain with Consumer Returns. Discret. Dyn. Nat. Soc. 2014, 11, 1–10. [Google Scholar] [CrossRef] [Green Version]
  14. Dan, B.; Lou, Y.; Ma, S. The pricing and optimization strategy of product service supply chain for service promotion. Manag. Rev. 2017, 29, 211–222. [Google Scholar]
  15. Rand, D. Exotic phenomena in games and duopoly models. J. Math. Econ. 1978, 5, 173–184. [Google Scholar] [CrossRef]
  16. Ranjbar, S.; Sahebi, H.; Ashayeri, J.; Teymouri, A. A competitive dual recycling channel in a three-level closed loop supply chain under different power structures: Pricing and collecting decisions. J. Clean. Prod. 2020, 272, 122623. [Google Scholar] [CrossRef]
  17. He, Q.; Wang, N.; Yang, Z.; He, Z.; Jiang, B.; Yang, Z. Competitive collection under channel inconvenience in closed-loop supply chain. Eur. J. Oper. Res. 2019, 275, 155–166. [Google Scholar] [CrossRef]
  18. Zhao, J.; Wei, J.; Li, M. Collecting channel choice and optimal decisions on pricing and collecting in a remanufacturing supply chain. J. Clean. Prod. 2017, 167, 530–544. [Google Scholar] [CrossRef]
  19. Liu, R.; Du, B.; Yuan, W.; Li, G. Competing Reverse Channels’ Performance with Sustainable Recycle Innovation Input. Appl. Sci. 2020, 10, 5429. [Google Scholar] [CrossRef]
  20. Dumrongsiri, A.; Fan, M.; Jain, A.; Moinzadeh, K. A supply chain model with direct and retail channels. Eur. J. Oper. Res. 2008, 187, 691–718. [Google Scholar] [CrossRef]
  21. Zhou, Y.; Guo, J.; Zhou, W. Pricing/service strategies for a dual-channel supply chain with free riding and service-cost sharing. Int. J. Prod. Econ. 2018, 196, 198–210. [Google Scholar] [CrossRef]
  22. Yan, R.; Pei, Z. Retail services and firm profit in a dual-channel market. J. Retail. Consum. Serv. 2009, 16, 306–314. [Google Scholar] [CrossRef]
  23. Wu, D.; Chen, J.; Li, P.; Zhang, R. Contract coordination of dual channel reverse supply chain considering service level. J. Clean. Prod. 2020, 260, 121071. [Google Scholar] [CrossRef]
  24. Abbink, K.; Brandts, J. Pricing in Bertrand competition with increasing marginal costs. Games Econ. Behav. 2008, 63, 1–31. [Google Scholar] [CrossRef] [Green Version]
  25. Agiza, H.N.; Elsadany, A.A. Nonlinear dynamics in the Cournot duopoly game with heterogeneous players. Physica A 2003, 320, 512–524. [Google Scholar] [CrossRef] [Green Version]
  26. Agiza, H.N.; Hegazi, A.S.; Elsadany, A.A. Complex dynamics and synchronization of a duopoly game with bounded rationality. Math. Comput. Simul. 2002, 58, 133–146. [Google Scholar] [CrossRef]
  27. Elsadany, A.A. Dynamics of a delayed duopoly game with bounded rationality. Math. Comput. Model. 2010, 52, 1479–1489. [Google Scholar] [CrossRef]
  28. Xue, R.; Zhang, F.; Tian, F. A System Dynamics Model to Evaluate Effects of Retailer-Led Recycling Based on Dual Chains Competition: A Case of e-Waste in China. Sustainability 2018, 10, 3391. [Google Scholar] [CrossRef] [Green Version]
  29. Li, T.; Yan, D.; Ma, X. Stability Analysis and Chaos Control of Recycling Price Game Model for Manufacturers and Retailers. Complexity 2019, 2019, 3157407. [Google Scholar] [CrossRef]
  30. Li, T.; Yan, D.; Sui, S. Research on the Complexity of Game Model about Recovery Pricing in Reverse Supply Chain considering Fairness Concerns. Complexity 2020, 2020, 9621782. [Google Scholar]
  31. Kong, L.; Liu, Z.; Pan, Y.; Xie, J.; Yang, G. Pricing and service decision of dual-channel operations in an O2O closed-loop supply chain. Ind. Manag. Data Syst. 2017, 117, 1567–1588. [Google Scholar] [CrossRef]
  32. Tsay, A.A.; Agrawal, N. Channel Dynamics under Price and Service Competition. Manuf. Serv. Oper. Manag. 2000, 2, 372–391. [Google Scholar] [CrossRef]
  33. Iyer, G. Coordinating Channels under Price and Nonprice Competition. Mark. Sci. 1998, 17, 338–355. [Google Scholar] [CrossRef] [Green Version]
  34. Li, Q.; Li, M.; Huang, Y. Dynamic Investigation in Green Supply Chain considering Channel Service. Complexity 2020, 2020, 1640724. [Google Scholar]
  35. Huang, Y.; Wang, Z. Closed-loop supply chain models with product take-back and hybrid remanufacturing under technology licensing. J. Clean. Prod. 2017, 142, 3917–3927. [Google Scholar] [CrossRef]
  36. Pu, X.; Gong, L.; Han, X. Consumer free riding: Coordinating sales effort in a dual-channel supply chain. Electron. Commer. Res. Appl. 2017, 22, 1–12. [Google Scholar] [CrossRef]
  37. Li, Q.; Li, B. Dual-channel supply chain equilibrium problems regarding retail services and fairness concerns. Appl. Math. Model. 2016, 40, 7349–7367. [Google Scholar] [CrossRef]
  38. Agiza, H.N.; Elsadany, A.A. Chaotic dynamics in nonlinear duopoly game with heterogeneous players. Appl. Math. Comput. 2004, 149, 843–860. [Google Scholar] [CrossRef]
  39. Zhang, J.; Da, Q.; Wang, Y. The dynamics of Bertrand model with bounded rationality. Chaos Solitons Fractals 2009, 39, 2048–2055. [Google Scholar] [CrossRef]
  40. Ma, J.; Guo, Z. The parameter basin and complex of dynamic game with estimation and two-stage consideration. Appl. Math. Comput. 2014, 248, 131–142. [Google Scholar] [CrossRef]
  41. Zhao, L. Nonlinear Complex Dynamics of Carbon Emission Reduction Cournot Game with Bounded Rationality. Complexity 2017, 2017, 8301630. [Google Scholar] [CrossRef]
  42. Puu, T. Attractors, Bifurcations and Chaos: Nonlinear Phenomena in Economics; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  43. Zhang, J.; Xi, X. Complexity Analysis of a 3-Player Game with Bounded Rationality Participating in Nitrogen Emission Reduction. Complexity 2020, 2020, 2069614. [Google Scholar] [CrossRef]
  44. Xi, X.; Zhang, J. Complexity analysis of a decision-making game concerning governments and heterogeneous agricultural enterprises with bounded rationality. Chaossolitons Fractals 2020, 140, 110220. [Google Scholar] [CrossRef]
  45. Bischi, G.I.; Chiarella, C. Global Analysis of Dynamic Models in Economics and Finance; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
Figure 1. Supply chain structure.
Figure 1. Supply chain structure.
Sustainability 12 09348 g001
Figure 2. The bifurcation diagram of the system (8) when g 1 changes: (a) P 1 and P 2 ; (b)   s .
Figure 2. The bifurcation diagram of the system (8) when g 1 changes: (a) P 1 and P 2 ; (b)   s .
Sustainability 12 09348 g002
Figure 3. The maximum Lyapunov exponent when g 1 changes.
Figure 3. The maximum Lyapunov exponent when g 1 changes.
Sustainability 12 09348 g003
Figure 4. The bifurcation diagram of system (8) when g 2 changes: (a) P 1 and P 2 ; (b)   s .
Figure 4. The bifurcation diagram of system (8) when g 2 changes: (a) P 1 and P 2 ; (b)   s .
Sustainability 12 09348 g004
Figure 5. The maximum Lyapunov exponent when g 2 changes.
Figure 5. The maximum Lyapunov exponent when g 2 changes.
Sustainability 12 09348 g005
Figure 6. The bifurcation diagram of system (8) when g 3 changes: (a) P 1 and P 2 ; (b)   s .
Figure 6. The bifurcation diagram of system (8) when g 3 changes: (a) P 1 and P 2 ; (b)   s .
Sustainability 12 09348 g006
Figure 7. The maximum Lyapunov exponent when g 3 changes.
Figure 7. The maximum Lyapunov exponent when g 3 changes.
Sustainability 12 09348 g007
Figure 8. Strange attractor: (a) g 1 = 0.1 , g 2 = 0.01 , g 3 = 0.01 ; (b) g 1 = 0.14 , g 2 = 0.01 , g 3 = 0.01 ; (c) g 1 = 0.18 , g 2 = 0.01 , g 3 = 0.01 .
Figure 8. Strange attractor: (a) g 1 = 0.1 , g 2 = 0.01 , g 3 = 0.01 ; (b) g 1 = 0.14 , g 2 = 0.01 , g 3 = 0.01 ; (c) g 1 = 0.18 , g 2 = 0.01 , g 3 = 0.01 .
Sustainability 12 09348 g008
Figure 9. Sensitivity of initial values in the system (8) when in stable state: (a) s = 0.7 ,   P 2 = 15 ; (b) P 1 = 11 ,   P 2 = 15 ; (c) P 1 = 11 ,   s = 0.7 .
Figure 9. Sensitivity of initial values in the system (8) when in stable state: (a) s = 0.7 ,   P 2 = 15 ; (b) P 1 = 11 ,   P 2 = 15 ; (c) P 1 = 11 ,   s = 0.7 .
Sustainability 12 09348 g009
Figure 10. Sensitivity of initial values in the system (8) when in chaotic state: (a) s = 0.7 ,   P 2 = 15 ; (b) P 1 = 11 , P 2 = 15 ; (c) P 1 = 11 , s = 0.7 .
Figure 10. Sensitivity of initial values in the system (8) when in chaotic state: (a) s = 0.7 ,   P 2 = 15 ; (b) P 1 = 11 , P 2 = 15 ; (c) P 1 = 11 , s = 0.7 .
Sustainability 12 09348 g010
Figure 11. 3D stable region for g 1 , g 2 , and g 3 under different values of γ : (a) γ = 1 , (b) γ = 3 , (c) γ = 5 .
Figure 11. 3D stable region for g 1 , g 2 , and g 3 under different values of γ : (a) γ = 1 , (b) γ = 3 , (c) γ = 5 .
Sustainability 12 09348 g011
Figure 12. 3D stable region for g 1 , g 2 , and g 3 under different values of σ : (a) σ = 50 , (b) σ = 100 .
Figure 12. 3D stable region for g 1 , g 2 , and g 3 under different values of σ : (a) σ = 50 , (b) σ = 100 .
Sustainability 12 09348 g012
Figure 13. The bifurcation diagrams of players’ profit when g 1 changes: (a,b)   g 2 = 0.01 , g 3 = 0.01 ; (c,d)   g 2 = 0.02 , g 2 = 0.01 .
Figure 13. The bifurcation diagrams of players’ profit when g 1 changes: (a,b)   g 2 = 0.01 , g 3 = 0.01 ; (c,d)   g 2 = 0.02 , g 2 = 0.01 .
Sustainability 12 09348 g013aSustainability 12 09348 g013b
Figure 14. The bifurcation diagrams of players’ profit when g 3 changes: (a,b) g 1 = 0.1 , g 2 = 0.01 ; (c,d)   g 1 = 0.1 ,   g 2 = 0.02 .
Figure 14. The bifurcation diagrams of players’ profit when g 3 changes: (a,b) g 1 = 0.1 , g 2 = 0.01 ; (c,d)   g 1 = 0.1 ,   g 2 = 0.02 .
Sustainability 12 09348 g014aSustainability 12 09348 g014b
Figure 15. The bifurcation diagram of system (14) when L changes: (a)   P 1 and P 2 ; (b)   s .
Figure 15. The bifurcation diagram of system (14) when L changes: (a)   P 1 and P 2 ; (b)   s .
Sustainability 12 09348 g015
Figure 16. Strange attractor of the system (13): (a)   L = 0.1 ; (b)   L = 0.4 .
Figure 16. Strange attractor of the system (13): (a)   L = 0.1 ; (b)   L = 0.4 .
Sustainability 12 09348 g016
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Zhu, C. Research on the Dynamic Pricing and Service Decisions in the Reverse Supply Chain considering Consumers’ Service Sensitivity. Sustainability 2020, 12, 9348. https://doi.org/10.3390/su12229348

AMA Style

Zhang J, Zhu C. Research on the Dynamic Pricing and Service Decisions in the Reverse Supply Chain considering Consumers’ Service Sensitivity. Sustainability. 2020; 12(22):9348. https://doi.org/10.3390/su12229348

Chicago/Turabian Style

Zhang, Jixiang, and Chen Zhu. 2020. "Research on the Dynamic Pricing and Service Decisions in the Reverse Supply Chain considering Consumers’ Service Sensitivity" Sustainability 12, no. 22: 9348. https://doi.org/10.3390/su12229348

APA Style

Zhang, J., & Zhu, C. (2020). Research on the Dynamic Pricing and Service Decisions in the Reverse Supply Chain considering Consumers’ Service Sensitivity. Sustainability, 12(22), 9348. https://doi.org/10.3390/su12229348

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop