• <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頻段的覆蓋能力差異
    久久精品夜色国产| 国产欧美亚洲国产| 永久免费av网站大全| 精品一区二区三区四区五区乱码 | 国产深夜福利视频在线观看| 黄色配什么色好看| 日本色播在线视频| 久久久久久久精品精品| 99久国产av精品国产电影| 国产精品人妻久久久影院| 国产欧美日韩综合在线一区二区| 国产精品 国内视频| 欧美激情极品国产一区二区三区 | 成人影院久久| 亚洲欧洲日产国产| 51国产日韩欧美| 久久国产精品男人的天堂亚洲 | 狂野欧美激情性bbbbbb| 精品一区二区三卡| 欧美xxⅹ黑人| 午夜久久久在线观看| 老司机亚洲免费影院| 晚上一个人看的免费电影| 国产又色又爽无遮挡免| 亚洲av中文av极速乱| 亚洲内射少妇av| 老司机亚洲免费影院| 国产色爽女视频免费观看| 看十八女毛片水多多多| 美女主播在线视频| 中文字幕av电影在线播放| 蜜桃国产av成人99| 国产在线免费精品| 日韩中字成人| 亚洲精品久久午夜乱码| 999精品在线视频| a级片在线免费高清观看视频| a级毛片在线看网站| 亚洲av福利一区| 国内精品宾馆在线| 精品一区二区三卡| 男女啪啪激烈高潮av片| 91久久精品国产一区二区三区| 最近手机中文字幕大全| 亚洲欧美一区二区三区国产| 亚洲熟女精品中文字幕| 亚洲性久久影院| 成年女人在线观看亚洲视频| 边亲边吃奶的免费视频| 国产亚洲最大av| 久久人人爽人人片av| 欧美变态另类bdsm刘玥| 2022亚洲国产成人精品| 色视频在线一区二区三区| 另类精品久久| 国产白丝娇喘喷水9色精品| 欧美人与性动交α欧美精品济南到 | 日本av手机在线免费观看| 考比视频在线观看| 午夜精品国产一区二区电影| 免费黄网站久久成人精品| 国产欧美日韩一区二区三区在线| 丝袜在线中文字幕| 久久人人爽人人片av| 国产 精品1| 国产成人aa在线观看| av在线老鸭窝| 婷婷色麻豆天堂久久| av一本久久久久| 99热全是精品| 久久久久国产网址| av免费观看日本| 一区二区三区四区激情视频| 男人舔女人的私密视频| 全区人妻精品视频| 亚洲欧洲日产国产| 欧美3d第一页| 免费高清在线观看日韩| 777米奇影视久久| 五月天丁香电影| 男人添女人高潮全过程视频| 蜜臀久久99精品久久宅男| 国产精品一区二区在线不卡| 女的被弄到高潮叫床怎么办| 老司机影院毛片| 色婷婷av一区二区三区视频| 9热在线视频观看99| 久久久精品免费免费高清| 亚洲精品乱码久久久久久按摩| 欧美少妇被猛烈插入视频| 最近2019中文字幕mv第一页| 免费在线观看完整版高清| 国产一区二区三区av在线| av在线播放精品| 午夜精品国产一区二区电影| 老司机影院毛片| 中国三级夫妇交换| 亚洲精品第二区| 欧美激情 高清一区二区三区| 欧美精品一区二区免费开放| 丰满少妇做爰视频| 免费少妇av软件| 男的添女的下面高潮视频| 91久久精品国产一区二区三区| 在线免费观看不下载黄p国产| 国产精品人妻久久久久久| 街头女战士在线观看网站| 久久久久久久精品精品| 亚洲伊人色综图| 香蕉丝袜av| 日韩大片免费观看网站| 亚洲精品日本国产第一区| 久久久精品免费免费高清| 亚洲激情五月婷婷啪啪| 亚洲色图综合在线观看| 午夜福利视频在线观看免费| 美女大奶头黄色视频| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三区在线| 精品国产露脸久久av麻豆| xxx大片免费视频| 亚洲人与动物交配视频| www.av在线官网国产| 国产成人91sexporn| 老司机影院成人| 国产成人午夜福利电影在线观看| 日韩精品有码人妻一区| 老司机亚洲免费影院| 少妇的逼水好多| 这个男人来自地球电影免费观看 | 国产精品三级大全| 久久久精品94久久精品| 色吧在线观看| 一级毛片电影观看| 男男h啪啪无遮挡| 热re99久久精品国产66热6| 亚洲国产毛片av蜜桃av| 夫妻性生交免费视频一级片| 91成人精品电影| 秋霞伦理黄片| 纵有疾风起免费观看全集完整版| 日本vs欧美在线观看视频| 中文字幕av电影在线播放| 韩国高清视频一区二区三区| 亚洲欧美成人精品一区二区| 亚洲国产av新网站| 最近的中文字幕免费完整| 久久精品国产亚洲av涩爱| 中文欧美无线码| 免费久久久久久久精品成人欧美视频 | 一级a做视频免费观看| 久久韩国三级中文字幕| a 毛片基地| 极品少妇高潮喷水抽搐| 波野结衣二区三区在线| 岛国毛片在线播放| 国产欧美日韩一区二区三区在线| www.av在线官网国产| 欧美精品人与动牲交sv欧美| √禁漫天堂资源中文www| 在现免费观看毛片| 飞空精品影院首页| 中国美白少妇内射xxxbb| 国产白丝娇喘喷水9色精品| 五月伊人婷婷丁香| 中国国产av一级| 婷婷色综合www| 国产日韩欧美亚洲二区| 亚洲经典国产精华液单| 少妇被粗大猛烈的视频| 国产综合精华液| 女性被躁到高潮视频| 国产又色又爽无遮挡免| 草草在线视频免费看| 99久国产av精品国产电影| 18禁国产床啪视频网站| 成人国语在线视频| 国产成人欧美| 国产成人精品久久久久久| 久久影院123| 黑人巨大精品欧美一区二区蜜桃 | 一二三四中文在线观看免费高清| 91成人精品电影| 大香蕉久久网| 免费人妻精品一区二区三区视频| 成人国语在线视频| 天堂8中文在线网| 色婷婷久久久亚洲欧美| 日韩一区二区三区影片| videos熟女内射| 国产欧美亚洲国产| 亚洲熟女精品中文字幕| 日本爱情动作片www.在线观看| 啦啦啦啦在线视频资源| 国产免费一级a男人的天堂| 综合色丁香网| 在线观看免费视频网站a站| 亚洲天堂av无毛| av不卡在线播放| kizo精华| 国产成人精品婷婷| 亚洲图色成人| 国产精品无大码| 一本大道久久a久久精品| 国产成人午夜福利电影在线观看| 一区二区av电影网| 亚洲精品成人av观看孕妇| 久久久国产精品麻豆| 日本猛色少妇xxxxx猛交久久| 欧美成人精品欧美一级黄| 国产免费福利视频在线观看| 久久 成人 亚洲| 成年女人在线观看亚洲视频| 久久久久久伊人网av| 秋霞在线观看毛片| 日本-黄色视频高清免费观看| 欧美性感艳星| 在线观看人妻少妇| 久久精品久久久久久噜噜老黄| 亚洲国产色片| 国产av精品麻豆| 搡老乐熟女国产| 久久精品国产亚洲av天美| 亚洲精品乱码久久久久久按摩| 99久久综合免费| 午夜免费男女啪啪视频观看| 亚洲国产欧美在线一区| av线在线观看网站| av免费观看日本| a级毛片黄视频| 精品酒店卫生间| 国产精品欧美亚洲77777| 内地一区二区视频在线| av播播在线观看一区| 婷婷色综合www| 欧美人与善性xxx| 大码成人一级视频| 免费黄网站久久成人精品| 日韩av不卡免费在线播放| 亚洲高清免费不卡视频| 久久久久久久久久成人| 91久久精品国产一区二区三区| 日日啪夜夜爽| 国产欧美日韩综合在线一区二区| 国产精品国产三级国产专区5o| 欧美日韩一区二区视频在线观看视频在线| 看免费av毛片| 亚洲欧美成人精品一区二区| 亚洲内射少妇av| 男人舔女人的私密视频| 只有这里有精品99| 亚洲人与动物交配视频| 国产熟女欧美一区二区| 91精品伊人久久大香线蕉| 国产又色又爽无遮挡免| 亚洲精品乱久久久久久| 老熟女久久久| 国产一区二区在线观看日韩| 欧美3d第一页| 亚洲精品久久久久久婷婷小说| 丰满少妇做爰视频| 国产白丝娇喘喷水9色精品| 一边亲一边摸免费视频| 女性被躁到高潮视频| 国产日韩欧美亚洲二区| 国产探花极品一区二区| 成人手机av| 国产亚洲最大av| 韩国高清视频一区二区三区| 不卡视频在线观看欧美| 亚洲av成人精品一二三区| 日韩欧美精品免费久久| 青春草视频在线免费观看| 精品国产乱码久久久久久小说| 久久久久精品性色| 寂寞人妻少妇视频99o| 国产精品久久久久久av不卡| 亚洲性久久影院| 一边亲一边摸免费视频| 女的被弄到高潮叫床怎么办| 亚洲精品日韩在线中文字幕| 97人妻天天添夜夜摸| 久久国产精品大桥未久av| 观看美女的网站| 中文字幕制服av| 欧美日韩综合久久久久久| 日韩大片免费观看网站| a 毛片基地| 熟女av电影| 久久精品久久久久久噜噜老黄| 大香蕉97超碰在线| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| av国产精品久久久久影院| 伊人久久国产一区二区| 纵有疾风起免费观看全集完整版| 视频中文字幕在线观看| 日韩欧美一区视频在线观看| 如日韩欧美国产精品一区二区三区| 中文字幕av电影在线播放| 久久精品久久久久久久性| videosex国产| 日韩三级伦理在线观看| 黄网站色视频无遮挡免费观看| 国产色爽女视频免费观看| 人成视频在线观看免费观看| 久久综合国产亚洲精品| 精品99又大又爽又粗少妇毛片| 国产精品久久久久久久电影| 欧美国产精品一级二级三级| 国产成人欧美| 日本-黄色视频高清免费观看| 亚洲欧洲精品一区二区精品久久久 | 一级片免费观看大全| 精品国产一区二区三区久久久樱花| 99香蕉大伊视频| 国语对白做爰xxxⅹ性视频网站| 美女中出高潮动态图| 91午夜精品亚洲一区二区三区| av国产精品久久久久影院| 日本91视频免费播放| 亚洲精品日韩在线中文字幕| 精品少妇内射三级| h视频一区二区三区| 一本色道久久久久久精品综合| 中文字幕亚洲精品专区| 最新中文字幕久久久久| 如何舔出高潮| 最新中文字幕久久久久| 欧美成人精品欧美一级黄| 中文欧美无线码| 麻豆精品久久久久久蜜桃| 久久久亚洲精品成人影院| 免费观看无遮挡的男女| 亚洲国产精品国产精品| 日韩一本色道免费dvd| 国产亚洲午夜精品一区二区久久| 三级国产精品片| 国产探花极品一区二区| 一区二区av电影网| 国产亚洲精品久久久com| 国产精品人妻久久久久久| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 捣出白浆h1v1| 一级片'在线观看视频| 国产深夜福利视频在线观看| 精品一区二区免费观看| 熟妇人妻不卡中文字幕| 极品人妻少妇av视频| 丝瓜视频免费看黄片| 欧美bdsm另类| 亚洲精品一区蜜桃| 色吧在线观看| 国产xxxxx性猛交| av又黄又爽大尺度在线免费看| 一个人免费看片子| 啦啦啦在线观看免费高清www| 天天影视国产精品| 免费久久久久久久精品成人欧美视频 | 精品一区二区三卡| √禁漫天堂资源中文www| 在线观看免费视频网站a站| 国产男女内射视频| 午夜免费男女啪啪视频观看| 久久人人爽av亚洲精品天堂| 日韩,欧美,国产一区二区三区| 天美传媒精品一区二区| 少妇高潮的动态图| 成人无遮挡网站| 99国产综合亚洲精品| 久久精品国产综合久久久 | 日韩 亚洲 欧美在线| 22中文网久久字幕| 麻豆乱淫一区二区| 波多野结衣一区麻豆| 黄色 视频免费看| 国产熟女午夜一区二区三区| 自线自在国产av| 日本wwww免费看| 日日啪夜夜爽| 777米奇影视久久| 亚洲国产看品久久| 激情五月婷婷亚洲| 亚洲成av片中文字幕在线观看 | 又黄又粗又硬又大视频| 成人毛片60女人毛片免费| 下体分泌物呈黄色| 日产精品乱码卡一卡2卡三| 建设人人有责人人尽责人人享有的| 日韩av不卡免费在线播放| 亚洲精品第二区| 久久99精品国语久久久| 国产精品蜜桃在线观看| 亚洲伊人久久精品综合| 巨乳人妻的诱惑在线观看| 国产亚洲欧美精品永久| 80岁老熟妇乱子伦牲交| 天堂8中文在线网| 日韩,欧美,国产一区二区三区| 啦啦啦视频在线资源免费观看| 久久久久久久久久久免费av| 咕卡用的链子| 国产女主播在线喷水免费视频网站| 黄色视频在线播放观看不卡| 最新的欧美精品一区二区| 天堂俺去俺来也www色官网| 制服诱惑二区| 久久精品国产a三级三级三级| 女人精品久久久久毛片| 日韩电影二区| 成人漫画全彩无遮挡| 久久久久精品人妻al黑| videossex国产| 免费观看在线日韩| 少妇 在线观看| 国产精品成人在线| 高清毛片免费看| 十分钟在线观看高清视频www| 午夜影院在线不卡| 亚洲精品国产色婷婷电影| 欧美老熟妇乱子伦牲交| 久久 成人 亚洲| 大陆偷拍与自拍| 热99国产精品久久久久久7| 男的添女的下面高潮视频| 国产日韩欧美视频二区| 日韩伦理黄色片| av黄色大香蕉| 欧美bdsm另类| 一二三四在线观看免费中文在 | av在线观看视频网站免费| 亚洲国产最新在线播放| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看| 高清毛片免费看| 69精品国产乱码久久久| 一边摸一边做爽爽视频免费| 亚洲精品aⅴ在线观看| 中文字幕精品免费在线观看视频 | 精品国产国语对白av| 在线天堂最新版资源| 曰老女人黄片| 国产一区二区三区av在线| 免费黄网站久久成人精品| 国产精品国产三级国产av玫瑰| 精品一品国产午夜福利视频| 色视频在线一区二区三区| 久久精品人人爽人人爽视色| www日本在线高清视频| 午夜日本视频在线| 如何舔出高潮| 18禁裸乳无遮挡动漫免费视频| 午夜福利乱码中文字幕| 国产精品嫩草影院av在线观看| 日韩电影二区| 少妇的丰满在线观看| 高清在线视频一区二区三区| 免费在线观看黄色视频的| 亚洲成人一二三区av| 亚洲国产欧美日韩在线播放| 高清毛片免费看| 日本与韩国留学比较| 久久久久国产精品人妻一区二区| 99国产综合亚洲精品| 日日爽夜夜爽网站| 不卡视频在线观看欧美| 国产精品一区二区在线不卡| 男女无遮挡免费网站观看| 亚洲精品色激情综合| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 久久韩国三级中文字幕| 丰满饥渴人妻一区二区三| 制服丝袜香蕉在线| 日本欧美国产在线视频| 九九在线视频观看精品| 91精品国产国语对白视频| 蜜桃在线观看..| 久久久国产一区二区| 天天影视国产精品| 哪个播放器可以免费观看大片| av国产久精品久网站免费入址| 亚洲久久久国产精品| 人人澡人人妻人| 欧美激情 高清一区二区三区| 制服人妻中文乱码| 成人免费观看视频高清| 一二三四中文在线观看免费高清| 伦理电影大哥的女人| 国产亚洲午夜精品一区二区久久| 精品一区在线观看国产| 丁香六月天网| 亚洲欧美成人综合另类久久久| 建设人人有责人人尽责人人享有的| 满18在线观看网站| 99久久精品国产国产毛片| 亚洲成人av在线免费| 美女福利国产在线| 国产在线视频一区二区| 欧美国产精品一级二级三级| h视频一区二区三区| 有码 亚洲区| 国产欧美另类精品又又久久亚洲欧美| 99热6这里只有精品| 制服丝袜香蕉在线| 大片免费播放器 马上看| 国产白丝娇喘喷水9色精品| 狂野欧美激情性xxxx在线观看| 国产亚洲av片在线观看秒播厂| 亚洲美女视频黄频| 18禁观看日本| 国产精品无大码| 一边摸一边做爽爽视频免费| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影小说| 永久免费av网站大全| 国产精品久久久久久精品电影小说| 男女免费视频国产| 日本欧美国产在线视频| 成年美女黄网站色视频大全免费| 亚洲av免费高清在线观看| 精品国产露脸久久av麻豆| 一级毛片我不卡| 中文字幕人妻丝袜制服| 午夜老司机福利剧场| 亚洲人成77777在线视频| videossex国产| 免费大片黄手机在线观看| 久久精品久久精品一区二区三区| 精品人妻熟女毛片av久久网站| 在线观看美女被高潮喷水网站| 免费在线观看完整版高清| 亚洲成人手机| 国产黄色免费在线视频| 亚洲少妇的诱惑av| 国产日韩一区二区三区精品不卡| 香蕉丝袜av| 欧美丝袜亚洲另类| 国产精品人妻久久久影院| 久久精品国产亚洲av涩爱| 免费日韩欧美在线观看| 最黄视频免费看| 成人综合一区亚洲| 日日爽夜夜爽网站| 久久久亚洲精品成人影院| 乱人伦中国视频| 国产男女内射视频| 97超碰精品成人国产| 黄色怎么调成土黄色| 在线观看免费高清a一片| 亚洲成色77777| 丰满少妇做爰视频| 夫妻性生交免费视频一级片| av电影中文网址| 肉色欧美久久久久久久蜜桃| 国产在视频线精品| 日本黄大片高清| 99久国产av精品国产电影| 午夜福利网站1000一区二区三区| av国产久精品久网站免费入址| 成人亚洲欧美一区二区av| 精品久久久精品久久久| 日韩av不卡免费在线播放| 伊人久久国产一区二区| 久久热在线av| 满18在线观看网站| 九九爱精品视频在线观看| 成人影院久久| 免费观看在线日韩| 纯流量卡能插随身wifi吗| 一级a做视频免费观看| 欧美最新免费一区二区三区| 韩国高清视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美亚洲国产| 少妇 在线观看| 中国三级夫妇交换| 久久久久国产网址| 亚洲精品第二区| 18禁裸乳无遮挡动漫免费视频| 日本黄色日本黄色录像| 国产黄色视频一区二区在线观看| 桃花免费在线播放| 国产一区二区三区综合在线观看 | av.在线天堂| 久久久久久久大尺度免费视频| 丰满乱子伦码专区| 中文字幕另类日韩欧美亚洲嫩草| 香蕉精品网在线| www.熟女人妻精品国产 | 三级国产精品片| 五月玫瑰六月丁香| 久久精品人人爽人人爽视色| 久久99热6这里只有精品| 午夜激情久久久久久久| 波野结衣二区三区在线| 欧美3d第一页| 男女边吃奶边做爰视频| 国产有黄有色有爽视频| 久久久久久久大尺度免费视频| 国产亚洲最大av| 亚洲欧美一区二区三区黑人 | 欧美日韩综合久久久久久| 蜜桃国产av成人99| 一级毛片电影观看| 日本-黄色视频高清免费观看| 免费黄网站久久成人精品| 国产精品女同一区二区软件| 边亲边吃奶的免费视频| 新久久久久国产一级毛片| 亚洲色图综合在线观看| 18在线观看网站| 精品国产一区二区三区四区第35| 精品国产国语对白av|