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

    基于流形分離技術(shù)的DOA和極化參數(shù)聯(lián)合估計方法*

    2018-10-11 12:25:30胡顯智戴旭初
    測試技術(shù)學報 2018年5期
    關鍵詞:流形信號源參數(shù)估計

    胡顯智, 王 寧, 戴旭初

    (1. 中國科學技術(shù)大學 信息科學技術(shù)學院, 安徽 合肥 230027;2. 南京電子技術(shù)研究所, 江蘇 南京 210039)

    0 引 言

    信號多維參數(shù)估計是陣列信號處理的重要任務之一. 相較于由增益特性完全相同的多個天線組成的標量陣列, 極化敏感陣列(polarization sensitive array)利用內(nèi)部陣元極化選擇特性的多元化獲取信號的極化信息, 被證明可有效提高參數(shù)估計的性能[1,2]. 針對DOA和極化參數(shù)的聯(lián)合估計問題, 諸多有效算法被陸續(xù)提出: 文獻[3]提出了pencil-MUSIC算法, 可得到高精度的二維DOA和極化參數(shù)的聯(lián)合估計, 但計算量較大; 文獻[4,5]將ESPRIT算法推廣到極化敏感陣列, 得到了低復雜度的參數(shù)估計算法, 但其只適用于特定結(jié)構(gòu)的陣列. 此外, 針對相干信號源的DOA和極化參數(shù)聯(lián)合估計, 雖然也出現(xiàn)了一些算法, 如空域平滑方法[6]、 極化域平滑方法[7]等, 但這些方法只對部分具有特殊結(jié)構(gòu)的陣列有效; 實現(xiàn)一般結(jié)構(gòu)的極化敏感陣列下相干信號源的參數(shù)估計仍然是一個難題.

    流形分離技術(shù)(MST)是一種對陣列接收數(shù)據(jù)建模的新方法, 該技術(shù)源于波形域建模(wavefield modeling[8,9])思想, 即將陣列的接收數(shù)據(jù)中與信號相關的部分表示為陣列對入射信號波場的采樣. 基于這一思想, Belloni等人提出了流形分離的概念, 并利用MST將信號導向矢量分解為采樣矩陣和一個范德蒙德結(jié)構(gòu)矢量的乘積, 進而將root-MUSIC算法推廣到任意結(jié)構(gòu)陣列[10]. 隨后, Costa等人進一步考慮了方位角和俯仰角聯(lián)合估計的情況, 利用MST得到了適用于任意結(jié)構(gòu)陣列的低復雜度二維DOA估計方法[11].

    本文進一步將流形分離技術(shù)拓展到極化敏感陣列中, 研究任意結(jié)構(gòu)的極化敏感陣列下DOA和極化參數(shù)的聯(lián)合估計方法. 同時兼顧信號源部分相關和完全相干的情況. 特別是本文算法在未知陣列結(jié)構(gòu)的確切參數(shù), 但是能夠獲得實測的陣列標定數(shù)據(jù)(array calibration measurements[10])時仍然有效, 具有重要的實際應用價值.

    1 信號和陣列模型

    X=Aθ,φ,γ,ηS+N,

    (1)

    式中:S∈CK×L和N∈CM×L分別為信號采樣矩陣和噪聲采樣矩陣. 陣列流形矩陣A=[aθ1,φ1,γ1,η1,…,aθK,φK,γK,ηk]∈CM×K的第k列為第k個來波對應的信號導向矢量, 其形式為[12]

    aθk,φk,γk,ηk=Uθk,φkβΞθk,φkhγk,ηk,

    (2)

    式中: 對角陣U為信號的空域相位矩陣;β被稱為極化敏感陣列的極化敏感矩陣, 每行代表相應陣元對信號波場的響應特性;Ξ(θk,φk)=[εh,εv]由水平方向矢量εh=[-sinθk,cosθk,0]T和垂直方向矢量εp=-[sinφkcosθk,sinφksinθk,cosφk]T構(gòu)成;h=[cosγk,sinγkejηk]T為電磁波的極化矢量, 包含來波的極化信息.

    綜合上述概念, 可進一步將信號導向矢量寫為(為方便表述, 略去了下標k)

    aθ,φ,γ,η=[ah(θ,φ),av(θ,φ)]h(γ,η),

    (3)

    式中:ah(θ,φ)=Uθ,φβεh和av(θ,φ)=Uθ,φβεv被稱為水平分量導向矢量和垂直分量導向矢量, 分別表示極化敏感陣列對來自(θ,φ)方向上的水平極化電磁波和垂直極化電磁波的響應.

    由于信號源個數(shù)的估計方法較多, 也較為成熟, 本文不對信號源個數(shù)的估計問題進行討論. 不失一般性, 假設信號源的個數(shù)K和信號子空間維度Ks(即信號協(xié)方差矩陣RS的秩)已知.

    2 流形分離技術(shù)

    在陣列信號處理中, 若信號導向矢量具有范德蒙德結(jié)構(gòu), 會帶來諸多的便利, 如root-MUSIC, ESPRIT等高效的快速算法在參數(shù)估計問題中的應用; 但這一般要求陣列具有特定的結(jié)構(gòu)(如ULA等). 流形分離技術(shù)可將任意結(jié)構(gòu)陣列的信號導向矢量分解為采樣矩陣和一個范德蒙德結(jié)構(gòu)矢量的乘積, 其中: 采樣矩陣僅與陣列結(jié)構(gòu)有關, 與信號參數(shù)無關; 范德蒙德結(jié)構(gòu)矢量僅包含信號DOA參數(shù), 與陣列無關. 這使得在任意結(jié)構(gòu)的陣列下實現(xiàn)參數(shù)估計的快速算法成為可能.

    具體到極化敏感陣列, 式(3)中的水平/垂直分量導向矢量ah(θ,φ)和av(θ,φ)同時包含信號DOA參數(shù)和陣列結(jié)構(gòu)參數(shù); 流形分離技術(shù)實質(zhì)上是尋找一組正交基對信號導向矢量進行正交展開, 從而將導向矢量中與陣列相關的部分和與信號相關的部分相互分離, 即

    ah(θ,φ)=Ghd(θ,φ),av(θ,φ)=Gvd(θ,φ),

    (4)

    式中: 矩陣Gh∈CM×MaMe和Gv∈CM×MaMe即為前面所說的采樣矩陣, 這里分別稱之為水平分量采樣矩陣和垂直分量采樣矩陣, 其完全由陣列結(jié)構(gòu)決定; 正交基采用二維傅里葉基矢量d(θ,φ)∈CMaMe×1, 即有

    d(θ,φ)=d(θ)?d(φ),

    (5)

    式中: ?為克羅內(nèi)積. 正整數(shù)Ma和Me是所截取的模式數(shù): 事實上式(4)中信號導向矢量的正交展開具有無窮多諧波分量(即式(5)中ma→∞,me→∞), 但可證明當ma和me足夠大時, 后續(xù)諧波分量基本可以忽略不計[13], 從而可截取有限的模式數(shù)對信號導向矢量進行建模且保持足夠的精度, 便于工程應用(實際應用時, 可根據(jù)陣列的有效孔徑選取合適的Ma和Me).

    結(jié)合式(3)和式(4), 得到基于MST的任意結(jié)構(gòu)極化敏感陣列的信號導向矢量模型為

    aθ,φ,γ,η=[Ghd(θ,φ),Gvd(θ,φ)]h(γ,η)=G(I2?d(θ,φ))h(γ,η),

    (6)

    式中: 采樣矩陣G=[Gh,Gv]∈CM×2MaMe;I2為2×2的單位陣;h(γ,η)為前文提到的信號極化矢量.

    現(xiàn)在的問題是采樣矩陣Gh和Gv如何求解. 文獻[10]利用一種被稱為有效孔徑分布函數(shù)(Effective Aperture Distribution Function)的方法來確定采樣矩陣. 具體而言, 分別在θ∈[-π,π]和φ∈[-π,π]內(nèi)均勻選取Qa(Qa>Ma)和Qe(Qe>Me)個標定點, 得到對應的陣列標定數(shù)據(jù)

    CM×QaMe,

    (7)

    3 DOA和極化參數(shù)聯(lián)合估計

    利用2小節(jié)討論的MST, 可將式(1)中極化敏感陣列的接收數(shù)據(jù)采樣矩陣X重新表示為

    X=G(I2?D(θ,φ))V(γ,η)S+N,

    (8)

    式中: 矩陣D(θ,φ)=[d(θ1,φ1),…,d(θK,φK)]∈CMaMe×K包含來波信號的DOA參數(shù); 矩陣V(γ,η)=[Vh(γ),Vv(γ,η]T包含信號的極化參數(shù), 有

    Vh(γ)=diag{cos(γ1),…,cos(γK)}∈RK×K,

    Vv(γ,η)=diag{sin(γ1)ejη1,…,sin(γK)ejηK}∈CK×K.

    (9)

    基于式(8)所示的陣列接收數(shù)據(jù)模型, 下面討論兩種適用于任意結(jié)構(gòu)極化敏感陣列的DOA和極化參數(shù)聯(lián)合估計方法.

    3.1 基于MST的極化FFT-MUSIC算法

    在零均值圓空-時-極化白噪聲的假設下, 陣列接收數(shù)據(jù)協(xié)方差矩陣可表示為

    (10)

    式中:Es為信號子空間, 由RX的Ks個主特征矢量組成;En為噪聲子空間, 由剩余M-Ks個次特征矢量組成, 與Es相互正交. 當滿足Ks=K時(即信號源非相干), 可通過搜索式(11)的極小值點獲得來波信號參數(shù)

    f(θ,φ,γ,η)=hH(γ,η)M(θ,φ)h(γ,η),

    (11)

    式中: 矩陣M(θ,φ)∈C2×2的形式為

    (12)

    直接對式(11)進行搜索, 需要對信號的角度和極化參數(shù)進行4-D聯(lián)合搜索, 計算量很大. 為減小計算量, 利用相關矩陣理論, 可將其分兩步進行:

    (13)

    式中:Gmin{·}表示矩陣的最小特征值;Gmin{·}表示矩陣最小特征值對應的特征矢量. 由此, 整個4-D搜索過程被分解成一個2-D角度參數(shù)搜索過程和一個求解矩陣最小特征值及特征矢量的問題.

    以上即為傳統(tǒng)的極化MUSIC算法[12]的基本思路, 其中第二步的矩陣M(θ,φ)的最小特征值求解具有簡單的閉式表達, 計算量主要集中在第一步中對角度參數(shù)的二維搜索中, 其搜索間隔必須取得足夠小以滿足精度要求, 而相應需要搜索的點就會很多, 由此帶來計算量的負擔.

    利用式(8)所示的陣列輸出數(shù)據(jù)模型, 可得到基于MST的參數(shù)估計快速方法: 極化FFT- MUSIC算法. 具體而言, 結(jié)合式(6)和式(12), 可將矩陣M(θ,φ)重新表示為

    (14)

    進一步有

    (15)

    式中:d(θ)C(2Ma-1)和d(φ)∈C(2Me-1)為一維傅里葉基矢量(見式(5)); 系數(shù)矩陣Chh,Chv,Cvh和Cvv的形式為

    (16)

    式中:a?h,b?v; 而∑diag{Π,k}表示矩陣Π第k條(塊)對角線上(塊)對角元之和, 矩陣左下角元計為k=1.

    觀察式(15), 2×2矩陣M(θ,φ)的每一項實質(zhì)都是一個標準的2-D DFT形式, 從而可利用2-D FFT算法快速計算得到

    (17)

    3.2 基于MST的極化MVP-SSF算法

    對于信號源完全相干的情況, 即Ks

    考慮式(8)~式(10)的陣列模型, 但允許Ks≤K, 此時可利用信號子空間擬合技術(shù)得到來波信號的參數(shù)估計[14]

    (18)

    目標問題是一個非線性最小二乘問題, 文獻[14]采用一種被稱為MVP(Modified Variable Projection)的算法對其進行求解. MVP算法實質(zhì)是對高斯-牛頓法的一種改進, 其通過迭代完成求解

    ξi=ξi-1-μH-1V′,

    (19)

    式中:ξi是第i次迭代的結(jié)果; 梯度向量V=f(ξ)∈C4K×1, 矩陣H∈C4K×4K是Hessian矩陣2f(ξ)的近似;μ為搜索步長, 保證算法的收斂性. 關于V′和H的具體表達式, 文獻[14]中做了詳盡的推導, 將其略作擴展即可得到極化敏感陣列下的形式. 然而, 上述求解需要知道陣列流形的顯式表達, 在實際中有時可能無法獲得陣列結(jié)構(gòu)的確切參數(shù)而只能得到實測的陣列標定數(shù)據(jù), 此時無法直接使用MVP算法. 此外, 即便可獲得陣列的所有參數(shù), 當陣列結(jié)構(gòu)較為復雜時, 計算陣列流形的導數(shù)也較為繁瑣.

    利用式(8)所示基于MST的極化敏感陣列模型, 可以得到V′和H的簡單閉式表達, 且其在僅知道實測的陣列標定數(shù)據(jù)時仍然有效, 這就是極化MVP-SSF算法的基本思想. 具體而言, 利用式(8)的數(shù)據(jù)模型, 結(jié)合文獻[14]的討論, 得到V′和H閉式表達為

    (20)

    式中: 矩陣Z=[IK,IK,IK,IK]∈RK×4K, ·表示Hardmard積;JA(ξ)∈CM×4K為陣列流形矩陣A(ξ)對ξ的求導, 利用MST有閉式表達

    (21)

    (22)

    設模式數(shù)Ma=Me=M. 極化MVP-SSF算法單次迭代的計算復雜度約為O(M2MK), 主要體現(xiàn)在JA(ξ)的計算上(GΨ可提前計算好, 不需要參與迭代).

    作為一個迭代算法, 極化MVP-SSF算法需要對信號參數(shù)進行初步估計作為迭代的初始值. 信號非相干(即Ks=K)時, 極化FFT-MUSIC算法是一個理想的選擇, 其估計快速且結(jié)果可靠. 對于信號完全相干(即Ks

    4 實驗及結(jié)果分析

    圖 1 8陣元極化敏感陣均勻圓陣Fig.1 Eight-element polarimetric uniform circular array

    實驗1

    考慮如圖 1 所示的極化敏感陣列, 它由8個指向不盡相同且均勻分布在圓上的短偶極子天線組成(各陣元皆指向圓心, 其輸出電信號正比于與之平行的電場分量), 半徑為λ. 兩個等功率的非相關信號分別自(θ1=20°,φ1=40°)和(θ2=30°,φ2=50°)的方向入射至陣列; 極化參數(shù)分別為(γ1=60°,η1=30°)和(γ2=20°,η2=70°). 快拍數(shù)L=200, 噪聲為加性高斯白噪聲. 獨立實驗500次, 分別采用傳統(tǒng)的極化MUSIC算法和本文的極化FFT-MUSIC算法(Ma=Me=21)進行信號參數(shù)估計. 在信噪比SNR=10 dB下, 兩算法其中一次典型實現(xiàn)的運行時間及參數(shù)估計結(jié)果的比較如表 1 所示(搜索間隔為0.25°).

    表 1 極化FFT-MUSIC算法與傳統(tǒng)MUSIC方法運行時間的比較

    圖 2 則反映了在不同信噪比下, 兩個算法的參數(shù)估計根均方誤差的變化情況.

    圖 2 參數(shù)估計隨信噪比變化的根均方誤差Fig.2 RMSE as a function of the SNR

    由圖 2 可知, 在整個信噪比變化范圍內(nèi)(0~30 dB), 極化FFT-MUSIC算法的統(tǒng)計性能與傳統(tǒng)的極化MUSIC算法相比, 幾無差別. 這說明即便信噪比達到30 dB, 由于截取有限個(Ma=Me=21)模式所造成的誤差與噪聲相比仍然足夠小, 可以忽略不計. 而根據(jù)表1可知, 在同一臺機器下, 極化FFT-MUSIC算法進行參數(shù)估計所需要的時間遠遠小于傳統(tǒng)MUSIC算法. 綜上所述, 極化FFT-MUSIC算法在不損失原極化MUSIC算法估計精度的同時, 有效降低了計算復雜度, 具有更強的實時性.

    圖 3 參數(shù)估計隨信號源相關系數(shù)變化的根均方誤差Fig.3 RMSE as a function of correlation coefficient between the two sources

    實驗2

    考慮半徑為2λ的極化敏感均勻圓陣, 陣元數(shù)為20個, 每個陣元仍指向圓心. 信號的DOA和極化參數(shù)同實驗1. 快拍數(shù)L=200, 信噪比SNR=10 dB. 獨立實驗1 000次, 圖 3 反映了在不同的信號源相關系數(shù)下, 極化FFT-MUSIC算法和極化MVP-SSF算法的參數(shù)估計根均方誤差的變化情況. 實驗中, 兩個算法所截取的模式數(shù)皆取Ma=Me=41.

    由圖 3 可知, 在信號源相關度較低時, 極化FFT-MUSIC算法和極化MVP-SSF算法的精度基本相當; 而當相關系數(shù)大于0.6后, 極化MVP-SSF算法的參數(shù)估計性能要明顯優(yōu)于極化FFT-MUSIC算法. 特別是當信號源完全相干時(相干系數(shù)為1), 極化FFT-MUSIC算法已然失效, 而極化MVP-SSF算法仍然具有良好的性能. 由此可見, 極化MVP-SSF算法可以很好地解決強相關源乃至相干信號源的參數(shù)估計問題.

    5 結(jié) 論

    本文將流形分離技術(shù)(MST)拓展到極化敏感陣列信號處理中, 研究了兩種DOA和極化參數(shù)的聯(lián)合估計方法: 極化FFT-MUSIC算法和極化MVP-SSF算法, 它們皆適用于任意結(jié)構(gòu)的極化敏感陣列. 實驗結(jié)果表明: 極化FFT-MUSIC算法具有與傳統(tǒng)極化MUSIC算法基本相同的性能, 而計算復雜度有效降低; 極化MVP-SSF算法則可以實現(xiàn)一般結(jié)構(gòu)的極化敏感陣列下強相關源乃至相干信號源的參數(shù)估計.

    猜你喜歡
    流形信號源參數(shù)估計
    一種基于可編程邏輯器件的多功能信號源設計
    基于新型DFrFT的LFM信號參數(shù)估計算法
    緊流形上的Schr?dinger算子的譜間隙估計
    迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
    Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
    聚焦4K視頻播放展望未來信號源發(fā)展
    Logistic回歸模型的幾乎無偏兩參數(shù)估計
    基于向前方程的平穩(wěn)分布參數(shù)估計
    基于競爭失效數(shù)據(jù)的Lindley分布參數(shù)估計
    發(fā)射機信號源的自動處理和控制系統(tǒng)
    亚洲精品久久久久久婷婷小说| 国产av国产精品国产| 天天添夜夜摸| 久久免费观看电影| 精品国产超薄肉色丝袜足j| 亚洲伊人久久精品综合| 巨乳人妻的诱惑在线观看| 超碰成人久久| 日韩熟女老妇一区二区性免费视频| 一本色道久久久久久精品综合| 国产成人一区二区在线| 老司机深夜福利视频在线观看 | 国产色婷婷99| av卡一久久| 麻豆乱淫一区二区| 国产黄色免费在线视频| 妹子高潮喷水视频| 色视频在线一区二区三区| 亚洲,欧美精品.| 欧美乱码精品一区二区三区| 日韩人妻精品一区2区三区| 黄网站色视频无遮挡免费观看| 亚洲一码二码三码区别大吗| 女性被躁到高潮视频| 国产成人啪精品午夜网站| 三上悠亚av全集在线观看| 大片电影免费在线观看免费| 九草在线视频观看| 黄色视频在线播放观看不卡| 国产一区二区在线观看av| 亚洲精品美女久久久久99蜜臀 | 日日啪夜夜爽| 亚洲国产精品一区三区| 在线天堂最新版资源| 嫩草影视91久久| 亚洲精品视频女| 韩国av在线不卡| 成年女人毛片免费观看观看9 | 国产精品熟女久久久久浪| 日韩伦理黄色片| www.自偷自拍.com| 少妇人妻 视频| 精品亚洲乱码少妇综合久久| 9色porny在线观看| 国产极品粉嫩免费观看在线| 国产免费现黄频在线看| 国产精品嫩草影院av在线观看| 国产97色在线日韩免费| 国产淫语在线视频| 99久久精品国产亚洲精品| 女人高潮潮喷娇喘18禁视频| 青春草亚洲视频在线观看| 亚洲国产av影院在线观看| 亚洲精品美女久久久久99蜜臀 | 亚洲成人手机| 大陆偷拍与自拍| 国产亚洲最大av| 人体艺术视频欧美日本| 国产深夜福利视频在线观看| 制服人妻中文乱码| 99香蕉大伊视频| 国产女主播在线喷水免费视频网站| 18禁国产床啪视频网站| 亚洲av成人不卡在线观看播放网 | 亚洲成av片中文字幕在线观看| 丝袜脚勾引网站| 十分钟在线观看高清视频www| 一级,二级,三级黄色视频| av网站免费在线观看视频| 欧美日韩亚洲综合一区二区三区_| 捣出白浆h1v1| 国产一区二区三区综合在线观看| 亚洲成人一二三区av| 最近手机中文字幕大全| 精品一区二区三区av网在线观看 | 另类精品久久| 久久精品亚洲熟妇少妇任你| 久久韩国三级中文字幕| 久久99热这里只频精品6学生| 一区二区三区精品91| 亚洲欧美一区二区三区黑人| 久久精品人人爽人人爽视色| 韩国精品一区二区三区| 亚洲欧美激情在线| 亚洲av电影在线观看一区二区三区| 亚洲国产av新网站| 一区二区三区精品91| 男女边摸边吃奶| 久久精品熟女亚洲av麻豆精品| 久久久久久人人人人人| 久久韩国三级中文字幕| 午夜久久久在线观看| 欧美日韩一级在线毛片| 永久免费av网站大全| 国产一区二区激情短视频 | 国产欧美日韩综合在线一区二区| 国产福利在线免费观看视频| 看免费成人av毛片| 久久毛片免费看一区二区三区| 新久久久久国产一级毛片| 亚洲av综合色区一区| a级毛片黄视频| 大片免费播放器 马上看| 9色porny在线观看| 熟女少妇亚洲综合色aaa.| 亚洲伊人久久精品综合| 男人操女人黄网站| 欧美精品人与动牲交sv欧美| 久久精品人人爽人人爽视色| 人人妻人人爽人人添夜夜欢视频| 精品国产一区二区三区四区第35| 国产爽快片一区二区三区| 亚洲精品成人av观看孕妇| 无限看片的www在线观看| 久久久久视频综合| 国产精品国产av在线观看| 丝袜喷水一区| 国产精品久久久人人做人人爽| 久热这里只有精品99| 日韩电影二区| 女的被弄到高潮叫床怎么办| 自线自在国产av| 精品国产一区二区三区四区第35| 嫩草影视91久久| 女人被躁到高潮嗷嗷叫费观| 天天躁夜夜躁狠狠躁躁| 国产亚洲av片在线观看秒播厂| tube8黄色片| 黄色视频在线播放观看不卡| 午夜91福利影院| 精品少妇内射三级| 在现免费观看毛片| 侵犯人妻中文字幕一二三四区| 亚洲精品国产av成人精品| 老司机深夜福利视频在线观看 | 美女中出高潮动态图| 晚上一个人看的免费电影| 丝袜脚勾引网站| 亚洲 欧美一区二区三区| 亚洲伊人色综图| 乱人伦中国视频| 91成人精品电影| 黄色 视频免费看| 久久久精品区二区三区| 9热在线视频观看99| 热99国产精品久久久久久7| 国产免费一区二区三区四区乱码| 黄色怎么调成土黄色| 久久久久人妻精品一区果冻| 亚洲综合色网址| 大码成人一级视频| 日韩大片免费观看网站| 少妇人妻久久综合中文| 亚洲av成人精品一二三区| 男人操女人黄网站| 国产精品一区二区在线观看99| 十八禁网站网址无遮挡| 色网站视频免费| 男的添女的下面高潮视频| 国产精品无大码| 亚洲av男天堂| 欧美国产精品一级二级三级| 亚洲国产av影院在线观看| av免费观看日本| 成人国产av品久久久| 亚洲一区中文字幕在线| 女人爽到高潮嗷嗷叫在线视频| 国产成人a∨麻豆精品| 王馨瑶露胸无遮挡在线观看| 啦啦啦中文免费视频观看日本| 热re99久久国产66热| 黄色视频在线播放观看不卡| 久久久久久人人人人人| 欧美成人精品欧美一级黄| 999精品在线视频| 999久久久国产精品视频| 悠悠久久av| 久久精品久久久久久久性| 黑人巨大精品欧美一区二区蜜桃| 巨乳人妻的诱惑在线观看| 亚洲精品中文字幕在线视频| 最近2019中文字幕mv第一页| 国产精品免费视频内射| 日韩 亚洲 欧美在线| 黄色视频在线播放观看不卡| 国产精品99久久99久久久不卡 | 大香蕉久久网| av有码第一页| 午夜福利免费观看在线| 国产男女超爽视频在线观看| 人体艺术视频欧美日本| 中文字幕另类日韩欧美亚洲嫩草| 免费观看人在逋| 免费女性裸体啪啪无遮挡网站| 香蕉丝袜av| 久久人妻熟女aⅴ| 天天操日日干夜夜撸| 久久久久久久久久久久大奶| 午夜日本视频在线| 青草久久国产| 99香蕉大伊视频| 亚洲国产av影院在线观看| 18禁裸乳无遮挡动漫免费视频| 精品国产一区二区久久| 欧美日韩一级在线毛片| 爱豆传媒免费全集在线观看| 9色porny在线观看| 久久久久网色| 精品亚洲乱码少妇综合久久| 亚洲第一青青草原| 中文天堂在线官网| 精品久久久久久电影网| 亚洲免费av在线视频| 国产国语露脸激情在线看| 只有这里有精品99| 亚洲欧洲国产日韩| 午夜激情久久久久久久| 在线观看国产h片| av不卡在线播放| 国产高清国产精品国产三级| 自拍欧美九色日韩亚洲蝌蚪91| 欧美人与性动交α欧美精品济南到| 免费观看a级毛片全部| 高清av免费在线| 亚洲 欧美一区二区三区| 91精品伊人久久大香线蕉| 国产精品久久久久成人av| 亚洲av综合色区一区| 日韩 欧美 亚洲 中文字幕| 国产色婷婷99| 80岁老熟妇乱子伦牲交| 国产在线免费精品| 最近中文字幕2019免费版| 精品国产露脸久久av麻豆| 男的添女的下面高潮视频| 别揉我奶头~嗯~啊~动态视频 | 男女高潮啪啪啪动态图| 最黄视频免费看| 日日撸夜夜添| 欧美中文综合在线视频| 黄色怎么调成土黄色| 女人爽到高潮嗷嗷叫在线视频| 下体分泌物呈黄色| 热re99久久精品国产66热6| 亚洲av在线观看美女高潮| 一级爰片在线观看| 国产不卡av网站在线观看| av有码第一页| 成人毛片60女人毛片免费| 久久久亚洲精品成人影院| 久久性视频一级片| 少妇人妻 视频| 亚洲精品国产一区二区精华液| 久久热在线av| 9191精品国产免费久久| 国产片特级美女逼逼视频| 精品国产超薄肉色丝袜足j| 国产乱人偷精品视频| 欧美黑人精品巨大| 亚洲中文av在线| av网站免费在线观看视频| 亚洲熟女精品中文字幕| 亚洲成国产人片在线观看| 成人18禁高潮啪啪吃奶动态图| 好男人视频免费观看在线| 国产黄色视频一区二区在线观看| 亚洲国产成人一精品久久久| 两个人免费观看高清视频| 19禁男女啪啪无遮挡网站| 精品国产乱码久久久久久男人| 欧美亚洲 丝袜 人妻 在线| 夫妻午夜视频| 在线 av 中文字幕| 精品亚洲成国产av| 美女高潮到喷水免费观看| av在线观看视频网站免费| 久久久久人妻精品一区果冻| 亚洲欧美日韩另类电影网站| 国产无遮挡羞羞视频在线观看| 国产精品久久久久久精品古装| 秋霞伦理黄片| 亚洲七黄色美女视频| 日韩大片免费观看网站| 久久99热这里只频精品6学生| 久久精品久久精品一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产无遮挡羞羞视频在线观看| 日韩电影二区| 国产成人精品在线电影| 国产精品国产三级专区第一集| 亚洲av成人精品一二三区| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的| 亚洲精品久久成人aⅴ小说| 亚洲欧美成人精品一区二区| 欧美日韩一区二区视频在线观看视频在线| 亚洲熟女精品中文字幕| 女人精品久久久久毛片| 欧美中文综合在线视频| www日本在线高清视频| 午夜精品国产一区二区电影| 宅男免费午夜| 久久女婷五月综合色啪小说| 日韩大码丰满熟妇| 成年女人毛片免费观看观看9 | 欧美日韩av久久| 亚洲av成人不卡在线观看播放网 | 各种免费的搞黄视频| av国产精品久久久久影院| 国产精品一区二区精品视频观看| 自线自在国产av| 亚洲男人天堂网一区| 成人毛片60女人毛片免费| 午夜老司机福利片| 成人漫画全彩无遮挡| 成年人午夜在线观看视频| 亚洲精品国产av蜜桃| 国产又色又爽无遮挡免| 51午夜福利影视在线观看| 国产国语露脸激情在线看| 久久 成人 亚洲| 韩国精品一区二区三区| 亚洲国产日韩一区二区| 精品福利永久在线观看| 日日撸夜夜添| 亚洲av电影在线观看一区二区三区| 久久久精品94久久精品| 91精品伊人久久大香线蕉| www.自偷自拍.com| 亚洲国产精品国产精品| 国产日韩欧美亚洲二区| 不卡视频在线观看欧美| 亚洲成色77777| 最黄视频免费看| 国产成人免费无遮挡视频| 久久久久久久大尺度免费视频| 九草在线视频观看| bbb黄色大片| 老司机深夜福利视频在线观看 | 日韩,欧美,国产一区二区三区| 只有这里有精品99| 成人18禁高潮啪啪吃奶动态图| 亚洲三区欧美一区| 伊人久久大香线蕉亚洲五| 美女午夜性视频免费| xxxhd国产人妻xxx| 亚洲欧美日韩另类电影网站| 久热这里只有精品99| 老鸭窝网址在线观看| 久久久久精品国产欧美久久久 | 九色亚洲精品在线播放| h视频一区二区三区| 人体艺术视频欧美日本| 午夜av观看不卡| 国产人伦9x9x在线观看| 建设人人有责人人尽责人人享有的| 亚洲成人一二三区av| 午夜福利,免费看| 91aial.com中文字幕在线观看| 久久精品久久精品一区二区三区| 中文字幕制服av| 久久精品国产亚洲av涩爱| 亚洲av日韩在线播放| 国产97色在线日韩免费| 卡戴珊不雅视频在线播放| 黄色 视频免费看| 国产伦理片在线播放av一区| 一级,二级,三级黄色视频| 高清av免费在线| 亚洲精品日本国产第一区| 欧美精品高潮呻吟av久久| 老司机影院毛片| 亚洲国产av影院在线观看| 91精品国产国语对白视频| 久久久久久人人人人人| 亚洲第一青青草原| 免费高清在线观看视频在线观看| 国产片内射在线| 麻豆乱淫一区二区| 国产不卡av网站在线观看| 亚洲激情五月婷婷啪啪| 久久久精品94久久精品| 一级a爱视频在线免费观看| 成人午夜精彩视频在线观看| 日韩人妻精品一区2区三区| 免费久久久久久久精品成人欧美视频| 热99国产精品久久久久久7| 一级爰片在线观看| 国产精品av久久久久免费| 国产97色在线日韩免费| 高清欧美精品videossex| 无遮挡黄片免费观看| 国产伦人伦偷精品视频| 可以免费在线观看a视频的电影网站 | 亚洲精品国产一区二区精华液| 水蜜桃什么品种好| 另类亚洲欧美激情| 精品国产乱码久久久久久小说| 精品久久久精品久久久| 男人添女人高潮全过程视频| 日韩大码丰满熟妇| 日日摸夜夜添夜夜爱| 日本91视频免费播放| 日韩av不卡免费在线播放| 在线观看国产h片| 超碰97精品在线观看| 欧美黑人欧美精品刺激| 午夜老司机福利片| www日本在线高清视频| 国产一卡二卡三卡精品 | 一级毛片 在线播放| 啦啦啦视频在线资源免费观看| 自线自在国产av| 欧美激情高清一区二区三区 | 欧美av亚洲av综合av国产av | 日韩免费高清中文字幕av| 波多野结衣av一区二区av| 美女大奶头黄色视频| 我的亚洲天堂| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 精品一区在线观看国产| 久久久精品区二区三区| 捣出白浆h1v1| 一本色道久久久久久精品综合| 一本大道久久a久久精品| 国产精品香港三级国产av潘金莲 | 黄频高清免费视频| 欧美国产精品一级二级三级| 纯流量卡能插随身wifi吗| 亚洲成人av在线免费| 在线观看三级黄色| 侵犯人妻中文字幕一二三四区| 人人妻人人爽人人添夜夜欢视频| 国产熟女午夜一区二区三区| 色播在线永久视频| tube8黄色片| 电影成人av| 黄片播放在线免费| 韩国精品一区二区三区| 黄频高清免费视频| 精品一区在线观看国产| 老司机在亚洲福利影院| 97精品久久久久久久久久精品| 午夜福利,免费看| 制服诱惑二区| 丰满少妇做爰视频| 18禁国产床啪视频网站| 亚洲av男天堂| 亚洲熟女毛片儿| 丝袜美腿诱惑在线| 久久综合国产亚洲精品| av在线app专区| 超色免费av| 老司机亚洲免费影院| 精品少妇黑人巨大在线播放| 午夜av观看不卡| 午夜日韩欧美国产| 18禁观看日本| 天天操日日干夜夜撸| 日日爽夜夜爽网站| 欧美黑人精品巨大| 国产一区二区在线观看av| 我要看黄色一级片免费的| 成年人午夜在线观看视频| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 亚洲欧美色中文字幕在线| 日本vs欧美在线观看视频| 国产爽快片一区二区三区| 亚洲天堂av无毛| 丰满饥渴人妻一区二区三| 中文天堂在线官网| 色婷婷久久久亚洲欧美| 欧美激情高清一区二区三区 | 午夜福利,免费看| 国产亚洲最大av| 日韩欧美一区视频在线观看| 国产精品人妻久久久影院| videos熟女内射| 日韩制服骚丝袜av| 熟妇人妻不卡中文字幕| 最黄视频免费看| 悠悠久久av| 午夜av观看不卡| 久久精品久久久久久久性| 超色免费av| 麻豆精品久久久久久蜜桃| 日韩av在线免费看完整版不卡| 黄色视频在线播放观看不卡| 亚洲精品aⅴ在线观看| 女性被躁到高潮视频| 黑丝袜美女国产一区| 午夜日韩欧美国产| 90打野战视频偷拍视频| 王馨瑶露胸无遮挡在线观看| 人人妻人人爽人人添夜夜欢视频| 最近中文字幕高清免费大全6| 丝袜美腿诱惑在线| 国产日韩欧美视频二区| 少妇人妻 视频| 久久精品熟女亚洲av麻豆精品| 精品一区二区三区av网在线观看 | 我的亚洲天堂| 香蕉丝袜av| 麻豆乱淫一区二区| 少妇人妻久久综合中文| 国产又色又爽无遮挡免| 丝袜美足系列| 日韩av不卡免费在线播放| 在线亚洲精品国产二区图片欧美| 亚洲综合色网址| 99国产精品免费福利视频| 国产在视频线精品| 只有这里有精品99| 精品亚洲成a人片在线观看| 激情五月婷婷亚洲| 国产成人精品在线电影| 最近的中文字幕免费完整| 久久久久精品久久久久真实原创| 日韩电影二区| 日日爽夜夜爽网站| 欧美精品一区二区大全| 中文字幕亚洲精品专区| 一边亲一边摸免费视频| 岛国毛片在线播放| 欧美日韩精品网址| 成人国产麻豆网| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站| 国产精品久久久久久人妻精品电影 | 国产免费现黄频在线看| 最近2019中文字幕mv第一页| 国产国语露脸激情在线看| 在线观看免费视频网站a站| 国产av国产精品国产| 男女下面插进去视频免费观看| 老司机靠b影院| 中文精品一卡2卡3卡4更新| 久久久久精品人妻al黑| 国产片特级美女逼逼视频| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看一区二区三区激情| 午夜免费鲁丝| 满18在线观看网站| 三上悠亚av全集在线观看| 国产男女内射视频| 色94色欧美一区二区| 最近的中文字幕免费完整| 国产一区二区在线观看av| 久久毛片免费看一区二区三区| 王馨瑶露胸无遮挡在线观看| 亚洲精品久久久久久婷婷小说| av在线app专区| 国产日韩一区二区三区精品不卡| 大陆偷拍与自拍| 国产精品一二三区在线看| 国产精品亚洲av一区麻豆 | 美国免费a级毛片| 91国产中文字幕| 国产日韩欧美在线精品| 亚洲,一卡二卡三卡| 男的添女的下面高潮视频| 自线自在国产av| 成人影院久久| av.在线天堂| 国产一级毛片在线| 大片电影免费在线观看免费| 美女高潮到喷水免费观看| 亚洲欧美一区二区三区国产| 极品人妻少妇av视频| 熟女av电影| 亚洲四区av| 少妇被粗大猛烈的视频| 欧美变态另类bdsm刘玥| 国语对白做爰xxxⅹ性视频网站| 国产精品成人在线| 欧美在线黄色| 国产av国产精品国产| 日韩人妻精品一区2区三区| 精品亚洲乱码少妇综合久久| 国产精品女同一区二区软件| 满18在线观看网站| 自线自在国产av| av天堂久久9| av.在线天堂| 欧美日韩亚洲高清精品| 一级毛片电影观看| 香蕉丝袜av| 飞空精品影院首页| 丰满乱子伦码专区| 免费观看人在逋| 18禁动态无遮挡网站| 激情五月婷婷亚洲| 欧美日韩一区二区视频在线观看视频在线| 美女视频免费永久观看网站| 国产亚洲午夜精品一区二区久久| 亚洲成人av在线免费| 国产精品久久久久久精品古装| 免费久久久久久久精品成人欧美视频| 国产成人精品在线电影| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| 国产精品一区二区在线观看99| 建设人人有责人人尽责人人享有的| 大码成人一级视频| 久久精品人人爽人人爽视色| 一区二区三区激情视频| 高清av免费在线| 国产高清不卡午夜福利| 国产一级毛片在线| 久久久精品94久久精品| 免费观看av网站的网址| 国产爽快片一区二区三区| 丰满迷人的少妇在线观看|