非线性需求下多寡头博弈模型的动力学分析论文

非线性需求下多寡头博弈模型的动力学分析

路正玉，于欢欢，王文静，赵 娜

(兰州交通大学数理学院，甘肃兰州 730070)

摘 要： 博弈论可以为寡头企业定价提供正确的指导．本文利用博弈论的相关理论分析了多寡头博弈模型，并通过分岔图、最大Lyapunov 指数图和混沌吸引子研究了系统在不同参数范围内的动力学特性．结果表明，随着调整速度、敏感度以及边际成本的增大，系统呈现分岔现象甚至混沌态，即当相关参数超过一定范围，市场处于混乱状态．据此可制定对企业最有利的价格策略．

关键词： 非线性需求；价格博弈；复杂性；价格决策

在社会和科技不断发展的同时，经济也在不断发展，像银行、保险等．伴随着经济的发展，出现了一些严重的社会问题，价格竞争^[1]就是其中一个．价格竞争的出现源自于各大厂商都想着扩大自己的市场，想通过调整自己产品的价格和其它厂家竞争，于是就出现了寡头垄断现象^[2]．市场往往被几个寡头垄断，这是市场普遍存在的一种形式．López 等^[3]基于经济学家Bertrand 的研究提出了一个双重寡头垄断模型，这个模型是以通过价格竞争来赢得市场为目的的一种模型．消费者喜欢选择价格低的产品，寡头垄断者就通过降低他们产品价格来赢得更多的消费者．Agiza等^[4-5]对多寡头博弈下市场稳定存在时商品的价格变化区间进行了研究，在该区间发现了混沌现象的存在．Ahmed 等^[6-10]通过研究市场上各厂家为了吸引消费者所作的广告，得出了一个新的模型：广告竞争模型，并且还得出了Lyapunov 指数．通过这些研究结论能够分析出价格的复杂变化．杨树等^[11-14]研究了房地产市场的一些垄断情况，发现真正的纳什均衡可以在有界有理不完全信息的动态重复博弈下进行．

崔国祯等[3]通过对斑马鱼肠下静脉（SIVc）血管直径、内皮细胞数及节间血管（ISVs）的观察，发现不同浓度的丹参素衍生物（ADTM）对于正常、损伤型模型均能起到增加与保护作用。可以激活胰岛素信号通路，发挥血管新生作用。

企业之间的行为都是相关的，它们之间互相影响．企业主要通过价格竞争在市场上相互制约，商品合理的定价直接影响到价格竞争的成功与否．在当今的寡头市场上，产品的性质基本上是相同的，每一家企业的定价决策和目的都很容易被其它企业觉察到．可以看出，现在的寡头企业的市场份额之间存在着一种促进和抑制的关系，所以，寡头企业的生存与发展的着重点就在定价的策略上．寡头企业定制价格，不仅要有丰富经验的先驱决策者，还要有先进的科学的指导理论，博弈论就是能够为寡头企业定价策略提供正确指导的科学理论．本文将利用博弈论的相关理论建立模型，分析企业在同样的产品性质下的价格定制策略．

1 模型的建立

假设市场上有两个生产同类产品的企业，企业1 和企业2 生产的产品具有一定的替代关系，生产的产品分别记为A 和B．一个寡头企业不仅影响自己的决策，而且影响其它寡头企业的决策，甚至改变整个市场的运行状态．产品价格、销售量分别记为p_i ,q_i (i =1, 2)，则寡头企业1 和企业2 的线性需求函数可设为：

本文利用WRF V3.5模式中的双参MOR方案和MY方案分别对2010年7月20—21日辽宁地区的一次强降水过程进行数值模拟,并在文中从地面累积降水量、降水强度、模拟强降水各种水凝物粒子随高度的分布以及最后对雨滴的一个细致的分析来对两方案进行了详细的对比。对比结果可以总结为以下几点:

假设企业i (i =1, 2)成本函数形式为c_ip_i ，其中c_i (i =1, 2)为企业技术水平对价格的敏感系数，则寡头企业i 的利润函数表达式为：

整理可得：

根据判断系统状态的最大的Lyapunov 指数，得到相对于图1（a）的最大Lyapunov 指数如图1（b）所示．最大Lyapunov 指数的不同表示系统的状态不相同，当最大Lyapunov 指数小于0时，系统处于稳定状态，最大Lyapunov 指数大于0时，系统处于混沌状态．

