Complexity
Volume 2018 (2018), Article ID 7812784, 13 pages
https://doi.org/10.1155/2018/7812784

Research Article

Research on Supply Chain Stability Driven by Consumer’s Channel Preference Based on Complexity Theory

Yi Tian, Junhai Ma, and Wandong Lou

College of Management and Economics, Tianjin University, Tianjin 300072, China

Correspondence should be addressed to Junhai Ma; mjhtju@aliyun.com and Wandong Lou; wdl160211@126.com

Received 22 January 2018; Accepted 3 May 2018; Published 7 August 2018

Academic Editor: Borja Sanz-Altamira

Copyright © 2018 Yi Tian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

At present, most of the manufacturers are increasingly laying emphasis on dual-channel selling. The level of consumer’s preference is diverse with different channels. This study discusses a supply chain where one common manufacturer uses two different channels to sell product. The bounded rationality expectation has been applied to explore the decision-making mechanism in finalizing the wholesale prices and sales commissions for the manufacturer and retailers. In addition, the dynamic features of the system are simulated by 2D bifurcation diagram, the largest Lyapunov exponent, attractor variation, and time series. The simulation results suggest that if the adjustment speeds of the wholesale prices and sales commissions change drastically, the system would fall into a chaotic state. If the consumers prefer the online channel, the wholesale prices and sales commissions may be raised and vice versa. However, too strong preference of each channel will make the system appear to have a considerably periodic fluctuation, even chaos. In the meantime, the parameter adjustment method can help make the periodic or chaotic system go back to stability. Therefore, the significance of this study is with a practical meaning to make better pricing strategy for players in dual-channel selling supply chain system.

1. Introduction

The rapid development of Internet technology and social economy promoted a huge amount of manufacturers who had been transferred from traditional retailing to the combination of traditional and electronic retailing. For instance, according to the data provided by the China Academy of Information and Communications Technology (CAICT), sales of smart phone in China records nearly 390 million at year-end of 2016. Huawei took a 16.4 percent market share in China, while Apple, 9.6 percent. Smart phone manufacturers always distribute their products through both traditional and online channels. In this research, a pricing problem in supply chain with one manufacturer (Huawei) who supplies an identical product (Mate 9) to two retailers (China Resources Vanguard Co. Ltd. and JD.com) is brought to discussion. CR Vanguard is one of the world’s top five hundred enterprises, and it is China’s largest supermarket chains by a number of stores for years. JD.com is a typical Chinese e-commerce company which is one of the biggest e-malls in China with self-support storage and logistics system. CR Vanguard and JD.com are the typical representatives of traditional and online channels for smartphone sales, respectively. Different consumers will choose their own preferred buying channels. Therefore, the supply chain driven by the consumer’s channel preference is necessary to be studied.

Two-echelon supply chain has been widely paid attention to and discussed in academia. Ma, Ma and Bao [1, 2] investigated the bullwhip effect in a two-echelon supply chain. This study discusses a supply chain which consists of one common manufacturer and two retailers from traditional and online sales channels, respectively. The manufacturer determines the wholesale prices and the two retailers decide their own sales prices.

In the current studies [3–7], the consumer preference is frequently considered as one of the chief factors that influence the marketing. Especially, consumer’s channel preference is a hot topic in supply chain research. Khouja et al. [8] analyzed the manufacturer’s channel selection and pricing problem. They considered consumer’s preferences under different distribution strategies. Gao and Su [9] investigated a stylized model where a retailer operated the dual-channel. Customers strategically make channel choices that included both online and offline. Ke and Liu [10] discussed a dual-channel supply chain competition with channel preference under an uncertain environment. Their results indicated that the customer’s direct-channel preference had the opposite effect on the profit of supplier and retailer. Gan et al. [11] presented pricing decisions and profits of supply chain members in a closed-loop supply chain and considered the parameter of customer’s direct channel preference as a key factor.

