岳冬冬，王魯民，紀(jì)煒煒，阮 雯，王 茜，熊敏思，肖 黎，鄭 亮

(1. 農(nóng)業(yè)部遠洋與極地漁業(yè)創(chuàng)新重點實驗室，中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090; 2. 農(nóng)業(yè)部東海漁業(yè)資源開發(fā)利用重點實驗室，中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090)

基于協(xié)整理論和VAR模型的中國農(nóng)村居民人均純收入與水產(chǎn)品消費關(guān)系研究

岳冬冬1，王魯民2，紀(jì)煒煒1，阮 雯2，王 茜1，熊敏思2，肖 黎1，鄭 亮2

(1. 農(nóng)業(yè)部遠洋與極地漁業(yè)創(chuàng)新重點實驗室，中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090; 2. 農(nóng)業(yè)部東海漁業(yè)資源開發(fā)利用重點實驗室，中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090)

借助2000～2015年相關(guān)統(tǒng)計指標(biāo)數(shù)據(jù)，應(yīng)用協(xié)整理論和向量自回歸模型等計量經(jīng)濟方法，探討了農(nóng)村居民人均純收入與人均水產(chǎn)品消費量之間的動態(tài)變化關(guān)系。研究結(jié)果表明：農(nóng)村居民人均純收入與人均水產(chǎn)品消費量之間存在長期協(xié)整關(guān)系；農(nóng)村居民人均純收入增長變化是人均水產(chǎn)品消費量增長的單向Granger原因，即農(nóng)村居民收入增長有助于提高水產(chǎn)品消費量；農(nóng)村居民人均純收入對于人均水產(chǎn)品消費量的沖擊力度大于人均水產(chǎn)品消費量對于農(nóng)村居民人均純收入的沖擊力度，促進農(nóng)村居民收入水平增長是提高水產(chǎn)品消費的重要措施。

水產(chǎn)品消費量； 農(nóng)村居民人均純收入； 動態(tài)關(guān)系； 協(xié)整理論； VAR模型

新中國成立后，受當(dāng)時漁業(yè)生產(chǎn)條件與國民經(jīng)濟發(fā)展水平等條件限制，水產(chǎn)品長期處于“奢侈品”狀態(tài)，少有人能夠品嘗。改革開放以來，在國家各項積極政策的扶持下，漁業(yè)產(chǎn)業(yè)規(guī)模不斷發(fā)展壯大，水產(chǎn)品產(chǎn)量連年攀升，“吃魚難”問題較快得到了緩解和解決[1]，人均水產(chǎn)品消費量呈現(xiàn)出明顯的增長趨勢。尤其是進入新世紀(jì)以來，農(nóng)村居民人均水產(chǎn)品消費量增速高于城鎮(zhèn)居民約5.83%[2]，但人均消費量絕對數(shù)值仍低于城鎮(zhèn)居民[3]。在農(nóng)村居民收入實現(xiàn)較快增長的背景下，提高農(nóng)村居民水產(chǎn)品消費量具有潛在的積極拉動力，并對調(diào)整水產(chǎn)品市場供給結(jié)構(gòu)和漁業(yè)產(chǎn)業(yè)轉(zhuǎn)型升級等具有聯(lián)動作用，因此，研究農(nóng)村居民水產(chǎn)品消費特征及影響因素成為重要課題。

