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

    一種預(yù)測(cè)頻段占用度的時(shí)間序列分析方法

    2011-09-26 02:00:12段洪濤曾繁聲李景春
    無(wú)線電工程 2011年7期
    關(guān)鍵詞:頻段頻譜檢驗(yàn)

    段洪濤,曾繁聲,李景春

    (國(guó)家無(wú)線電監(jiān)測(cè)中心,北京100037)

    0 引言

    無(wú)線電頻譜資源是稀缺資源,頻譜占用度測(cè)量數(shù)據(jù)對(duì)于無(wú)線電管理工作至關(guān)重要。在國(guó)際電聯(lián)(ITU)出版的文獻(xiàn)中,頻譜占用度的作用可以歸納為:反映頻段內(nèi)的頻率使用情況;為頻率的分配和指配提供參考依據(jù);反映一個(gè)地區(qū)特定頻段的頻率使用變化趨勢(shì)。自2007年開始,無(wú)線電管理部門重點(diǎn)開展了頻譜占用度監(jiān)測(cè)月報(bào)工作,對(duì)測(cè)試方法和統(tǒng)計(jì)方法都進(jìn)行了規(guī)范,上報(bào)數(shù)據(jù)的分析和挖掘工作也凸顯重要性。

    目前對(duì)于頻譜占用情況的掌握都是直接采用測(cè)量結(jié)果,這種測(cè)量一般都是根據(jù)實(shí)際工作目標(biāo)的短期測(cè)量行為,實(shí)際上是將測(cè)量結(jié)果作為現(xiàn)在和未來(lái)一段時(shí)間內(nèi)的頻譜占用情況,前提是認(rèn)為頻譜占用度情況平穩(wěn)變化[1,2]。如果長(zhǎng)期開展測(cè)量工作,則可以得到更為精確的預(yù)測(cè)結(jié)果和變化趨勢(shì)。

    1 頻譜占用度及時(shí)間序列

    1.1 頻率(頻譜)占用度測(cè)量

    信道占用度:使用監(jiān)測(cè)接收機(jī)或頻譜分析儀對(duì)特定的信道進(jìn)行測(cè)量,信號(hào)大于某一門限電平值的時(shí)間與總測(cè)量時(shí)間的百分比。

    頻段占用度:在測(cè)量時(shí)間內(nèi),用監(jiān)測(cè)接收機(jī)或頻譜分析儀對(duì)某一頻段用固定步長(zhǎng)(信道)進(jìn)行順序測(cè)量,大于某一門限電平值的信道數(shù)與總信道數(shù)的百分比[3]。

    在頻譜占用度測(cè)量工作中有2個(gè)方面的問(wèn)題需要深入研究:①測(cè)量本身的問(wèn)題,即如何測(cè)量頻譜占用度,在此方面應(yīng)保證測(cè)試過(guò)程的合理性、測(cè)試數(shù)據(jù)的真實(shí)性,測(cè)量結(jié)果的準(zhǔn)確性和可靠性,涉及信號(hào)檢測(cè)、儀器原理以及統(tǒng)計(jì)學(xué)的相關(guān)知識(shí);②頻段占用度數(shù)據(jù)分析的問(wèn)題,即怎樣分析這些占用度數(shù)據(jù),能得到什么結(jié)果,這個(gè)方面涉及數(shù)據(jù)分析和挖掘的相關(guān)知識(shí)。

    占用度數(shù)據(jù)的第1個(gè)作用通過(guò)測(cè)量本身體現(xiàn),第2個(gè)和第3個(gè)作用從嚴(yán)格意義上來(lái)講,都屬于占用度的預(yù)測(cè)和分析范疇。舉個(gè)例子:為了在某個(gè)頻段指配若干個(gè)頻率,進(jìn)行了一段時(shí)間的專項(xiàng)監(jiān)測(cè),得到了頻段占用度數(shù)據(jù),然后根據(jù)占用度數(shù)據(jù)進(jìn)行了頻率指配。要知道占用度數(shù)據(jù)僅是針對(duì)過(guò)去特定時(shí)間的統(tǒng)計(jì)值,如果根據(jù)其進(jìn)行頻率指配,則隱含了預(yù)測(cè)分析的內(nèi)容。

    頻率占用度測(cè)量的一般過(guò)程是:用接收機(jī)順序掃描頻段內(nèi)的若干個(gè)信道,在大量掃描后對(duì)每個(gè)信道都有一組測(cè)試數(shù)據(jù)作為分析基礎(chǔ)數(shù)據(jù),通常情況下這些測(cè)試數(shù)據(jù)在時(shí)間上都是等間隔分布的。這些基礎(chǔ)數(shù)據(jù)經(jīng)過(guò)統(tǒng)計(jì)分析形成了每個(gè)信道以及整個(gè)頻段的占用度數(shù)據(jù),這些占用度數(shù)據(jù)可以是測(cè)試時(shí)間內(nèi)的,通常是一天的;也可以是每周的,或者是每月的占用度統(tǒng)計(jì)數(shù)據(jù)[4]。

    1.2 時(shí)間序列

    時(shí)間序列是按時(shí)間次序排列的隨機(jī)變量序列。在實(shí)際問(wèn)題中所能得到的數(shù)據(jù)只是時(shí)間序列的有限樣本。時(shí)間序列分析的主要任務(wù)就是根據(jù)觀測(cè)數(shù)據(jù)的特點(diǎn)為數(shù)據(jù)建立盡可能合理的統(tǒng)計(jì)模型,然后利用模型的統(tǒng)計(jì)特性去解釋數(shù)據(jù)的統(tǒng)計(jì)規(guī)律,以期達(dá)到控制或預(yù)報(bào)目的。

    在頻率占用度的測(cè)量中,無(wú)論從基礎(chǔ)數(shù)據(jù)獲得的微觀時(shí)間角度,還是從每月上報(bào)數(shù)據(jù)的宏觀時(shí)間角度,這些占用度數(shù)據(jù)都自然構(gòu)成了時(shí)間序列[5]。需要說(shuō)明的是時(shí)間序列分析對(duì)數(shù)據(jù)也有一定的要求,具體地說(shuō)就是數(shù)據(jù)的數(shù)量和質(zhì)量,數(shù)據(jù)數(shù)量太少則不能完全反映數(shù)據(jù)之間的聯(lián)系,一般來(lái)說(shuō)常見的時(shí)間序列分析數(shù)據(jù)數(shù)量都在20個(gè)以上。數(shù)據(jù)的質(zhì)量指的是數(shù)據(jù)的獲得應(yīng)該在同一條件下。對(duì)于頻譜占用度的時(shí)序分析來(lái)說(shuō),這一點(diǎn)尤其重要,因?yàn)闇y(cè)試數(shù)據(jù)由于測(cè)試設(shè)備、門限電平、統(tǒng)計(jì)方法、甚至是測(cè)試人員的不同而有較大的變化。計(jì)算所使用的占用度數(shù)據(jù)是在連續(xù)的一段時(shí)間按內(nèi)由同一位測(cè)試人員使用同一套設(shè)備獲得的,基本可以認(rèn)為測(cè)試條件相同。2009年北京某站上報(bào)的針對(duì)450~470MHz測(cè)試的頻段占用度上報(bào)數(shù)據(jù)如表1所示。

    表1 450~470 MHz頻段占用度數(shù)據(jù)

    任何時(shí)間序列經(jīng)過(guò)合理的函數(shù)變換后都可以認(rèn)為是由3個(gè)部分疊加而成,這3個(gè)部分是趨勢(shì)項(xiàng)部分、周期項(xiàng)部分和隨機(jī)噪聲項(xiàng)部分,表示為:

    式中,{Tt}為趨勢(shì)項(xiàng);{St}為季節(jié)項(xiàng);{Rt}為隨機(jī)項(xiàng)。時(shí)間序列{Xt}是這3項(xiàng)的疊加。時(shí)間序列分析的首要任務(wù)是通過(guò)對(duì)觀測(cè)樣本的觀察分析,把時(shí)間序列的趨勢(shì)項(xiàng)、季節(jié)項(xiàng)和隨機(jī)項(xiàng)分解出來(lái)。對(duì)于趨勢(shì)項(xiàng)和季節(jié)項(xiàng),一般可以用非隨機(jī)的函數(shù)進(jìn)行刻畫,比較簡(jiǎn)單。當(dāng)時(shí)間序列在適當(dāng)?shù)娜サ糈厔?shì)項(xiàng)和季節(jié)項(xiàng)后,剩下的隨機(jī)部分通常會(huì)有某種平穩(wěn)性,而平穩(wěn)時(shí)間序列的歷史X1,X2,…Xn中往往含有Xn+1的信息,這就使得利用歷史樣本X1,X2,…Xn預(yù)測(cè)將來(lái)Xn+1成為可能,這也是時(shí)間序列進(jìn)行預(yù)測(cè)的本質(zhì)[5]。

    2 ARIMA模型分析

    ARIMA模型有3種基本類型:自回歸(AR)模型、移動(dòng)平均(MA)模型和自回歸移動(dòng)平均(ARIMA)模型。

    2.1 自回歸模型

    如果時(shí)間序列yt是它的前期值和隨機(jī)項(xiàng)的線性函數(shù),即可表示為:

    則稱該時(shí)間序列是自回歸序列,p為自回歸模型的階數(shù),記為AR(p),隨機(jī)項(xiàng)ut是相互獨(dú)立的白噪聲序列,記 Bk為k步滯后算子,即

    則模型可以表示為:

    模型可以簡(jiǎn)寫為:

    AR(p)過(guò)程平穩(wěn)的條件是滯后多項(xiàng)式 φ(B)的根均在單位圓外,即 φ(B)=0的根大于1。

    2.2 移動(dòng)平均模型

    如果時(shí)間序列 yt是它的當(dāng)期和前期的隨機(jī)誤差項(xiàng)的線性函數(shù),即可表示為:

    則稱該時(shí)間序列是移動(dòng)平均序列,q為移動(dòng)平均模型的階數(shù),記作MA(q)。

    引入滯后算子,并令

    yt可以簡(jiǎn)寫為:

    移動(dòng)平均過(guò)程無(wú)條件平穩(wěn)要求,但通常希望AR過(guò)程與MA過(guò)程能相互表出,即為可逆過(guò)程。因此要求滯后多項(xiàng)式θ(B)的根都在單位圓外。

    2.3 自回歸移動(dòng)平均模型

    如果時(shí)間序列yt是它的當(dāng)期和前期的隨機(jī)誤差項(xiàng)以及前期值的線性函數(shù),即可表示為:

    則稱該時(shí)間序列是自回歸移動(dòng)平均序列,記作ARIMA(p,q),該模型是由博克思(Box)和詹金斯(Jenkins)于20世紀(jì)70年代初提出的一著名時(shí)間序列預(yù)測(cè)方法,所以又稱為box-jenkins模型、博克思-詹金斯法。這種方法具有較高的預(yù)測(cè)精度,適用于短期預(yù)測(cè)。其中ARIMA(p,d,q)稱為差分自回歸移動(dòng)平均模型,AR為自回歸,p為自回歸項(xiàng);MA為移動(dòng)平均,q為移動(dòng)平均項(xiàng)數(shù),d為時(shí)間序列成為平穩(wěn)時(shí)所做的差分次數(shù)[6]。

    3 性能分析

    下面將分成以下幾個(gè)步驟對(duì)上面所提模型進(jìn)行性能分析:①時(shí)間序列的平穩(wěn)性識(shí)別;②對(duì)非平穩(wěn)的時(shí)間序列數(shù)據(jù)進(jìn)行平穩(wěn)化處理;③選取適用的模型及參數(shù);④模型檢驗(yàn);⑤利用已通過(guò)檢驗(yàn)的模型進(jìn)行預(yù)測(cè)。

    為了檢驗(yàn)預(yù)測(cè)效果,將第10個(gè)觀測(cè)值留出,作為評(píng)價(jià)預(yù)測(cè)精度的參照對(duì)象。建模的樣本期為2009年3月—2009年11月,在具體應(yīng)用模型過(guò)程中可以使用EViews軟件幫助計(jì)算和分析。

    3.1 平穩(wěn)性識(shí)別

    平穩(wěn)時(shí)間序列的一般概念是時(shí)間序列的統(tǒng)計(jì)特征不隨時(shí)間變化而變化,寬平穩(wěn)的定義是:序列的均值和方差都是與時(shí)間 t無(wú)關(guān)的常數(shù),協(xié)方差僅與時(shí)間間隔有關(guān),與t無(wú)關(guān)。關(guān)于頻譜占用度的預(yù)測(cè)問(wèn)題首先要從判定時(shí)間序列的平穩(wěn)性開始,單從定義來(lái)看序列的平穩(wěn)性比較抽象,不容易判斷,其實(shí)在實(shí)際工作中時(shí)間序列的平穩(wěn)性有完備的判斷理論和簡(jiǎn)單方法。

    最簡(jiǎn)單的方法是通過(guò)時(shí)間序列的時(shí)序圖來(lái)判斷時(shí)間序列是否存在周期性或趨勢(shì)性。這種方法的優(yōu)點(diǎn)是簡(jiǎn)單直觀,對(duì)于那些明顯是非平穩(wěn)的時(shí)間序列很容易判斷。缺點(diǎn)是對(duì)于一般的時(shí)間序列,特別是不帶有明顯特征的時(shí)間序列,不容易判斷準(zhǔn)確。占用度隨時(shí)間變換的時(shí)序圖如圖1所示。

    從圖1可以明顯看出占用度呈現(xiàn)緩慢上升的趨勢(shì),應(yīng)該不是平穩(wěn)序列。

    圖1 時(shí)序圖

    判斷時(shí)間序列平穩(wěn)性比較理論化的方法是單位根檢驗(yàn)(Unit Root Test,或稱為迪基——福勒檢驗(yàn))。這種方法在實(shí)際中比較常用,但是數(shù)據(jù)計(jì)算量較大??梢岳脮r(shí)間序列分析中常用的分析軟件(Eviews)對(duì)表1的數(shù)據(jù)進(jìn)行單位根的計(jì)算。

    經(jīng)過(guò)計(jì)算:時(shí)間序列檢驗(yàn)統(tǒng)計(jì)值為-0.983 937,大于顯著水平為10%的臨界值-2.801 384,所以不能拒絕原假設(shè),序列存在單位根,是非平穩(wěn)的。

    3.2 平穩(wěn)化處理

    通過(guò)時(shí)序圖和單位根檢驗(yàn),基本可判斷原序列不平穩(wěn),下面就是對(duì)序列進(jìn)行平穩(wěn)化處理。在時(shí)間序列分析中,如果一個(gè)時(shí)間序列是非平穩(wěn)的,更一般的方法是通過(guò)取差分的方法而形成平穩(wěn)序列。

    經(jīng)過(guò)計(jì)算:一階差分后的時(shí)間序列在1%的顯著水平下拒絕原假設(shè),即99%的置信水平下不存在單位根,該序列是平穩(wěn)的。

    3.3 模型識(shí)別及參數(shù)估計(jì)

    所謂隨機(jī)時(shí)間序列模型的識(shí)別,就是對(duì)于一個(gè)平穩(wěn)的隨機(jī)時(shí)間序列,找出生成它的合適的隨機(jī)過(guò)程或模型,即判斷該時(shí)間序列是遵循一純AR過(guò)程、還是遵循一純MA過(guò)程或ARIMA過(guò)程。

    所使用的工具主要是時(shí)間序列的自相關(guān)函數(shù)(Autocorrelation Function,ACF)及偏自相關(guān)函數(shù)(Partial Autocorrelation Function,PACF)。ARIMA(p,q)的自相關(guān)函數(shù)可以看作MA(q)的自相關(guān)函數(shù)和AR(p)的自相關(guān)函數(shù)的混合物。

    當(dāng)p=0時(shí),具有截尾性質(zhì);

    當(dāng)q=0時(shí),具有拖尾性質(zhì);

    當(dāng)p、q都不為0時(shí),具有拖尾性質(zhì)。

    從識(shí)別上看,通常ARIMA(p,q)過(guò)程的偏自相關(guān)函數(shù)(PACF)可能在p階滯后前有幾項(xiàng)明顯的尖柱(spikes),但從 p階滯后項(xiàng)開始逐漸趨向于零;而它的自相關(guān)函數(shù)(ACF)則是在q階滯后前有幾項(xiàng)明顯的尖柱,從q階滯后項(xiàng)開始逐漸趨向于零。一階差分后的時(shí)間序列自相關(guān)與偏相關(guān)計(jì)算結(jié)果如表2所示。

    表2 一階差分時(shí)間序列相關(guān)性計(jì)算結(jié)果

    模型的參數(shù)估計(jì)和檢驗(yàn)結(jié)果如表3所示。由表3可以看出自相關(guān)函數(shù)在滯后2期和3期都明顯不為零,之后逐漸表現(xiàn)為拖尾性,偏相關(guān)函數(shù)也在滯后2期后表現(xiàn)為拖尾特性。所以對(duì)于差分后的時(shí)間序列可初步選用 ARIMA(2,1)、ARIMA(2,2)、ARIMA(3,1)和ARIMA(3,2)模型。

    表3 模型的參數(shù)估計(jì)和檢驗(yàn)結(jié)果

    從表3可以看出,ARIMA(3,1)模型的多項(xiàng)式倒數(shù)根都落在單位圓內(nèi),滿足過(guò)程平穩(wěn)的基本要求,其他模型的倒數(shù)跟都出現(xiàn)了在單位圓外的情況,所以可以認(rèn)為ARIMA(3,1)模型是合適的。

    3.4 模型檢驗(yàn)

    模型檢驗(yàn)主要是指對(duì)模型的合適性進(jìn)行驗(yàn)證,通常是對(duì)模型的殘差序列et進(jìn)行白噪聲檢驗(yàn)。若殘差序列不是白噪聲序列,意味著殘差序列還存在有用信息沒(méi)被提取,還需要進(jìn)一步改進(jìn)模型。通常側(cè)重于檢驗(yàn)殘差序列的隨機(jī)性。殘差序列樣本的自相關(guān)系數(shù)應(yīng)該近似為零。殘差自相關(guān)檢驗(yàn)結(jié)果如表4所示。表4表明殘差近于白噪聲,模型已經(jīng)提取了有規(guī)律的信息,模型的擬合效果較好。

    表4 殘差自相關(guān)檢驗(yàn)結(jié)果

    3.5 模型預(yù)測(cè)

    根據(jù)ARIMA(3,1)模型對(duì)頻段占用度進(jìn)行預(yù)測(cè)結(jié)果如圖2所示。實(shí)線表示預(yù)測(cè)值,虛線表示預(yù)測(cè)值2倍標(biāo)準(zhǔn)誤差帶。經(jīng)過(guò)計(jì)算,均方根誤差為0.188 120,平均絕對(duì)誤差為0.162 667,平均絕對(duì)百分誤差為59.639 49,表明預(yù)測(cè)精度較好。預(yù)測(cè)的12月份頻段占用度為14.094 276,而實(shí)際的測(cè)試值為14.13,可以看出預(yù)測(cè)值很準(zhǔn)確。

    圖2 占用度預(yù)測(cè)圖

    4 結(jié)束語(yǔ)

    通過(guò)實(shí)際運(yùn)用可以看出:對(duì)于頻段占用度測(cè)量工作,測(cè)試及統(tǒng)計(jì)條件的一致性能是后續(xù)分析的基礎(chǔ),也是發(fā)現(xiàn)頻段占用度精確變化規(guī)律的前提條件,沒(méi)有準(zhǔn)確的采集數(shù)據(jù)也就沒(méi)有后續(xù)的有價(jià)值的分析。就拿本例來(lái)說(shuō),該站后來(lái)更換了靈敏度更高的設(shè)備,450~470MHz頻段的占用度測(cè)量結(jié)果大幅度增加到50%左右,這種變化不是由于頻率使用量激增的結(jié)果,而是由于能夠更精確地區(qū)分小信號(hào)和噪聲的結(jié)果。

    長(zhǎng)時(shí)間積累數(shù)據(jù)很有價(jià)值。對(duì)于本例的分析,原想找到更多的數(shù)據(jù)進(jìn)行分析,因?yàn)檎加枚葴y(cè)量工作具有一定的階段性,每次任務(wù)時(shí)間一般是幾個(gè)月到半年,數(shù)據(jù)樣本點(diǎn)比較少,通常間隔一段時(shí)間后占用度數(shù)據(jù)還會(huì)有比較明顯的變化,經(jīng)過(guò)咨詢多數(shù)情況是測(cè)試人員變化造成的,因?yàn)殚T限電平設(shè)置、信道被占用判斷都是根據(jù)經(jīng)驗(yàn),可以考慮用軟件固定設(shè)置值解決此類問(wèn)題。

    利用EViews分析軟件對(duì)監(jiān)測(cè)站上報(bào)的頻段占用度監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析,重點(diǎn)介紹了ARIMA模型在分析和預(yù)測(cè)頻段占用度數(shù)據(jù)方面的應(yīng)用過(guò)程。關(guān)于頻率占用度數(shù)據(jù)的分析、挖掘有很多東西需要研究,特別是結(jié)合頻段劃分和業(yè)務(wù)劃分的研究,如果有可能結(jié)合臺(tái)站數(shù)據(jù)庫(kù)進(jìn)行分析比對(duì)可以更好地預(yù)測(cè)占用度變化情況,體現(xiàn)日常監(jiān)測(cè)工作的價(jià)值,對(duì)無(wú)線電管理工作發(fā)揮更大的作用。

    [1]MCHENRY MD.ROBERSON M,D,MACDONALDJ.Chicago Spectrum Occupancy Measurements[R].Report for the WIL of IIT,2006.

    [2]ELLINGSON S W.Spectral occupancy at VHF:Implication for Frequency-agile Cognitive Radios[C].In Proc.of IEEE Vehicular Technology Conference(VTC),2005:1397-1382.

    [3]SPAULDING A D,HAGN,G H.On the Definition and Estimation of Spectrum Occupancy[J].IEEE Trans.on EMC,1977,19(3):269-280.

    [4]CEPT/ERC Recommendation 01-10 E.Frequency Channel Occupancy Measurements[S],2000.

    [5]朱照宣.關(guān)于時(shí)間序列分析[J].力學(xué)實(shí)踐,1989,11(1):22-26.

    [6]何書元.應(yīng)用時(shí)間序列分析[M].北京:北京大學(xué)出版社,2009.

    猜你喜歡
    頻段頻譜檢驗(yàn)
    序貫Lq似然比型檢驗(yàn)
    gPhone重力儀的面波頻段響應(yīng)實(shí)測(cè)研究
    地震研究(2021年1期)2021-04-13 01:04:56
    一種用于深空探測(cè)的Chirp變換頻譜分析儀設(shè)計(jì)與實(shí)現(xiàn)
    2021年《理化檢驗(yàn)-化學(xué)分冊(cè)》征訂啟事
    對(duì)起重機(jī)“制動(dòng)下滑量”相關(guān)檢驗(yàn)要求的探討
    一種基于稀疏度估計(jì)的自適應(yīng)壓縮頻譜感知算法
    關(guān)于鍋爐檢驗(yàn)的探討
    推擠的5GHz頻段
    CHIP新電腦(2016年3期)2016-03-10 14:07:52
    認(rèn)知無(wú)線電頻譜感知技術(shù)綜述
    TD—LTE在D頻段和F頻段的覆蓋能力差異
    在线观看午夜福利视频| www.精华液| 久久人妻熟女aⅴ| 国产精品九九99| 亚洲精品美女久久av网站| 婷婷六月久久综合丁香| 久久亚洲精品不卡| 免费av毛片视频| 一级a爱片免费观看的视频| 丁香欧美五月| 黄色女人牲交| 国产又爽黄色视频| 欧美午夜高清在线| 99在线人妻在线中文字幕| 一二三四社区在线视频社区8| 88av欧美| 亚洲激情在线av| 亚洲第一青青草原| 国产精品自产拍在线观看55亚洲| 老司机午夜十八禁免费视频| 夜夜躁狠狠躁天天躁| 成人精品一区二区免费| 黄色毛片三级朝国网站| 精品久久久久久久久久免费视频| 久9热在线精品视频| 国产成年人精品一区二区| 久久国产乱子伦精品免费另类| 亚洲色图av天堂| 午夜福利影视在线免费观看| 此物有八面人人有两片| 大陆偷拍与自拍| 国产精品久久视频播放| 大香蕉久久成人网| 国产三级黄色录像| 国产又爽黄色视频| 国产亚洲精品久久久久久毛片| 国产熟女xx| 国产亚洲欧美精品永久| 91大片在线观看| 亚洲色图av天堂| av在线播放免费不卡| 精品免费久久久久久久清纯| 亚洲欧美精品综合久久99| 咕卡用的链子| 亚洲熟女毛片儿| 激情视频va一区二区三区| 亚洲欧美日韩高清在线视频| 性少妇av在线| 在线观看午夜福利视频| 日本精品一区二区三区蜜桃| 亚洲中文字幕日韩| 日本一区二区免费在线视频| 黄片小视频在线播放| 久久久久国产一级毛片高清牌| 99在线视频只有这里精品首页| 91大片在线观看| 免费在线观看日本一区| 精品国产乱子伦一区二区三区| 一本综合久久免费| 亚洲av五月六月丁香网| 国产精品99久久99久久久不卡| 亚洲 欧美 日韩 在线 免费| 97碰自拍视频| 国产午夜精品久久久久久| 真人一进一出gif抽搐免费| 亚洲国产精品久久男人天堂| 高清在线国产一区| 淫妇啪啪啪对白视频| 亚洲自偷自拍图片 自拍| 黑人巨大精品欧美一区二区mp4| 久久久久久人人人人人| 精品国产美女av久久久久小说| 18禁国产床啪视频网站| av福利片在线| 咕卡用的链子| 国产主播在线观看一区二区| 亚洲五月婷婷丁香| 国产激情欧美一区二区| 久久人人97超碰香蕉20202| av有码第一页| 成年女人毛片免费观看观看9| 好看av亚洲va欧美ⅴa在| 日本五十路高清| 国产三级黄色录像| 亚洲av日韩精品久久久久久密| 国产色视频综合| 最近最新中文字幕大全电影3 | 国内精品久久久久久久电影| 老司机深夜福利视频在线观看| 色在线成人网| 免费观看人在逋| 少妇的丰满在线观看| 国产av一区在线观看免费| 91老司机精品| 此物有八面人人有两片| 欧美在线一区亚洲| 欧美大码av| 国产视频一区二区在线看| 免费在线观看完整版高清| 成人手机av| 精品一区二区三区av网在线观看| 日韩精品青青久久久久久| 99国产综合亚洲精品| 亚洲免费av在线视频| 中文字幕精品免费在线观看视频| 国产又色又爽无遮挡免费看| 日韩高清综合在线| 国内精品久久久久精免费| 国产欧美日韩一区二区三区在线| 亚洲无线在线观看| 在线av久久热| 亚洲av成人一区二区三| 黄网站色视频无遮挡免费观看| 亚洲,欧美精品.| 久久性视频一级片| 搡老熟女国产l中国老女人| 国产成人系列免费观看| 久久精品成人免费网站| 亚洲精品av麻豆狂野| 亚洲 国产 在线| 免费高清在线观看日韩| 国产激情久久老熟女| 久久国产精品影院| 国产亚洲欧美在线一区二区| 色综合欧美亚洲国产小说| 欧美日韩精品网址| 久久久水蜜桃国产精品网| 无人区码免费观看不卡| 老司机在亚洲福利影院| 久久久水蜜桃国产精品网| 精品一区二区三区四区五区乱码| 97人妻天天添夜夜摸| 欧美老熟妇乱子伦牲交| 老司机午夜十八禁免费视频| 97人妻天天添夜夜摸| 久久中文看片网| 欧美人与性动交α欧美精品济南到| 精品久久久久久久人妻蜜臀av | 两性夫妻黄色片| 国产99久久九九免费精品| 欧美+亚洲+日韩+国产| 免费看美女性在线毛片视频| 日韩欧美国产在线观看| 亚洲 国产 在线| 午夜亚洲福利在线播放| 国产成人av激情在线播放| 999久久久国产精品视频| 91av网站免费观看| 国产日韩一区二区三区精品不卡| 亚洲伊人色综图| 日日干狠狠操夜夜爽| 欧美国产日韩亚洲一区| 成人国产一区最新在线观看| 日本 欧美在线| 人人澡人人妻人| 亚洲 国产 在线| а√天堂www在线а√下载| www.熟女人妻精品国产| 国产精品久久久久久亚洲av鲁大| 久久久久国内视频| 一本综合久久免费| 日本精品一区二区三区蜜桃| 亚洲人成伊人成综合网2020| 人人澡人人妻人| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人精品中文字幕电影| 男人的好看免费观看在线视频 | 法律面前人人平等表现在哪些方面| 一区二区三区激情视频| 黑丝袜美女国产一区| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| www.www免费av| 国产1区2区3区精品| 亚洲成人精品中文字幕电影| 后天国语完整版免费观看| 久久久久久久久中文| 亚洲国产毛片av蜜桃av| 一本久久中文字幕| 日本在线视频免费播放| 国产私拍福利视频在线观看| 嫩草影视91久久| 岛国视频午夜一区免费看| 国产高清有码在线观看视频 | 欧美精品亚洲一区二区| 女人高潮潮喷娇喘18禁视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲天堂国产精品一区在线| 国产精品久久久av美女十八| 国产色视频综合| 精品高清国产在线一区| 亚洲第一青青草原| 91麻豆av在线| 欧美日韩黄片免| 久久中文字幕一级| 50天的宝宝边吃奶边哭怎么回事| 精品人妻1区二区| 亚洲国产精品成人综合色| 黄片播放在线免费| 国产不卡一卡二| 久久午夜亚洲精品久久| 日韩精品中文字幕看吧| 久久久精品欧美日韩精品| 丰满人妻熟妇乱又伦精品不卡| 国产欧美日韩一区二区三区在线| 丝袜人妻中文字幕| 色综合婷婷激情| 国产av一区二区精品久久| 国产国语露脸激情在线看| 日韩国内少妇激情av| 国产精品二区激情视频| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| 性色av乱码一区二区三区2| 两个人看的免费小视频| 欧美在线黄色| 免费在线观看亚洲国产| 丝袜在线中文字幕| 一级,二级,三级黄色视频| 真人做人爱边吃奶动态| 国产一卡二卡三卡精品| 成人三级做爰电影| 国产三级在线视频| 一本久久中文字幕| 国产乱人伦免费视频| 美女高潮到喷水免费观看| 免费久久久久久久精品成人欧美视频| 日韩精品青青久久久久久| 午夜福利在线观看吧| 亚洲av片天天在线观看| 亚洲成av片中文字幕在线观看| 十分钟在线观看高清视频www| 动漫黄色视频在线观看| 欧美日韩乱码在线| 男男h啪啪无遮挡| 美女大奶头视频| 两个人视频免费观看高清| 校园春色视频在线观看| 两个人视频免费观看高清| 精品午夜福利视频在线观看一区| 超碰成人久久| 色综合亚洲欧美另类图片| 成人手机av| 9191精品国产免费久久| 亚洲精品在线观看二区| av视频在线观看入口| 国产激情久久老熟女| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 精品久久蜜臀av无| 一二三四在线观看免费中文在| 欧美午夜高清在线| 精品无人区乱码1区二区| 一本大道久久a久久精品| 久久人人精品亚洲av| 色精品久久人妻99蜜桃| 夜夜躁狠狠躁天天躁| 成人av一区二区三区在线看| 51午夜福利影视在线观看| 97超级碰碰碰精品色视频在线观看| 一边摸一边抽搐一进一出视频| 日本精品一区二区三区蜜桃| 69精品国产乱码久久久| 亚洲午夜理论影院| 人人妻人人澡欧美一区二区 | 女性生殖器流出的白浆| 精品久久蜜臀av无| 免费在线观看视频国产中文字幕亚洲| 日韩精品免费视频一区二区三区| 免费看十八禁软件| 亚洲 欧美一区二区三区| av天堂在线播放| 两个人视频免费观看高清| √禁漫天堂资源中文www| 国产成人欧美| 嫩草影院精品99| 老司机午夜福利在线观看视频| 黄色成人免费大全| 中文字幕人成人乱码亚洲影| 国产成人影院久久av| 看免费av毛片| 91字幕亚洲| 精品久久久久久久毛片微露脸| 亚洲一区高清亚洲精品| 亚洲国产欧美网| 亚洲国产高清在线一区二区三 | 国产aⅴ精品一区二区三区波| 亚洲一区二区三区不卡视频| 91成人精品电影| 国产成人精品无人区| 久久久久国产一级毛片高清牌| 中文字幕精品免费在线观看视频| av视频免费观看在线观看| 欧美成狂野欧美在线观看| 成人国产综合亚洲| 亚洲av熟女| 精品卡一卡二卡四卡免费| 一二三四社区在线视频社区8| 成人18禁在线播放| 久久精品国产亚洲av香蕉五月| 很黄的视频免费| 亚洲成人久久性| 国产精品久久视频播放| 黄色成人免费大全| 久久精品91无色码中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美一区二区三区黑人| 中文字幕精品免费在线观看视频| 91av网站免费观看| 色播在线永久视频| 亚洲狠狠婷婷综合久久图片| 99国产精品99久久久久| 久久香蕉国产精品| 成人欧美大片| 国产一区二区激情短视频| xxx96com| 男女床上黄色一级片免费看| 精品一区二区三区av网在线观看| 男女床上黄色一级片免费看| 我的亚洲天堂| 99国产精品免费福利视频| 国产精品永久免费网站| 日韩欧美一区二区三区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 男女做爰动态图高潮gif福利片 | 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 最近最新中文字幕大全电影3 | 国产亚洲av高清不卡| 免费女性裸体啪啪无遮挡网站| 日日爽夜夜爽网站| 欧美激情 高清一区二区三区| 日本三级黄在线观看| 欧美在线黄色| 免费久久久久久久精品成人欧美视频| 在线观看免费视频网站a站| 中出人妻视频一区二区| 美女午夜性视频免费| 人人妻人人澡人人看| 午夜视频精品福利| 在线观看日韩欧美| 久久久久久免费高清国产稀缺| 91成年电影在线观看| 天堂影院成人在线观看| 欧美av亚洲av综合av国产av| 母亲3免费完整高清在线观看| 淫妇啪啪啪对白视频| 动漫黄色视频在线观看| cao死你这个sao货| 国产精品一区二区精品视频观看| 精品国产乱码久久久久久男人| 亚洲精品国产色婷婷电影| 国产一区在线观看成人免费| 色av中文字幕| 男人舔女人下体高潮全视频| 欧美另类亚洲清纯唯美| 日韩 欧美 亚洲 中文字幕| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 在线av久久热| 国产一区在线观看成人免费| 欧洲精品卡2卡3卡4卡5卡区| 一进一出好大好爽视频| 9色porny在线观看| 99国产精品免费福利视频| 亚洲在线自拍视频| 黄色视频不卡| 一级毛片精品| 亚洲无线在线观看| av电影中文网址| 国内毛片毛片毛片毛片毛片| 一本大道久久a久久精品| 香蕉丝袜av| 一本大道久久a久久精品| 制服丝袜大香蕉在线| 国产免费男女视频| 欧美成人性av电影在线观看| 亚洲午夜精品一区,二区,三区| 69av精品久久久久久| 亚洲最大成人中文| 99国产精品一区二区蜜桃av| 精品福利观看| 嫁个100分男人电影在线观看| 中文亚洲av片在线观看爽| 久久久久国内视频| 亚洲精品中文字幕在线视频| 女警被强在线播放| 亚洲伊人色综图| 99久久精品国产亚洲精品| 国产亚洲精品第一综合不卡| 亚洲免费av在线视频| 久久久国产欧美日韩av| 色av中文字幕| 国产欧美日韩一区二区三| 午夜久久久久精精品| 午夜免费观看网址| 免费在线观看日本一区| 亚洲av电影在线进入| 亚洲国产欧美网| 极品教师在线免费播放| 两性午夜刺激爽爽歪歪视频在线观看 | 免费在线观看黄色视频的| 男人舔女人的私密视频| 亚洲人成伊人成综合网2020| 美女大奶头视频| 国产亚洲欧美精品永久| 久久婷婷人人爽人人干人人爱 | 视频区欧美日本亚洲| 国产亚洲精品av在线| 欧美激情高清一区二区三区| 亚洲av五月六月丁香网| 男人的好看免费观看在线视频 | 国产精品99久久99久久久不卡| 亚洲五月天丁香| 午夜两性在线视频| 欧美日韩福利视频一区二区| 天天躁夜夜躁狠狠躁躁| 日本撒尿小便嘘嘘汇集6| 涩涩av久久男人的天堂| 成人三级做爰电影| 1024香蕉在线观看| 中文字幕av电影在线播放| 亚洲av日韩精品久久久久久密| 在线国产一区二区在线| 国产高清有码在线观看视频 | 久久久久久久久久久久大奶| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 午夜激情av网站| 久久久久国内视频| 国产高清视频在线播放一区| 亚洲成a人片在线一区二区| 搡老岳熟女国产| 日韩欧美在线二视频| 欧美激情 高清一区二区三区| 怎么达到女性高潮| 成在线人永久免费视频| 亚洲最大成人中文| 亚洲人成电影观看| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | svipshipincom国产片| 黄片小视频在线播放| 国产激情欧美一区二区| 91av网站免费观看| 在线观看66精品国产| 欧美乱妇无乱码| 好男人电影高清在线观看| 最新在线观看一区二区三区| 亚洲片人在线观看| 精品国产乱子伦一区二区三区| 女性生殖器流出的白浆| 亚洲最大成人中文| 亚洲五月婷婷丁香| 日韩欧美国产一区二区入口| 嫩草影视91久久| 18美女黄网站色大片免费观看| 精品国产超薄肉色丝袜足j| 国产极品粉嫩免费观看在线| 无遮挡黄片免费观看| 亚洲av电影不卡..在线观看| 人人妻人人爽人人添夜夜欢视频| 91精品国产国语对白视频| 制服丝袜大香蕉在线| 97人妻天天添夜夜摸| 久久性视频一级片| 国产av在哪里看| 精品欧美一区二区三区在线| 在线视频色国产色| 久久久久久亚洲精品国产蜜桃av| www.自偷自拍.com| 国产亚洲欧美精品永久| 男女床上黄色一级片免费看| 欧美日本视频| 麻豆av在线久日| 久久久久久人人人人人| 久久久久国产一级毛片高清牌| 免费av毛片视频| 亚洲欧美日韩高清在线视频| 黄色视频,在线免费观看| 免费无遮挡裸体视频| 青草久久国产| 国产三级在线视频| 正在播放国产对白刺激| 久久精品人人爽人人爽视色| 手机成人av网站| 国产色视频综合| 99国产极品粉嫩在线观看| 人妻久久中文字幕网| 亚洲男人天堂网一区| 亚洲av成人一区二区三| 亚洲精华国产精华精| 每晚都被弄得嗷嗷叫到高潮| 久久久久久亚洲精品国产蜜桃av| 一进一出好大好爽视频| 中文字幕av电影在线播放| 国产精品日韩av在线免费观看 | 欧美黑人精品巨大| 十分钟在线观看高清视频www| 日韩欧美国产一区二区入口| 欧美黑人欧美精品刺激| 制服诱惑二区| 久久香蕉精品热| 69av精品久久久久久| 国产精品乱码一区二三区的特点 | 露出奶头的视频| 最新美女视频免费是黄的| 国产在线观看jvid| 国产精品秋霞免费鲁丝片| 好男人电影高清在线观看| 欧美人与性动交α欧美精品济南到| 夜夜看夜夜爽夜夜摸| 久久精品91无色码中文字幕| 曰老女人黄片| 99国产精品一区二区蜜桃av| 欧美绝顶高潮抽搐喷水| 欧美日韩亚洲综合一区二区三区_| 看片在线看免费视频| 国产午夜精品久久久久久| 国产精品电影一区二区三区| 人人妻人人爽人人添夜夜欢视频| 91大片在线观看| 黄片小视频在线播放| 国内久久婷婷六月综合欲色啪| 久久国产乱子伦精品免费另类| 黄色片一级片一级黄色片| 亚洲 国产 在线| 人人妻人人澡人人看| 12—13女人毛片做爰片一| 一区二区三区精品91| 麻豆成人av在线观看| 一区二区三区精品91| 不卡一级毛片| 亚洲成国产人片在线观看| 岛国视频午夜一区免费看| 一本综合久久免费| 两性夫妻黄色片| 美女大奶头视频| 精品国产乱子伦一区二区三区| 久久人妻熟女aⅴ| 亚洲中文字幕一区二区三区有码在线看 | 亚洲最大成人中文| 久久久国产成人免费| 日本欧美视频一区| 亚洲 欧美一区二区三区| 国产成人精品无人区| 一区二区三区高清视频在线| 免费久久久久久久精品成人欧美视频| 成人免费观看视频高清| 在线十欧美十亚洲十日本专区| 人成视频在线观看免费观看| 精品不卡国产一区二区三区| 99在线视频只有这里精品首页| 国产一区在线观看成人免费| 亚洲精品av麻豆狂野| 久久久久久久久久久久大奶| 此物有八面人人有两片| 亚洲天堂国产精品一区在线| 精品国产乱码久久久久久男人| 在线观看www视频免费| 老司机午夜福利在线观看视频| 亚洲av成人av| 免费无遮挡裸体视频| 国产成人一区二区三区免费视频网站| 免费av毛片视频| 亚洲一区二区三区不卡视频| 久久午夜亚洲精品久久| 国产又色又爽无遮挡免费看| 一区二区三区国产精品乱码| 日韩高清综合在线| 亚洲精品一区av在线观看| 久久久久久久久免费视频了| 国产伦人伦偷精品视频| 少妇粗大呻吟视频| 成人18禁在线播放| 国内毛片毛片毛片毛片毛片| 欧美日韩黄片免| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片 | 亚洲av成人不卡在线观看播放网| 欧美大码av| 日本a在线网址| 亚洲一区二区三区色噜噜| 亚洲欧美精品综合久久99| 国产亚洲欧美在线一区二区| 国产亚洲av高清不卡| 狂野欧美激情性xxxx| 国产黄a三级三级三级人| 一级a爱视频在线免费观看| 久久久精品国产亚洲av高清涩受| 国产精品免费一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 一区在线观看完整版| 午夜a级毛片| 欧美日本视频| 国产欧美日韩一区二区精品| 国产av一区二区精品久久| 一本综合久久免费| 99香蕉大伊视频| 他把我摸到了高潮在线观看| 免费在线观看视频国产中文字幕亚洲| 精品熟女少妇八av免费久了| 亚洲视频免费观看视频| 亚洲va日本ⅴa欧美va伊人久久| 国产极品粉嫩免费观看在线| 国产91精品成人一区二区三区| 欧美国产精品va在线观看不卡| 99香蕉大伊视频| 久久精品国产清高在天天线| 欧美精品啪啪一区二区三区| 免费在线观看黄色视频的| 国产精品,欧美在线| 国产精品乱码一区二三区的特点 | 精品一品国产午夜福利视频| 韩国av一区二区三区四区|