同质企业研发竞争行为的动力学性态分析

天津工业大学学报

ＪＯＵＲＮＡＬＯＦＴＩＡＮＪＩＮＰＯＬＹＴＥＣＨＮＩＣＵＮＩＶＥＲＳＩＴＹ

Vol.34No.3June2015

DOI：10.3969/j.issn.1671-024x.2015.03.016

同质企业研发竞争行为的动力学性态分析

黄东卫，付玉霞，于沈新

（天津工业大学理学院，天津

300387）

摘要：应用非线性动力系统中的分岔和混沌理论及相关的经济学理论，建立了同质三寡头R&D（研究与开发）竞

争在有限理性预期下的动力学模型，分析了该系统的均衡点及其稳定性，并以投入调整速度及技术溢出水平为控制参数，分析了系统产生分岔及混沌状态的条件，进一步得到了寡头间博弈引起的利润变化性态；结合Mathematics对不同参数条件下的系统进行了数值模拟，发现当企业的R&D投入调整速度过大时，系统由稳定状态进入倍周期分岔进而进入混沌状态；随着企业的技术溢出水平增大，其自身的利润下降，竞争对手的利润上升.模拟结果表明企业的利润在稳定状态下是稳定变化的，在混沌状态下则是连续波动的；当企业的投入调整速度过大时，其自身的平均利润呈下降趋势，竞争对手的平均利润呈上升趋势.

关键词：同质企业；竞争行为；动力学性态分析；研究与开发；有限理性；分岔；均衡点；Lyapunov指数；技术溢出中图分类号：F274；F224

文献标志码：A

文章编号：１６７１－０２４Ｘ（２０15）０3－００73－０5

DynamicsanalysisofstateofhomogenouscorporateR&Dcompetition

HUANGDong-wei，FUYu-xia，YUShen-xin

（SchoolofScience，TianjinPolytechnicUniversity，Tianjin300387，China）

Abstract：Weapplythetheoryofbifurcationandchaosofnonlineardynamicalsystemsandrelatedtheoryofeconomics，

andthenestablishaR&Dcompetitiondynamicmodelunderthelimitedrationalexpectationsofthreehomogeneousoligarchsandanalyzethissystem′sequilibriumsandtheirstability.WealsotakeadjustingspeedofR&Dinputandtechnicalspilloverslevelasthecontrolparameterandanalyzetheconditionofthesystembifurcationandchaos，thenwegetthechangingbehaviorofprofitascompetitionamongoligarchs；UndertheconditionofdifferentparameterswesimulatethesystembyMathematicsandillustratetheresult.WefindthatwhentheadjustingspeedofR&Dinputoftheenterpriseistoolarge，thenthesystembecomesfromstablestateintoperiod-doublingbifurcationandchaoticstate.Asthetechnicalspilloversleveloftheenterpriseincreases，itsprofitfallsandtheopponent′sprofitrises.Inaddition，theresultsofsimulationshowthatenterprise′sprofitschangessteadylyunderthestablestateandchangesvalatily，underthechaoticstate.Whentheadjustingspeedofinputoftheenterpriseistoolarge，itsaverageprofitfallsandtheopponents′rises.

Keywords：homogenouscorporate；competition；dynamicsanalysis；researchanddevelopment（R&D）；boundedratio鄄

nality；bifurcation；equilibriums；Lyapunovindex；technicalspillovers

研究和开发（R&D）已经成为企业发展的主力，是

企业核心竞争力和科技提升的重要科技手段.企业开展R&D活动可以获取更低的生产成本，更高的生产效率，从而获取更大的利润.然而，企业的R&D成果具有公品的部分属性，这就意味着在进行R&D活动的过程中，一定存在着技术溢出.技术溢出效应是经济市场中不可避免的现象，所以，这个因素也是不可忽略的.中外学者对此做了很多研究.J．A.SchumPeter创新的技术理论推动了R&D的相关研究的发

收稿日期：2014-12-11

