# Complexity Dynamic Character Analysis of Retailers Based on the Share of Stochastic Demand and Service

Authors: Junhai Ma; Weiya Di; Hao Ren
Journal: Complexity
 (2017)

---

## Abstract

Apart from the price fluctuation, the retailers’ service level becomes another key factor that affects the market demand. This paper depicts a modified price and demand game model based on the stochastic demand and the retailer’s service level which influences the market demand decided by customers’ preference, while the market demand is stochastic in this model. We explore how the price adjustment speed affects the stability of the supply chain system with respect to service level and stochastic demand. The dynamic behavior of the system is researched by simulation and the stability domain and the bifurcation phenomenon are shown clearly. The largest Lyapunov exponent and the chaotic attractor are also given to confirm the chaotic characteristic of the system. The simulation results indicate that relatively small price adjustment speed may maintain the system at stable state. With the price adjustment speed gradually increasing, the price system gets unstable and finally becomes chaotic. This chaotic phenomenon will perturb the product market and this phenomenon should be controlled to keep the system stay in the stable region. So the chaos control is done and the chaos can be controlled completely. The conclusion makes significant contribution to the system referring to the price fluctuation based on the service level and stochastic demand.


---

## Body

## 1. Introduction

Bertrand duopoly model is frequently used when the competitive behaviors are discussed in proper research. Hence, the dynamic behaviors of the whole supply chain system have already been probed. However, the retail service level plays an important role in influencing the retail price which indirectly changes the dynamic behaviors of the whole supply chain system that contains one supplier and two retailers.The two retailers are duopoly in the product market. Xiao and Yang [1] established a price-service competition model with the circumstance whereby the demand is uncertainty in two competing supply chains considering the risk sensitivity. The outcome showed the higher service investment efficiency of one retailer may lead to lower optimal price and service level of his rival. Yan and Pei [2] focused on the retail service in a dual-channel supply chain, while the result suggested that the improved retail services enhanced the supply chain performance effectively in a dual-channel supply chain. Hall et al. [3] engaged in the dynamic retail price and ordering decisions in the category management setting, and this research found brand’s own price and the cross-price effect had an effect on the brand-by-brand market. Ghosh and Mitra [4] examined the different and opposite results to the ones obtained from a similar comparison in the standard setting through comparing the Bertrand and Cournot model in mixed markets with public and private enterprises considering maximum welfare and maximum profit separately.C. H. Tremblay and V. J. Tremblay [5] introduced a Cournot-Bertrand duopoly model in a Cournot-type firm and a Bertrand-type firm based on the degree of product differentiation. The Cournot-Bertrand equilibrium was calculated and the analysis of the equilibrium was also conducted. Sieke et al. [6] investigated the service level based supply contracts including flat penalty contract and unit penalty contract in a coordinated supply chain. The optimal values of the contract parameters were decided and service level measures were connected to traditional service level measures. Ma and Zhang [7] analyzed the property of three oligarchs which have different rationalities in insurance market based on the price game model and variable feedback control method. Naimzada and Tramontana [8] studied the liner demand and cost functions with product differentiation, and a dynamic model was established to probe into the stability properties in the mixed Cournot-Bertrand duopoly model. Elsadany [9] developed a triopoly game to indicate the characteristic of a dynamic Cournot game through three bounded rational players. Jha and Shanker [10] introduced the circumstance where a supplier produced products and supplied them to a set of buyers. The lead time can be reduced with an added cash cost and a service level constraint was also included in the research. Wang and Ma [11] centered on the limited information in a Cournot-Bertrand mixed duopoly game model. The bounded rational principles were used to help the two firms make management decisions. He discussed the Nash equilibrium point and the local stability of the game. Mahmoodi and Eshghi [12] added the price competition and the stochastic demand into the duopoly supply chains; the effect of price and the uncertainty demand on the system were discussed.Ahmed et al. [13] pointed out that the gradient adjustment mechanism of price was applied in a dynamic Bertrand duopoly in which the goods were differentiated. Kawabata and Takarada [14] probed into the effect of free trade agreements on the welfare in Bertrand and Cournot competition model under differentiated products. Brianzoni et al. [15] studied the mathematical properties and dynamics of a Bertrand duopoly with horizontal differentiated products and nonlinear costs, and the new evidence has been offered. Liu and Wang [16] investigated the quality control game that discussed the different risk attitudes affecting quality control game of supply chain especially in logistics service supply chain, instead of considering one member’s risk attitude and ignoring the combination of two members’ risk attitudes. Fang and Shou [17] managed supply uncertainty under supply chain Cournot competition between two supply chains consisting of a retailer and an exclusive supplier, respectively.Esmaeilzadeh and Taleizadeh [18] conducted the optimal pricing decisions for two complementary products in a two-echelon supply chain under different market powers with game theory approaches. Ma and Guo [19] examined the impacts of information on the dynamical price game in two Bertrand game models, whereby the player obtains information of his rival before making decisions. Q.-H. Li and B. Li [20] designed the dual-channel supply chain; thus channel competition became inevitable. Hence, value-added services provided by retailers were considered. Equilibrium problems regarding retail services and fairness concerns were analyzed. Protopappa-Sieke et al. [21] developed the optimal two-period inventory allocation policies under multiple service level contracts in view of the fact that optimal inventory allocation had a significant impact on profits in the retail industry. How a manufacturer responds to a service level contract was analyzed as considering to minimize the expected costs.The structure of the paper is as follows. Firstly, the basic Bertrand duopoly model considering the service level and stochastic demand is established in Section2. Secondly, the analysis of the model mentioned in the above section has been conducted in Section 3. Then, we simulate the stability domain and the dynamical behaviors, and the related discussions have been made in Section 4. Lastly, the conclusion will be given in Section 5.


## 2. Demand Model

