摘要：基于1978—2016年畜禽產(chǎn)品年度產(chǎn)量數(shù)據(jù)，采用HP濾波技術和非線性MS-AR模型，分析中國畜牧業(yè)生產(chǎn)波動的特征，從畜牧業(yè)生產(chǎn)區(qū)制轉換的視角剖析畜牧業(yè)生產(chǎn)波動的具體路徑，探討影響畜牧業(yè)生產(chǎn)波動和路徑轉換的原因。結果表明：從生產(chǎn)波動變化趨勢上看，豬肉生產(chǎn)增長高峰期出現(xiàn)最早，其次是牛肉、羊肉、禽肉和禽蛋生產(chǎn)，牛奶生產(chǎn)增長高峰期出現(xiàn)時間最晚；從生產(chǎn)波動比較上看，豬肉生產(chǎn)波動幅度最小，其他畜牧產(chǎn)業(yè)生產(chǎn)波動幅度相對較大，各畜牧產(chǎn)業(yè)生產(chǎn)波動趨勢存在差異，但近年來總體趨于收斂。通過模型識別將畜牧業(yè)生產(chǎn)劃分為低速、中速和快速增長3個區(qū)制，低速和中速增長區(qū)制下各畜禽業(yè)產(chǎn)品產(chǎn)量平均增長率在3%和5%左右，而快速增長區(qū)制下平均增長率在10%以上；豬肉生產(chǎn)傾向于保持中速增長，牛肉、羊肉和牛奶生產(chǎn)容易出現(xiàn)低中速增長往復徘徊的現(xiàn)象，禽肉和禽蛋生產(chǎn)在各區(qū)制下均較為穩(wěn)定且區(qū)制變遷較為平穩(wěn)；改革開放以來，除牛奶產(chǎn)業(yè)以外，其他畜牧產(chǎn)業(yè)生產(chǎn)基本在低速或中速增長區(qū)制時間最長。基于此，提出針對不同畜禽產(chǎn)業(yè)在不同發(fā)展階段相機制定不同調控政策，推動畜禽產(chǎn)業(yè)規(guī)模化養(yǎng)殖和產(chǎn)業(yè)化經(jīng)營水平，提高畜禽產(chǎn)業(yè)抗沖擊能力等建議。

　　關鍵詞：畜牧業(yè)；生產(chǎn)波動；路徑識別；HP濾波；MS-AR模型

　　1模型構建與數(shù)據(jù)來源

　　1.1模型構建MS-AR

　　模型由Hamilton[14-15]提出，本質上屬于變參數(shù)模型，使用該模型可以從數(shù)據(jù)中提取過去狀態(tài)的轉換信息，一致而有效的估計參數(shù)，可以有效的捕捉經(jīng)濟變量中非線性動態(tài)特征[16]。考慮K維時間序列經(jīng)濟變量yt=y1,y2,…,yt，t=1,2,…,T，可以構成以下的AR(p)模型：

