供應鏈倉儲物流金融服務定價模型研究*

汪 洋,邱紅星

(華南理工大學 工商管理學院,廣東 廣州 510640)

一、引言

港口群作為交通運輸樞紐協調經濟腹地和國際貿易物流,是區域經濟乃至國民經濟發展的重要支撐力量。全球金融風暴引發港口群貨物吞吐量和集裝箱吞吐量的增長整體下降,直接導致倉儲物流積滯,供應鏈系統物流形成大量倉儲,企業流動資金被占用以至出現資金短缺。[1]倉儲物流金融作為商業銀行新型融資業務,能夠有效緩解供應鏈系統物流融資倉儲的倉儲物流積壓與資金短缺困境。[2]

物流金融的過往研究主要集中在其概念、模式及系統風險方面。物流金融是商業銀行為降低交易成本與風險,攜同物流企業依據供應鏈系統物流管理而進行的金融活動。[3]有學者運用層次模糊綜合評價法構建基于三方經營主體的定量與定性指標相結合的風險評價模型,為商業銀行物流金融倉儲融資服務提供有益的決策參考。[4]有關研究從供應鏈系統的角度探討物流金融倉儲融資業務,通過構建倉儲融資物流資金流量結構邏輯模型及風險/收益優化模型,論證物流倉儲融資能促進大型交通樞紐發展并使商業銀行在較低風險下獲得較高收益；[5-6]運用博弈論方法論證倉儲融資各參與主體融資行為的選擇是博弈的必然結果。[7]

本文試圖構建供應鏈倉儲物流金融服務定價模型,采用非線性動力學方法分析其定價機制及價格穩定性,科學地制訂倉儲物流金融服務價格,協調雙方利益關系,將商業銀行倉儲物流金融服務與供應鏈系統物流倉儲業務有機結合,拓展融資空間,盤活倉儲物流積壓資金,促進國民經濟發展。

二、研究邊界

供應鏈系統是供應商、制造商、分銷商和零售商等多個物流環節構成,由物流、信息流和資金流連接起來的系統,具有開放性、動態性、網絡復雜性及各個環節交互推動作用等特點,具有開放型的復雜系統特征,屬于復雜性系統的研究范疇。[8]供應鏈倉儲物流金融服務,是商業銀行以供應鏈倉儲物流質押為授信條件,利用倉儲物流所有權融資而提供的綜合金融服務。[9]

供應鏈倉儲物流金融服務定價的研究范圍,主要是探討供應鏈系統中的倉儲物流,商業銀行金融服務的倉儲物流價格制定機制和定價模型、價格穩定性及其條件等。研究主體涉及商業銀行、倉儲監管方和供應鏈系統倉儲物流所有權者。[10]倉儲物流金融服務定價的影響變量主要探討倉儲物流金融服務業務成本(C)、業務利潤率(RP)、業務需求量(Q)等。

三、供應鏈倉儲物流金融服務定價機制探討

供應鏈系統倉儲物流金融服務定價主要由商業銀行、供應鏈系統倉儲物流所有權者及倉儲監管方三個業務主體的影響因素,決定其間特定的利益與權責關系。供應鏈系統倉儲物流所有權者以其一定價值的動產質押,獲得商業銀行貸款認可,兩者間形成借貸關系與質押關系。[11]出于降低成本與專業化管理考慮,商業銀行將其對質押品的監管外包,與倉儲監管方形成委托代理關系。倉儲監管方與倉儲物流所有權者形成監管關系。

(一)倉儲物流金融服務定價因素分析

供應鏈倉儲物流金融服務系統的三方參與主體的協同關系可以考慮為： 倉儲物流所有權者是倉儲物流金融服務的受用方,支付倉儲物流金融服務費用,即倉儲物流金融服務價格的接受方；商業銀行與倉儲監管方是倉儲物流金融服務的提供方,收取倉儲物流金融服務中間業務費,是服務價格的制定方。商業銀行在倉儲物流金融服務中處于中心地位,有專業的貸款定價資源,主導倉儲物流金融服務價格的制定。倉儲監管方提供定價過程所需業務倉儲信息及倉儲物流數據,為倉儲物流金融服務價格的接受方和制定方,協商價格乃至最后確定價格,提供管理服務并從中收取服務費。

