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

    基于互信息的深證股票復(fù)雜網(wǎng)絡(luò)拓?fù)湫再|(zhì)分析

    2016-11-28 02:08:26胡振華覃子龍楊燕
    統(tǒng)計(jì)與決策 2016年20期
    關(guān)鍵詞:概率密度函數(shù)互信息信息熵

    胡振華,覃子龍,楊燕

    (中南大學(xué)商學(xué)院,長(zhǎng)沙410083)

    基于互信息的深證股票復(fù)雜網(wǎng)絡(luò)拓?fù)湫再|(zhì)分析

    胡振華,覃子龍,楊燕

    (中南大學(xué)商學(xué)院,長(zhǎng)沙410083)

    文章針對(duì)以往在非線性系統(tǒng)中不當(dāng)使用線性系統(tǒng)統(tǒng)計(jì)指標(biāo)——Pearson相關(guān)系數(shù)來(lái)構(gòu)建復(fù)雜網(wǎng)絡(luò)的方法,以深證股票為研究對(duì)象,使用在非線性系統(tǒng)中衡量隨機(jī)變量之間相關(guān)性的指標(biāo)——互信息來(lái)構(gòu)建復(fù)雜網(wǎng)絡(luò)并分析了該網(wǎng)絡(luò)的拓?fù)湫再|(zhì),認(rèn)為深證股票復(fù)雜網(wǎng)絡(luò)具有小世界特征和無(wú)標(biāo)度特征,深證股票緊密程度高,市場(chǎng)信息傳遞效率也較高,社團(tuán)結(jié)構(gòu)呈現(xiàn)第二產(chǎn)業(yè)集聚、第三產(chǎn)業(yè)分散的特點(diǎn),股票收益率的波動(dòng)與產(chǎn)業(yè)鏈具有更大的相關(guān)性,網(wǎng)絡(luò)對(duì)隨機(jī)故障具有較高的魯棒性,而對(duì)蓄意攻擊具有高度的脆弱性。

    非線性系統(tǒng);信息熵;互信息;復(fù)雜網(wǎng)絡(luò);拓?fù)湫再|(zhì)

    0 引言

    國(guó)內(nèi)外大量的研究表明,金融市場(chǎng)具有參與要素眾多、內(nèi)部結(jié)構(gòu)復(fù)雜等特點(diǎn),在本質(zhì)上是一個(gè)典型的非線性動(dòng)力學(xué)系統(tǒng),如收益率“尖峰厚尾”的分布特征、收益率波動(dòng)的集聚現(xiàn)象等[1,2]。非線性理論對(duì)此提供了新的研究范式,其中網(wǎng)絡(luò)科學(xué)把金融市場(chǎng)視為一個(gè)系統(tǒng),認(rèn)為系統(tǒng)內(nèi)的個(gè)體是非勻質(zhì)的,以異質(zhì)的個(gè)體為節(jié)點(diǎn),以個(gè)體之間存在的一定關(guān)系為邊,構(gòu)建出相應(yīng)的復(fù)雜網(wǎng)絡(luò)模型[3,4]。構(gòu)建股票復(fù)雜網(wǎng)絡(luò)模型主要采用兩種方法:閾值法[5]和拓?fù)浼s束法[6],這兩種方法均以Pearson相關(guān)系數(shù)為基礎(chǔ)。在計(jì)算Pearson相關(guān)系數(shù)時(shí),除了要求兩個(gè)變量符合聯(lián)合正態(tài)分布外,更要求兩個(gè)變量之間呈直線相關(guān)關(guān)系。而股票收益率時(shí)間序列數(shù)據(jù)大多呈現(xiàn)出“尖峰厚尾”的特征,并不服從正態(tài)分布的假設(shè);并且在上述方法中,以線性系統(tǒng)中的統(tǒng)計(jì)指標(biāo)為基礎(chǔ)來(lái)建立非線性系統(tǒng)模型,其適用性也值得商榷,可能產(chǎn)生潛在的估計(jì)錯(cuò)誤風(fēng)險(xiǎn)。事實(shí)上,起源于信息論的互信息概念,可用于衡量非線性系統(tǒng)中兩個(gè)或多個(gè)變量之間的相關(guān)關(guān)系,目前已經(jīng)被廣泛應(yīng)用于生物基因、水文、經(jīng)濟(jì)等諸多領(lǐng)域。本文以信息熵和互信息為基礎(chǔ),以深證100指數(shù)為研究對(duì)象,闡明股票復(fù)雜網(wǎng)絡(luò)模型的構(gòu)建過(guò)程,并分析該復(fù)雜網(wǎng)絡(luò)的拓?fù)湫再|(zhì)。

    1 信息熵和互信息

    1.1 信息熵

    信息熵的概念隸屬于信息論領(lǐng)域,由熱力學(xué)熵的概念泛化而來(lái),Shannon首先定義了信息熵作為某件事情不確定性的度量。對(duì)于離散型隨機(jī)變量X,其信息熵H(X)的定義式為:

    其中,p(x)為x每種取值的概率,對(duì)于底數(shù)b的不同取值意味著信息熵H(X)有不同的量綱,以下不再贅述。對(duì)于連續(xù)型隨機(jī)變量X,其信息熵H(X)的定義式為:

    二維推廣,對(duì)于離散型、連續(xù)型的隨機(jī)變量X、y,分別有:

    其中,p(x,y)為X、y的聯(lián)合概率,f(x,y)為X、y的聯(lián)合概率密度函數(shù)。

    1.2 互信息

    互信息用于度量?jī)蓚€(gè)或多個(gè)變量所共享的信息量,互信息越大則變量間的相關(guān)性越強(qiáng),反之則反然。對(duì)于離散型、連續(xù)型的二維隨機(jī)變量X、y,其定義式分別為:

    其中,pX(x)、py(y)分別為X、y的邊緣概率, fX(x,y)、fy(x,y)分別為X、y的邊緣概率密度函數(shù)。

    根據(jù)上述公式,易得:

    1.3 信息熵和互信息的估計(jì)

    對(duì)信息熵和互信息的估計(jì)主要采用兩種方法:參數(shù)方法和非參數(shù)方法。使用參數(shù)方法時(shí),首先需假定概率密度函數(shù)的形式,再根據(jù)樣本數(shù)據(jù)估計(jì)出其中的各項(xiàng)參數(shù),進(jìn)而求得信息熵和互信息。但在實(shí)際應(yīng)用中,或者概率密度函數(shù)通常難以得知,或者參數(shù)方法無(wú)法處理非線性不規(guī)則的聯(lián)合分布,使得參數(shù)方法的應(yīng)用受到極大限制,而發(fā)展出的非參數(shù)方法得到更為廣泛的應(yīng)用。

    常見(jiàn)的非參數(shù)方法有直方圖法[7]和核函數(shù)法[8]。直方圖法是把區(qū)間或者區(qū)域劃分成若干個(gè)bin,進(jìn)而估計(jì)概率密度函數(shù),求得信息熵和互信息的方法。該方法簡(jiǎn)單易行,但精度較低,只適用于低維數(shù)據(jù),計(jì)算誤差主要來(lái)源于bin的數(shù)量的選擇,換而言之,與bin的區(qū)間寬度h有關(guān)。核函數(shù)法是在每個(gè)樣本點(diǎn)上疊置一個(gè)核函數(shù),將所有核函數(shù)疊加而作為概率密度函數(shù)的一個(gè)估計(jì),求得信息熵和互信息的方法。該方法計(jì)算復(fù)雜度高,精度也高,但只適用于小樣本的情形,當(dāng)樣本量較大時(shí),核參數(shù)的估計(jì)將會(huì)變得十分困難。值得一提的是,為了避免聯(lián)合概率密度函數(shù)的估計(jì),可以較為方便的通過(guò)Copula熵來(lái)估計(jì)互信息[9,10]。

    2 基于互信息的深證股票復(fù)雜網(wǎng)絡(luò)模型構(gòu)建

    鑒于時(shí)間序列的樣本數(shù)據(jù)規(guī)模通常較大,本文舍棄核函數(shù)法估計(jì)信息熵和互信息;在將信息熵和互信息轉(zhuǎn)換為歐氏空間距離時(shí),為了保持分子、分母計(jì)算口徑的一致性,本文舍棄了Copula熵法估計(jì)互信息,統(tǒng)一選用直方圖法估計(jì)信息熵和互信息,進(jìn)而構(gòu)建相應(yīng)的復(fù)雜網(wǎng)絡(luò)。本文所有的數(shù)據(jù)處理和分析通過(guò)Matlab 7.11.0、Mathematica 9.0和Ucinet6實(shí)現(xiàn)。

    2.1 數(shù)據(jù)來(lái)源與預(yù)處理

    本文選取深證100指數(shù)作為分析對(duì)象,選取指標(biāo)為每日收盤(pán)價(jià)。其中,深證100指數(shù)來(lái)自深圳證券交易所網(wǎng)站(http://www.szse.cn)發(fā)布的成分股,數(shù)據(jù)下載時(shí)間為2014年6月23日。股價(jià)數(shù)據(jù)的時(shí)間跨度為2011年1月4日至2014年3月31日,共783個(gè)觀測(cè)日,74,895個(gè)觀測(cè)數(shù)據(jù),數(shù)據(jù)來(lái)源為銳思數(shù)據(jù)庫(kù)(http://www.resset.cn)。對(duì)因上市時(shí)間較晚、停牌時(shí)間較長(zhǎng)以及其他原因使得數(shù)據(jù)缺失較多的10支成分股數(shù)據(jù)予以剔除,對(duì)剩下90支股票的缺失數(shù)據(jù)采用線性插值法予以補(bǔ)齊。

    2.2 構(gòu)建股票復(fù)雜網(wǎng)絡(luò)

    (1)計(jì)算股票的對(duì)數(shù)收益率

    其中,Rx(t)表示股票X第t期的對(duì)數(shù)收益率,px(t)表示股票X第t期的收盤(pán)價(jià)。

    (2)采用等距法計(jì)算兩兩股票X、y的信息熵和互信息

    ①計(jì)算bin的最優(yōu)寬度h*[11]

    對(duì)于股票X,有:

    其中,n為樣本規(guī)模,σx為Rx的樣本標(biāo)準(zhǔn)差,同理可得

    ②估計(jì)概率密度函數(shù)

    構(gòu)造Rx的M個(gè)bin區(qū)間設(shè)Rx落在區(qū)間上的樣本點(diǎn)數(shù)量為ui,則有Rx的概率密度函數(shù):

    其中,N為Ry的bin區(qū)間數(shù)量,vi為Ry落在區(qū)間上的樣本點(diǎn)數(shù)量,wij為Rx落在區(qū)間、Ry落在區(qū)間[ymin+(j-1)上的樣本點(diǎn)數(shù)量。

    ③計(jì)算信息熵和互信息

    依據(jù)式(2)、式(4)分別計(jì)算H(X)、H(y)和H(X,y),依據(jù)式(7)計(jì)算I(X,y)。

    (3)計(jì)算歐氏距離

    由于互信息不完全滿足三角不等式,還需將其轉(zhuǎn)換為歐氏空間距離d(X,y),計(jì)算公式為[12]:

    (4)構(gòu)建復(fù)雜網(wǎng)絡(luò)

    由于閾值法具有一定的主觀性,并且閾值過(guò)低起不到良好的過(guò)濾冗余信息的效果;最小生成樹(shù)的生成方法又存在信息過(guò)濾過(guò)度的問(wèn)題;而平面最大過(guò)濾圖不僅保留了最小生成樹(shù)所包含的網(wǎng)絡(luò)層級(jí)結(jié)構(gòu),還包含了網(wǎng)絡(luò)內(nèi)部結(jié)構(gòu)的其他大量信息[13],因此本文使用平面最大過(guò)濾圖來(lái)生成相應(yīng)的復(fù)雜網(wǎng)絡(luò)。

    以股票為節(jié)點(diǎn),以歐氏空間距離為邊,設(shè)復(fù)雜網(wǎng)絡(luò)共包含V個(gè)節(jié)點(diǎn),以平面最大過(guò)濾圖的方法來(lái)構(gòu)建相應(yīng)的復(fù)雜網(wǎng)絡(luò)的步驟為[13]:初始網(wǎng)絡(luò)不含任何連邊,將d(X,y)按降序排列,逐次將d(X,y)所對(duì)應(yīng)的邊加入網(wǎng)絡(luò),當(dāng)且僅當(dāng)生成的網(wǎng)絡(luò)仍是可平面圖時(shí)加入該邊,否則取消這條邊的加入,直至網(wǎng)絡(luò)包含3(V-2)條邊,生成的深證股票平面最大過(guò)濾圖PMFG如圖1所示。

    圖1 深證股票復(fù)雜網(wǎng)絡(luò)PMFG

    3 深證股票復(fù)雜網(wǎng)絡(luò)的拓?fù)湫再|(zhì)分析

    3.1 小世界特性

    較多復(fù)雜網(wǎng)絡(luò)都具有與相同規(guī)模隨機(jī)網(wǎng)絡(luò)相當(dāng)?shù)妮^小的平均路徑長(zhǎng)度,同時(shí)又具有與相同規(guī)模規(guī)則網(wǎng)絡(luò)相當(dāng)?shù)妮^大的平均聚類系數(shù),復(fù)雜網(wǎng)絡(luò)所表現(xiàn)出的這種統(tǒng)計(jì)特性就是小世界特性。根據(jù)平面最大過(guò)濾圖的生成規(guī)則, PMFG的平均度必然為以此為標(biāo)準(zhǔn),生成相同規(guī)模的ER隨機(jī)圖和最近鄰耦合網(wǎng)絡(luò),PMFG與它們的拓?fù)湫再|(zhì)對(duì)比見(jiàn)表1。可以看到, PMFG既具有與ER隨機(jī)圖相近的平均路徑長(zhǎng)度,又具有與最近鄰耦合網(wǎng)絡(luò)相近的平均聚類系數(shù),小世界特征十分明顯。平均路徑長(zhǎng)度刻畫(huà)了某個(gè)節(jié)點(diǎn)平均需要經(jīng)過(guò)多少個(gè)節(jié)點(diǎn)即可到達(dá)另一個(gè)節(jié)點(diǎn),PMFG較小的平均路徑長(zhǎng)度說(shuō)明深證股票市場(chǎng)的信息傳遞效率較高;平均聚類系數(shù)刻畫(huà)的是某個(gè)節(jié)點(diǎn)的任意兩個(gè)鄰居節(jié)點(diǎn)也互為鄰居的可能性,PMFG較大的平均聚類系數(shù)說(shuō)明深證股票的緊密程度較高。

    表1 復(fù)雜網(wǎng)絡(luò)PMFG和隨機(jī)網(wǎng)絡(luò)、規(guī)則網(wǎng)絡(luò)的拓?fù)湫再|(zhì)對(duì)比

    3.2 無(wú)標(biāo)度特性

    不同于隨機(jī)變量常見(jiàn)的鐘形分布,大部分取值位于均值/峰值的附近,或者說(shuō)落在以特征標(biāo)度為中心的一個(gè)相對(duì)較小的區(qū)間里?,F(xiàn)實(shí)中的很多復(fù)雜網(wǎng)絡(luò)的大部分節(jié)點(diǎn)度都比較小,但也存在少數(shù)較大的節(jié)點(diǎn)度,即服從“長(zhǎng)尾”分布,不存在單一的特征標(biāo)度,這一統(tǒng)計(jì)特性就是復(fù)雜網(wǎng)絡(luò)的無(wú)標(biāo)度特性。如圖2(a)所示,復(fù)雜網(wǎng)絡(luò)PMFG的度分布很不均勻,大多數(shù)節(jié)點(diǎn)的度數(shù)都很小,一半以上股票的度數(shù)為2;但也有少量度數(shù)較大的Hub節(jié)點(diǎn),節(jié)點(diǎn)度最高的股票是樂(lè)視網(wǎng),度數(shù)達(dá)到66,也就是說(shuō),該股票與網(wǎng)絡(luò)中近75%的股票有連接,存在相關(guān)關(guān)系。如圖2(b)所示,在雙對(duì)數(shù)直角坐標(biāo)系下,其度分布大致在一條斜率為-1.1708的直線附近,亦即PMFG的度分布服從參數(shù)為γ=1.1708的冪律分布。

    圖2 復(fù)雜網(wǎng)絡(luò)PMFG的度分布

    3.3 社團(tuán)結(jié)構(gòu)

    很多實(shí)際網(wǎng)絡(luò)都具有社團(tuán)結(jié)構(gòu),即整個(gè)復(fù)雜網(wǎng)絡(luò)由若干個(gè)社團(tuán)構(gòu)成,每個(gè)社團(tuán)內(nèi)部的節(jié)點(diǎn)之間連接相對(duì)緊密,社團(tuán)之間的節(jié)點(diǎn)之間連接相對(duì)稀疏。模塊度Q是常用來(lái)衡量社團(tuán)劃分質(zhì)量的指標(biāo)[14],其上限為1,通常認(rèn)為只要Q值大于0.3就說(shuō)明網(wǎng)絡(luò)存在明顯的社團(tuán)結(jié)構(gòu)[15]。Ucinet6在劃分社團(tuán)結(jié)構(gòu)時(shí)采用的是Girvan-Newman算法,該算法是以邊介數(shù)為社團(tuán)間瓶頸的社團(tuán)發(fā)現(xiàn)算法[16]。當(dāng)PMFG劃分為8個(gè)社團(tuán)時(shí),模塊度Q取得最大值0.427,相應(yīng)的社團(tuán)結(jié)構(gòu)及其行業(yè)分布如圖3、表2所示,其中行業(yè)分類依據(jù)為中國(guó)證券監(jiān)督管理委員會(huì)2012年頒布的《上市公司行業(yè)分類指引》。

    圖3 復(fù)雜網(wǎng)絡(luò)PMFG的社團(tuán)結(jié)構(gòu)

    表2 復(fù)雜網(wǎng)絡(luò)PMFG社團(tuán)的行業(yè)分布

    從表2中可以看到,除社團(tuán)1和社團(tuán)2的節(jié)點(diǎn)成員較少外,其余社團(tuán)的行業(yè)分布呈現(xiàn)如下兩個(gè)特點(diǎn):①第二產(chǎn)業(yè)呈現(xiàn)集聚的特點(diǎn),其中制造業(yè)(C)基本占了各個(gè)社團(tuán)節(jié)點(diǎn)成員數(shù)量的一半,甚至絕大多數(shù);②第三產(chǎn)業(yè)呈現(xiàn)分散的特點(diǎn),尤其是金融業(yè)(J)和房地產(chǎn)業(yè)(K)幾乎出現(xiàn)在各個(gè)社團(tuán)中。以上社團(tuán)結(jié)構(gòu)特點(diǎn)說(shuō)明:其一,第二產(chǎn)業(yè)在國(guó)民經(jīng)濟(jì)中仍占據(jù)著舉足輕重的地位;其二,第三產(chǎn)業(yè)在廣度、深度上都取得了長(zhǎng)足的發(fā)展,在國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展中的地位日益突出;其三,相較于行業(yè),股票收益率的波動(dòng)應(yīng)與產(chǎn)業(yè)鏈具有更大的相關(guān)性。股票網(wǎng)絡(luò)的社團(tuán)結(jié)構(gòu)劃分對(duì)于組合投資具有重要意義,投資者應(yīng)當(dāng)盡量將資金投資于不同社團(tuán)的股票以分散投資風(fēng)險(xiǎn),增強(qiáng)投資的穩(wěn)定性,僅僅依賴股票的行業(yè)歸屬等直觀信息進(jìn)行決策是不可靠的。

    3.4 魯棒性與脆弱性

    當(dāng)復(fù)雜網(wǎng)絡(luò)受到攻擊時(shí),網(wǎng)絡(luò)的連通程度就會(huì)受到破壞,傳輸能力也會(huì)受到影響[17]。為了分析復(fù)雜網(wǎng)絡(luò)PMFG的魯棒性與脆弱性,考慮兩種節(jié)點(diǎn)去除策略[18]:①蓄意攻擊策略,即從去除網(wǎng)絡(luò)中度最高的節(jié)點(diǎn)開(kāi)始,依次去除部分度最高的節(jié)點(diǎn);②隨機(jī)故障策略,即完全隨機(jī)地去除網(wǎng)絡(luò)中的部分節(jié)點(diǎn)。f為去除的節(jié)點(diǎn)數(shù)占復(fù)雜網(wǎng)絡(luò)PMFG總節(jié)點(diǎn)數(shù)的比例,S為最大連通子圖的相對(duì)大小,L為平均路徑長(zhǎng)度,各指標(biāo)的變化趨勢(shì)如圖4所示,其中對(duì)于隨機(jī)故障策略取50次獨(dú)立實(shí)驗(yàn)的平均值??梢钥吹?由于復(fù)雜網(wǎng)絡(luò)PMFG具有無(wú)標(biāo)度特征,對(duì)隨機(jī)故障具有較高的魯棒性,最大連通子圖的相對(duì)大小下降得較為緩慢,平均路徑長(zhǎng)度的變化則并不顯著;但是,對(duì)于蓄意攻擊具有高度的脆弱性,去除網(wǎng)絡(luò)中少量度最高的節(jié)點(diǎn)就可破壞網(wǎng)絡(luò)的連通性,在去除40%左右的節(jié)點(diǎn)時(shí),S就趨緊于0了。由于股市對(duì)經(jīng)濟(jì)的波動(dòng)是異常敏感的,因此搜集關(guān)鍵股票節(jié)點(diǎn)的信息,掌握經(jīng)濟(jì)沖擊股市的方式和范圍,對(duì)維持股市的穩(wěn)定、防止遭受惡意破壞、保持國(guó)家經(jīng)濟(jì)的穩(wěn)健發(fā)展都有著十分重要的意義。

    圖4 復(fù)雜網(wǎng)絡(luò)PMFG的魯棒性與脆弱性

    4 結(jié)論

    復(fù)雜網(wǎng)絡(luò)的建模方法對(duì)于非線性系統(tǒng)的分析提供了一種新的研究方法,但以往以線性系統(tǒng)中的Pearson相關(guān)系數(shù)為基礎(chǔ)構(gòu)建復(fù)雜網(wǎng)絡(luò)的方法與非線性系統(tǒng)的建模初衷是不相吻合的。本文以深證股票為研究對(duì)象,以非線性系統(tǒng)中衡量?jī)蓚€(gè)變量相關(guān)關(guān)系的互信息指標(biāo)為基礎(chǔ),構(gòu)建了復(fù)雜網(wǎng)絡(luò)PMFG。

    通過(guò)對(duì)其拓?fù)湫再|(zhì)進(jìn)行分析,認(rèn)為復(fù)雜網(wǎng)絡(luò)PMFG具有小世界特征,說(shuō)明深證股票市場(chǎng)的信息傳遞效率較高,股票的緊密程度也較高;網(wǎng)絡(luò)中大多數(shù)的節(jié)點(diǎn)度數(shù)很小,但也存在少量Hub節(jié)點(diǎn),具有無(wú)標(biāo)度特征;網(wǎng)絡(luò)具有明顯的社團(tuán)結(jié)構(gòu),呈現(xiàn)第二產(chǎn)業(yè)集聚、第三產(chǎn)業(yè)分散的特點(diǎn),股票收益率的波動(dòng)與產(chǎn)業(yè)鏈具有更大的相關(guān)性;網(wǎng)絡(luò)對(duì)隨機(jī)故障具有較高的魯棒性,而對(duì)蓄意攻擊具有高度的脆弱性。

    在以往對(duì)以Pearson相關(guān)系數(shù)為基礎(chǔ)構(gòu)建的復(fù)雜網(wǎng)絡(luò)的研究中,復(fù)雜網(wǎng)絡(luò)同樣具有明顯的社團(tuán)結(jié)構(gòu),但顯示的是股票收益率的波動(dòng)與行業(yè)具有更大的相關(guān)性[19]。本文的結(jié)論與此有所區(qū)別,該差異可能來(lái)自于以下幾個(gè)原因:研究方法的差異導(dǎo)致結(jié)論的差異、我國(guó)股票市場(chǎng)的有效程度不如美國(guó)等發(fā)達(dá)國(guó)家、我國(guó)的產(chǎn)業(yè)結(jié)構(gòu)仍處在不斷調(diào)整的過(guò)程中。因此,就我國(guó)股票網(wǎng)絡(luò)的社團(tuán)結(jié)構(gòu)及其內(nèi)部構(gòu)成還有待進(jìn)一步的分析。

    [1]黃登仕.金融市場(chǎng)的標(biāo)度理論[J].管理科學(xué)學(xué)報(bào),2000,3(2).

    [2]Patton A J.Copula一Based Models for Financial Time Series[A].Andersen T G,Davis R A,KreiB J一P,et al.Handbook of Financial Time Series[M].Heidelberg:Springer,2009.

    [3]馬源源,莊新田,李凌軒.股市中危機(jī)傳播的SIR模型及其仿真[J].管理科學(xué)學(xué)報(bào),2013,16(7).

    [4]張鼎,莊新田,盧文娟,曲紅濤.基于復(fù)雜網(wǎng)絡(luò)理論的期貨指數(shù)網(wǎng)絡(luò)的研究[J].系統(tǒng)管理學(xué)報(bào),2014,23(1).

    [5]張鼎,莊新田,盧文娟.我國(guó)滬深股市同配性和網(wǎng)絡(luò)導(dǎo)航現(xiàn)象研究——基于復(fù)雜網(wǎng)絡(luò)理論[J].技術(shù)經(jīng)濟(jì),2013,32(7).

    [6]黃瑋強(qiáng),莊新田,姚爽.中國(guó)股票關(guān)聯(lián)網(wǎng)絡(luò)拓?fù)湫再|(zhì)與聚類結(jié)構(gòu)分析[J].管理科學(xué),2008,21(3).

    [7]Meyer PE,Lafitte F,BontempiG.Minet:A R/Bioconductor Package for Inferring Large TranscriptionalNetworksUsingMutual Information [J].BMCBioinformatics,2008,9.

    [8]Moon y I,Rajagopalan B,Lall U.Estimation ofMutual Information Using KernelDensity Estimators[J].PhysicalReview E,1995,52(3).

    [9]Zeng X,Durrani TS.Estimation ofMutual Information Using Copula Density Function[J].Electronics Letters,2011,47(8).

    [10]韓敏,劉曉欣.基于Copula熵的互信息估計(jì)方法[J].控制理論與應(yīng)用,2013,30(7).

    [11]龔偉.基于信息熵和互信息的流域水文模型不確定性分析[D].北京:清華大學(xué)碩士論文,2012.

    [12]Kraskov A,Stogbauer H,Andrzejak RG,et al.Hierarchical Clustering Using Mutual Information[J].Europhysics Letters,2005,70 (2).

    [13]Tumminello M,Aste T,DiM T,et al.A Tool for Filtering Information in Complex Systems[J].Proc.Natl.Acad.Sci.USA,2005,102 (30).

    [14]Newman M E J,Girvan M.Finding and Evaluating Community Structure in Networks[J].PhysicalReview E,2004,69(2).

    [15]Newman M E J.Detecting Community Structure in Networks[J].European Physical Journal B,2004,38(2).

    [16]Girvan M,Newman M E J.Community Structure in Socialand BiologicalNetworks[J].Proc.Natl.Acad.Sci.USA,2002,99(12).

    [17]KurantM,Thiran P,Hagmann P.Error and Attack Tolerance of Layered Complex Networks[J].PhysicalReview E,2007,76(2).

    [18]Albert R,Jeong H,BarabásiA L.Error and Attack Tolerance ofComplex Networks[J].Nature,2000,406.

    [19]Jiang X F,Chen TT,Zheng B.Structure of Local Interactions in Complex FinancialDynamics[J].Scientific Reports,2014,(4).

    (責(zé)任編輯/亦民)

    F830.91

    A

    1002-6487(2016)20-0160-04

    胡振華(1962—),男,湖南邵陽(yáng)人,教授,研究方向:金融計(jì)量、貨幣政策。覃子龍(1982—),男,湖南張家界人,博士研究生,研究方向:復(fù)雜網(wǎng)絡(luò)、金融計(jì)量。

    猜你喜歡
    概率密度函數(shù)互信息信息熵
    基于信息熵可信度的測(cè)試點(diǎn)選擇方法研究
    冪分布的有效估計(jì)*
    已知f(x)如何求F(x)
    基于信息熵的實(shí)驗(yàn)教學(xué)量化研究
    一種基于信息熵的雷達(dá)動(dòng)態(tài)自適應(yīng)選擇跟蹤方法
    基于互信息的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)
    聯(lián)合互信息水下目標(biāo)特征選擇算法
    基于信息熵的IITFN多屬性決策方法
    改進(jìn)的互信息最小化非線性盲源分離算法
    基于增量式互信息的圖像快速匹配方法
    每晚都被弄得嗷嗷叫到高潮| 欧美性长视频在线观看| 亚洲精品国产一区二区精华液| 中国美女看黄片| 国内毛片毛片毛片毛片毛片| 国产精品久久久人人做人人爽| 嫩草影视91久久| 三级男女做爰猛烈吃奶摸视频| 脱女人内裤的视频| 色尼玛亚洲综合影院| 99精品久久久久人妻精品| 90打野战视频偷拍视频| www日本黄色视频网| 少妇粗大呻吟视频| av免费在线观看网站| 在线观看舔阴道视频| 黑人欧美特级aaaaaa片| 色av中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久久久人妻精品电影| 久久久国产精品麻豆| 丰满人妻一区二区三区视频av | 亚洲精品在线美女| 日韩精品中文字幕看吧| 日日爽夜夜爽网站| 欧美成人免费av一区二区三区| 亚洲人与动物交配视频| 精品久久久久久久末码| 亚洲精品中文字幕一二三四区| 少妇裸体淫交视频免费看高清 | 亚洲aⅴ乱码一区二区在线播放 | 丰满人妻熟妇乱又伦精品不卡| 在线看三级毛片| 亚洲欧美日韩高清专用| 色精品久久人妻99蜜桃| 国产精品爽爽va在线观看网站| 亚洲成人精品中文字幕电影| 欧美+亚洲+日韩+国产| 亚洲av成人一区二区三| 国产亚洲欧美在线一区二区| 在线免费观看的www视频| 男人舔女人下体高潮全视频| 日本五十路高清| ponron亚洲| 精品少妇一区二区三区视频日本电影| 午夜日韩欧美国产| 两人在一起打扑克的视频| 人成视频在线观看免费观看| 丝袜美腿诱惑在线| 亚洲人成77777在线视频| 老汉色av国产亚洲站长工具| 国产精品,欧美在线| 久久婷婷人人爽人人干人人爱| 97人妻精品一区二区三区麻豆| 国产精品影院久久| 成人av一区二区三区在线看| 欧美不卡视频在线免费观看 | 久久天堂一区二区三区四区| 婷婷六月久久综合丁香| 一级a爱片免费观看的视频| 欧美av亚洲av综合av国产av| 婷婷精品国产亚洲av在线| 91字幕亚洲| 久久久久亚洲av毛片大全| 日本黄大片高清| 国产精品乱码一区二三区的特点| 日日爽夜夜爽网站| 日日摸夜夜添夜夜添小说| 精品国内亚洲2022精品成人| 高潮久久久久久久久久久不卡| 脱女人内裤的视频| 国产激情偷乱视频一区二区| 真人做人爱边吃奶动态| 国产激情欧美一区二区| 久久久精品大字幕| 国产av麻豆久久久久久久| 窝窝影院91人妻| 黄色毛片三级朝国网站| 亚洲中文日韩欧美视频| 在线观看舔阴道视频| 一本综合久久免费| 国产免费男女视频| 欧美绝顶高潮抽搐喷水| 少妇熟女aⅴ在线视频| 国产精品电影一区二区三区| 成人欧美大片| 国产精品野战在线观看| 亚洲国产欧美网| 久久精品国产综合久久久| 国产亚洲精品av在线| 亚洲国产日韩欧美精品在线观看 | 国产精品一及| 美女黄网站色视频| 99久久精品热视频| 波多野结衣高清作品| 国产av一区二区精品久久| 亚洲片人在线观看| 欧美在线黄色| 两个人免费观看高清视频| 岛国视频午夜一区免费看| 天堂动漫精品| 九九热线精品视视频播放| 90打野战视频偷拍视频| 日本免费a在线| 亚洲欧美一区二区三区黑人| 69av精品久久久久久| 国模一区二区三区四区视频 | 亚洲精品粉嫩美女一区| 国产亚洲欧美在线一区二区| videosex国产| 俄罗斯特黄特色一大片| 免费看美女性在线毛片视频| 一夜夜www| 俺也久久电影网| 99riav亚洲国产免费| 毛片女人毛片| 午夜老司机福利片| 免费看日本二区| 亚洲国产中文字幕在线视频| 法律面前人人平等表现在哪些方面| 成人18禁高潮啪啪吃奶动态图| 欧美最黄视频在线播放免费| 欧美乱码精品一区二区三区| 欧美不卡视频在线免费观看 | 免费电影在线观看免费观看| 美女免费视频网站| 啦啦啦免费观看视频1| 国产主播在线观看一区二区| 午夜免费观看网址| 色综合站精品国产| а√天堂www在线а√下载| 久久久水蜜桃国产精品网| 国产私拍福利视频在线观看| 长腿黑丝高跟| 国产精品久久久人人做人人爽| 亚洲成人久久性| 两个人的视频大全免费| 好男人电影高清在线观看| 久久婷婷人人爽人人干人人爱| 国产精品免费一区二区三区在线| 国产精品日韩av在线免费观看| 精品国产乱码久久久久久男人| 最近最新中文字幕大全免费视频| 亚洲五月天丁香| 后天国语完整版免费观看| 国产亚洲欧美98| 黄色丝袜av网址大全| 日韩大尺度精品在线看网址| 国产伦一二天堂av在线观看| 欧美日韩精品网址| 日韩欧美 国产精品| 好男人电影高清在线观看| 成人三级黄色视频| 久久婷婷人人爽人人干人人爱| 波多野结衣高清作品| av在线天堂中文字幕| 久久精品国产亚洲av高清一级| 国产亚洲av嫩草精品影院| 久久久国产成人免费| 亚洲免费av在线视频| xxxwww97欧美| 91麻豆av在线| 50天的宝宝边吃奶边哭怎么回事| 欧美乱码精品一区二区三区| 欧美三级亚洲精品| 久久人人精品亚洲av| 久久精品夜夜夜夜夜久久蜜豆 | 老汉色av国产亚洲站长工具| 九九热线精品视视频播放| 夜夜爽天天搞| 真人一进一出gif抽搐免费| 欧美人与性动交α欧美精品济南到| 三级毛片av免费| 久久草成人影院| 日本一本二区三区精品| 欧美日本视频| 婷婷精品国产亚洲av| 亚洲 国产 在线| avwww免费| 十八禁网站免费在线| 国产精品爽爽va在线观看网站| 亚洲一码二码三码区别大吗| 好男人在线观看高清免费视频| 天堂av国产一区二区熟女人妻 | 窝窝影院91人妻| 国产亚洲精品综合一区在线观看 | 狠狠狠狠99中文字幕| 高清毛片免费观看视频网站| 成人国语在线视频| 中文字幕久久专区| 9191精品国产免费久久| 日日爽夜夜爽网站| 亚洲自偷自拍图片 自拍| 一级黄色大片毛片| 人人妻,人人澡人人爽秒播| 久久这里只有精品中国| 九九热线精品视视频播放| 日本 欧美在线| 久久午夜亚洲精品久久| x7x7x7水蜜桃| 正在播放国产对白刺激| 国产伦一二天堂av在线观看| 精品国产亚洲在线| 国产亚洲精品第一综合不卡| 88av欧美| 最近最新免费中文字幕在线| aaaaa片日本免费| 国产成人欧美在线观看| 日日夜夜操网爽| 亚洲va日本ⅴa欧美va伊人久久| 国产精品99久久99久久久不卡| 老司机靠b影院| 99re在线观看精品视频| 精品久久久久久久毛片微露脸| 色播亚洲综合网| 在线观看免费日韩欧美大片| 国产97色在线日韩免费| 中文字幕最新亚洲高清| 搞女人的毛片| 两个人的视频大全免费| 欧美黑人精品巨大| 露出奶头的视频| 18美女黄网站色大片免费观看| 色综合欧美亚洲国产小说| 国产99久久九九免费精品| 精品高清国产在线一区| 婷婷亚洲欧美| 伊人久久大香线蕉亚洲五| avwww免费| or卡值多少钱| 久久久国产成人免费| 午夜成年电影在线免费观看| 一本综合久久免费| 人妻夜夜爽99麻豆av| 亚洲成人中文字幕在线播放| 国产一区二区激情短视频| 日韩欧美免费精品| 国产黄色小视频在线观看| 十八禁人妻一区二区| 中文字幕精品亚洲无线码一区| 观看免费一级毛片| 成人永久免费在线观看视频| 曰老女人黄片| 999精品在线视频| 欧美黑人精品巨大| 久久中文字幕一级| 禁无遮挡网站| 熟女少妇亚洲综合色aaa.| 男男h啪啪无遮挡| 全区人妻精品视频| av欧美777| 色综合亚洲欧美另类图片| 亚洲av美国av| 狠狠狠狠99中文字幕| 非洲黑人性xxxx精品又粗又长| 欧美成狂野欧美在线观看| 老司机午夜福利在线观看视频| xxxwww97欧美| 99精品在免费线老司机午夜| 午夜两性在线视频| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 国产成人aa在线观看| 国产精品一区二区免费欧美| 搞女人的毛片| 日本五十路高清| 女生性感内裤真人,穿戴方法视频| 18美女黄网站色大片免费观看| 亚洲第一欧美日韩一区二区三区| 天堂av国产一区二区熟女人妻 | 色在线成人网| 久久久精品欧美日韩精品| 最近在线观看免费完整版| 久久国产精品人妻蜜桃| 在线观看66精品国产| 美女免费视频网站| 天天添夜夜摸| 久久这里只有精品19| 亚洲 欧美一区二区三区| 又紧又爽又黄一区二区| 美女午夜性视频免费| 精品国产乱码久久久久久男人| 精品熟女少妇八av免费久了| 亚洲人成网站高清观看| 精品少妇一区二区三区视频日本电影| 亚洲国产精品sss在线观看| 精品久久久久久久毛片微露脸| 在线观看免费午夜福利视频| 久久久久久久久免费视频了| 国产精品99久久99久久久不卡| 国产一区二区在线观看日韩 | 18禁黄网站禁片午夜丰满| 黄色a级毛片大全视频| 国产一区二区三区视频了| 国产精品一及| 1024视频免费在线观看| 蜜桃久久精品国产亚洲av| 日韩精品中文字幕看吧| 亚洲成人精品中文字幕电影| 久久久精品欧美日韩精品| 美女 人体艺术 gogo| 国产成人啪精品午夜网站| 不卡av一区二区三区| 97人妻精品一区二区三区麻豆| 香蕉丝袜av| 久久久国产成人精品二区| 人人妻人人看人人澡| 色哟哟哟哟哟哟| 久久中文看片网| 国模一区二区三区四区视频 | 可以免费在线观看a视频的电影网站| 成人欧美大片| 国产高清视频在线播放一区| 男女床上黄色一级片免费看| 欧美成狂野欧美在线观看| 亚洲av成人精品一区久久| 亚洲中文字幕日韩| 1024香蕉在线观看| 男插女下体视频免费在线播放| 久久午夜综合久久蜜桃| 国产视频一区二区在线看| 久久天堂一区二区三区四区| 人成视频在线观看免费观看| 蜜桃久久精品国产亚洲av| 国产精品国产高清国产av| 日本 av在线| 麻豆av在线久日| 天堂影院成人在线观看| 国产精品久久久久久久电影 | 亚洲第一欧美日韩一区二区三区| 熟妇人妻久久中文字幕3abv| 男女午夜视频在线观看| 国产精品自产拍在线观看55亚洲| 老司机午夜十八禁免费视频| 亚洲国产日韩欧美精品在线观看 | 无人区码免费观看不卡| 观看免费一级毛片| 亚洲 国产 在线| 日韩欧美三级三区| 欧美黑人巨大hd| 一本精品99久久精品77| 动漫黄色视频在线观看| 久久这里只有精品中国| 午夜福利在线观看吧| 两个人看的免费小视频| 999精品在线视频| 国产欧美日韩精品亚洲av| 欧美色视频一区免费| 欧美最黄视频在线播放免费| 999久久久国产精品视频| 琪琪午夜伦伦电影理论片6080| 亚洲欧美日韩高清专用| 淫妇啪啪啪对白视频| 97碰自拍视频| 天堂√8在线中文| 深夜精品福利| 中出人妻视频一区二区| 激情在线观看视频在线高清| 神马国产精品三级电影在线观看 | 国产精品免费视频内射| 国产又色又爽无遮挡免费看| 亚洲欧美一区二区三区黑人| 国产精品久久电影中文字幕| 国产精品免费视频内射| 巨乳人妻的诱惑在线观看| 亚洲av熟女| 国产精品久久电影中文字幕| 麻豆成人av在线观看| 免费在线观看完整版高清| 午夜精品在线福利| 脱女人内裤的视频| 国产精品久久久久久亚洲av鲁大| 国产一区二区三区在线臀色熟女| 夜夜夜夜夜久久久久| 成人高潮视频无遮挡免费网站| 国产精品久久久av美女十八| 欧美成人性av电影在线观看| 国产精品 欧美亚洲| 亚洲中文字幕一区二区三区有码在线看 | 亚洲七黄色美女视频| 国产视频内射| 欧美激情久久久久久爽电影| 久久这里只有精品中国| 日本一二三区视频观看| 国产99久久九九免费精品| 日韩精品青青久久久久久| 欧美日韩乱码在线| 亚洲avbb在线观看| 久热爱精品视频在线9| 88av欧美| 非洲黑人性xxxx精品又粗又长| 国产精品亚洲美女久久久| 啦啦啦免费观看视频1| 母亲3免费完整高清在线观看| 国产在线精品亚洲第一网站| 日本免费一区二区三区高清不卡| 国产精品自产拍在线观看55亚洲| 夜夜夜夜夜久久久久| 国产黄片美女视频| 五月玫瑰六月丁香| 亚洲人成电影免费在线| 一本综合久久免费| 丰满的人妻完整版| 亚洲人成电影免费在线| 一本大道久久a久久精品| 99久久精品热视频| 一级作爱视频免费观看| 欧美 亚洲 国产 日韩一| 香蕉av资源在线| 无人区码免费观看不卡| 成人高潮视频无遮挡免费网站| 国产精品一及| 精品无人区乱码1区二区| 手机成人av网站| 国产黄a三级三级三级人| 久久久久性生活片| 99精品在免费线老司机午夜| 日韩成人在线观看一区二区三区| 一级片免费观看大全| 91老司机精品| 亚洲国产精品999在线| 高清毛片免费观看视频网站| 国产成人精品久久二区二区免费| 国产一区二区在线观看日韩 | 国产人伦9x9x在线观看| 三级国产精品欧美在线观看 | 国产一区二区在线av高清观看| 久久草成人影院| 动漫黄色视频在线观看| 欧美性猛交╳xxx乱大交人| 国产精品 欧美亚洲| 我的老师免费观看完整版| 18禁美女被吸乳视频| 亚洲人成网站高清观看| 亚洲乱码一区二区免费版| 又大又爽又粗| 亚洲午夜理论影院| 国产精品久久久久久精品电影| 亚洲人与动物交配视频| 日韩国内少妇激情av| 宅男免费午夜| 日本撒尿小便嘘嘘汇集6| 国产乱人伦免费视频| 久久香蕉国产精品| 美女扒开内裤让男人捅视频| 亚洲 欧美 日韩 在线 免费| 极品教师在线免费播放| 欧美大码av| 无遮挡黄片免费观看| 一本久久中文字幕| 在线观看舔阴道视频| 免费电影在线观看免费观看| 夜夜夜夜夜久久久久| 国产视频一区二区在线看| 亚洲人成77777在线视频| 日韩欧美国产在线观看| 999久久久精品免费观看国产| 国产不卡一卡二| 亚洲人成电影免费在线| 白带黄色成豆腐渣| videosex国产| 国产片内射在线| 国产精品久久久久久亚洲av鲁大| 久久久久久久精品吃奶| 国产av麻豆久久久久久久| avwww免费| 99久久国产精品久久久| 黄色成人免费大全| videosex国产| 中文在线观看免费www的网站 | 久久九九热精品免费| 一本综合久久免费| 免费电影在线观看免费观看| 婷婷丁香在线五月| 少妇粗大呻吟视频| cao死你这个sao货| 天堂av国产一区二区熟女人妻 | 国产成人精品无人区| 又黄又粗又硬又大视频| 亚洲性夜色夜夜综合| 免费在线观看日本一区| 少妇被粗大的猛进出69影院| 亚洲人成伊人成综合网2020| 国产亚洲av高清不卡| 午夜久久久久精精品| 伊人久久大香线蕉亚洲五| 久久九九热精品免费| 9191精品国产免费久久| 亚洲精品美女久久av网站| 国产精品亚洲av一区麻豆| 老汉色∧v一级毛片| 日日干狠狠操夜夜爽| 国产探花在线观看一区二区| 男人舔女人下体高潮全视频| 国产精品美女特级片免费视频播放器 | 十八禁人妻一区二区| 国产精品乱码一区二三区的特点| 久久久水蜜桃国产精品网| 亚洲男人的天堂狠狠| 嫩草影视91久久| 1024视频免费在线观看| 无限看片的www在线观看| 国产黄a三级三级三级人| 国产av一区在线观看免费| 嫩草影视91久久| 在线观看免费日韩欧美大片| 日本熟妇午夜| 丝袜美腿诱惑在线| 一级毛片精品| 长腿黑丝高跟| 久久精品国产99精品国产亚洲性色| 动漫黄色视频在线观看| 久久亚洲真实| 亚洲va日本ⅴa欧美va伊人久久| 日韩三级视频一区二区三区| 成熟少妇高潮喷水视频| 好男人电影高清在线观看| 香蕉国产在线看| 日本成人三级电影网站| 免费人成视频x8x8入口观看| 99久久99久久久精品蜜桃| 久久香蕉国产精品| 九色国产91popny在线| 国产av一区二区精品久久| 亚洲第一欧美日韩一区二区三区| 神马国产精品三级电影在线观看 | 白带黄色成豆腐渣| 欧美日韩乱码在线| 女人爽到高潮嗷嗷叫在线视频| 国产三级黄色录像| 日韩大尺度精品在线看网址| 久久精品国产清高在天天线| 蜜桃久久精品国产亚洲av| 美女免费视频网站| 在线看三级毛片| 18禁观看日本| 别揉我奶头~嗯~啊~动态视频| 国产午夜福利久久久久久| 丝袜人妻中文字幕| 国产精品一区二区三区四区免费观看 | 婷婷精品国产亚洲av在线| 真人一进一出gif抽搐免费| 亚洲熟女毛片儿| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| 搡老妇女老女人老熟妇| 中文资源天堂在线| a级毛片a级免费在线| 中文字幕久久专区| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三| www.精华液| www日本黄色视频网| 性欧美人与动物交配| 亚洲美女黄片视频| 热99re8久久精品国产| 嫩草影视91久久| 在线观看免费午夜福利视频| 日韩欧美 国产精品| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 一本一本综合久久| 国产一区二区三区视频了| 青草久久国产| 国产午夜福利久久久久久| 午夜精品在线福利| 欧美一区二区精品小视频在线| 精品一区二区三区四区五区乱码| av天堂在线播放| 最近在线观看免费完整版| 亚洲一区二区三区色噜噜| 久久九九热精品免费| 在线观看免费视频日本深夜| 久久精品91无色码中文字幕| 白带黄色成豆腐渣| cao死你这个sao货| 老熟妇乱子伦视频在线观看| 99久久无色码亚洲精品果冻| 757午夜福利合集在线观看| 久久人妻福利社区极品人妻图片| 天堂av国产一区二区熟女人妻 | 国产精品久久视频播放| 国语自产精品视频在线第100页| 亚洲黑人精品在线| 777久久人妻少妇嫩草av网站| av超薄肉色丝袜交足视频| 精品福利观看| 最近最新免费中文字幕在线| 亚洲午夜精品一区,二区,三区| 日本a在线网址| 国产精品乱码一区二三区的特点| 51午夜福利影视在线观看| 亚洲欧美日韩东京热| av视频在线观看入口| 日韩精品中文字幕看吧| 大型av网站在线播放| 精品福利观看| 国产精品久久视频播放| 亚洲专区国产一区二区| 久久精品成人免费网站| 午夜两性在线视频| tocl精华| 欧美一区二区精品小视频在线| 亚洲成人久久爱视频| www.自偷自拍.com| 亚洲成a人片在线一区二区| 午夜福利18| 99久久综合精品五月天人人| av在线播放免费不卡| 久久伊人香网站| 国产精品,欧美在线| 夜夜躁狠狠躁天天躁| 国产三级黄色录像| 99久久99久久久精品蜜桃|