Since oligopoly game model has already been researched in previous literature, while the oligopoly situation has been discussed, we consider that the other retailer comes into the monopoly market. Due to the entrance of the other retailer, exclusive monopoly was broken. What is more, the marketing competition is more vehement. Retailer’s demand is influenced by the other retailer’s price. In addition to retail price, we also consider that the market demand is associated with the service level of themselves and the competitors. Meanwhile, we assume the products that two retailers sell are homogeneous. Hence, the demand function of the two retailers is shown as follows:(1)d1=a1-b1p1+b2p2+k1v1-k2v2,d2=a2-b1p2+b2p1+k1v2-k2v1,where ki (i=1,2) shows the ratio of customers buying the product, which represents the service discrepancy level of the products. bi (i=1,2) represents the price sensitive coefficient, which expresses the degree of market demand influenced by the price of their products and competition’s product price. Due to the difference of the price and service level between two products, customers will make their own choice according to the price and service level. We must be clear that the higher the service level of the products is, the more the consumers will choose to buy the products in condition of the same price level facing the homogeneous product.In the expression above, we find that the market demand is affected by both retail price and the service level. The higher service level will bring about the larger cost. The cost of service level is shown in(2)ci=ηivi22,i=1,2.In (2), vi means the service level of two retailers. ηi is the influence coefficient, which is used to signify the influence degree of the service level to cost.While the product market is complex and stochastic, the final demand of the customers is uncertain. We set the demand as stochastic variableai, whereby the stochastic demand is the sum of the fixed demand ai,0 and the potential stochastic part with market uncertainty e, that is, ai=ai,0+e. And the expectation and variance of e are 0 and r, respectively. We assume that retailers can predict the stochastic demand of their products. Retailers can help to make decisions with the market information forecast, when facing the uncertain factors in the product market. The forecast value of the stochastic demand is f, and f is composed with stochastic variable ai and the error term ε, that is, fi=ai+ε. The expectation and variance of ε are 0 and s, separately. We have to point out that e and s are independent. Just like the former research that in Li [22], we know(3)Eaif=siri+siai,0+riri+sifi,Ef-ai,02=r+s.The profit function of the manufacturer is depicted in (4). In this equation, w represents the wholesale prices that manufacturer sells the products to retailer:(4)πm=d1+d2w.The retailers’ operating goal is to maximize profits. The retailer’s profit function is shown as follows:(5)πR,i=pi-w-cidi,i=1,2.When retailers and manufacturer in the market share the information, the expected profits decision models of manufacturer and retailers are shown as follows, respectively:(6)max⁡Eπmis=d1+d2w,max⁡EπR,1is=p1-w-c1d1f1,max⁡EπR,2is=p2-w-c2d2f2.We can deduce the marginal profits of two retailers through calculating the first-order partial derivative of (6):(7)∂EπR,1NI∂p1=1-g1a0,1+g1f1-2b1p1+b2p2+k1v1-k2v2+b1w+b1η1v122,∂EπR,2NI∂p2=1-g2a0,2+g2f2-2b1p2+b2p1+k1v2-k2v1+b1w+b1η2v222.For the sake of simplicity, we makegi=ri/ri+si  (i=1,2) in (7).


### 2.1. Optimal Decision

If there is only one retailer in the market, the retailer regards profit maximization as the goal to give the optimal decision in the circumstance of discrete game. When there are two retailers in the market, the optimal decision of retailer 1 is influenced by retailer 2. Meanwhile, retailer 2 optimal decision is also affected by retailer 1 decision. Hence, we assume the information in the market can be shared between two retailers; that is, the optimal decisions of two retailers are known to each other. When the marginal profits reduce to zero, we have the optimal decision for retailers. The retail price and the service level are the decision variables.For retailer 1, we calculate the first-order derivative of its profit. When the marginal profit gets to zero, the optimal retail price and optimal service level will be derived in equation. The detailed process is as follows:(8)∂EπR,1NI∂p1=0,∂EπR,1NI∂V1=0,(9)p1is=12b11-g1a0,1+g1f1+b2p2-k2v2+b1w+3k122η1b1,v1is=k1b1η1.Analogously, the optimal decision of retailer 2 will be conducted. Equation (10) has the same meaning as (8). The optimal decision of retailer 2 is expressed in (11):(10)∂EπR,2NI∂p2=0,∂EπR,2NI∂V2=0,(11)p2is=12b11-g2a0,2+g2f2+b2p1-k2v1+b1w+3k122η2b1,v2is=k1b1η2.When two retailers make their own optimal decision, we may further find the second derivatives of the retail prices and service level are less than zero. Therefore, the optimal decision of the system will be derived, which means the system has the equilibrium point:p1∗,p2∗,v1∗,v2∗. The calculation will be expressed as follows:(12)∂2EπR,1NI∂2p1<0,∂2EπR,1NI∂2V1<0,∂2EπR,2NI∂2p2<0,∂2EπR,2NI∂2V2<0.Through jointing (9) and (11), the optimal solution of the system will be solved:(13)p1∗=2b14b12-b22R1+b22b1R2+b22b1B+A,p2∗=2b14b12-b22R2+b22b1R1+b22b1A+B,v1∗=k1b1η1,v2∗=k1b1η2.In (13), we have R1=(1-g1)a0,1+g1f1, R2=(1-g2)a0,2+g2f2, A=b1w+k1(3k1η2-k2η1)/b1η1η2, B=b1w+k1(3k1η1-k2η2)/b1η1η2.


### 2.2. The Complex Dynamics of the System