供應鏈倉儲物流金融服務定價的運作流程和定價因素如圖1所示。商業銀行業務審核部門對物流融資方(倉儲物流所有權者)與倉儲監管方(倉儲物流監管方)進行資信審查和信息采集,主要包括歷史倉儲信息、資信記錄、質押品所有權與價值等,倉儲物流融資信息及數據輸入商業銀行數據處理中心,其業務拓展部門根據物流融資方及倉儲監管方的信息和數據資料核定適宜的業務種類。

圖1 供應鏈倉儲物流金融服務定價因素

商業銀行業務風險分析系統根據倉儲物流金融業務種類與數據,確定業務的法律風險系數、信用風險系數等,將數據輸入數據處理中心。經過審核并確定業務類型后,商業銀行分析具體的因素影響及程度,并錄入商業銀行倉儲物流金融服務定價數據庫,繼而進入數據處理中心。倉儲物流金融服務定價的影響因素主要有物流金融影響因素、物流倉儲影響因素、物流財務影響因素及宏觀經濟影響因素。這四大影響因素可細分為二級影響因素。物流財務影響因素可細分為倉儲物流總價值、利率、倉儲物流融資量、應付賬款等；物流倉儲影響因素的二級影響因素為倉儲監管成本、目標利潤率、人力資源成本等；物流金融影響因素的二級影響因素為營銷成本、貸款折扣率、融資利率等；宏觀經濟影響因素的二級影響因素為倉儲物流金融服務的競爭價格、通貨膨脹率等。倉儲物流金融服務定價分析系統分析價格與定價因素之間的關系,制定系統倉儲物流金融服務參考價格,通過信息管理中心,知曉商業銀行和倉儲物流所有權者。定價主體雙方將審定的參考價格反饋倉儲物流金融服務定價模型,最終制定并輸出倉儲物流金融服務價格。

(二)倉儲物流金融服務定價機制

供應鏈倉儲物流金融服務定價機制(圖2),主要由三大模塊構成： 分析模塊、數據模塊與定價模塊。分析模塊研究、分析倉儲物流金融服務定價內部影響因素與外部影響因素,為倉儲物流金融服務價格的最終確定提供決策支持。

定價模塊主要完成定價數理模型構建。數據模塊收集倉儲物流金融服務系統內、外影響因素的信息與數據,信息與數據經處理后傳輸至定價模塊和分析模塊,為其提供數據支持。倉儲物流金融服務參考價格與最終后定價格傳給三方參與主體,最終確定供應鏈倉儲物流金融服務價格。

圖2 供應鏈倉儲物流金融服務定價機制

四、供應鏈倉儲物流金融服務定價模型研究

商業銀行各項業務均需滿足盈利性、安全性及流動性的要求。依據這三個原則,商業銀行融資服務定價的影響因素主要有貸款成本、信貸風險、貸款額度與期限、客戶關系、商業銀行收益目標、借貸市場競爭程度等。[12]

(四)倉儲物流金融服務定價模型變量

供應鏈倉儲物流金融服務定價的影響變量可分為三類： 倉儲物流金融服務定價的成本因素變量、利潤率因素變量、需求量因素變量。供應鏈倉儲物流金融服務定價的影響變量及變量關系如表1所示。

表1 供應鏈倉儲物流金融服務定價的影響變量及關系

倉儲物流金融服務的成本因素變量主要包括： 倉儲物流金融服務的資金成本、倉儲監管成本、人力資源成本、營銷成本；倉儲物流金融服務的利潤率因素變量主要包括： 倉儲物流金融服務目標利潤率、風險補償、通貨膨脹率、商業銀行其他業務利潤率；倉儲物流金融服務的需求量因素變量主要包括： 倉儲物流總價值、利潤額、倉儲物流金融服務的競爭價格。

(二)供應鏈倉儲物流金融服務定價模型變量的關系分析

倉儲物流金融服務定價影響變量,通過倉儲物流金融服務的成本、利潤率及需求量三大類要素的變動,影響倉儲物流金融服務價格的穩定性。供應鏈系統倉儲物流金融服務定價變量關系如圖3所示。

圖3 供應鏈倉儲物流金融服務價格影響變量關系圖

1. 倉儲物流金融服務定價成本(C)