水產(chǎn)品是優(yōu)質(zhì)動物蛋白質(zhì)、脂溶性維生素和礦物質(zhì)的良好來源，同時還具有蛋白質(zhì)含量較高、脂肪含量較低、不飽和脂肪酸含量較高等“健康食品”的特征，越來越受到消費者的青睞[4]。關(guān)于農(nóng)村居民水產(chǎn)品消費特征已有大量研究。例如：王恩濤利用固定彈性需求函數(shù)模型對主要食物的收入彈性進行了分析，結(jié)果發(fā)現(xiàn)農(nóng)村居民水產(chǎn)品收入彈性較高[5]；馬冠生等對城鄉(xiāng)居民水產(chǎn)品消費特征的研究顯示，在食用率、食用頻率指標(biāo)方面，城鎮(zhèn)居民均高于農(nóng)村居民[6]；張玉梅等利用QUAIDS 方法模擬水產(chǎn)品消費，認(rèn)為農(nóng)村居民收入水平對水產(chǎn)品消費需求具有顯著制約性[7]；關(guān)潔認(rèn)為經(jīng)濟收入、消費習(xí)慣、水產(chǎn)品的加工程度等是影響我國農(nóng)村居民水產(chǎn)品消費的主要因素，并以山東省調(diào)查樣本為例，實證分析發(fā)現(xiàn)受訪者農(nóng)業(yè)純收入等對其魚類產(chǎn)品的消費有顯著影響[8]。從上述研究來看，農(nóng)村居民收入狀況對于水產(chǎn)品消費具有重要影響作用，但具體的影響程度，不同學(xué)者采用了不同的方法、數(shù)據(jù)，在結(jié)論方面也存在一定差異，尤其是對于農(nóng)村居民收入與水產(chǎn)品消費量之間的動態(tài)關(guān)系鮮見有論述?；诖?，本文借助2000～2015年農(nóng)村居民人均純收入、農(nóng)村居民人均水產(chǎn)品消費量指標(biāo)，利用協(xié)整分析和向量自回歸模型，探討農(nóng)村居民收入水平對水產(chǎn)品消費量的動態(tài)影響作用，以期進一步豐富農(nóng)村居民水產(chǎn)品消費研究。

1 數(shù)據(jù)來源與處理

1.1指標(biāo)選擇與數(shù)據(jù)來源

本文選取2000～2015年作為研究周期，以實證反映進入新世紀(jì)以來農(nóng)村居民水產(chǎn)品消費特征。在具體衡量指標(biāo)選擇方面，本文的“經(jīng)濟收入指標(biāo)”采用農(nóng)村居民人均純收入(IR)，“水產(chǎn)品消費量”采用農(nóng)村居民人均水產(chǎn)品消費量(CR)指標(biāo)，上述指標(biāo)數(shù)據(jù)均來自《中國統(tǒng)計年鑒》[9-10]。

1.2數(shù)據(jù)預(yù)處理

為消除價格因素對農(nóng)村居民收入指標(biāo)的影響，以及可能造成的研究結(jié)論失真問題，本文利用《中國統(tǒng)計年鑒》歷年農(nóng)村居民消費價格指數(shù)，以2000年為基年進行平減，以便消除價格因素的擾動。平減公式為：

IRt=IR2000×(價格指數(shù)t/價格指數(shù)2000) (1)

式中IRt表示t年農(nóng)村居民人均純收入平減后數(shù)值，IR2000表示基年農(nóng)村居民人均純收入，價格指數(shù)t表示t年農(nóng)村居民消費價格指數(shù)，價格指數(shù)2000表示2000年農(nóng)村居民消費價格指數(shù)。

為消除平減后數(shù)據(jù)IR與CR指標(biāo)可能存在的異方差現(xiàn)象，對IR和CR指標(biāo)分別取自然對數(shù)，并記為LNIR和LNCR。數(shù)據(jù)如表1所示。

1.2研究方法

本文采用計量軟件EViews6.0分析農(nóng)村居民收入水平與人均水產(chǎn)品消費量之間的關(guān)系。其中，長期影響關(guān)系分析運用了協(xié)整理論，并構(gòu)建協(xié)整方程進行定量分析；而短期影響關(guān)系則利用向量自回歸模型、脈沖響應(yīng)函數(shù)以及方差分解等進行量化分析。

表1 農(nóng)村居民收入和水產(chǎn)品消費量處理后數(shù)據(jù)Tab.1 Processed data of per capita net income of rural households and per capita consumption of aquatic products

2 結(jié)果分析

2.1平穩(wěn)性檢驗

協(xié)整理論是由Engle和Granger 于1978年首次提出的，主要用于研究兩個或多個非平穩(wěn)變量之間的均衡關(guān)系，并要求進行協(xié)整分析首先要對變量進行平穩(wěn)性檢驗。本文選擇較為常用的ADF檢驗法來進行變量平穩(wěn)性檢驗[11]。利用EViews6.0分別繪制LNCR和LNIR的趨勢圖及其一階差分序列圖，LNCR和LNIR的一階差分別用D(LNCR)和D(LNIR)表示，具體如圖1和圖2所示。其中，從圖1可以大致看出LNCR和LNIR為非平穩(wěn)序列，兩者均表現(xiàn)出一定的增長趨勢；從圖2可以看出，D(LNCR)和D(LNIR)序列則圍繞均衡位置上下波動。