展.Dssgupta根据创新收益的不同结构，把R&D分为竞争和非竞争两种类型.D’Aspremont和Jac-quemin最早研究了R&D合作问题，提出了R&D两阶段理论[1].Ziss（1944）和Bischi（2002）建立了具有技术溢出水平的R&D博弈模型[2].Trshler[3]和Luckraz[4]从外生溢出的角度讨论了R&D的竞争和合作.国内不少学者对R&D的决策问题也做了一些研究.霍沛军等[5]研究了双寡头R&D合作与非合作时最优溢出问题.盛昭瀚等[6]运用非线性动力系统理论对R&D投入动态竞争

基金项目：国家自然科学基金资助项目（11102132）

通信作者：黄东卫（1966—），男，教授，博士，研究方向为随机动力系统理论及应用.E-mail：tjhuangdw@163.com

—74—

天津工业大学学报第34卷

系统进行了复杂性分析.胡荣等[7]对异质产品双寡头企业R&D最佳技术含量的博弈进行了分析.然而，多数文献主要研究了双寡头R&D的决策问题，对于三维、四维等更高维寡头之间的竞争却研究较少.而真实的市场中往往都是多于两个寡头的竞争.本文对不同溢出水平下有限理性同质三寡头的两阶段动态竞争模型作了复杂性分析，并用数值模拟分析了投入速度、技术溢出水平对企业利润及平均利润的影响.

1模型的建立

假设市场上竞争的三寡头企业生产的产品是同

质的，竞争中存在技术溢出.由D’Aspremont和Jacquemin的经典R&D两个阶段研究：R&D活动分为

R&D投入和产量决策两个阶段.同时，假设三企业的R&D投入都是有限理性的，即企业第二阶段的决策是根据边际利润效应做出判断的.

设第i家企业（i=1，2，3）在t时期的价格为p（it），q（it

）表示企业（ii=1，2，3）在t时期生产的产品产量，市场的总产出Q（t）=q（1t）+q（2t）+q（3t）.设三企业在t时期的逆需求函数为p（it

）=a-bQ（a，b＞0）.其中：a为市场最高价；b为市场需求弹性.设第i家企业（i=1，2，3）在t时期的第一阶段的研究投入量为xit）（≥0）.博弈过程中存在着技术溢出，为了讨论方便，设企业i对企业j和企业k的技术溢出水平都为α（i0≤α≤1；i，j，k=1，2，3；i≠j≠k

），则企业i通过企业j和企业k的技术溢出使本企业的投入增加量为αix（jt）+αkx（kt）.设未进行R&D活动时的三企业的单位产品成本均为c，开展R&D活动后企业i单位产品的成本减少量为εijk（t）=β（x（it）+αjx（jt）+αkx（kt））（β≥0），则成本变为C（it

）=c-εijk（t）.本文假设企业i的研究投入的成本是非线性的：λ（ix（it））22（i=1，2，3）.于是三企业在t时期第二阶段的利润函数为π（it）=（p（it）-C（it））q（it

）-λ2

ixi（t）

2（i=1，2，3）（1）

根据利润最大化原则，分别令鄣π

（it）鄣q（=0（i=it）1，2，3）可以得到t时期的Nash-Cournot均衡产量为

q*

ijk=（[a-c）+x（it）β（3-2αi）+x（jt）β（2αj-1）+x（kt）β（2αk-1）]/4b（i，j，k=1，2，3；i≠j≠k）

（2）

把（2）代入（1）得在t时期第一阶段企业i选择投

入量为x（it

）（i=1，2，3）时利润最大，为2

π*

）=b（q*

）2-λixi（t）i（tijk（t）2（i=1，2，3）（3）

又知3个企业都采用有限理性预期，企业的研究投入不是通过精确计算得出的使利润最大化的决策，而是根据边际利润做出判断，边际利润大于零时，企业进行投入，边际利润小于零时，企业不投入，于是企业在t+1期的投入量为

x（it+1）=x（it

）+*

γ鄣π（x1，x2，x3）ix（it

）i鄣x（i=1，2，3）（4）

i

式中：γi＞0（i=1，2，3）为企业i在博弈过程中R&D投入调整速度.由式（1）—式（4）得到同质三寡头R&D博弈动态模型为：

*

x（qit+1）=x（it）+γix（i

t）（ijk（t）β（3-2αi）