Different from the above mentioned literatures, this paper discusses the decision-making of wholesale prices and sales commissions for manufacturer and retailers based on the bounded rationality expectation by employing complexity theory to explore the effect of the adjustment speed of the above variables on the system stability. More and more researchers used bounded rational theory to analyze the economic model. Ma and Wang [12] applied two noncooperative game models to study the bounded rationality expectation. Song and Zhao [13] pointed out a decision-making problem of newsvendor system with strategic customers of bounded rationality. Ma and Si [14] established a continuous Bertrand duopoly game model with two-stage delay to investigate the influence of delay and weight on the complex dynamic characteristics of the system. Ma and Xie [15] analyzed the comparison and complexity on a dual-channel supply chain under the different channel power structures and uncertain demand. Ma and Xie [16] discussed the impact of loss sensitivity on a mobile phone supply chain system stability based on the chaos theory. In conclusion, it has received a considerable amount of attention to apply bounded rationality and complexity theory to study the management problems in the supply chain.

In reality, the exploration of the supply chain driven by consumer’s channel preference based on chaos theory has a great significance. In this paper, besides the adjusted basic parameters such as the wholesale prices and sales commissions, the consumer’s channel preference is also discussed. The results indicated that the preference plays an important role in status transformation from stability to chaos. In this analysis, a supply chain model which depicts one common manufacturer using two different channels to sell its product is built in the first place. Then, the bounded rationality expectation is applied to study the complex dynamics of the system. Furthermore, the dynamic characteristics of the system are discussed by 2D bifurcation diagram, the largest Lyapunov exponent, attractor variation, and time series. The results indicate that not only the drastic adjustment of traditional parameters such as the wholesale prices and sales commissions can make the system go chaotic but also the consumer’s channel preference can also lead to system considerably periodic fluctuation, even to chaos, which is the primary finding in this article. Finally, the parameter adjustment control method is performed on the system to return the periodic or the chaotic state to stability. The conclusions of this paper can give constructive suggestions for managers of both manufacturers and retailers.

The structure of this paper is as follows. In Section 2, the basic model considering the influence of consumer’s channel preference is established. In this part, the optimal decision of the complex dynamic system is analyzed. Following Section 2, this study further applies numerical simulation to investigate the stability aspect and the dynamical behaviors. 2D bifurcation diagram, the largest Lyapunov exponent, attractor, and time series have been utilized to study the influence on the supply chain. In Section 4, we introduce the method to control the chaos system. Finally, the conclusion of this research is provided.

2. Demand Model

A supply chain system consists of one manufacturer and two retailers, one of which is a traditional retailer (retailer 1) and the other is an online retailer (retailer 2), which has been restricted in this study. The manufacturer provides an identical product at a constant unit cost () and wholesale prices () to retailer 1 and retailer 2. Then, the retailers distribute the product to consumers in a common market by traditional and online channels, respectively. The retailers decide their unit sales commissions . Retailer 1 and retailer 2 have a competitive relationship between each other as shown in Figure 1.

Figure 1: The model of supply chain.

Hence, the sales prices of this supply chain can be described as . The demand function is simplified as the following linear form:

This simplified linear demand function has been used by Anderson and Bao [17], Wu et al. [18], Yang and Zhou [19], and Choi [20]. Here, is the market base of channel , is the sales price of retailer , is the price sensitivity of the product, and denotes the substitutability of the retailers. In order to discuss the influence of consumer’s channel preference and pricing strategies on the performance of the supply chain members without loss of generality, the method proposed by Mcguire and Staelin [21] has been applied to rescale the demand functions in (1). The scaled model can conduct some analytical comparisons instead of numerical experiments reported in Xiao et al. [22]. Let ; the above linear function will be changed as

The substitutability of the two retailers in terms of changes in prices has been denoted as parameter . Let , and , and define . The parameter depicts the consumer’s channel preference.

represents the scenario in which the submarket is monopolized by one retailer. represents the scenario in which there exist two completely substitutable retailers. Then, the prices can be standardized as , , and .

The demand function can be formulated as follows:

can be obtained by using the same method:

So we can get the demands as below:

But in reality, the market requirements are often influenced by some other factors; here, is used to represent these disturbance factors in research.

follows a normal distribution such that , . The expected profit decision models of the manufacturer and two retailers are shown as

In order to deduce the marginal profits of the manufacturer and two retailers, the first-order partial derivative of (7) can be calculated as follows:

2.1. Optimal Decision