圖1 LNCR和LNIR趨勢圖Fig.1 Trend of LNCR and LNIR

運用EViews6.0對LNCR和LNIR兩個變量序列進行ADF檢驗，結(jié)果顯示兩個變量序列的ADF單位根檢驗值均大于臨界值(5%顯著水平)，因此LNCR和LNIR為非平穩(wěn)序列。進一步對D(LNCR)和D(LNIR)兩個序列進行ADF檢驗，結(jié)果表明一階差分序列的ADF檢驗值小于臨界值(5%顯著水平)，表明D(LNCR)和D(LNIR)是平穩(wěn)的，說明LNCR和LNIR是一階單整，可以進行協(xié)整檢驗。結(jié)果如表2所示。

圖2 LNCR和LNIR一階差分序列趨勢圖Fig.2 First Order Differential Sequence Trends of LNCR and LNIR

表2 序列ADF檢驗結(jié)果Tab.2 Result of ADF test of the variable

注：檢驗形式(C,T,L)中的C,T,L分別表示截距、時間趨勢和滯后期

Note: C,T,L in the test form(C,T,L) denotes intercept, time trend and largest lag length.

2.2協(xié)整模型構(gòu)建

Engle-Granger兩步法和Johansen協(xié)整檢驗方法是協(xié)整檢驗常用的兩種方法，本文采用前者對LNCR和LNIR之間協(xié)整關(guān)系的存在性進行檢驗[11]。

2.2.1 建立變量的最小二乘回歸

利用Eviews6.0對LNCR和LNIR進行OLS回歸，具體模型表達如下：

LNCR=-6.956 1+1.088 2LNIR

(2)

(-6.721 4) (8.309 5)

R2=0.831 4，R2(Adjusted)=0.819 4，Log likelihood=20.282 7，F(xiàn)=69.048 2，括號內(nèi)為T檢驗結(jié)果。

2.2.2 殘差序列平穩(wěn)性檢驗

運用EViews6.0提取式(2)的殘差，記為e，并進行ADF檢驗。由表3可知，ADF檢驗值-2.040 7小于顯著水平5%的臨界值-1.966 3。因此，殘差序列e是平穩(wěn)的。上述2個步驟的結(jié)果表明， LNCR和LNIR之間存在長期均衡關(guān)系，其中式(2)即為兩者之間動態(tài)關(guān)系的數(shù)學(xué)模型。從式(2)的表達式可以看出，農(nóng)村居民人均純收入對水產(chǎn)品消費量具有正向促進作用，并且農(nóng)村居民人均純收入每增加一個單位，人均水產(chǎn)品消費量會相應(yīng)增加1.088 2個單位。

表3 殘差序列e的ADF檢驗結(jié)果Tab.3 ADF test result of series e

2.3格蘭杰(Granger)因果關(guān)系檢驗

協(xié)整檢驗結(jié)論說明農(nóng)村居民人均純收入增長與水產(chǎn)品消費之間存在長期均衡關(guān)系，但要判斷兩者之間是否存在因果關(guān)系，則需要通過格蘭杰因果關(guān)系檢驗進行驗證。通過不同滯后期的選擇，對LNCR和LNIR的格蘭杰因果關(guān)系檢驗，結(jié)果如表4所示。

從表4可以看出，滯后期為1時，在0.098 7的概率下，拒絕假設(shè)“LNIR does not Granger Cause LNCR”，即LNIR是LNCR的格蘭杰原因；在0.454 6的概率下，接受假設(shè)“LNCR does not Granger Cause LNIR”，即LNCR不是LNIR的格蘭杰原因。檢驗結(jié)果表明，短期內(nèi)，農(nóng)村居民人均純收入增長是水產(chǎn)品消費量增長的原因。

表4 格蘭杰因果關(guān)系檢驗Tab.4 Granger causality test

2.4誤差修正模型建立

由于LNCR和LNIR之間具有協(xié)整關(guān)系，因此可以建立誤差修正模型，用以反映農(nóng)村居民人均純收入與水產(chǎn)品消費之間的長期均衡關(guān)系及短期偏離向長期均衡修正的機制。首先通過EViews6.0判斷滯后期，結(jié)果如表5所示。