2-λix（it

））（i=1，2，3）（5）

2不动点及局部稳定性分析

为了研究系统的不动点的局部稳定性，首先计算

系统（5）的雅克比矩阵

J（x1，x2，x3）=

3JJ11

JJ21J12J31J22J13J233233

󰀦其中

*

J（3-2αii=1+γ（i

β2x（i3-2α）i

28b+qijkβi）2-xiλi）（i=1，2，3）

Jij=β2

xiγ（

i2αi-3）（2αj-1）8b（i，j=1，2，3，i≠j）（6）

令x（it+1）=x（it），求得系统的解析解比较繁琐，为了简便，给某些参数赋值.取α1=0.3，α2=0.2，α3=0.25，β=1.2，λ1=2，λ2=2.3，λ3=2.2，a=30，b=2，c=4得到系统的8个不动点E（10，0，0）E（20，0，2.9771）E30，2.9972，0）E（40，2.8073，2.7457）E（53.1586，0，0）E（62.9107，2.8361，0）E（72.9538，0，2.8148）E（82.7356，2.6655，2.6070）.显然，这些不动点和γ（ii=1，2，3）无关；其中Ei（i=1，2，3，4，5，6，7）称有界均衡点，E8是唯一的Nash均衡点[8-9].

定理有界均衡点E（ii=1，2，3，4，5，6，7）是不稳

定均衡点.

证明

将E（10，0，0

）代入雅克比矩阵得（（第3期

黄东卫，等：同质企业研发竞争行为的动力学性态分析

—75—

3\"\"\"1+4.68γ1

0

0$%%%J（E\"%1）=

\"\"\"01+5.07γ2

0

%%%\"\"\"#

0

01+4.875γ%3

%%&

则它的特征值分别是1+4.68γ1，1+5.07γ2，1+4.875γ3，又因为参数γi＞0（i=1，2，3），所以这3个特征值的模都是大于1的，所以该不动点是不稳定的.

同理，将E（ii=1，2，3，4，5，6，7）分别代入雅克比矩阵计算出它们各自的特征值（不再一一列举出来），结果发现特征值的模都至少有一个大于1[10]，即E（ii=1，2，3，4，5，6，7）是不稳定的.

下面讨论Nash均衡点E8的稳定性情况，将E8代入雅克比矩阵得

31-4.0530γ$\"1-0.3545γ1-0.2954γ1

%\"%\"%J（E\"%8）=\"-0.2495γ21-4.50γ2-0.3119γ2

%\"%\"%\"-0.2346γ%\"\"3-0.3519γ31-4.2690γ3%%#

&

它的特征方程为（fλ）=λ3+Aλ3+Bλ+C，其中A=4.0530γ1+4.50γ2+4.2690γ3

B=-8.1060γ1-9.0179γ2-8.5379γ3+

18.1862γ1γ2+17.2328γ1γ3+19.1387γ2γ3+3C=4.053γ1+4.50γ2+4.2690γ3-18.1862γ1γ2-17.2328γ1γ3-19.1387γ2γ3+76.9308γ1γ2γ3-1

根据Routh-Hurwitz稳定性判据，系统在E8渐近稳定的充要条件是特征方程的所有零点落在复平面的单位圆内，必须满足以下4点[11]

＞))（f1）=A+B+C+1＞0))))-f（-1）=-A+B-C+1＞0()))C2-1＜0

（7）

（)))*

1-C2）2-（B-AC）2＞0

根据式（7）解出来的γ（ii=1，2，3）就是系统的稳定区域.稳定的Nash均衡点表明：尽管三寡头初始的R&D投入量不同，但是经过重复博弈，企业每次对初始投入量的调整最终将稳定在Nash均衡点处的稳定的投入量水平，一旦企业的调整速度超出了稳定区域，系统就会进入倍周期分岔甚至混沌状态.

3数值模拟及模型的复杂动力学特性分析

为了更好地了解系统（5）的动态演化过程，识别

系统的混沌行为[12-13]，本文采用数值模拟来分析.给出部分参数的值，作出R&D投入量的分岔图和李雅普诺夫指数谱.取数值模拟的初始参数值为α1=0.3，α2=0.2，α3=0.25，β=1.2，λ1=2，λ2=2.3，λ3=2.2，a=