In the market, the optimal decisions of the retailers and manufacturer are influenced by each other. Therefore, this study assumes that they share overlapping information. In this research, when the marginal profits reduce to zero, the optimal decision for retailers and manufacturer can be obtained. The wholesale prices and sales commissions are the decision variables. The research calculates the first-order derivative of the retailers’ and manufacturer’s profit. When the marginal profit equals to zero, the optimal wholesale prices and sales commissions can be derived. The detailed process is as follows:

The optimal solution of the system will be solved.

2.2. The Complex Dynamics of the System

The real market is a changeable and complex system. The manufacturer and retailers cannot know the complete market information, and predict other decision-makers’ information accurately. Hence, the bounded rationality expectations have been universally adopted by market decision-makers. The decision of period is as follows:

In the paper, the price adjustment speeds consist of two wholesale prices and sales commissions. Here, ; and are the adjustment speed of the sales commissions with retailer 1 and retailer 2, respectively. and represent the adjustment speed of the wholesale prices with the manufacturer. Based on the above equation, the four-dimensional discrete dynamic system can be acquired and expressed in (12). The wholesale prices and the sales commissions in the next period depend on the estimation of these period marginal profits. Base on that, the adjustment of bounded rationality is completed. If the entity’s current marginal profit is positive, the price will increase during the next period for higher profits. If the entity’s current marginal profit is negative, the price should be reduced during the next period. The dynamic adjustment processes are as follows:

After a long period of operation, when and , the system will become steady. Under the equilibrium state, all parties in the game cannot increase profits for themselves by changing the decision variables. Therefore, when the decision variables take the values of the equilibrium state, every player in the game can obtain the maximized profit. There are sixteen equilibriums of system (12):

For these above equilibriums, it can be found that is the only point which is completely nonzero. The others are boundary points. is the unique Nash equilibrium point. The Jacobian matrix of system (12) can be calculated as follows:

The Jacobian matrix and Jury criterion have been applied to discuss the local stability of the Nash equilibrium in this paper.

The characteristic polynomial of the Jacobian matrix is as follows:

In order to guarantee the local stability of the Nash equilibrium, must satisfy the following conditions:

The details of , and are as below:

Because the conditions are very complex, the process to solve (12) is a tedious task. If the Nash equilibrium satisfies (12), it can be considered as locally stable. In the next section, we will utilize numerical simulation to further analyze the dynamic characteristic on the system.

3. Dynamic Characteristics of the System

Recently, the numerical simulation method has been widely applied to represent complex dynamic system in literatures. In this study, is the sales commissions’ adjustment speed of retailer . represents the adjustment speed of the wholesale prices () of the manufacturer. The values of the system parameters are defined as , , , , , , , , , and . The stable region, bifurcation, and chaos will be studied as follows.

3.1. The Stable Region

The decision variables of the system are wholesale prices and sales commissions. The system will change to be either stable or unstable, if the wholesale prices and sales commissions obtain different values. The stability region of the deterministic system is shown as the blue and green part in Figures 2 and 3, respectively. The white areas are the unstable range. If the applied values come from the blue and green region of Figures 2 and 3, after the iteration of period by period, the system will tend to reach the equilibrium state. If the values landed in the white areas, the system will exhibit unstable phenomena.

Figure 2: The 2D stable region of and .

Figure 3: The 2D stable region of and .

In Figure 2, when increases, the tendency of should be reduced. In the stable region, it has the opposite effect on and , so does Figure 3. Through Figures 2 and 3, we found that the stable range of and is almost twice as much as and , when the value of is larger than others’ price adjustment parameters. Period doubling and chaos should be avoided. The supply chain system should be kept in the stable state as long as possible.

3.2. The Bifurcation and Chaos Behavior

The parameter basin has been used to describe 2D bifurcation effect. The results are shown in Figures 4 and 5. Different colors represent different cycles. Blue represents the stable region in Figure 4, and green represents the stable region in Figure 5. Orange, pink, purple, indigo, brown, yellow, and black regions represent the stable cycles of the periods from 2 to 8, respectively. The gray region is non-convergence, and the white region is divergence in both Figures 4 and 5. The system can remain stable, when the adjustment parameter is relatively small. As and increase, the system will go into the stable cycles of period 2. If and continue to increase, the system will go into the stable cycles of period 4, which is shown in Figure 4. However, the system continually changes to the stable cycles of periods 2, 4, 8, and so on. If and keep on increasing, the system will go into chaos state. This process is depicted in Figure 5. In this situation, the market will become chaotic and the market participants cannot make good decisions to reach the best profits.