由于大型渠道混凝土衬砌施工作业多在1∶2或更陡的坡面上进行切缝施工，且自动行走，能使其平稳地行驶在坡面上，配重调整过重过轻、配重的位置都是影响实际应用的主要因素，所以机具的配重尤为重要。

其中p_i (t + 1)和p_i (t )分别表示第t + 1期和第t 期的价格．第t 期，原寡头企业的边际利润函数为：

令∂π ₁∂p ₁=0，∂π ₂∂p ₂=0时就能得到企业1 和企业2 的反应方程：

其中a ₁＞0和a ₂＞0分别表示产品的市场容量；α ₁＞0和α ₂＞0分别表示企业产品A 和B 对自身价格的敏感度；β ₁＞0和β ₂＞0分别表示企业产品A 和B 对对方产品价格的敏感度．在两个企业同时决策的情况下，假定两企业的边际成本分别为c ₁和c ₂，没有固定的成本．

依据以上分析，得出如下基于不同预期的双寡头投资动态调整模型：

2 模型分析

产品A 产品B 存在着一定的替代关系，但不是相同性质的产品，此时α ₁,α ₂,β ₁,β ₂均为有限的正实数．考虑反应方程：

由隐函数定理得，当D ≠ 0时，反应方程组能够确定唯一的一个函数：p ₁(α ₁,α ₂,β ₁,β ₂)，p ₂(α ₁,α ₂,β ₁,β ₂)，所以，可以确定唯一的一组均衡价格p ₁，p ₂．

3 数值模拟

为了更好地了解系统的动力学特性，在这部分，主要通过模拟系统在不同参数范围下的复杂动力学行为来定性地分析系统的动力学特性．

3.1 调整速度v ₁和v ₂对系统稳定性的影响

下面讨论寡头厂商A（B）所生产的产品对寡头厂商B（A）生产的产品的替代率β ₁（β ₂）对系统稳定性的影响．

因为市场上各企业占据了不同的市场信息，市场信息表现为不完全性，对于寡头垄断者的企业也是如此，在有限的信息下，他们在市场竞争中只能对价格的调整策略作出合理的决策．当两个企业的价格调整策略差不多时，认为他们是采用了有限的决策方式，主要就是通过市场上的利润反馈信息来对价格进行调整．下一期的价格调整信息是由这一期的市场利润和价格信息分析得出的．两个企业采用上述规则进行决策时，也就是说如果边际利润为正，市场上产品的需求量就会增加，相反，如果边际利润为负，市场上产品的需求就会减少．下列式子中v ₁和v ₂为价格调整速度，价格的调整策略可以表述为：

通过高中数学课程的学习，学生能够在具体的真实情境中积累经验，形成抽象思维，把握事情的本质，运用数学抽象的思维方式思考并解决实际问题．在2018年的高考数学命题中，对学生数学抽象核心素养的考查已经落实到具体的题目当中，凸显高考数学命题鲜明的素养导向和育人导向．

图2 是系统寡头企业商品售价p ₁和p ₂关于参数v ₂的分岔图和相应的最大Lyapunov 指数图．通过不断调整其它参数，固定一组参数不变：a ₁=0.2，a ₂=0.2，α ₁=0.8，α ₂=1,β ₁=0.3，β ₂=0.2，v ₁=0.1，c ₁=0.2，c ₂=0.3，绘制系统（7）关于v ₂的单参数分岔图，如图2（a）所示．从图2 中可以看出，当参数0＜v ₂＜1.36时，系统处于稳定周期，市场表现为稳定有序盈利，有利于社会的稳定发展，企业之间的相对竞争力度较小，各企业相对稳定运营；当参数1.36＜v ₂＜1.71时，就变成了混沌状态．图2（a）的最大Lyapunov 指数如图2（b）所示．

图1 系统（7）关于参数v ₁的分岔图和相应最大Lyapunov 指数图
Fig 1 Bifurcation Diagram and Corresponding Maximum Lyapunov Exponent Diagram about Parametersv ₁of System (7)

图2 系统（7）关于参数v ₂的分岔图和相应最大Lyapunov 指数图
Fig 2 Bifurcation Diagram and Corresponding Maximum Lyapunov Exponent Diagram about Parametersv ₂of System (7)

3.2 敏感度₁β 和₂β 对系统稳定性的影响

