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

    FDA-MIMO雷達(dá)目標(biāo)距離、角度及未知互耦參數(shù)估計(jì)

    2021-11-11 06:14:46喬理?yè)P(yáng)崔忠馬楊赟秀
    關(guān)鍵詞:參數(shù)估計(jì)天線雷達(dá)

    喬理?yè)P(yáng), 崔忠馬, 楊赟秀, 舒 勤,*

    (1. 四川大學(xué)電氣工程學(xué)院, 四川 成都 610065; 2. 北京遙感設(shè)備研究所, 北京 100854;3. 西南技術(shù)物理研究所, 四川 成都 610045)

    0 引 言

    多輸入多輸出(multiple input multiple output,MIMO)技術(shù)是一種具有許多優(yōu)勢(shì)的靈活技術(shù)[1-3]。MIMO雷達(dá)通過(guò)使用多個(gè)天線同時(shí)發(fā)射相互正交的信號(hào),并使用多個(gè)天線接收回波。這使得MIMO雷達(dá)具有比傳統(tǒng)的相控陣?yán)走_(dá)更大的自由度,導(dǎo)致MIMO雷達(dá)的可識(shí)別性和分辨率顯著提高。頻控陣(frequency diverse array,FDA)與相控陣具有相似的結(jié)構(gòu),作為一種新型陣列雷達(dá)在2006年被首次提出[4]。FDA通過(guò)控制每個(gè)天線的載波頻率來(lái)控制具有不同距離-角度關(guān)系的波束形式[5]。正是這種距離-角度依賴(lài)的特性,使得FDA有了在距離和角度兩個(gè)維度上定位目標(biāo)的能力。MIMO技術(shù)與FDA在文獻(xiàn)[6]中被首次合并,自此針對(duì)FDA-MIMO雷達(dá)的研究成了熱點(diǎn)問(wèn)題。

    陣列的參數(shù)估計(jì)研究有悠久的歷史,許多優(yōu)秀的估計(jì)方法被提出。例如基于旋轉(zhuǎn)不變技術(shù)的信號(hào)參數(shù)估計(jì)技術(shù)(estimation of signal parameters via rotational invariance technique,ESPRIT)和多重信號(hào)分類(lèi)(multiple signal classification,MUSIC)技術(shù)[7-8]。目前,大多數(shù)的估算方法都假定天線陣列是理想的。但實(shí)際中,天線陣列始終存在未知的誤差。典型的天線誤差包括陣列間的相互耦合,這是由天線元件的輻射效應(yīng)引起的必然存在的誤差[9-10]。天線陣列互耦合會(huì)使得信號(hào)發(fā)射與回波參數(shù)估計(jì)帶來(lái)障礙。發(fā)射天線陣列互耦會(huì)影響電壓駐波比,可能導(dǎo)致能量輻射失敗[11-12];接收陣列互耦合可能導(dǎo)致目標(biāo)參數(shù)估計(jì)誤差增大。因此,消除陣列互耦合對(duì)提高陣列與估計(jì)算法性能是非常必要的。本文重點(diǎn)研究接收陣列存在互耦合對(duì)目標(biāo)距離和角度估計(jì)的影響,以及消除互耦合的參數(shù)估計(jì)方法。

    傳統(tǒng)MIMO雷達(dá)的參數(shù)估計(jì)方法中,有許多針對(duì)存在陣列耦合的情況。文獻(xiàn)[13-14]提出了均勻線性陣列(uniform linear array,ULA)校準(zhǔn)的方法,并且總結(jié)了ULA的信號(hào)模型與互耦合矩陣模型。文獻(xiàn)[15]中提出了用貝葉斯學(xué)習(xí)來(lái)研究到達(dá)角(direction of arrival,DOA)估計(jì)性能和陣列自校準(zhǔn)的問(wèn)題。文獻(xiàn)[16]提出一種DOA和互耦合系數(shù)聯(lián)合估計(jì)的方法,從此陣列參數(shù)和互耦合系數(shù)的聯(lián)合估計(jì)成為了新的研究方向。文獻(xiàn)[17]從信號(hào)角度分析信號(hào)的額外信息,從而進(jìn)行參數(shù)估計(jì)。文獻(xiàn)[18]中從信號(hào)矩陣結(jié)構(gòu)角度思考,提出了一種數(shù)據(jù)轉(zhuǎn)換策略方法來(lái)進(jìn)行參數(shù)估計(jì),然后通過(guò)迭代求解互耦合系數(shù),迭代過(guò)程浪費(fèi)了巨大的計(jì)算資源。文獻(xiàn)[19]提出了一種遞進(jìn)方法,首先避免耦合影響并估計(jì)DOA,然后使用估計(jì)得到的DOA通過(guò)解決半定規(guī)劃問(wèn)題估計(jì)耦合系數(shù),迭代收斂的過(guò)程耗時(shí)較多,效率低下。文獻(xiàn)[20]改進(jìn)了迭代的方法,通過(guò)建立拉格朗日乘數(shù)法解決二次優(yōu)化問(wèn)題,計(jì)算效率提升。

    上述互耦合系數(shù)估計(jì)方法均適用于傳統(tǒng)MIMO雷達(dá)的信號(hào)結(jié)構(gòu)。目前,鮮見(jiàn)有研究FDA-MIMO雷達(dá)互耦合校準(zhǔn)方法,由于天線陣陣列間電磁輻射引起的互耦合效應(yīng)不可避免,因而對(duì)FDA-MIMO雷達(dá)互耦合系數(shù)估計(jì)及影響分析缺失,對(duì)于FDA-MIMO雷達(dá)陣列性能評(píng)估是不充分的。本文針對(duì)FDA-MIMO雷達(dá)陣列間存在未知互耦合問(wèn)題,提出一種天線間存在耦合情況下的距離和角度估計(jì)技術(shù)和耦合系數(shù)估計(jì)方法。

    本文首先將傳統(tǒng)陣列耦合校準(zhǔn)技術(shù)拓展到FDA-MIMO雷達(dá),得到存在互耦的FDA-MIMO雷達(dá)接收信號(hào)模型,對(duì)此模型分解得到信號(hào)子空間和噪聲子空間,利用信號(hào)子空間和噪聲子空間的正交性得到譜函數(shù),通過(guò)一種數(shù)據(jù)轉(zhuǎn)換方法,成功得到了距離和角度的估計(jì)值。再通過(guò)類(lèi)似線性約束最小方差重構(gòu)優(yōu)化問(wèn)題,求解未知耦合系數(shù)。仿真結(jié)果證明了方法的有效性。

    1 信號(hào)模型

    考慮一個(gè)具有M個(gè)發(fā)射陣元和N個(gè)接收陣元的單基地均勻線性FDA-MIMO雷達(dá),則第m個(gè)陣元發(fā)送的復(fù)信號(hào)可以表示為

    sm(t)=φm(t)ej2πfmt, 0≤t≤T;m=1,2,…,M

    (1)

    式中:T為雷達(dá)脈沖持續(xù)時(shí)間;fm第m個(gè)天線上面增加的頻率偏移,并且有fm=f1+Δfm,f1為參考載波頻率,Δfm為第m(m>1)個(gè)陣元的頻率偏移,該頻率偏移相對(duì)載波頻率來(lái)說(shuō)可以忽略不計(jì);φm(t)表示單位能量波形,且波形滿(mǎn)足正交條件。

    電磁波在空間任意傳播,假設(shè)雷達(dá)的遠(yuǎn)場(chǎng)空間中存在n個(gè)目標(biāo),每個(gè)目標(biāo)均會(huì)反射來(lái)自發(fā)射天線的信號(hào),則接收陣列中,第n個(gè)接收陣元接收來(lái)自k個(gè)目標(biāo)的反射信號(hào)表示為

    (2)

    (3)

    式中:rk為第θk個(gè)目標(biāo)的徑向距離;θk為第k個(gè)目標(biāo)的角度;c為光速;dT和dR分別為發(fā)射陣列陣元間距和接收陣列陣元間距。

    在接收端,通過(guò)下變頻和匹配濾波處理接收到的信號(hào),輸出信號(hào)的第l個(gè)快拍可以表示為

    y(l)=As(l)+n(l)

    (4)

    式中:A=[aR(θ1)?aT(r1,θ1),aR(θ2)?aT(r2,θ2),…,aR(θk)?aT(rk,θk)]∈CMN×K為導(dǎo)向矢量矩陣,?為Kronecker積;s(l)=[s1(l),s2(l),…,sK(l)]T∈CK×1為包含各目標(biāo)復(fù)反射系數(shù)的信號(hào)矢量;n(l)∈CMN×1為獨(dú)立0均值復(fù)高斯分布的噪聲向量;aT(rk,θk)∈CM×1為發(fā)射陣列的導(dǎo)向矢量,aR(θk)∈CN×1為接收陣列的導(dǎo)向矢量,其中,

    (5)

    (6)

    在發(fā)射和接收陣列陣元間存在耦合的情況下,陣列發(fā)射和接收的信號(hào)將會(huì)產(chǎn)生列相關(guān)效應(yīng),且aT和aR均會(huì)被互耦合矩陣影響[20]。

    具有M×M個(gè)陣列天線的天線間互耦合效應(yīng)可以通過(guò)M×M的互耦合矩陣來(lái)描述,其矩陣為可以表示為

    (7)

    在互耦合矩陣C中,元素ci, j為陣列中第i個(gè)陣元與第j個(gè)陣元之間的耦合系數(shù)。ci, j的大小與陣元之間的距離成反比,隨著陣元間距超過(guò)某一定值,耦合將可忽略不計(jì)[13]。ULA的耦合矩陣具有對(duì)稱(chēng)Toeplitz結(jié)構(gòu)[21],設(shè)CT和CR分別為發(fā)射陣列和接收陣列的互耦合系數(shù)矩陣,并且假設(shè)發(fā)射陣列與接收陣列的每行各有P和Q個(gè)非零互耦合系數(shù),則

    (8)

    cT=[cT1,cT2,…,cTP]T和cR=[cR1,cR2,…,cRQ]T為互耦合系數(shù)向量。

    因此,L個(gè)快拍信號(hào)模型寫(xiě)為矩陣形式為

    Y=(CRAR)?(CTAT)S+N=CAS+N

    (9)

    式中:C=CR?CT∈CMN×MN為耦合矩陣;S∈CK×L為信號(hào)矩陣;N為高斯噪聲矩陣。

    2 所提算法

    2.1 距離與角度估計(jì)

    根據(jù)子空間算法,首先求得Y的協(xié)方差矩陣:

    R=E{YYH}=CARsAHCH+σ2IMN

    (10)

    式中:Rs是信號(hào)協(xié)方差矩陣;σ2是噪聲方差;IMN為單位矩陣。

    對(duì)式(10)進(jìn)行特征值分解,得到:

    (11)

    可以使用MUSIC的方法,遍歷搜索區(qū)域中的距離Z和角度Θ,可以用以下譜函數(shù)的K個(gè)峰估算所有目標(biāo)的距離和角度[22]:

    (12)

    式中:F(Z,Θ)的K個(gè)最大峰值將對(duì)應(yīng)估計(jì)的源的角度和距離。

    但是,上述譜函數(shù)因?yàn)殛嚵形粗ヱ詈闲?yīng)的存在,譜函數(shù)會(huì)受未知互耦合矩陣C影響,形成嚴(yán)重的相關(guān)效應(yīng),估計(jì)性能急劇下降。為了消除陣列互耦合的影響,可以使用下述方法。

    首先,令:

    (13)

    式中:ΦT∈CM×P與ΦR∈CN×P均為與導(dǎo)向矢量相關(guān)的轉(zhuǎn)換矩陣;CT和CR為影響發(fā)射陣列和接收陣列的耦合矩陣,設(shè)φ=1,2,…,P,φ=1,2,…,Q。并且,令

    (14)

    (15)

    經(jīng)過(guò)上述轉(zhuǎn)換,譜函數(shù)可以表達(dá)為

    (16)

    (17)

    則距離和角度的估計(jì)值可表示為

    (18)

    2.2 互耦合系數(shù)估計(jì)

    在估計(jì)距離和角度的過(guò)程中需要避免c為0的特殊情況,通過(guò)類(lèi)似線性約束最小方差來(lái)重構(gòu)此問(wèn)題。在這里定義eTc=1。其中,e=[1,0,…,0]T∈CM×1。

    因此,式(18)可以表達(dá)為

    min(cHΠc)

    s.t.eTc=1

    (19)

    構(gòu)建代價(jià)函數(shù):

    L(Z,Θ)=cHΠc-λ(eTc-1)

    (20)

    式中:λ是常數(shù)。對(duì)c求偏導(dǎo),令?/?cL(Z,Θ)=0得

    2Πc+λe=0

    (21)

    從式(21)中獲得

    c=μΠ-1e

    (22)

    其中,μ為常數(shù):

    (23)

    由式(22)和式(23),可以得到互耦合系數(shù)的估計(jì)值:

    (24)

    3 仿真分析

    在仿真中,M和N均設(shè)置為20,dT和dR均設(shè)置為0.5c/fM+B,B≈Δf2,f1設(shè)置為10 GHz,對(duì)每個(gè)天線的頻率偏移采用對(duì)數(shù)法則[23],即Δf=lnmΔf,其中Δf=20 kHz。發(fā)射陣列和接收陣列的互耦系數(shù)分別設(shè)置為cR=[1,0.8+0.5j,0.1+0.2j]T,cT=[1,0.6+0.3j,0.2+0.1j]T??紤]空間中有兩個(gè)待測(cè)目標(biāo)其距離和角度分別為(80 km,20°)和(130 km,55°),信源發(fā)射信號(hào)快拍數(shù)為200,如圖1所示。

    圖1 發(fā)射陣列與接收陣列示意圖Fig.1 Schematic diagram of transmitting array and receiving array

    仿真 1所提算法對(duì)目標(biāo)進(jìn)行估計(jì)的結(jié)果

    采用所提算法對(duì)目標(biāo)參數(shù)進(jìn)行估計(jì)的結(jié)果如圖2所示。角度的搜索范圍為[20 km,200 km],間隔為[20 km,200 km];距離的搜索范圍為[20 km,200 km],間隔為1 m,信噪比(signal to noise ratio, SNR)設(shè)置為10 dB。可以看出,在收發(fā)陣列存在互耦合影響的情況下,所提算法在不需要預(yù)知互耦合系數(shù)的前提下,就可以較為精確地獲得目標(biāo)距離和角度的估計(jì)值,這是本文算法的優(yōu)勢(shì)。由于FDA-MIMO結(jié)構(gòu)本身不能消除互耦合,因而直接使用傳統(tǒng)估計(jì)算法結(jié)果是不理想的。圖2為參數(shù)估計(jì)的譜結(jié)果三維圖,可以清晰看到目標(biāo)的距離和角度參數(shù)估計(jì)結(jié)果,圖中強(qiáng)峰為每個(gè)目標(biāo)的距離和角度參數(shù)。

    圖2 所提算法的譜圖(三維)Fig.2 Spectra of the proposed algorithm (three dimensional)

    圖3是圖2譜結(jié)果圖在距離角度平面上的投影,其中橫坐標(biāo)為角度,縱坐標(biāo)為距離,圖中紅框位置為譜峰,譜峰即為距離和角度的估計(jì)值。

    圖3 所提算法的譜圖(二維)Fig.3 Spectra of the proposed algorithm (two dimensional)

    仿真 2算法估計(jì)距離和角度的統(tǒng)計(jì)性能

    定義距離和角度估計(jì)結(jié)果的均方根誤差(root mean aquare error, RMSE)分別為

    (25)

    式中:K為目標(biāo)個(gè)數(shù);W為蒙特卡羅實(shí)驗(yàn)次數(shù);仿真中W=500。SNR與角度和距離結(jié)果的RMSE的關(guān)系如圖4所示。

    圖4 角度和距離的RMSE與SNR的關(guān)系Fig.4 Relationship between RMSE and SNR of angle and range

    圖4顯示,在陣列互耦合存在的情況下,所提算法在一定SNR范圍內(nèi)均具備較好的估計(jì)性能。文獻(xiàn)[18]和文獻(xiàn)[19]中算法同樣可以進(jìn)行參數(shù)估計(jì)但是誤差相對(duì)較大,傳統(tǒng)MUSIC方法則在互耦合情況下性能喪失,所提方法在保證準(zhǔn)確估計(jì)的前提下,提升了估計(jì)準(zhǔn)度。

    仿真 3算法估計(jì)互耦合系數(shù)的統(tǒng)計(jì)性能

    定義互耦合系數(shù)估計(jì)結(jié)果的RMSE為

    (26)

    根據(jù)所提算法分別對(duì)發(fā)射陣列耦合系數(shù)cR和接收陣列耦合系數(shù)cT進(jìn)行估計(jì)。圖5給出了耦合系數(shù)估計(jì)時(shí)的SNR與RMSE的關(guān)系。

    圖5 耦合系數(shù)估計(jì)的RMSE與SNR的關(guān)系Fig.5 Relationship between RMSE and SNR of mutual coupling coefficient estimation

    由圖5可以看到,在一定SNR范圍內(nèi),耦合系數(shù)估計(jì)的RMSE隨SNR的增加而減小,所提方法與對(duì)比算具有相似的估計(jì)性能,所提方法體現(xiàn)更優(yōu)的估計(jì)結(jié)果。

    4 結(jié) 論

    本文提出一種FDA-MIMO雷達(dá)陣列間存在互耦合效應(yīng)時(shí)的參數(shù)估計(jì)方法。該方法首先構(gòu)造了存在未知互耦合陣列的接收信號(hào)模型,通過(guò)子空間分解的方法得到用于參數(shù)估計(jì)的譜函數(shù);并且采用一種數(shù)據(jù)轉(zhuǎn)換策略,解決了譜函數(shù)由于未知耦合矩陣存在而失真的問(wèn)題,最終得到了FDA距離和角度參數(shù)的估計(jì)值;通過(guò)類(lèi)似線性約束最小方差重構(gòu)優(yōu)化問(wèn)題,估計(jì)出耦合系數(shù)矩陣。仿真結(jié)果顯示,本文提出的方法解決了發(fā)射陣和接收陣共同存在耦合情況的參數(shù)估計(jì)問(wèn)題,并且避免了迭代方法造成的資源浪費(fèi),提升了計(jì)算效率,并且估計(jì)性能均優(yōu)于現(xiàn)有算法。

    猜你喜歡
    參數(shù)估計(jì)天線雷達(dá)
    有雷達(dá)
    大自然探索(2023年7期)2023-08-15 00:48:21
    基于新型DFrFT的LFM信號(hào)參數(shù)估計(jì)算法
    雷達(dá)
    ETC相控陣天線與普通天線應(yīng)用對(duì)比分析
    Logistic回歸模型的幾乎無(wú)偏兩參數(shù)估計(jì)
    ALLESS轉(zhuǎn)動(dòng)天線射頻旋轉(zhuǎn)維護(hù)與改造
    電子制作(2016年1期)2016-11-07 08:42:54
    基于向前方程的平穩(wěn)分布參數(shù)估計(jì)
    理論宣講要上接天線、下接地氣
    基于競(jìng)爭(zhēng)失效數(shù)據(jù)的Lindley分布參數(shù)估計(jì)
    基于空時(shí)二維隨機(jī)輻射場(chǎng)的彈載雷達(dá)前視成像
    秋霞伦理黄片| 亚洲欧美日韩无卡精品| 久久人人爽人人爽人人片va| 天天躁日日操中文字幕| 精品99又大又爽又粗少妇毛片| 日韩欧美一区视频在线观看 | 不卡视频在线观看欧美| 免费高清在线观看视频在线观看| 一本—道久久a久久精品蜜桃钙片| 综合色丁香网| av在线老鸭窝| 欧美最新免费一区二区三区| 国内揄拍国产精品人妻在线| 亚洲,一卡二卡三卡| 国产视频内射| 精品少妇黑人巨大在线播放| 亚洲成人中文字幕在线播放| 啦啦啦视频在线资源免费观看| 成年免费大片在线观看| 久久久久久久久久成人| 哪个播放器可以免费观看大片| 国产视频首页在线观看| 天堂俺去俺来也www色官网| 人妻 亚洲 视频| 免费观看在线日韩| 日韩欧美 国产精品| 亚洲精品日韩av片在线观看| 国国产精品蜜臀av免费| av卡一久久| 亚洲av不卡在线观看| 麻豆成人av视频| 边亲边吃奶的免费视频| 一级毛片电影观看| 免费在线观看成人毛片| 国产伦在线观看视频一区| 亚洲一区二区三区欧美精品| 多毛熟女@视频| 久久人人爽av亚洲精品天堂 | 国产成人精品一,二区| 99re6热这里在线精品视频| 久久久久久久久久久免费av| 亚洲成人av在线免费| av国产免费在线观看| 亚洲国产精品专区欧美| 亚洲av.av天堂| 高清av免费在线| 啦啦啦中文免费视频观看日本| 欧美日韩国产mv在线观看视频 | 免费观看在线日韩| 一区在线观看完整版| 色哟哟·www| 嫩草影院入口| videossex国产| 五月天丁香电影| 亚洲av不卡在线观看| 欧美高清成人免费视频www| 人人妻人人澡人人爽人人夜夜| 全区人妻精品视频| 91久久精品电影网| 亚洲怡红院男人天堂| av又黄又爽大尺度在线免费看| 亚洲伊人久久精品综合| 美女国产视频在线观看| av国产免费在线观看| 欧美精品人与动牲交sv欧美| 美女脱内裤让男人舔精品视频| 五月开心婷婷网| 黄色一级大片看看| 欧美激情极品国产一区二区三区 | 亚洲av欧美aⅴ国产| 交换朋友夫妻互换小说| 国产极品天堂在线| 国产精品av视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩亚洲高清精品| 欧美日韩一区二区视频在线观看视频在线| 午夜激情久久久久久久| 久久久久性生活片| www.av在线官网国产| 最黄视频免费看| 最近中文字幕高清免费大全6| 一级av片app| 久久久久久久久久成人| 亚洲综合色惰| 国产美女午夜福利| 99re6热这里在线精品视频| 3wmmmm亚洲av在线观看| 性色avwww在线观看| 中文在线观看免费www的网站| 美女脱内裤让男人舔精品视频| 各种免费的搞黄视频| 2021少妇久久久久久久久久久| 99热这里只有是精品50| 汤姆久久久久久久影院中文字幕| 深爱激情五月婷婷| 蜜桃久久精品国产亚洲av| 成人亚洲精品一区在线观看 | 国产精品伦人一区二区| 老司机影院毛片| av免费在线看不卡| 亚洲欧美一区二区三区国产| 亚洲人成网站在线观看播放| 美女主播在线视频| 国产精品久久久久成人av| 在线观看一区二区三区激情| av在线观看视频网站免费| 久久综合国产亚洲精品| 伊人久久国产一区二区| 日韩人妻高清精品专区| 91久久精品国产一区二区成人| 丝瓜视频免费看黄片| 久久国产精品大桥未久av | 国产精品久久久久成人av| 精品一区二区三卡| 国产精品无大码| 亚洲,一卡二卡三卡| 伊人久久国产一区二区| 我的老师免费观看完整版| 麻豆国产97在线/欧美| 秋霞在线观看毛片| 亚洲经典国产精华液单| 日本猛色少妇xxxxx猛交久久| 免费黄色在线免费观看| 女人久久www免费人成看片| 亚洲伊人久久精品综合| 99热全是精品| 不卡视频在线观看欧美| 亚洲av男天堂| 亚洲欧美中文字幕日韩二区| 亚洲欧美日韩东京热| 亚洲美女视频黄频| 青青草视频在线视频观看| 成人二区视频| 免费观看av网站的网址| 中国国产av一级| 婷婷色av中文字幕| 免费观看av网站的网址| 搡老乐熟女国产| 建设人人有责人人尽责人人享有的 | 一个人看的www免费观看视频| 七月丁香在线播放| 最近2019中文字幕mv第一页| 女性生殖器流出的白浆| 2018国产大陆天天弄谢| 欧美日韩亚洲高清精品| 国内少妇人妻偷人精品xxx网站| 我要看日韩黄色一级片| 国产爽快片一区二区三区| 高清午夜精品一区二区三区| 国产精品女同一区二区软件| 欧美高清性xxxxhd video| 大陆偷拍与自拍| 亚洲自偷自拍三级| 国产69精品久久久久777片| 精品酒店卫生间| 最近最新中文字幕免费大全7| 国产淫片久久久久久久久| 亚洲精华国产精华液的使用体验| 国产成人91sexporn| 国产精品一区www在线观看| 中国美白少妇内射xxxbb| 3wmmmm亚洲av在线观看| 水蜜桃什么品种好| 成人亚洲欧美一区二区av| 男女边吃奶边做爰视频| 老熟女久久久| 成人毛片60女人毛片免费| 在线观看国产h片| 欧美极品一区二区三区四区| 大话2 男鬼变身卡| 久久6这里有精品| 国产在线一区二区三区精| 三级国产精品欧美在线观看| 91在线精品国自产拍蜜月| 乱码一卡2卡4卡精品| 午夜福利高清视频| 毛片一级片免费看久久久久| 少妇的逼水好多| 美女主播在线视频| 国产黄色视频一区二区在线观看| 三级经典国产精品| 激情五月婷婷亚洲| 精品人妻一区二区三区麻豆| 午夜日本视频在线| 18禁裸乳无遮挡免费网站照片| 日韩一区二区视频免费看| 欧美日韩视频精品一区| 男人爽女人下面视频在线观看| 久久亚洲国产成人精品v| 又大又黄又爽视频免费| 秋霞在线观看毛片| 日本黄大片高清| 黄色怎么调成土黄色| 老师上课跳d突然被开到最大视频| 22中文网久久字幕| 亚洲欧美清纯卡通| 国产探花极品一区二区| 国产男女内射视频| 国产日韩欧美亚洲二区| 国产一区二区三区av在线| 午夜老司机福利剧场| av视频免费观看在线观看| 看非洲黑人一级黄片| 热99国产精品久久久久久7| 久久av网站| 亚洲国产精品专区欧美| 免费看av在线观看网站| 男女免费视频国产| 久久韩国三级中文字幕| 久久毛片免费看一区二区三区| 日本wwww免费看| 国产 精品1| 欧美激情极品国产一区二区三区 | 99国产精品免费福利视频| 国产精品一及| 妹子高潮喷水视频| 国产爱豆传媒在线观看| 能在线免费看毛片的网站| 久久精品熟女亚洲av麻豆精品| 亚洲激情五月婷婷啪啪| 最新中文字幕久久久久| 性高湖久久久久久久久免费观看| 男女免费视频国产| 麻豆国产97在线/欧美| 国产免费一区二区三区四区乱码| 国产乱来视频区| 99久久精品一区二区三区| 亚洲欧美精品自产自拍| 亚洲欧洲日产国产| 日本av手机在线免费观看| 99热这里只有是精品在线观看| 高清视频免费观看一区二区| 人体艺术视频欧美日本| 一区二区av电影网| 成人18禁高潮啪啪吃奶动态图 | 久久婷婷青草| 国产成人freesex在线| 国产精品国产三级专区第一集| 久久99热这里只有精品18| 一级毛片我不卡| 精品99又大又爽又粗少妇毛片| 极品少妇高潮喷水抽搐| 80岁老熟妇乱子伦牲交| 成人亚洲欧美一区二区av| 一区在线观看完整版| 欧美国产精品一级二级三级 | 热99国产精品久久久久久7| 免费黄频网站在线观看国产| 中国美白少妇内射xxxbb| 免费av不卡在线播放| 在线观看av片永久免费下载| 亚洲自偷自拍三级| 国产在视频线精品| 乱码一卡2卡4卡精品| 永久免费av网站大全| 国产av国产精品国产| 久久毛片免费看一区二区三区| 中国三级夫妇交换| 亚洲不卡免费看| 亚洲真实伦在线观看| 亚洲综合精品二区| h日本视频在线播放| 美女xxoo啪啪120秒动态图| 国产一区有黄有色的免费视频| 精品国产乱码久久久久久小说| 国产 精品1| 美女脱内裤让男人舔精品视频| 日韩成人av中文字幕在线观看| av专区在线播放| 日本wwww免费看| 国产片特级美女逼逼视频| 91精品一卡2卡3卡4卡| 爱豆传媒免费全集在线观看| 99热这里只有精品一区| 美女主播在线视频| 美女视频免费永久观看网站| 男女边摸边吃奶| 毛片女人毛片| 欧美日韩综合久久久久久| 欧美一区二区亚洲| av专区在线播放| 日韩大片免费观看网站| 久久久久久久大尺度免费视频| 中文天堂在线官网| av女优亚洲男人天堂| 极品少妇高潮喷水抽搐| 草草在线视频免费看| 成人二区视频| 亚洲av国产av综合av卡| 欧美xxxx性猛交bbbb| 国产有黄有色有爽视频| 性色avwww在线观看| 国产成人精品福利久久| 国产综合精华液| 99久久精品一区二区三区| 色吧在线观看| 国产 一区 欧美 日韩| 免费播放大片免费观看视频在线观看| 中文字幕亚洲精品专区| 3wmmmm亚洲av在线观看| 欧美成人精品欧美一级黄| 亚洲精品亚洲一区二区| 亚洲精品日本国产第一区| 色网站视频免费| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 国产成人精品婷婷| 联通29元200g的流量卡| 亚洲av福利一区| 亚洲av在线观看美女高潮| 亚洲欧美精品专区久久| 麻豆精品久久久久久蜜桃| 男的添女的下面高潮视频| 2021少妇久久久久久久久久久| 国产高清不卡午夜福利| 六月丁香七月| 亚洲经典国产精华液单| 99久久综合免费| 亚洲国产日韩一区二区| 久久久久久久精品精品| 久久久久久久久久成人| 校园人妻丝袜中文字幕| 99国产精品免费福利视频| 国产欧美另类精品又又久久亚洲欧美| 欧美xxxx性猛交bbbb| 欧美日韩亚洲高清精品| 高清视频免费观看一区二区| 亚洲精品视频女| 欧美少妇被猛烈插入视频| 亚洲欧美中文字幕日韩二区| 少妇的逼水好多| 亚洲欧美一区二区三区黑人 | 我的女老师完整版在线观看| 嫩草影院入口| 建设人人有责人人尽责人人享有的 | 成人二区视频| 成人高潮视频无遮挡免费网站| xxx大片免费视频| 丝瓜视频免费看黄片| 26uuu在线亚洲综合色| 国产v大片淫在线免费观看| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 国产黄色免费在线视频| 秋霞在线观看毛片| 亚洲欧美精品专区久久| 久热久热在线精品观看| 精品国产乱码久久久久久小说| 国产精品伦人一区二区| 看非洲黑人一级黄片| 欧美成人精品欧美一级黄| 亚洲熟女精品中文字幕| 国产精品蜜桃在线观看| 精品久久久久久久久亚洲| 免费看av在线观看网站| 黄片wwwwww| 国产一级毛片在线| 美女主播在线视频| 国产久久久一区二区三区| 成人黄色视频免费在线看| 亚洲成人一二三区av| 亚洲精品,欧美精品| 韩国高清视频一区二区三区| 一区二区三区精品91| 99热6这里只有精品| 亚洲精品自拍成人| 性色av一级| 国产精品偷伦视频观看了| 免费在线观看成人毛片| 午夜日本视频在线| 妹子高潮喷水视频| 26uuu在线亚洲综合色| 久久久久久久国产电影| 亚洲伊人久久精品综合| 久久久久久久亚洲中文字幕| 免费不卡的大黄色大毛片视频在线观看| 97在线人人人人妻| 美女国产视频在线观看| 成人美女网站在线观看视频| 高清av免费在线| 日韩电影二区| 狂野欧美激情性xxxx在线观看| 国产av一区二区精品久久 | 国产在线一区二区三区精| 精品久久久久久电影网| 国产成人精品一,二区| 中文精品一卡2卡3卡4更新| av黄色大香蕉| 亚洲精品色激情综合| 国产色爽女视频免费观看| 国产精品国产三级国产av玫瑰| 一级爰片在线观看| 色哟哟·www| 国产欧美日韩精品一区二区| 伦理电影免费视频| 人妻一区二区av| 久久综合国产亚洲精品| 女人久久www免费人成看片| 99久久综合免费| 大陆偷拍与自拍| 中文字幕精品免费在线观看视频 | 韩国av在线不卡| 男人添女人高潮全过程视频| 亚洲精品自拍成人| 日韩欧美精品免费久久| 少妇的逼水好多| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 高清不卡的av网站| 色网站视频免费| 国产精品久久久久成人av| 在线看a的网站| 色婷婷久久久亚洲欧美| 亚洲,欧美,日韩| 3wmmmm亚洲av在线观看| 亚洲成人av在线免费| 精品久久久久久久久av| 香蕉精品网在线| 欧美精品一区二区大全| 久久午夜福利片| 性高湖久久久久久久久免费观看| 免费观看性生交大片5| 欧美人与善性xxx| 国产国拍精品亚洲av在线观看| 精品久久久噜噜| 高清日韩中文字幕在线| 国产精品久久久久久久久免| 亚洲一区二区三区欧美精品| 99久久精品国产国产毛片| 中国美白少妇内射xxxbb| 亚洲av.av天堂| 国产精品蜜桃在线观看| 欧美国产精品一级二级三级 | 一级黄片播放器| 国产女主播在线喷水免费视频网站| 国产成人a∨麻豆精品| 亚洲在久久综合| 在线观看国产h片| 人妻 亚洲 视频| 色吧在线观看| 高清毛片免费看| 日本午夜av视频| 日韩一区二区三区影片| 中文天堂在线官网| 日韩亚洲欧美综合| av.在线天堂| av国产久精品久网站免费入址| 亚洲欧美日韩东京热| 久久综合国产亚洲精品| 久久婷婷青草| 国产欧美另类精品又又久久亚洲欧美| 香蕉精品网在线| 精品国产乱码久久久久久小说| 国产久久久一区二区三区| 亚洲精品色激情综合| 小蜜桃在线观看免费完整版高清| 午夜免费鲁丝| 夫妻性生交免费视频一级片| 国产精品嫩草影院av在线观看| 国产永久视频网站| 亚洲av综合色区一区| 日韩精品有码人妻一区| 国产精品麻豆人妻色哟哟久久| 免费av中文字幕在线| 麻豆成人av视频| 美女内射精品一级片tv| 永久免费av网站大全| 性色av一级| 亚洲国产av新网站| 男人添女人高潮全过程视频| 综合色丁香网| 在线天堂最新版资源| 黄色视频在线播放观看不卡| 哪个播放器可以免费观看大片| 国产日韩欧美亚洲二区| 精品一区在线观看国产| 国产白丝娇喘喷水9色精品| 国内精品宾馆在线| 亚洲真实伦在线观看| 高清毛片免费看| 国产亚洲欧美精品永久| 人妻一区二区av| 18禁在线无遮挡免费观看视频| 永久免费av网站大全| 热re99久久精品国产66热6| 精品国产三级普通话版| 久久婷婷青草| 国精品久久久久久国模美| 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 精品久久久久久久久亚洲| 又黄又爽又刺激的免费视频.| 国产精品欧美亚洲77777| 人妻制服诱惑在线中文字幕| 久久久色成人| 亚洲天堂av无毛| tube8黄色片| 久久久久人妻精品一区果冻| av天堂中文字幕网| 人妻少妇偷人精品九色| 哪个播放器可以免费观看大片| 观看免费一级毛片| 偷拍熟女少妇极品色| 黑丝袜美女国产一区| 亚洲欧美一区二区三区国产| 免费播放大片免费观看视频在线观看| 成年女人在线观看亚洲视频| 午夜免费鲁丝| 制服丝袜香蕉在线| 美女视频免费永久观看网站| 在线观看三级黄色| 国产在线视频一区二区| 免费观看的影片在线观看| 免费人妻精品一区二区三区视频| 直男gayav资源| 日本猛色少妇xxxxx猛交久久| 日韩,欧美,国产一区二区三区| 午夜福利影视在线免费观看| 在线观看免费日韩欧美大片 | 男人狂女人下面高潮的视频| 97超视频在线观看视频| 欧美激情国产日韩精品一区| 大香蕉97超碰在线| av播播在线观看一区| 国产 一区 欧美 日韩| 久久影院123| .国产精品久久| 亚洲精品第二区| 精品一区二区免费观看| 蜜桃在线观看..| 六月丁香七月| 2022亚洲国产成人精品| 精品人妻一区二区三区麻豆| 成人漫画全彩无遮挡| 成人亚洲欧美一区二区av| 韩国av在线不卡| 国产成人免费观看mmmm| 97超碰精品成人国产| 亚洲国产最新在线播放| 日本与韩国留学比较| 国产色爽女视频免费观看| 久久人人爽av亚洲精品天堂 | 天堂俺去俺来也www色官网| av免费在线看不卡| 亚洲,欧美,日韩| 男人和女人高潮做爰伦理| av国产免费在线观看| 精品视频人人做人人爽| 一区二区三区四区激情视频| av网站免费在线观看视频| 纯流量卡能插随身wifi吗| 国产国拍精品亚洲av在线观看| 岛国毛片在线播放| 国产 精品1| 亚洲国产av新网站| 欧美日本视频| 日韩免费高清中文字幕av| 国产精品无大码| 欧美97在线视频| 免费观看性生交大片5| 欧美精品一区二区大全| 丝瓜视频免费看黄片| 激情 狠狠 欧美| 日本免费在线观看一区| 妹子高潮喷水视频| 一级片'在线观看视频| 亚洲最大成人中文| 水蜜桃什么品种好| 美女国产视频在线观看| 国产精品99久久久久久久久| av不卡在线播放| 国产 一区精品| 亚洲丝袜综合中文字幕| 精品久久久久久久末码| 国产精品国产三级专区第一集| 一个人免费看片子| 免费大片18禁| 亚洲精品aⅴ在线观看| 国产91av在线免费观看| 噜噜噜噜噜久久久久久91| 久久女婷五月综合色啪小说| 韩国高清视频一区二区三区| 国产男人的电影天堂91| h日本视频在线播放| a级毛色黄片| 自拍欧美九色日韩亚洲蝌蚪91 | www.色视频.com| 伦精品一区二区三区| 性色avwww在线观看| 精品人妻偷拍中文字幕| 久久精品久久久久久久性| av又黄又爽大尺度在线免费看| 大香蕉久久网| 日本欧美国产在线视频| 国产亚洲5aaaaa淫片| 国产久久久一区二区三区| 国产精品一区二区在线不卡| 国产免费福利视频在线观看| 内地一区二区视频在线| 多毛熟女@视频| 欧美日韩一区二区视频在线观看视频在线| 婷婷色av中文字幕| 亚洲欧洲日产国产| 在线天堂最新版资源| 亚洲国产成人一精品久久久| 国产精品爽爽va在线观看网站| 午夜福利在线观看免费完整高清在| 国产爽快片一区二区三区| 成人漫画全彩无遮挡| 纵有疾风起免费观看全集完整版| 亚洲国产精品国产精品| 久久久久久久精品精品| 亚洲av中文av极速乱| tube8黄色片| 天天躁日日操中文字幕|