Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 2014, Article ID 749769, 14 pages http://dx.doi.org/10.1155/2014/749769

Research Article

The Complex Dynamics of Bertrand-Stackelberg Pricing Models in a Risk-Averse Supply Chain

Junhai Ma and Qiuxiang Li

College of Management and Economics, Tianjin University, Tianjin 300072, China Correspondence should be addressed to Qiuxiang Li; lqxkycg@163.com Received 22 September 2013; Accepted 17 February 2014; Published 18 May 2014 Academic Editor: Mingshu Peng

Copyright © 2014 J. Ma and Q. Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We construct dynamic Bertrand-Stackelberg pricing models including two manufacturers and a common retailer in a risk-averse supply chain with the uncertain demand. The risk-averse supply chain follows these strategies: Bertrand game between the two manufacturers and Stackelberg game between the manufacturer and the retailer. We study the effect of the price adjustment speed, the risk preference, and the uncertain demand on the stability of the risk-averse supply chain using bifurcation, power spectrum, attractor, and so forth. It is observed that there exists slip bifurcation when the price adjustment speed across some critical value, the stable region, and total profit of the risk-averse supply chain will increase with increase of RMi and decrease with increase of a. The profit of the supply chain and the two manufacturers will decrease and the weaker (retailer) is a beneficiary when the supply chain is in chaos. The fluctuation in the supply chain can be gradually controlled by the control of the price adjustment speed.

1. Introduction

With the development of economic globalization, the relationship among the supply chain members becomes more and more complex under the different environment. An enterprise which is involved in the middle of multiple supply chains has all kinds of complicated relationship when the parameter of the upstream and downstream enterprises is changed, such as market uncertainty and risk preference factor; the enterprise's decision-making behaviors become more complicated and hard to predict.

We all know that the price is always a sensitive topic, which can affect customer needs and wants, distribution of the products and services among the supply chains. Scholars at home and abroad have done a lot of research on this aspect. Wei et al. [1] studied pricing decisions in a supply chain with two manufacturers and one common retailer and constructed five pricing models under decentralized decision cases with consideration of different market power structures. Mukhopadhyay et al. [2] considered two separate firms, which had private forecast information about market uncertainties and offered complement goods in a leader-follower type, and devised a "simple to implement" information sharing scheme under which both firms and the total

system are better off. These literatures analyzed and compared the optimal solution under different market power structures, but they did not consider the effect of risk reference of the participants on the optimal solution.

There are many literatures taking risk preference into account. Caliskan-Demirag et al. [3] constructed models of the supply chain with a risk-averse retailer by adopting the conditional-value-at-risk (CVaR) decision criterion. Luo and Huang [4] explored the impact of risk preference on the strategies of the supply chain by taking the different attention of high profit and low profit as retailers risk measurement. These literatures only consider the unilateral risk, but, in the uncertainty environment, there exists bilateral risk among the participants which corresponds to the actual situation. In this paper, we will consider a supply chain under bilateral risk with two risk-averse manufacturers and a risk-averse retailer, which make the supply chain more complex.

Many literatures study the complexity of supply chain. Huang and Chen [5] studied the sale-surety contract option of supply chain with effort dependent demand and risk preference and got some meaningful conclusions. Guan and Zhou [6] researched the integrated optimization problem

of three-level supply chain consisting of suppliers, distributors, and retailers under decision-makers having different risk attitude. Huang and Yang [7] studied a two-echelon supply chain model with one supplier and one retailer in a newsvendor problem; the supplier with different risk attitude has great influence on the retailer's optimal order quantity; the operation efficiency of the supply chain will be underperformed when the supplier is much too risk-averse. These literatures have studied the participant's behaviors of the risk supply chain, but they did not present the dynamic complex features of the risk supply chain.

Research on dynamical complexity of a system has been of concern to scholars. Puu [8] found that the Cournot three oligopoly model could appear strange attractors with fractal dimension, and he studied the situation of the duopoly game. Many researchers improved the classical Cournot model and found that certain dynamical behaviors of the system occurred in the course of repeated games with three or four duopolies. Many experts have also studied this field under different conditions, such as different expectations and incomplete information select it Ma and Sun [9] established a decentralized pricing game model and studied its complex dynamic characteristics of triopoly under different decision-making rule; the result showed that the process of game would tend to a Nash equilibrium at a lower price adjustment speed, and, with the increase of the value of adjustment speed, the system would appear to be unstable and gradually gone into a chaos state. Ma and Bangura [10] studied the dynamic complexity of financial and economic system under the condition of three parameters changing.