Figure 4: 2D bifurcation diagram of and .

Figure 5: 2D bifurcation diagram of and .

When the adjustment parameter changes in a small range from 0, the variables will be stable. However, if the decision variables adjust sharply, the variables will be unstable. Especially, in chaotic state, the system will become sensitive and unpredictable. In economic terms, when the adjustment parameter exceeds a certain limit, the intense competition and price drastic fluctuation will occur.

The bifurcation diagram is an intuitive way to show the dynamic characteristics of the system. One of the parameters is fixed and the other changed within a certain range. In the following bifurcation diagrams, the blue line represents and the red line . The green line represents and the purple line .

The largest Lyapunov exponent (LLE) is another efficient way to make clear the complex dynamic characteristics of the system, which also is corresponding to the variation of the bifurcation diagram when discussing the same parameter. The dramatic variation of LLE is accompanied by the decision variables changing. If the LLE < 0, it means the system is in a stable state. When the LLE = 0 for the first time, the system stays at a critical point and usually the bifurcation appears. If the LLE > 0, it indicates that the system runs into the unstable state. The manufacturer and retailers should be cautious to prevent system chaos.

The following part gives the bifurcation diagram and LLE diagram with respect to the parameters , , , and . The bifurcation diagrams of , , , and , with respect to the adjustment speeds of the sales commissions, are shown in Figures 6 and 7. In the same condition, Figures 8 and 9 depict the variation of LLE with respect to and .

Figure 6: Bifurcation diagram of , , , when .

Figure 7: Bifurcation diagram of , , , when .

Figure 8: The largest Lyapunov exponent of the system with respect to .

Figure 9: The largest Lyapunov exponent of the system with respect to .

When and with the other parameters settled down, the conditions of , , , and are shown in Figure 6. Because the dots of and gather together, we further enlarge critical parts in sub-windows in order to describe the bifurcation conditions of and . The upper small picture is the blowup of the bifurcation diagram for and , and the other one is the blowup of the bifurcation diagram for . When , , , , and are all in the stable state. When grows to 2.339, , , , and bifurcate at the first time and enter into the stable cycles of period 2. If continues to increase, the system will enter into the stable cycles of period 4.

Figure 8 reveals the variation of LLE along with increasing, which expresses the pattern of the complex dynamic behavior caused by the growth of the price adjustment speed. Figure 8 explores the LLE of the system with respect to . When , the LLE is below zero, which indicates that the system is stable. When increases to 2.339, the LLE is equal to zero and the bifurcation occurs in the system. The bifurcation point in Figure 8 is the same as that in Figure 6.

Figures 7 and 9 also have the same bifurcation point. , , , and bifurcate at 2.156 as varies from 0 to 2.98. If continues increasing, the system runs into the cycles of period 4. In order to describe the changing processes of the system from the stable state to the chaotic region more clearly, and start from 1.5 in Figures 6–9.

The bifurcation diagrams of , , , and , with respect to the adjustment speeds of the wholesale prices, are shown in Figures 10 and 11. In the same condition, Figures 12 and 13 depict the variation of LLE with respect to and .

Figure 10: Bifurcation diagram of , , , when .

Figure 11: Bifurcation diagram of , , , when .

Figure 12: The largest Lyapunov exponent of the system with respect to .

Figure 13: The largest Lyapunov exponent of the system with respect to .

, , , and bifurcate at 1.096 as varies from 0 to 1.70. , , , and bifurcate at 1.107 as varies from 0 to 1.73. If and continue increasing, the system runs into cycles of periods 4, 8, and so on. Next, the system will fall into chaos gradually. In this part, and start from 0.7 in Figures 10–13.

Above all, Figures 10–13 sufficiently indicate that if the wholesale prices of the manufacturer are adjusted drastically, the system will fall into chaos. On the other side, as shown in Figures 6–9, if the sales commissions of the retailers are adjusted radically, the bifurcation will appear but the chaos will not occur. From an economic point of view, if the adjustment parameter rapidly changes, the market will fall into the chaotic state. Therefore, the manufacturer and retailers should keep the price adjustments within a reasonable range.