倉儲物流金融服務定價最基本的要求是利率(價格)可以覆蓋業務開展產生的成本和合理的利潤。成本與倉儲物流金融服務的價格成正向關系,成本越高,價格越高。倉儲物流金融服務的業務成本來自商業銀行與倉儲監管方兩個參與主體,其影響變量主要包括資金成本,倉儲監管成本,人力資源成本,營銷成本。

倉儲物流金融服務成本因素(C)受資金成本(F)、倉儲監管成本(W)、人力資源成本(H)及銷售成本(M)的影響。倉儲物流金融業務中,商業銀行讓渡資金的使用權,存在資金成本(F)。倉儲監管方為質押貨物提供倉儲監管服務,倉庫租賃與監管費用構成倉儲監管成本(W)。商業銀行對倉儲物流所有權者授信調查,業務審核,倉單真實性及質押物狀況跟蹤調查等,其人力費用為其人力資源成本(H)。倉儲物流金融業務宣傳推廣,構成銷售成本(M)。這四類變量與倉儲物流金融服務成本均為正相關。

2. 倉儲物流金融服務定價的利潤率(RP)

貸款是商業銀行資金運用的主要業務,貸款收益率直接影響到商業銀行總體盈利目標的實現。[13]因此,倉儲物流金融服務利潤率成為倉儲物流金融服務定價的重要影響因素。倉儲物流金融業務利潤率與其價格正相關,利潤率越高,價格越高。倉儲物流金融服務利潤率的確定會考慮管理層要求的目標利潤率,同時業務開展所面臨風險、宏觀經濟變量以及商業銀行其他業務利潤率也會對倉儲物流金融服務利潤率的確定產生影響。

倉儲物流金融服務利潤率(RP)主要受目標利潤率(RT)、通貨膨脹率(RI)、其他業務利潤率(RO)、風險補償(CR)的影響。商業銀行倉儲物流金融服務存在業務風險,要求一定的風險補償(CR)。通貨膨脹會導致倉儲物流所有權者償還貸款金額的實際價值減小,繼而倉儲物流金融業務的利潤減少。目標利潤率(RT)、通貨膨脹率(RI)、風險補償(CR)與倉儲物流金融服務的利潤率(RP)均為正相關。商業銀行其它業務利潤率(RO)反映倉儲物流金融服務的機會成本,對倉儲物流金融服務利潤率產生一定影響。

3. 倉儲物流金融服務定價需求量(Q)

現階段,我國物流公司在法律上不具備開展物流金融服務的資格,物流金融服務只由少數商業銀行開展,市場存在較大的進入障礙,可將物流金融服務的市場看作寡頭壟斷的市場類型。寡頭壟斷市場介于完全競爭與完全壟斷市場之間,市場價格確定不僅需要考慮供求關系,同時還要考慮競爭對手的行為。在本文中,將競爭對手行為對價格確定的影響歸入其對服務需求的影響,通過服務需求影響價格。因而本文可以直接考慮供求關系對價格的確定。由供求原理可知,供應鏈倉儲物流金融服務價格是與倉儲物流融資需求成正相關的,物流融資需求越大,倉儲物流金融服務價格越高。

倉儲物流金融服務需求量(Q)受到融資方倉儲物流總價值(S)、利潤額(E)和倉儲物流金融服務競爭價格(Pc)的影響。融資方倉儲物流價值絕對額(S),會影響倉儲物流金融服務的融資需求,倉儲物流絕對額越大,物流所有權者被占用資金越多,融資需求越強。倉儲物流總價值(S)與倉儲物流金融服務的需求量(Q)正相關。倉儲物流所有權者利潤額(E)越高,其以所得利潤填補被占用資金的意愿越強。因而,倉儲物流所有權者利潤額(E)與需求量(Q)負相關。商業銀行競爭對手倉儲物流金融服務的價格(PC)會對業務的需求量(Q)產生一定的影響。

(三)供應鏈倉儲物流金融服務定價模型探討

倉儲物流金融是金融業與物流業相互合作而衍生的新型金融業務,屬于商業銀行貸款業務。商業銀行貸款定價是貸款利率的確定,倉儲物流金融服務定價是確定倉儲物流所有權者向商業銀行融資所需要支付的利率。商業銀行傳統的定價模型主要有成本加成定價模型、價格先導定價模型、客戶盈利分析模式等,[14]根據貸款業務的具體特點,選取適合的定價模型制定價格?；谏虡I銀行全面管理理念與對優勢客戶的爭奪,單一類型的定價方法不適用于倉儲物流金融服務定價。本文擬構建供應鏈倉儲物流金融服務的綜合貸款定價模型。