图1 是系统寡头企业的商品售价p ₁和p ₂关于参数v ₁的分岔图和相应的最大Lyapunov 指数图．通过不断调整其它参数，尝试着用Matlab 绘制一些图形．取一组参数不变：a ₁=0.4，a ₂=0.3，α ₁=0.9，α ₂=0.9，β ₁=0.05，β ₂=0.7，c ₁=0.4，c ₂=0.5，v ₂=0.2，绘制系统（7）关于v ₁的单参数分岔图，如图1（a）所示．从图1（a）中看出，当参数0＜v ₁＜2.42时，系统是很稳定的，当系统处于稳定的区间时，市场表现为稳定有序，有利于社会的稳定发展，企业之间的相对竞争力度较小，各企业相对稳定存在；当参数2.42＜v ₁＜2.9时，系统就出现了混沌现象，当混沌现象发生时市场就会进入到一种混乱无序的杂乱状态，在这种情况下，无论是企业自己还是他们的竞争者都无法避免受混乱状态的影响，尤其对市场上的一些根基并不稳固的小厂家来说，生存是很困难的，市场的不稳定影响的就不只是市场了，整个社会都会受其影响，这个时候，寡头企业1 的管理者就应该控制价格的相关调整策略和他们产品的价格调整速度，尽量不要因为眼前的利益发动价格战，进而产生市场的混乱．

窖泥中微生物总数和耐热芽孢杆菌数对窖池发酵有着重要影响，但并不是微生物数量越多，窖泥的质量就越好[13]。窖泥微生物在生产过程中不断被驯化，优胜劣汰，并非所有优势微生物都对白酒生产有促进作用[14-15]，如同白酒中含量最多的风味物质并非对白酒风格起决定性作用[16]。因此窖泥中微生物指标作为一个参考，并没有赋予较高的权重。

图3 系统（7）关于参数 ₁β 的分岔图和相应最大Lyapunov 指数图
Fig 3 Bifurcation Diagram and Corresponding Maximum Lyapunov Exponent Diagram about Parameters ₁β of System (7)

图4 是系统寡头厂商商品售价p ₁和p ₂关于参数β ₂的分岔图和相应的最大Lyapunov 指数图．通过不断调整其它参数，固定一组参数不变：a ₁=0.6，a ₂=0.3，α ₁=0.5，α ₂=1.5，β ₁=0.5，v ₁=0.2，v ₂=0.3，c ₁=0.1，c ₂=0.1，绘制系统（7）关于β ₂的单参数分岔图，如图4（a）所示，图形的形状与图3 类似．从图中可以看出当参数0＜β ₂＜4.17时，系统是很稳定的，当系统处于稳定的区间，市场表现为稳定有序；当参数4.17＜β ₂＜4.22时，系统表现为混沌状态，当混沌现象发生时，市场就会进入到一种混乱无序的杂乱状态，在这种情况下，无论是厂商自己还是他们的竞争者都无法避免受混乱状态的影响，尤其对市场上的一些根基并不稳固的小厂家来说生存很是困难，市场的不稳定影响的就不只是市场了，整个社会都会受其影响，在这个时候，寡头厂商B 管理者应尽力控制产品替代率，尽量不要因为眼前利益就发动价格战，进而产生市场的混乱．图4（a）的最大Lyapunov 指数如图4（b）所示．

图3 是系统寡头厂商商品售价p ₁和p ₂关于参数β ₁的分岔图和相应的最大Lyapunov 指数图．通过不断调整其它参数，取一组参数不变：a ₁=1 ，a ₂=0.5，α ₁=2 ，α ₂=0.5，β ₂=0.5，v ₁=0.2，v ₂=0.2，c ₁=0.2，c ₂=0.1，绘制系统（7）关于β ₁的单参数分岔图，如图3（a）所示．可以看出，当参数0＜β ₁＜6.42时，系统是很稳定的，当系统处于稳定的区间，市场表现为稳定有序；当参数6.42＜β ₁＜6.79时，就变成了混沌状态，当混沌现象发生时市场就会进入到一种混乱无序的杂乱状态，在这种情况下，无论是厂商自己还是他们的竞争者都无法避免受混乱状态的影响，尤其对市场上的一些根基并不稳固的小厂家来说，生存是很困难的，市场的不稳定影响的就不只是市场了，整个社会都会受其影响，在这个时候，寡头厂商A 管理者应尽力控制产品替代率，尽量不要因为眼前利益就发动价格战，进而造成市场混乱．根据判断系统状态的最大Lyapunov 指数，得到相对于图3（a）的最大Lyapunov 指数，如图3（b）所示．最大Lyapunov指数的不同表示系统的状态不相同，当最大Lyapunov 指数小于0时，系统处于稳定状态，最大Lyapunov 指数大于0时，系统处于混沌状态．