The substitutability of the retailers will be discussed in the following part. represents the demands which are absolutely independent. Each retailer is a monopolist in its own area. represents the two retailers who are totally substitutable. When and the other parameters are kept constant, the situations of , , , and are shown in Figure 14. When grows to 0.872, the system begins the first bifurcation which is the stable cycles of period 2. If continues to increase, the system runs into the stable cycles of periods 4, 8, and so on. Then, the system will fall into chaos gradually. The small picture in Figure 14 is the blowup of the bifurcation region for , , , and , when .

Figure 14: Bifurcation diagram of , , , when .

Figure 15 explores the LLE of the system with respect to . When , the LLE is less than zero, and the system is in a stable state. While grows to 0.872, LLE changes into zero for the first time, and the system falls into the bifurcation state, the completely same bifurcation points in Figures 15 and 14.

Figure 15: The largest Lyapunov exponent of the system with respect to.

In Figure 14, the manufacturer is more sensitive to the substitutability of the retailer, which leads to a rise in the wholesale price. However, the retailer’s profits are robust to the changes in the retailer’s substitutability. Under this condition, fierce competition may lead to the dominance of the manufacturer. Nevertheless, bifurcation and chaos will still occur in the market after the fierce competition, once a certain extent of substitutability remains among the retailers of different channels.

The impacts of consumer’s channel preference on the dynamic system will be discussed in the following part. When , , , , , and , the situations of , , , and are shown in Figure 16. The lines in the bifurcation diagram indicate that other parameters are stabilized separately, when is changing in the range of [0, 3.80]. When grows to 2.558, the system begins the first bifurcation which is the stable cycles of period 2. If continues to increase, the system runs into the stable cycles of periods 4, 8, and so on. After that, the system will fall into chaos gradually. Figure 17 proves the same conclusion as Figure 16 does.

Figure 16: Bifurcation diagram of , , , when .

Figure 17: The largest Lyapunov exponent of the system with respect to .

The curves of the traditional retailers and online retailers are shown in Figure 16. When , it can be found that the channel preference of consumers will convert from one channel to another. If one channel has more consumer preference, its wholesale prices and sales commission will rise; and are displayed in Figure 16. However, if the channel preference obviously comes to one channel, it is possible that bifurcation and chaos appear in the system.

In this paper, the chaos attractor has been investigated to the dynamic characteristics of the system. Strange attractor is a set of numbers of indefinite points. The larger the space of the strange attractor is, the more complex the structure will be. No matter how instable the system is, the points are always attracted along a certain trajectory.

Figure 18 shows the chaotic attractor of and , when , , , , , , , , , and , the system becomes complex and even moves into chaotic state. The dynamic characteristics of the system are obvious. Based on parameter studies, changes in a larger range than when the adjustment speeds of wholesale price change more intensively than .

Figure 18: The chaotic attractor of and .

The time series of different decision variables are shown in Figures 19 and 20. In the chaotic state, the system is more sensitive to the initial value. The phenomenon can be called as the initial value sensitivity. In Figure 19, to better show the character of the chaotic system, we assume the retailer 1 to make decisions under the condition that and show the time series of after 100 iterations. The blue line represents the initial value with , and the red line represents the initial value with , when the dynamic system state is in a chaotic state. We also set the other parameters’ values as follows: , , , , , , , and . As shown in Figure 19, although the initial value of sales commissions only changes from 0.1 to 0.101, after 100 iterations, a large gap between the two tracks will occur.

Figure 19: Time series of in a chaotic state.

Figure 20: Time series of , , , in a stable state.

Figure 20 shows the time series of , , , and in a stable state. We set the initial values of the wholesale prices and sales commissions as , , , and in Figure 20. We set the initial value of and from a low value, while we set the initial value of and from a high value. We set the same initial values of , , , and with those above. , , , and would gradually flatten after 12 iterations shown in Figure 20.

In short, Figures 19 and 20 have different meanings. The adjustment speeds of the wholesale prices and sales commissions depend on the price fluctuation levels. As such, the manufacturer and retailers should keep price adjustment in reasonable ranges.

