摘 要 從生態學視角來研究企業競爭的動態演化已成為一個研究熱點.針對現有研究文獻中只討論二維系統的不足,以三個企業之間的競爭為例,構造了企業間競爭的三維模型,運用微分方程穩定性理論分析其穩定性,并由此揭示出企業間競爭的動力學機制.數值仿真結果表明,該三維模型能有效地模擬企業間競爭的動態演化規律.
關鍵詞 生態競爭方程;企業競爭;穩定性分析;動力學機制
中圖分類號 F062.2 文獻標識碼 A
Dynamic Evolution of Enterprises’ Competition Based on Ecological Competition Equations
YAO Yi
(School of Mathematical Sciences, Nanjing Normal University, Jiangsu, Nanjing 210046,China)
Abstract
The research on dynamic evolution of enterprises competition from the perspective of ecology has attracted much attention. Aiming atthe drawbacks of the existing studies, a three-dimensional competitive model for three enterprises was formulated. The stability of this model was analyzed using the qualitative theory of differential equations, which reveals the dynamical mechanism of enterprises competition. The numerical simulations suggest that this three-dimensional model can reflect the law of dynamic evolution of enterprises competition effectively.
Keywords ecological competition equations;enterprise competition; stability analysis; dynamic mechanism
1 引 言
競爭是現代經濟發展的動力,企業要在激烈的競爭中生存與發展,就必須了解所處的市場,深刻理解競爭的內在發展規律.考慮到經濟系統中企業間競爭與生態系統中種群間競爭的相似性,國內外一些學者嘗試利用生態學理論來研究企業競爭的動態演化規律,并取得了豐碩的成果[1-8].
文獻[1]利用修改的Logistic模型分析了IT企業集群的共生性分析,認為集群內存在若干個核心企業, 它們形成互補或替代關系, 眾多的衛星企業寄生在各自的核心企業下,核心企業和衛星企業在集群內能夠自動維持著這種分工協作的穩定性.文獻[2]利用Logistic模型,對中小企業與大型企業之間的惡性競爭、寄生關系、捕食關系和互惠共生四種互動模式進行數學描述,分析各種模式的穩定條件,并據此提出了相應的中小企業競爭策略的選擇.依據生態系統共生模型,文
獻[3]建立了一個包含正負作用的一般企業的競合模型,通過模型的分析說明了企業集群的發展取決于市場容納量與固定成長率,而這兩者的提高與企業集群初始狀態、企業集群環境與市場環境密切相關,同時,合作能增加企業集群的穩定性.文獻[4]在此前研究的基礎上,增加了企業產出水平的下臨界點,進行了企業競爭模型的局部穩定點分析和全局漸近性態分析,并通過數值仿真實驗具體說明了企業競爭的各種情況和穩定性條件.
已有文獻的研究是深入且有效的,是對企業企業競爭理論的完善和充實.然而,目前這些研究主要集中于利用二維系統研究企業競爭機制,其不足之處在于只能單純刻畫企業之間的競爭或合作模式,無法揭示企業之間合作并與其他企業競爭的過程.本文將討論企業競爭的更一般情形,以三維系統為研究對象,運用微分方程穩定性理論討論它們在完全合作市場、完全競爭市場和混合市場下競爭的動力學機制,并通過數值仿真實驗呈現企業競爭的動態演化規律.
2 生態競爭模型
假設x1(t),x2(t),x3(t)分別表示企業1、2、3在t時刻的產出水平,這三個量均為時間t的函數.為了分析的方便,這里的時間t不僅包含一般時間的含義,而且還包含技術、信息、成本、專業化和分工等影響產出水平的因素.
Lotka-Volterra競爭模型描述了這三個企業間的競爭:
dx1dt=α1x1(1-x1N1-γ12x2N2-γ13x3N3),
dx2dt=α2x2(1-x2N2-γ21x1N1-γ23x3N3),
dx3dt=α3x3(1-x3N3-γ31x1N1-γ32x2N2), (1)
其中,α1,α2,α3是三個企業的產出水平在獨立生存狀態下的固定增長率,N1、N2、N3是其對應的最大產出水平,γij表示企業i與j之間的競爭系數,這里假設γij=γji.
本文將市場定義為以下三類:
Ⅰ.完全合作市場:當γij<1時,所有的企業之間的競爭相對較弱,每個企業都能取得一定的發展.需要說明的是,這里所謂的“合作”只是意味著較弱的競爭,并不考慮真正的合作現象.
Ⅱ.完全競爭市場:當γij>1時,所有的企業之間存在強的競爭,只有一個企業能在競爭中獲勝,最終取得全部市場.
Ⅲ.混合市場:在這種市場中,兩個企業之間競爭較弱,他們與第三個企業之間存在強的競爭.
下面,本文將分析在這三種市場中企業競爭的動力學機制和動態演化過程.
3 企業競爭動態演化的分析
采用微分方程定性理論[9]分析系統(1)中的平衡點.由定態方程,得到模型的平衡點分別為:P0=(0,0,0);P1=(N1,0,0),P2=(0,N2,0),P3=(0,0,N3),P1、P2、P3表示只有一個企業能獲勝;P12=(N11+γ12,N21+γ12,0),P13=(N11+γ13,0,N31+γ13),P23=(0,N21+γ23,N31+γ23),P12,P13,P23表示兩個“合作”企業將在競爭中獲勝;P123表示三個企業各自占有一份市場的情形.為了計算上的方便,這里假設完全對稱的條件:γ12=γ13=γ23=γ,α1=α2=α3=α,N1=N2=N3=N,P123=(N2γ+1,N2γ+1,N2γ+1).
為了進行平衡點的穩定性分析,計算Jacobian矩陣,
α1-2α1N1y1-α1γ12N2y2-α1γ13N3y3-α1γ12N2y1-α1γ13N3y1-α2γ12N1y2α2-2α2N2y2-α2γ12N1y1-α2γ23N3y3-α2γ23N3y2-α3γ13N1y3-α3γ23N2y3α3-2α3N3y3-α3γ13N1y1-α3γ23N2y2,(2)
其中的y1、y2、y3分別表示平衡點中的對應分量,并計算相應矩陣的特征值,結果如下:
①對于P0=(0,0,0),Jacobian矩陣的特征值為ζ1=α1,ζ2=α2,ζ3=α3,所有的特征值大于0,P0是不穩定的,三個企業都將從零開始發展壯大.
②P1=(N1,0,0),Jacobian矩陣的特征值為ζ1=α3-α3γ13,ζ2=-α1,ζ3=α2-α2γ12.
P2=(0,N2,0),對應的特征值為ζ1=-α2,ζ2=α1-α1γ12,ζ3=α3-α3γ23;
P3=(0,0,N3),對應的特征值為ζ1=α1-α1γ13,ζ2=α2-α2γ23,ζ3=-α3.
以P3為例進行分析.當γ13>1,γ23>1,即企業3與企業1、2同時都存在強的競爭時,P3是穩定的.完全競爭市場競爭模型的仿真實驗(見圖1,
其中α1=α2=α3=1,N1=N2=N3=1,γ12=γ13=γ23=1.5)表明,若初始值的位置靠近P3,則軌線會趨于平衡點P3.因此,在完全競爭市場,競爭的初始實力決定了競爭的最終結果.
圖1 完全競爭市場競爭模型
同時,由于參數γ12沒有影響P3的穩定性,在混合市場(企業1和企業2“合作”), P3也是穩定的.混合市場競爭模型的仿真實驗(見圖2, 其中α1=α2=α3=1,N1=N2=N3=1,γ12=0.5,γ13=γ23=1.5)表明,當初始值接近P3時,軌線也將收斂到平衡點P3.
③對于P12=(N11+γ12,N21+γ12,0),特征值為
ζ1=α3(1+γ12-γ13-γ23)1+γ12,
ζ2=α2+α1-α22-2α1α2+α21+4α2γ212α12(1+γ12),
ζ3=α2+α1+α22-2α1α2+α21+4α2γ212α12(1+γ12).
以P12為例進行說明.當γ12<1,γ13+γ23>1+γ12,即企業1和企業2“合作”,都與企業3有強的競爭時,P12是穩定的.圖2中也表明,當初始值接近P12時,軌線將收斂到平衡點P12.
圖2 混合市場競爭模型
由此可見,初始實力在混合市場的競爭中依然是決定因素.
P13=(N11+γ13,0,N31+γ13),其特征值為
ζ1=α2(1+γ13-γ12-γ23)1+γ13,
ζ2=-α3+α1-α23-2α1α3+α21+4α3γ213α12(1+γ13),
ζ3=-α3+α1+α23-2α1α3+α21+4α3γ213α12(1+γ13);
P23=(0,N21+γ23,N31+γ23),其特征值為
ζ1=α1(1+γ23-γ12-γ13)1+γ23,
ζ2=-α3+α3-α23-2α2α3+α22+4α2γ223α32(1+γ23),
ζ2=-α3+α3+α23-2α2α3+α22+4α2γ223α32(1+γ23).
P13與P23的分析與P12類似.
④對于P123=(N2γ+1,N2γ+1,N2γ+1),特征值為
ζ1=(γ-1)α2γ+1,
ζ2=(γ-1)α2γ+1,
ζ3=-α,
當γ<1,所有企業都處在弱的競爭狀態,即在完全合作市場,P123是穩定的.完全合作市場競爭模型的仿真實驗(見圖3,其中α1=α2=α3=1,N1=N2=N3=1,γ12=γ13=γ23=0.5
)表明,在弱競爭狀態下,三個企業的發展水平將收斂到P123.
圖3完全合作市場競爭模型
4數值仿真實驗
為了更清楚地理解企業之間的競爭,以下將通過數值仿真實驗來模擬企業競爭的動態演化過程.
1)在完全競爭市場,三個企業之間都存在強的競爭,此時競爭將導致實力最強的企業將戰勝其他實力較弱的企業,并最終能達到自己的最大產出水平(見圖4, 其中α1=α2=α3=1,N1=N2=N3=1,γ12=γ13=γ23=1.5,初始條件:x1(0)=0.2,x2(0)=0.3,x3(0)=0.5).
時間t
圖4 企業3將贏得競爭
2)在混合市場中,企業的初始水平將決定到底是“合作”企業獲勝還是與之對抗的企業獲勝.當“合作”企業的初始水平與對抗企業相差不大時,“合作”企業將贏得競爭.同時,由于獲勝的“合作”企業之間競爭較弱,使得雙方的產出能力都不能達到各自最大的水平.當“合作”企業的初始水平遠低于對抗企業時,“合作”企業在競爭中失敗.此時,對抗企業卻能通過發展達到其最大產出水平(見圖5,其中α1=α2=α3=1,N1=N2=N3=1,γ12=0.5,γ13=γ23=1.5,初始條件:x1(0)=0.2,x2(0)=0.3,x3(0)=0.5和圖6, 其中α1=α2=α3=1,N1=N2=N3=1,γ12=0.5,γ13=γ23=1.5,初始條件:x1(0)=0.2,x2(0)=0.3,x3(0)=0.8).在實驗中也發現,勝負的結果與“合作”企業之間的競爭水平關系不大(見圖7,其中α1=α2=α3=1,N1=N2=N3=1,γ12=0.8,γ13=γ23=1.5,初始條件:x1(0)=0.2,x2(0)=0.3,x3(0)=0.5
).
時間t
圖5 “合作”企業1和2贏得競爭
3)在完全合作市場中,由于所有企業之間的競爭較弱,影響了所有企業的產出水平,雖然最終的產出水平都趨于穩定,但遠小于企業所能達到的最大產出水平(見圖8,其中α1=α2=α3=1,N1=N2=N3=1,γ12=γ13=γ23=0.5,初始條件:x1(0)=0.2,x2(0)=0.3,x3(0)=0.5).
5 結 論
本文以三個企業為例分析企業間競爭的動態演化規律,通過微分方程定性理論分析了企業間競爭可能出現的不同情況以及穩定性條件,通過數值仿真實驗動態的展示了競爭的演化過程:在完全競爭市場,實力強的企業將戰勝其他企業;在完全合作市場,所有的企業都能發展,但產出水平不能達到各自
的最大水平;在混合市場,實力也將決定競爭的勝負,當“合作”企業的實力與對抗企業之間相差不大,“合作”企業獲勝,然而,“合作”企業均無法達到各自最大的產出水平,對抗企業的實力若遠大于“合作”企業,“合作”企業將失敗,此時,對抗企業最終能獲得自己的最大產出水平.由于本文的分析是在三個企業的競爭基礎上展開,并將競爭與合作的情況全面考慮,因而,所得到的結果更具有一般意義,是對傳統二維模型的有效擴展.
參考文獻
[1] 陶長琪. IT企業集群的共生性和穩定性研究[J].科技管理研究,2004,(5):62-64.
[2] 沈運紅,王恒山.中小企業與大型企業的互動策略選擇[J].系統工程理論方法與應用,2006,15(5):417-420.
[3] 喻小軍,譚建.企業集群的生態系統競合模型[J].系統工程,2007,25(7):108-111.
[4] 張睿,錢省三,高臻.基于生態學理論的企業競爭模型[J].系統工程,2008,26(2):116-119.
[5] LUIS López, MIGUEL A FSANJUN. Defining strategies to win in the Internet market [J]. Physica A:Statistical Mechanics and its Application, 2001,301(1/4): 512-534.
[6] 王子龍,譚清美,許簫迪.企業集群共生演化模型及實證研究[J].中國管理科學,2006,14(2):141-148.
[7] 周浩. 企業集群的共生模型及穩定性分析[J].系統工程,2003,21(4):32-37.
[8] 李艷會, 朱思銘. 一類電子商務網站競爭模型分析[J]. 中山大學學報:自然科學版,2003,42(5):6-10.
[9] J K HALE. Theory of Functional Differential Equations [M], Springer, NewYork, 1977.
注:本文中所涉及到的圖表、注解、公式等內容請以PDF格式閱讀原文