In recent years, many experts apply the dynamical complexity to study the dynamic change process of supply chain. Hwarng and Xie [11] found that there existed the chaotic enlargement phenomenon among the members of the supply chain which enriched the connotation of the bullwhip effect. J. Wang and X. Wang [12] established nonlinear supply chain inventory system models with forbidden returning and limited supply capacity; numerable simulations showed that the supply chain inventory system had complex dynamic behaviors under certain parameter settings; they gave some suggestions to eliminate the complexity of the dynamic supply chain. Ma and Feng [13] presented investigation simulations of retailer's demand and stock; the behaviors of the system exhibited deterministic chaos with consideration of system constraints. These literatures researched the dynamic complexity of the supply chain but did not consider the influence of the decision-maker's risk behaviors on the supply chain management. In this paper, we will study the dynamic complexity of a risk-averse supply chain with two manufacturers and a common retailer under uncertain demand. Considering the change of parameters in the dynamic risk supply chain, such as the price adjustment speed, risk preference, and uncertain demand, we can study the influence of parameters on the price and stable region of the two manufacturers and retailer.

The remainder of this paper is organized as follows. In Section 2, we describe the supply chain problem, make assumptions of the system model, and discuss the system model. In Section 3, we construct a Bertrand-Stackelberg

dynamic pricing model which consists of two manufactures and one retailer with risk-averse attitude. Analysis is made under different variable conditions in Section 4. In Section 5, the variable feedback control method will be used to control chaos in the system. In Section 6, we outline some conclusions and hence relevant recommendations for future research.

2. Model

2.1. Description of the Problem. In this section, we construct a dynamic pricing game model in a risk-averse supply chain which consists of two manufacturers (M1 and M2)and acom-mon retailer (R). The two manufacturers are competitive and sell respective products to the common retailer, and the common retailer sells two kinds of products to consumers directly. The customers' demand is stochastic. We consider the supply chain following these strategies: Bertrand game between the two manufacturers and Stackelberg game between the manufacturer and the retailer. In these strategies, the two manufacturers and the retailer make their own decisions, respectively, for maximizing their profit; the decision process is as follows: the two manufacturers, as the Stackelberg leader, determine the respective wholesale price (wt) (i = 1,2); the retailer as the follower sets his own optimal retail price (pt) (i = 1,2) based on the manufacturer's decisions.

Furthermore, in order to capture the uncertain demand which is affected by the change of economic and business conditions and prediction errors, we assume the market demand random variable a is as follows: a = a + e, where a is the primary demand level and e follows a normal distribution such as E(e) = 0, Var(e) = a2. However, the normality assumption has been used extensively in the literature (e.g., Gal-Or [14]; Raju and Roy [15]; Vives [16]). The two manufacturers and the retailer know the distribution of the uncertain demand and determine their behaviors, respectively.

Because the customer demand is stochastic, there is financial risk to the two manufacturers and the retailer. Therefore, we should consider the effect of the riskpreference of the two manufacturers and the retailer on pricing decision. The preference theory provides the framework which incorporates the participators' financial risk preference into their decision process. The valuation measure we use is known as the certainty equivalent in the preference theory and is defined as certain value that a participator is just willing to accept an uncertain event (Kunstman [17]).

One form of the utility function in both theoretical and applied work in areas of decision theory and finance is the exponential utility function which can be expressed as (n) = -e-*i/R> (i = M1, M2, R), where Rt is the risk tolerance level of the two manufacturers and retailer, nt is the profit, and e is the exponential constant. When Rt < >x>, it implies that the decision-maker has risk-averse behavior, and Rt approaches >x> which implies the decision-maker is risk-neutral (Walls [18]). If the decision-maker is risk-averse, n follows a normal distribution, and expected utility is

E(U) = E(n) - (Var(n)/2R), where E(n) is the mean of n and Var(n) is the variance of n.

2.2. Assumption of the System

(1) Customer demand is always satisfied, demand function is linear, and the two manufacturers and the retailer make decentralized decision.

(2) We consider two partly substitutable products coming from a competitive market in which consumers can buy any one of them.

(3) The consumer demand is stochastic; the two manufacturers and the retailer are all risk aversion.

(4) CM is marginal cost of M^ CM is marginal cost of

E(Um1 ) = E(nMi )- J" Ml ^

