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

    門限CARR模型的實(shí)證研究
    ——以上證市場(chǎng)為例

    2015-02-20 05:23:56鄒文俊
    關(guān)鍵詞:極差門限限值

    鄒文俊

    門限CARR模型的實(shí)證研究
    ——以上證市場(chǎng)為例

    鄒文俊

    (天津大學(xué) 管理與經(jīng)濟(jì)學(xué)部,天津 300072)

    金融數(shù)據(jù)的波動(dòng)往往存在著非線性的特征,為了更好地刻畫上證市場(chǎng)股價(jià)指數(shù)的波動(dòng)特征,將門限模型與CARR模型相結(jié)合,構(gòu)建了門限CARR模型。首先借鑒門限自回歸條件持續(xù)期模型(TACD)的建模方法,構(gòu)建門限CARR模型;其次選取2002年1月4日至2012年2月15日的上證綜合股價(jià)指數(shù)的日數(shù)據(jù)進(jìn)行研究,通過非線性檢驗(yàn)和門限值識(shí)別的方法找到樣本序列的門限值;最后建立了擬合上證股票市場(chǎng)波動(dòng)的非線性特征的門限CARR模型,得到上證股票市場(chǎng)在高、低波動(dòng)下的不同波動(dòng)特征。

    非線性;門限模型;CARR模型;波動(dòng)

    針對(duì)金融數(shù)據(jù)的非線性特征的研究是金融時(shí)間序列分析中的重要研究領(lǐng)域,對(duì)于非線性時(shí)間序列模型的研究也逐漸成為學(xué)術(shù)研究的熱點(diǎn)。常見的非線性模型有非線性自回歸模型(NLAR)、廣義自回歸模型(GAR)、雙線性模型(BL)、門限自回歸模型(TAR)、雙門限模型、平滑轉(zhuǎn)移模型(STAR)、向量平滑轉(zhuǎn)移模型(STVAR模型)、人工神經(jīng)網(wǎng)絡(luò)模型(ANN)、馬爾科夫體制轉(zhuǎn)換模型等。

    在以上幾種常見的非線性模型中被廣泛使用的主要有兩種類型,一類是門限模型;另一類是馬爾科夫機(jī)制轉(zhuǎn)換模型。Hamilton(1989)[1]提出了馬爾科夫機(jī)制轉(zhuǎn)換模型(Markov regime switching model),此后馬爾科夫機(jī)制轉(zhuǎn)換模型在時(shí)間序列分析上被廣泛應(yīng)用。馬爾科夫機(jī)制轉(zhuǎn)換模型只能推斷出變量在各機(jī)制之間的轉(zhuǎn)換概率,而無法描繪出機(jī)制轉(zhuǎn)換的非線性形式。門限模型(Threshold models)由Tong(1978)[2]提出,與馬爾科夫機(jī)制轉(zhuǎn)換模型相比,門限模型以變量本身作為結(jié)構(gòu)轉(zhuǎn)折點(diǎn),利用類似分段建模的形式采用門限空間來改進(jìn)線性逼近,能夠更加清晰地描述時(shí)間序列的機(jī)制轉(zhuǎn)換運(yùn)動(dòng)狀態(tài)。

    而對(duì)于時(shí)間序列波動(dòng)性的刻畫,Chou(2005)[3]提出了CARR(Conditional Auto-Regressive Range)模型,即自回歸條件極差模型,將極差與GARCH模型的思想結(jié)合起來,有效刻畫極差的動(dòng)態(tài)結(jié)構(gòu)。使得極差在波動(dòng)性的預(yù)測(cè)上,比已被廣泛接受的GARCH模型更具有優(yōu)勢(shì)。丁忠民等(2005)[4]以上海股市作為研究對(duì)象,分別運(yùn)用CARR模型和GARCH模型進(jìn)行波動(dòng)性預(yù)測(cè),進(jìn)而對(duì)兩種方法的預(yù)測(cè)能力進(jìn)行比較,實(shí)證結(jié)果表明CARR模型在擬合波動(dòng)性方面優(yōu)于GARCH模型。程細(xì)玉等(2009)[5]以2006年美元對(duì)港幣5分鐘匯率比價(jià)和日元對(duì)美元15分鐘匯率比價(jià)的極差和收益率作為樣本數(shù)據(jù)進(jìn)行實(shí)證分析,研究結(jié)果表明CARR模型及其衍生模型在高頻金融時(shí)間序列的價(jià)格波動(dòng)性捕捉方面表現(xiàn)突出。

    然而,傳統(tǒng)CARR模型和GARCH模型無法很好地?cái)M合存在變結(jié)構(gòu)現(xiàn)象的非線性金融時(shí)間序列。對(duì)于存在不同波動(dòng)體制的波動(dòng)序列還需要引進(jìn)門限模型的思想。本文將門限模型與CARR模型相結(jié)合建模,構(gòu)建門限CARR模型,選取2002年1月4日至2012年2月15日的上證綜合股價(jià)指數(shù)的日數(shù)據(jù)進(jìn)行研究,分別為上證市場(chǎng)股價(jià)指數(shù)的高波動(dòng)和低波動(dòng)建模,刻畫上證市場(chǎng)在不同波動(dòng)體制下的波動(dòng)特征。

    1 門限CARR模型的構(gòu)建

    1.1 CARR模型簡(jiǎn)介

    傳統(tǒng)的GARCH模型能夠很好地解決金融時(shí)間序列中出現(xiàn)的異方差問題,對(duì)于金融資產(chǎn)的波動(dòng)持續(xù)性、聚集性和非對(duì)稱性能夠較好地描述, 所以GARCH模型方法在經(jīng)濟(jì)金融領(lǐng)域被廣泛接受。在此基礎(chǔ)上,Chou(2005)[3]提出了CARR(Conditional Auto-Regressive Range)模型,將極差與GARCH模型的思想結(jié)合起來,有效刻畫極差的動(dòng)態(tài)結(jié)構(gòu),使得極差在波動(dòng)性的預(yù)測(cè)上,比已被廣泛接受的GARCH模型更具有優(yōu)勢(shì)。GARCH模型通過收益率數(shù)據(jù)來估計(jì)波動(dòng)性,CARR模型則是通過極差數(shù)據(jù)來估計(jì)波動(dòng)率,相對(duì)于收益率所包含的資產(chǎn)價(jià)格變動(dòng)信息不足的弱點(diǎn),極差在反映資產(chǎn)價(jià)格變動(dòng)方面占有優(yōu)勢(shì)。極差利用了單位交易區(qū)間中資產(chǎn)的最高價(jià)和最低價(jià)兩個(gè)時(shí)間點(diǎn)上的信息,更能反映單位交易區(qū)間內(nèi)的波動(dòng)情況。

    就模型的形式及性質(zhì)而言,CARR模型類似于GARCH模型,同時(shí)也類似于Engle等(1998)[6]所提出的ACD(Autoregressive Conditional Duration)模型。CARR(p,q)模型如下式所示:

    Rt=λtεt

    (1)

    (2)

    εt~iidf(·)

    (3)

    (4)

    式中:Rt表示第t期股票價(jià)格Vt取自然對(duì)數(shù)后的極差;λt表示在t期之前的所有信息集合It-1已知的情況下,對(duì)第t期極差Rt進(jìn)行預(yù)測(cè)的條件期望,λt=E(Rt/It-1),λt≥0;干擾項(xiàng)εt是標(biāo)準(zhǔn)化極差εt=Rt/λt,不同時(shí)期的干擾項(xiàng)相互獨(dú)立,它的分布被假定為服從一個(gè)有單位均值的密度函數(shù)f(·)的分布;ω表示極差中的不確定性,也可以表示極差的初始水平,ω >0;αi是極差的滯后期系數(shù),可以表示極差條件均值的短期影響,且αi≥0,i=1,2,…,p;βj是極差條件均值的滯后期系數(shù),可以表示極差條件均值的長(zhǎng)期影響,且βj≥0,j=1,2,…,q。

    Chou提出根據(jù)CARR模型干擾項(xiàng)服從不同分布,CARR模型可分為不同種類模型。當(dāng)干擾項(xiàng)εt服從單位指數(shù)分布時(shí),CARR模型的對(duì)數(shù)似然函數(shù)形式為

    (5)

    干擾項(xiàng)服從上述指數(shù)分布的CARR模型稱為ECARR模型。

    當(dāng)干擾項(xiàng)εt服從Weibull分布,即εt~Gamma(γ,1/γ)時(shí),CARR模型的對(duì)數(shù)似然函數(shù)形式為:

    (6)

    干擾項(xiàng)服從上述Weibull分布的CARR模型稱為WCARR模型。

    1.2 構(gòu)建門限CARR模型

    CARR模型較好地?cái)M合了金融資產(chǎn)的波動(dòng)性特征,而在現(xiàn)實(shí)情況中,金融資產(chǎn)的波動(dòng)往往會(huì)出現(xiàn)不同波動(dòng)體制的變結(jié)構(gòu)現(xiàn)象,因此構(gòu)建非線性的CARR模型十分必要。本文構(gòu)建門限CARR模型的思路借鑒了門限自回歸條件持續(xù)期模型(TACD)的建模方法。Zhang等(2001)[7]提出了依靠門限變量來描述復(fù)雜經(jīng)濟(jì)過程的ACD模型,成為TACD模型。首先定義

    Rj∈[rj-1,rj),j=1,2,…,J

    (7)

    式中:J為正整數(shù);-∞=r0

    門限變量定義為

    Zi=h(xi,…,x1;yi,…,y1)

    (8)

    式中{yi}為一組與{xi}相關(guān)的經(jīng)濟(jì)變量。

    當(dāng)Zi-1∈Rj

    (9)

    依據(jù)TACD模型的構(gòu)建方法,本文將構(gòu)建門限條件自回歸極差模型即TCARR模型。將極差作為研究變量,門限變量為極差Rt的一階滯后項(xiàng),定義門限空間為Sj,則TCARR模型的基本形式為

    (10)

    (11)

    兩體制TCARR模型TCARR(2,p,q)可以較好地反映極差波動(dòng)的非線性特征,同時(shí)將金融資產(chǎn)波動(dòng)分為高波動(dòng)和低波動(dòng)兩種體制,能夠反映出更多的市場(chǎng)微觀結(jié)構(gòu)信息。

    2 非線性檢驗(yàn)

    門限模型是非線性時(shí)間序列模型,所以在構(gòu)建門限模型之前,首先要對(duì)樣本進(jìn)行檢驗(yàn),考察其是否具有非線性的特征。金融時(shí)間序列的非線性特性檢驗(yàn)方法有Hurst指數(shù)法、BDS檢驗(yàn)法、替代數(shù)據(jù)法和基于Ljung-Box Q統(tǒng)計(jì)量的Mcleod-Li檢驗(yàn)法。BDS檢驗(yàn)?zāi)苡行У貦z測(cè)出其他統(tǒng)計(jì)檢驗(yàn)所忽視的非線性特征,因而在非線性檢驗(yàn)中具有很強(qiáng)的功效性,應(yīng)用較為廣泛。本文簡(jiǎn)要介紹BDS檢驗(yàn)方法。

    BDS檢驗(yàn)是由Brock等(1996)[8]正式提出的。BDS檢驗(yàn)運(yùn)用了關(guān)聯(lián)積分法,零假設(shè)為所檢驗(yàn)的時(shí)間序列是獨(dú)立同分布的。在進(jìn)行檢驗(yàn)之前,先確定區(qū)間大小ε。如果時(shí)刻s,t的觀察值為ys和yt,則所有的觀察值(ys,yt)按對(duì)構(gòu)建為

    {(ys,yt),(ys+1,yt+1),…,(ys+m,yt+m)}

    (12)

    式中m為嵌入?yún)^(qū)間。每對(duì)觀察值滿足ε條件的共有概率被定義為概率cm(ε)。在獨(dú)立同分布的零假設(shè)下共有概率cm(ε)是每對(duì)概率c1(ε),c2(ε),…的簡(jiǎn)化:

    (13)

    如果需要觀察的有n個(gè)樣本,那么就要通過滿足ε條件的數(shù)對(duì)的數(shù)量和所有被觀察數(shù)對(duì)來估計(jì):

    (14)式中Iε為指示函數(shù):

    (15)

    這里的cm,n也可以作為關(guān)聯(lián)積分項(xiàng)來定義。通過關(guān)聯(lián)積分項(xiàng),BDS檢驗(yàn)統(tǒng)計(jì)量被定義為

    (16)

    在零假設(shè)為“時(shí)間序列yt為獨(dú)立同分布”時(shí),BDS檢驗(yàn)統(tǒng)計(jì)值近似于一個(gè)標(biāo)準(zhǔn)正態(tài)分布:BDS(n,m,ε)→N(0,1)

    如果BDS檢驗(yàn)統(tǒng)計(jì)量是正值,即cm,n(ε)-[c1,n-m+1(ε)]m>0,說明“積聚效應(yīng)”在價(jià)格波動(dòng)中要比純粹的隨機(jī)過程中發(fā)生的頻率要高。一個(gè)正的BDS檢驗(yàn)統(tǒng)計(jì)量拒絕了殘差為獨(dú)立正態(tài)同分布的假設(shè)。在實(shí)際操作中還應(yīng)該注意ε和“嵌入?yún)^(qū)間”m的選擇范圍。Scheinkman等(1989)[9]指出,當(dāng)ε取值過小,可能會(huì)導(dǎo)致接受零假設(shè);而如果ε過大,那么就會(huì)存在對(duì)數(shù)據(jù)過高估計(jì)的危險(xiǎn)。Hsieh(1989)[10]則建議“嵌入?yún)^(qū)間”m應(yīng)該從2到10的區(qū)間里選取。m越大,相關(guān)函數(shù)被高估的風(fēng)險(xiǎn)就越大。文中根據(jù)樣本個(gè)數(shù),m取值為6.

    3 門限值識(shí)別方法

    在建立門限CARR(TCARR)模型時(shí),本文選用極差的一階滯后項(xiàng)Rt-1作為模型的門限變量。接下來就要設(shè)定門限值,由于門限值是未知的,需要人為地進(jìn)行尋找。Chan(1993)[11]提出了一種獲得門限值的超一致估計(jì)量的方法。本文運(yùn)用這種方法尋找門限CARR模型的門限值。其估計(jì)方法為:

    首先,對(duì)樣本中所有的極差觀察值從小到大進(jìn)行排序,即:

    R1

    對(duì)于每個(gè)Rj,令τ=Rj,首先將其作為門限值進(jìn)行門限CARR模型的參數(shù)估計(jì)。門限值的一致估計(jì)量就是該模型的殘差最小時(shí)所采用的樣本值。

    而在實(shí)際實(shí)證分析中,通常不考慮樣本中最大的15%和最小的15%的{Rj},使門限的兩側(cè)都存在足夠多的觀察值。將排序后得到的新序列(記其樣本共有T′個(gè)變量)的第一個(gè)值作為門限值帶入模型進(jìn)行門限模型參數(shù)的估計(jì),同時(shí)需要求出此時(shí)的模型的殘差平方和,記做resid(1);然后,再將排序后序列中的每個(gè)樣本值作為門限值進(jìn)行門限CARR模型的估計(jì),得到殘差平方和序列{resid(1)},i=1,2,…,t′。該序列中的最小值resid(k)(k≤t′)所對(duì)應(yīng)的極差序列樣本值Rk即為門限CARR模型的門限值。

    4 實(shí)證研究

    4.1 數(shù)據(jù)描述

    本文選取2002年1月4日至2012年2月15日的上證綜合股價(jià)指數(shù)的日數(shù)據(jù)進(jìn)行研究,經(jīng)過整理之后總共有2369個(gè)極差觀察值。本文首先利用上證綜合股價(jià)指數(shù)極差序列所擬合的CARR模型的殘差序列進(jìn)行BDS檢驗(yàn),以判斷上證股票市場(chǎng)的極差序列是否具有非線性特征。

    表1 上證綜合指數(shù)的BDS檢驗(yàn)

    4.2 門限值的設(shè)定

    由4.1得出上證綜合指數(shù)的極差序列具有非線性特征的結(jié)論,在此基礎(chǔ)上,本文建立門限CARR模型。首先要確定門限值。按照前面介紹的門限值識(shí)別方法,本文將原極差序列樣本進(jìn)行從小到大的排序,并排除前15%和后15%的樣本得到樣本容量為1658個(gè)的新序列。然后將新序列中的每個(gè)樣本值一次帶入門限CARR(2,1,1)模型中并計(jì)算其殘差平方和,得到殘差平方和序列,見圖1.

    圖1 上證綜合指數(shù)殘差平方和序列

    從殘差平方和序列分布圖并結(jié)合樣本序列,得出上海綜合指數(shù)極差序列的第1211個(gè)樣本擬合的門限CARR(2,1,1)模型的殘差平方和最小,得到其模型的門限值為2.744642603。

    4.3 建立兩體制門限CARR(2,1,1)模型

    由4.2得出了擬合上海綜合指數(shù)極差序列的兩體制門限CARR(2,1,1)模型的門限值,將Rt-1作為門限變量,建立擬合上證股票市場(chǎng)波動(dòng)的非線性特征的門限CARR模型。為了更好地說明門限模型在擬合非線性金融數(shù)據(jù)中的作用,本文將新構(gòu)建的兩體制CARR(2,1,1)模型與傳統(tǒng)CARR(1,1)模型分別擬合樣本數(shù)據(jù),得到的實(shí)證結(jié)果如下:

    描述上證綜合指數(shù)極差序列的兩體制門限CARR(2,1,1)模型(τ=2.744642603)數(shù)據(jù)見表2。

    表2 兩體制門限CARR(2,1,1)模型數(shù)據(jù)

    描述上證綜合指數(shù)極差序列的CARR(1,1)模型數(shù)據(jù)見表3.

    表3 CARR(1,1)模型擬合數(shù)據(jù)

    表2中,各均值方程各參數(shù)值旁的方括號(hào)和圓括號(hào)內(nèi)的數(shù)據(jù)為其參數(shù)的T統(tǒng)計(jì)量及P值。

    在95%的置信水平下,T(40)=1.684,而T(2369)

    表3中,ω、β和γ的P值在95%的置信水平下不能推翻參數(shù)為零的假設(shè),故參數(shù)不顯著。

    從以上的實(shí)驗(yàn)結(jié)果可以看出,兩體制門限CARR(2,1,1)模型相對(duì)于傳統(tǒng)CARR(1,1)模型,對(duì)于上海股票市場(chǎng)數(shù)據(jù)擬合程度更好,并且對(duì)于上海股票市場(chǎng)在低波動(dòng)和高波動(dòng)體制下的波動(dòng)特征進(jìn)行了有效地刻畫,其對(duì)數(shù)似然函數(shù)值達(dá)到了-2924.18841599,較好地體現(xiàn)了上海股票市場(chǎng)波動(dòng)的非線性特征。

    圖2為上證綜合指數(shù)波動(dòng)在不同體制下的模型的殘差項(xiàng)序列

    圖2 上證綜合指數(shù)波動(dòng)處于兩種體制下的

    從圖2中可以明顯看出,處于不同體制下的波動(dòng)序列存在明顯差別,在高波動(dòng)體制下波動(dòng)特征更加明顯,這充分說明了兩體制門限CARR(2,1,1)模型可以非常明顯地刻畫出上證股票市場(chǎng)在不同波動(dòng)體制下的波動(dòng)特征。

    從上述實(shí)證結(jié)果可看出,我國(guó)上證股票市場(chǎng)存在不同波動(dòng)體制,這與現(xiàn)實(shí)情況十分契合。在現(xiàn)實(shí)情況中,由于金融危機(jī),信息不對(duì)稱,莊家操縱股價(jià),投資者不理性等問題,股市會(huì)產(chǎn)生劇烈波動(dòng),這時(shí)如果仍用單一波動(dòng)體制來刻畫金融時(shí)間序列就會(huì)刻畫不充分。將變結(jié)構(gòu)的非線性模型與傳統(tǒng)刻畫波動(dòng)的模型相結(jié)合,能夠更加準(zhǔn)確地刻畫劇烈波動(dòng)的金融市場(chǎng)。

    5 結(jié)束語

    本文將門限模型與CARR模型相結(jié)合構(gòu)建門限CARR模型,利用構(gòu)建的門限CARR模型對(duì)上證市場(chǎng)股價(jià)指數(shù)的波動(dòng)進(jìn)行了研究。在實(shí)證研究中,兩體制門限CARR(2,1,1)模型很好地描述了上證市場(chǎng)股價(jià)指數(shù)在高波動(dòng)和低波動(dòng)時(shí)的不同特征。本文嘗試用非線性模型中的門限模型與刻畫波動(dòng)的CARR模型相結(jié)合來擬合金融市場(chǎng)數(shù)據(jù),在未來的研究中,可以擴(kuò)展到更加廣泛的非線性模型的應(yīng)用,比如馬爾科夫體制轉(zhuǎn)換模型等。

    [1] Hamilton J D,Susmel R.Autoregressive conditional heteroskedasticity and changes in regime[J].Journal of Econometrics,1994,(64):307-333.

    [2]Tong H.On a Threshold Model[M].Pattern Recognition and Signal Processing,Amsterdan:Sijthoff&Noordhoff,1978:101-141.

    [3]Chou R Y.Forecasting financial volatilities with extreme values:the conditional autoregressive range (CARR) model[J].Journal of Money’ Credit and Banking,2005,37(3):561-582.

    [4]丁忠明,夏萬軍.中國(guó)股市波動(dòng)的CARR模型分析[J].商業(yè)經(jīng)濟(jì)與管理,2005,(12):41-45.

    [5]程細(xì)玉,夏天.金融市場(chǎng)波動(dòng)性CARR類模型與GARCH類模型的比較研究[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2009,(13):12-18.

    [6]Engle R F,Russell J F.Autoregressive Conditional Duration:a New Model for Irregular Spaced Transaction Data[J].Econometrica,1998,66(5):1127-1162.

    [7]Zhang M Y,Russell J,Tsay R.A nonlinear autoregressive conditonal duration model with applications to financial transactions data[J].Journal of Econometrics,2001,(104):179-207.

    [8]Brock W A,Dechert W D.LeBaron Scheinkman J.A.A Test for Independence Based on the Correlation Dimension [J].Economic Review,1996,(15):197-235.

    [9]Scheinkman J A,LeBaron B.Nonlinear Dynamics and Stock Returns [J].The Journal of Business,1989,62 (3):311-337.

    [10]Hsieh D A.Modeling Heteroscedasticity in Daily Foreign-Exchange Rates[J].Journal of Business & Economic Statistics,1989,7(3):307-317.

    [11]Chan K S.Consistency and limiting distribution of the least squares estimator of a threshold autoregressive model [J].Ann.Statist ,1993,(21):520-533.

    (責(zé)任編輯:馬金發(fā))

    Empirical Study on the Threshold CARR Model:An Example from Shanghai Stock Market

    ZOU Wenjun

    (College of Management and Economics,Tianjin University,Tianjin 300072,China)

    The volatility of financial data often has nonlinear characteristics.In order to describe the volatility characteristics of Shanghai stock market,we combined the threshold model and CARR (Conditional Auto-Regressive Range) model to construct the threshold CARR model.First,we borrowed the lessons from the threshold autoregressive conditional duration model (TACD) to construct the threshold CARR model;Secondly,we used the daily data of Shanghai stock market from January 4th,2002 to February 15th,2012 as the sample and find the threshold value through the nonlinear test and threshold recognition method;At last,we established the threshold CARR model which fit the nonlinear characteristics of the volatility of Shanghai stock market and describe the different volatility characteristics of Shanghai stock market.

    nonlinear characteristics;threshold model;CARR model;volatility

    2014-09-23

    國(guó)家自然科學(xué)基金資助項(xiàng)目(70901055);國(guó)家社會(huì)科學(xué)基金資助項(xiàng)目(14CTJ012);教育部博士點(diǎn)新教師基金資助項(xiàng)目(20090032120031);天津大學(xué)自主創(chuàng)新基金項(xiàng)目

    鄒文俊(1990—),女,碩士研究生,研究方向:金融時(shí)間序列分析.

    1003-1251(2015)02-0078-06

    F830

    A

    猜你喜歡
    極差門限限值
    基于規(guī)則的HEV邏輯門限控制策略
    地方債對(duì)經(jīng)濟(jì)增長(zhǎng)的門限效應(yīng)及地區(qū)差異研究
    隨機(jī)失效門限下指數(shù)退化軌道模型的分析與應(yīng)用
    二次函數(shù)極差的規(guī)律與應(yīng)用
    關(guān)于廢水排放特別限值的思考
    遼寧省遼河流域石油煉制排放限值的制定
    在分位數(shù)回歸中結(jié)構(gòu)突變的經(jīng)驗(yàn)極差檢驗(yàn)
    中美煉鋼行業(yè)污染物排放限值研究
    生產(chǎn)性服務(wù)業(yè)集聚與工業(yè)集聚的非線性效應(yīng)——基于門限回歸模型的分析
    湖湘論壇(2015年3期)2015-12-01 04:20:17
    環(huán)境保護(hù)部解讀新發(fā)布的大氣污染物特別排放限值
    久久国产精品影院| 肉色欧美久久久久久久蜜桃| 一区在线观看完整版| av一本久久久久| 亚洲欧美一区二区三区久久| 老熟女久久久| 黄色视频不卡| 久久午夜亚洲精品久久| 精品第一国产精品| 999精品在线视频| 国产精品久久久人人做人人爽| 肉色欧美久久久久久久蜜桃| 777米奇影视久久| 亚洲av成人一区二区三| 新久久久久国产一级毛片| 日韩三级视频一区二区三区| 国产极品粉嫩免费观看在线| 国精品久久久久久国模美| 黄网站色视频无遮挡免费观看| 国产不卡一卡二| 日韩视频一区二区在线观看| e午夜精品久久久久久久| 真人做人爱边吃奶动态| 亚洲欧美精品综合一区二区三区| 久久精品国产亚洲av香蕉五月 | 十分钟在线观看高清视频www| 亚洲精品av麻豆狂野| 亚洲精品中文字幕在线视频| 99在线人妻在线中文字幕 | 久久热在线av| 黄色 视频免费看| 99精品欧美一区二区三区四区| 欧美日本中文国产一区发布| 欧美日韩视频精品一区| 91精品国产国语对白视频| 婷婷成人精品国产| 岛国毛片在线播放| av福利片在线| 日日摸夜夜添夜夜添小说| 另类精品久久| 亚洲精品国产精品久久久不卡| 两性夫妻黄色片| 国产精品一区二区在线不卡| 少妇 在线观看| 成年动漫av网址| 亚洲av第一区精品v没综合| 夜夜夜夜夜久久久久| 亚洲中文日韩欧美视频| tube8黄色片| 欧美激情高清一区二区三区| 一本—道久久a久久精品蜜桃钙片| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| 国产精品av久久久久免费| av片东京热男人的天堂| 99九九在线精品视频| 黑人巨大精品欧美一区二区蜜桃| 俄罗斯特黄特色一大片| 2018国产大陆天天弄谢| 久久国产精品影院| 欧美日本中文国产一区发布| 精品国产乱码久久久久久小说| 国精品久久久久久国模美| 日韩免费av在线播放| 亚洲熟妇熟女久久| 99久久精品国产亚洲精品| 亚洲av片天天在线观看| 女人高潮潮喷娇喘18禁视频| 纵有疾风起免费观看全集完整版| 久久热在线av| 99国产极品粉嫩在线观看| 精品久久久久久久毛片微露脸| 亚洲国产av新网站| 欧美精品啪啪一区二区三区| 蜜桃在线观看..| 大香蕉久久成人网| 黄色毛片三级朝国网站| 亚洲欧美色中文字幕在线| 国产伦人伦偷精品视频| 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 国产91精品成人一区二区三区 | 国产精品影院久久| 99国产精品一区二区蜜桃av | 欧美日韩av久久| av欧美777| 最近最新免费中文字幕在线| 久久久精品94久久精品| 久久国产精品影院| 国产精品.久久久| 视频区图区小说| 黑人操中国人逼视频| 97在线人人人人妻| 久久久久久久精品吃奶| 黑人欧美特级aaaaaa片| 女人高潮潮喷娇喘18禁视频| 久久免费观看电影| 一级毛片女人18水好多| videos熟女内射| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 国产日韩欧美视频二区| 久久国产精品大桥未久av| 在线av久久热| 一级黄色大片毛片| 久久久久精品人妻al黑| 国产成人一区二区三区免费视频网站| 9热在线视频观看99| 一区二区日韩欧美中文字幕| 男人操女人黄网站| 美国免费a级毛片| 欧美精品高潮呻吟av久久| 高清黄色对白视频在线免费看| 成年女人毛片免费观看观看9 | 少妇粗大呻吟视频| 视频区欧美日本亚洲| 国产精品美女特级片免费视频播放器 | 婷婷成人精品国产| 成人亚洲精品一区在线观看| 欧美国产精品一级二级三级| 亚洲国产成人一精品久久久| 免费高清在线观看日韩| svipshipincom国产片| 精品一区二区三区视频在线观看免费 | 午夜日韩欧美国产| 国产精品影院久久| 久久久久国内视频| 大片免费播放器 马上看| 精品亚洲成国产av| 免费观看av网站的网址| 最近最新中文字幕大全电影3 | 亚洲美女黄片视频| a级毛片在线看网站| 久久人妻av系列| 国产高清视频在线播放一区| 亚洲av第一区精品v没综合| 欧美国产精品一级二级三级| 天天躁夜夜躁狠狠躁躁| 日韩人妻精品一区2区三区| 国产精品一区二区在线不卡| 色在线成人网| 中国美女看黄片| 国产99久久九九免费精品| 成人国产一区最新在线观看| 亚洲欧美一区二区三区黑人| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 久久亚洲真实| 日日摸夜夜添夜夜添小说| 波多野结衣av一区二区av| 电影成人av| 精品少妇久久久久久888优播| 91字幕亚洲| 天天添夜夜摸| svipshipincom国产片| 悠悠久久av| 黑丝袜美女国产一区| av又黄又爽大尺度在线免费看| 一级黄色大片毛片| av有码第一页| 国产成人啪精品午夜网站| 动漫黄色视频在线观看| 国产亚洲欧美在线一区二区| 法律面前人人平等表现在哪些方面| 一区二区三区激情视频| 国产成人精品在线电影| 多毛熟女@视频| 成人av一区二区三区在线看| 亚洲av日韩在线播放| 丰满人妻熟妇乱又伦精品不卡| 少妇粗大呻吟视频| 亚洲精华国产精华精| 国产成人av激情在线播放| 99久久人妻综合| 天天操日日干夜夜撸| 两个人看的免费小视频| 日韩成人在线观看一区二区三区| tocl精华| netflix在线观看网站| 美女午夜性视频免费| 亚洲人成电影观看| 看免费av毛片| 久久精品aⅴ一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 日韩免费高清中文字幕av| 蜜桃在线观看..| 国产成人一区二区三区免费视频网站| 国产精品一区二区精品视频观看| 性高湖久久久久久久久免费观看| 国产aⅴ精品一区二区三区波| 男人舔女人的私密视频| 午夜福利在线观看吧| 国产99久久九九免费精品| 国产成人系列免费观看| 久久久国产成人免费| 美国免费a级毛片| 亚洲成人国产一区在线观看| 日本wwww免费看| 国产精品久久久久久精品电影小说| 亚洲男人天堂网一区| aaaaa片日本免费| 黄色视频在线播放观看不卡| 亚洲av日韩精品久久久久久密| 午夜激情久久久久久久| 精品人妻1区二区| 又大又爽又粗| 超碰成人久久| 高清视频免费观看一区二区| 亚洲国产中文字幕在线视频| 99国产极品粉嫩在线观看| 1024视频免费在线观看| av电影中文网址| 国产精品国产av在线观看| 99香蕉大伊视频| 不卡av一区二区三区| 老汉色av国产亚洲站长工具| 在线播放国产精品三级| 国产精品久久久av美女十八| 在线观看免费午夜福利视频| 亚洲欧美激情在线| 日韩欧美免费精品| 亚洲欧美日韩高清在线视频 | 国产人伦9x9x在线观看| 国产xxxxx性猛交| 久久午夜亚洲精品久久| 国产精品免费视频内射| 黄频高清免费视频| 丁香六月欧美| 久久久精品免费免费高清| 亚洲国产av新网站| 999精品在线视频| cao死你这个sao货| 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频| 成人国语在线视频| 久久久久久人人人人人| 夜夜爽天天搞| 欧美日韩国产mv在线观看视频| 黄片小视频在线播放| 久久精品国产99精品国产亚洲性色 | 精品乱码久久久久久99久播| 国产精品久久电影中文字幕 | 他把我摸到了高潮在线观看 | 1024香蕉在线观看| 精品午夜福利视频在线观看一区 | 汤姆久久久久久久影院中文字幕| 国产老妇伦熟女老妇高清| 国产亚洲欧美在线一区二区| 色视频在线一区二区三区| 国产成人系列免费观看| 日本vs欧美在线观看视频| 日韩中文字幕欧美一区二区| 欧美在线一区亚洲| 香蕉久久夜色| 午夜福利在线免费观看网站| 最近最新中文字幕大全电影3 | 精品少妇黑人巨大在线播放| 久久久精品94久久精品| av在线播放免费不卡| 国产麻豆69| 亚洲欧美激情在线| 久久天躁狠狠躁夜夜2o2o| 香蕉久久夜色| 亚洲人成77777在线视频| av天堂久久9| 亚洲欧美色中文字幕在线| 精品国产乱码久久久久久小说| 精品少妇一区二区三区视频日本电影| av网站免费在线观看视频| 99九九在线精品视频| 精品少妇黑人巨大在线播放| 国产一区二区 视频在线| 午夜福利欧美成人| bbb黄色大片| 中文亚洲av片在线观看爽 | 亚洲国产中文字幕在线视频| 日本欧美视频一区| 色婷婷av一区二区三区视频| 嫩草影视91久久| 欧美日韩成人在线一区二区| 精品少妇一区二区三区视频日本电影| 99国产精品一区二区蜜桃av | 午夜福利视频在线观看免费| 国产一区二区在线观看av| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av高清一级| 亚洲成人手机| 人人妻人人爽人人添夜夜欢视频| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 日韩大片免费观看网站| 久久久久久久大尺度免费视频| 欧美精品一区二区大全| 色老头精品视频在线观看| 国产av精品麻豆| 亚洲色图综合在线观看| 国产91精品成人一区二区三区 | 国产一区二区三区综合在线观看| 人人澡人人妻人| 在线观看免费视频网站a站| 女人精品久久久久毛片| 日韩欧美国产一区二区入口| 如日韩欧美国产精品一区二区三区| 午夜激情av网站| 99国产综合亚洲精品| 日日摸夜夜添夜夜添小说| 欧美黄色片欧美黄色片| 两个人看的免费小视频| 人妻一区二区av| www.精华液| 国产又爽黄色视频| 777米奇影视久久| 丁香六月天网| 菩萨蛮人人尽说江南好唐韦庄| 国产成人欧美| 国产熟女午夜一区二区三区| 免费在线观看影片大全网站| 制服诱惑二区| 精品一区二区三区视频在线观看免费 | 国产精品久久久久久精品古装| 热re99久久国产66热| 两个人免费观看高清视频| 人人妻人人澡人人看| 亚洲精品国产区一区二| 成人黄色视频免费在线看| 亚洲欧美精品综合一区二区三区| 色精品久久人妻99蜜桃| 国产精品99久久99久久久不卡| 久久精品国产亚洲av高清一级| 丰满少妇做爰视频| 日韩视频一区二区在线观看| 天天躁日日躁夜夜躁夜夜| 黄色视频不卡| 午夜福利视频精品| 91成年电影在线观看| xxxhd国产人妻xxx| 亚洲精品一二三| 婷婷丁香在线五月| 免费少妇av软件| 亚洲成a人片在线一区二区| 天天躁夜夜躁狠狠躁躁| 国产麻豆69| 伊人久久大香线蕉亚洲五| 美女主播在线视频| 天天躁夜夜躁狠狠躁躁| 精品一区二区三区四区五区乱码| 亚洲av成人不卡在线观看播放网| 国产成人av教育| 亚洲精品中文字幕在线视频| 如日韩欧美国产精品一区二区三区| 亚洲成人免费电影在线观看| 欧美成狂野欧美在线观看| 黄网站色视频无遮挡免费观看| 国产黄色免费在线视频| 97人妻天天添夜夜摸| 性色av乱码一区二区三区2| 丰满少妇做爰视频| 黄频高清免费视频| 一边摸一边抽搐一进一出视频| 叶爱在线成人免费视频播放| 我要看黄色一级片免费的| 丝袜美腿诱惑在线| 在线十欧美十亚洲十日本专区| 欧美激情高清一区二区三区| 老司机亚洲免费影院| 亚洲av第一区精品v没综合| 欧美人与性动交α欧美软件| 国产伦人伦偷精品视频| 美女国产高潮福利片在线看| 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 麻豆乱淫一区二区| 成年人黄色毛片网站| 一区二区三区激情视频| 久久久国产成人免费| 成人av一区二区三区在线看| 国产精品秋霞免费鲁丝片| 999久久久国产精品视频| 欧美日韩一级在线毛片| 伊人久久大香线蕉亚洲五| 国产麻豆69| 十八禁高潮呻吟视频| 亚洲人成伊人成综合网2020| 色在线成人网| 免费久久久久久久精品成人欧美视频| 久久久欧美国产精品| 欧美大码av| 1024香蕉在线观看| 怎么达到女性高潮| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 国产精品九九99| 国产1区2区3区精品| av有码第一页| 一本综合久久免费| 亚洲精品中文字幕一二三四区 | 香蕉国产在线看| 亚洲精品中文字幕在线视频| 热99re8久久精品国产| 黑人巨大精品欧美一区二区蜜桃| 在线av久久热| 午夜福利乱码中文字幕| 久久影院123| 黄片小视频在线播放| 午夜免费鲁丝| 女人久久www免费人成看片| 亚洲欧美一区二区三区久久| 在线看a的网站| 97人妻天天添夜夜摸| 视频在线观看一区二区三区| 天堂俺去俺来也www色官网| 久久性视频一级片| 成人av一区二区三区在线看| 操美女的视频在线观看| 天天影视国产精品| av一本久久久久| 深夜精品福利| 精品少妇内射三级| 免费观看人在逋| 热99国产精品久久久久久7| cao死你这个sao货| 久久久精品区二区三区| 99久久国产精品久久久| 亚洲欧美色中文字幕在线| 国产精品免费大片| 高清视频免费观看一区二区| 热re99久久国产66热| 欧美激情高清一区二区三区| 国产精品1区2区在线观看. | 十分钟在线观看高清视频www| 国产xxxxx性猛交| svipshipincom国产片| 日韩欧美三级三区| 亚洲精品av麻豆狂野| 99精国产麻豆久久婷婷| 激情在线观看视频在线高清 | 成人影院久久| 黄色视频不卡| 国产精品久久久久成人av| 亚洲精品中文字幕一二三四区 | bbb黄色大片| 亚洲成人免费av在线播放| 老司机午夜福利在线观看视频 | 可以免费在线观看a视频的电影网站| 在线观看免费高清a一片| 成人亚洲精品一区在线观看| 国产精品熟女久久久久浪| 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 中文字幕人妻熟女乱码| 免费黄频网站在线观看国产| 在线播放国产精品三级| 国产午夜精品久久久久久| 欧美乱码精品一区二区三区| 久久免费观看电影| 老汉色∧v一级毛片| 热99re8久久精品国产| 国产在线观看jvid| 亚洲欧洲精品一区二区精品久久久| 精品人妻熟女毛片av久久网站| 国产1区2区3区精品| 久久人妻熟女aⅴ| 水蜜桃什么品种好| 香蕉国产在线看| 亚洲伊人色综图| av又黄又爽大尺度在线免费看| 少妇被粗大的猛进出69影院| 国产av国产精品国产| 国产有黄有色有爽视频| 一级黄色大片毛片| 日本一区二区免费在线视频| 国产一区有黄有色的免费视频| 成人18禁在线播放| 捣出白浆h1v1| 欧美人与性动交α欧美软件| 十八禁人妻一区二区| 飞空精品影院首页| 老司机午夜福利在线观看视频 | 大片电影免费在线观看免费| 美女高潮喷水抽搐中文字幕| 色播在线永久视频| 久久久久久久精品吃奶| 999精品在线视频| 免费不卡黄色视频| 久久亚洲精品不卡| 日日摸夜夜添夜夜添小说| 男女之事视频高清在线观看| 一区二区三区激情视频| 99久久精品国产亚洲精品| 久久久精品免费免费高清| 99精品在免费线老司机午夜| 男女免费视频国产| 一区二区三区激情视频| 亚洲成a人片在线一区二区| 国产av一区二区精品久久| 窝窝影院91人妻| 黄频高清免费视频| 在线观看www视频免费| 视频区图区小说| av在线播放免费不卡| 精品国产乱子伦一区二区三区| 久久人妻福利社区极品人妻图片| 一边摸一边做爽爽视频免费| aaaaa片日本免费| 考比视频在线观看| 亚洲成人手机| 久久精品91无色码中文字幕| 国产欧美日韩一区二区三区在线| 高清av免费在线| 日本黄色日本黄色录像| 欧美性长视频在线观看| 如日韩欧美国产精品一区二区三区| 涩涩av久久男人的天堂| 在线播放国产精品三级| 一边摸一边抽搐一进一出视频| 欧美精品一区二区免费开放| 无遮挡黄片免费观看| 最新美女视频免费是黄的| 成人亚洲精品一区在线观看| 亚洲专区国产一区二区| 精品一区二区三区四区五区乱码| 欧美变态另类bdsm刘玥| 男男h啪啪无遮挡| 一级片免费观看大全| 狠狠精品人妻久久久久久综合| 久久久国产一区二区| 成年人免费黄色播放视频| 不卡av一区二区三区| 日本av免费视频播放| 天天躁狠狠躁夜夜躁狠狠躁| 我的亚洲天堂| 在线看a的网站| 免费在线观看影片大全网站| 大片电影免费在线观看免费| 亚洲精品国产精品久久久不卡| a级毛片黄视频| 免费观看a级毛片全部| 精品一区二区三卡| 久久 成人 亚洲| av超薄肉色丝袜交足视频| 久久久久久久大尺度免费视频| 欧美另类亚洲清纯唯美| 亚洲精品粉嫩美女一区| 欧美日韩av久久| 99久久国产精品久久久| 热99久久久久精品小说推荐| 欧美日韩国产mv在线观看视频| 亚洲人成电影观看| tube8黄色片| xxxhd国产人妻xxx| 欧美在线一区亚洲| 韩国精品一区二区三区| svipshipincom国产片| 亚洲av电影在线进入| 欧美成人午夜精品| 欧美黄色淫秽网站| 国产欧美日韩精品亚洲av| 老司机亚洲免费影院| 中文字幕人妻熟女乱码| 久久精品国产亚洲av香蕉五月 | 午夜福利视频在线观看免费| 淫妇啪啪啪对白视频| 满18在线观看网站| 91字幕亚洲| 欧美精品av麻豆av| 久久中文字幕一级| 在线观看www视频免费| 欧美一级毛片孕妇| 国产精品.久久久| 侵犯人妻中文字幕一二三四区| 国产麻豆69| 国产精品免费一区二区三区在线 | 国产亚洲av高清不卡| 欧美黄色淫秽网站| 亚洲熟妇熟女久久| 一二三四在线观看免费中文在| 色婷婷av一区二区三区视频| 深夜精品福利| 亚洲午夜精品一区,二区,三区| 亚洲人成电影免费在线| 久久毛片免费看一区二区三区| 国产亚洲精品一区二区www | 国产免费视频播放在线视频| 国产黄色免费在线视频| 欧美精品高潮呻吟av久久| 中文字幕av电影在线播放| 久久久久久久久久久久大奶| 国产男靠女视频免费网站| 日本av手机在线免费观看| 欧美精品人与动牲交sv欧美| 99国产极品粉嫩在线观看| 精品一区二区三区av网在线观看 | 热re99久久国产66热| 高清毛片免费观看视频网站 | 国产精品99久久99久久久不卡| 免费不卡黄色视频| 久久久久精品人妻al黑| 一二三四在线观看免费中文在| 成人国语在线视频| 欧美激情极品国产一区二区三区| 亚洲精品成人av观看孕妇| 精品久久蜜臀av无| 国产精品.久久久| 丰满饥渴人妻一区二区三| 国产亚洲欧美精品永久| 久久精品人人爽人人爽视色| 动漫黄色视频在线观看| 免费久久久久久久精品成人欧美视频| 黑丝袜美女国产一区| 老汉色∧v一级毛片| www.熟女人妻精品国产| 免费女性裸体啪啪无遮挡网站| 亚洲专区中文字幕在线| 中文字幕最新亚洲高清|