30，b=2，c=4，γ2=0.3，γ3=0.4，3个寡头的初始R&D投入量分别是x1=1.7，x2=2，x3=2.1.

3.1R&D投入调整速度对系统状态的影响

为了研究企业R&D的调整速度对系统（5）的状态影响，本文选取企业1的R&D调整速度作为研究对象，企业2和企业3的R&D的调整速度对系统（5）的影响可以用同样的方法做.企业1的R&D投入量的分岔图如图1所示.

43企业1

企业2

量入投2企业3

D&R100.0

0.2

0.40.60.8

γ1

图1R&D投入量的分岔图

Fig.1R&Dinputbifurcationdiagram

观察图1可知，系统（5）随着企业1的R&D投入调整速度γ1的连续变化呈现稳定、倍周期分岔、混沌状态.γ1较小时，系统稳定；γ1较大时，系统出现分岔进而出现混沌状态.同时也可以用系统的Lyapunov指数谱（图2）进一步阐述该过程的变化趋势.Lyapunov指数是用来表征相轨道间分离快慢的指标，如果Lya－punov指数为正，则系统为混沌状态，如果Lyapunov指数为负，则系统处于稳定状态.

0.4数0.2指von0.0upayL-0.2大最-0.4

0.1

0.2

0.3

0.4

0.5

0.6

0.7

γ1

图2系统的最大Lyapunov指数谱Fig.2

MaximalLyapunovexponentofsystem

从图2可以看到，企业1的R&D投入调整速度

γ1较小时，系统的最大Lyapunov指数小于零，系统处于稳定状态；γ1=0.49时，系统的最大Lyapunov指数

等于零，系统进入倍周期分岔状态；当γ1增大到0.62时，系统的最大Lyapunov指数开始大于零，系统进入混沌状态.这与图1观察到的结果是一致的.3.2R&D投入调整速度对企业利润的影响

仍然以企业1的R&D投入调整速度为研究对

—76—

天津工业大学学报第34卷

象，模拟三寡头企业的利润在不同的γ1取值下，随博弈周期变化趋势图如图3所示.选取博弈周期T=100，分别取γ1=0.32和γ1=0.66，作出企业利润在这两个取值下的变化趋势图.由图1和图2知，γ1=0.32对应系统（5）的稳定状态，γ1=0.66对应系统（5）的混沌状态.

2624企业2

润利22的业企业1

企20企业3

1816

020

406080100

博弈周期T

（a）γ1󰀃=󰀃0.32时三企业的利润变化趋势图

2624润利22的业企20

1816

020

4060

80100

博弈周期T

（b）γ1=0.66时三企业的利润变化趋势图

21.7638

321.7638π21.763721.763721.7636

23.521.0

23.0π2

π1

21.5

22.5

（c）γ1=0.32时对应的三企业的利润

233π222521

2419

23

π2

π201

21

22

22

（d）γ1=0.66时对应的三企业的利润

图3

三企业的利润变化趋势图及相图

Fig.3

Changingprofitofthreeenterprisesandphasediagram

由图3（a）可知，γ1=0.32即系统处于稳定状态

时，三企业的利润随博弈周期稳定变化；由图3（b）可

知，γ1=0.32系统处于混沌状态时，三企业的利润连续波动.图3（c）和图3（d）是对应的三企业的利润的相图，γ1=0.32时相图基本是一条直线，说明利润稳定变化，γ1=0.32时相图是混沌吸引子，说明利润是波动的，这和图3（a）、图3（b）的结论的是一致的.

为了进一步说明企业1的R&D投入调整速度γ1、γ2和γ3对企业利润的影响，模拟企业平均利润在100个博弈周期内随着γ1、γ2和γ3的变化趋势图，如图4所示.

24企业2

润利22均平的20企业1

企业3业企1816

0.10.20.3

0.40.50.60.7

γ1

（a）平均利润随γ1的变化趋势图

24企业2

润利22均平1

企业3

的20企业业企1816

0.1

0.2

0.3

0.4

0.5

0.6