假設供應鏈倉儲物流所有權者根據其資金需求對價值為SF的倉儲物流融資,商業銀行根據其融資倉儲物流價值(SF)確定貸款折扣率(ρ),向倉儲物流所有權者發放貸款(Q)。同時,商業銀行根據倉儲物流金融服務的成本,利潤率及需求量的分析制定倉儲物流金融服務價格(P)。

在倉儲物流融資發生的時點,商業銀行倉儲物流金融服務的成本率與利潤率為其價格的直接構成因素,為線性疊加關系,同時假設需求量因素即融資額(Q)對價格(P)產生π的影響,則商業銀行倉儲物流金融服務的時點價格為：P=C+RP+πQ；式中,π為物流融資額對倉儲物流金融服務價格的影響系數。

1. 成本模型

倉儲物流金融服務的資金成本(F)、銷售成本(M)、人力資源成本(H)、倉儲監管成本(W)構成了倉儲物流金融服務成本(C)。當融資倉儲物流量變化時,倉儲監管單位成本會隨著融資倉儲物流量的增大而遞減(絕對值增加),銷售成本與人力資源成本變化不大,因而可假設銷售成本與人力資源成本是固定的,則商業銀行倉儲物流金融服務的時點成本可表示為：

C=(M+W+H)/Q+F

W=μSF；式中,μ為倉儲費用率。

2. 利潤率模型

倉儲物流金融服務利潤率的主要構成成分是目標利潤率與風險補償率,通貨膨脹率與其他業務利潤率亦會對其產生一定程度影響。具體關系為：RP=RT+CR+αRI+βRO

式中,α和β分別為通貨膨脹率與其他業務利潤率對倉儲物流金融服務利潤率的影響系數。由于商業銀行的貸款資金來源一般可分解為兩個方面： 舉債和產權資本,因此商業銀行可以根據既定的產權資本目標(RZ),客戶貸款額度(Q)及貸款的資本金支持率(RC)來確定時點目標利潤率,[12]因此：RT=RZ·Q·RC

3. 需求量模型

倉儲物流所有權者利潤額,倉儲物流金融服務競爭價格,均會影響倉儲物流所有權者倉儲融資的意愿；商業銀行根據倉儲物流所有權者的融資倉儲物流以貸款折扣率ρ發放貸款,倉儲物流所有權者的融資倉儲物流時點價值(SF)及時點融資額度(Q)為：

SF=S·(ω1·E+ω2·PC)；

Q=SF·ρ；式中,ω1、ω2分別為倉儲物流所有權者利潤額,倉儲物流金融服務競爭價格對倉儲物流所有權者融資意愿的影響系數。

綜合以上各式,可得商業銀行倉儲物流金融服務時點定價公式為：

P=[(M+H)/φ]·S-1+(RZ·RC·φ+π·φ)·S+ζ；

上式中,φ=ρ(ω1·E+ω2·PC)；

ζ=(CR+αRI+βRO)+(μ/ρ+F) 。

商業銀行采集倉儲物流金融服務影響變量信息及數據后,將數據輸入上述定價模型,即可求得倉儲物流金融服務時點價格。

隨著時間的推移,物流金融服務時點定價公式中,H、F、PC、S等都隨時間變化,價格也在波動的。下面對倉儲物流金融服務價格穩定性進行討論。

因為林露白跟家人住，加上兩人都剛畢業，年紀也還小，沒想過結婚的事。直到戀愛第三年的時候，才跟家里講，雙方見了父母，彼此都很滿意，初步定下了結婚的日期。但是，林露白的父母都是教師，有些傳統的觀念，對她唯一的要求是不允許夜不歸宿，她只能周末找個借口出去逛街的時候，跟魏舟去酒店開個房間溫存一會兒，之后再匆忙地趕回家。

(四)倉儲物流金融服務價格穩定性分析