Due to the complexity of the product market, the competition between enterprises in the market is increasingly fierce. Retailers do not know the perfect information of the market, while the market is complex and customers’ demands are uncertain. In addition, the more the transactions are, the greater the uncertainty will be, the more incomplete the information will be. Hence, the market demand cannot be forecasted accurately. The bounded rationality expectations solve this problem faultlessly. The retail price of next period will be decided according to the retail price in this period and the marginal profits, which makes up the discrete dynamic system of the product market.(14)p1t+1=p1t+θ1p1t1-g1a0,1+g1f1-2b1p1+b2p2+k1v1-k2v2+b1w+b1η1v122,p2t+1=p2t+θ2p2t1-g2a0,2+g2f2-2b1p2+b2p1+k1v2-k2v1+b1w+b1η2v222.In the equation above,θi  (i=1,2) is the price adjustment speed of retailer i, and the price adjustment speed is positive, which is set to make sure that retail price of next period is practically significant. When the marginal profit takes positive value, the retail price of next period will increase. However, when the marginal profit takes negative value, the retail price of next period will decrease.The system (14) is nonlinear; the system is running under constant interference. When the system is disturbed, whether it can maintain a predetermined trajectory, the stability of the system becomes the primary question to study. What is different from the linear system is that the nonlinear system has multiple equilibriums. We obtain four equilibriums of system (14): E0(0,0), E1(0,R2+B/2b1), E2(R1+A/2b1,0), E∗((2b1/4b12-b22)(R1+b2/2b1R2+b2/2b1B+A),(2b1/(4b12-b22))(R2+(b2/2b1)R1+b2/2b1A+B)).We may find thatE0, E1, E2 are boundary points, and E∗ is within the system. We guess that E∗ is the unique Nash equilibrium. In order to verify the guess and analyze the stability of the system, we calculate the Jacobian matrix of system (14):(15)JE=J11J12J21J22.In the equation above,(16)J11=1+θ1R1+θ1b2p2+θ1A-4θ1b1p1,J12=θ1b2p1,J21=θ2b2p2,J22=1+θ2R2+θ1b2p1+θ2B-4θ2b1p.E
∗ is the unique Nash equilibrium; the best profit of the system can be obtained when the best decision is in E∗. We all know that when the margin profit is more than zero, the retailers will adjust their retail prices, until the margin profit goes up to zero. At this time, the system reaches to equilibrium state. Retailers will not change their own retail price any more, and the equilibrium state will be maintained by these two game players. Therefore, among all the equilibrium points, only the Nash equilibrium E∗ has practical significance. Then, this part will explore the local stability of the Nash equilibrium by using the Jacobian matrix and the July criterion.The characteristic polynomial of the Jacobian matrix takes the following form:(17)Fλ=λ2-J11+J22λ-J11J22-J12J21.In order to guarantee that the Nash equilibrium is local stability,Ji,j  (i,j=1,2) must satisfy the following limitations:(18)-J11+J22+J11J22-J12J21+1>0,J11+J22+J11J22-J12J21+1>0,J11J22-J12J21<1,J11+J22>0.Due to the fact that limitations are so complex, the process to solve inequality (18) is very complicated. If the Nash equilibrium satisfies inequality (18), we may insure that Nash equilibrium is locally stable. We will probe into the dynamic characteristic of the system through numerical simulation in next section.


## 2.1. Optimal Decision

If there is only one retailer in the market, the retailer regards profit maximization as the goal to give the optimal decision in the circumstance of discrete game. When there are two retailers in the market, the optimal decision of retailer 1 is influenced by retailer 2. Meanwhile, retailer 2 optimal decision is also affected by retailer 1 decision. Hence, we assume the information in the market can be shared between two retailers; that is, the optimal decisions of two retailers are known to each other. When the marginal profits reduce to zero, we have the optimal decision for retailers. The retail price and the service level are the decision variables.For retailer 1, we calculate the first-order derivative of its profit. When the marginal profit gets to zero, the optimal retail price and optimal service level will be derived in equation. The detailed process is as follows:(8)∂EπR,1NI∂p1=0,∂EπR,1NI∂V1=0,(9)p1is=12b11-g1a0,1+g1f1+b2p2-k2v2+b1w+3k122η1b1,v1is=k1b1η1.Analogously, the optimal decision of retailer 2 will be conducted. Equation (10) has the same meaning as (8). The optimal decision of retailer 2 is expressed in (11):(10)∂EπR,2NI∂p2=0,∂EπR,2NI∂V2=0,(11)p2is=12b11-g2a0,2+g2f2+b2p1-k2v1+b1w+3k122η2b1,v2is=k1b1η2.When two retailers make their own optimal decision, we may further find the second derivatives of the retail prices and service level are less than zero. Therefore, the optimal decision of the system will be derived, which means the system has the equilibrium point:p1∗,p2∗,v1∗,v2∗. The calculation will be expressed as follows:(12)∂2EπR,1NI∂2p1<0,∂2EπR,1NI∂2V1<0,∂2EπR,2NI∂2p2<0,∂2EπR,2NI∂2V2<0.Through jointing (9) and (11), the optimal solution of the system will be solved:(13)p1∗=2b14b12-b22R1+b22b1R2+b22b1B+A,p2∗=2b14b12-b22R2+b22b1R1+b22b1A+B,v1∗=k1b1η1,v2∗=k1b1η2.In (13), we have R1=(1-g1)a0,1+g1f1, R2=(1-g2)a0,2+g2f2, A=b1w+k1(3k1η2-k2η1)/b1η1η2, B=b1w+k1(3k1η1-k2η2)/b1η1η2.


## 2.2. The Complex Dynamics of the System