表5 VEC模型滯后期選擇Tab.5 Lag order selection criteria of VEC

從表5可以看出，LR、FPE、AIC、SC、HQ的最佳滯后期為1，因此將LNCR和LNIR的誤差修正模型最佳滯后期定為1。運用EViews6.0建立誤差修正模型，結(jié)果如下：

DLNCR=0.045 3-0.288 1DLNIR-0.391 3ecmt-1

(3)

(1.8 285) (0.3 844) (-1.6 766)

R2=0.218 3 R2(adjusted)= 0.088 1 D.W.= 2.230 3

式(3)中，括號內(nèi)為T檢驗結(jié)果，ecmt-1表示非均衡誤差，其系數(shù)大小反映了對偏離長期均衡的調(diào)整力度。式中ecmt-1的系數(shù)為-0.391 3，符合誤差反向修正機制，說明當(dāng)農(nóng)村居民人均純收入、水產(chǎn)品消費量短期波動偏離長期均衡關(guān)系后，將以0.391 3的調(diào)整力度將非均衡狀態(tài)拉回到均衡狀態(tài)。

2.5脈沖響應(yīng)函數(shù)

為分析農(nóng)村居民人均純收入增長與水產(chǎn)品消費之間的相互影響程度，利用EViews6.0建立VAR(1)模型，進行脈沖響應(yīng)函數(shù)分析，結(jié)果如圖3所示。圖3橫軸表示隨機擾動項沖擊作用的滯后期數(shù)(本文取10期)，縱軸為響應(yīng)指數(shù)，數(shù)值越大表示對擾動的響應(yīng)程度越高。其中LNIR對LNCR的響應(yīng)過程為，在第一期出現(xiàn)負(fù)響應(yīng)，之后呈現(xiàn)一定的正向響應(yīng)，幅度不大，說明農(nóng)村居民人均純收入的短期增長對水產(chǎn)品消費會產(chǎn)生影響，長期持續(xù)增長的話，對于水產(chǎn)品消費量增長有促進作用，但強度不大。LNCR對LNIR的響應(yīng)過程為，在第一期的響應(yīng)為零，之后逐漸呈現(xiàn)正向響應(yīng)過程，說明水產(chǎn)品消費量對農(nóng)村居民人均純收入具有正向促進作用，且與LNIR對LNCR的響應(yīng)相比，響應(yīng)力度較大。

圖3 LNCR和LNIR對一個標(biāo)準(zhǔn)誤擾動的響應(yīng)結(jié)果Fig.3 Response of LNCR and LNIR to a S.E.

2.6方差分解

方差分解可以分析每一個結(jié)構(gòu)沖擊對內(nèi)生變量影響的貢獻度，從而評價不同結(jié)構(gòu)沖擊的重要性，結(jié)果如圖4所示。在水產(chǎn)品消費的方差分解中，LNCR受自身沖擊影響較大，但影響力度不斷減弱；相反，LNIR的變化對其產(chǎn)生的沖擊則不斷加強，到第10期時分別占60%和40%。在農(nóng)村居民人均純收入增長的方差分解中，LNIR受自身沖擊的影響較大，且保持平穩(wěn)，到第10期保持在近93%，而LNCR的沖擊影響則較小，僅為6%。由此可見，農(nóng)村居民人均純收入增長對于促進水產(chǎn)品消費的貢獻非常顯著，而水產(chǎn)品消費對于農(nóng)村居民人均純收入增長的貢獻則并不顯著。

圖4 LNCR和LNIR的方差分解結(jié)果Fig. 4 Result of variance decomposition of LNCR and LNIR

3 結(jié)論與討論

本文利用計量經(jīng)濟方法對中國2000～2015年農(nóng)村居民人均純收入與人均水產(chǎn)品消費量關(guān)系進行了定量研究，結(jié)論如下：

(1)農(nóng)村居民人均純收入與人均水產(chǎn)品消費量之間具有協(xié)整關(guān)系，兩者之間呈現(xiàn)穩(wěn)定的長期正相關(guān)關(guān)系，其中農(nóng)村居民人均純收入每增加一個單位，人均水產(chǎn)品消費量會相應(yīng)增加1.088 2個單位。