4. Chaos Control

The excessive competition will cause the chaotic phenomenon of a market, which poises negative effect on all competitive subjects in the market and impedes the economic operation. Therefore, some chaos control methods have been made use of to propel the system to return to stable equilibrium. There are a lot of chaos control methods such as continuous feedback control, self-adaptation control, and OGY. In this research, parameter adjustment control method has been used. The manager can control the chaotic behavior through introducing the control parameter . Huang [23] conducted the chaos control with setting the control parameter as . The values of the other parameters have been set in the sections above. The stable range and the bifurcation diagram will be given to account for the state of the system. Under the parameter adjustment control, the system can be expressed as follows:

Under the process of chaos control, three parameters, including wholesale price, consumer’s channel preference, and the substitutability of retailers, will be discussed in this part. If equals 1.65, the system will be chaotic. When the control parameter varies from 0 to 0.7, the changes of the system are shown in Figure 21. Besides, if equals 0.892, the system will be chaotic. When the control parameter varies from 0 to 0.3, the changes of the system are shown in Figure 22. Finally, if we set , , , , and , the system will be chaotic. When the control parameter varies from 0 to 0.4, the changes of system are shown in Figure 23.

Figure 21: Bifurcation diagram when varies from 0 to 0.7 and .

Figure 22: Bifurcation diagram when varies from 0 to 0.3 and .

Figure 23: Bifurcation diagram when varies from 0 to 0.4 and .

The above figures indicate that the chaos of the system has been controlled when introducing the controlling parameter . With increasing, the system changes from the chaotic region to the bifurcation region. When is big enough, the system will fall into the steady state. For the manufacturer and retailer, they should try their best to control the system into the stable state and maximize their profits.

5. Conclusion

A supply chain system consisted of one manufacturer and two retailers have been discussed in this study. One is a traditional retailer and the other an online retailer. The complex dynamics of the system is studied by numerical simulation. The dynamic characteristics are investigated through 2D bifurcation diagram, the largest Lyapunov exponent, chaotic attractor, and time series. Finally, parameter adjustment control method is used to control the chaos.

In this paper, a pricing problem in the supply chain with one manufacturer who supplies an identical product to two retailers has been researched. For example, Huawei supplies Mate 9 to China Resources Vanguard Co. Ltd. and JD.com. CR Vanguard and JD.com are the typical representatives of traditional and online channels, respectively. These research results indicate that if the retailers compete with each other, the substitutability between the two sales channels will be generated. The competitions between retailers are helpful for the manufacturers. However, when the competition is fierce, the traditional and online channel can be replaced by each other, and the system will be reduced into a chaotic state. This article has demonstrated that consumer’s channel preference can affect the wholesale prices and sales commissions. If the consumers prefer any channel of the two, its wholesale prices and sales commissions may be raised. When the preference is too intense, the system will be unstable and even lead to chaotic phenomena. The parameter adjustment control method is performed on the parameters of wholesale price; the substitution of retailers and consumer’s channel preference make the system stable by adjusting the parameter . For all these works, the achievements of this study have both theoretical and practical significances for the manufacturer and retailers.

The conclusions of this paper can be helpful to the managers; however, there still exist some limitations. On the one hand, it considered a model of two-stage and dual-channel supply chain, which might be simple in practice. On the other hand, we only did the numerical simulation but did not make an empirical analysis of influence factors as an example. Hence, in the future, research based on large amounts of real data is significant. We may do further study in multi-stage and multi-channel supply chains in a practical market.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest regarding the publication of this paper.

Acknowledgments

The research was supported by the National Natural Science Foundation of China (71571131).

References

J. Ma and X. Ma, “Measure of the bullwhip effect considering the market competition between two retailers,” International Journal of Production Research, vol. 55, no. 2, pp. 313–326, 2016.

J. Ma and B. Bao, “Research on bullwhip effect in energy-efficient air conditioning supply chain,” Journal of Cleaner Production, vol. 143, no. 2, pp. 854–865, 2017.

T. Shu, H. Liao, S. Chen, S. Wang, K. K. Lai, and L. Gan, “Analysing remanufacturing decisions of supply chain members in uncertainty of consumer preferences,” Applied Economics, vol. 48, no. 34, pp. 3208–3227, 2016.