　　y=v+Ay+∑μ（1）式中：yt為時間序列經(jīng)濟變量，v為截距項，Ai(i=1,2,…,p)為系數(shù)向量，p為滯后階數(shù)μt～NID(0,ε)。（1）式即為具有截距的線性AR(p)模型，將其轉變?yōu)榫敌问降腁R(p)模型，即：-()ptit-iti=1y=Ay-+μμμ∑（2）式中：yt為時間序列經(jīng)濟變量，Ai(i=1,2,…,p)為系數(shù)向量，μt～NID(0,…)，μ表示yt的均值。在上述的AR(p)模型均為線性模型，在AR(p)中引入了馬爾可夫鏈，假定模型中存在著不可觀測的M種狀態(tài)，用St(t=1,2,…,M)來表示，并假定截距項v和誤差項u均具有狀態(tài)轉換特征，St是一組服從馬爾可夫鏈的離散隨機變量，其轉移概率矩陣形式為：111212121212.....................MMMMMMPPPPPPPPPP????=??（3）1)P=P(sjsiijrtt==-（4）1∑Pij=（5）如果截距項v具有狀態(tài)轉換特征，即在（1）式基礎上可以構建變參數(shù)v(St)的MSI(M)-AR(p)模型。依據(jù)模型中自回歸系A1,A2,…,AP，和誤差項μt的異方差ε是否引入?yún)^(qū)制狀態(tài)變量St，可將模型劃分為MSI(M)-AR(p)模型、MSIA(M)-AR(p)模型、MSIH(M)-AR(p)模型和MSIAH(M)-AR(p)模型等。如果均值項μt具有狀態(tài)轉換特征，即在（2）式基礎上可以構建變參數(shù)μ(St)的MSM(M)-AR(p)模型，依據(jù)模型中自回歸系數(shù)A1,A2,…,AP和誤差項μt的異方差ε是否引入?yún)^(qū)制狀態(tài)變量St，可將模型劃分為MSM(M)-AR(p)模型、MSMA(M)-AR(p)模型、MSMH(M)-AR(p)模型和MSMAH(M)-AR(p)模型等。MS(M)-AR(p)模型參數(shù)可通過Hamilton[15]提出的期望最大化（expectationmaximization，EM）算法和Krolzig[17]提出的極大似然估計（maximumlikelihood，ML）得到。詳細的推導過程可參考Krolzig[17]的研究，在此不再贅述。具體模型形式的選擇上，可根據(jù)AIC、SC和HQ值最小化原則確定最優(yōu)模型形式。

　　1.2數(shù)據(jù)來源

　　本研究涵蓋主要畜禽產(chǎn)品，具體包括豬肉、牛肉、羊肉、禽肉、禽蛋和牛奶等6大類畜禽產(chǎn)品，研究數(shù)據(jù)均為各畜禽產(chǎn)品產(chǎn)量數(shù)據(jù)，樣本區(qū)間為1978—2016年，其中1978—1999年產(chǎn)量數(shù)據(jù)來源于聯(lián)合國糧農組織統(tǒng)計數(shù)據(jù)庫（FAOSTAT），2000—2016年產(chǎn)量數(shù)據(jù)來源于農業(yè)部畜牧業(yè)司、全國畜牧總站編撰的《中國畜牧業(yè)統(tǒng)計》。

　　2結果與分析

　　2.1中國畜牧業(yè)生產(chǎn)波動的基本特征

　　改革開放以來，中國畜牧業(yè)各類產(chǎn)品產(chǎn)量增長率的動態(tài)變動軌跡（圖1）。首先，可以看出改革開放以來各畜禽產(chǎn)品產(chǎn)量增長率波動起伏，不同時期畜禽產(chǎn)品產(chǎn)量增長率存在一定差異；其次，可將各畜禽產(chǎn)品產(chǎn)量增長率大致劃分為兩個階段，一是1979—2000年，該時期畜禽產(chǎn)品產(chǎn)量增長率相對較高，增長率波動幅度相對較大，二是2000—2016年，該時期畜禽產(chǎn)品產(chǎn)量增長率相對較低，增長率波動幅度相對較小；最后，通過觀察2012—2016年畜禽產(chǎn)品產(chǎn)量增長率可以看出，多個畜禽產(chǎn)品產(chǎn)量增長率出現(xiàn)明顯的負向增長態(tài)勢。為了更清晰直觀地剖析中國畜牧業(yè)不同畜禽產(chǎn)品產(chǎn)量增長率的演變軌跡，下文將基于HP濾波技術和MS-AR模型對各畜禽產(chǎn)品產(chǎn)量增長率進行具體分析。

　　2.2中國畜牧業(yè)生產(chǎn)波動的演變軌跡

