This paper analyses the dynamics of dual-channel energy supply chain model with heterogeneous retailers (as regards the type of expectations’ formation). On the basis of analyzing the stabilities of four fixed points in the three-dimensional dynamic system, local stable regions of Nash equilibrium are obtained. Effects of S on the stable regions and profit are studied. Simulation results show that the adjustment of price speed has an obvious impact on the complexity of competition. The performances of the model in different period are measured by using the index of average profit. The results show that unstable behavior in economic system is often an unfavorable outcome. So this paper discusses the application of parameters control method when the model is in chaos and then allows the oligarchs to eliminate the negative effects.
1. Introduction
A number of industrial and governmental statistical reports show that commerce on the Internet is growing at an attractive rate. It means that consumers can get product from dual-channel which is the traditional retail channel and online direct channel in the market. With the emergence and development of the online channel, much previous research on dual-channel management tends to focus on pricing strategies. Balasubramanian [1] models competition in a dual-channel environment from a strategic perspective and shows that the level of market coverage can be used as a mechanism to control competition. Chiang et al. [2] argue that the vertically integrated direct channel allows a manufacturer to constrain the partner retailer’s pricing behavior. Hua et al. [3] show that the difference between the demand transfer ratios in the two channels with respect to lead time and direct sale price and customer acceptance of the direct channel have great effects on the lead time and pricing decisions. Furthermore, some researchers have studied pricing strategies under information asymmetry [4–6].
However, most of the discussion grounds are focused on well-established and sophisticated microeconomics of retailers, which prevalently either is in a static context or assumes retailers’ perfect foresight. So the above papers are unable to explain the dynamic behaviors of the economic system. The dual-channel market describes the structure of oligarch competition model. It is the main characteristic of interdependence that makes oligarch consider the reaction of the competitor. In an oligarch competition, the firms are either output or price setters; refer to Cournot [7] and Bertrand [8]. Rand [9] may be the first man who suggested that the Cournot adjustment process may also fail to converge to a Nash equilibrium and may also exhibit cyclical and even chaotic dynamics. Puu [10–12] found that the Cournot duopoly model can appear as a strange attractor with fractal dimension, and he studied the strategy of the duopoly game. Agiza et al. [13–15] have studied the complex dynamics of abounded rationality duopoly game with a nonlinear demand function. The modification of the duopoly game depends on the strategy that the firms use, such as homogeneous and heterogeneous, and the expectations of the output the firms have to maximize, such as bounded rationality, naive expectation, and adaptive expectation. Bischi et al. [16–18] gave the general formula of the oligopoly model with a form of bounded rationality. They discussed the global and local stability of the duopoly game with a particular form of bounded rationality. Ma and Pu [19] studied complex behaviors of a Cournot-Bertrand mixed duopoly model with the application of nonlinear dynamics theory. Wang and Ma [20], based on the players with bounded rationality, proposed a Cournot-Bertrand mixed game model and discussed the stability of the system. Aims of all the modifications previously mentioned or not are to make oligarch competition model become economically more justified in the world. In recent years, several scholars have considered more realistic questions, such as Ma et al. [21, 22] who studied bounded rational oligarch game from insurance and R&D input competition market; they found the chaos phenomenon and proposed some control methods.
In our paper, we considered one online direct retailer and one dual-channel retailer in the market. We assumed that the dual-channel retailer made price decision through both network and traditional channel. Based on the analysis of this situation, we established the corresponding dynamic price game model. Local stable regions of Nash equilibrium and effects of S on the stable regions and profit are studied. Adjustment of price speed has an obvious impact on the complexity of competition. The results show that unstable behavior in economic system is often an unfavorable outcome.
So the rest of this paper is as follows: In Section 2, we describe the model of dual-channel Bertrand Game with heterogeneous expectations. In Section 3, we study the local stability of equilibrium points in the model. In Section 4, we do simulation. In Section 5, the control measure is adopted. Finally, some conclusions are drawn in Section 6.
2. Model Framework2.1. Assumptions
The model is based on the following assumptions:
The dual-channel retailer and online direct retailer sell similar products from different channels on the basis of price competition.
The two retailers are with heterogeneous expectations. The dual-channel retailer is bounded rationality and the online direct retailer is thinking with adaptive expectations.
The dual-channel retailers and online direct retailer have the same cost.
2.2. Model Construction
We consider two oligopoly retailers in supply chain: one retailer sells products to consumers through network and the traditional store, while the other sells products to consumers only through the Internet. The structure of supply chain is as shown in Figure 1. Different from the classic Bertrand model which is a special game model with price as its decision varies, the dual-channel retailer makes price decision both in the tradition channel and network channel. The online direct retailer makes price decision in network channel.
The structure of dual-channel energy supply chain.
The demand function is given by(1)qd0=sa-bp0+β1p1+β2p2,qd1=1-sa-bp1+β0p0+β2p2,qn=1-sa-bp2+β1p1+β0p0.
Here subscripts d, n represent the dual-channel retailer and the online direct retailer, respectively. a(a>0) represents the possible largest demand.
s(0<s<1) represents the percentage share of the demand going to traditional channel. b is the coefficient of price elasticity of demand. βi is cross-price coefficient which represents the retail price sensitivities of the demands in the traditional channel and direct channel. In reality, βi is related to channel loyalty. To make the discussion easier to follow, we let β0=β1=β2=β.
The cost function is linear, so the profit function of dual-channel retailer is(2)πd=p0qd0+p1qd1-cqd0+qd1.
The profit function of online retailer is(3)πn=p2-cqn.
Based on first-order conditions of the profit function ∂π/p=0, we can calculate the optimal reaction function of each retailer, and the intersection of reaction function (4) is the Bertrand Nash equilibrium point:(4)p0=2βp1+βp2+as+cb-β2b,p1=2βp0+βp2+a1-s+cb-β2b,p2=bc+βp0+βp1+a1-s2b.
As the price of the game among retailers is ongoing, price game process is repeated in the complex process for a long term. In oligopoly, game players can choose simple expectation rules such as naive or complicated as adaptive expectations and bounded rationality. The players can use the same strategy (homogeneous expectations) [23, 24] or can use different strategy (heterogeneous expectations). In this study we consider heterogeneous retailers such that each retailer thinks with different strategy to maximize his outputs.
We assume the dual-channel is bounded rationality player [17, 25, 26]. In reality, he could only get part of market information; hence they try to use local information based on estimation of marginal profit in the current period. That is, if the current marginal profit is positive, the retailer will raise his price in the next period. Otherwise, he will reduce it. The price adjustment process is(5)pit+1=pit+αipit∂πip0,p1,p2∂pit,0≤αi≤1,i=0,1,where αi is a positive parameter which represents the relative speed of price adjustment.
The online direct retailer can revise their beliefs according to the adaptive expectations rules. If the online direct retailer thinks with adaptive expectations [27, 28], its outputs are obtained with weights between last period’s price decision and his reaction function (4). Hence the dynamic equation of the adaptive expectation retailer has the form(6)p2t+1=p2t+vp2t∗-p2t,0≤v≤1,where v is the speed of adjustment parameter of an adaptive player. If v=1 then (6) of adaptive expectations is identical to naive expectation. This implies that naive expectation is a special case of adaptive expectations behavior.
Then we get a discrete dynamic system of the dual-channel game model as follows:(7)p0t+1=p0t+α0p0t2βp1t-2bp0t+βp2t+as+cb-β,p1t+1=p1t+α1p1t2βp0t-2bp1t+βp2t+a1-s+cb-β,p2t+1=p2t+vbc+βp0t+βp1t+a1-s2b-p2t.
3. Model Analysis3.1. Equilibrium Solutions
In system (7), let pi(t+1)=pi(t); then(8)p0t2βp1t-2bp0t+βp2t+as+cb-β=0,p1t2βp0t-2bp1t+βp2t+a1-s+cb-β=0,bc+βp0t+βp1t+a1-s2b-p2t=0.
We can get the equilibrium solutions. Before we solve for the equilibrium solutions of system (7), firstly, we assign some values to parameters by considering actual competition. Let a=4, b=0.7, s=0.6, c=1, β=0.3. Then we get the four equilibrium solutions of system, which are(9)E1=0,0,1.643,E2=2.465,0,2.171,E3=0,1.866,2.043,E4=4.509,4.109,3.489.Obviously, E1, E2, and E3 are boundary equilibrium solutions, while E4 is the only Nash equilibrium solution.
3.2. Stability of Equilibrium Solutions
To study the stability of the equilibrium solutions, we will calculate the Jacobian matrix of system (7):(10)JE=U12βα0p0βα0p02βα1p1U2βα1p1βv2bβv2b1-v,where (11)U1=1+α02βp1-2bp0+βp2+as+cb-β-2bα0p0,U2=1+α12βp0-2bp1+βp2+a1-s+cb-β-2bα1p1.
Then we put the value of E1, E2, and E3 into (10); according to the conditions for the stability of equilibrium solutions, the moduli of all characteristic roots should be less than 1 [29]. So we can see that E1, E2, and E3 are unstable equilibrium solutions.
At the Nash equilibrium point E4=(4.509,4.109,3.489), its Jacobian matrix is(12)JE4=1-6.313α02.705α01.353α02.465α11-5.753α11.233α10.2143v0.2143v1-v.And the characteristic equation of J(E4) is λ3-Aλ2-Bλ-C=0, where(13)A=6.31α0+5.75α1+v-3,B=6+6.02vα0-2v-5.49vα1-12.63α0-11.50α0+29.65α0α1,C=6.31α0+5.75α1-6.02vα0-5.49vα1+24.89vα0α1+v-1-29.65α0α1.
According to the Jury stability criterion, the necessary and sufficient condition of asymptotic stabilization at E4 is that all zero points of the characteristic polynomial are inside the unit circle of the complex plane. So the following four conditions must be met:(14)f1=1-A-B-C>0,-f-1=-A+B-C+1>0,1-C2>0,1-C22-B+AC2>0.
Conditions (14) define a region of system stability in the place of the adjustment parameters which can be shown in Figure 2(a). Equations (14) are conditions for local stability of the Nash equilibrium E4 which is asymptotically stable with the values of (α0,α1,v). In the stable region, the price will reach the Nash equilibrium point by modulating limited times with random initial price. If the values of (α0,α1,v) are out of the stable region, the market will lose stability or even fall into chaos.
As can be seen in Figure 2, the system is stable in v∈[0,1], which means v almost has no influence on the stability of system. The market is stable when v∈[0,1] and the values of αi are small, but the market may be unstable when αi increases. In the stable region shown in Figure 2, ranges of α0 and α1 are almost the same.
The economic meaning of the stable region is that if the value of (α0,α1,v) is in the stable region, prices of three firms, p1, p2, and p3, will achieve the Nash equilibrium price after a number of games. Firm 1 and firm 2 should control their price adjustment speed to maintain the market stability.
3.3. Effect of S on the Stable Regions and Profit
In order to show the effects of the percentage share of the demand in traditional channel on the stable regions, we set S=0.3, 0.5, 0.7, respectively; according to the parameters above, we can get Figures 2(b), 2(c), and 2(d).
As can be seen from the comparison of Figures 2(a)–2(d), with the increase of the percentage share of the demand in traditional channel S, ranges of α1 expand while ranges of α0 narrow. From the perspective of economics, in order to keep the market stable, the smaller the percentage share of the demand in traditional channel (S), the larger the rang of the price adjustment speed of the traditional channel firm (α0) and the smaller the rang of the price adjustment speed of the network channel firm (α1).
When we set S=0.3, 0.5, 0.6, 0.7, respectively; the Nash equilibrium E4 is also different which can be shown in Table 1.
Nash equilibrium and Nash profit equilibrium with different S.
S
0.3
0.5
0.6
0.7
Nash equilibrium
p0=4.2915;p1=5.0915;p2=4.5106
p0=4.4362;p1=4.4362;p2=3.8298;
p0=4.5085;p1=4.1085;p2=3.4894
p0=4.5809;p1=3.7809;p2=3.1489;
Nash profit equilibrium
πd=11.22;πn=8.63
πd=9.45;πn=5.61
πd=8.84;πn=4.34
πd=8.41;πn=3.23
According to the profit function of dual-channel retailer and the profit function of online retailer, πd=p0qd0+p1qd1-c(qd0+qd1) and πn=(p2-c)qn. The Nash profit equilibrium is also shown in Table 1.
As can be seen from the comparison of Figures 2(a)–2(d), with the increase of the percentage share of the demand in traditional channel (S), Nash equilibrium price (p0) in traditional channel increases; Nash equilibrium prices (p1,p2) in network channel decrease.
With the increase of the percentage share of the demand in network channel (1-S), Nash equilibrium profit of all the retailers increases.
3.4. Parameter Basin Plots
Figures 3, 4, and 5 present the parameter basin plots in the (α1,α2) with different value of parameter v, in which different colors represent different states. The light blue indicates stable steady state and red is for cycles of period 2, green for cycles of period 3, yellow for cycles of period 4, black for cycles of period 5, pink for cycles of period 8, dark blue for chaos, and white for divergence (which means one of the players will be out of the market in economics). We can find that when the parameters (α1,α2) from the light blue area pass through yellow, pink, and dark blue areas in turn, system (7) enters into chaos through flip bifurcation. When (α1,α2) from the black go to the dark blue area directly, system (7) enters into chaos through Neimark-Sacker bifurcation.
The parameter basin when v=0.2.
The parameter basin when v=0.5.
The parameter basin when v=0.8.
4. Numerical Simulations for System
The main goal of this section is to show that the complicated dynamic features of the dual supply model with heterogeneous retailers. To provide some numerical evidence for the chaotic behavior of system (7), we present various numerical results here to show the chaoticity, including its bifurcations diagrams, Lyapunov exponents, strange attractors, and sensitive dependence on initial conditions. In order to study the properties of the dynamic model, it is convenient to take the parameters values assigned above.
4.1. The Effect of the Price Adjustment Speed
Figure 6 describes the Lyapunov exponents and price bifurcation with change of α0 when α1=0.2, v=0.2. As shown in Figure 3, the Nash equilibrium point is stable at α0<0.2508. With α0 increasing, the stability of equilibrium point changes; the system experiences period-doubling bifurcation and falls into chaos eventually. The positive Lyapunov exponent is used to mark the chaos. The system is in chaos when most of the Lyapunov exponents are positive. Figure 7 shows a chaos attractor of system (7), which is another characteristic of the system. The sensitive dependence on initial conditions is one of the important features of chaos. Figure 8 shows the changes of retailers different price separately with process of games when (p0,p1,p2)=(1,1,1) and (p0,p1,p2)=(1.001,1,1). From Figure 8 we can see that when p0 changes a little, the price decision will have distinct changes. Now we have studied the dynamic behavior in system (7) with change of v. Figure 9 shows that the changes of parameter v have no effect on system (7) when α0=0.2, α1=0.2; this means that the Nash equilibrium point is stable for 0≤v≤1.
Lyapunov exponent and bifurcation diagram with α0 when α1=0.2,v=0.2.
Chaotic attractor of the system when α0=0.42, α1=0.2.
Sensitivity dependence to initial condition when system in chaos period.
Bifurcation diagram with v when α1=0.2, α2=0.2.
4.2. The Effect of Parameters Changes on the Profit
Figure 10 shows the profit of the two retailers in 50 games when the system is in a stationary period, two period-doubling bifurcations, and a chaotic period. Figure 11 shows the average profits of both retailers with respect to the change of α0. From the figure, we see that when the price adjustment speed increases, the system goes through different stages. So we can make conclusions that the profit of the two retailers in a chaotic period is less than that in other periods. These explain why chaos is undesirable in the real market.
The profits of retailers when system in (a) stable period; (b) two period-doubling bifurcations; (c) chaotic period.
The change of average profits with respect to α0.
5. Chaos Control
Our study found that once the market is in chaos, the total profit of the system is less than that in the equilibrium state, so this kind of state is not expected to appear. However, the current situation is that the oligarchs often maximize profit by any kind of means in the process of marketization. So the market will be out of order and finally fall into chaos. It is particularly important for each party that some control measures need to be adopted timely to make the system return to its stable equilibrium. Therefore, we propose the adaptive adjustment mechanism [30, 31] to control system (7). Assume system (7) is generally described by pi(t+1)=fi(p0(t),p1(t),p2(t)). Then the controlled system is(15)p0t+1=1-τg0p0t,p1t,p2t+τp0t,p1t+1=1-τg1p0t,p1t,p2t+τp1t,p2t+1=1-τg2p0t,p1t,p2t+τp2t,where τ is the adjustment parameter. Figure 12 shows the bifurcation diagram of price against the adaptive parameter τ. From Figure 12, with control parameter τ increasing, chaotic system was gradually controlled at Nash equilibrium solutions. Figure 13 shows the cases of control parameters τ=0.2 and τ=0.35, under which the system rapidly converges to aperiodic-2 orbits and the Nash equilibrium solutions, respectively. In real market, we can consider τ a learning ability or adaptability of market. For instance, the retailers analyzed the information in the past and then adopted proper price adjustment speed.
Bifurcation diagram with change of the control parameter τ.
Time history of prices of the retailers (a) τ=0.2 and (b) τ=0.35.
6. Conclusions
In this study, a dual-channel supply chain model with a traditional retailer and an online direct retailer is proposed.
Highlights of innovation included in the paper are the following 2 points:
Performances of the model in different period are measured by using the index of average profit.
Unstable behavior in economic system often brings an unfavorable outcome.
The bifurcations, chaos, and other complex phenomena were found and the chaotic features are justified numerically via computing Lyapunov exponents, showing the chaotic attractor and sensitive dependence on initial conditions. When chaos occurs, it breaks the whole system, and the market will become abnormal, irregular, and unpredictable. The system’s performance was measured by the profits. With the increase of the percentage share of the demand in network channel (1-S), Nash equilibrium profit of all the retailers increases. It can be found from our study that the optimal strategy is in the stable region of the Nash equilibrium. Once the market is in chaos, the profits of system are less than that of the system in equilibrium state. At last, we use adaptive adjustment mechanism to control the system, showing that the controlled system can delay chaos or even eliminate chaos completely. So it is particularly important for each retailer to adopt some control measures timely in order to make the system keep in stable state and avoid the loss of system profits.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research was supported by the National Natural Science Foundation of China (no. 61273231) and Doctoral Fund of Ministry of Education of China (Grant no. 20130032110073) and supported by The Foundation for Outstanding Young Scientist in Shandong Province (Fund no. BS2014SF021).
BalasubramanianS.Mail versus mall: a strategic analysis of competition between direct marketers and conventional retailers199817318119510.1287/mksc.17.3.1812-s2.0-0032218131ChiangW.-Y. K.ChhajedD.HessJ. D.Direct marketing, indirect profits: a strategic analysis of dual-channel supply-chain design200349112010.1287/mnsc.49.1.1.127492-s2.0-0037287242HuaG.WangS.ChengT. C. E.Price and lead time decisions in dual-channel supply chains2010205111312610.1016/j.ejor.2009.12.0122-s2.0-75849127286YaoD.-Q.LiuJ. J.Competitive pricing of mixed retail and e-tail distribution channels200533323524710.1016/j.omega.2004.04.0072-s2.0-11144253598YueX.LiuJ.Demand forecast sharing in a dual-channel supply chain2006174164666710.1016/j.ejor.2004.12.0202-s2.0-33845269555YanR.GhoseS.Forecast information and traditional retailer performance in a dual-channel competitive market2010631778310.1016/j.jbusres.2009.02.0172-s2.0-70349833151CournotA.1838Paris, FranceHachetteBertrandJ.Review of Theorie mathematique de la richesse sociale and of Recherches sur les principles mathematiques de la theorie des richesses188367499508RandD.Exotic phenomena in games and duopoly models19785217318410.1016/0304-4068(78)90022-8MR5114862-s2.0-0000374060PuuT.The chaotic duopolists revisited1998333-438539410.1016/S0167-2681(97)00064-42-s2.0-0032517943PuuT.Complex dynamics with three oligopolists19967122075208110.1016/S0960-0779(96)00073-22-s2.0-0001327121PuuT.On the stability of Cournot equilibrium when the number of competitors increases2008663-444545610.1016/j.jebo.2006.06.0102-s2.0-42649095987AgizaH. N.HegaziA. S.ElsadanyA. A.Complex dynamics and synchronization of a duopoly game with bounded rationality200258213314610.1016/s0378-4754(01)00347-0MR18740612-s2.0-0036133248AgizaH. N.ElsadanyA. A.Chaotic dynamics in nonlinear duopoly game with heterogeneous players2004149384386010.1016/S0096-3003(03)00190-5MR2033167ZBL1064.910272-s2.0-0345867023AgizaH. N.ElsadanyA. A.El-DessokyM. M.On a new Cournot Duopoly game20132013548780310.1155/2013/487803BischiG. I.GardiniL.KopelM.Analysis of global bifurcations in a market share attraction model2000245–785587910.1016/s0165-1889(99)00028-7MR1748840BischiG. I.NaimzadazA.Global analysis of a dynamic duopoly game with bounded rationality20005Boston, Mass, USABirkhäuser361385Annals of the International Society of Dynamic Games10.1007/978-1-4612-1336-9_20BischiG. I.SbragiaL.SzidarovszkyF.Learning the demand function in a repeated Cournot oligopoly game200839440341910.1080/00207720701792131MR23986782-s2.0-40249113324MaJ.PuX.The research on Cournot-Bertrand duopoly model with heterogeneous goods and its complex characteristics201372489590310.1007/s11071-013-0761-7MR30645082-s2.0-84878560876WangH.MaJ.Complexity analysis of a Cournot-Bertrand duopoly game model with limited information20132013628737110.1155/2013/287371MR3037742JunhaiM.JunlingZ.Research on the price game and the application of delayed decision in oligopoly insurance market20127042327234110.1007/s11071-012-0566-0MR29981652-s2.0-84870859412LiT.MaJ.The complex dynamics of R&D competition models of three oligarchs with heterogeneous players2013741-2455410.1007/s11071-013-0947-z2-s2.0-84883762111Dubiel-TeleszynskiT.Complex dynamics in Bertrand duopoly games with homogeneous and heterogeneous players20092295116FurthD.Anything goes with heterogeneous, but not always with homogeneous oligopoly200933118320310.1016/j.jedc.2008.04.012MR24776832-s2.0-56949103265AgizaH. N.HegaziA. S.ElsadanyA. A.The dynamics of Bowley's model with bounded rationality20011291705171710.1016/S0960-0779(00)00021-7ZBL1036.910042-s2.0-0035400910ElsadanyA. A.Competition analysis of a triopoly game with bounded rationality201245111343134810.1016/j.chaos.2012.07.0032-s2.0-84866500820BischiG. I.KopelM.Equilibrium selection in a nonlinear duopoly game with adaptive expectations20014617310010.1016/S0167-2681(01)00188-32-s2.0-0001902879DingJ.MeiQ.YaoH.Dynamics and adaptive control of a duopoly advertising model based on heterogeneous expectations201267112913810.1007/s11071-011-9964-yMR28687882-s2.0-82255164400LiC.PengJ.HuangQ.Controlling synchronization of chaos and hyperchaos systems20013334347HuangW.Stabilizing nonlinear dynamical systems by an adaptive adjustment mechanism2000612R1012R10152-s2.0-0000173216HuangW.Theory of adaptive adjustment20005424726310.1155/S102602260000056X