Due to the complexity of the product market, the competition between enterprises in the market is increasingly fierce. Retailers do not know the perfect information of the market, while the market is complex and customers’ demands are uncertain. In addition, the more the transactions are, the greater the uncertainty will be, the more incomplete the information will be. Hence, the market demand cannot be forecasted accurately. The bounded rationality expectations solve this problem faultlessly. The retail price of next period will be decided according to the retail price in this period and the marginal profits, which makes up the discrete dynamic system of the product market.(14)p1t+1=p1t+θ1p1t1-g1a0,1+g1f1-2b1p1+b2p2+k1v1-k2v2+b1w+b1η1v122,p2t+1=p2t+θ2p2t1-g2a0,2+g2f2-2b1p2+b2p1+k1v2-k2v1+b1w+b1η2v222.In the equation above,θi  (i=1,2) is the price adjustment speed of retailer i, and the price adjustment speed is positive, which is set to make sure that retail price of next period is practically significant. When the marginal profit takes positive value, the retail price of next period will increase. However, when the marginal profit takes negative value, the retail price of next period will decrease.The system (14) is nonlinear; the system is running under constant interference. When the system is disturbed, whether it can maintain a predetermined trajectory, the stability of the system becomes the primary question to study. What is different from the linear system is that the nonlinear system has multiple equilibriums. We obtain four equilibriums of system (14): E0(0,0), E1(0,R2+B/2b1), E2(R1+A/2b1,0), E∗((2b1/4b12-b22)(R1+b2/2b1R2+b2/2b1B+A),(2b1/(4b12-b22))(R2+(b2/2b1)R1+b2/2b1A+B)).We may find thatE0, E1, E2 are boundary points, and E∗ is within the system. We guess that E∗ is the unique Nash equilibrium. In order to verify the guess and analyze the stability of the system, we calculate the Jacobian matrix of system (14):(15)JE=J11J12J21J22.In the equation above,(16)J11=1+θ1R1+θ1b2p2+θ1A-4θ1b1p1,J12=θ1b2p1,J21=θ2b2p2,J22=1+θ2R2+θ1b2p1+θ2B-4θ2b1p.E
∗ is the unique Nash equilibrium; the best profit of the system can be obtained when the best decision is in E∗. We all know that when the margin profit is more than zero, the retailers will adjust their retail prices, until the margin profit goes up to zero. At this time, the system reaches to equilibrium state. Retailers will not change their own retail price any more, and the equilibrium state will be maintained by these two game players. Therefore, among all the equilibrium points, only the Nash equilibrium E∗ has practical significance. Then, this part will explore the local stability of the Nash equilibrium by using the Jacobian matrix and the July criterion.The characteristic polynomial of the Jacobian matrix takes the following form:(17)Fλ=λ2-J11+J22λ-J11J22-J12J21.In order to guarantee that the Nash equilibrium is local stability,Ji,j  (i,j=1,2) must satisfy the following limitations:(18)-J11+J22+J11J22-J12J21+1>0,J11+J22+J11J22-J12J21+1>0,J11J22-J12J21<1,J11+J22>0.Due to the fact that limitations are so complex, the process to solve inequality (18) is very complicated. If the Nash equilibrium satisfies inequality (18), we may insure that Nash equilibrium is locally stable. We will probe into the dynamic characteristic of the system through numerical simulation in next section.


## 3. Dynamic Characteristics of the System

In consideration of the complexity of system (14), the numerical simulation method will be employed in this paper as this method has been widely used in most of the literatures. The stable region, bifurcation, and chaos will be investigated; meanwhile, the largest Lyapunov exponent, entropy, chaotic attractor, and the time domain response will be discussed to verify the dynamic characteristics of this system. At last, the sensitivity of initial value will be analyzed. We take the parameter values as follows: g1=0.7, g2=0.6, a0,1=0.5, a0,2=0.55, f1=0.6, f2=0.65, b1=1, b2=0.5, k1=1.1, k2=1.2, v1=0.7, v2=0.6, w=0.1, η1=0.2, η2=0.3.


### 3.1. The Stable Region

The decision variables of the system are retail prices of two retailers. As retail prices obtain different values, the system will change to either the stable range or the unstable range. Figure1 depicts the stable region of the system with blue range and the unstable of the system with white region. When employed values are in the blue range, the system will tend to the equilibrium state after iteration period by period. Then the system will stay at the Nash equilibrium state unless some factors outside break the equilibrium state.Figure 1
The stable range of system (14).


### 3.2. The Bifurcation and Chaos Behavior