　　通過HP濾波技術得到不同畜禽產(chǎn)品產(chǎn)量增長率的趨勢成分和波動成分變動軌跡（圖2）。改革開放以來豬肉產(chǎn)量增長率總體呈不斷下降態(tài)勢，尤其是2010年以來，下滑趨勢較為明顯；從波動成分上看，豬肉產(chǎn)量增長率改革開放初期波動幅度較大，自90年代后期開始波動幅度有所收斂，但2012年以來再次表現(xiàn)出較大的波動態(tài)勢。牛肉產(chǎn)量增長率經(jīng)歷了先上升后下降的過程，在90年代初期達到頂峰，之后開始下降，2010年基本達到谷底并趨于平穩(wěn)；從波動成分上看，牛肉產(chǎn)量增長率在改革開放后出現(xiàn)了較大的波動起伏現(xiàn)象，但隨著時間推移波動幅度呈下降態(tài)勢，2012年以來波動基本趨于0。對于羊肉產(chǎn)量增長率趨勢成分變化趨勢而言，改革開放后，羊肉產(chǎn)量增長率緩慢下降，之后開始緩慢回升，但自90年代中期開始出現(xiàn)較大幅度下滑，近幾年有緩慢上升趨勢；從波動成分上看，同豬肉和牛肉產(chǎn)量增長率波動相似，前期波動幅度較大，之后波動開始收斂。從禽肉產(chǎn)量增長率趨勢成分變化趨勢上看，其變化趨勢與牛肉產(chǎn)量增長率趨勢成分變化趨勢相似，經(jīng)歷了先上升后下降的過程，在90年代初期達到頂峰，之后開始下降，2005年之后基本趨于平穩(wěn)；從波動成分上看，禽肉產(chǎn)量增長率波動經(jīng)歷了先上升后下降的過程，近年來，除個別年份由于禽流感等疫情影響出現(xiàn)大幅波動外，波動幅度均相對較小。禽蛋產(chǎn)量增長率趨勢成分經(jīng)歷了先緩慢上升又快速下降的過程，進入新世紀后基本保持較為平穩(wěn)；從波動成分上看，改革開放后波動幅度較大，但自90年代后期逐漸收斂，進入新世紀以來，除個別年份外，基本呈現(xiàn)小幅波動態(tài)勢。牛奶產(chǎn)量增長率趨勢成分變化趨勢不同于其他畜禽產(chǎn)品產(chǎn)量增長率變化趨勢，其經(jīng)歷了先緩慢上升、又緩慢下降之后又快速上升、再快速下降的“M”型變化過程，增長率高峰期與其他畜禽產(chǎn)品相比，出現(xiàn)相對滯后，在2004年左右達到峰值，之后出現(xiàn)了“過山車式”快速下降；從波動成分上看，在改革開放后的前20年，基本呈現(xiàn)小幅波動態(tài)勢，進入新世紀后波動幅度開始增大。

　　3結論與政策啟示

　　3.1結論