供應鏈倉儲物流金融服務系統是一個復雜的動態開放性系統。系統內部的物流、信息流與資金流在不斷地發生變化。因而影響供應鏈倉儲物流金融服務的物流微觀變量和宏觀變量隨時間變化,供應鏈物流動態系統的倉儲物流金融服務因素影響具有時間延遲、信息放大、動態性和非線性的特征。[15]在時間t0到t1序列內,倉儲物流金融服務定價影響變量在變動,商業銀行倉儲物流金融服務價格也在波動,用x表示隨時間變化的倉儲物流金融服務的價格。從物流金融服務時點定價公式進行分析,主要是融資倉儲物流量和倉儲物流金融服務成本這兩類的宏觀變量與微觀變量隨時間在變化,用y和z表示倉儲物流所有權者融資倉儲物流量的變化程度和倉儲物流金融服務成本的變化程度。

同理,貨幣供應量M2,CPI,PPI,工資率,倉儲租賃費率,存款利率都會隨時間變化,并導致倉儲物流金融服務成本的變動。例如,其他影響變量不變,M2的增加導致市場流動性增加,資金借貸成本下降,倉儲物流金融服務成本亦相應下降。在t0-t1,變量波動方向具有不確定性,倉儲物流金融服務成本z的變化可看為其影響的非線性集成。有：

供應鏈倉儲物流金融服務價格、成本與融資倉儲物流量的影響變量隨時間變化的不確定性,導致倉儲物流金融服務價格、成本與融資倉儲物流量呈現復雜的非線性波動,且三者相互影響的變化速率具有非線性變化關系。[16]宏觀經濟中匯率等變量波動將影響融資倉儲物流量y的變化,導致價格的波動；當融資倉儲物流量增加時,倉儲物流金融服務的需求增長,其價格上漲趨勢增強。倉儲物流金融服務價格與其成本成線性關系,因此有：

(1)

倉儲物流金融服務的價格變化反作用于融資倉儲物流量。服務價格上升,倉儲物流所有權者考慮從商業銀行的競爭對手獲取貸款,或以發行債券、股票等方式融資,倉儲物流金融服務價格的變化將導致融資倉儲物流量加速的變化。倉儲物流金融服務成本上升,會使倉儲物流金融服務價格上升,最終致使融資倉儲物流量加速下降,經濟增長速度σ與融資倉儲物流量成負向線性關系：

(2)

(3)

聯立(1)、(2)、(3)三式,供應鏈倉儲物流金融服務價格波動的非線性特征為：

令Df(θ)=0,分析此非線性系統在平衡點θ的穩定性,有：

Df(θ)-λE=

可得其特征方程為：

-λ3+[(m-an)-σ]λ2+

(4)

當以上特征方程的特征值均具有負實部時,平衡點θ是穩定的[17]。本文不直接對此方程組求解。根據羅斯-霍維茲判據,此三維非線性系統定點漸近穩定的充要條件是[18]：

由于Δ1=(m-an)-σ>0,因而可消去Δ2>0中因子(m-an)-σ。又有Δ3=Δ1·Δ2,因此不等式方程組只要滿足Δ1>0和Δ2>0即可。故上述不等式可簡化為：

由于倉儲物流金融服務為新拓展業務,市場需求旺盛,商業銀行具有較高的定價自主權,如業務需求出現波動,即調整其倉儲物流金融服務價格。因而,a>1,故有： 2ab+(2a-1)(m-an)>0。再由b+m-an>0,則上式可進一步化簡化為：

則σ的取值范圍：

由以上分析可知,供應鏈倉儲物流金融服務價格的穩定性與經濟發展速度有密切的關系。當經濟發展速度維持在一定的水平,即當：

五、結論

供應鏈倉儲物流金融服務定價系統是多層次的非線性復雜系統。其定價機制受倉儲物流金融服務定價的業務成本因素、業務利潤率因素和業務需求量因素的影響。這些變量參數在供應鏈系統中不斷發生變化,變量參數的廣泛性與波動性使定價系統形成具有自組織功能的復雜網絡結構。

通過供應鏈系統倉儲物流金融服務定價的影響因素： 成本因素、利潤率因素、市場需求因素及其數據變量的影響程度分析,構建倉儲物流金融服務時點定價穩定性模型,以供倉儲物流融資參與主體對金融服務定價的確認,商業銀行具有較高的定價自主權。