(2)利用誤差修正模型分析隨機擾動干擾時，農(nóng)村居民人均純收入與人均水產(chǎn)品消費量的長期關(guān)系對短期關(guān)系的波動具有一定影響，而當(dāng)農(nóng)村居民人均純收入、人均水產(chǎn)品消費量短期波動偏離長期均衡關(guān)系后，將以0.391 3的調(diào)整力度將非均衡狀態(tài)拉回到均衡狀態(tài)。

(3)脈沖響應(yīng)函數(shù)與方差分解的結(jié)果顯示，農(nóng)村居民人均純收入與人均水產(chǎn)品消費量之間的相互作用不同，其中農(nóng)村居民人均純收入對于人均水產(chǎn)品消費量的沖擊力度大于后者對前者的沖擊力度。

(4)農(nóng)村居民人均純收入增長對于人均水產(chǎn)品消費量具有較為明顯的促進作用，努力確保農(nóng)村居民人均純收入持續(xù)穩(wěn)定增長，有助于實現(xiàn)農(nóng)村居民人均水產(chǎn)品消費量保持增長趨勢。

長期以來，中國居民水產(chǎn)品消費量雖然有了較大幅度提高，但與水產(chǎn)品消費大國的平均水平仍有較大差距[12-13]，同時城鄉(xiāng)差距、區(qū)域消費差距特征明顯[14]。從個體消費角度看，收入水平是重要的影響因素，但現(xiàn)有關(guān)于農(nóng)村居民人均純收入與水產(chǎn)品消費之間關(guān)系的定量研究較少，主要集中在收入彈性方面[7,15]。本文利用協(xié)整理論和向量自回歸模型對農(nóng)村居民人均純收入和水產(chǎn)品消費量之間的動態(tài)關(guān)系進行了探討，結(jié)果表明：農(nóng)村居民人均純收入能夠促進水產(chǎn)品消費，但水產(chǎn)品消費對促進農(nóng)村居民人均純收入增長的影響則較為微弱。這一結(jié)論符合現(xiàn)有關(guān)于農(nóng)村居民水產(chǎn)品消費行為特征研究結(jié)論，即收入是影響水產(chǎn)品消費的重要因素[5,7-8]；在收入彈性方面，食物消費升級模式與糧食安全政策分析評估課題組的結(jié)論為0.900 2[15]，而本文結(jié)論為1.088 2，造成差異的原因可能是時間序列數(shù)據(jù)的調(diào)整以及消費特征的改變等。從人均消費量角度看，目前農(nóng)村居民人均水產(chǎn)品消費量僅相當(dāng)于城鎮(zhèn)居民的一半[16]，造成這一差距的原因包括收入因素、文化因素以及水產(chǎn)品加工業(yè)的發(fā)展程度，隨著農(nóng)村居民收入水平的逐步提高、電子商務(wù)平臺在農(nóng)村的普及以及對水產(chǎn)品營養(yǎng)特征的進一步宣傳等因素的影響，廣大農(nóng)村有可能成為水產(chǎn)品消費增長的重要增長動力源。通過倡導(dǎo)健康消費和營養(yǎng)價值理念，創(chuàng)新水產(chǎn)品營銷策略，有助于進一步開發(fā)農(nóng)村水產(chǎn)品消費市場，進而促進漁業(yè)生產(chǎn)、加工及流通環(huán)節(jié)不斷升級，實現(xiàn)漁業(yè)產(chǎn)業(yè)鏈延長和價值鏈提升。

[1] 岳冬冬,王莎,王魯民,等.我國水產(chǎn)品產(chǎn)量波動特征的多維度分析[J]. 南方農(nóng)業(yè)學(xué)報, 2013，44(2):344-349.

[2] 孫煒琳,劉佩.近10年來我國水產(chǎn)品供求變化及趨勢分析[J]. 中國農(nóng)業(yè)信息, 2014(3):61-64.

[3] 董曉霞.我國生鮮農(nóng)產(chǎn)品居民消費差異分析[J].中國食物與營養(yǎng),2005(8):25-27.

[4] WANG S, LAY S, YU H,etal. Dietary guidelines for Chinese residents(2016):comments and comparisons[J]. Journal of Zhejiang University Science, 2016,17(9):649-656.