( \f- U ,r \ (wi -CmI )2°2

= (wi - CM, )(a-biPi + dip2)---

e(um2) = e(vm2)- 2[ Ml) 2rm2

= {W2 -Cu2 )(a-b2p2 + d2pi)

(w2 -CM2 fo2

From formulas (2), (3), and (4), we obtain the revenue functions of the two manufacturers and the retailer which is more in accordance with the actual situation using the utility function. When a manufacturer changes the value of parameter, how to adjust value of parameter and what impact it will have on the other manufacturers and the retailer are the main innovation points of this paper.

2.3. Revenue Function of the System. In this study, wt and pt (i = 1,2) are decision variables and other variables are exogenous variables. As is known in the supply chain, we assume that pt > wt (i = 1,2); this inequality ensures that each participant can obtain a positive profit.

Extend the demand function in Banker et al. [19]. We assume that the primary demand function in this paper is decided by p; (i = 1,2) as follows:

Di =a-bipi +dip2, D2 = a-b2p2 +d2pi,

where a represent the base demand level of the product, b (i = 1,2) are price sensitive coefficient of demand, the cross-price sensitive coefficient dt (i = 1,2) reflects the substitution degree of the products, and a,bt,dt > 0. We can obtain the expected utility functions of the two manufacturers and the retailer as follows:

E(UR) = E(nR)-

= (Pi - Wi) (a - biPi + dip2) + (P2 - W2)

(Pi + p2 -Wi - W2)2°2

x (a-b2p2 + d2pi) -

3. Bertrand-Stackelberg Model

Suppose that the two manufacturers and the retailer have principal and subordinate relationship, the two manufacturers are Stackelberg leaders, the retailer is follower, and there is Bertrand competition between the two manufacturers. Then, the manufacturers and the retailer process sequential dynamic game; the game equilibrium is Stackelberg equilibrium. In this game, the two manufacturers make decisions for wholesale price according to the market information; the retailer makes decisions according to the two manufacturers. Using backward induction, we first find the response functions of the second stage from the game model. The optimal marginal utility of the retailer can be obtained by the firstorder conditions of formula (2); the calculation results are as follows:

dE(UR) (Wl +w2)a2

= a + blwl - a2w2 +

0 +2bi)pi + [di + d2 )P2,

dE(UR)_ (wi +w2)o1

—--= a + b2w2 - diwi +--

0 +2b2)p2 +(di +d2 )pi.

The retailer's reaction functions are as follows by solving formula (5):

pi = *

AE + BC DE-C2 ' AC + BD DE-C2 '

where A = a + b1w1 - d2w2 + ((w1 + w2)a2/Rr),B = a + b2w2 -d1w1 + ((w1 + w2)a2/RR), C = d1 + d2 - (a2/RR), D = (o2/Rr) + 2b1,E = (o2/Rr) + 2b2.

Formula (6) is the optimal decision making of the retailer on the premise of w1,w2; the retailer can obtain the decision after it observes the manufacturer's behavior. Substitute formula (6) into formulas (3) and (4); the optimal wholesale price of the two manufacturers can be obtained by the firstorder conditions of formulas (3) and (4):

dE(UMi) dw1

= a - bipi + dip2 -

(wi -cu1 )a

RMj (w1 - CMJ )

(a2/RR + 2b1) (g2/Rr + 2b2) - (d1 + d2 - o2/Rr)2

+ \ d, + d2 -

+ d2 -

+ \T + 2b1

dE(UM2)

= a - b2p2 + d2p1 -

(w2 -cm2)a

(w2 - CU2 )

(a2/RR + 2b1) (g2/Rr + 2b2) - (d1 + d2 - o2/Rr)2

+ \d1+d2 (b2 +a

41 + d-TJ\TR-"2

+ ( +2b1 )(b2 + a

Let dE(UMj )/dw1 = 0 and dE(UM2 )/dw2 = 0; we can find out the equilibrium solutions of the two manufacturers to the retailer at some stage; the equilibrium solutions express optimal decision of the two manufacturers in various possible situations in a game stage.

In the actual decision process, the economic behavior of the decision-maker always shows limited rational characteristics, such as risk-averse behavior. Then, decision results are different from the one which is perfectly rational. Optimal solution is an optimal state when other parameters in the system change. Although it is very much rare for the system to go into an optimal state, how to determine the adjustment speed and pricing orientation is the focus of this research. In this paper, the two manufacturers make decisions based on limited rational expectations; they adjust the game process on the basis of last period marginal utilities. If the marginal utilities in period t are positive, they will continue their output adjustment strategy in period t +1. The process can be modeled as follows:

dE(UM )

w1 (t+1) = w1 (t) + k1w1 (t)- v j/

P1 (t) =

w2 (t+1)=w2 (t)+k2w2 (t)

a + b1w1 (t) - d2w2 (t)

+ (w1 (t)+w2 (t))a2 Rr

a + b2w2 (t) - d1w1 (t)

+ (w1 (t)+w2 (t))a2' R

R.+2t\

d, + di -

dw1 dE(UM2)

+ d2 -

p2 (t) =

a+ b1w1 (t) - d2w2 (t)

(w(t) + w2 (t))o2 ' Rr .

a + b2w2 (t) - d1w1 (t) (w1 (t)+w2 (t))a2] I

1 + d2 -

xi +2b1 )(i+2b2

+ d2 -

45 40 35 30 25 20 15 10 5 0

X: 0.00102 Y: 31.22

X: 0.00103 Y: 26.33

■■ U/j-iVC

■ " V- 'jt

44 42 40 38 <p 36 34 32 30 28

0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5

0.2 0.4 0.6

0.8 k1

1 1.2 1.4 1.6 x10-

Figure 1: Price bifurcation and Lyapunov of two manufacturers and retailer with change of kl when k2 = 0.001: (a) wholesale price; (b) retail price; (c) Lyapunov.

where k1 is the adjustment coefficient of w1 and k2 is the adjustment coefficient of u>2. According to the dynamic adjustment process, we can see that the wholesale price of the manufacturer is related to the price adjustment speed, retail price, the mean and variance of the base demand level, and risk tolerance level ofmanufacturers. Similarly, the retail price of the retailer is related to the adjustment speed coefficient, wholesale price, the mean and variance of the base demand level, and risk tolerance level of the retailer.

4.1. The Fixed Point. In system (8), letting + 1) = u>;(i)(i = 1,2), we can get the fixed points of system (8). Before we solve the fixed points of system (8), we first assign some parameters considering the actual competition: a = 1000, bx = 1.2, fc2 = =

