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

    基于Sobel算子的雙基地MIMO雷達(dá)陣列診斷方法

    2019-12-23 09:01:18陳金立卓齊剛竇思鈺翟介新董雨晴
    關(guān)鍵詞:信噪比算子灰度

    陳金立,卓齊剛,竇思鈺,翟介新,董雨晴

    (1. 南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,南京 210044;2. 南京信息工程大學(xué),南京 210044)

    0 引 言

    多輸入多輸出(Multiple-Input Multiple-Output,MIMO)雷達(dá)是在數(shù)字陣列雷達(dá)、多基地雷達(dá)以及MIMO通信技術(shù)的基礎(chǔ)上發(fā)展起來的一種新體制雷達(dá)[1-2]。與傳統(tǒng)的相控陣?yán)走_(dá)相比,MIMO雷達(dá)能靈活地根據(jù)需求設(shè)計(jì)發(fā)射波形,還可通過虛擬陣元的形成來擴(kuò)展收發(fā)陣列的孔徑,以提升目標(biāo)可識(shí)別數(shù)目的上限,并提高目標(biāo)分辨力以及參數(shù)估計(jì)性能[3-5]。天線陣列單元是雷達(dá)的關(guān)鍵組成部分,能夠?qū)崿F(xiàn)天線波束快速掃描、空間功率合成和多波束形成等功能,但長期工作于高頻高功率下的收發(fā)組件容易損壞,從而導(dǎo)致相應(yīng)的陣元失效。陣元失效導(dǎo)致其輸出信號(hào)僅為熱噪聲,而不包含空間目標(biāo)信息,從而導(dǎo)致MIMO雷達(dá)探測性能的明顯下降。MIMO雷達(dá)探測性能的恢復(fù)需已知失效陣元的位置,因此研究有效的天線陣列診斷方法以實(shí)現(xiàn)失效陣元位置的高成功概率診斷具有重要的理論與實(shí)際意義[6-7]。

    針對天線陣列失效單元的診斷問題,一些學(xué)者提出了一類采用波導(dǎo)探頭測量天線陣列的遠(yuǎn)場或近場數(shù)據(jù)的陣列診斷方法。文獻(xiàn)[8]提出了一種反向傳播算法,該方法利用波導(dǎo)探頭采樣被測陣列的近場數(shù)據(jù),再對近場數(shù)據(jù)作傅里葉變換來獲得陣列的口徑場分布特性,從而診斷出失效陣元的位置。文獻(xiàn)[9]提出了一種矩陣算法,該算法首先建立探頭電壓與陣列激勵(lì)系數(shù)的線性關(guān)系,然后利用矩陣求逆算法求解該線性關(guān)系以獲得陣列激勵(lì)系數(shù),最后通過激勵(lì)系數(shù)獲得陣列的遠(yuǎn)場方向圖從而檢測出失效陣元的位置。文獻(xiàn)[8-9]所提的陣列診斷方法要求使用探頭以半波長為步進(jìn)間隔采集陣列近場數(shù)據(jù),同時(shí)要求采樣數(shù)量不少于陣元數(shù)目,因此隨著陣列規(guī)模的增大將導(dǎo)致上述方法出現(xiàn)采樣數(shù)據(jù)大、測量時(shí)間長等問題。壓縮感知 (Compressive Sensing, CS) 理論的出現(xiàn)為陣列診斷提供了全新的思路,該方法能夠利用較少數(shù)量的采樣數(shù)據(jù)檢測失效陣元,逐漸成為陣列診斷領(lǐng)域研究的熱點(diǎn)。文獻(xiàn)[10]將貝葉斯壓縮感知理論應(yīng)用于失效陣元位置的診斷,克服了傳統(tǒng)CS算法的測量矩陣必須滿足約束等距性質(zhì)的限制。文獻(xiàn)[11]提出了一種基于lp范數(shù)最小化的均勻線陣的陣列診斷方法,能進(jìn)一步地減少遠(yuǎn)場數(shù)據(jù)的采樣數(shù)量。然而當(dāng)失效陣元的數(shù)量較多時(shí),失效陣元的激勵(lì)系數(shù)不再滿足稀疏性,導(dǎo)致基于CS理論的陣列診斷方法失效。波導(dǎo)探頭的使用會(huì)增加硬件成本并干擾待測場,文獻(xiàn)[12]直接對雷達(dá)工作期間的陣列接收數(shù)據(jù)進(jìn)行差分處理來診斷失效陣元的位置,無需使用額外的波導(dǎo)探頭,但需要遍歷所有失效陣元組合情況下各差分位置上虛擬陣元的重復(fù)次數(shù)來進(jìn)行失效陣元位置的搜索,且任意差分位置的虛擬陣元重復(fù)次數(shù)的錯(cuò)誤估計(jì)均會(huì)導(dǎo)致陣列診斷失敗,因此該方法穩(wěn)定性較差,同時(shí)當(dāng)陣列中的失效陣元位置使同一差分位置的虛擬陣元同時(shí)失效時(shí),該方法無法準(zhǔn)確診斷出失效陣元的位置。

    針對上述問題及研究現(xiàn)狀,本文提出了一種基于Sobel算子的雙基地MIMO雷達(dá)陣列診斷方法,該方法利用失效陣元與正常陣元的圖像邊緣強(qiáng)度的差異特性,將失效陣元位置的診斷問題轉(zhuǎn)化為圖像中較低邊緣強(qiáng)度位置的診斷問題。本文方法能夠在無需使用波導(dǎo)探頭的前提下快速有效地診斷出雙基地MIMO雷達(dá)失效陣元的位置。

    1 陣元失效MIMO雷達(dá)信號(hào)模型

    假設(shè)雙基地MIMO雷達(dá)系統(tǒng)由M個(gè)發(fā)射陣元和N個(gè)接收陣元組成,發(fā)射陣列和接收陣列都為均勻線性分布陣列,發(fā)射和接收陣元的間距分別為dt和dr,假設(shè)空間中存在L個(gè)遠(yuǎn)場窄帶目標(biāo),第l(l=1,…,L)個(gè)目標(biāo)相對于發(fā)射陣列和接收陣列的發(fā)射角(DOD)和接收角(DOA)分別為(φl,θl)。各發(fā)射陣元同時(shí)發(fā)射相同載頻和帶寬的周期性正交信號(hào)B=[b1,…,bM]T,bm為第m個(gè)發(fā)射陣元的發(fā)射信號(hào),(·)T表示轉(zhuǎn)置。則k時(shí)刻發(fā)射信號(hào)經(jīng)L個(gè)目標(biāo)反射由接收陣列接收的信號(hào)為

    (1)

    式中,βl為第l個(gè)目標(biāo)的反射系數(shù);ar(θl)=[1,e-j2πdrsin(θl)/λ,…,e-j2π(N-1)drsin(θl)/λ]T∈N×1為接收陣列的導(dǎo)向矢量;at(φl)=[1,e-j2πdtsin(φl)/λ,…,e-j2π(M-1)dtsin(φl)/λ]T∈M×1為發(fā)射陣列的導(dǎo)向矢量;為均值為零的高斯白噪聲且與信號(hào)彼此獨(dú)立。

    X=AS+W

    (2)

    式中,A=[ar(θ1)?at(φ1),…,ar(θL)?at(φL)],?表示Kronecker積;S=[S(1),…,S(K)],S(k)=[s1,…,sL]T為在k時(shí)刻的目標(biāo)反射系數(shù)向量;W=[w(1),…,w(K)]為高斯白噪聲矩陣。

    在實(shí)際應(yīng)用中,隨著雙基地MIMO雷達(dá)天線陣列的不斷增多,以及受高低溫變化或振動(dòng)等較惡劣的外界環(huán)境和硬件老化等因素的影響,陣列中時(shí)常會(huì)出現(xiàn)陣元失效的情況。定義Q為失效接收陣元的位置集合,則陣元失效雙基地MIMO雷達(dá)的虛擬陣列輸出信號(hào)矩陣為

    (3)

    其中,diag(·)表示構(gòu)造對角矩陣;AT=[at(φ1),…,at(φL)]為發(fā)射陣列流形矩陣;0L×L為大小為L×L的全零矩陣。

    2 常用邊緣檢測算子

    在數(shù)字圖像中,邊緣是圖像的一項(xiàng)基本特征,是指圖像中局部特性的不連續(xù)導(dǎo)致灰度發(fā)生急劇變化的區(qū)域,邊緣主要存在于目標(biāo)與目標(biāo)和目標(biāo)與背景之間。圖像邊緣包含了大量的重要信息,利用圖像的邊緣能夠描繪出目標(biāo)的輪廓,因此邊緣檢測的研究獲得了眾多學(xué)者的關(guān)注。邊緣檢測的目的是找出圖像中灰度急劇變化的區(qū)域,圖像灰度的梯度能夠表示圖像中局部灰度的變化信息,圖像中(x,y)位置上的像素的圖像梯度可定義為

    (4)

    (1)Roberts算子

    Roberts邊緣檢測算子利用斜向偏差分來計(jì)算梯度,其方向垂直于圖像邊緣,其大小代表圖像邊緣的強(qiáng)度,計(jì)算時(shí)采用對角線方向相鄰兩像素之差,則Roberts的梯度幅度值為:

    G(i,j)=|f(i,j)-f(i+1,j+1)|+|f(i+1,j)-

    f(i,j+1)|

    (5)

    式(5)可用卷積模板來表示

    G(i,j)=|Gx|+|Gy|

    (6)

    (2)Sobel算子

    Sobel算子是從不同的方向檢測邊緣,采用的是中心像素上下左右鄰點(diǎn)灰度的加權(quán),該算子計(jì)算邊緣強(qiáng)度的公式表示如下:

    G(i,j)=|Gx|+|Gy|

    (7)

    (3)Laplace算子

    Laplace算子是一種二階微分算子,對于數(shù)字圖像f(x,y),Laplace算子在像素點(diǎn)位置(x,y)的定義為

    (8)

    3 基于Sobel算子的陣列診斷方法

    (9)

    其中,round(·)表示四舍五入取整;max(·)表示返回最大值。

    圖1為雙基地MIMO雷達(dá)虛擬陣列輸出信號(hào)的灰度圖像,其中雙基地MIMO雷達(dá)的發(fā)射陣元為5個(gè),接收陣元為15個(gè),其中第4、6、9、14個(gè)接收陣元失效。由圖1可以看出,正常接收陣元對應(yīng)的圖像塊矩陣具有較多圖像邊緣,而失效陣元對應(yīng)的圖像塊矩陣中存在較少圖像邊緣,因此正常接收陣元的圖像邊緣多于失效接收陣元的圖像邊緣。根據(jù)邊緣檢測理論可知,正常接收陣元的邊緣強(qiáng)度高于失效接收陣元的邊緣強(qiáng)度,因此可以將評(píng)價(jià)圖像中邊緣多少的邊緣檢測理論應(yīng)用于雙基地MIMO雷達(dá)的陣列診斷中。

    圖1 虛擬陣列輸出信號(hào)的灰度圖像

    在眾多邊緣檢測算子中,Sobel算子在漏檢邊緣、假檢邊緣、邊緣定位精度和抗圖像噪聲能力等方面具有明顯的優(yōu)勢[15]。本文選取Sobel邊緣檢測算子計(jì)算各接收陣元對應(yīng)圖像塊的邊緣強(qiáng)度,由式(7)可得圖像塊yn,n=1,…,N中的各像素的邊緣強(qiáng)度為

    Gn(i,j)=|Gx|+|Gy|,i=1,…,M,j=1,…,K

    (10)

    式中,Gx和Gy分別代表經(jīng)橫向及縱向邊緣檢測的圖像,(i,j)表示圖像塊yn中各像素的位置;Gx=yn(i-1,j+1)+2yn(i,j+1)+yn(i+1,j+1)-yn(i-1,j-1)-2yn(i,j-1)-yn(i+1,j-1);Gy=yn(i+1,j-1)+2yn(i+1,j)+yn(i+1,j+1)-yn(i-1,j-1)-2yn(i-1,j)-yn(i-1,j+1)。

    根據(jù)式(10),依次對各圖像塊yn進(jìn)行邊緣強(qiáng)度計(jì)算,從而獲得各圖像塊yn的邊緣強(qiáng)度總值,即

    (11)

    由于失效接收陣元所對應(yīng)的圖像塊的邊緣強(qiáng)度較低,而正常接收陣元所對應(yīng)的圖像塊的邊緣強(qiáng)度較高,因此通過設(shè)定合適的邊緣強(qiáng)度門限Gt=τ·mean(G)來診斷失效接收陣元的位置,其中τ為調(diào)節(jié)因子,mean(·)表示取均值,即

    (12)

    式中,Z(n)=1時(shí)表示第n個(gè)接收陣元失效,Z(n)=0時(shí)表示第n個(gè)接收陣元正常工作。

    4 仿真與分析

    為了驗(yàn)證本文方法在雙基地MIMO雷達(dá)陣列診斷方面的有效性,設(shè)計(jì)了幾組仿真實(shí)驗(yàn),將Zhu的方法[12]和本文方法應(yīng)用于陣列診斷時(shí)的性能進(jìn)行對比。在以下仿真中,雙基地MIMO雷達(dá)的發(fā)射陣元數(shù)M=5,接收陣元數(shù)N=15,收發(fā)陣元間距均為半波長,假設(shè)發(fā)射陣列各陣元發(fā)射相互正交的Hadamard編碼信號(hào),在每個(gè)重復(fù)周期內(nèi)的相位編碼個(gè)數(shù)為256??臻g中存在3個(gè)遠(yuǎn)場目標(biāo),方位角度分別(20°,8°),(-10°,23°),(0°,35°),回波噪聲選取均值為零的加性高斯白噪聲,回波信噪比定義為

    (13)

    定義失效陣元診斷的成功概率為

    (14)

    其中,MT為蒙特卡洛實(shí)驗(yàn)次數(shù);Tr為失效陣元診斷成功的次數(shù)。在本文方法中,調(diào)節(jié)因子τ=0.7。

    仿真實(shí)驗(yàn)1:診斷成功概率隨信噪比的變化關(guān)系

    假設(shè)接收陣列中存在4個(gè)位置隨機(jī)的陣元失效,快拍數(shù)為K=100,信噪比變化范圍為-30~0 dB,進(jìn)行100次蒙特卡洛實(shí)驗(yàn),陣元失效的情況分為兩種:(a)接收陣列中首尾陣元正常工作,每次實(shí)驗(yàn)的失效陣元在其余位置隨機(jī)選取;(b)每次實(shí)驗(yàn)假設(shè)首尾陣元至少存在一個(gè)陣元發(fā)生失效。由圖2(a)時(shí)的仿真結(jié)果可知,在陣列首尾處不出現(xiàn)陣元失效時(shí),隨著信噪比的不斷增加,本文方法和Zhu的方法診斷失效陣元的成功概率均逐漸提高,但本文方法在低信噪比下仍然能準(zhǔn)確診斷失效陣元的位置,因此其診斷性能優(yōu)于Zhu的方法。由于當(dāng)陣列中的失效陣元位置使同一差分位置的虛擬陣元同時(shí)失效時(shí),Zhu的方法會(huì)因匹配搜索失敗而失效,因此即使在高信噪比時(shí)Zhu的方法在個(gè)別實(shí)驗(yàn)中仍然會(huì)出現(xiàn)診斷失敗的情況。由圖2(b)的仿真結(jié)果可知,在陣列首尾存在失效陣元時(shí),Zhu的方法完全失效,而本文方法依然能夠準(zhǔn)確地診斷出失效陣元的位置,因此本文方法適用于任意位置下的失效陣元診斷問題,具有良好的穩(wěn)定性。

    圖2 診斷成功概率隨信噪比的變化關(guān)系

    仿真實(shí)驗(yàn)2: 診斷成功概率隨快拍數(shù)的變化關(guān)系

    圖3 診斷成功概率隨快拍數(shù)的變化關(guān)系

    假設(shè)接收陣列中存在4個(gè)隨機(jī)的陣元失效,在陣列首尾處不出現(xiàn)失效陣元以保證Zhu的方法工作正常,信噪比為-10 dB,令快拍數(shù)由20~100變化,進(jìn)行100次蒙特卡洛實(shí)驗(yàn),圖3為診斷成功概率隨快拍數(shù)的變化關(guān)系。由圖3可知,在不同的快拍數(shù)下,Zhu的方法始終在個(gè)別實(shí)驗(yàn)中存在診斷失敗的情況,而本文方法均能準(zhǔn)確的診斷出失效陣元位置,無診斷失敗的情況出現(xiàn),因此進(jìn)一步驗(yàn)證了本文方法的穩(wěn)定性。

    仿真實(shí)驗(yàn)3:診斷成功概率隨失效接收陣元數(shù)的變化關(guān)系

    為了驗(yàn)證本文方法在不同失效接收陣元數(shù)下的穩(wěn)健性,假設(shè)接收陣列中的隨機(jī)失效陣元數(shù)由1~11依次增加,在陣列首尾處不出現(xiàn)失效陣元以保證Zhu的方法工作正常,快拍數(shù)為K=100,信噪比為-10 dB,進(jìn)行100次蒙特卡洛實(shí)驗(yàn),仿真結(jié)果如圖4所示。由圖4可知,在失效接收陣元數(shù)不斷增加時(shí),Zhu的方法的診斷性能逐漸惡化甚至失效,而本文方法始終保持著100%的診斷成功概率,具有良好的穩(wěn)健性。

    圖4 診斷成功概率隨失效接收陣元數(shù)的變化關(guān)系

    5 結(jié) 語

    針對現(xiàn)有陣列診斷方法使用波導(dǎo)探頭會(huì)增加硬件成本并干擾待測場等問題,本文提出了一種無需波導(dǎo)探頭采樣的雙基地MIMO雷達(dá)陣列診斷方法。將雙基地MIMO雷達(dá)虛擬陣列輸出信號(hào)經(jīng)取模處理后的數(shù)據(jù)矩陣等效為灰度圖像數(shù)據(jù)矩陣,由于正常接收陣元對應(yīng)的圖像塊中目標(biāo)像素與背景像素互相重疊,而失效接收陣元對應(yīng)的圖像塊中均為背景,因此兩者的圖像邊緣強(qiáng)度存在明顯差異,因此利用Sobel算子計(jì)算各接收陣元對應(yīng)的圖像邊緣強(qiáng)度,從而實(shí)現(xiàn)雙基地MIMO雷達(dá)的失效接收陣元的診斷。理論分析與仿真實(shí)驗(yàn)表明,本文所提陣列診斷方法在無需使用波導(dǎo)探頭同時(shí)具有良好的診斷性能。

    猜你喜歡
    信噪比算子灰度
    采用改進(jìn)導(dǎo)重法的拓?fù)浣Y(jié)構(gòu)灰度單元過濾技術(shù)
    擬微分算子在Hp(ω)上的有界性
    基于灰度拉伸的圖像水位識(shí)別方法研究
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    基于深度學(xué)習(xí)的無人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    一類Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫
    低信噪比下LFMCW信號(hào)調(diào)頻參數(shù)估計(jì)
    電子測試(2018年11期)2018-06-26 05:56:02
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    基于最大加權(quán)投影求解的彩色圖像灰度化對比度保留算法
    Roper-Suffridge延拓算子與Loewner鏈
    一级av片app| 18禁在线播放成人免费| 欧美区成人在线视频| 国产欧美亚洲国产| 久久久精品免费免费高清| 高清在线视频一区二区三区| 91精品国产九色| 下体分泌物呈黄色| 精品国产国语对白av| 日日爽夜夜爽网站| 午夜免费鲁丝| 又粗又硬又长又爽又黄的视频| 色吧在线观看| 女人精品久久久久毛片| 69精品国产乱码久久久| 国产精品国产三级专区第一集| 在线天堂最新版资源| 亚洲国产日韩一区二区| 日日摸夜夜添夜夜添av毛片| 精品视频人人做人人爽| 久久狼人影院| 男人添女人高潮全过程视频| 中文在线观看免费www的网站| 在线亚洲精品国产二区图片欧美 | 国国产精品蜜臀av免费| 亚洲国产欧美在线一区| xxx大片免费视频| 18禁裸乳无遮挡动漫免费视频| 97超碰精品成人国产| 久久狼人影院| 午夜影院在线不卡| 美女中出高潮动态图| 久久久a久久爽久久v久久| 亚洲精品一区蜜桃| 国产乱人偷精品视频| 乱人伦中国视频| 纯流量卡能插随身wifi吗| 国产一级毛片在线| 一级毛片电影观看| 高清av免费在线| 久久久久视频综合| 亚洲av.av天堂| 又粗又硬又长又爽又黄的视频| 国产亚洲91精品色在线| xxx大片免费视频| 爱豆传媒免费全集在线观看| 亚洲综合色惰| 国产精品秋霞免费鲁丝片| 国产一区二区在线观看av| 日韩不卡一区二区三区视频在线| 国产熟女午夜一区二区三区 | 22中文网久久字幕| 寂寞人妻少妇视频99o| 欧美三级亚洲精品| 你懂的网址亚洲精品在线观看| av天堂中文字幕网| 日本免费在线观看一区| 免费观看在线日韩| 亚洲欧美一区二区三区黑人 | 草草在线视频免费看| 日韩伦理黄色片| 日韩人妻高清精品专区| 91久久精品电影网| 欧美精品一区二区大全| 日韩一本色道免费dvd| 亚洲精品久久久久久婷婷小说| 中文欧美无线码| 蜜臀久久99精品久久宅男| 2022亚洲国产成人精品| 亚洲av中文av极速乱| 亚洲美女黄色视频免费看| av免费观看日本| 看非洲黑人一级黄片| 最黄视频免费看| 日韩制服骚丝袜av| av免费在线看不卡| 最近最新中文字幕免费大全7| 日韩中文字幕视频在线看片| 爱豆传媒免费全集在线观看| 97超视频在线观看视频| 七月丁香在线播放| 又粗又硬又长又爽又黄的视频| 国产视频首页在线观看| 偷拍熟女少妇极品色| 久久综合国产亚洲精品| 噜噜噜噜噜久久久久久91| 免费观看性生交大片5| 啦啦啦中文免费视频观看日本| 亚洲av电影在线观看一区二区三区| 五月玫瑰六月丁香| 又粗又硬又长又爽又黄的视频| 黑人高潮一二区| 男女边吃奶边做爰视频| 国产精品一区二区在线不卡| 综合色丁香网| 啦啦啦视频在线资源免费观看| 免费观看无遮挡的男女| 精品国产国语对白av| 日韩免费高清中文字幕av| 男人爽女人下面视频在线观看| 下体分泌物呈黄色| 超碰97精品在线观看| 成年av动漫网址| 精品亚洲成a人片在线观看| 国产精品蜜桃在线观看| 熟妇人妻不卡中文字幕| 新久久久久国产一级毛片| 成人美女网站在线观看视频| 国产免费又黄又爽又色| 亚洲av成人精品一区久久| 久久精品国产亚洲av天美| 十八禁高潮呻吟视频 | 亚洲精品视频女| 国产精品无大码| 极品人妻少妇av视频| 日韩中字成人| www.色视频.com| 少妇丰满av| 天天躁夜夜躁狠狠久久av| 最新的欧美精品一区二区| 欧美精品一区二区大全| 亚洲av国产av综合av卡| 婷婷色综合大香蕉| 丰满人妻一区二区三区视频av| 一级片'在线观看视频| 国产高清国产精品国产三级| 激情五月婷婷亚洲| av国产精品久久久久影院| 王馨瑶露胸无遮挡在线观看| 亚洲欧美一区二区三区国产| 妹子高潮喷水视频| 亚洲精品国产av蜜桃| av黄色大香蕉| 国产一区有黄有色的免费视频| 在线精品无人区一区二区三| 亚洲电影在线观看av| 亚洲av中文av极速乱| 国产极品天堂在线| 国产伦精品一区二区三区视频9| 女人精品久久久久毛片| 搡老乐熟女国产| 18禁在线无遮挡免费观看视频| 丝袜脚勾引网站| 免费看日本二区| 丝袜喷水一区| 亚洲国产精品999| 97精品久久久久久久久久精品| 亚洲色图综合在线观看| 春色校园在线视频观看| 老司机影院成人| 亚洲国产成人一精品久久久| 亚洲国产毛片av蜜桃av| 另类亚洲欧美激情| 久久久久久久久久久免费av| 永久免费av网站大全| 69精品国产乱码久久久| 精品久久久久久久久亚洲| 夫妻性生交免费视频一级片| 国产 精品1| 国产高清不卡午夜福利| 丝瓜视频免费看黄片| 一级av片app| 久久99热这里只频精品6学生| 免费看光身美女| 高清毛片免费看| 亚洲国产毛片av蜜桃av| 高清视频免费观看一区二区| 欧美日韩一区二区视频在线观看视频在线| 在线观看免费日韩欧美大片 | 亚洲真实伦在线观看| 美女主播在线视频| 国产成人精品婷婷| 色视频在线一区二区三区| 一本大道久久a久久精品| 男人爽女人下面视频在线观看| 麻豆乱淫一区二区| 黑人巨大精品欧美一区二区蜜桃 | 91久久精品国产一区二区三区| 亚洲高清免费不卡视频| 国产老妇伦熟女老妇高清| 丰满饥渴人妻一区二区三| 国产在线一区二区三区精| 少妇的逼水好多| 国产黄片美女视频| 中文天堂在线官网| 99久久精品热视频| av一本久久久久| 丰满少妇做爰视频| 秋霞在线观看毛片| 丁香六月天网| 国产精品国产av在线观看| 久久久国产精品麻豆| 只有这里有精品99| 国产成人免费观看mmmm| 伦理电影免费视频| 一区二区三区免费毛片| 日韩免费高清中文字幕av| 国产精品久久久久久精品古装| 久久精品国产亚洲av天美| 午夜福利,免费看| 最近中文字幕2019免费版| 欧美精品一区二区免费开放| 国产男人的电影天堂91| 免费观看无遮挡的男女| 黄色视频在线播放观看不卡| 亚洲美女搞黄在线观看| kizo精华| 精品国产露脸久久av麻豆| 国产精品福利在线免费观看| 亚洲一级一片aⅴ在线观看| 亚洲情色 制服丝袜| 亚洲精品国产av成人精品| 国产乱来视频区| 成人美女网站在线观看视频| 亚洲性久久影院| 国产免费一级a男人的天堂| 免费观看a级毛片全部| 久久久久人妻精品一区果冻| 性色avwww在线观看| 国产91av在线免费观看| 午夜福利,免费看| 国产成人精品一,二区| 精品少妇黑人巨大在线播放| 亚洲第一区二区三区不卡| 观看av在线不卡| 日韩av免费高清视频| 91久久精品国产一区二区成人| 久久人人爽人人片av| 少妇精品久久久久久久| 国产成人freesex在线| 22中文网久久字幕| 久久久久久伊人网av| 青青草视频在线视频观看| 51国产日韩欧美| 观看免费一级毛片| 亚洲av二区三区四区| 亚洲天堂av无毛| 亚洲精品国产色婷婷电影| 久久久久网色| 乱人伦中国视频| 看十八女毛片水多多多| 亚洲精品,欧美精品| 亚洲国产欧美在线一区| 一本久久精品| 成年av动漫网址| 99国产精品免费福利视频| 日本-黄色视频高清免费观看| 蜜臀久久99精品久久宅男| 一边亲一边摸免费视频| 国产精品熟女久久久久浪| 新久久久久国产一级毛片| 久久女婷五月综合色啪小说| 免费观看a级毛片全部| a级毛片免费高清观看在线播放| 春色校园在线视频观看| 亚洲国产成人一精品久久久| 亚洲,欧美,日韩| 只有这里有精品99| 亚洲欧洲日产国产| 春色校园在线视频观看| 亚洲精品国产av蜜桃| 久久精品国产亚洲av涩爱| 国产精品成人在线| 伦理电影免费视频| 在现免费观看毛片| 97在线人人人人妻| 看十八女毛片水多多多| av播播在线观看一区| 黑人猛操日本美女一级片| 欧美日韩国产mv在线观看视频| 如何舔出高潮| 伊人久久精品亚洲午夜| 中文欧美无线码| 少妇裸体淫交视频免费看高清| 美女视频免费永久观看网站| 久久久国产欧美日韩av| 亚洲国产日韩一区二区| 亚洲美女黄色视频免费看| 成人特级av手机在线观看| 大香蕉97超碰在线| 国产高清三级在线| 精品国产露脸久久av麻豆| 亚洲av.av天堂| 亚洲婷婷狠狠爱综合网| 内地一区二区视频在线| 成人免费观看视频高清| 亚洲,一卡二卡三卡| 夜夜骑夜夜射夜夜干| av线在线观看网站| 成人亚洲精品一区在线观看| 国产极品粉嫩免费观看在线 | 欧美日韩在线观看h| 成人18禁高潮啪啪吃奶动态图 | 青春草亚洲视频在线观看| 国产69精品久久久久777片| 尾随美女入室| 亚洲精品久久久久久婷婷小说| 国产精品一区二区在线不卡| 超碰97精品在线观看| 国产白丝娇喘喷水9色精品| 一区二区av电影网| 久久青草综合色| 丝瓜视频免费看黄片| 九九在线视频观看精品| 亚洲精品乱久久久久久| 亚洲综合精品二区| 亚洲av成人精品一区久久| 国产成人精品无人区| 免费久久久久久久精品成人欧美视频 | 91aial.com中文字幕在线观看| av天堂久久9| 亚洲欧美日韩另类电影网站| 蜜桃在线观看..| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| 精品人妻熟女av久视频| 日韩欧美一区视频在线观看 | 精品久久久精品久久久| 中文在线观看免费www的网站| 久久久久久久久久人人人人人人| 菩萨蛮人人尽说江南好唐韦庄| 亚洲综合精品二区| 一个人免费看片子| 中文天堂在线官网| 91久久精品国产一区二区三区| 亚洲精品中文字幕在线视频 | 久久久午夜欧美精品| 一级二级三级毛片免费看| 亚洲av综合色区一区| 人妻制服诱惑在线中文字幕| 久久精品国产自在天天线| 国产成人a∨麻豆精品| 国产日韩欧美亚洲二区| 欧美日韩一区二区视频在线观看视频在线| 成人国产av品久久久| av国产久精品久网站免费入址| 内地一区二区视频在线| 欧美日韩综合久久久久久| 黑人高潮一二区| 亚洲欧美一区二区三区国产| 国产精品不卡视频一区二区| 国产黄频视频在线观看| 亚洲,一卡二卡三卡| 亚洲国产精品国产精品| 国产欧美亚洲国产| 久久久久久久久大av| 中文乱码字字幕精品一区二区三区| 另类亚洲欧美激情| 一区在线观看完整版| 国模一区二区三区四区视频| 久久久a久久爽久久v久久| 国产女主播在线喷水免费视频网站| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 伊人久久精品亚洲午夜| 日韩制服骚丝袜av| 国产精品嫩草影院av在线观看| 午夜免费观看性视频| 日本av免费视频播放| 日日摸夜夜添夜夜爱| 少妇精品久久久久久久| 在线观看av片永久免费下载| 日本av手机在线免费观看| av线在线观看网站| 高清黄色对白视频在线免费看 | 免费少妇av软件| 亚洲精品亚洲一区二区| 午夜免费鲁丝| 亚洲国产精品国产精品| 欧美精品一区二区大全| 免费大片18禁| 伊人久久国产一区二区| 两个人的视频大全免费| 午夜久久久在线观看| 国产一级毛片在线| 中文字幕人妻丝袜制服| 一级a做视频免费观看| 国产无遮挡羞羞视频在线观看| 黑丝袜美女国产一区| 亚洲内射少妇av| 亚洲天堂av无毛| 一区二区三区免费毛片| 午夜影院在线不卡| 久久国产精品大桥未久av | 精品人妻一区二区三区麻豆| 国产女主播在线喷水免费视频网站| 久久人人爽av亚洲精品天堂| 91成人精品电影| 伊人亚洲综合成人网| 卡戴珊不雅视频在线播放| 亚洲图色成人| 国产免费福利视频在线观看| 亚洲欧洲日产国产| 女的被弄到高潮叫床怎么办| 夫妻性生交免费视频一级片| 在线观看免费日韩欧美大片 | 桃花免费在线播放| 老司机影院成人| 欧美高清成人免费视频www| 日韩不卡一区二区三区视频在线| 国产精品久久久久成人av| 少妇被粗大的猛进出69影院 | 黄色毛片三级朝国网站 | av免费在线看不卡| 老司机影院成人| 久久久欧美国产精品| 又大又黄又爽视频免费| 最近手机中文字幕大全| 精品一品国产午夜福利视频| 精华霜和精华液先用哪个| 日韩亚洲欧美综合| 亚洲无线观看免费| 我的女老师完整版在线观看| 国产黄色视频一区二区在线观看| 国产在线视频一区二区| 久久久久久久久大av| 国产日韩欧美视频二区| 精品国产一区二区三区久久久樱花| 国产精品欧美亚洲77777| 亚洲av不卡在线观看| 大话2 男鬼变身卡| 久久青草综合色| 久久精品国产a三级三级三级| 边亲边吃奶的免费视频| 天美传媒精品一区二区| 国产成人精品一,二区| 国产欧美日韩综合在线一区二区 | 中文字幕人妻熟人妻熟丝袜美| 成年美女黄网站色视频大全免费 | 国精品久久久久久国模美| 中文字幕久久专区| 看非洲黑人一级黄片| 在现免费观看毛片| 内射极品少妇av片p| 日韩三级伦理在线观看| 中文字幕亚洲精品专区| 亚洲一区二区三区欧美精品| 亚洲成人手机| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 18禁动态无遮挡网站| 中文天堂在线官网| 亚洲欧美日韩卡通动漫| 在线观看人妻少妇| av福利片在线| 97精品久久久久久久久久精品| 高清午夜精品一区二区三区| 嫩草影院新地址| 成人午夜精彩视频在线观看| 十分钟在线观看高清视频www | 人人妻人人添人人爽欧美一区卜| 免费黄频网站在线观看国产| 自拍偷自拍亚洲精品老妇| 亚洲va在线va天堂va国产| 99re6热这里在线精品视频| 国产男人的电影天堂91| 极品教师在线视频| 国产女主播在线喷水免费视频网站| 日日爽夜夜爽网站| 99久久精品一区二区三区| 亚洲综合色惰| 妹子高潮喷水视频| 肉色欧美久久久久久久蜜桃| 嫩草影院新地址| 久久毛片免费看一区二区三区| 日韩一本色道免费dvd| 一区二区三区乱码不卡18| 蜜桃在线观看..| 久久久久视频综合| 久久久久久久国产电影| 日韩人妻高清精品专区| 一区二区三区免费毛片| 91精品一卡2卡3卡4卡| 大片电影免费在线观看免费| 中文字幕制服av| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 久久鲁丝午夜福利片| 免费av中文字幕在线| 久久精品国产a三级三级三级| 国产精品成人在线| 97精品久久久久久久久久精品| 99精国产麻豆久久婷婷| 久久人人爽av亚洲精品天堂| 免费少妇av软件| 日韩不卡一区二区三区视频在线| 久久精品国产亚洲网站| 永久网站在线| 国产成人免费无遮挡视频| 亚洲精华国产精华液的使用体验| av又黄又爽大尺度在线免费看| 欧美成人午夜免费资源| 欧美日韩视频精品一区| 亚洲国产精品成人久久小说| 大陆偷拍与自拍| 高清不卡的av网站| 80岁老熟妇乱子伦牲交| 国产成人91sexporn| 性色av一级| 91午夜精品亚洲一区二区三区| 性色avwww在线观看| av有码第一页| 一边亲一边摸免费视频| 国产高清三级在线| 精品久久久久久电影网| 亚洲国产欧美日韩在线播放 | 卡戴珊不雅视频在线播放| 国产精品人妻久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 免费高清在线观看视频在线观看| 日本欧美国产在线视频| 人体艺术视频欧美日本| 久久韩国三级中文字幕| 看免费成人av毛片| 国产黄片视频在线免费观看| 18+在线观看网站| 在线观看三级黄色| 观看免费一级毛片| 国产欧美日韩综合在线一区二区 | 国产午夜精品一二区理论片| 美女主播在线视频| 秋霞在线观看毛片| 国产精品国产三级专区第一集| 狂野欧美激情性bbbbbb| 国产熟女欧美一区二区| 午夜91福利影院| 国产精品.久久久| 青青草视频在线视频观看| 尾随美女入室| 成人18禁高潮啪啪吃奶动态图 | 建设人人有责人人尽责人人享有的| 午夜福利网站1000一区二区三区| 国产老妇伦熟女老妇高清| 一级毛片黄色毛片免费观看视频| 18禁在线无遮挡免费观看视频| 国产在线免费精品| av国产精品久久久久影院| 看免费成人av毛片| 人人妻人人添人人爽欧美一区卜| 蜜桃在线观看..| 国产视频首页在线观看| 一级毛片aaaaaa免费看小| 一本—道久久a久久精品蜜桃钙片| 日本91视频免费播放| 九九久久精品国产亚洲av麻豆| 自线自在国产av| 日韩伦理黄色片| 99热这里只有精品一区| 日本av免费视频播放| 男人狂女人下面高潮的视频| 日韩av不卡免费在线播放| 亚洲欧美日韩另类电影网站| 亚洲第一区二区三区不卡| 日本欧美国产在线视频| 天天操日日干夜夜撸| 乱系列少妇在线播放| 久久免费观看电影| 午夜激情福利司机影院| 婷婷色麻豆天堂久久| 少妇人妻精品综合一区二区| freevideosex欧美| 一级黄片播放器| 只有这里有精品99| 男人狂女人下面高潮的视频| 精华霜和精华液先用哪个| 久久久久久久国产电影| 欧美日韩精品成人综合77777| 熟妇人妻不卡中文字幕| 欧美bdsm另类| 熟妇人妻不卡中文字幕| 美女国产视频在线观看| 久久久久久久亚洲中文字幕| 九色成人免费人妻av| 少妇被粗大的猛进出69影院 | 黑人巨大精品欧美一区二区蜜桃 | 十八禁高潮呻吟视频 | 亚洲精品一区蜜桃| 日本vs欧美在线观看视频 | 美女国产视频在线观看| 又黄又爽又刺激的免费视频.| 中文精品一卡2卡3卡4更新| 男女无遮挡免费网站观看| 国产男人的电影天堂91| h视频一区二区三区| 亚洲精品国产成人久久av| 亚洲欧美成人精品一区二区| 校园人妻丝袜中文字幕| 性高湖久久久久久久久免费观看| 亚洲成人手机| 肉色欧美久久久久久久蜜桃| 亚洲欧美中文字幕日韩二区| 国内少妇人妻偷人精品xxx网站| 熟女人妻精品中文字幕| 免费高清在线观看视频在线观看| 另类精品久久| 亚洲欧美成人综合另类久久久| av在线app专区| av又黄又爽大尺度在线免费看| 亚洲丝袜综合中文字幕| 成人特级av手机在线观看| 国产精品99久久久久久久久| 一二三四中文在线观看免费高清| 夫妻性生交免费视频一级片| 欧美 亚洲 国产 日韩一| 天天躁夜夜躁狠狠久久av| 人妻系列 视频| 精品久久久久久久久亚洲| 国产乱人偷精品视频| 亚洲国产欧美在线一区| 另类亚洲欧美激情| 男女边吃奶边做爰视频| 伦精品一区二区三区| 六月丁香七月| 寂寞人妻少妇视频99o| 精品一区二区三区视频在线|