The fixed point is stable when all the eigenvalues of the Jacobian are in the unit circle. Whenθ1 and θ2 satisfy the critical value, the system appears as saddle point. If one characteristic root λ2=-1 and all the other eigenvalues are within the unit circle, the fixed point will be unstable and produce period-doubling bifurcation. From Figure 2, we can see that system (14) will appear as period-doubling bifurcation.Figure 2
The characteristic root of system (14).We discuss the bifurcation and chaos behaviors of the system in this section with parameter basin plots and the bifurcation diagram. We set the initial values of the system asp1(0)=0.25 and p2(0)=0.2, which represent the retail price of two retailers at the first moment.Figure3 describes the parameter basin plots of system (14) with respect to the price adjustment speeds θ1 and θ2. Different colors in this figure represent the stable region, stable cycles of periods, and chaotic states separately: the light blue range is the stable region, pink range means the stable cycles of period 2, orange range is the stable cycles of period 3, yellow range is the stable cycles of period 4, green range is the stable cycles of period 5, dark blue range is the stable cycles of period 6, cyan range is the stable cycles of period 7, purple range is the stable cycles of period 8, red range represents chaos, and white means divergence. If the retail prices are in the white region, one of the retailers is outside the product market, and the market will be oligopolistic market. Through the observation of Figure 3, we find pink, yellow, dark blue, and purple constitute the road from even period cycles to chaos, the road being called the flip bifurcation in discrete dynamical system; meanwhile, orange, green, and cyan constitute the road from odd period cycles to chaos.Figure 3
The parameter basin plots of system (14).By observing this figure, we may find the retail price changes from the stable region to the even period cycles and next it gets to chaotic region. This phenomenon is the flip bifurcation, whereby the retailers adjust their own retail price, the system will go from the even period cycles to the chaotic region directly.Figures4(a) and 4(b) give the bifurcation diagram of retail prices with respect to θ1 when θ2=1.2. This figure can explain the complex dynamic behavior of the system more intuitively, whereby one of the price adjustment speeds is fixed and the other adjustment speed changes within a certain range. This paper investigates the adjustment speed of retailers’ prices firstly by keeping θ2 fixed at 1.2 and θ1varying from zero to 3. In Figure 4(a), retailer 1’s retail price is stable at 0.49 as θ1 vary from zero to 1.95. Just at this time, the system begins the first bifurcation, while the system runs into stable cycles of period 2. As θ1 increases to 2.47, the system gets to stable cycles of period 4. When θ1 gets the value of 2.58, retail 1’s price goes into stable cycles of period 8. After that, the system falls into a chaotic state.Figure 4
(a)
Bifurcation diagram ofp1 with respect to θ1 when θ2=1.2
(b)
Bifurcation diagram ofp2 with respect to θ1 when θ2=1.2Figure4(b) depicts the bifurcation diagram of p2 with respect to θ1 when θ2=1.2. The bifurcation points of this figure are completely the same as the bifurcation points in Figure 4(a). However, what is different from Figure 4(a) is the value of retail price. Retailer 2’s retail price is stable at 0.415 as θ1 vary from zero to 1.95. Next, p2 gets to cycles of periods 2, 4, and 8 in sequence. In the end, p2 also falls into chaos.Figures5(a) and 5(b) describe the bifurcation behavior of retail prices when θ2 increases gradually and θ1 keeps the fixed value 1.2, which is the same as Figure 4; this shows the complex dynamic behavior of retail price considering the price adjustment speed of retailer 2 varying from zero to 3.5. In Figure 5(a), as θ2 is small and less than 2.22; the retail price of retailer 1 keeps a fixed number of 0.488. After that, the retail price falls into the stable cycles of period 2; when θ2increases to 2.85, the retail price falls into the stable cycles of period 4; then when θ2 goes to 2.98, the retail price falls into the stable cycles of period 8; and so forth. After the system passes through multiple iterations, it will turn into chaos, which will not be expected to appear.Figure 5
(a)
Bifurcation diagram ofp1 with respect to θ2 when θ1=1.2
(b)
Bifurcation diagram ofp2 with respect to θ2 when θ1=1.2Analogously, the complex dynamic behavior of retail price 2 considering the price adjustment speed of retailer 2 varying from zero to 3.5 will be shown in Figure5(b). Firstly, retailer price 2 is stable at 0.414 when θ2 is at [0, 2.22]. Then the system turns into bifurcation and chaos state orderly. It is worth noting that the stable price of retailer 2 is less than that of retailer 1, which means retailer 2 may possess more market share. This phenomenon is not conducive to retailer 1. And retailer 1 must take measures to decrease its own retail price to win more market demand.Through the analysis of Figures4 and 5, we find that relatively small price adjustment speed may maintain the system at stable state. Figure 6 shows the stability of the system from a three-dimensional perspective. With the price adjustment speed gradually increasing, the price system gets unstable and finally becomes chaotic. Hence, from the point of view of management, the retailer should choose a relatively reasonable adjustment speed to alter the retail price in each decision unit.Figure 6
3D bifurcation diagram ofp1 and p2.


### 3.3. The Largest Lyapunov Exponent

The largest Lyapunov exponent (LLE) is another way to make clear the complex dynamic characteristics of the system. And the LLE represents the extreme sensitivity to initial value of the system vividly that two of the same initial value generated in the trajectory separate according to index methods with the passage of time. As we want to judge whether the chaotic phenomenon exists in the system, we just need to observe the LLE intuitively. If the LLE is more than zero, it means that the differences of the initial value increase with exponential form after periods of iteration; that is, the system falls into chaotic region. On the contrary, if the LLE is less than zero, it means that the system is in stable state, while the system is not initially sensitive and keeps a fixed value all the time. When the LLE equals zero, the system is periodic and quasi-periodic. The expression of LLE is as follows:(19)λ=limn→∞⁡1n∑n=0n-1ln⁡dfxn,μdx.We used MATLAB to simulate Figures7(a) and 7(b). Figures 7(a) and 7(b) explore the LLE with the increase of price adjustment speed. The system acts out the complex dynamic behavior as the price adjustment speed changes. In Figure 7(a), when θ1 is in [0, 1.95], the LLE is less than zero and the system is in stable state; while θ1 is in [1.95, 2.6], the LLE is equal to zero and the system is in bifurcation region; when θ1 is more than 2.6, the LLE is more than zero and the system is in chaotic region. Figure 7(b) shows the LLE with respect to θ2, and it has the same changing process with θ1. Figure 7 accounts for the state of the dynamic system more intuitively; meanwhile, the bifurcation point and the chaos point in Figure 7 are in consistent with the bifurcation diagrams.Figure 7
(a)
The largest Lyapunov exponent of system with respect toθ1 when θ2=1.2
(b)
The largest Lyapunov exponent of system with respect toθ2 when θ1=1.2


### 3.4. The Entropy

Entropy is an important quantity to characterize chaotic systems. In different types of dynamic systems, the entropy values are different. Entropy values can be used to distinguish between regular motion, chaotic motion, and random motion. In the random motion system, the entropy is unbounded; in the regular motion system, the entropy is zero; in the chaotic motion system, the entropy is greater than zero. The greater the entropy, the greater the rate of loss of information, the greater the degree of chaos in the system, or the more complex the system.It can be seen from Figure8 that as the system progressively enters chaos, the entropy of the system increases rapidly from zero. The law of change is also the same as the change of bifurcation and maximum Lyapunov exponent.Figure 8
(a)
The entropy of system with respect toθ1 when θ2=1.2
(b)
The entropy of system with respect toθ2 when θ1=1.2


### 3.5. The Chaotic Attractor