L. Feng, K. Govindan, and C. Li, “Strategic planning: design and coordination for dual-recycling channel reverse supply chain considering consumer behavior,” European Journal of Operational Research, vol. 260, no. 2, pp. 601–612, 2017.

S. Du, J. Zhu, H. Jiao, and W. Ye, “Game-theoretical analysis for supply chain with consumer preference to low carbon,” International Journal of Production Research, vol. 53, no. 12, pp. 3753–3768, 2014.

R. T. Frambach, H. C. A. Roest, and T. V. Krishnan, “The impact of consumer internet experience on channel preference and usage intentions across the different stages of the buying process,” Journal of Interactive Marketing, vol. 21, no. 2, pp. 26–41, 2007.

M. Okimoto, “International price competition among food industries: the role of income, population and biased consumer preference,” Economic Modelling, vol. 47, pp. 327–339, 2015.

M. Khouja, S. Park, and G. Cai, “Channel selection and pricing in the presence of retail-captive consumers,” International Journal of Production Economics, vol. 125, no. 1, pp. 84–95, 2010.

F. Gao and X. Su, “Omnichannel retail operations with buy-online-and-pick-up-in-store,” Management Science, vol. 63, no. 8, pp. 2478–2492, 2017.

H. Ke and J. Liu, “Dual-channel supply chain competition with channel preference and sales effort under uncertain environment,” Journal of Ambient Intelligence and Humanized Computing, vol. 8, no. 5, pp. 781–795, 2017.

S. S. Gan, I. N. Pujawan, Suparno, and B. Widodo, “Pricing decision for new and remanufactured product in a closed-loop supply chain with separate sales-channel,” International Journal of Production Economics, vol. 190, pp. 120–132, 2017.

J. Ma and H. Wang, “Complexity analysis of dynamic non-cooperative game models for closed-loop supply chain with product recovery,” Applied Mathematical Modelling, vol. 38, no. 23, pp. 5562–5572, 2014.

Y. Song and X. Zhao, “A newsvendor problem with boundedly rational strategic customers,” International Journal of Production Research, vol. 55, no. 1, pp. 228–243, 2016.

J. Ma and F. Si, “Complex dynamics of a continuous Bertrand duopoly game model with two-stage delay,” Entropy, vol. 18, no. 7, p. 266, 2016.

J. Ma and L. Xie, “The comparison and complex analysis on dual-channel supply chain under different channel power structures and uncertain demand,” Nonlinear Dynamics, vol. 83, no. 3, pp. 1379–1393, 2016.

J. Ma and L. Xie, “The impact of loss sensitivity on a mobile phone supply chain system stability based on the chaos theory,” Communications in Nonlinear Science and Numerical Simulation, vol. 55, pp. 194–205, 2018.

E. J. Anderson and Y. Bao, “Price competition with integrated and decentralized supply chains,” European Journal of Operational Research, vol. 200, no. 1, pp. 227–234, 2010.

C. H. Wu, C. W. Chen, and C. C. Hsieh, “Competitive pricing decisions in a two-echelon supply chain with horizontal and vertical competition,” International Journal of Production Economics, vol. 135, no. 1, pp. 265–274, 2012.

S. L. Yang and Y. W. Zhou, “Two-echelon supply chain models: considering duopolistic retailers different competitive behaviors,” International Journal of Production Economics, vol. 103, no. 1, pp. 104–116, 2006.

S. C. Choi, “Price competition in a channel structure with a common retailer,” Marketing Science, vol. 10, no. 4, pp. 271–296, 1991.

T. W. Mcguire and R. Staelin, “An industry equilibrium analysis of downstream vertical integration,” Marketing Science, vol. 2, no. 2, pp. 161–191, 1983.

T. Xiao, T. M. Choi, and T. C. E. Cheng, “Product variety and channel structure strategy for a retailer-Stackelberg supply chain,” European Journal of Operational Research, vol. 233, no. 1, pp. 114–124, 2014.

W. Huang, “Stabilizing nonlinear dynamical systems by an adaptive adjustment mechanism,” Physical Review E, vol. 61, no. 2, pp. R1012–R1015, 2000.