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

    基于似然函數(shù)的雙曲調(diào)頻信號參數(shù)估計快速算法

    2021-05-30 07:22:34馬碧云元達鵬劉嬌蛟
    電子與信息學(xué)報 2021年5期
    關(guān)鍵詞:下界均方參數(shù)估計

    馬碧云 元達鵬 劉嬌蛟*

    ①(華南理工大學(xué)電子與信息學(xué)院 廣州510640)

    ②(自然資源部海洋環(huán)境探測技術(shù)與應(yīng)用重點實驗室 廣州510300)

    1 引言

    線性調(diào)頻(Liner Frequency Modulation,LFM)信號被廣泛地應(yīng)用于雷達偵查、地震勘測、水聲通信與仿生聲吶等領(lǐng)域,而能否準(zhǔn)確地估計出接收信號的參數(shù),將直接影響探測精度,所以對LFM信號的參數(shù)估計是一個近年來比較重要的研究方向,國內(nèi)外學(xué)者也作了深入廣泛的研究[1–5]。然而在水聲環(huán)境等復(fù)雜的應(yīng)用場景中,物體間的相對運動會造成多普勒效應(yīng),這種多普勒效應(yīng)不僅會影響接收信號的脈沖壓縮結(jié)果,也會影響接收信號的參數(shù)估計結(jié)果,因此LFM信號難以適應(yīng)此類動態(tài)場景的需求。有鑒于此,研究者引入了雙曲調(diào)頻(Hyperbolic Frequency Modulation,HFM)信號。HFM信號不僅具有良好的脈沖壓縮性能,還具有多普勒不變性[6],即經(jīng)過多普勒效應(yīng)影響的接收信號仍可與匹配濾波器匹配,因此HFM信號近年來在水聲等多普勒影響嚴(yán)重的應(yīng)用領(lǐng)域中得到了廣泛的研究,但是目前針對HFM信號參數(shù)估計問題的研究較為有限[7,8]。在參數(shù)估計問題上,最大似然估計由于其自身估計算法的精度極高,受到廣大研究學(xué)者的推崇,然而最大似然估計涉及多維搜索的過程,存在因計算量大而無法直接在實際工程中使用的問題。因此本文針對HFM信號,提出一種基于似然函數(shù)的HFM信號參數(shù)估計快速方法。

    本文安排如下:第1節(jié)敘述HFM信號以及其參數(shù)估計問題的研究背景以及應(yīng)用意義;第2,3節(jié)將介紹HFM信號的特征,并推導(dǎo)HFM信號參數(shù)估計的Cramer-Rao下界,作為本文參數(shù)估計方法的性能評估標(biāo)準(zhǔn);第4,5節(jié)將基于噪聲的高斯隨機過程,構(gòu)建HFM信號的似然函數(shù),利用GABC算法對似然函數(shù)進行全局最優(yōu)化搜索,并針對算法的優(yōu)化過程,提出一種改進的適應(yīng)度函數(shù);第6節(jié)將通過蒙特卡洛試驗驗證本文所述方法的有效性,結(jié)果表明本方法在保證估計精度的同時能提高算法收斂速度。

    2 HFM信號

    可見式(2)是關(guān)于時間t的雙曲函數(shù),因此稱為雙曲調(diào)頻信號。

    另外,HFM信號的多普勒不變性表現(xiàn)為:經(jīng)過多普勒效應(yīng)影響的HFM信號除了存在一個時域上的時間偏移 ?t外,并不會產(chǎn)生多普勒頻移,且其瞬時頻率仍然保持相同的調(diào)頻率k與初始頻率f0,因此仍可與匹配濾波器完全匹配[9],上述結(jié)果可表示為

    式中,fH(t)為 原HFM信號的瞬時頻率,fDH(t)為HFM信號經(jīng)過多普勒效應(yīng)影響后的瞬時頻率,時間偏移?t=(1?α)/αkf0,其中α為多普勒尺度因子。

    因此相比于LFM信號,HFM信號具有多普勒不變性,經(jīng)多普勒尺度變化后的信號參數(shù)仍保持不變,在復(fù)雜水聲環(huán)境中的信號參數(shù)估計問題上具有一定優(yōu)勢。

    3 高斯白噪聲下的Cramer-Rao下界

    在參數(shù)估計問題中,Cramer-Rao下界為任意無偏估計量的方差確定下界,即無偏估計量的方差只能無限逼近Cramer-Rao下界,但不會小于該下界[10]。因此,本節(jié)通過推導(dǎo)Cramer-Rao下界,為任意無偏估計量的性能比較提供了一個參考標(biāo)準(zhǔn)。

    在高斯白噪聲環(huán)境下,采樣所得的HFM信號服從高斯隨機過程

    4 HFM信號的似然函數(shù)

    HFM信號參數(shù)的最大似然估計為漸近無偏估計,在理論上最大似然估計可以逼近Cramer-Rao下界,具有極高的估計精度。

    其中,fs為采樣頻率;w(n),n∈{1 2···M}為零均值的復(fù)高斯白噪聲隨機過程,其實部與虛部方差均為常量σ2,且相互獨立;其中調(diào)頻率k與初始頻率f0的最大似然估計值為HFM信號參數(shù)估計的最終目的。

    為了方便化簡,式(12)可表示為向量形式

    5 似然函數(shù)優(yōu)化問題

    傳統(tǒng)的最優(yōu)化算法如牛頓迭代法、最速下降法、共軛梯度法等,由于對目標(biāo)函數(shù)有著較高要求,比如要求目標(biāo)函數(shù)必須可導(dǎo),因此在大多數(shù)實際問題中無法使用。近年來,仿生智能優(yōu)化算法因其對復(fù)雜問題具有高效的優(yōu)化性能且無需目標(biāo)函數(shù)導(dǎo)數(shù)信息等優(yōu)點,受到研究學(xué)者的廣泛關(guān)注。在最優(yōu)化算法的搜索過程中,均需要通過反復(fù)計算適應(yīng)度函數(shù)值進行個體篩選,因此適應(yīng)度函數(shù)的選取將直接決定算法的收斂速度以及能否尋得最優(yōu)解。有鑒于此,本文根據(jù)最優(yōu)化算法的搜索過程改進了其中的適應(yīng)度函數(shù),使得最優(yōu)化算法在保證準(zhǔn)確性的同時降低一定的計算量。

    5.1 最優(yōu)化算法

    隨著1975年遺傳算法(Genetic Algorithms,GA)的提出,經(jīng)過幾十年的發(fā)展,研究學(xué)者針對不同的應(yīng)用問題提出了多種仿生智能算法,各種方法因適用范圍不同,都有一定的優(yōu)缺點,需要根據(jù)實際的問題進行選擇應(yīng)用。

    Karaboga等人[12]于2007年提出人工蜂群(Artificial Bee Colony,ABC)算法用于解決多變量函數(shù)優(yōu)化問題,后來Zhu等人[13]提出全局最優(yōu)引導(dǎo)的人工蜂群算法(Global best-guided Artificial Bee Colony algorithm,GABC),它的基本思想是將每個蜜源位置表示為待求問題的一個可行解,并通過模仿蜜蜂采蜜行為來進行最優(yōu)解搜索。GABC算法在多維函數(shù)優(yōu)化問題上表現(xiàn)優(yōu)秀,且具有魯棒性強、控制參數(shù)較少等優(yōu)點。相比于ABC算法,GABC算法通過跟隨當(dāng)前最優(yōu)值進行搜索,優(yōu)化結(jié)果精度更高。因此針對似然函數(shù)的優(yōu)化問題,本文選擇GABC算法作為下文研究HFM信號參數(shù)估計問題的最優(yōu)化算法。

    5.2 時間復(fù)雜度分析

    在優(yōu)化過程中,最優(yōu)化算法無一例外地需要在每一次迭代中,對適應(yīng)度函數(shù)進行與個體數(shù)量同等次數(shù)的計算,以得到所在位置參數(shù)下的適應(yīng)度大小,因此不同最優(yōu)化算法的時間復(fù)雜度如表1所示。

    其中最大似然估計中的D為待估計參數(shù)的總數(shù),u j與lj分別是第j個參數(shù)的搜索范圍上界與下界,p j為第j個 參數(shù)的遍歷精度;CGA,CABC,CGABC分別為GA算法、ABC算法、GABC算法的個體數(shù)量,NGA,NABC,NGABC分別為GA算法,ABC算法,GABC算法的最大迭代次數(shù),Tfit為適應(yīng)度函數(shù)的時間復(fù)雜度。

    由表1可以看出,適應(yīng)度函數(shù)的復(fù)雜度是最優(yōu)化算法整體復(fù)雜度的主要組成部分,因此適應(yīng)度函數(shù)的設(shè)計需要盡可能簡單,使得整體算法的時間復(fù)雜度最小。

    5.3 適應(yīng)度函數(shù)的優(yōu)化

    在最優(yōu)化算法中,優(yōu)化過程正是對解成員進行“優(yōu)勝劣汰”的過程,而解成員的好壞由它的適應(yīng)度函數(shù)值所決定,因此適應(yīng)度函數(shù)的選取將直接決定算法的收斂速度以及能否尋得最優(yōu)解。

    表1 不同方法的時間復(fù)雜度

    對于HFM信號的參數(shù)估計問題,傳統(tǒng)的適應(yīng)度函數(shù)為壓縮概率函數(shù)[14],即蜜源個體mi的適應(yīng)度fiti(m i)為

    然而若采用式(24)作為適應(yīng)度函數(shù),每次計算都將涉及矩陣取逆的高復(fù)雜度操作;并且由于參數(shù)(k i,f0,i)在優(yōu)化過程中是隨機選取,因此存在SH(k i,f0,i)S(k i,f0,i)矩陣不可逆的可能性;此外在部分最優(yōu)化算法(如ABC算法、GABC算法)中,涉及同時對多個個體進行適應(yīng)度計算,而式(24)的每次計算只能得到單個個體的適應(yīng)度,缺乏并行計算能力。

    使用式(25)作為最優(yōu)化算法的適應(yīng)度函數(shù)計算公式,不僅避免了式(24)復(fù)雜的矩陣求逆運算,以及存在矩陣不可逆而導(dǎo)致的精度下降,且避免通過低效率的循環(huán)去計算多個個體[m1m2···m C]的適應(yīng)度值,從而提升了最優(yōu)算法的整體收斂速度。

    6 實驗驗證

    實驗1:適應(yīng)度函數(shù)的有效性驗證

    由于最優(yōu)化算法是通過不斷比較適應(yīng)度函數(shù)值的大小進行優(yōu)化搜索,因此適應(yīng)度函數(shù)的有效性需要首先得到驗證,即通過對適應(yīng)度函數(shù)進行優(yōu)化搜索能得到正確的參數(shù)估計值。圖1分別為適應(yīng)度函數(shù)式(24)與優(yōu)化后的適應(yīng)度函數(shù)式(25)在不同參數(shù)組合下的2維遍歷,其中HFM信號的初始頻率f0為1.5 MHz,調(diào)頻率k為–0.077802。

    如圖1(a)所示,當(dāng)適應(yīng)度函數(shù)取得最大值時,橫縱坐標(biāo)所代表的參數(shù)估計值與初始設(shè)定值相同。其中縱軸為待估計HFM信號的初始頻率,橫軸表示待估計HFM信號的調(diào)頻率。而相比于圖1(a),圖1(b)雖在部分有一定的隆起,但在該適應(yīng)度函數(shù)取得最大值時,橫縱坐標(biāo)所代表的參數(shù)估計值也與初始設(shè)定值相同,滿足作為適應(yīng)度函數(shù)的條件。

    實驗2:估計值的均方誤差與Cramer-Rao下界對比

    本文通過1000次蒙特卡洛仿真試驗,對比了多種參數(shù)估計方法在HFM信號參數(shù)估計問題中的性能。在優(yōu)化算法的參數(shù)設(shè)置如下:在GA算法中,令種群大小為40,最大遺傳代數(shù)為1000,個體長度為20,代溝為0.95,交叉概率為0.7,變異概率為0.01;在ABC算法中,蜂群群落大小為60,最大迭代次數(shù)為1000;在GABC算法中,蜂群群落大小為80,最大迭代次數(shù)為100;在最大似然估計法中,f0的遍歷精度為50 Hz;k 的遍歷精度為 1 ×10?5。

    圖1 不同適應(yīng)度函數(shù)的2維遍歷結(jié)果

    圖2 不同方法下估計參數(shù)的均方誤差與Cramer-Rao下界

    圖2為在–10~10 d B信噪比內(nèi)多種參數(shù)估計方法的性能比較。如圖2所示,5種估計方法的均方誤差都隨著信噪比的增加,總體呈現(xiàn)下降的趨勢:其中最大似然估計的性能最優(yōu),均方誤差可隨著搜索精度的提升而逼近Cramer-Rao下界;對于GA算法,在信噪比達到0 d B以上時,受信噪比的影響將減小,均方誤差的變化將趨于平緩;與ABC算法而言,由于問題針對性較強,因此其估計結(jié)果的均方誤差小于GA算法,但仍與最大似然估計的性能有明顯差距;相比之下,GABC算法由于引入了全局最優(yōu)引導(dǎo)項,因此其估計結(jié)果的精度更高,均方誤差能隨著信噪比的增加而逼近Cramer-Rao下界。

    另外從圖中可知,標(biāo)準(zhǔn)GABC算法相比于采用優(yōu)化后的適應(yīng)度函數(shù)式(25)的GABC算法,前者存在矩陣不可逆的可能性,因此兩曲線在圖中沒有完全重合,但總體性能相差不大;在信噪比為4 dB以上時,后者的均方誤差更接近于最大似然估計的性能。而在信噪比為10 d B時,對于參數(shù)k和f0,與最大似然估計的均方誤差相比,前者相差分別為4.44 d B與3.22 d B,而后者分別相差2.09 d B與0.68 d B,后者更接近最大似然估計的性能。在水下應(yīng)用中由于噪聲水平不高,信噪比大致保持在5 dB以上,因此文中所述方法能被有效地用于HFM信號的參數(shù)估計問題中。

    實驗3:收斂速度對比

    在完成有效性的驗證后,下面對比GABC算法、與采用優(yōu)化后的適應(yīng)度函數(shù)式(25)的GABC算法在HFM信號參數(shù)估計問題中的收斂速度以及算法耗時,其中兩種GABC算法的參數(shù)設(shè)置均相同,蜂群群落大小為20,最大迭代次數(shù)為80。

    如圖3所示,GABC算法和采用優(yōu)化后的適應(yīng)度函數(shù)式(25)的GABC算法均能在50次迭代內(nèi)尋得最優(yōu)值,雖然后者收斂速度稍快,但總體相差不大。因此可以預(yù)見的是,在相同的迭代次數(shù)內(nèi),適應(yīng)度函數(shù)的時間復(fù)雜度越低,GABC算法的實現(xiàn)效率就越高。

    圖3 基于不同適應(yīng)度函數(shù)的GABC算法收斂速度

    因此本文通過Matlab中的探查器(Profile)可得算法的時間消耗,即最大似然估計的時間消耗為13738.8 s,GABC算法的時間消耗為1.785 s,適應(yīng)度函數(shù)優(yōu)化后的GABC算法的時間消耗為0.715 s。可見針對HFM信號的參數(shù)估計問題,采用本文提出的適應(yīng)度函數(shù)式(25)的GABC算法在具有較高精準(zhǔn)度的同時,也具有較快的收斂速度。

    7 總結(jié)

    本文首先推導(dǎo)出HFM信號參數(shù)估計的Cramer-Rao下界,作為本文參數(shù)估計性能的評估標(biāo)準(zhǔn)。然后基于高斯隨機噪聲,構(gòu)建了HFM信號的似然函數(shù),并結(jié)合數(shù)據(jù)向量化的特點提出一種改進的適應(yīng)度函數(shù),最后利用GABC算法對該適應(yīng)度函數(shù)進行極值尋優(yōu),從而實現(xiàn)HFM信號的參數(shù)估計;通過蒙特卡洛仿真證明了本方法在信噪比為3 d B以上時,HFM信號的參數(shù)估計結(jié)果的均方誤差更逼近Cramer-Rao下界,且運算量約是原來的1/3。

    猜你喜歡
    下界均方參數(shù)估計
    一類隨機積分微分方程的均方漸近概周期解
    基于新型DFrFT的LFM信號參數(shù)估計算法
    Beidou, le système de navigation par satellite compatible et interopérable
    Lower bound estimation of the maximum allowable initial error and its numerical calculation
    Logistic回歸模型的幾乎無偏兩參數(shù)估計
    基于向前方程的平穩(wěn)分布參數(shù)估計
    基于競爭失效數(shù)據(jù)的Lindley分布參數(shù)估計
    基于抗差最小均方估計的輸電線路參數(shù)辨識
    矩陣Hadamard積的上下界序列
    最大度為10的邊染色臨界圖邊數(shù)的新下界
    久久青草综合色| 少妇 在线观看| 日韩免费高清中文字幕av| 精品久久久精品久久久| 老司机深夜福利视频在线观看| 久久精品亚洲精品国产色婷小说| 欧美 亚洲 国产 日韩一| 美女福利国产在线| 亚洲精品av麻豆狂野| 国精品久久久久久国模美| 少妇 在线观看| 午夜福利欧美成人| 久久中文字幕人妻熟女| 亚洲精品乱久久久久久| 亚洲精品美女久久久久99蜜臀| 大陆偷拍与自拍| 美女主播在线视频| 9色porny在线观看| 国产精品国产高清国产av | 婷婷成人精品国产| 国产伦理片在线播放av一区| 久久性视频一级片| 高潮久久久久久久久久久不卡| 亚洲欧美日韩另类电影网站| 国产av一区二区精品久久| 18禁美女被吸乳视频| 国产极品粉嫩免费观看在线| 精品少妇内射三级| 欧美亚洲日本最大视频资源| 国产一区二区激情短视频| 久久午夜亚洲精品久久| av国产精品久久久久影院| 国产亚洲一区二区精品| 国产精品久久久人人做人人爽| 亚洲av第一区精品v没综合| 十分钟在线观看高清视频www| 中文字幕高清在线视频| 日本五十路高清| 国产亚洲欧美在线一区二区| 亚洲国产av影院在线观看| 亚洲色图综合在线观看| 激情视频va一区二区三区| 在线永久观看黄色视频| 色精品久久人妻99蜜桃| 两个人看的免费小视频| 国产亚洲av高清不卡| 午夜福利一区二区在线看| 精品福利永久在线观看| av一本久久久久| 老司机在亚洲福利影院| 一区福利在线观看| 国产精品 国内视频| 色综合欧美亚洲国产小说| 不卡av一区二区三区| 久久久国产欧美日韩av| 国产福利在线免费观看视频| 午夜激情久久久久久久| 一级黄色大片毛片| 国产欧美日韩一区二区三| 亚洲七黄色美女视频| 亚洲精品乱久久久久久| 欧美人与性动交α欧美软件| 午夜福利视频在线观看免费| 丝袜美腿诱惑在线| 成人国产av品久久久| kizo精华| 欧美人与性动交α欧美精品济南到| 天堂8中文在线网| 丰满饥渴人妻一区二区三| 国产91精品成人一区二区三区 | 欧美 日韩 精品 国产| 免费在线观看完整版高清| 一区二区三区精品91| 成人影院久久| 国产成人欧美| 午夜福利在线免费观看网站| 精品国产一区二区久久| 免费在线观看视频国产中文字幕亚洲| 在线观看免费视频网站a站| 纵有疾风起免费观看全集完整版| 午夜激情久久久久久久| 国产日韩欧美视频二区| 19禁男女啪啪无遮挡网站| 亚洲一卡2卡3卡4卡5卡精品中文| 色综合婷婷激情| 好男人电影高清在线观看| 视频区欧美日本亚洲| 18禁黄网站禁片午夜丰满| av国产精品久久久久影院| 老熟妇仑乱视频hdxx| 亚洲精品国产精品久久久不卡| 亚洲人成电影免费在线| 亚洲自偷自拍图片 自拍| 天天躁狠狠躁夜夜躁狠狠躁| 中亚洲国语对白在线视频| 成人18禁在线播放| 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| 久久精品亚洲av国产电影网| 脱女人内裤的视频| 精品人妻1区二区| 叶爱在线成人免费视频播放| 中文字幕色久视频| 欧美成人午夜精品| 男女床上黄色一级片免费看| av福利片在线| 超碰成人久久| 精品亚洲成国产av| 丝袜美腿诱惑在线| 国产人伦9x9x在线观看| 成人国产一区最新在线观看| 18禁美女被吸乳视频| 亚洲avbb在线观看| 久久久久久免费高清国产稀缺| 两个人免费观看高清视频| 亚洲国产看品久久| 正在播放国产对白刺激| 久久久精品94久久精品| 人人妻人人添人人爽欧美一区卜| 脱女人内裤的视频| 女同久久另类99精品国产91| 黑人巨大精品欧美一区二区mp4| 在线观看www视频免费| 国产亚洲精品第一综合不卡| 国产精品二区激情视频| 亚洲av日韩精品久久久久久密| 最新的欧美精品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 国产成人啪精品午夜网站| 中文字幕高清在线视频| 日本黄色视频三级网站网址 | 日日夜夜操网爽| 久久性视频一级片| 国产成人啪精品午夜网站| 97人妻天天添夜夜摸| av在线播放免费不卡| 久久人妻熟女aⅴ| 精品熟女少妇八av免费久了| 免费在线观看日本一区| 欧美日韩黄片免| 国产在线视频一区二区| 香蕉国产在线看| 久久国产精品人妻蜜桃| 亚洲专区国产一区二区| 韩国精品一区二区三区| 国产精品久久久人人做人人爽| 亚洲国产欧美一区二区综合| 欧美精品av麻豆av| 国产精品1区2区在线观看. | 国产老妇伦熟女老妇高清| 日本五十路高清| 深夜精品福利| 久久中文字幕一级| 午夜精品久久久久久毛片777| 亚洲国产av新网站| 日韩精品免费视频一区二区三区| 精品国产一区二区三区四区第35| 国产亚洲欧美在线一区二区| 久久午夜亚洲精品久久| 国产成人免费观看mmmm| 一级毛片女人18水好多| 777久久人妻少妇嫩草av网站| 日韩制服丝袜自拍偷拍| 宅男免费午夜| 高清黄色对白视频在线免费看| 精品国产乱码久久久久久小说| 夫妻午夜视频| av不卡在线播放| 乱人伦中国视频| 99在线人妻在线中文字幕 | 国产精品av久久久久免费| 最黄视频免费看| 夜夜夜夜夜久久久久| 多毛熟女@视频| 一区二区三区乱码不卡18| 精品亚洲成国产av| 90打野战视频偷拍视频| 国产精品av久久久久免费| 2018国产大陆天天弄谢| 亚洲va日本ⅴa欧美va伊人久久| 十八禁高潮呻吟视频| 纵有疾风起免费观看全集完整版| 99国产精品99久久久久| 国产欧美日韩一区二区三| 精品国产乱码久久久久久小说| 超碰成人久久| 欧美老熟妇乱子伦牲交| 日韩中文字幕视频在线看片| 亚洲五月色婷婷综合| 国产三级黄色录像| 伦理电影免费视频| 日韩成人在线观看一区二区三区| 久久香蕉激情| 少妇 在线观看| 五月开心婷婷网| 777米奇影视久久| 精品国产超薄肉色丝袜足j| videosex国产| 高清av免费在线| 精品一区二区三区视频在线观看免费 | 18禁黄网站禁片午夜丰满| 色婷婷久久久亚洲欧美| 国产精品九九99| 91麻豆av在线| 自线自在国产av| 国产成人精品无人区| 国产精品一区二区免费欧美| 丝袜喷水一区| 操出白浆在线播放| 国产精品亚洲一级av第二区| 久久精品国产亚洲av高清一级| 久久热在线av| av欧美777| 免费在线观看黄色视频的| 欧美精品啪啪一区二区三区| 久久av网站| 成人国产av品久久久| 久久久久久久久久久久大奶| 午夜福利视频精品| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| www.自偷自拍.com| 青草久久国产| 水蜜桃什么品种好| 免费女性裸体啪啪无遮挡网站| 欧美久久黑人一区二区| h视频一区二区三区| 99香蕉大伊视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品免费一区二区三区在线 | 大片免费播放器 马上看| 在线观看免费视频网站a站| 欧美日韩一级在线毛片| 国产精品亚洲一级av第二区| 99国产精品一区二区蜜桃av | 每晚都被弄得嗷嗷叫到高潮| 亚洲精华国产精华精| 成人av一区二区三区在线看| 性高湖久久久久久久久免费观看| 久久国产精品人妻蜜桃| 欧美黑人精品巨大| 亚洲免费av在线视频| 精品国产亚洲在线| 电影成人av| 国产免费视频播放在线视频| 精品乱码久久久久久99久播| www.精华液| 自线自在国产av| 午夜福利在线观看吧| 99久久99久久久精品蜜桃| 日韩精品免费视频一区二区三区| 俄罗斯特黄特色一大片| 欧美黑人精品巨大| 亚洲七黄色美女视频| av网站在线播放免费| 午夜两性在线视频| 19禁男女啪啪无遮挡网站| 国产97色在线日韩免费| 久久精品国产综合久久久| 新久久久久国产一级毛片| 久久久国产精品麻豆| 亚洲人成电影免费在线| 久久青草综合色| 欧美+亚洲+日韩+国产| 欧美日韩中文字幕国产精品一区二区三区 | 成人av一区二区三区在线看| 大香蕉久久网| 精品国产乱码久久久久久小说| 日日爽夜夜爽网站| 久久婷婷成人综合色麻豆| 欧美日本中文国产一区发布| 人人妻人人澡人人看| 99九九在线精品视频| 桃红色精品国产亚洲av| 成人手机av| 精品一区二区三区av网在线观看 | 中文字幕制服av| 午夜福利免费观看在线| 最新在线观看一区二区三区| 超碰成人久久| 国产精品二区激情视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 国产欧美日韩一区二区三区在线| 成在线人永久免费视频| 婷婷丁香在线五月| 又紧又爽又黄一区二区| 男女床上黄色一级片免费看| 女警被强在线播放| 午夜91福利影院| 91成人精品电影| 男女免费视频国产| 一级毛片电影观看| 国产一区二区三区视频了| 亚洲九九香蕉| 久久久久久久久免费视频了| 精品亚洲成a人片在线观看| 亚洲精品中文字幕一二三四区 | 免费在线观看完整版高清| 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 国产在线视频一区二区| 欧美日韩亚洲国产一区二区在线观看 | 国产欧美日韩综合在线一区二区| 一本综合久久免费| 国产91精品成人一区二区三区 | 亚洲第一av免费看| 啪啪无遮挡十八禁网站| 国产99久久九九免费精品| 亚洲精品久久成人aⅴ小说| 国产熟女午夜一区二区三区| 99国产精品99久久久久| 考比视频在线观看| 性少妇av在线| 久久国产精品影院| 人人妻人人澡人人爽人人夜夜| 国产男女内射视频| 精品一区二区三区视频在线观看免费 | 精品国产乱码久久久久久小说| 色老头精品视频在线观看| 精品少妇黑人巨大在线播放| 精品久久蜜臀av无| 亚洲欧美激情在线| 国产精品自产拍在线观看55亚洲 | av国产精品久久久久影院| tube8黄色片| 一区二区三区精品91| 夜夜骑夜夜射夜夜干| 在线天堂中文资源库| 欧美成人免费av一区二区三区 | 精品福利观看| 欧美激情极品国产一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产精品98久久久久久宅男小说| 国产精品.久久久| 久久久精品94久久精品| 大型黄色视频在线免费观看| 国产精品 国内视频| 啦啦啦视频在线资源免费观看| 午夜福利视频在线观看免费| 亚洲一区二区三区欧美精品| 操出白浆在线播放| 一级黄色大片毛片| 露出奶头的视频| 最近最新中文字幕大全电影3 | 亚洲伊人久久精品综合| 亚洲色图av天堂| 视频区欧美日本亚洲| 桃花免费在线播放| 在线观看免费视频日本深夜| 久久精品人人爽人人爽视色| 亚洲免费av在线视频| 成年人午夜在线观看视频| 一进一出好大好爽视频| 中文欧美无线码| 91九色精品人成在线观看| 国产精品免费大片| 欧美日韩亚洲高清精品| 亚洲九九香蕉| 久久影院123| 丝瓜视频免费看黄片| 高清黄色对白视频在线免费看| 国产精品一区二区免费欧美| 国产在线一区二区三区精| 精品免费久久久久久久清纯 | 看免费av毛片| h视频一区二区三区| 老司机午夜十八禁免费视频| 欧美成人午夜精品| 99精国产麻豆久久婷婷| 制服人妻中文乱码| 男女午夜视频在线观看| 欧美日韩av久久| 下体分泌物呈黄色| 90打野战视频偷拍视频| 最新美女视频免费是黄的| 精品亚洲成国产av| 久久久国产成人免费| 黄色视频在线播放观看不卡| 日韩欧美免费精品| 黄片播放在线免费| 日本欧美视频一区| 最新的欧美精品一区二区| 国产成人精品久久二区二区免费| 丁香六月天网| a级毛片在线看网站| 高潮久久久久久久久久久不卡| 日韩视频在线欧美| 免费观看a级毛片全部| 国产成人av激情在线播放| 高清av免费在线| 麻豆乱淫一区二区| 色94色欧美一区二区| 中文字幕色久视频| 国产一区二区三区视频了| 久久精品91无色码中文字幕| 精品高清国产在线一区| 精品国产一区二区三区四区第35| 国产熟女午夜一区二区三区| av电影中文网址| 一本一本久久a久久精品综合妖精| 一区二区av电影网| 国产精品免费视频内射| 国产一区二区三区在线臀色熟女 | 国产高清视频在线播放一区| cao死你这个sao货| 亚洲av片天天在线观看| 久久av网站| 国产精品免费大片| av天堂在线播放| 欧美乱妇无乱码| 热99国产精品久久久久久7| 悠悠久久av| 欧美 亚洲 国产 日韩一| 露出奶头的视频| 久久 成人 亚洲| 91麻豆精品激情在线观看国产 | 久热爱精品视频在线9| 叶爱在线成人免费视频播放| 男男h啪啪无遮挡| 黄色视频不卡| 欧美久久黑人一区二区| 久久精品aⅴ一区二区三区四区| 精品久久久久久电影网| 丝袜在线中文字幕| 国产精品成人在线| 一边摸一边做爽爽视频免费| 久久久久久久久久久久大奶| 欧美国产精品va在线观看不卡| 精品久久久久久久毛片微露脸| 黄色怎么调成土黄色| 91麻豆精品激情在线观看国产 | 国产福利在线免费观看视频| 亚洲av片天天在线观看| 国产成人精品在线电影| 精品国内亚洲2022精品成人 | 亚洲精品在线观看二区| 中文欧美无线码| 超碰成人久久| 五月开心婷婷网| 色视频在线一区二区三区| 色婷婷av一区二区三区视频| 亚洲精品国产色婷婷电影| 国产欧美日韩一区二区三区在线| 十八禁高潮呻吟视频| 亚洲欧美一区二区三区黑人| 欧美乱码精品一区二区三区| 国产日韩欧美亚洲二区| 国产欧美日韩一区二区精品| 日韩中文字幕欧美一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产av国产精品国产| 成人18禁在线播放| 久久久精品区二区三区| 国产一区二区三区在线臀色熟女 | 啦啦啦免费观看视频1| 日本一区二区免费在线视频| 一区二区三区激情视频| av网站免费在线观看视频| 男女下面插进去视频免费观看| 伦理电影免费视频| 亚洲精品久久成人aⅴ小说| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三区在线| 又大又爽又粗| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 久久久精品免费免费高清| 一进一出抽搐动态| 久久人人97超碰香蕉20202| 国产伦人伦偷精品视频| 丁香六月欧美| 精品乱码久久久久久99久播| 成人国产av品久久久| 另类亚洲欧美激情| 亚洲avbb在线观看| 一个人免费看片子| 免费不卡黄色视频| 中文字幕制服av| 久久久国产成人免费| 中文字幕制服av| 亚洲五月婷婷丁香| 免费日韩欧美在线观看| 欧美激情高清一区二区三区| 亚洲精品国产色婷婷电影| 欧美老熟妇乱子伦牲交| 亚洲精品乱久久久久久| av线在线观看网站| 99热国产这里只有精品6| 久热爱精品视频在线9| 国产精品av久久久久免费| 又黄又粗又硬又大视频| 中文字幕另类日韩欧美亚洲嫩草| av电影中文网址| 免费高清在线观看日韩| 国产免费视频播放在线视频| 人人妻人人爽人人添夜夜欢视频| svipshipincom国产片| 久久狼人影院| 国产精品久久电影中文字幕 | 又紧又爽又黄一区二区| a级毛片在线看网站| 侵犯人妻中文字幕一二三四区| 欧美精品人与动牲交sv欧美| 欧美 亚洲 国产 日韩一| 91av网站免费观看| 精品国产国语对白av| 国产午夜精品久久久久久| 国产精品一区二区免费欧美| 18禁国产床啪视频网站| 亚洲精品国产色婷婷电影| 日韩欧美一区二区三区在线观看 | 男女免费视频国产| 淫妇啪啪啪对白视频| 波多野结衣av一区二区av| 99精品久久久久人妻精品| 亚洲五月色婷婷综合| 欧美在线一区亚洲| 欧美 亚洲 国产 日韩一| 99国产综合亚洲精品| 亚洲国产中文字幕在线视频| 人人妻,人人澡人人爽秒播| 黄色片一级片一级黄色片| 又紧又爽又黄一区二区| 黑丝袜美女国产一区| 午夜日韩欧美国产| 俄罗斯特黄特色一大片| 一本久久精品| 国产男女超爽视频在线观看| 国产成人av激情在线播放| 免费在线观看完整版高清| 男女高潮啪啪啪动态图| 日韩精品免费视频一区二区三区| 多毛熟女@视频| 亚洲视频免费观看视频| 久久ye,这里只有精品| 国产欧美日韩一区二区三| 国产激情久久老熟女| 国产欧美日韩一区二区三| 午夜成年电影在线免费观看| 亚洲国产看品久久| 美女扒开内裤让男人捅视频| 午夜福利欧美成人| 精品久久蜜臀av无| 无遮挡黄片免费观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美乱妇无乱码| 麻豆乱淫一区二区| 在线看a的网站| 国产亚洲精品一区二区www | 天堂中文最新版在线下载| 久久久久网色| 满18在线观看网站| 色老头精品视频在线观看| 久久久久国内视频| 国产无遮挡羞羞视频在线观看| 两性夫妻黄色片| 18禁观看日本| 免费在线观看黄色视频的| 国产男靠女视频免费网站| 亚洲综合色网址| 精品国产一区二区三区四区第35| 国产精品免费大片| 国产主播在线观看一区二区| 亚洲avbb在线观看| 两个人免费观看高清视频| 午夜久久久在线观看| 久久亚洲精品不卡| 国产高清国产精品国产三级| 欧美黑人欧美精品刺激| 午夜久久久在线观看| 91麻豆精品激情在线观看国产 | 日韩熟女老妇一区二区性免费视频| 日日摸夜夜添夜夜添小说| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲黑人精品在线| 国产精品熟女久久久久浪| 亚洲国产欧美一区二区综合| 国产精品99久久99久久久不卡| 一级毛片女人18水好多| 久久亚洲精品不卡| 91大片在线观看| www.999成人在线观看| 亚洲av片天天在线观看| 免费女性裸体啪啪无遮挡网站| 在线 av 中文字幕| 视频在线观看一区二区三区| 老汉色∧v一级毛片| 桃花免费在线播放| 日韩精品免费视频一区二区三区| 正在播放国产对白刺激| 韩国精品一区二区三区| 老司机午夜十八禁免费视频| 欧美 日韩 精品 国产| 女警被强在线播放| 女人高潮潮喷娇喘18禁视频| av网站在线播放免费| 国产伦理片在线播放av一区| 人人澡人人妻人| a级片在线免费高清观看视频| 人人澡人人妻人| 日韩熟女老妇一区二区性免费视频| 午夜91福利影院| 人成视频在线观看免费观看| 97人妻天天添夜夜摸| 高清毛片免费观看视频网站 | 久久中文看片网| 国产精品一区二区免费欧美| 视频区欧美日本亚洲| 午夜福利影视在线免费观看| 天天操日日干夜夜撸| 亚洲精品国产区一区二| 人人妻人人澡人人看| 精品视频人人做人人爽|