The chaotic attractor investigates the chaotic characteristics of the system and the chaotic attractors will be given in Figures9(a), 9(b), and 9(c), while the price adjustment speeds are in stable region and chaotic region, respectively. Strange attractor is the joint action’s product of overall stability and local instability. All the movement outside the attractor will finally converge to the chaotic attractor. Meanwhile the movement inside the attractor is instable. If the attractor is a fixed point, the system is in stable state. On the contrary, when the attractor is not a fixed point and period cycle, the system is in chaos.Figure 9
(a)
The chaotic attractor of system whenθ1=θ2=1.2
(b)
The chaotic attractor of system whenθ1=2.8, θ2=1.2
(c)
The chaotic attractor of system whenθ1=1.2, θ2=3.2In Figure9(a), we judged on the parameter basin plots, bifurcation diagram, the largest Lyapunov exponent, and entropy of system; when θ1 and θ2 equal 1.2, they are all in stable states, and at this moment the attractor is a fixed point.From the parameter basin plots, bifurcation diagram, the largest Lyapunov exponent, and entropy of system, we may see that the systems are in chaotic states whileθ1=2.8, θ2=1.2 and θ1=1.2, θ2=3.2. Hence, the attractors expand to a big area, and the system begins to do limit movement in Figures 9(b) and 9(c). Then the attractors turn into chaotic attractors. Just now, the system becomes complex and the dynamic characteristics of the system are more obviously.


### 3.6. Time Domain Response

The time domain response depicts the changing process of the price data with the increase of time, which is arranged by the whole time index. There are different shapes for the time domain response. Figures10(a), 10(b), and 10(c) describe the time domain responses of the system when the price adjustment speeds are in stable region and chaotic region, respectively. In Figure 10, blue line represents the retail price of retailer 1, and red line represents the retailer price of retailer 2.Figure 10
(a)
The time domain response of system whenθ1=θ2=1.2
(b)
The time domain response of system whenθ1=2.8, θ2=1.2
(c)
The time domain response of system whenθ1=1.2, θ2=3.2In Figure10(a), θ1 and θ2 equal 1.2. The price adjustment speeds are all in stable states. Then the two retail price fix at the equilibrium point E∗(0.488,0.414), and two retail price will not change as the time changes. It is the same as Figure 9; we set θ1=2.8, θ2=1.2, and θ1=1.2, θ2=3.2 in Figures 10(b) and 10(c). In Figure 10(b), while θ1 is in stable region and θ2 is in chaotic area, the fluctuation of blue line is obviously big compared to that of the red line, and the red line just tends to the value 0.36. Uniformly in Figure 10(c), the fluctuation of blue line is obviously small compared to that of the red line, and the blue line just tends to the value 0.45 as θ2 is in stable region and θ1 is in chaotic region. This phenomenon shows that when the price adjustment speed is in stable state, the time domain response is smooth. On the contrary, the time domain response will be in big fluctuation while the price adjustment speed is in chaotic state.


## 3.1. The Stable Region

The decision variables of the system are retail prices of two retailers. As retail prices obtain different values, the system will change to either the stable range or the unstable range. Figure1 depicts the stable region of the system with blue range and the unstable of the system with white region. When employed values are in the blue range, the system will tend to the equilibrium state after iteration period by period. Then the system will stay at the Nash equilibrium state unless some factors outside break the equilibrium state.Figure 1
The stable range of system (14).


## 3.2. The Bifurcation and Chaos Behavior

The fixed point is stable when all the eigenvalues of the Jacobian are in the unit circle. Whenθ1 and θ2 satisfy the critical value, the system appears as saddle point. If one characteristic root λ2=-1 and all the other eigenvalues are within the unit circle, the fixed point will be unstable and produce period-doubling bifurcation. From Figure 2, we can see that system (14) will appear as period-doubling bifurcation.Figure 2
The characteristic root of system (14).We discuss the bifurcation and chaos behaviors of the system in this section with parameter basin plots and the bifurcation diagram. We set the initial values of the system asp1(0)=0.25 and p2(0)=0.2, which represent the retail price of two retailers at the first moment.Figure3 describes the parameter basin plots of system (14) with respect to the price adjustment speeds θ1 and θ2. Different colors in this figure represent the stable region, stable cycles of periods, and chaotic states separately: the light blue range is the stable region, pink range means the stable cycles of period 2, orange range is the stable cycles of period 3, yellow range is the stable cycles of period 4, green range is the stable cycles of period 5, dark blue range is the stable cycles of period 6, cyan range is the stable cycles of period 7, purple range is the stable cycles of period 8, red range represents chaos, and white means divergence. If the retail prices are in the white region, one of the retailers is outside the product market, and the market will be oligopolistic market. Through the observation of Figure 3, we find pink, yellow, dark blue, and purple constitute the road from even period cycles to chaos, the road being called the flip bifurcation in discrete dynamical system; meanwhile, orange, green, and cyan constitute the road from odd period cycles to chaos.Figure 3
The parameter basin plots of system (14).By observing this figure, we may find the retail price changes from the stable region to the even period cycles and next it gets to chaotic region. This phenomenon is the flip bifurcation, whereby the retailers adjust their own retail price, the system will go from the even period cycles to the chaotic region directly.Figures4(a) and 4(b) give the bifurcation diagram of retail prices with respect to θ1 when θ2=1.2. This figure can explain the complex dynamic behavior of the system more intuitively, whereby one of the price adjustment speeds is fixed and the other adjustment speed changes within a certain range. This paper investigates the adjustment speed of retailers’ prices firstly by keeping θ2 fixed at 1.2 and θ1varying from zero to 3. In Figure 4(a), retailer 1’s retail price is stable at 0.49 as θ1 vary from zero to 1.95. Just at this time, the system begins the first bifurcation, while the system runs into stable cycles of period 2. As θ1 increases to 2.47, the system gets to stable cycles of period 4. When θ1 gets the value of 2.58, retail 1’s price goes into stable cycles of period 8. After that, the system falls into a chaotic state.Figure 4
(a)
Bifurcation diagram ofp1 with respect to θ1 when θ2=1.2
(b)
Bifurcation diagram ofp2 with respect to θ1 when θ2=1.2Figure4(b) depicts the bifurcation diagram of p2 with respect to θ1 when θ2=1.2. The bifurcation points of this figure are completely the same as the bifurcation points in Figure 4(a). However, what is different from Figure 4(a) is the value of retail price. Retailer 2’s retail price is stable at 0.415 as θ1 vary from zero to 1.95. Next, p2 gets to cycles of periods 2, 4, and 8 in sequence. In the end, p2 also falls into chaos.Figures5(a) and 5(b) describe the bifurcation behavior of retail prices when θ2 increases gradually and θ1 keeps the fixed value 1.2, which is the same as Figure 4; this shows the complex dynamic behavior of retail price considering the price adjustment speed of retailer 2 varying from zero to 3.5. In Figure 5(a), as θ2 is small and less than 2.22; the retail price of retailer 1 keeps a fixed number of 0.488. After that, the retail price falls into the stable cycles of period 2; when θ2increases to 2.85, the retail price falls into the stable cycles of period 4; then when θ2 goes to 2.98, the retail price falls into the stable cycles of period 8; and so forth. After the system passes through multiple iterations, it will turn into chaos, which will not be expected to appear.Figure 5
(a)
Bifurcation diagram ofp1 with respect to θ2 when θ1=1.2
(b)
Bifurcation diagram ofp2 with respect to θ2 when θ1=1.2Analogously, the complex dynamic behavior of retail price 2 considering the price adjustment speed of retailer 2 varying from zero to 3.5 will be shown in Figure5(b). Firstly, retailer price 2 is stable at 0.414 when θ2 is at [0, 2.22]. Then the system turns into bifurcation and chaos state orderly. It is worth noting that the stable price of retailer 2 is less than that of retailer 1, which means retailer 2 may possess more market share. This phenomenon is not conducive to retailer 1. And retailer 1 must take measures to decrease its own retail price to win more market demand.Through the analysis of Figures4 and 5, we find that relatively small price adjustment speed may maintain the system at stable state. Figure 6 shows the stability of the system from a three-dimensional perspective. With the price adjustment speed gradually increasing, the price system gets unstable and finally becomes chaotic. Hence, from the point of view of management, the retailer should choose a relatively reasonable adjustment speed to alter the retail price in each decision unit.Figure 6
3D bifurcation diagram ofp1 and p2.