0.8, d2 = 0.7, a = 60, Rr = 60, £

Mj = 15, and CM = 10. We will calculate all the fixed points and only consider the Nash equilibrium point = 31.22, w2 = 26.33, p1 = 37.37, and p2 = 36.59).

Jacobian matrix of (8) in the Nash equilibrium point is

4. The Complex Dynamic Behavior

The ultimate goal of the supply chain is to pursue profit maximization for each of the participants and to achieve optimum overall. Therefore, they should adjust price based on their marginal profit of last period.

/1 + 91895.38^ 146744.3fc1 . ' = ' 202024.66fc2 1 + 119860.3fc2 (9)

Characteristic polynomial of (9) is/(A) = A2 -AA + £, where A = 91895.38fc1 + 119860.3fc2 + 2 and B = 91895.38fc1 + 119860.3fc2 + 10718148142.27^ + 1.

45 40 35 30 25 20 15 10 5 0

50 45 40 35 p 30 25 20 15 10

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 ki and ^2

X: 0.001 Y: 37.37

X: 0.001 Y: 36.59

'<Z2i IS

■ -s-S"*

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 k1 and k2 x1°-3

Figure 2: Price bifurcation of two manufacturers and retailer with change of kl and k2: (a) wholesale price; (b) retail price.

50 45 -40 35 30 4

X: 31.22 Y: 37.37 Z: 36.59

P 35 p1 30

50 40 -30 -20 -10

25 25 1 0 0

(a) (b)

Figure 3: The chaos attractor of the supply chain system: (a) kl = k2 = 0.001; (b) kl =k2 = 0.0015.

45 40 35 30 25 20 15 10 5 0

X: 0.00102 Y : 45.09

X : 0.00102 Y : 45

Figure 4: Bifurcation diagram of wl,w2,pl, and p2 with change of kx when k2 = 0.001,Rr = 120: (a) w1tw2; (b) pt,p2.

60 55 50 45 40 35 30 25 20 15

Figure 5: Bifurcation diagram of wl,w2,pl, and p2 with change of k1 when k2 = 0.001 ,RMl = 120: (a) w1,w2; (b) p1,p2.