图4 系统（7）关于参数 ₂β 的分岔图和相应最大Lyapunov 指数图
Fig 4 Bifurcation Diagram and Corresponding Maximum Lyapunov Exponent Diagram about Parameters ₂β of System (7)

3.3 边际成本c ₁和c ₂对系统稳定性的影响

图5 是系统寡头厂商商品售价p ₁和p ₂关于参数c ₁的分岔图和相应的最大Lyapunov 指数图．通过不断调整其它参数，取一组参数不变：a ₁=3 ，a ₂=3 ，α ₁=4 ，α ₂=4 ，β ₁=0.8，β ₂=0.8，v ₁=0.4，v ₂=0.4，c ₂=0.2，绘制系统（7）关于c ₁的单参数分岔图，如图5（a）所示．从图中可以看出当参数0＜c ₁＜0.21时，系统是很稳定的，当系统处于稳定的区间，市场也就稳定了；当参数0.21＜c ₁＜0.318时，系统处于混沌状态，在混沌状态下，市场也就变得混乱了，各厂商都会受其影响，这个时候厂商A 就应该控制敏感系数c ₁，避免市场发生混乱．图5（a）的最大Lyapunov 指数如图5（b）所示．

接下来讨论寡头厂商A（B）的价格敏感系数c ₁（c ₂）对系统稳定性的影响．

资源实现共享，必定涉及到多方利益。这些所涉及到的利益包括：高校教师和学生，高校自身以及相关的政府部门。在资源使用的过程中，要秉承“谁使用，谁负责”的原则。协调好各方的利益是关键，如果出现利益受损的一方，则受损方应得到相应的补偿。因为在资源共享的过程中受益人或受益单位有支付相应成本的义务，正所谓“谁受益，谁负担”。对于教育资源使用的各方，只有在教育资源不受损、可以利用资源共同收益的状况下，才能提高各部门以及个人的积极性。

图5 系统（7）关于参数c ₁的分岔图和相应最大Lyapunov 指数图
Fig 5 Bifurcation Diagram and Corresponding Maximum Lyapunov Exponent Diagram about Parametersc ₁of System (7)

图6 是系统寡头厂商商品售价p ₁和p ₂关于参数c ₂的分岔图和相应的最大Lyapunov 指数图．取其它参数不变：a ₁=3 ，a ₂=3 ，α ₁=4 ，α 2=4 ，β 1=0.8，β ₂=0.8，v ₁=0.4，v ₂=0.4，c ₁=0.2，图形的形状与图6 类似．

从图6 中可以看出，当参数0＜c ₂＜0.21时系统是很稳定的，当系统处于稳定的区间时，市场也处于稳定状态；当参数0.21＜c ₂＜0.318时，系统处于混沌状态，在混沌状态下，市场也就变得混乱了，各厂商都会受其影响，这个时候厂商B 就应该控制敏感系数c ₂，以避免市场发生混乱．图6（a）的最大Lyapunov 指数如图6（b）所示．

图6 系统（7）关于参数c ₂的分岔图和相应最大Lyapunov 指数图
Fig 6 Bifurcation Diagram and Corresponding Maximum Lyapunov Exponent Diagram about Parametersc ₂of System (7)

3.4 对系统（7）混沌吸引子和相图的模拟

当参数v ₁的取值为4.5、4.4、4.3、4.2、4.1、4.0、3.9、3.8 时，系统（7）的吸引子和相图见图7．从图7 中得出，系统（7）表现混沌状态时，对应图（a）和图（b）中的混沌吸引子．

4.高温高湿是本病的一个重要诱因，因此，夏天霉雨季节(4～8月)应对10～15日龄雏鸡进行药物预防，以免暴发本病。

4 结 语

本文建立了一个寡头垄断博弈模型，分析了价格调整对寡头厂商利益以及整个市场的影响．数值模拟结果表明，一旦价格调整率超过稳定范围，整个市场将陷入混乱，寡头厂商自身利润也将减少．寡头垄断厂商或当前寡头垄断厂商应采取较低的价格调整率，以实现利润最大化．从上文可以得出，在寡头市场价格调整时会出现分岔甚至混沌现象，这对社会发展非常不利，有关部门应该规范市场，从而让市场和社会健康持续地发展下去．