## 3.3. The Largest Lyapunov Exponent

The largest Lyapunov exponent (LLE) is another way to make clear the complex dynamic characteristics of the system. And the LLE represents the extreme sensitivity to initial value of the system vividly that two of the same initial value generated in the trajectory separate according to index methods with the passage of time. As we want to judge whether the chaotic phenomenon exists in the system, we just need to observe the LLE intuitively. If the LLE is more than zero, it means that the differences of the initial value increase with exponential form after periods of iteration; that is, the system falls into chaotic region. On the contrary, if the LLE is less than zero, it means that the system is in stable state, while the system is not initially sensitive and keeps a fixed value all the time. When the LLE equals zero, the system is periodic and quasi-periodic. The expression of LLE is as follows:(19)λ=limn→∞⁡1n∑n=0n-1ln⁡dfxn,μdx.We used MATLAB to simulate Figures7(a) and 7(b). Figures 7(a) and 7(b) explore the LLE with the increase of price adjustment speed. The system acts out the complex dynamic behavior as the price adjustment speed changes. In Figure 7(a), when θ1 is in [0, 1.95], the LLE is less than zero and the system is in stable state; while θ1 is in [1.95, 2.6], the LLE is equal to zero and the system is in bifurcation region; when θ1 is more than 2.6, the LLE is more than zero and the system is in chaotic region. Figure 7(b) shows the LLE with respect to θ2, and it has the same changing process with θ1. Figure 7 accounts for the state of the dynamic system more intuitively; meanwhile, the bifurcation point and the chaos point in Figure 7 are in consistent with the bifurcation diagrams.Figure 7
(a)
The largest Lyapunov exponent of system with respect toθ1 when θ2=1.2
(b)
The largest Lyapunov exponent of system with respect toθ2 when θ1=1.2


## 3.4. The Entropy

Entropy is an important quantity to characterize chaotic systems. In different types of dynamic systems, the entropy values are different. Entropy values can be used to distinguish between regular motion, chaotic motion, and random motion. In the random motion system, the entropy is unbounded; in the regular motion system, the entropy is zero; in the chaotic motion system, the entropy is greater than zero. The greater the entropy, the greater the rate of loss of information, the greater the degree of chaos in the system, or the more complex the system.It can be seen from Figure8 that as the system progressively enters chaos, the entropy of the system increases rapidly from zero. The law of change is also the same as the change of bifurcation and maximum Lyapunov exponent.Figure 8
(a)
The entropy of system with respect toθ1 when θ2=1.2
(b)
The entropy of system with respect toθ2 when θ1=1.2


## 3.5. The Chaotic Attractor

The chaotic attractor investigates the chaotic characteristics of the system and the chaotic attractors will be given in Figures9(a), 9(b), and 9(c), while the price adjustment speeds are in stable region and chaotic region, respectively. Strange attractor is the joint action’s product of overall stability and local instability. All the movement outside the attractor will finally converge to the chaotic attractor. Meanwhile the movement inside the attractor is instable. If the attractor is a fixed point, the system is in stable state. On the contrary, when the attractor is not a fixed point and period cycle, the system is in chaos.Figure 9
(a)
The chaotic attractor of system whenθ1=θ2=1.2
(b)
The chaotic attractor of system whenθ1=2.8, θ2=1.2
(c)
The chaotic attractor of system whenθ1=1.2, θ2=3.2In Figure9(a), we judged on the parameter basin plots, bifurcation diagram, the largest Lyapunov exponent, and entropy of system; when θ1 and θ2 equal 1.2, they are all in stable states, and at this moment the attractor is a fixed point.From the parameter basin plots, bifurcation diagram, the largest Lyapunov exponent, and entropy of system, we may see that the systems are in chaotic states whileθ1=2.8, θ2=1.2 and θ1=1.2, θ2=3.2. Hence, the attractors expand to a big area, and the system begins to do limit movement in Figures 9(b) and 9(c). Then the attractors turn into chaotic attractors. Just now, the system becomes complex and the dynamic characteristics of the system are more obviously.


## 3.6. Time Domain Response