　　研究表明，從產(chǎn)量增長率變化趨勢上看，豬肉產(chǎn)量增長率高峰期出現(xiàn)最早，其次是牛肉、羊肉、禽肉和禽蛋產(chǎn)量增長率，而牛奶產(chǎn)量增長率高峰期出現(xiàn)時間最晚；從不同畜禽產(chǎn)品產(chǎn)量波動比較上看，豬肉產(chǎn)量增長率波動幅度最小，但近幾年有增大的趨勢；從產(chǎn)量增長率的總趨勢來看，近幾年，牛肉產(chǎn)量增長率最為平穩(wěn)，其次為羊肉和禽蛋產(chǎn)量增長率，雞肉產(chǎn)量增長率在最近一些年份也出現(xiàn)較大波動，而牛奶產(chǎn)量增長率處于波動調整期，并未趨于平穩(wěn)。通過模型識別將畜禽產(chǎn)品生產(chǎn)劃分為“低速增長區(qū)制”、“中速增長區(qū)制”和“快速增長區(qū)制”3個區(qū)制，其中低速增長區(qū)制下畜禽產(chǎn)品產(chǎn)量平均增長率基本在3%左右，中速增長區(qū)制中平均增長率在5%左右，快速增長區(qū)制中各畜禽產(chǎn)品產(chǎn)量平均增長率基本在10%以上；豬肉生產(chǎn)傾向于保持中速增長，牛肉、羊肉和牛奶生產(chǎn)容易出現(xiàn)低速增長和中速增長往復徘徊的現(xiàn)象，禽肉和禽蛋產(chǎn)量增長率在各區(qū)制下均較為穩(wěn)定。改革開放以來，除牛奶以外，其他畜牧產(chǎn)業(yè)基本在低速或中速增長區(qū)制時間最長，從平均持續(xù)期來看，除禽蛋產(chǎn)業(yè)以外，其他畜牧產(chǎn)業(yè)均在快速增長區(qū)制下平均持續(xù)期最長。從區(qū)制變遷上看，豬肉、牛肉、羊肉和牛奶產(chǎn)量增長率區(qū)制變遷較為頻繁，而禽肉和禽蛋產(chǎn)量增長率區(qū)制變遷較為平穩(wěn)。造成畜禽產(chǎn)品生產(chǎn)區(qū)制轉移的原因主要有供需差異、價格波動、進口沖擊和畜禽疫情等，除此之外，規(guī)模化程度和政府調控政策也會對其造成影響。預計未來一段時期，豬肉、牛肉、禽蛋和牛奶生產(chǎn)將以低速增長為主，羊肉生產(chǎn)主要處于中速增長區(qū)制，禽肉生產(chǎn)仍將處于中低速發(fā)展狀態(tài)。

　　3.2政策啟示

　　首先，針對不同畜禽產(chǎn)業(yè)發(fā)展特征，制定不同產(chǎn)業(yè)調控政策，穩(wěn)定重要畜禽產(chǎn)品生產(chǎn)，加大對弱勢畜禽產(chǎn)業(yè)扶持，因地制宜推動產(chǎn)業(yè)平穩(wěn)發(fā)展。其次，當產(chǎn)業(yè)發(fā)展進入快速發(fā)展階段時，應提高產(chǎn)業(yè)抗風險和外部沖擊能力，一方面健全和完善獸醫(yī)防疫體系，提高畜禽產(chǎn)業(yè)疫病防控能力，加大重大動物疫病的防控，強化動物防疫監(jiān)督執(zhí)法，維護畜禽養(yǎng)殖業(yè)安全發(fā)展；另一方面提高產(chǎn)業(yè)進口及價格波動監(jiān)測預警，防止外部沖擊對產(chǎn)業(yè)造成大幅波動。

　　參考文獻：

　　[1]張新煥,閻新華,楊德剛,等.中國畜牧業(yè)發(fā)展的時空演變及趨勢分析[J].干旱區(qū)地理,2003,26(3):281-285.ZhangXH,YanXH,YangDG,etal.AnalysisonspatialandtemporalevolutionandtrendofanimalhusbandryinChina[J].AridLandGeography,2003,26(3):281-285.

　　[2]徐雪高,陳潔,李靖.中國畜牧業(yè)發(fā)展的歷程與特征[J].中國畜牧雜志,2011,47(20):14-17.XuXG,ChenJ,LiJ.HistoryandcharacteristicsofChina’slivestock[J].ChineseJournalofAnimalScience,2011,47(20):14-17.

　　[3]楊子剛,毛文坤,于海朋.中國畜牧業(yè)生產(chǎn)區(qū)域布局演變分析[J].中國畜牧雜志,2012,48(10):11-14.YangZG,MaoWK,YuHP.AnalysisonevolutionofanimalhusbandryproductionregionaldistributioninChina[J].ChineseJournalofAnimalScience,2012,48(10):11-14.

　　[4]隋建利,蔡琪瑤.中國農業(yè)經(jīng)濟周期的路徑演化識別——改革開放以來的實踐與經(jīng)驗[J].中國農村經(jīng)濟,2016(9):30-43.SuiJL,CaiQY.IdentificationoftheagriculturaleconomiccycleinChina:Asummarysincethereformandopening[J].ChineseRuralEconomy,2016(9):30-43.