[5] 王恩濤.中國農(nóng)村居民食物消費與營養(yǎng)狀況研究[D].北京：中國農(nóng)業(yè)科學(xué)院, 2007.

[6] 馬冠生,胡小琪,崔朝輝,等.我國居民食物消費現(xiàn)狀與特點[J].中國食物與營養(yǎng),2007(7): 4-6.

[7] 張玉梅,喻聞,李志強.中國農(nóng)村居民食物消費需求彈性研究[J].江西農(nóng)業(yè)大學(xué)學(xué)報(社會科學(xué)版),2012,11(2):7-13.

[8] 關(guān)潔.我國農(nóng)村居民水產(chǎn)品消費影響因素探析——以山東農(nóng)村居民魚類產(chǎn)品消費為例[J].生產(chǎn)力研究,2013(7): 25-26.

[9] 中華人民共和國國家統(tǒng)計局.2016中國統(tǒng)計年鑒[M].北京:中國統(tǒng)計出版社, 2016.

[10] 中華人民共和國國家統(tǒng)計局.1999中國統(tǒng)計年鑒[M].北京:中國統(tǒng)計出版社, 1999.

[11] 高鐵梅.計量經(jīng)濟分析方法與建模:EViews應(yīng)用及實例(第2版)[M].北京:清華大學(xué)出版社,2009.

[12] 王國華.日本水產(chǎn)品消費的變動與啟示[J]. 世界農(nóng)業(yè), 2012(2):66-69.

[13] 張琳. 世界人均水產(chǎn)品消費走勢分析[J]. 漁業(yè)致富指南,2006(16):4-4.

[14] 小遠.中國水產(chǎn)品消費動向統(tǒng)計分析[J].漁業(yè)致富指南,2010(12):17-19.

[15] 食物消費升級模式與糧食安全政策分析評估課題組.中國農(nóng)民食物消費研究[M].北京:中國農(nóng)業(yè)出版社,2007.

[16] 岳冬冬,耿瑞,趙文武.西北地區(qū)中青年群體水產(chǎn)品消費特征研究[J].農(nóng)產(chǎn)品質(zhì)量與安全,2017(3):85-90.

OntherelationshipbetweenpercapitanetincomeofruralresidentsandpercapitaconsumptionofaquaticproductsinChinabasedontheVARmodel

YUE Dong-dong1, WANG Lu-min2, JI Wei-wei1, RUAN Wen2, WANG Qian1, XIONG Min-si2, XIAO Li1, ZHENG Liang2

(1.KeyLaboratoryofOceanicandPolarFisheries,MinistryofAgriculture;EastChinaSeaFisheriesResearchInstitute,ChineseAcademyofFisherySciences,Shanghai200090,China; 2.KeyLaboratoryofEastChinaSeaFisheryResourcesExploitation,MinistryofAgriculture;EastChinaSeaFisheriesResearchInstitute,ChineseAcademyofFisherySciences,Shanghai200090,China)

This paper analyzed the relationship between the per capita net income of China’s rural residents and their per capita consumption of aquatic products from 2000 to 2015 by using the method of econometric analysis such as cointegration theory and vector autoregression model. The results show that there is a long-term cointegration relationship between per capita net income of rural residents and per capita consumption of aquatic products. The change in per capita net income of rural households is the one-way Granger cause of per capita consumption of aquatic products, that is to say, the increase in per capita net income of rural households will help to increase the consumption of aquatic products. The impact of per capita net income of rural residents on per capita consumption of aquatic products is greater than the impact of per capita consumption of aquatic products on per capita net income of rural residents, and promoting per capita net income of rural households growth will be an important measure to increase consumption of aquatic products.

consumption of aquatic products; per capita net income of rural households; dynamic relationship; cointegration theory; VAR model

2095-3666(2017)03-0161-07

10.13233/j.cnki.fishis.2017.03.001

2017-5-20

：2017-07-19

農(nóng)業(yè)部委托項目“水產(chǎn)品消費量與人均GDP關(guān)系研究”

岳冬冬(1982-)，博士，副研究員，研究方向為漁業(yè)經(jīng)濟與發(fā)展戰(zhàn)略。E-mail:yuedong1982@126.com

F 326.4

：A