The time domain response depicts the changing process of the price data with the increase of time, which is arranged by the whole time index. There are different shapes for the time domain response. Figures10(a), 10(b), and 10(c) describe the time domain responses of the system when the price adjustment speeds are in stable region and chaotic region, respectively. In Figure 10, blue line represents the retail price of retailer 1, and red line represents the retailer price of retailer 2.Figure 10
(a)
The time domain response of system whenθ1=θ2=1.2
(b)
The time domain response of system whenθ1=2.8, θ2=1.2
(c)
The time domain response of system whenθ1=1.2, θ2=3.2In Figure10(a), θ1 and θ2 equal 1.2. The price adjustment speeds are all in stable states. Then the two retail price fix at the equilibrium point E∗(0.488,0.414), and two retail price will not change as the time changes. It is the same as Figure 9; we set θ1=2.8, θ2=1.2, and θ1=1.2, θ2=3.2 in Figures 10(b) and 10(c). In Figure 10(b), while θ1 is in stable region and θ2 is in chaotic area, the fluctuation of blue line is obviously big compared to that of the red line, and the red line just tends to the value 0.36. Uniformly in Figure 10(c), the fluctuation of blue line is obviously small compared to that of the red line, and the blue line just tends to the value 0.45 as θ2 is in stable region and θ1 is in chaotic region. This phenomenon shows that when the price adjustment speed is in stable state, the time domain response is smooth. On the contrary, the time domain response will be in big fluctuation while the price adjustment speed is in chaotic state.


## 4. Chaos Control

With the findings of sections above, we come to the conclusion that if the adjustment speeds of the retail prices are in chaotic region, the system will be complex, and the manager cannot make decision in time that phenomenon is not expected to appear in the product market. Therefore, we should take measure to expand the stable region and control the chaos further.From the management point of view, the market is sharply fluctuating while the adjustment speeds of the retail prices expand from the stable region to chaotic region; hence the adjustment speeds of the retail prices must be in stable region so as to keep the market smooth. Then the chaos control method becomes necessary.Here, we employ the widely used control method, that is, parameter adjustment control method to control the chaotic phenomenon. Ma and Li [23] adopted this kind of method to control the complex dynamic behavior of the pricing game models in a risk-averse supply chain, and the chaos has been controlled perfectly. Hence, we also adopt parameter adjustment control method so as to solve this problem in this section. The manager may control the chaotic behavior through introducing the control parameter u, and we set u=0.6 as Ma et al. [24] have conducted the chaos control with setting the control parameter as 0.6. The other parameters employ the values set in the sections above. After that, the stable range and the bifurcation diagram will be given to account for the state of the system. With the parameter adjustment control, the system will be(20)p1t+1=1-up1t+θ1p1t1-g1a0,1+g1f1-2b1p1+b2p2+k1v1-k2v2+b1w+b1η1v122+up1t,p2t+1=1-up2t+θ2p2t1-g2a0,2+g2f2-2b1p2+b2p1+k1v2-k2v1+b1w+b1η2v222+up2t.Figure11 depicts the stable region after having been controlled by the control parameter. Compared to the stable region of system (14), the stable region of system (20) increases obviously, which means a relative large price adjustment speed will not make the system run into chaos. This is helpful for managers to adjust the retail price with a more range.Figure 11
The stable range of system (20).The complex dynamic behavior of system (20) has been investigated in Figures 12 and 13. The first bifurcation point of the system is big compared to that before controlling, the value of which is 5, and it is bigger than value (1.95) of first bifurcation point in system (14) in Figure 9 depicting the bifurcation diagram of p with respect to θ1 when θ2=1.2. Figure 10 shows the bifurcation diagram of p with respect to θ2 when θ1=1.2. It is similar to the phenomenon in Figure 12 that the first bifurcation point of the system is big compared to that before controlling. The first bifurcation point expands to 5.92 that the stable region of system (20) is enlarged. Hence, the managers will make decisions relatively simple.Figure 12
Bifurcation diagram ofp with respect to θ1 when θ2=1.2.Figure 13
Bifurcation diagram ofp with respect to θ2 when θ1=1.2.Bifurcation diagram ofp when the control parameter u varies from 0 to 1 is given in Figure 14. When u takes the value of zero, system (20) will be the same as system (14), and then the system is chaotic. With the increase of u, the system changes from the chaotic region to the bifurcation region and finally falls into the stable region, which means the chaos behavior has been controlled completely.Figure 14
Bifurcation diagram ofp when u varies from 0 to 1.


## 5. Conclusion

Previous researches have explored the price fluctuation with price and demand game model in the supply chain management. However the retailers’ service level has not been involved in the basic demand model. This paper depicts a modified Bertrand duopoly model based on the stochastic demand and the retailer’s service level which influences the market demand decided by customers’ preference, while the market demand is stochastic in this model. How the price adjustment speed affects the stability of the supply chain system with respect to service level and stochastic demand has been studied in this paper. The dynamic behavior of the system is researched by simulation and the stability domain and the bifurcation phenomenon are shown clearly. The largest Lyapunov exponent, the chaotic attractor, and the time series are also given to confirm the chaotic characteristic of the system.We come to the conclusion that (1) with the price adjustment speed gradually increasing, the price system gets unstable and finally becomes chaotic; (2) a relatively small price adjustment speed may maintain the system at stable state; (3) chaos can be controlled completely with the parameter control method. This chaotic phenomenon will perturb the product market. Hence, the conclusion makes significant contribution to the system, referring to the price fluctuation based on the service level and stochastic demand.In spite of the contribution that this paper has offered to the managers, some limitations still exist in this paper. Firstly, the supply chain only contains two retailers; in order for this study to be more close to the reality and be used in more scales of enterprises, we may study the dynamic game behavior in the three-oligarch product market. Secondly, this paper only does the numerical simulation, while there are no real data to support the conclusion of the paper. Hence, some real data should be collected to identify the conclusion.



---

*Source: 1382689-2017-09-24.xml*