γ2

（b）平均利润随γ2的变化趋势图

24企业2

润利22均平的20企业1

企业3

业企1816

0.10.20.30.40.50.6

γ3

（c）平均利润随γ3的变化趋势图

图4平均利润随γ1、γ2、γ3时变化趋势图Fig.4Averageprofitwithchangeofγ1，γ2，γ3

从图4（a）可知，企业1R&D投入调整速度γ1较

小时，系统处于稳定状态，三企业的平均利润稳定变化；随着γ1的增大，系统进入倍周期分岔甚至混沌状态，企业1的平均利润呈下降趋势，但企业2与企业3的平均利润呈上升趋势.这说明企业的投入调整速度过大时对企业自身是不利的，对竞争对手是有利的，图4（b）与图4（c）也说明了这个结果.

第3期黄东卫，等：同质企业研发竞争行为的动力学性态分析

—77—

3.3

技术溢出水平对企业平均利润的影响

以企业1的技术溢出水平α1为研究对象，分别选

取α2=0.2和α2=0.6，α3=0.25和，α3=0.7，模拟三企业的平均利润随着α1的变化趋势图，其他参数不变，选择博弈周期T=100，见图5.

40润30企业3利企业2

均平20企业1

的业企100

0.0

0.20.4

0.60.81.0

α1

（a）α2=0.2，α3=0.25三企业的平均利润变化趋势图

403

润30企业1

企业利企业2

均平的20业企100

0.0

0.2

0.4

0.6

0.8

1.0

α1

（b）α2=0.6，α3=0.25三企业的平均利润变化趋图

40润30企业1

企业2

利企业3

均平的20业企100

0.0

0.20.4

0.60.81.0

α1

（c）α2=0.2，α3=0.7三企业的平均利润变化趋势图

40

企业1

润30企业2

利均平的20企业3

业企100

0.0

0.20.4

0.60.81.0

α1

（d）α2=0.6，α3=0.7三企业的平均利润变化趋势图

图5

三企业的平均利润随α1的变化趋势图

Fig.5

Averageprofitofthreeenterpriseswithchangeofα1

比较图5（a）和图5（b）可以看出，当企业3的技术溢出水平固定为α3=0.25，企业2的技术溢出水平分

别为α2=0.2和α2=0.6时，企业1的平均利润随着它自身的技术溢出水平α1的增加呈下降趋势；企业2和企业3的平均利润呈上升趋势；图5（b）中三企业的大体趋势和图5（a）是一样的，但是企业1和企业3的平均利润整体上高于图5（a）中企业1和企业3的平均利润，企业2的平均利润整体上则相对降低，同时企业2与企业3之间的平均利润差距变大.这说明，企业的技术溢出使自身利润下降，则对竞争对手是有利的.反过来，如果企业降低自己的技术溢出水平，对手的利润就会降低，由于竞争者同样的心理，对手也会降低自身的技术溢出水平来降低对方的利润，最终各企业将在技术溢出水平相等时达到博弈均衡.比较图5（a）和图5（c），图5（b）和图5（d）亦能得到同样的结论.

由此可见，企业之间是相互影响又相互依存的，只有企业间相互合作，达到知识共享，资源共享，技术

共享，优势互补，才能共赢，才能为社会创造更多的财富和价值.

4结束语

本文建立了不同溢出水平下同质三寡头在有限理性预期下R&D竞争模型，对此模型作了复杂的动力学分析，发现当企业的R&D投入调整速度过大时，系统由稳定状态进入倍周期分岔进而进入混沌状态；随着企业的技术溢出水平增大，它自身的利润下降，对手的利润上升，对手的技术溢出水平变大时，它的利润也会相对变大.同时用数值模拟分析了企业的R&D投入调整速度及技术溢出水平对企业利润和平均利润等的影响.结果表明：企业的利润在稳定状态下是稳定变化的，在混沌状态下则是连续波动的；企业的平均利润与自身的调整速度密切相关，当某个企业调整速度过大时，自身的平均利润呈下降趋势，竞争对手呈上升趋势.企业在竞争的过程中，企业的技术溢出水平对自身和对对手的影响都是利大于弊的.企业间只有合作，形成资源共享，知识共享，技术共享才能共赢.本文仅就同质三寡头R&D之间的竞争做