由于國民經濟宏觀變量依時間序列變動,供應鏈倉儲物流金融服務定價影響變量數據亦隨時序變化,使倉儲物流金融服務定價處于動態變化的開系統中。運用非線性動力學方法,構建供應鏈倉儲物流金融服務價格動態模型；由模型推導可知,倉儲物流金融服務價格與經濟發展速度存在密切關系。當經濟發展速度維持在某區間時,倉儲物流金融服務價格保持穩定狀態；若經濟發展速度出現異常(危機/泡沫),系統倉儲物流金融服務定價將出現分岔的不穩定態。

本課題還有待深入研究： 倉儲物流金融服務參與各方利益分配的定量分析,金融工程技術在倉儲物流融資管理上的應用,倉儲物流金融服務參與各方所面臨風險的定量分析,倉儲物流金融服務信用與評估體系問題等。

參考文獻：

[1] WANG Yang,MA Xingrui,CHEN Yuan zhi. Study on the Harmonious Operation of Logistics Finance Management of Supply Chain in Ports-Hinterland Based on System Dynamics [M]//I CRMEM 2008. The 2nd International Conference on Risk Management and Engineering Management, Beijing, 2008:82-86.

[2] 汪洋,黃雯雯. 港口群/腹地區域供應鏈系統物流財務管理的國民經濟動態性分析[J]. 工業工程,2010(01)： 36-41.

[3] 宋焱,李偉杰. 物流金融： 出現動因、發展模式與風險防范[J]. 南方金融,2009(12)： 41-46.

[4] 李電生,員麗芬. 基于物流金融模式的中小型企業融資風險評價研究[J]. 技術經濟與管理研究,2010(01)： 62-64.

[5] 汪洋,馬興銳. 港口群/腹地區域供應鏈系統物流倉儲融資的動態性分析[J]. 現代財經,2010(2)： 71-75.

[6] 汪洋,馬興銳. 金融危機態勢下的港口群/腹地區域物流系統倉儲融資機理研究[J]. 工業技術經濟,2009(11)： 96-101.

[7] Wang Yang,Ma Xingrui. Study on Warehousing Financing Risk-return of Supply Chain System in Chinese Ports-Hinterland [M]//ICIM2009 . Proceedings of the 5THInternational Conference on Innovation and Management. Wuhan： Published by Wuhan University of Technology Press,2009： 1551-1559.

[8] WANG Yang,PENG Yulian,LUO Yuanjing. Study on Complexity and Non-linear of Logistics Finance Management of Supply Chain System Based on System Dynamics [M]//ICME2008 . Proceedings of 2008 International Conference on Management Science & Engineering. Long Beach, USA： Copyright,2008IEEE,2008： 1203-1213.

[9] WANG Yang,YU Donglin,LUO Yuanjing,HAO Hongyu. Study on Features of Logistics Finance of Supply Chain Systemand Pledge/Factoring Model[EB/OL].(2011-09-11) [2011-09-307]http： //ieeexplore. ieee. org/stamp/stamp. jsp?tp=&arnumber=4679363.

[10] Wang Yang,Ma Xingrui. Economic Analysis on Ports-Hinterland Supply Chain System Warehousing Financing [M]//ECBI2009. Electronic Commerce Business Intelligence. Tokyo： Published by 2009 IEEE(CPS),2009： 202-206.

[11] 陳曉蕾. CS公司倉儲物流融資研究[D]. 廣州：華南理工大學工商管理學院, 2009.

[12] 舒寧. 下限加點： 商業銀行貸款利率的定價模型[J]. 武漢金融,2009(2)： 63-64.

[13] 吳占權. 新型農村金融機構的貸款定價問題探討[J]. 農村經濟,2009(10)： 69-72.

[14] 陳忠. 商業銀行貸款定價理論與實踐[J]. 金融理論與實踐,2007(11)： 28-31.

[15] Wang Yang,Peng Yulian,Ma Xingrui. A Dynamics Research on Ports-Hinterland Supply Chain System Warehousing Financing Based on Financial Logistics [M]//ICIM2008 . Proceedings of the 5THInternational Conference on Innovation and Management. Wuhan： Published by Wuhan University of Technology Press,2008： 1402-1407.

[16] 汪洋,羅遠菁,彭玉蓮. 供應鏈系統物流財務管理的復雜性和非線性研究[J]. 現代財經,2008(4)： 59-63.

[17] 楊萬利,王鐵寧. 非線性動力學理論方法及應用[M]. 北京： 國防工業出版社,2007.

[18] 劉秉正,彭建華. 非線性動力學[M]. 北京： 高等教育出版社,2004.