• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    中國現(xiàn)階段環(huán)境約束下的適度經(jīng)濟增長率研究——基于廢氣排放量視角的分析

    2014-05-12 10:22:20李寶瑜
    統(tǒng)計與信息論壇 2014年4期
    關(guān)鍵詞:城鎮(zhèn)化率約束條件增長率

    李寶瑜,劉 洋

    (山西財經(jīng)大學(xué) 統(tǒng)計學(xué)院,山西 太原 030006)

    中國現(xiàn)階段環(huán)境約束下的適度經(jīng)濟增長率研究
    ——基于廢氣排放量視角的分析

    李寶瑜,劉 洋

    (山西財經(jīng)大學(xué) 統(tǒng)計學(xué)院,山西 太原 030006)

    經(jīng)濟發(fā)展與環(huán)境污染的關(guān)系歷來都是一對不可回避的矛盾,也是學(xué)術(shù)界研究的熱點,國內(nèi)外大多數(shù)學(xué)者都在考察人均GDP與環(huán)境污染的庫茲涅茨曲線,但對于環(huán)境與經(jīng)濟增長率的關(guān)系卻很少有人研究。首先建立環(huán)境變量與經(jīng)濟增長率關(guān)系的變參數(shù)狀態(tài)空間模型,并將建立的狀態(tài)空間模型轉(zhuǎn)化成線性規(guī)劃中的決策方程,之后應(yīng)用參數(shù)區(qū)間估計和集合論兩種方法,通過約束相應(yīng)的變量及變量的系數(shù)來估計中國現(xiàn)階段環(huán)境約束下的適度經(jīng)濟增長率。

    適度經(jīng)濟增長;環(huán)境約束;狀態(tài)空間模型;線性規(guī)劃模型

    一、引 言

    改革開放以來,中國的經(jīng)濟增速舉世矚目,同時,高增速所帶來的環(huán)境污染也越來越嚴重。因此,我們必須思考環(huán)境和經(jīng)濟之間到底是怎樣的關(guān)系,經(jīng)濟加速發(fā)展是否必然帶來環(huán)境污染,先污染后治理的規(guī)律能否打破。

    關(guān)于經(jīng)濟發(fā)展和環(huán)境污染的關(guān)系,主要圍繞環(huán)境庫茲涅茨曲線展開。該理論認為,環(huán)境污染與經(jīng)濟發(fā)展呈倒“U”型曲線(環(huán)境隨著經(jīng)濟的發(fā)展會經(jīng)歷一個先惡化后改善的過程),通常采用的變量是人均國民收入或相應(yīng)經(jīng)濟總量。然而,還有另一條研究路徑常常被人忽視,這就是經(jīng)濟增長率與環(huán)境污染的關(guān)系研究。從不同國家的數(shù)據(jù)可以看出:增長率與不同的環(huán)境變量呈現(xiàn)的是一種具有不同斜率的線性關(guān)系,絕不是隨著經(jīng)濟增長率的提高,環(huán)境污染呈現(xiàn)先走高、后平穩(wěn)、然后再下降的倒“U”型曲線關(guān)系。既然環(huán)境污染與增長率之間呈線性關(guān)系,就必然出現(xiàn)要加速經(jīng)濟發(fā)展就必須付出環(huán)境代價的問題,經(jīng)濟與環(huán)境就會陷入一種矛盾之中。那么,是否存在一個能夠使二者互相妥協(xié),共生共進的適度經(jīng)濟增長率呢?能否找到一個污染最小條件下的適度經(jīng)濟增長率區(qū)間?本文將通過建立模型,尋找環(huán)境約束下的適度經(jīng)濟增長率。

    二、文獻綜述

    在環(huán)境與經(jīng)濟關(guān)系這個大命題下,有著眾多不同的研究路徑。1991年美國經(jīng)濟學(xué)家Grossman和Krueger提出經(jīng)濟發(fā)展與環(huán)境發(fā)展的關(guān)系呈倒“U”型曲線關(guān)系的理論,圍繞這一現(xiàn)象是否具有普遍意義,經(jīng)濟增長是否可以避免環(huán)境污染的問題就出現(xiàn)了很大的爭論。一些學(xué)者認為環(huán)境庫茲涅茨曲線(EKC)不可避免,而另一些學(xué)者則認為可以打破。支持環(huán)境庫茲涅茨曲線存在的學(xué)者Markus等,從經(jīng)濟規(guī)模、結(jié)構(gòu)、技術(shù)、收入差異、國際貿(mào)易與環(huán)境關(guān)系的理論角度證明環(huán)境倒“U”型曲線的存在[1]。另一些不支持環(huán)境庫茲涅茨曲線的學(xué)者如Dasgupta等,通過實證研究發(fā)現(xiàn)某些污染物與人均收入水平?jīng)]有倒“U”型關(guān)系,而是線性甚至“N”型關(guān)系[2]。

    在中國,這些年來經(jīng)濟增速創(chuàng)造了世界奇跡,但是否進入了倒“U”型軌道?先污染后治理的規(guī)律是否能夠被打破?不同的人對此也有不同的看法。中國一些學(xué)者注意到EKC曲線是經(jīng)驗數(shù)據(jù)的描述,并不能說明人均收入和環(huán)境之間存在必然聯(lián)系。不同的省份和環(huán)境變量會呈現(xiàn)不同的曲線。他們通過理論與實證研究認為EKC曲線并不存在[3]。

    環(huán)境與經(jīng)濟增長的另一路徑是將環(huán)境變量引入內(nèi)生經(jīng)濟增長模型所做的研究。自從20世紀(jì)80年代中后期以來,以Romer、Lucas等人為代表的內(nèi)生增長模型的出現(xiàn),經(jīng)濟學(xué)家們開始將環(huán)境變量引入生產(chǎn)函數(shù),把環(huán)境質(zhì)量引入效用函數(shù),進而在內(nèi)生增長模型框架下討論環(huán)境惡化與經(jīng)濟可持續(xù)發(fā)展問題。代表模型有:Bovenberg等在Romer模型的基礎(chǔ)上將環(huán)境變量引入生產(chǎn)函數(shù)的模型,Ligthart通過將“AK”模型進行擴展來研究環(huán)境污染外部性與經(jīng)濟持續(xù)增長問題,Grimaud等將環(huán)境污染和有限的不可再生資源引入新熊彼特模型的研究。上述敘述均參考Panayotou的綜述[4]。

    而國內(nèi)對于該問題的研究較少,大多是在國外研究成果的基礎(chǔ)上進行模型改進。梅國平等將環(huán)境質(zhì)量內(nèi)生化,構(gòu)建Hamilton動態(tài)優(yōu)化模型尋求均衡狀態(tài)[5];岳書敬等使用三種不同的模型,考察了在環(huán)境因素的約束下,工業(yè)的增長效率及其影響因素[6];也有使用面板數(shù)據(jù)進行省際間的實證研究。

    關(guān)于環(huán)境約束下的經(jīng)濟增長問題,國內(nèi)外的研究主要有兩種研究思路:是研究環(huán)境庫茲涅茨曲線是否可以避免;二是在內(nèi)生經(jīng)濟增長模型的基礎(chǔ)上,使用動態(tài)優(yōu)化方法進行求解,得出在相應(yīng)條件下經(jīng)濟能夠持續(xù)增長的結(jié)論。相應(yīng)的研究都是在求解過程中使用經(jīng)濟總量為研究對象,得出了固定的經(jīng)濟增長率,并沒有直接對經(jīng)濟增長率進行研究。本文擬直接對經(jīng)濟增長率與環(huán)境變量的關(guān)系進行研究。

    為了求證二者關(guān)系,首先需要找到一個環(huán)境典型變量和一個經(jīng)濟增長變量。經(jīng)濟增長變量可以用目前通行的GDP增長率來度量,環(huán)境變量則比較復(fù)雜。環(huán)境污染有多個方面,每個方面都是以物質(zhì)計量單位統(tǒng)計的,這些方面目前還難以貨幣化計量,也缺乏一個綜合指標(biāo),所以我們選擇用某個典型變量作為環(huán)境代表變量,本文選用廢氣排放量進行實際建模和分析。而對于環(huán)境的其他變量,諸如廢水排放量或者其它某個環(huán)境指標(biāo)與經(jīng)濟增長率的關(guān)系,也可參照本文的方法進行建模和分析,原理類似。

    本文設(shè)計了一套廢氣排放約束下的經(jīng)濟增長模型,完整的方案包括兩個模型和兩種方法。具體思路為:建立狀態(tài)空間計量模型估計環(huán)境(廢氣排放)與經(jīng)濟增長的關(guān)系;將估計得到的環(huán)境模型作為目標(biāo)函數(shù)建立線性規(guī)劃模型,給出約束條件;在各種細分約束條件下,求解多維約束和單項約束下的適度經(jīng)濟增長率;使用區(qū)間估計和集合求交集法,分別求解多維約束下和多個單項約束下的適度經(jīng)濟增長率區(qū)間。

    三、用狀態(tài)空間模型估計環(huán)境與經(jīng)濟增長率的關(guān)系

    (一)變量與模型選擇

    經(jīng)驗證明,如果以年度作為研究對象,經(jīng)濟增長率對環(huán)境有顯著的影響,然而環(huán)境對經(jīng)濟增長率的影響盡管存在,但它卻是一個緩慢和復(fù)雜的過程,不一定表現(xiàn)在年度對應(yīng)關(guān)系上,不過環(huán)境變化自身卻可能在時間上有連續(xù)性,即本期的環(huán)境受前期影響較大。基于這種考慮,我們選擇能夠反應(yīng)環(huán)境連續(xù)變化的以環(huán)境變化作為被解釋變量的可變參數(shù)狀態(tài)空間模型。

    環(huán)境變化特別是空氣質(zhì)量變化不僅受經(jīng)濟增長率的影響,同時還受其它因素的影響,如城鎮(zhèn)化速度、能源消耗水平、經(jīng)濟增長對環(huán)境的彈性系數(shù)等,以環(huán)境變量作為被解釋變量時,還必須選擇相關(guān)解釋變量。

    我們以工業(yè)廢氣排放總量對數(shù)(lnFQ)作為環(huán)境變量,經(jīng)濟增長率(ZZL)、城鎮(zhèn)化率(CZHL)、能源強度(NYQD)、lnFQ的滯后一期作為解釋變量,分別從增長速度、城鎮(zhèn)化水平和生產(chǎn)技術(shù)水平三個方面來解釋經(jīng)濟發(fā)展造成的污染和影響,建立如下狀態(tài)空間模型:

    其中式(1)是量測方程,反映了各要素與廢氣排放總量的關(guān)系,α是截距項,β、γ、η分別表示廢氣排放總量對經(jīng)濟增長率、能源強度和城鎮(zhèn)化率變動的敏感程度,εt是均值為0、協(xié)方差連續(xù)不相關(guān)的擾動項;φ為可變參數(shù),反映各個時點上廢氣排放總量受前一期影響的程度。式(2)為狀態(tài)方程,描述了可變參數(shù)φ的生成過程。

    (二)數(shù)據(jù)的預(yù)處理

    本文選擇的數(shù)據(jù)為年度數(shù)據(jù),相應(yīng)的原始數(shù)據(jù)均可從歷年《中國統(tǒng)計年鑒》獲得,在建立模型之前需要對所選擇的指標(biāo)數(shù)據(jù)進行一定的預(yù)處理。

    工業(yè)廢氣排放總量對數(shù)(lnFQ):中國從1985年開始公布環(huán)境污染的相關(guān)數(shù)據(jù),我們可以從《中國統(tǒng)計年鑒》和《環(huán)境統(tǒng)計年鑒》獲得工業(yè)廢氣排放總量(億標(biāo)立方米)1985-2010年的數(shù)據(jù),但是從2011年開始中國不再公布工業(yè)廢氣排放總量,改為公布廢氣中主要污染物的排放情況,包括二氧化硫、氮氧化物和煙(粉)塵,單位為萬噸。數(shù)據(jù)公布口徑的變化使得工業(yè)廢氣排放總量這一指標(biāo)2011年的數(shù)據(jù)缺失,本文以2011年中國工業(yè)的增長率代替工業(yè)廢氣排放總量的增長率,計算出2011年的工業(yè)廢氣排放總量來彌補數(shù)據(jù)的缺失。對廢氣排放總量序列取對數(shù),得到lnFQ序列。

    經(jīng)濟增長率(ZZL):本文選用的增長率指標(biāo)為全國GDP剔除價格因素的環(huán)比經(jīng)濟增長率。

    城鎮(zhèn)化率(CZHL):按照慣例,城鎮(zhèn)化率這一指標(biāo)用城鎮(zhèn)人口占總?cè)丝诘谋戎貋肀硎?。中國城?zhèn)化率的變化比較穩(wěn)定,但在1995年存在一個拐點,1995年之前基本以每年0.5%的速度增長,在此之后,每年上升約1.4%。這種數(shù)據(jù)上的異??赡芨y(tǒng)計口徑的變化有關(guān),并不符合實際情況。因此,本文以1995年的數(shù)據(jù)為基礎(chǔ),用1995年至2011年中國城鎮(zhèn)化率的平均變化速度代替之前的變化速度,得到了消除拐點的城鎮(zhèn)化率序列。

    能源強度(NYQD):能源強度即單位GDP能耗,單位是噸標(biāo)準(zhǔn)煤/萬元。中國只公布了各省份的萬元GDP耗能,并沒有全國的能源強度水平,因此需要用全國的GDP和能源消耗總量來計算。為了保證各年能源強度的可比性,首先對GDP進行價格平減,將各年的GDP換算為2011年的價格水平,再用能源消耗總量除以GDP計算能源強度的值。

    (三)模型的估計

    1.變量的平穩(wěn)性及協(xié)整關(guān)系檢驗

    狀態(tài)空間方法構(gòu)造的變參數(shù)模型要求方程中的變量是平穩(wěn)的,或變量之間存在協(xié)整關(guān)系,因此,需要對變量進行單位根檢驗,檢驗結(jié)果見表1。

    表1 變量的單位根檢驗結(jié)果表

    由上表可知,lnFQ、ZZL、CZHL和NYQD四個變量均為一階單整,符合建立協(xié)整關(guān)系的前提條件。另外,還需對變量間的協(xié)整關(guān)系進行Johansen檢驗,結(jié)果見表2。

    表2 變量的協(xié)整檢驗結(jié)果表

    協(xié)整檢驗結(jié)果表明變量間存在協(xié)整關(guān)系,這樣就避免偽回歸,可以建立狀態(tài)空間模型,進一步分析變量間的動態(tài)關(guān)系。

    2.模型的估計結(jié)果與檢驗

    根據(jù)式(1)和式(2)式中模型的設(shè)定,運用Eviews6.0對模型進行求解,得到狀態(tài)空間模型的估計結(jié)果。

    經(jīng)檢驗,該模型R2=0.997,殘差的ADF單位根檢驗為平穩(wěn),不存在一階自相關(guān)。

    四、用線性規(guī)劃模型求解環(huán)境約束下的經(jīng)濟增長率集合

    (一)目標(biāo)函數(shù)和約束條件

    1.目標(biāo)函數(shù)的建立

    前面通過狀態(tài)空間模型估計了廢氣與經(jīng)濟增長率的關(guān)系,根據(jù)計量模型的假定,在殘差平穩(wěn)的條件下,殘差的均值是0,所以可以省略誤差項,將估計的關(guān)系式作為確定性關(guān)系式應(yīng)用,將其作為線性規(guī)劃模型的目標(biāo)函數(shù)。因為要求解的是最優(yōu)經(jīng)濟增長率,所以將原狀態(tài)空間方程經(jīng)濟增長率變量放到目標(biāo)函數(shù)的左邊,廢氣變量放到方程的右邊,這樣,就把原來的環(huán)境污染最小化問題,轉(zhuǎn)化成經(jīng)濟增長率最大化問題。線性規(guī)劃模型為:

    2.約束條件的確定

    根據(jù)狀態(tài)空間模型的變系數(shù)方程,可以約束方程中的系數(shù),也可以約束目標(biāo)方程中的變量。設(shè)計約束條件的指導(dǎo)思想是:國家有規(guī)劃目標(biāo)的,按國家規(guī)劃目標(biāo)測算;國家沒有規(guī)劃目標(biāo)的,按理論標(biāo)準(zhǔn)或?qū)嶋H數(shù)據(jù)趨勢劃定約束界限,其中有的要設(shè)定上限值,有的則約束下限值。在本文的目標(biāo)規(guī)劃方程中,與經(jīng)濟增長率有關(guān)的變量除了廢氣排放量,還有能源強度、城鎮(zhèn)化率以及經(jīng)濟增長率對廢氣排放量的彈性。下文對于廢氣排放量的約束是借鑒“十二五”規(guī)劃目標(biāo)值,其他變量也需要用與廢氣發(fā)展?fàn)顟B(tài)相對應(yīng)的約束值才能求得環(huán)境約束下的最大經(jīng)濟增長率。因此,在該線性規(guī)劃中,除了要對廢氣排放量進行約束,還需要對其他的有關(guān)變量進行約束。

    廢氣排放約束:根據(jù)國家出臺的《十二五規(guī)劃》等一系列相關(guān)經(jīng)濟、環(huán)境政策文件,給出“十二五”期間具體的目標(biāo)值,見表3。

    表3 “十二五”規(guī)劃目標(biāo)值表

    《綱要》上并沒有給出廢氣排放量約束的具體指標(biāo),而是給出了廢氣中二氧化硫、氮氧化物的約束指標(biāo),本文將全國二氧化硫的約束指標(biāo)10%作為廢氣排放量下降的上限,將全國二氧化氮下降7%的約束指標(biāo)作為廢氣排放量下降的下限,則現(xiàn)階段廢氣排放量降幅的約束區(qū)間為[7%,10%],對應(yīng)具體廢氣排放量的區(qū)間為[467 251.2,482 849.1]。

    城鎮(zhèn)化率約束:《綱要》中要求到2015年達到51.5,考慮到2011年的中國城鎮(zhèn)化率已經(jīng)為51.3,通過由中國公布的數(shù)據(jù)計算得到近10年中國城鎮(zhèn)化率的增長速度在2%~3%的范圍,按照這個速度,到“十二五”末期,中國城鎮(zhèn)化率的區(qū)間約束為[56,58]。

    單位國內(nèi)生產(chǎn)總值能源消耗約束:《綱要》中要求在2015年單位GDP能耗降低16%,按照2010年的GDP能耗下降16%得到0.63噸/萬元,因此本文對GDP能耗進行約束時,要求其小于0.63噸/萬元。考慮GDP能耗現(xiàn)階段不可能為0,因此其下限一定大于0。

    廢氣排放與經(jīng)濟增長率的彈性系數(shù)約束:彈性系數(shù)的求解,根據(jù)方程關(guān)系,可以推導(dǎo)出彈性系數(shù)公式為:

    根據(jù)公式,可以算出經(jīng)濟增長率對工業(yè)廢氣排放總量增長率的彈性系數(shù)。計算結(jié)果表明:中國經(jīng)濟增長率與環(huán)境污染呈斜率不同的線性變動關(guān)系,并且該斜率的變化有一定的波動范圍,取值區(qū)間為[0.18,0.27],平均值為0.22。在2000年以前,彈性系數(shù)主要在0.20~0.22區(qū)間變動,2001年開始,彈性系數(shù)隨時間呈線性增長趨勢(除2008年經(jīng)濟危機有所降低外)??紤]到現(xiàn)階段國家對于環(huán)保投資加大,淘汰高能耗產(chǎn)業(yè)以及科技進步等因素,該系數(shù)的上限現(xiàn)階段不會超過0.18,同時考慮到其下降的速度不會太快,故本文設(shè)定其現(xiàn)階段的下限為0.16。由此,確定了約束條件,對應(yīng)到對數(shù)模型中經(jīng)濟增長率前的系數(shù)約束區(qū)間為[0.011 324,0.012 484]。

    廢氣滯后系數(shù)平均化處理。由于方程中含有廢氣的滯后一期項,因此其系數(shù)原則上不能超過1,對應(yīng)到對數(shù)模型中,其系數(shù)就不能大于0,雖然本文不對其做約束,但為了簡化計算,對其做了平均化處理??紤]到這個系數(shù)變化并不是很大,取其之前所有系數(shù)的平均值(-0.004 54)作為目標(biāo)方程廢氣排放量滯后一期項前的系數(shù)。

    (二)適度經(jīng)濟增長率求解

    1.求解多維約束條件下的最大經(jīng)濟增長率

    依據(jù)不同的問題,線性規(guī)劃的算法有所不同。通常情況下,線性規(guī)劃的算法有圖解法、單純形算法、基線法(由于它們都是在可行集的邊界上取得最優(yōu)值,故合稱為邊界點算法)、內(nèi)點法、橢球法及其各種改進算法。單純形算法的基本原理是:先找出一個基可行解,判斷其是否為最優(yōu)解,如為否,則轉(zhuǎn)換到相鄰的基可行解,并使目標(biāo)函數(shù)值不斷增大,一直找到最優(yōu)解為止。內(nèi)點算法是針對單形法的“邊界趨近”觀念而改采“內(nèi)部逼近”的路線,相對于只沿著可行域的邊沿進行移動的單純形算法,內(nèi)點算法能夠在可行域內(nèi)移動,本文使用內(nèi)點法進行求解。具體作法為:同時將各種約束的各種可能取值帶入方程,得到各種組合條件下的經(jīng)濟增長率。

    考慮到每一個約束條件都屬于時點數(shù)據(jù),并不連續(xù),因此在求解時,只能選取每個條件中的一個數(shù)值,根據(jù)每個條件的區(qū)間范圍,從每個條件中均選取10個點,在選取每個點時采取下列方法:在增長率前的系數(shù)β、城鎮(zhèn)化率和廢氣排放量的選擇上,都將其約束區(qū)間10等分,取其區(qū)間的上限點或者下限點。在能源強度的選擇上,以0.63為起點,按照1%的下降速度分別求出9個點。因為廢氣排放量的滯后一期的數(shù)值與廢氣排放量是一一對應(yīng)的,根據(jù)之前的廢氣排放量的數(shù)據(jù),就可以求出其滯后一期的數(shù)值。最終選取符合各個條件的數(shù)值見表4。

    表4 各變量細分值表

    上述四個變量都有10個細分點,每個細分點都帶入目標(biāo)方程中計算一次,最后共得到=10 000個經(jīng)濟增長率的樣本點。在這10 000個經(jīng)濟增長率樣本點中最小值為89.72,最大值為120.03。

    2.求解單項約束條件下的最大經(jīng)濟增長率

    另一種方法是對四個約束條件分別進行約束,在其他條件固定的情況下考察單項約束下的經(jīng)濟增長率。最終分別求得四個約束條件下可能的經(jīng)濟增長率的區(qū)間。設(shè)Xi為某一約束變量,a與b分別為其約束的下限和上限,Y1、Y2、Y3分別為其余三個約束條件,k、h、f分別為三個約束取值范圍中的一個固定點,則模型可以表達為:

    單項約束的具體做法是在研究某一約束條件時將其他條件固定不變,可以得出該約束條件的每一個取值所對應(yīng)的經(jīng)濟增長率。例如,廢氣排放單項約束時,僅讓廢氣約束變化,能源強度、城鎮(zhèn)化率和彈性系數(shù)都固定在某一給定的點。對于其它約束的固定取值,確定為各自約束條件下接近于現(xiàn)狀的上限或下限值。四個約束條件的固定取值見表5。

    表5 四個約束條件的固定取值表

    對此模型求解,可以得出每一個單項約束下可能的經(jīng)濟增長率集合與最大經(jīng)濟增長率,具體數(shù)值見表6。

    表6 單項約束下可能的經(jīng)濟增長率表

    從表6可以看到,廢氣排放、能源強度、城鎮(zhèn)化率、彈性系數(shù)各自約束下的最大經(jīng)濟增長率分別為105.13、107.73、107.50、112.94。

    五、適度經(jīng)濟增長率的估計與選擇

    在現(xiàn)階段,與中國的實際情況相結(jié)合,上文所求得的最大經(jīng)濟增長率要么太高,要么太低,并不是適合中國國情的環(huán)境約束下的適度經(jīng)濟增長率。在上文所求得的所有經(jīng)濟增長率集合中,我們使用兩種方法來進行適度經(jīng)濟增長率區(qū)間的選取。

    (一)多維約束下的適度經(jīng)濟增長率求解

    多維約束下的適度經(jīng)濟增長率,是指在多維約束條件下求得的所有經(jīng)濟增長率的取值點中,樣本密集分布區(qū)間的經(jīng)濟增長率。

    上文在多維約束條件下,使用線性規(guī)劃的方法求出了10 000個經(jīng)濟增長率的可能取值點并獲取了最大值。但最大值不一定等于適度值,因為最大值和最小值,都是在各個約束變量取極端值時產(chǎn)生的,因此需要找到一種具有普遍適用性的規(guī)律。根據(jù)統(tǒng)計學(xué)大數(shù)定律,可以將數(shù)據(jù)密集區(qū)域,即圍繞在均值周圍的區(qū)域看做是符合普遍性的適度區(qū)域。圖1展示了前面多維約束條件下產(chǎn)生的10 000個結(jié)果的統(tǒng)計分布。

    圖1 多項約束下可能的經(jīng)濟增長率的概率分布

    從分布圖上直觀可以看出,多維約束下得到的可能的經(jīng)濟增長率大致屬于正態(tài)分布,所以可以在給出一定置信度的條件下,用正態(tài)分布的區(qū)間估計方法求出樣本均值及其浮動范圍作為適度增長率。為了判定其分布,需要計算其統(tǒng)計量,具體結(jié)果見表7。

    表7 多項約束下經(jīng)濟增長速度分布的主要統(tǒng)計量表

    從表7可以看出,樣本的偏度為0.087接近于0,峰度為2.554接近3,考慮其樣本量已經(jīng)達到10 000,屬于大樣本,因此可以認為該總體已經(jīng)近似服從正態(tài)分布。對于總體服從正態(tài)分布,均值和方差都未知的大樣本,使用下述公式進行均值區(qū)間估計:

    (二)單項約束下的最優(yōu)經(jīng)濟增長率求解

    現(xiàn)在換一種思路來求解適度經(jīng)濟增長率區(qū)間。在前面單項約束條件下,各自產(chǎn)生了10個經(jīng)濟增長率點,如果把這些點都匯集到一起,是否會有一些公共的因素呢?如果有公共交集部分,就可以將其看做同時符合各單項約束的適度經(jīng)濟增長率。這樣就把問題轉(zhuǎn)換成了將四個變量所對應(yīng)的增長率所組成的集合求其交集問題。

    通過上述計算可以得到,四個單項約束條件下的經(jīng)濟增長率的共同區(qū)域為[103.66,105.13],見圖2。

    圖2 單項約束下的共同經(jīng)濟增長率區(qū)間

    六、研究結(jié)論

    本文通過建立狀態(tài)空間和目標(biāo)規(guī)劃模型,使用區(qū)間估計和集合求交集的方法,研究了環(huán)境約束下的適度經(jīng)濟增長率,得出如下結(jié)論:

    1.根據(jù)狀態(tài)空間計量模型估計結(jié)果,證明了環(huán)境與經(jīng)濟增長率是一種正向因果關(guān)系,經(jīng)濟增長率越高,環(huán)境污染將越大。這也證明了我們所提出的二者相克、相濟共生的理論觀點。

    2.在環(huán)境約束下,存在一個經(jīng)濟增長率的極大值,在兩種不同的方法測算下有不同的結(jié)果,多維約束條件下和單項條件約束下的極大值分別約為20%和13%。

    3.環(huán)境約束下的極大值不一定是適度值,分別用正態(tài)分布區(qū)間估計和集合論求交集的兩種方法測算的適度經(jīng)濟增長率都集中在4%左右,非常一致,將兩種方法的共同信息都考慮進去,可以確定中國現(xiàn)階段環(huán)境約束下的適度經(jīng)濟增長率區(qū)間在3%~5%范圍。超過這個范圍,可能就需要以犧牲環(huán)境為代價。

    4.現(xiàn)階段,中國不能單純追求高經(jīng)濟增長率,而必須考慮環(huán)境約束,這些約束主要包括環(huán)境對經(jīng)濟增長的彈性系數(shù)約束、單位GDP能源消耗約束、城鎮(zhèn)化速度約束和污染物排放約束。

    [1] Markus P.Technical Progress,Structural Change,and the Environmental Kuznets Curve[J].Ecological Economics,2002,42(3).

    [2] Dasgupta S,Laplante B,Wang H,et al.Confronting the Environmental Kuznets Curve[J].The Journal of Economic Perspectives,2002,16(1).

    [3] 胡宗義,劉亦文,唐偉杰.低碳經(jīng)濟背景下碳排放的庫茲涅茨曲線研究[J].統(tǒng)計研究,2013(2).

    [4] Panayotou T.Economic Growth and the Environment[J].Economic Survey of Europe,2003(2).

    [5] 梅國平,萬建香.環(huán)境約束下我國經(jīng)濟增長的內(nèi)生機理——基于CDE與FBA的數(shù)理分析與數(shù)字校正[J].經(jīng)濟管理,2012(3).

    [6] 岳書敬,劉富華.環(huán)境約束下的經(jīng)濟增長效率及其影響因素[J].?dāng)?shù)量經(jīng)濟技術(shù)經(jīng)濟研究,2009(5).

    Researching on Moderate Economic Growth under the Environmental Restriction in China:The Analysis Based on Exhaust Emissions

    LI Bao-yu,LIU Yang
    (School of Statistics,Shanxi University of Finance & Economics,Taiyuan 030006,China)

    The relationship between economic development and environmental pollution has always been one pair of unavoidable conflicts,which is also the academic research focus,the majority of domestic and foreign scholars have been paid their attention on the study of the environmental Kuznets curve,but the relationship between environment and economic growth always be ignored.In this paper,we firstly establish the time-varying parameter state space model of the relationship between environment and economic growth,and transformed the time-varying parameter state space model into a decision-making equation of linear programming,and then we apply the methods of interval estimation of parameters and set theory to estimate China's current optimal economic growth under the environmental constraints by constraining the variables and coefficient.At last we get the range of optimal economic growth which is between 3%and 5%.

    moderate economic growth;environmental restriction;state space model;linear programming model

    F205

    A

    1007-3116(2014)04-0023-07

    2013-11-15;修復(fù)日期:2014-01-17

    李寶瑜,男,山西文水人,教授,博士生導(dǎo)師,研究方向:國民經(jīng)濟核算與宏觀經(jīng)濟統(tǒng)計分析;

    劉 洋,女,山西介休人,碩士生,研究方向:國民經(jīng)濟核算與宏觀經(jīng)濟統(tǒng)計分析。

    (責(zé)任編輯:張治國)

    猜你喜歡
    城鎮(zhèn)化率約束條件增長率
    31 省份最新城鎮(zhèn)化率:9 省份超70%,這10 個省份城鎮(zhèn)人口最多
    基于一種改進AZSVPWM的滿調(diào)制度死區(qū)約束條件分析
    2020年河北省固定資產(chǎn)投資增長率
    2019年河北省固定資產(chǎn)投資增長率
    A literature review of research exploring the experiences of overseas nurses in the United Kingdom (2002–2017)
    國內(nèi)生產(chǎn)總值及其增長率
    線性規(guī)劃的八大妙用
    貨幣供應(yīng)量同比增長率
    安徽省生態(tài)足跡與城鎮(zhèn)化率關(guān)系實證研究
    2011年全省城鎮(zhèn)化率達45.6%
    如何舔出高潮| 嫩草影院入口| 久久精品国产亚洲av涩爱| 黄片播放在线免费| 亚洲欧美色中文字幕在线| 久久影院123| 免费观看av网站的网址| 婷婷色综合www| 国产在线免费精品| 亚洲,一卡二卡三卡| 黄色视频不卡| 另类精品久久| 99热网站在线观看| 久久久久久久久久久免费av| 性色av一级| 国产片内射在线| 免费观看人在逋| 亚洲伊人色综图| av在线播放精品| 这个男人来自地球电影免费观看 | 精品国产一区二区三区四区第35| 午夜日韩欧美国产| 国产黄色免费在线视频| 美女午夜性视频免费| 欧美人与性动交α欧美软件| 青草久久国产| 久久天堂一区二区三区四区| 亚洲四区av| 制服丝袜香蕉在线| 久久久国产欧美日韩av| 国产av国产精品国产| 美女主播在线视频| 国产一区二区三区综合在线观看| 999精品在线视频| 午夜久久久在线观看| 一级片免费观看大全| 欧美日韩av久久| 国产精品三级大全| 欧美精品一区二区大全| 极品少妇高潮喷水抽搐| 我的亚洲天堂| av在线观看视频网站免费| 老熟女久久久| 99久久综合免费| 亚洲av成人不卡在线观看播放网 | 亚洲成人免费av在线播放| 亚洲色图综合在线观看| 一级片免费观看大全| 蜜桃在线观看..| 中文字幕最新亚洲高清| 久久久久久久久久久久大奶| 999久久久国产精品视频| 高清av免费在线| 免费高清在线观看日韩| 嫩草影视91久久| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜美腿诱惑在线| 18禁观看日本| 免费高清在线观看视频在线观看| 亚洲人成网站在线观看播放| 国产精品一区二区精品视频观看| 国产精品 欧美亚洲| 国产极品天堂在线| 亚洲欧美日韩另类电影网站| 国产一区二区 视频在线| 欧美中文综合在线视频| av免费观看日本| 黄色一级大片看看| 激情视频va一区二区三区| 色网站视频免费| 97人妻天天添夜夜摸| www.熟女人妻精品国产| 中文字幕另类日韩欧美亚洲嫩草| 国产xxxxx性猛交| 精品国产国语对白av| 日韩人妻精品一区2区三区| 精品亚洲成a人片在线观看| 久久精品久久精品一区二区三区| 波多野结衣一区麻豆| 亚洲国产精品999| 一区二区三区激情视频| 国产精品免费大片| videosex国产| 超碰97精品在线观看| 国产精品女同一区二区软件| 久久精品熟女亚洲av麻豆精品| av福利片在线| 久久精品国产a三级三级三级| 丁香六月欧美| 青草久久国产| 男女床上黄色一级片免费看| 另类精品久久| 亚洲成人手机| av又黄又爽大尺度在线免费看| 你懂的网址亚洲精品在线观看| 国产亚洲精品第一综合不卡| 亚洲av电影在线观看一区二区三区| 一区福利在线观看| 亚洲精品国产区一区二| 丁香六月欧美| 亚洲av福利一区| 久久人妻熟女aⅴ| 亚洲欧美成人精品一区二区| 国产精品蜜桃在线观看| 9191精品国产免费久久| av又黄又爽大尺度在线免费看| 美女脱内裤让男人舔精品视频| 一区二区日韩欧美中文字幕| 午夜激情av网站| 黄片播放在线免费| 少妇的丰满在线观看| 天天操日日干夜夜撸| 欧美少妇被猛烈插入视频| 18在线观看网站| 国产精品.久久久| 国产无遮挡羞羞视频在线观看| 日日爽夜夜爽网站| 国产日韩欧美在线精品| 国产精品免费视频内射| 久久人人爽av亚洲精品天堂| 熟女少妇亚洲综合色aaa.| 香蕉丝袜av| 一个人免费看片子| 下体分泌物呈黄色| 啦啦啦在线免费观看视频4| 国产精品无大码| 九草在线视频观看| 999久久久国产精品视频| 美女福利国产在线| 日韩制服骚丝袜av| 免费黄网站久久成人精品| 好男人视频免费观看在线| 人成视频在线观看免费观看| 狂野欧美激情性bbbbbb| 久久久久网色| 欧美另类一区| 亚洲一码二码三码区别大吗| 亚洲一码二码三码区别大吗| 在线天堂中文资源库| 国产又爽黄色视频| 晚上一个人看的免费电影| 欧美亚洲日本最大视频资源| 在线观看免费高清a一片| 色综合欧美亚洲国产小说| 国产黄色免费在线视频| 国产成人精品久久久久久| 国产精品av久久久久免费| 少妇被粗大的猛进出69影院| 亚洲精品成人av观看孕妇| 久久ye,这里只有精品| 日本色播在线视频| 国产熟女欧美一区二区| 免费在线观看完整版高清| 国产深夜福利视频在线观看| 老司机亚洲免费影院| 国产一区二区在线观看av| 大片电影免费在线观看免费| 亚洲国产欧美在线一区| 免费观看a级毛片全部| 亚洲精品日韩在线中文字幕| 这个男人来自地球电影免费观看 | 老汉色av国产亚洲站长工具| 精品视频人人做人人爽| 亚洲欧美色中文字幕在线| 精品视频人人做人人爽| 熟妇人妻不卡中文字幕| 嫩草影视91久久| 精品一品国产午夜福利视频| avwww免费| 深夜精品福利| 在线看a的网站| 日韩不卡一区二区三区视频在线| 最黄视频免费看| 51午夜福利影视在线观看| 国产精品无大码| 国产亚洲欧美精品永久| 成人影院久久| 日韩av在线免费看完整版不卡| 爱豆传媒免费全集在线观看| 丝袜美足系列| 欧美日韩精品网址| 精品一区二区三卡| 久久毛片免费看一区二区三区| 岛国毛片在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 日日撸夜夜添| 久久久久久久精品精品| 国产精品一二三区在线看| 看免费成人av毛片| 免费观看av网站的网址| 操出白浆在线播放| 日韩成人av中文字幕在线观看| 亚洲精品日本国产第一区| 亚洲av成人不卡在线观看播放网 | 天天躁夜夜躁狠狠久久av| 日韩大片免费观看网站| 久久热在线av| 男人爽女人下面视频在线观看| 人人妻,人人澡人人爽秒播 | 日韩不卡一区二区三区视频在线| 国产精品.久久久| 18在线观看网站| av一本久久久久| 国产成人a∨麻豆精品| 女人被躁到高潮嗷嗷叫费观| 另类亚洲欧美激情| 啦啦啦 在线观看视频| 波多野结衣一区麻豆| 一级片'在线观看视频| 综合色丁香网| 2021少妇久久久久久久久久久| 亚洲精品在线美女| 一区福利在线观看| 在线天堂最新版资源| 国产成人精品无人区| 欧美日韩福利视频一区二区| 免费av中文字幕在线| 精品国产国语对白av| 亚洲四区av| 亚洲,欧美,日韩| 日日撸夜夜添| 久久精品久久久久久久性| 久久女婷五月综合色啪小说| 我的亚洲天堂| 热99久久久久精品小说推荐| 色综合欧美亚洲国产小说| 如何舔出高潮| 亚洲精品日韩在线中文字幕| 热99国产精品久久久久久7| 中文字幕人妻丝袜一区二区 | 国产日韩欧美在线精品| 亚洲国产av影院在线观看| 99热全是精品| 亚洲成色77777| 国产男女内射视频| 少妇精品久久久久久久| 99香蕉大伊视频| 下体分泌物呈黄色| 亚洲第一区二区三区不卡| 欧美日韩视频高清一区二区三区二| 搡老岳熟女国产| 一区二区三区精品91| 人妻人人澡人人爽人人| 国产探花极品一区二区| 成人影院久久| 一区二区日韩欧美中文字幕| 777米奇影视久久| 午夜精品国产一区二区电影| 又大又爽又粗| 国产精品久久久久成人av| 久久久久精品人妻al黑| 久久久精品94久久精品| videos熟女内射| 夜夜骑夜夜射夜夜干| 欧美变态另类bdsm刘玥| 国产97色在线日韩免费| 老熟女久久久| 观看av在线不卡| 亚洲欧美一区二区三区久久| 免费在线观看完整版高清| 国产成人欧美在线观看 | 亚洲伊人色综图| 国产极品天堂在线| 亚洲欧美清纯卡通| 欧美日韩综合久久久久久| 国产一区二区 视频在线| 国产色婷婷99| 日韩大片免费观看网站| 亚洲国产精品国产精品| 哪个播放器可以免费观看大片| 亚洲,一卡二卡三卡| 成年动漫av网址| 国产麻豆69| 免费看av在线观看网站| 美女大奶头黄色视频| 丰满少妇做爰视频| 亚洲av电影在线进入| 亚洲国产精品一区三区| 十八禁高潮呻吟视频| av免费观看日本| 亚洲精品第二区| 你懂的网址亚洲精品在线观看| 国精品久久久久久国模美| 1024香蕉在线观看| 久久久久国产一级毛片高清牌| 老汉色∧v一级毛片| 伊人亚洲综合成人网| 欧美日韩国产mv在线观看视频| 黄片小视频在线播放| 美女扒开内裤让男人捅视频| 国产色婷婷99| 777米奇影视久久| 亚洲国产精品国产精品| 午夜老司机福利片| 中文字幕亚洲精品专区| 99精品久久久久人妻精品| 青春草国产在线视频| 久久99一区二区三区| 久久精品国产综合久久久| 深夜精品福利| 人人澡人人妻人| 男人添女人高潮全过程视频| 可以免费在线观看a视频的电影网站 | 欧美成人午夜精品| 交换朋友夫妻互换小说| av不卡在线播放| av有码第一页| 国产精品一区二区在线不卡| 亚洲人成网站在线观看播放| 777米奇影视久久| 97在线人人人人妻| 国产精品秋霞免费鲁丝片| 老汉色∧v一级毛片| 十八禁人妻一区二区| 老司机在亚洲福利影院| 母亲3免费完整高清在线观看| 久久精品国产综合久久久| 岛国毛片在线播放| 国产成人av激情在线播放| 午夜免费男女啪啪视频观看| www.熟女人妻精品国产| 国产精品国产av在线观看| 国产精品人妻久久久影院| 极品少妇高潮喷水抽搐| 国产成人免费观看mmmm| 涩涩av久久男人的天堂| 最近最新中文字幕免费大全7| 天天躁夜夜躁狠狠躁躁| 亚洲欧洲国产日韩| 国产精品蜜桃在线观看| 免费在线观看黄色视频的| 国产精品免费视频内射| 国产女主播在线喷水免费视频网站| 一本久久精品| 狠狠婷婷综合久久久久久88av| 中文字幕av电影在线播放| 人妻一区二区av| 国产精品一区二区在线不卡| 尾随美女入室| 9色porny在线观看| 一本一本久久a久久精品综合妖精| 国产精品 欧美亚洲| avwww免费| 成人三级做爰电影| 欧美97在线视频| 十分钟在线观看高清视频www| 人人妻,人人澡人人爽秒播 | 亚洲一级一片aⅴ在线观看| 国产在线一区二区三区精| 少妇人妻 视频| 韩国高清视频一区二区三区| 如日韩欧美国产精品一区二区三区| 亚洲精品第二区| 久久免费观看电影| 97人妻天天添夜夜摸| 国产精品久久久久久人妻精品电影 | 巨乳人妻的诱惑在线观看| 视频区图区小说| 纯流量卡能插随身wifi吗| 亚洲一卡2卡3卡4卡5卡精品中文| 天天躁夜夜躁狠狠躁躁| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久人妻| 两性夫妻黄色片| 热re99久久精品国产66热6| 国产黄频视频在线观看| 老司机靠b影院| 亚洲天堂av无毛| 久久精品久久久久久久性| 大香蕉久久成人网| 久久精品久久久久久久性| 青草久久国产| 中文字幕制服av| 免费在线观看黄色视频的| 中文字幕制服av| 精品国产乱码久久久久久小说| 妹子高潮喷水视频| 日韩伦理黄色片| 久久久久久人妻| 久久久精品国产亚洲av高清涩受| av线在线观看网站| 国产成人欧美在线观看 | 深夜精品福利| 日日爽夜夜爽网站| 午夜福利在线免费观看网站| 午夜福利影视在线免费观看| 亚洲第一区二区三区不卡| 我的亚洲天堂| 看十八女毛片水多多多| 午夜福利影视在线免费观看| 亚洲一码二码三码区别大吗| 97精品久久久久久久久久精品| 欧美在线一区亚洲| 99热全是精品| 哪个播放器可以免费观看大片| 中文字幕亚洲精品专区| 国产一区二区在线观看av| 亚洲少妇的诱惑av| 亚洲综合色网址| 日韩电影二区| 国产精品欧美亚洲77777| 国产精品久久久人人做人人爽| svipshipincom国产片| 国产精品国产av在线观看| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 久久综合国产亚洲精品| 日韩,欧美,国产一区二区三区| 日日爽夜夜爽网站| 亚洲成人av在线免费| 中文字幕精品免费在线观看视频| av.在线天堂| 超碰成人久久| 1024视频免费在线观看| 欧美亚洲日本最大视频资源| 日本wwww免费看| 国产日韩欧美亚洲二区| 777米奇影视久久| 欧美黑人精品巨大| 青青草视频在线视频观看| 看免费成人av毛片| 天天躁日日躁夜夜躁夜夜| 亚洲第一区二区三区不卡| 9热在线视频观看99| 亚洲av中文av极速乱| 热re99久久精品国产66热6| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 自拍欧美九色日韩亚洲蝌蚪91| 日韩欧美一区视频在线观看| 国产高清国产精品国产三级| 十八禁网站网址无遮挡| 午夜91福利影院| 国产成人精品久久久久久| 水蜜桃什么品种好| 久久精品亚洲av国产电影网| 欧美 日韩 精品 国产| 久久这里只有精品19| 涩涩av久久男人的天堂| 午夜久久久在线观看| 久久亚洲国产成人精品v| 久久青草综合色| 十八禁网站网址无遮挡| 精品久久久精品久久久| 天天躁夜夜躁狠狠久久av| 久久精品久久久久久久性| 国产欧美日韩综合在线一区二区| 一个人免费看片子| 国产片内射在线| 精品人妻熟女毛片av久久网站| 中文字幕高清在线视频| 色网站视频免费| 男人爽女人下面视频在线观看| 国产黄色免费在线视频| 亚洲精品国产av蜜桃| 99香蕉大伊视频| 新久久久久国产一级毛片| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩视频精品一区| 国产精品蜜桃在线观看| 成人国语在线视频| 国产成人午夜福利电影在线观看| 看免费av毛片| 午夜福利免费观看在线| 亚洲国产精品成人久久小说| 1024视频免费在线观看| 国产欧美日韩一区二区三区在线| 成年女人毛片免费观看观看9 | 男女床上黄色一级片免费看| 亚洲国产欧美网| 亚洲av日韩在线播放| 国产日韩欧美视频二区| 亚洲欧美成人综合另类久久久| 精品亚洲成国产av| 亚洲欧洲精品一区二区精品久久久 | 交换朋友夫妻互换小说| 亚洲av综合色区一区| 少妇人妻久久综合中文| 国产淫语在线视频| 国产一区二区 视频在线| videosex国产| 亚洲欧美中文字幕日韩二区| 青草久久国产| 亚洲少妇的诱惑av| av一本久久久久| 国产成人精品福利久久| 久久久久久久国产电影| 国产97色在线日韩免费| 热99久久久久精品小说推荐| 美女福利国产在线| 国产精品久久久久久人妻精品电影 | 国产熟女欧美一区二区| 亚洲精品第二区| 国产欧美日韩综合在线一区二区| 美女午夜性视频免费| 久久综合国产亚洲精品| 精品国产露脸久久av麻豆| 国产精品熟女久久久久浪| 只有这里有精品99| 伊人久久国产一区二区| 国产精品三级大全| 精品酒店卫生间| 亚洲欧美一区二区三区久久| 亚洲激情五月婷婷啪啪| 中文字幕亚洲精品专区| 久久久国产一区二区| 99re6热这里在线精品视频| 丝袜美足系列| 日韩av不卡免费在线播放| 中文字幕制服av| 色综合欧美亚洲国产小说| 欧美黑人欧美精品刺激| 欧美乱码精品一区二区三区| 丝袜人妻中文字幕| 高清在线视频一区二区三区| 国产亚洲午夜精品一区二区久久| 岛国毛片在线播放| 曰老女人黄片| av免费观看日本| 在线观看一区二区三区激情| √禁漫天堂资源中文www| 不卡视频在线观看欧美| 日本欧美国产在线视频| 午夜av观看不卡| 欧美av亚洲av综合av国产av | 男女床上黄色一级片免费看| 中文字幕精品免费在线观看视频| 国产精品成人在线| av在线老鸭窝| 我的亚洲天堂| 51午夜福利影视在线观看| 久久久久国产精品人妻一区二区| 99香蕉大伊视频| 亚洲国产看品久久| 亚洲,一卡二卡三卡| 国产女主播在线喷水免费视频网站| 国产精品久久久久久久久免| 色吧在线观看| 久热这里只有精品99| 男男h啪啪无遮挡| 捣出白浆h1v1| 久久久久网色| 国产欧美日韩综合在线一区二区| 日韩视频在线欧美| 伦理电影免费视频| 亚洲久久久国产精品| 人妻人人澡人人爽人人| 一区二区三区精品91| 一个人免费看片子| 校园人妻丝袜中文字幕| 成人毛片60女人毛片免费| 亚洲精品中文字幕在线视频| 爱豆传媒免费全集在线观看| 老汉色av国产亚洲站长工具| 亚洲国产精品成人久久小说| 日韩精品有码人妻一区| 国产精品国产av在线观看| av国产久精品久网站免费入址| 日韩av不卡免费在线播放| 伊人久久大香线蕉亚洲五| 欧美少妇被猛烈插入视频| 夫妻性生交免费视频一级片| 亚洲欧美激情在线| 国产97色在线日韩免费| 国产精品二区激情视频| 99久国产av精品国产电影| 69精品国产乱码久久久| 日韩一区二区视频免费看| 嫩草影院入口| 久久99精品国语久久久| 丁香六月天网| 精品少妇久久久久久888优播| 国产福利在线免费观看视频| 欧美日韩亚洲高清精品| 国产成人啪精品午夜网站| 欧美变态另类bdsm刘玥| 一级片'在线观看视频| kizo精华| 亚洲av男天堂| 成人18禁高潮啪啪吃奶动态图| 午夜免费男女啪啪视频观看| 精品第一国产精品| 亚洲成人国产一区在线观看 | 99精品久久久久人妻精品| 免费黄网站久久成人精品| 桃花免费在线播放| 国产精品久久久av美女十八| 国产精品免费大片| 亚洲精华国产精华液的使用体验| 国产精品人妻久久久影院| 日韩欧美一区视频在线观看| 亚洲国产精品一区三区| 啦啦啦 在线观看视频| 国产成人a∨麻豆精品| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边抽搐一进一出视频| 亚洲精品日韩在线中文字幕| 丝袜在线中文字幕| 久久久久网色| 1024香蕉在线观看| 亚洲成国产人片在线观看| 日韩电影二区| 国产精品女同一区二区软件| 色94色欧美一区二区| www日本在线高清视频| 欧美亚洲 丝袜 人妻 在线| 中文字幕亚洲精品专区| 日韩大片免费观看网站| 欧美黑人精品巨大| 国产成人a∨麻豆精品| 欧美亚洲 丝袜 人妻 在线| 国产xxxxx性猛交| 午夜日韩欧美国产| 国产一区二区三区av在线| 人妻 亚洲 视频|