了一些分析，对于产品有差异性企业间的竞争，不同预期条件下企业间的竞争及更高维寡头之间的竞争没有作分析，这些问题还需要进一步的研究.参考文献：

[1]

D′ASPREMONTC，JACQUEMINA.Cooperativeandnon-

（下转第84页）

第3期张芳，等：双渠道供应链博弈模型的复杂性分析

—83—

󰀁󰀁󰀁󰀁p（mt+1）=p（mt）+β2p（mt

）[am+b2ar/2b1+󰀡b22/b

1

-2b1󰀢p（mt

）+b2ω（t）+（2b2

-b2

12-b1b2）c/2b1]

（18）

其中，零倍商决策为p（rt）=[ar+b2p（mt

）+b1ω（t）]/2b1式中：β1为ω的调整系数；β2为pm的调整系数.

由系统（18）可见，制造商首先决策ω和pm，其调整参数为β1、β2，零售商看到制造商的决策后再来决策pr.3.2.2动态模型分析

系统（18）有4个均衡点：E0=（0，0），此时，Pr=ar2b1

E1=（ar+（b1-b2）c，0），P=3ar+（b1-b2）c2br4b，11

b2

2

E1am+A+（2b1-b1b2-b2=（0，2）c

2F），22

P+b1b2am+b2A+b（22br=ar1-b1b2-b2）c2b1

4b，

1F

E*

=（ω*

，p*

*

m），Pr=pr，

式中：ω*=B+12C；p**

m=A2F+12C；pr2Fr=a4b+

1

B2F+b1+b24bC.1

类似集中决策下均衡点的证明，可以得出分散决策下E0、E1、E2是不稳定的，E*是Nash均衡解.

设a=400，θ=0.45，b1=5，b2=3，c=30，此时均

衡解为ω*=63.75，p**

m=66.25，pr=69.75，由Jury条件解出稳定域为0<β1<8/1275，0<β2<8/2173，则系统18）关于参数β1、β2的稳定域如图7所示.

0.0050.0040.0032β0.0020.001

0.000

0.0000.0010.0020.0030.0040.0050.0060.007

β1

图7分散决策下Nash均衡点稳定域

Fig.7Nashequilibriumregionofstabilityunderdecentralized

decision

固定β2=0.003，可以得到价格ω、pm、pr随β1变化的分岔图以及利润πm、πr、πc随β1变化的分岔图，如图8所示.

由图8可以看出，当β1=0.00318时，价格ω、pm、

1009080X：0.00318

Y：66.59

70pr

rp，m60pω

pm

，ω5040302010

00.51.01.5

2.02.53.03.54.04.5

β1/×10-3

（a）ω，pm，pr随β1变化的分岔图

60005000πc

4000πm

X：0.00318Y：4757

cπ3000，rπ，2000mπ1000

πr

0-1000-20000

1

2

34

5

6

β1/×10-3

（b）πm，πr，πc随β1变化的分岔图

图8价格ω、pm、pr和利润πm、πr、πc随β1变化的分岔图Fig.8

Bifurcationdiagramsofω，pm，prandπm，πr，πcwithrespecttoβ1

pr和利润πm、πr、πc均进入倍周期分岔.

固定β1=0.001，可得价格ω、pm、pr随β2变化的分岔图以及利润πm、πr、πc随β2变化的分岔图，如图9所示.

9080X：0.0027prpm

70Y：70.03

rp60ω

，mp，50ω403020100

0.5

1.0

1.5

2.02.53.03.5

4.0

4.5

β2/×10-3

（a）ω、pm、pr随β2变化的分岔图

60005000πc

4000πm

X：0.00357cπ3000Y：4777

，rπ，2000mπ1000

πr0-1000-20000

1

2

34

5

6

β2/×10-3

（b）πm、πr、πc随β2变化的分岔图

图9价格ω、pm、pr和利润πm、πr、πc随β2变化的分岔图Fig.9

Bifurcationdiagramsofω，pm，prandπm，πr，πcwithrespecttoβ2

（