32 30 28 26 24 w 22 (-■ 20 18 16 14

X: Q.QQQ85 Y: 26.98

X: Q.QQQ87 Y: 22.Q6

•il!:»". :' ■ Sill«,..* *

lilllliihijlijl

36 34 32 p 30 28 26 24

X: Q.QQQ85 Y: 31.14

Q.QQQ85 29.99

- !.!!: ./I!

.■¡ii

*i!'i I.!!,

•■• .-.raw ■■■■. ■■ ■;•■: "te . " v, »

8 9 10 11 12 13 4 5 6 7 8 9 10 11 12 13

k1 x1Q-4 k1 x1Q-4

(a) (b)

Figure 6: Bifurcation diagram of wl,w2,pl, and p2 with change of kx when k2 = 0.001, a = 70: (a) wx ,w2; (b) px,p2-

The local stability of Nash equilibrium can be gained according to Routh-Hurwiz's condition:

/(1) = 1 - A + B > 0, F (-1) = 1+A + £>0, (10)

B - 1 < 0.

Condition (10) gives the necessary and sufficient conditions of stable region of the Nash equilibrium point. Economic meaning of the stable region is that, whatever initial price is chosen by the two manufacturers in local stable region, they will eventually achieve Nash equilibrium price after finite games. It is important to notice that the two manufacturers may accelerate the price adjustment speed in order to increase their profit. Price adjustment parameter does not change

Nash equilibrium point. Once one manufacturer adjusts price too fast and pushes ku k2 out of the stable region, the system tends to become unstable and falls into chaos. On the basis of the parameters assigned above, we use numerical simulation to describe the dynamic behaviors of system (8).

4.2. The Effect ofPrice Adjustment Speed on the System

(1) The Influence Which the Price Adjustment Speed Has on the Behaviors of the Two Manufacturers and the Retailer. Since the two manufacturers' behaviors are similar, we only discuss the influence on system behaviors when the parameter of M1 is changed.

First, we can get the price trajectory diagrams of the two manufacturers and retailer with change of k1 when

Figure 7: Power spectrum of w1 ,w2,p1, and p2 when k1 = k2 = 0.0005.

k2 = 0.001, w1 = 20, and w2 = 15, as shown in Figure 1(a) and Figure 1(b). We can obtain that w2 is less affected by change of k1 and w1,w2, p1, and p2 change from the stable period, period-doubling bifurcation to the chaos in three trajectories. When k1 e [0, 0.00102], w1, w2,p1, and p2 are stable. When k1 = 0.00102, the first bifurcation appears in w1,w2,p1, and p2, and the price of the two manufacturers and retailer vibrates in two points; after that the second bifurcation appears in the system; finally, the system goes into chaos, the price behaves more disorderly, and the market behaviors become unpredictable. Figure 1(c) shows corresponding change of the Lyapunov exponent. The positive Lyapunov exponent is used to mark the chaos; the bigger the Lyapunov exponent, the stronger the chaos. The system is in chaos when most of the Lyapunov exponents are positive. Figure 2 shows the price trajectory diagrams of the two manufacturers and retailer with the change of k1 and k2 simultaneously. We can observe that the stable regions of w1, w2 are smaller than the one with change of only k1 or only k2 and the w1,p1, and p2 change from the stable period and period-doubling bifurcation to the chaos in four trajectories; the dynamic characteristics of w2 in particular do not follow the period-doubling bifurcation. Because k1 and k2 change at the same time, the retailer's behaviors appear in complicated characteristics.

Figure 3(a) shows an attractor of the supply chain which is in stable state when k1 = k2 = 0.001. Figure 3(b) shows a chaos attractor of the supply chain which appears in a complex state when k1 = k2 = 0.0015; it is another chaos characteristic of the variables.

Proposition 1. The value of price adjustment speed determines whether the supply chain is stable or not; when we determine this value, each manufacturer should consider market reaction of the competitors and retailer. Only upstream and downstream enterprises keep stable, to ensure the stability of the supply chain and to maximize the enterprise's profit.

(2) The Influence Which Change of Risk Preference Has on the Behaviors of the Two Manufacturers and the Retailer. Let k2 = 0.001, Rr = 120; the value of other parameters is the same as in the previous assumption; we can obtain the price bifurcation diagrams of the two manufacturers and the retailer with change of k1, as shown in Figure 4. We can see that the stable region is not changed, wholesale price of the two manufacturers declines, and the retail price goes up with change of Rr. Figure 5 shows the price bifurcation diagrams of the two manufacturers and the retailer with change of k1 when k2 = 0.001 ,RMi = 120, the first bifurcation point is k1 = 0.00129, and the equilibrium value is (46.73, 26.4, 48.89, and 40.6). We can make a conclusion that the risk preference of the manufacturer can affect the stable region of system and change the equilibrium point of system, and the risk preference of the retailer cannot affect the stable region of system and change the equilibrium point of system.

(3) The Influence Which the Uncertain Demand Has on the System Variable Behaviors. Because the demand forecasting exists errors, the customer demand is always uncertain. Next we will observe the price change of the two manufacturers and retailer with change of k1 when a = 70, k2 = 0.001, and it is shown in Figure 6. We can see that, with a increasing, the first bifurcation point of the system is k1 = 0.00105, and w1,w2,p1, and p2 are all falling.

Proposition 2. Uncertain demand makes the supply chain access chaos quickly, and the wholesale price and the retail price decline. It is uncertain about their respective profit; thereby the competitiveness of the supply chain becomes weak.

(4) The Power Spectrum of Variables and the Sensitive Dependence on Initial Conditions. We adopt a cycle diagram method to estimate the power spectrum of variables. Next, we will observe the price change of the two manufacturers and retailer when k1 = k2 = 0.001 and the initial values are w1 = 20, w2 = 15; it is shown in Figure 7. We know that the supply chain is stable when k1 = k2 = 0.001, so the power spectrum of the variables is straight lines which conform to the attractor in Figure 3(a). When k1 = k2 = 0.0015, the supply chain is in chaos. From Figure 8, we can see that w1,w2,p1, and p2 vibrate with a frequency, w2 changes with approximate periodic motion, and range of movement of p1 is bigger than the one of p2.

Then, we will observe the price change of the two manufacturers and the retailer when k1 = k2 = 0.0015 and (w1 ,w2) = (20.01, 15) and (20, 15) which has smaller change in w1 and no change in o^they are shown in Figure 9. We can see that the price of the two manufacturers and retailer has distinct change. The sensitive dependence on initial conditions is another important characteristic of the chaotic system as it fully manifests the sensitive dependence on the initial conditions of the system (8).

Through the above analysis, we know that the risk-averse supply chain has been in chaos. At this state, the price of the two manufacturers and the retailer changes disorderly from the beginning of the Nash equilibrium.

(5) The Effect of Parameter Change on the Profit. Figure 10 shows the profit bifurcation of the two manufacturers and retailer with changes of k1. Obviously, the profit bifurcation is similar to price bifurcation including period-doubling bifurcation, four-period-doubling bifurcation, and chaos state. Figure 11 shows the profit of the two manufacturers and retailer in 50 games when the system is in a stationary period, two-period-doubling bifurcation, and a chaotic period. In the different period, the fluctuation range of price of M1 is larger than that of M2.Tables 1,2,and 3 give the profit data of the two manufacturers and retailer, respectively, in different periods with change of k1, RM^, and a.

Proposition 3. First, the total profit of the system in the chaotic period is less than that in other periods. Second, the total profit of the system and respective profit of the two manufacturers and the retailer will increase with increase ofRM^ and decrease with increase of a. Third, when the system is in chaos, the profit of the two manufacturers will decrease, but the profit of retailer will increase. Namely, when the system goes into chaos, the retailer is a beneficiary. This is why some participants set out contract to

Table 1: The profit of the two manufacturers and retailer in different periods when k1 = 0.001.

Different periods Stable period Two-period-doubling bifurcation Chaotic period

nR 8159.76 8144.99 8175.85

nMl 8074.74 6126.11 2403.59

nM2 8128.54 8149.99 8054.7

In 24363.04 22421.09 18634.14

avoid disordered competition in some situations and why some participants prefer chaotic market in some cases.

5. Chaos Control

Competitive manufacturers will certainly want to achieve maximum profit in the existing supply chain. Through the above analysis, we can see that the change of k1,k2, and a often causes disorder behaviors in the market which are

40 30 20 10

3 0 -10 -20 -30 -40

20 40 60 80

20 40 60 80 100

80 100

10 8 6 4 2 0 -2 -4 -6 -8 -10

80 100

Figure 9: Different change of price of the two manufacturers and the retailer with increase of game at k1 = k2 = 0.0015.

Table 2: The profit of the two manufacturers and retailer in different periods when RMi = 120.

Different periods Stable period Two-period-doubling bifurcation Chaotic period

nR 8182.74 8190.65 8203.86

15539.95 14331.78 12443.96

nM2 8144.84 8136.82 8125.35

In 31867.63 30659.25 28773.17

Table 3: The profit of the two manufacturers and retailer in different periods when a = 120.

Different periods Stable period Two period-doubling bifurcation Chaotic period

nR 6001.45 6009.92 6013.69

nMl 5918.76 4346.43 3573.6

nM2 5972.61 5968.58 5967.1

In 17892.81 16324.92 15554.38

of disadvantage to the stability of the supply chain and the development of the enterprise. However, the participants often maximize their own profit by any kind of means in the process of marketization. So the market will be out of order and finally falls into chaos. It is particularly important that each participant should make rational strategic decision timely for making the system return to the stable equilibrium.

Parameter adjustment and feedback control method of the system variable will be used to control the chaos of

the system (8). It is often used at the chaos control of general discrete dynamic system. We will analyze system chaos control based on the influence of k1 on the stability of the risk-averse supply chain. Assume system (8) is wt (t + 1) = fi(w1 (t), w2(t)); the system under control is as follows:

w1 (t + m) = (1 -u) f™ (w1 (t), w2 (t)) + uw1 (t), w2 (t + m) = (1 -u) f.2 (w1 (t), w2 (t)) + uw1 (t),

9000 8000 7000 6000 5000 4000 3000 2000 1000 0

X: 0.00104

Y: 8072

8220 8200 8180 8160 8140 8120 8100 8080 8060 8040 8020

Figure 10: The profit bifurcation with change of kl when k2 = 0.001.

Pi (t + m) = f3 (w1 (t), w2 (t)), Pi (t + m) = f4 (wi (t), w2 (t)).

When u = 0, the controlled system is the original system; thus, the controlled system has the same periodic orbit a the original system. When m = 1, each step iteration of the fixed point is in control. When m = 2, 4, and so forth, namely, to control the second cycle orbit, four cycles orbit. As long as we select an appropriate value for u, we can ensure delayed bifurcation with ki at fixed point and keep the supply chain being stable within large scope of ki. In addition, with appropriate adjustment value of m, we can realize stability control of higher cycle orbits for the chaotic attractor. Then, we discuss the stability control in the supply chain when m = 1 and try to stabilize the price on the fixed point. The controlled system can be represented as

( dE(UM )\ wi (t+1) = (1-u)(wi (t) + kiwi (t) —^-— ) + uwi (t),

w2 (t + 1) = (1 - u) \w2 (t) + k2w2 (t)

9E(uMi )

+ uw2 (t) ,

Pi (t)

a + biwi (t) - d2w2 (t)

(wi (t) + W2 (t))a2 R„

a + b2w2 (t) - diwi (t)

(wi (t) + W2 (t))g2 R„

— + 2bi ) ( — + 2b2

+ a-, -

— + 2b0 R»

+ do -

8200 8000 7800 7600 7400 2 7200

7000 6800 6600 6400 6200

8500 8000 7500 7000

« 6500

^ 6000 5500 5000 4500 4000

aaaaaaaaaaaaaaaaaaaaaaa^

I AaaJLAJLA/ [IIIIIIIIIIIIII

9000 8000 7000 6000

£ 5000

^ 4000 3000 2000 1000 0

Figure 11: The profit of the two manufacturers and retailer k2 = 0.001: (a) stable period; (b) two-period-doubling bifurcation; (c) chaotic

p2 (f)

a + b1w1 (t) - d2w2 (t)

x\ ( — + 2h)( — + 2b2

Rr 1 )\ R 2

(wi (t) + W2 (t)) a:

a + b2w2 (t) - d1w1 (t)

(wi (t) + W2 (t))a2 ' Ra

h +a2 -

+ do -

^ + 2b1 R*

Figures 12 and 13 show that the chaos system can be controlled gradually from the four-period-doubling bifurcation and two-period-doubling bifurcation to the fixed point with the

38 36 34 32 30 28 26 24 22 20

kl (a)

32 31 30 29 28 27 26 25 24 23

ki (b)

Figure 12: Control bifurcation diagram of w1 and w2 with change of k1 when k2 = 0.001: (a) u = 0.5; (b) u = 0.7.

50 45 40 35 p 30 25 20 15 10

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

(a) (b)

Figure 13: Control bifurcation diagram of w1, w2, p1 and p2 with change of u when k1 = 0.002, k2 = 0.001: (a) w1 and w2; (b) p1 and p.

45 40 35 30 25 20 15 10 5 0

control parameter u increasing. When u > 0.4825, the controlled system can be stabilized at the Nash equilibrium point. In a real market, we can consider u as the regulation on the price adjustment speed which could avoid market chaos when manufacturers pursue their maximum profit. We can also consider u as the learning ability or adaptability of the market. For instance, the manufacturer will adjust price according to the information in the past. Due to the complexity, disorder, and randomness in chaotic state, the manufacturer should have a clear forecast of chaos control. After adjusting parameter u, we can make the periodic orbit stable at the expected point.

6. Conclusions

Considering the randomness of customer's requirements and supply's risk preference, we study the supply chain in the market which consists of two manufacturers and a retailer. We

construct a dynamic Stackelberg and Bertrand pricing game model and find that bifurcation, chaos, and other complex phenomena occur when the price adjustment speed, the level of risk preference, and predict error change. When the chaos occurs, the stability of the whole supply system is broken, and the market becomes abnormal, irregular, and unpredictable. It is important to note that, when determining the value of the parameter, each manufacturer should consider the market reaction of competitors and retailer to ensure that the upstream and downstream enterprises keep stable and ensure that the supply chain keeps stable for their maximum profit. Finally, we use the parameter adjustment method to control the supply system. Then, we have obtained some conclusions from our research.

(1) The two manufacturers play the Bertrand game; when M1 changes its price adjustment speed, its wholesale price presents deterministic characteristics of chaos,

and the wholesale price of M2 is influenced less. The price of retailer is influenced largely. When price adjustment speed of the two manufacturers changes at the same time, it has no effect on the optimal pricing strategy, the stable region of the supply chain gets smaller, and the two manufacturers all show complicated behavior characteristics.

(2) Increasing the risk tolerance level of the retailer does not impact stable region, but retail price will increase and wholesale price of the two manufacturers will fall. If M1 increases his risk tolerance level, the stability between M1 and the retailer is enhanced, and the two manufacturers' wholesale price and products' retail price go up. Uncertain demand makes the wholesale prices and retailer prices access chaos quickly, two manufacturers' wholesale price and retail price decline, and total profit and their respective profit decrease.

(3) The total profit of the system in the chaotic period is less than that in other periods; it and respective profit of the two manufacturers and the retailer will increase with increase of and decrease with increase of a. When the system is in chaos, the profit of the two manufacturers will decrease, but the profit of retailer will increase. Namely, when the system goes into the chaos, the weaker (retailer) is a beneficiary.

This paper is a realistic guide for the risk-averse supply chain to formulate its parameter adjustment strategies to avoid the loss of the respective profit and the total profit. It is also a realistic reference for the managers to formulate relevant policies on macroeconomic control. The managers can adjust parameters to make the supply chain in good operating condition based on the operation condition of the supply chain; for example, when the variance of customer demand gets bigger, the manager should reduce the price adjustment speed for the stability of the supply chain.

There are several possible directions for the future study. First, one can study three-echelon supply chain with complexity behaviors under the different game strategies considering the delay time, such as the Stackelberg game between the two manufacturers. Second, one can adopt a different form of demand function and, finally, other control methods maybe applied in order to achieve different results.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank the reviewers for their careful reading and providing some pertinent suggestions. The research was supported by the National Natural Science Foundation of China (No. 61273231), Doctoral Fund of Ministry of Education of China (Grant No. 20130032110073) and supported by Tianjin University Innovation Fund.

References

[1] J. Wei, J. Zhao, and Y. Li, "Pricing decisions for complementary products with firms' different market powers," European Journal of Operational Research, vol. 224, no. 3, pp. 507-519, 2013.

[2] S. K. Mukhopadhyay, X. Yue, and X. Zhu, "A Stackelberg model of pricing of complementary goods under information asymmetry," International Journal of Production Economics, vol. 134, no. 2, pp. 424-433, 2011.

[3] O. Caliskan-Demirag, Y. Chen, and J. Li, "Customer and retailer rebates under risk aversion," International Journal of Production Economics, vol. 133, no. 2, pp. 736-750, 2011.

[4] C. Luo and J. Huang, "Pricing and ordering strategies in supply chain with different risk preference," Computer Integrated Manufacturing Systems, vol. 18, no. 4, pp. 859-866, 2012.

[5] B. Hu and X. Chen, "Supply chain sale-surety option contract with effort dependent demand and risk preference," Journal of Industrial Engineering: Engineering Management, vol. 26, no. 3, pp. 184-190, 2012.

[6] Z. Guan and B. Zhou, "Integrated optimization of three-level supply chain under different risk attitude," Operations Research and Management Science, vol. 20, no. 6, pp. 19-24, 2011.

[7] S. Huang and C. Yang, "Supply chain coordination model considering risk attitudes of agents with VaR constraints," Journal of Industrial Engineering: Engineering Management, vol. 25, no. 2, pp. 136-141, 2012.

[8] T. Puu, "The chaotic duopolists revisited," Journal of Economic Behavior and Organization, vol. 33, no. 3-4, pp. 385-394, 1998.

[9] J. Ma and Z. Sun, "The research on price game model and its complex characteristics oftriopoly in different decision-making rule," Nonlinear Dynamics, vol. 71, no. 1-2, pp. 35-53, 2013.

[10] J. Ma and H. I. Bangura, "Complexity analysis research of financial and economic system under the condition of three parameters' change circumstances," Nonlinear Dynamics, vol. 70, no. 4, pp. 2313-2326, 2012.

[11] H. B. Hwarng and N. Xie, "Understanding supply chain dynamics: a chaos perspective," European Journal of Operational Research, vol. 184, no. 3, pp. 1163-1178, 2008.

[12] J. Wang and X. Wang, "Complex dynamic behaviors of constrained supply chain systems," Systems Engineering-Theory & Practice, vol. 32, no. 4, pp. 746-751, 2012.

[13] J. Ma and Y. Feng, "The study of the chaotic behavior in retailer's demand model," Discrete Dynamics in Nature and Society, vol. 2008, Article ID 792031,12 pages, 2008.

[14] E. Gal-Or, "Information sharing in oligopoly," Econometrica, vol. 53, no. 2, pp. 329-343,1985.

[15] J. S. Raju and A. Roy, "Market information and firm performance," Management Science, vol. 46, no. 8, pp. 1075-1084, 2000.

[16] X. Vives, "Duopoly information equilibrium: cournot and Bertrand," Journal of Economic Theory, vol. 34, no. 1, pp. 71-94, 1984.

[17] A. Kunstman, "Decision making under uncertainty: models and choices," European Journal of Operational Research, vol. 8, no. 4, pp. 401-402,1981.

[18] M. R. Walls, "Combining decision analysis and portfolio management to improve project selection in the exploration and production firm," Journal of Petroleum Science and Engineering, vol. 44, no. 1-2, pp. 55-65, 2004.

[19] R. D. Banker, I. Khosla, and K. K. Sinha, "Quality and competition," Management Science, vol. 44, no. 9, pp. 1179-1192,1998.

Copyright of Discrete Dynamics in Nature & Society is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.