所有寡头垄断者初始状态的微小差异都会对分析结果产生很大影响．为了稳定市场，寡头垄断者应慎重选择博弈的初始条件．当市场处于混乱状态时，寡头垄断者可以使用延迟反馈控制来控制混沌．

在有限的理论下，博弈模型是建立在销售量仅与价格有关的条件下，实际上影响销量的因素很多，如产品质量、市场需求等，这些问题需要进一步分析．

图7 系统（7）关于参数v ₁变化时的混沌吸引子、相图
Fig 7 Chaotic Attractors and Phase Diagram about Parametersv ₁of System (7)

参考文献

[1]孙晓东，田澎，赵藜．消费者异质下基于质量-价格竞争的定价策略选择[J]．系统管理学报，2013，22（3）：349-357．

[2]王必好，黄浩杰．寡头垄断市场结构的技术创新效应研究：基于伯特兰和古诺均衡分析视角[J]．经济评论，2013（5）：13-21．

[3]López M C, Naylor R A. The Cournot-Bertrand profit differential: A reversal result in a differentiated duopoly with wage bargaining [J]. Eur Econ Rev, 2002, 48(3): 681-696.

[4]Agiza H N, Elsadany A A. Chaotic dynamics in nonlinear duopoly game with heterogeneous players [J]. Appl Math Comput, 2004, 149(3): 843-860.

[5]刘峰，李亚光，王宏兴．非线性需求下四寡头价格博弈模型及其复杂特性[J]．系统工程学报，2016（6）：719-728．

[6]Ahmed E, Agiza H N, Hassan S Z. On modelling advertisement in Cournot duopoly [J]. Chaos Soliton Fract, 1999,10(7): 1179-1184.

[7]罗宇．不同市场规则下发电企业投资决策分析[D]．长沙：长沙理工大学，2008：26-29．

[8]吉伟卓．寡头垄断电力市场产量博弈模型及其混沌复杂性研究[D]．天津：天津大学，2008：9-36．

[9]马军海，吴可菲．中国啤酒市场的四寡头价格博弈及其延迟决策[J]．系统工程学报，2013，28（6）：717-728．

[10]徐伟，马军海．保险市场中三寡头垄断的动态博弈模型研究[J]．复杂系统与复杂性科学，2013，10（2）：52-58．

[11]杨树，梁樑，熊立．Cournot 和Bertrand 竞争下均衡质量的比较[J]．系统管理学报，2008（1）：104-109．

[12]杨哲．非合作博弈在行为经济学思想下的均衡研究[D]．重庆：重庆大学，2012：21-48．

[13]黄伟，陈广南．推力中止过程一维不定常气相流动数值计算[J]．固体火箭技术，1996（3）：12-16．

[14]高忠科．两相流复杂网络非线性动力学特性研究[D]．天津：天津大学，2010：68-89．

Dynamic Analysis of Multi-oligopoly Game Model under Nonlinear Demand

LU Zhengyu, YU Huanhuan, WANG Wenjing, ZHAO Na
(College of Mathematics and Physics, Lanzhou Jiaotong University, Lanzhou, China 730070)

Abstract: Game theory is a scientific theory that can give guidance for oligarchic enterprises in price customization. In this paper, the multi-oligopoly game model is analyzed by using the game theory, and the dynamic characteristics of the system in the range of different parameters are studied through bifurcation graph, maximum Lyapunov exponent graph and chaotic attractor. The results show that with the increase of adjustment speed, sensitivity and marginal cost, the system presents bifurcation and even chaos, that is, when the relevant parameters exceed a certain range, the market will be in a state of chaos. In this light, the most favorable price strategy can be formulated for enterprises.

Key Words: Nonlinear Demand; Price Game; Complexity; Price Decision

中图分类号： O193

文献标志码： A

文章编号： 1674-3563(2019)03-0036-09

DOI： 10.3875/j.issn.1674-3563.2019.03.006 本文的PDF 文件可以从http://xuebao.wzu.edu.cn/获得

收稿日期： 2019-02-26

基金项目： 国家自然科学基金(61863022)；中国博士后科学基金(2017M623276)；甘肃省科技支撑计划项目(17JR5RA096)

作者简介： 路正玉(1992- )，男，甘肃靖远人，硕士研究生，研究方向：混沌经济学

(编辑：王一芳)

标签：非线性需求论文;  价格博弈论文;  复杂性论文;  价格决策论文;  兰州交通大学数理学院论文;