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

    基于輔助陣元法的非圓信號實值MUSIC算法

    2016-11-11 08:22:55鄭春紅賈潔民鄧歡歡
    關(guān)鍵詞:實值信源元法

    鄭春紅,賈潔民,鄧歡歡,楊 剛

    (西安電子科技大學(xué)電子工程學(xué)院,陜西 西安 710071)

    ?

    基于輔助陣元法的非圓信號實值MUSIC算法

    鄭春紅,賈潔民,鄧歡歡,楊剛

    (西安電子科技大學(xué)電子工程學(xué)院,陜西 西安 710071)

    針對非圓信號波達方向(direction of arrival,DOA)估計算法在工程應(yīng)用中受陣列誤差影響的問題,基于輔助陣元法的基本原理,并結(jié)合非圓信號的特征提出了一種基于輔助陣元法的非圓信號實值多重信號分類(multiple signal classification,MUSIC)算法。該算法利用非圓信號輸出陣列的實值擴展方式,既提高陣元輸出信號的使用率,也減少算法的計算量。同時詳細(xì)推導(dǎo)了獲取信號DOA估計的數(shù)學(xué)表達式,并對所提算法進行了仿真實驗。實驗結(jié)果驗證了該算法的有效性。

    波達方向估計; 輔助陣元法; 非圓信號; 實值

    0 引 言

    陣列信號處理中,波達方向(direction of arrival,DOA)估計作為陣列信號處理的一個基本問題,在近些年得到大量的研究[1-2]。多重信號分類法是一種最常用的子空間的算法[3-5],而非圓信號作為現(xiàn)代通信系統(tǒng)中的一種常用信號,在DOA估計當(dāng)中利用其旋轉(zhuǎn)不變性改進傳統(tǒng)算法[6],非圓信號偽協(xié)方差矩陣不為零的特點可將接收數(shù)據(jù)矩陣進行擴展,使其維數(shù)加倍,有陣列擴展作用,提高DOA估計的性能[7]。同時利用協(xié)方差和偽協(xié)方差矩陣可以提高檢測及參數(shù)估計性能[8-9]。針對非圓信號的特性,一些學(xué)者提出非圓信號DOA估計算法,例如文獻[10-11]提出非圓信號的多重信號分類(multiple signal classification algorithm for noncircular signals,NC-MUSIC)算法及求根(NC-root-MUSIC)算法,文獻[12]研究了將酉旋轉(zhuǎn)不變技術(shù)信號參數(shù)(estimation of signal parameters via rotational invariance techniques,ESPRIT)算法應(yīng)用于非圓信號的DOA估計[12]等。這些算法較傳統(tǒng)的DOA估計算法性能上有了較大的提高,然而在實際工程應(yīng)用由于陣列自身誤差的存在,算法很難達到實際的理論水平。并且在實現(xiàn)過程中要進行大量的復(fù)矩陣運算,要提高工程應(yīng)用中的估計精度,是解決問題的關(guān)鍵,也是本文要展開的主題方向。

    本文基于以上現(xiàn)狀,研究了文獻中提出的基于輔助陣元法的非圓信號的DOA估計算法[13-16],并且針對該算法提出了基于輔助陣元法的非圓信號的DOA估計實值算法,該算法能夠校正方位依賴陣元的幅相誤差,只需要參數(shù)的一維搜索,并且通過有效利用信號的非圓特征,充分利用信號的信息,提高DOA估計的精度,并且能夠在一定程度上解決基于輔助陣元法的非圓信號的DOA估計算法的計算復(fù)雜度高計算時間長的問題。

    1 基于輔助陣元法的DOA估計

    1.1方向依賴陣列信號模型及非圓信號的輔助陣元法

    1.1.1方向依賴陣列模型

    針對陣列接收天線,假設(shè)由N個等距線陣組成,相鄰陣元之間的間距為d??紤]D個遠(yuǎn)場的窄帶信號以平面波形式入射到空間某陣列上,這里假設(shè)陣元個數(shù)等于通道數(shù),即各陣元接收到信號后經(jīng)各自的傳輸通道送到處理器。在這里考慮的陣元都是存在方向依賴的陣元,陣元不滿足各向同性以及陣元之間存在誤差。于是對于整個陣列在t時刻的輸出[17-18]可以表示為

    (1)

    式中,X(t)表示為陣列接收的信號矢量;W(θ)=[w(θ1),w(θ2),…,w(θD)],其中w(θi)=diag[τ(θi)]a(θi),τ(θi)是陣列對入射角為θi信號的幅相校正矢量,diag[]表示由向量中元素組成的對角矩陣,a(θi)是在理想情況下的信號導(dǎo)向矢量,是由于陣元位置不同導(dǎo)致的信號波程差引起的信號相位;S(t)為入射信號復(fù)幅度向量;N(t)為K×1陣列噪聲矢量。信號和陣元噪聲相互獨立且噪聲為零均值的高斯白噪聲。

    1.1.2非圓信號的輔助陣元法

    輔助陣元法是通過精確校準(zhǔn)的輔助陣元,可在無需已知陣列誤差參數(shù)的條件下估計信源方位,基于信源方位的準(zhǔn)確估計又可對陣列誤差進行精確校正,從而實現(xiàn)陣列誤差參數(shù)和信源方位的“去耦合”估計。對于基于輔助陣元的非圓信號自校正算法,該算法利用最大非圓信號的非圓特性進行運算處理。假設(shè)由K個陣元組成的任何幾何陣列中,其中有M個陣元為精確校準(zhǔn)的陣元,N(N=K-M)個存在方向依賴的陣元。D個遠(yuǎn)場最大非圓信號以角度為[θ1,θ2,…,θD]的平面波入射到空間陣列上。陣元編碼以第一個精準(zhǔn)陣元為參考陣元。

    按照存在方向依賴陣列信號模型,則陣列中以θi角度入射信號的導(dǎo)向矩陣可表示為

    (2)

    式中,bdiag[]表示由矩陣或者矢量組成的塊對角矩陣;τ(θi)是陣列對入射角為θi信號的幅相校正矢量;a1(θi)和a2(θi)是a(θi)按輔助陣元和非校準(zhǔn)陣元的原始信源導(dǎo)向矢量組成的分塊矩陣。將w(θi)轉(zhuǎn)換表示為

    (3)

    即b(θi)是一個K×(N+1)的矩陣,d(θi)是一個(N+1)×1的向量,所以有w(θi)=b(θi)d(θi),則在t時刻,陣列的最終輸出信號可表示為

    (4)

    在對于非圓信號的輔助陣元法DOA估計中的具體做法是,首先將存在方向依賴的信號模型引入到非圓信號的DOA估計當(dāng)中,即陣列對最大非圓信號的輸出采用式(4)的形式,然后根據(jù)NC-MUSIC算法的思想,對方向依賴陣列輸出進行共軛擴展,最終求出算法的虛擬空間譜。

    1.2基于輔助陣元法的實值DOA估計算法

    對于存在方向依賴矩陣的非圓信號陣列模型,在上文中提到的關(guān)于非圓信號DOA估計的輔助陣元法在一定程度上減小了計算量,但是在計算過程中依然是一個復(fù)數(shù)矩陣,如果對復(fù)數(shù)矩陣進行特征分解計算量依然很大,不利于工程實現(xiàn)。為此,針對輔助陣元法,同時參考文獻[19-20],提出一種基于輔助陣元法的非圓信號實質(zhì)DOA估計算法,能夠?qū)?fù)數(shù)矩陣轉(zhuǎn)換成實矩陣,更適于工程計算。

    由最大非圓信號的實質(zhì)信號相移性質(zhì)可以得到存在方向依賴矩陣在t時刻的輸出

    (5)

    式中,Φ1/2=diag[ejφ1/2,ejφ2/2,…,ejφD/2],即是由信號的非圓相位構(gòu)成的流型對角矩陣。表示對式(5)進行最大非圓信號陣列輸出實值擴展,具體形式為

    (6)

    (7)

    式中,URS表示信號子空間;URn表示噪聲子空間;ΣRS表示信號空間對應(yīng)的特征向量;ΣRn表示噪聲空間對應(yīng)的特征向量。由信號的導(dǎo)向矢量與噪聲子空間正交的性質(zhì),可以得到

    (8)

    對于輔助陣元,將式(2)代入到WR可得入射角為θ的非圓信號在輔助陣元法的實值流形矢量為

    (9)

    (10)

    則WRE的導(dǎo)向矢量形式為

    (11)

    (12)

    將式(12)中WRE寫成列向分塊形式有

    (13)

    (14)

    為了簡化計算現(xiàn)在分析QE(θ)的結(jié)構(gòu):

    (15)

    (16)

    (17)

    2 仿真分析

    非圓信號輔助陣元法通過輔助陣元實現(xiàn)對陣元誤差的自校正,從而最大限度的提高的DOA估計的精度,以下仿真實驗將利用非圓信號輔助陣元法和其他相關(guān)算法對比,體現(xiàn)該算法的優(yōu)越性能。

    仿真實驗采用的陣列為7個陣元均勻線陣,陣元間距為入射信號波長的一半,其中前3個為精準(zhǔn)校正的陣元,入射信號為窄帶二相相移鍵控(binary phase shift keying,BPSK)信號,各個入射信號功率相等,入射到陣列的噪聲為零均值高斯白噪聲。

    實驗 1兩個入射信號,其非圓相位分別為π/6和π/3,入射角度分別為-15°和20°,分別用非圓信號MUSIC算法和非圓信號輔助陣元法算法來求取信號的空間譜,為了對比算法的性能,仿真實驗分別設(shè)定信噪比為70 dB、20 dB,快拍數(shù)為200次,實驗的算法仿真結(jié)果如圖1所示。

    圖1 不同信噪比對比結(jié)果圖Fig.1 Comparison results of different SNR

    對比以上結(jié)果,不難發(fā)現(xiàn),非圓信號的輔助陣元法實值算法在不同的信噪比下能保持良好的算法性能,NC-MUSIC算法就不同,對兩個空間譜的構(gòu)造明顯有一定的偏差,信號能量比差大,角度估計偏移量相對較大,對噪聲比較敏感。

    實驗 2兩個入射信號分別以0°和15°、0°和10°、0°和5°及0°和2°方式,在信噪比為20dB入射到陣列上,200次快拍采樣仿真結(jié)果如圖2所示。

    圖2 不同角度差下的結(jié)果對比Fig.2 Comparison results of different angle

    很明顯,隨著角度差的越來越小,NC-MUSIC算法的結(jié)果越來越差,在10°角度差以下就幾乎分辨不出結(jié)果,而非圓信號輔助陣元法實值算法則相對能保持比較理想的結(jié)果,在2°的角度差時仍然可以區(qū)分出兩個目標(biāo)。

    實驗 3兩個最大非圓信號分別以一定的角度入射到陣列上,利用非圓信號NC-MUSIC算法及其實值算法進行空間譜的求取,算法在一臺配置為酷睿i3的計算機上運行。對程序仿真結(jié)果如圖3所示。

    圖3 非圓信號輔助陣元法與實值算法對比Fig.3 Comparison result between NC-MUSIC and real value algorithm

    從圖3可以看到,對于同樣的入射信號,非圓信號輔助陣元法的分辨效果與實值算法的分辨效果相同,二者具有相同的算法性能。

    同時,對200組BPSK信號進行仿真空間譜求取,同時對非圓信號輔助陣元法及其MUSIC算法運行的平均時間(單位:s)進行記錄如下,算法分別在200次和500次快拍情況下進行,結(jié)果如表1所示。

    表1 算法時間對比

    從實驗結(jié)果中可以看出,非圓信號輔助陣元實值算法的運行時間明顯低于非圓信號輔助陣元法,實值算法的確把非圓信號輔助陣元法算法計算量減少了許多。

    實驗 4在本次實驗中設(shè)計4組實驗,每組實驗的非圓信號的功率相同,同時設(shè)置信噪比均為20 dB,采用相同的快拍數(shù),分別利用輔助陣元法和非圓信號輔助陣元法進行空間譜求取信號DOA。在實驗中,第1組為2個信源,第2組為3個信源,第3組為4個信源,第4組為5個信源,仿真結(jié)果如圖4所示。

    從以上的仿真結(jié)果可以看出,當(dāng)天線陣列有7個陣元時,針對輔助陣元法,當(dāng)信源數(shù)為2時,可以很好的分辨出信源,當(dāng)信源個數(shù)為3個及其以上時,該算法已經(jīng)失去了其性能。反觀非圓信號輔助陣元法,當(dāng)信源個數(shù)為4時,算法依然可以準(zhǔn)確分辨出信源,當(dāng)信源個數(shù)為5時,該算法也失去了其性能。從實驗仿真可以知道,非圓信號輔助陣元法算法能識別的目標(biāo)個數(shù)比傳統(tǒng)的輔助陣元增加了一倍。也就是說,在相同信源個數(shù)條件下,所需的輔助陣元的個數(shù)是原始算法中陣元個數(shù)的一半。

    圖4 信源個數(shù)判斷仿真Fig.4 Comparison result about different signlas

    3 結(jié) 論

    非圓信號的輔助陣元實值算法通過在陣元中添加輔助矩陣,在進行DOA估計時消去了陣列誤差對DOA估計的影響,從而提高了DOA估計的精度。同時利用非圓信號可擴展的特性,充分利用了陣列信號的利用率,使算法能識別的目標(biāo)個數(shù)比傳統(tǒng)的輔助陣元增加了一倍。并且通過利用實值進行計算大大減少了算法的計算量,在工程實現(xiàn)上具有一定的利用價值。

    [1] Yang G,Xin J,iia X,et al.A novel DOA estimation method for closely spaced multiple sources with large power differences[C]//Proc.of the Radar Conference,2015:1276-1279.

    [2] Reddy V V,Mubeen M,Boon P N.Reduced-complexity super-resolution DOA estimation with unknown number of sources[J].Signal Processing Letters,2015,22(6):772-776.

    [3] Shu C G,Liu Y M,Yu Z Y,et al.Research on modified root-MUSIC algorithm of DOA estimation based on covariance matrix reconstruction[J].Sensors & Transducers,2014,178(9):214-218.

    [4] Pasya I,Iwakiri N,Kobayashi T.Joint direction-of-departure and direction-of-arrival estimation in an ultra-wideband MIMO radar system[C]//Proc.of the Radio and Wireless Symposium,2014:52-54.

    [5] Oh D G,Ju Y H,Lee J H.Subspace-based auto-paired range and DOA estimation of dual-channel FMCW radar without joint diagonalisation[J].Electronics Letters,2014,50(18):1320-1322.

    [6] Wu C X,Zhang M,Wang H J,et al.Array extension method based on Noncircular signal[J].Journal of Projectiles,Rockets,Missiles and Guidance,2014,34(5):165-168.(吳晨曦,張旻,王紅軍,等.基于非圓信號的陣列擴展方法[J].彈箭與制導(dǎo)學(xué)報,2014,34(5):165-168.)

    [7] Song A M,Li Y,Liu J.DOA estimation of Noncircular signals with multistage Wiener filter and polynomial rooting[J].Journal of University of Electronic Science and Technology of China,2013,42(1):53-57.(宋愛民,李堰,劉劍.非圓信號多級維納濾波DOA估計求根算法[J].電子科技大學(xué)學(xué)報,2013,42(1):53-57.)

    [8] Wang Y L.The theory and algorithms of spatial spectrum estimation[M].Beijing:Tsinghua University Press,2004:438-441.(王永良.空間譜估計理論與算法[M].北京:清華大學(xué)出版社,2004:438-441.)

    [9] Bencheikh M L,Wang Y.Combined esprit-rootmusic for DOA-DOD estimation in polarimetric bistatic MIMO radar[J].Process in Electromagnetics Research Letters,2011(22):109-117.

    [10] Liu Z M,Huang Z T,Zhou Y Y,et al.Direction-of-arrival estimation of noncircular signals via sparse representation[J].Aerospace and Electronic Systems,2012,48(3):2690-2698.

    [11] Gounon P,Adnet C,Galy J.Localization angulaire de signaux non circulaires[J].Traitement du Signal,1998,15(1):17-23.

    [12] Harrdt M,Romer F.Enhancements of unitary ESPRIT for noncircular sources[C]//Proc.of the Acoustics,Speech,and Signal Processing Conference,2004:101-104.

    [13] Cao S,Xu D Y,Xu X,et al.DOA estimation for noncircular signals in the presence of mutual coupling[J].Signal Proces-sing,2014,105(12):12-16.

    [14] Zhong M L,Fan Z Y.Direction-of-arrival estimation for noncircular signals[C]//Proc.of the International Conference on Computer,Networks and Communication Engineering,2013:10-13.

    [15] Lin M,Cao L L,Ouyang J,et al.DOA estimation using virtual array technique for noncirlular signals[C]//Proc.of the Wireless Communications & Signal Processing Conference,2012:1-5.

    [16] Shen L,Liu Z W,Gou X M,et al.Polynomial-rooting based fourth-order MUSIC for direction-of-arrival estimation of noncircular signals[J].Journal of Systems Engineering and Electronics,2014,25(6):942-948.

    [17] Liu F L,Wang J K,Sun C Y.Spatial differencing method for DOA estimation under the coexistence of both uncorrelated and coherent signals[J].Antennas and Propagation,2014,60(4):2052-2062.

    [18] El-Barbary K A,Mohamed T S,Melad M S.High resolution direction of arrival estimation (coherent signal source DOA estimation)[J].International Journal of Engineering Research and Applications,2013,3(1):132-139.

    [19] Feng D Z,Zheng C D,Zhou Y.Real-value algorithm of MUSIC via exploitation of signal sources property[J].Chinese Journal of Radio Science,2007,22(2):331-334.(馮大政,鄭春弟,周祎.一種利用信號特點的實值MUSIC算法[J].電波科學(xué)學(xué)報,2007,22(2):331-334.)

    [20] Zheng C D,Xie C W,Li Y C.Root MUSIC algorithm based on real-valued eigenvalue decomposition[J].Journal of Data Acquisition and Processing,2010,25(2):154-157.(鄭春弟,解春維,李有才.基于實值特征值分解的求根MUSIC算法[J].數(shù)據(jù)采集與處理,2010,25(2):154-157.)

    MUSIC algorithm about real values for noncircular signals based on instrumental sensors

    ZHENG Chun-hong,JIA Jie-min,DENG Huan-huan,YANG Gang

    (School of Electronic Engineering,Xidian University,Xi’an 710071,China)

    As array errors affect the precision of direction of arrival (DOA)estimation,a multiple signal classification (MUSIC)algorithm for noncircular signals about real values for noncircular signals based on instrumental sensors is proposed,which combines the basic principle of the instrumental sensor method (ISM)and the characters of noncircular signals.This new algorithm that uses the real value of the non-circular signal output arrays extending mode not only improves utilization of the array output signals,but also decreases the calculation amount.A mathematical expressions about obtaining the signal DOA estimation are derived and a simulation is made.And the results of the algorithm simulation verify the superiority of the algorithm.

    direction of arrival (DOA)estimation; instrumental sensor method (ISM); noncircular signal; real value

    2015-08-13;

    2016-02-15;網(wǎng)絡(luò)優(yōu)先出版日期:2016-05-24。

    中央高?;究蒲袠I(yè)務(wù)費專項資金(JB140235);上海航天科技創(chuàng)新基金(SAST201453)資助課題

    TN 911.7

    ADOI:10.3969/j.issn.1001-506X.2016.11.03

    鄭春紅(1969-),女,副教授,主要研究方向為智能信號信息處理。

    E-mail:chzheng@xidian.edu.cn

    賈潔民(1990-),男,碩士研究生,主要研究方向為雷達陣列信號處理。

    E-mail:anchorely2008@126.com

    鄧歡歡(1989-),男,碩士研究生,主要研究方向為雷達陣列信號處理。

    E-mail:denghuanz@qq.com

    楊剛(1968-),男,教授,博士,主要研究方向為智能信號處理。

    E-mail:ygangxidian@gmail.com

    網(wǎng)絡(luò)優(yōu)先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20160524.1731.002.html

    猜你喜歡
    實值信源元法
    多粒度實值形式概念分析
    基于極化碼的分布式多信源信道聯(lián)合編碼
    無線電工程(2022年4期)2022-04-21 07:19:44
    換元法在解題中的運用
    基于離散元法的礦石對溜槽沖擊力的模擬研究
    重型機械(2019年3期)2019-08-27 00:58:46
    實值多變量維數(shù)約簡:綜述
    信源控制電路在功率容量測試系統(tǒng)中的應(yīng)用
    電子世界(2017年16期)2017-09-03 10:57:36
    換元法在解題中的應(yīng)用
    “微元法”在含電容器電路中的應(yīng)用
    信源自動切換裝置的設(shè)計及控制原理
    雙正交周期插值小波函數(shù)的實值對稱性
    国产成人免费无遮挡视频| 精品人妻熟女毛片av久久网站| 伊人久久国产一区二区| 大片免费播放器 马上看| 日韩成人av中文字幕在线观看| 日日撸夜夜添| 亚洲怡红院男人天堂| 又大又黄又爽视频免费| 国产成人免费观看mmmm| 日本vs欧美在线观看视频| 欧美日韩精品成人综合77777| 国产熟女午夜一区二区三区 | 大香蕉久久网| 成人综合一区亚洲| 免费黄网站久久成人精品| 韩国av在线不卡| 精品一区二区三区视频在线| 99热全是精品| 999精品在线视频| 又黄又爽又刺激的免费视频.| 久久久国产欧美日韩av| 亚洲av成人精品一区久久| 啦啦啦啦在线视频资源| 99热6这里只有精品| 99热这里只有精品一区| 久久久久国产精品人妻一区二区| 国产黄色视频一区二区在线观看| 黄色怎么调成土黄色| 男女高潮啪啪啪动态图| 免费不卡的大黄色大毛片视频在线观看| 国产熟女午夜一区二区三区 | 久久狼人影院| 99久久精品一区二区三区| www.av在线官网国产| 亚洲av欧美aⅴ国产| 只有这里有精品99| 亚洲av成人精品一二三区| 在线观看国产h片| 亚洲四区av| 在线亚洲精品国产二区图片欧美 | 日本色播在线视频| 五月玫瑰六月丁香| 国产日韩一区二区三区精品不卡 | 人妻系列 视频| 亚洲五月色婷婷综合| 日韩欧美精品免费久久| 亚洲欧美一区二区三区国产| 最近的中文字幕免费完整| 一本大道久久a久久精品| 国产 精品1| 韩国av在线不卡| 欧美成人精品欧美一级黄| 在线观看人妻少妇| 亚洲,欧美,日韩| 边亲边吃奶的免费视频| 国产黄频视频在线观看| 一级毛片 在线播放| 国产精品免费大片| 有码 亚洲区| av一本久久久久| 国产成人av激情在线播放 | 一本大道久久a久久精品| 国产亚洲精品久久久com| 久久韩国三级中文字幕| 黄色毛片三级朝国网站| 青春草亚洲视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av蜜桃| 日韩 亚洲 欧美在线| 成年人午夜在线观看视频| 亚洲欧美清纯卡通| 热re99久久国产66热| 91精品国产国语对白视频| 校园人妻丝袜中文字幕| 欧美日韩视频高清一区二区三区二| 久久久精品94久久精品| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久精品久久久久真实原创| 国产黄频视频在线观看| 91午夜精品亚洲一区二区三区| 亚洲精品aⅴ在线观看| 日本黄大片高清| 成年美女黄网站色视频大全免费 | av电影中文网址| 好男人视频免费观看在线| 日韩成人av中文字幕在线观看| 色视频在线一区二区三区| 国产高清有码在线观看视频| 一本大道久久a久久精品| 欧美 亚洲 国产 日韩一| 十分钟在线观看高清视频www| 亚洲精品日韩av片在线观看| 国产成人91sexporn| 搡女人真爽免费视频火全软件| 精品久久蜜臀av无| 久久综合国产亚洲精品| 国产成人免费无遮挡视频| 日韩欧美精品免费久久| 久久久精品区二区三区| 99热这里只有精品一区| av专区在线播放| 青春草视频在线免费观看| 午夜免费观看性视频| 免费少妇av软件| 日本91视频免费播放| 国产在线一区二区三区精| 18禁观看日本| 日本wwww免费看| 久久久久视频综合| 国产黄频视频在线观看| 视频在线观看一区二区三区| 亚洲av国产av综合av卡| 男女啪啪激烈高潮av片| 3wmmmm亚洲av在线观看| 欧美日韩视频精品一区| 满18在线观看网站| 最近中文字幕2019免费版| 国产精品.久久久| 插逼视频在线观看| 你懂的网址亚洲精品在线观看| 美女大奶头黄色视频| 久久热精品热| 最新的欧美精品一区二区| 午夜福利在线观看免费完整高清在| 免费高清在线观看日韩| 国产精品一区二区在线观看99| 久久久久国产网址| 久久鲁丝午夜福利片| 日韩,欧美,国产一区二区三区| 国产成人freesex在线| 亚洲在久久综合| 国产av一区二区精品久久| 国产免费又黄又爽又色| 涩涩av久久男人的天堂| 高清在线视频一区二区三区| 亚洲不卡免费看| 精品99又大又爽又粗少妇毛片| 久久久久久久国产电影| 自线自在国产av| 老熟女久久久| 在线观看人妻少妇| tube8黄色片| 一级爰片在线观看| 精品一区二区三卡| 少妇的逼水好多| 精品一品国产午夜福利视频| 日韩一区二区视频免费看| 久久毛片免费看一区二区三区| 欧美老熟妇乱子伦牲交| 国产白丝娇喘喷水9色精品| 色网站视频免费| 伊人亚洲综合成人网| 只有这里有精品99| 日本色播在线视频| 国产精品国产三级国产av玫瑰| 九色成人免费人妻av| 中文字幕最新亚洲高清| 黄色怎么调成土黄色| 亚洲欧美清纯卡通| 精品一区二区三区视频在线| 中文乱码字字幕精品一区二区三区| 亚洲精品成人av观看孕妇| 成年av动漫网址| 亚洲三级黄色毛片| 秋霞伦理黄片| 欧美日本中文国产一区发布| av有码第一页| 久久国产精品大桥未久av| 免费观看在线日韩| av在线观看视频网站免费| 一级片'在线观看视频| 9色porny在线观看| 亚洲高清免费不卡视频| av电影中文网址| 99热这里只有是精品在线观看| 欧美 日韩 精品 国产| 久久久久久久亚洲中文字幕| 日韩欧美一区视频在线观看| 国产白丝娇喘喷水9色精品| 男女免费视频国产| freevideosex欧美| 亚洲精品aⅴ在线观看| 满18在线观看网站| 成人毛片60女人毛片免费| 中文天堂在线官网| 国产精品成人在线| 美女主播在线视频| 97超视频在线观看视频| 国产亚洲欧美精品永久| 欧美另类一区| 国国产精品蜜臀av免费| 在线观看美女被高潮喷水网站| av在线播放精品| 国语对白做爰xxxⅹ性视频网站| av不卡在线播放| 妹子高潮喷水视频| 国产高清三级在线| 国产毛片在线视频| 大话2 男鬼变身卡| 这个男人来自地球电影免费观看 | 王馨瑶露胸无遮挡在线观看| av电影中文网址| 午夜激情av网站| 婷婷成人精品国产| 中文字幕人妻熟人妻熟丝袜美| 亚洲第一av免费看| 两个人免费观看高清视频| 亚洲性久久影院| 国产精品一区二区在线观看99| 亚洲精品第二区| 欧美亚洲日本最大视频资源| 高清在线视频一区二区三区| 成年女人在线观看亚洲视频| 久久av网站| 国产成人精品在线电影| 久久人妻熟女aⅴ| 欧美精品一区二区免费开放| 亚洲国产av影院在线观看| 国产成人精品婷婷| 天堂8中文在线网| 国产精品久久久久久av不卡| 一本大道久久a久久精品| 99热这里只有精品一区| 九草在线视频观看| 久久久久久久久久久久大奶| 亚洲一区二区三区欧美精品| 69精品国产乱码久久久| 十分钟在线观看高清视频www| 一级a做视频免费观看| 日本爱情动作片www.在线观看| 国产精品一区www在线观看| 婷婷成人精品国产| av在线老鸭窝| 在线播放无遮挡| 国产成人精品久久久久久| 一区二区三区精品91| 亚洲怡红院男人天堂| 免费av中文字幕在线| 亚洲情色 制服丝袜| 丰满少妇做爰视频| 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 国产欧美亚洲国产| 日韩人妻高清精品专区| 视频区图区小说| 久久影院123| 国内精品宾馆在线| 中文字幕免费在线视频6| av专区在线播放| 久久久国产精品麻豆| av网站免费在线观看视频| 国产精品无大码| av又黄又爽大尺度在线免费看| 欧美性感艳星| 国产亚洲精品久久久com| 久久久久国产精品人妻一区二区| 丝袜喷水一区| 亚洲,欧美,日韩| 亚洲欧美精品自产自拍| 国产精品国产三级专区第一集| 黑人猛操日本美女一级片| 99九九在线精品视频| 亚洲欧美日韩另类电影网站| 国产日韩欧美在线精品| 成人免费观看视频高清| 超色免费av| 美女xxoo啪啪120秒动态图| 国产黄色视频一区二区在线观看| 亚洲av在线观看美女高潮| 午夜视频国产福利| 一级毛片黄色毛片免费观看视频| 亚洲精品,欧美精品| 国产精品人妻久久久影院| 最近中文字幕高清免费大全6| 国产精品久久久久久精品电影小说| 纯流量卡能插随身wifi吗| 一区二区日韩欧美中文字幕 | 亚洲欧美清纯卡通| 婷婷色av中文字幕| 久久狼人影院| 国产成人精品在线电影| 亚洲精品,欧美精品| 在线观看三级黄色| 卡戴珊不雅视频在线播放| 少妇熟女欧美另类| 久久久a久久爽久久v久久| 黄色一级大片看看| 这个男人来自地球电影免费观看 | 国产欧美日韩综合在线一区二区| 国产精品三级大全| 免费观看无遮挡的男女| 亚洲精品成人av观看孕妇| 亚洲欧美色中文字幕在线| videosex国产| 成人手机av| 少妇的逼好多水| 午夜老司机福利剧场| 91精品国产九色| 91精品一卡2卡3卡4卡| 男女免费视频国产| 少妇人妻 视频| 日本爱情动作片www.在线观看| 国产精品99久久99久久久不卡 | 性色av一级| h视频一区二区三区| 夜夜看夜夜爽夜夜摸| 搡女人真爽免费视频火全软件| 久久久久久久精品精品| 99热全是精品| 夜夜爽夜夜爽视频| 制服丝袜香蕉在线| 最近最新中文字幕免费大全7| 成人18禁高潮啪啪吃奶动态图 | 久久久久久人妻| 尾随美女入室| 黄片无遮挡物在线观看| 女人精品久久久久毛片| 男男h啪啪无遮挡| 欧美 日韩 精品 国产| 亚洲av不卡在线观看| 少妇的逼水好多| 纵有疾风起免费观看全集完整版| 桃花免费在线播放| 热re99久久国产66热| 中文字幕免费在线视频6| 五月开心婷婷网| 日韩精品有码人妻一区| 最近中文字幕高清免费大全6| 久久精品夜色国产| 少妇的逼水好多| 永久免费av网站大全| 国产成人免费观看mmmm| 免费观看性生交大片5| 亚洲成人一二三区av| 色5月婷婷丁香| 最新中文字幕久久久久| av国产精品久久久久影院| 新久久久久国产一级毛片| 日韩熟女老妇一区二区性免费视频| 制服诱惑二区| 久久影院123| 亚洲欧美日韩另类电影网站| 成人午夜精彩视频在线观看| 91久久精品电影网| 欧美国产精品一级二级三级| 丝袜脚勾引网站| av在线老鸭窝| 国产熟女欧美一区二区| 国产极品粉嫩免费观看在线 | 亚洲国产精品成人久久小说| 国产成人精品一,二区| 99国产综合亚洲精品| 一本—道久久a久久精品蜜桃钙片| av卡一久久| 国产亚洲精品第一综合不卡 | 亚洲国产精品999| 亚洲精品视频女| 久久久精品94久久精品| 天堂俺去俺来也www色官网| 精品亚洲成国产av| 国产av码专区亚洲av| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,一卡二卡三卡| 三上悠亚av全集在线观看| 人体艺术视频欧美日本| 各种免费的搞黄视频| 亚洲精品国产av成人精品| 国产高清三级在线| 满18在线观看网站| 一级毛片电影观看| 亚洲精品第二区| 国产伦理片在线播放av一区| av播播在线观看一区| 校园人妻丝袜中文字幕| 特大巨黑吊av在线直播| 日韩熟女老妇一区二区性免费视频| 国产精品秋霞免费鲁丝片| 亚洲国产最新在线播放| av黄色大香蕉| 狠狠婷婷综合久久久久久88av| 七月丁香在线播放| 2021少妇久久久久久久久久久| 久久国产精品男人的天堂亚洲 | 成年人免费黄色播放视频| 日本午夜av视频| 国精品久久久久久国模美| 狠狠精品人妻久久久久久综合| videosex国产| 桃花免费在线播放| 日本欧美国产在线视频| 午夜免费观看性视频| 亚洲av国产av综合av卡| 欧美日韩视频高清一区二区三区二| 丝袜在线中文字幕| 五月开心婷婷网| 日日撸夜夜添| 国内精品宾馆在线| 一级毛片aaaaaa免费看小| 男女啪啪激烈高潮av片| 免费av中文字幕在线| 麻豆精品久久久久久蜜桃| 久久久久久久久久成人| 日韩中字成人| 91精品伊人久久大香线蕉| 久久精品国产亚洲av天美| 亚洲精品aⅴ在线观看| 成人影院久久| 日本黄色片子视频| 十八禁网站网址无遮挡| 国产一区二区在线观看av| 婷婷色综合www| 秋霞伦理黄片| 少妇被粗大猛烈的视频| 一区二区三区免费毛片| 国产高清三级在线| 人成视频在线观看免费观看| 中文字幕最新亚洲高清| 亚洲国产欧美日韩在线播放| 大话2 男鬼变身卡| 国产高清国产精品国产三级| 精品99又大又爽又粗少妇毛片| 日本-黄色视频高清免费观看| 欧美少妇被猛烈插入视频| 亚洲av不卡在线观看| 亚洲欧美精品自产自拍| 99九九线精品视频在线观看视频| 国产欧美日韩一区二区三区在线 | 久久综合国产亚洲精品| 内地一区二区视频在线| 色网站视频免费| 国产亚洲欧美精品永久| 一边摸一边做爽爽视频免费| 一个人看视频在线观看www免费| 国产欧美日韩综合在线一区二区| 嫩草影院入口| 亚洲精品一区蜜桃| 亚洲国产精品999| av免费观看日本| 天堂8中文在线网| 亚洲精品中文字幕在线视频| 波野结衣二区三区在线| 九色成人免费人妻av| 日本黄色片子视频| 成人毛片a级毛片在线播放| 国产片内射在线| 韩国高清视频一区二区三区| 制服丝袜香蕉在线| 黄片播放在线免费| 久久精品国产a三级三级三级| 国产成人精品久久久久久| 精品少妇久久久久久888优播| 国产色婷婷99| 狂野欧美激情性bbbbbb| 国产在线视频一区二区| 精品久久蜜臀av无| 欧美亚洲日本最大视频资源| 日本免费在线观看一区| 久久人人爽av亚洲精品天堂| 精品人妻一区二区三区麻豆| 久久免费观看电影| 久久韩国三级中文字幕| 日韩一区二区视频免费看| 精品久久久久久久久亚洲| 在线 av 中文字幕| 国产成人freesex在线| 精品久久蜜臀av无| 免费看av在线观看网站| 精品国产一区二区三区久久久樱花| 一级毛片电影观看| 熟女人妻精品中文字幕| 亚洲精品乱码久久久v下载方式| 满18在线观看网站| 色视频在线一区二区三区| 国产高清三级在线| av在线观看视频网站免费| 午夜福利视频精品| 热99久久久久精品小说推荐| 国产伦精品一区二区三区视频9| 国产 精品1| 久久久久久久久久久免费av| √禁漫天堂资源中文www| 亚洲精品乱码久久久v下载方式| 天天躁夜夜躁狠狠久久av| 少妇猛男粗大的猛烈进出视频| 午夜视频国产福利| 少妇人妻久久综合中文| 婷婷色av中文字幕| 这个男人来自地球电影免费观看 | 久久人人爽av亚洲精品天堂| 国产高清三级在线| 亚洲精品日韩av片在线观看| 国产片特级美女逼逼视频| videosex国产| 在线观看国产h片| 又黄又爽又刺激的免费视频.| av播播在线观看一区| 国产亚洲精品久久久com| 人妻一区二区av| 亚洲精品色激情综合| 久久精品国产亚洲av天美| 亚洲一区二区三区欧美精品| 大香蕉97超碰在线| 性高湖久久久久久久久免费观看| 国产片内射在线| 久久国产精品男人的天堂亚洲 | 亚洲天堂av无毛| 水蜜桃什么品种好| 大陆偷拍与自拍| 毛片一级片免费看久久久久| 欧美激情 高清一区二区三区| 国产成人精品婷婷| 丝袜喷水一区| 五月伊人婷婷丁香| 久久人人爽人人片av| 肉色欧美久久久久久久蜜桃| 久久久久久久精品精品| 日韩精品免费视频一区二区三区 | 国产又色又爽无遮挡免| 美女国产视频在线观看| 国产精品一二三区在线看| 亚洲欧美一区二区三区黑人 | 久久久久久久久久久免费av| 制服人妻中文乱码| 国产淫语在线视频| 桃花免费在线播放| 午夜精品国产一区二区电影| 亚洲精品日韩在线中文字幕| 国产黄色视频一区二区在线观看| 满18在线观看网站| 午夜老司机福利剧场| 国产免费现黄频在线看| xxxhd国产人妻xxx| 精品人妻熟女毛片av久久网站| 男女边吃奶边做爰视频| 欧美一级a爱片免费观看看| 日日撸夜夜添| 久久久久久久久久久丰满| 国产极品粉嫩免费观看在线 | 成年女人在线观看亚洲视频| 成人国产av品久久久| 久久久国产精品麻豆| 黄片播放在线免费| 国产成人午夜福利电影在线观看| 精品午夜福利在线看| 97在线人人人人妻| 99热国产这里只有精品6| 久久精品久久久久久久性| 成年美女黄网站色视频大全免费 | 亚洲国产精品999| 国国产精品蜜臀av免费| 欧美一级a爱片免费观看看| 久久久久久伊人网av| 在线亚洲精品国产二区图片欧美 | 下体分泌物呈黄色| 亚洲国产精品一区二区三区在线| 好男人视频免费观看在线| 欧美激情极品国产一区二区三区 | 久久人人爽av亚洲精品天堂| 少妇被粗大的猛进出69影院 | 欧美三级亚洲精品| 一区二区三区免费毛片| 久久久久国产网址| 日本爱情动作片www.在线观看| 午夜福利,免费看| 成人影院久久| 午夜视频国产福利| 精品久久久久久久久av| 一级爰片在线观看| 中国三级夫妇交换| 精品亚洲成a人片在线观看| 亚洲欧美成人综合另类久久久| 国产精品蜜桃在线观看| 午夜av观看不卡| 日本黄色片子视频| 欧美日韩成人在线一区二区| 亚洲精品成人av观看孕妇| 国产精品 国内视频| a级毛片黄视频| 交换朋友夫妻互换小说| 久久久久精品性色| 日韩av不卡免费在线播放| 夜夜骑夜夜射夜夜干| 色94色欧美一区二区| 18禁在线播放成人免费| 97超碰精品成人国产| 伦理电影免费视频| 国产午夜精品一二区理论片| 男女高潮啪啪啪动态图| 久久毛片免费看一区二区三区| 久久久精品94久久精品| 国产又色又爽无遮挡免| 欧美国产精品一级二级三级| 熟女av电影| 黄片播放在线免费| 欧美激情国产日韩精品一区| 久久综合国产亚洲精品| 欧美激情 高清一区二区三区| 搡老乐熟女国产| 亚洲性久久影院| 日日爽夜夜爽网站| 国产成人精品福利久久| 美女福利国产在线| videosex国产| 午夜91福利影院| 天堂中文最新版在线下载| 免费日韩欧美在线观看| 亚洲精品视频女| 中文字幕制服av| 七月丁香在线播放| 少妇精品久久久久久久| 丰满乱子伦码专区| kizo精华| 天堂8中文在线网| 国产精品一区二区三区四区免费观看| 午夜福利视频精品|