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

    基于局部散射中心的近、遠場微動回波時頻分布特性的解析表達

    2022-08-19 02:54:40占偉杰萬顯榮易建新
    電子與信息學報 2022年8期
    關(guān)鍵詞:扇葉遠場微動

    占偉杰 萬顯榮 易建新

    (武漢大學電子信息學院 武漢 430072)

    1 引言

    微多普勒效應(yīng)是由目標(或其部件)的轉(zhuǎn)動、振動、進動等微動引起的頻率調(diào)制現(xiàn)象[1,2]。它的特征表現(xiàn)與目標的幾何結(jié)構(gòu)和運動狀態(tài)密切相關(guān),是實現(xiàn)目標精細化建模的重要輔助信息,常被應(yīng)用于雷達目標狀態(tài)監(jiān)測、分類識別等領(lǐng)域[3–8]。微多普勒效應(yīng)源于目標微動,具有時變特性,因此對目標微動回波的時頻分布特性進行深入分析,有利于進一步挖掘它與目標參數(shù)之間的關(guān)系,拓展其應(yīng)用范圍。

    微動回波時頻分布特性的有效獲取與精確表示是微動研究領(lǐng)域的熱點。目前常利用短時傅里葉變換(Short-Time Fourier Transform, STFT)、短時分數(shù)階傅里葉變換(Short-Time Fractional Order Fourier Transform, STFRFT)、魏格納維爾分布(Wigner-Ville Distribution, WVD)、壓縮感知類等時頻變換方法獲取微動回波的時頻分布特性結(jié)果[9–11]。上述研究的目的是尋找能夠兼顧時間維、頻率維精度且交叉項小、聚焦性好的時頻分析方法,為后續(xù)微動回波時頻分布特性的深入分析提供有益的處理方法和工具。

    當微動部件的尺寸遠小于它與雷達之間的距離時,滿足遠場探測條件,入射波可近似看成平面波[12]。直升機、通航飛機、無人機等目標的旋翼尺寸較小,往往滿足遠場探測條件。目前已有大量研究展示了上述目標的微動回波時頻分布結(jié)果[13–17]。這些研究表明遠場微動回波的時頻圖中包含多種不同形狀的強能量條帶(稱為flash),且flash的位置、形狀等特征與扇葉數(shù)目、長度、轉(zhuǎn)速等參數(shù)存在一定關(guān)系。文獻[18]基于散射點散射系數(shù)和散射點在葉片上的分布情況,構(gòu)建了葉片回波的散射點模型,定性地解釋了上述flash的形成機理。

    近年來,基于微多普勒效應(yīng)的風電機組健康狀態(tài)監(jiān)測逐漸引起關(guān)注[19,20]。這些研究均是通過建立葉片微動回波的時頻分布特性與葉片狀態(tài)之間的關(guān)系來達到監(jiān)測目的的。風電機組葉片的尺寸較大,容易落入雷達的近場探測區(qū),此時其微動回波的時頻分布特性也呈現(xiàn)出更復(fù)雜的形態(tài)。因此有必要對近場微動回波的時頻分布特性展開研究。已有部分研究展示了近場風電機組葉片微動回波的時頻分布結(jié)果[21,22]。從這些結(jié)果中可以看出近場微動回波時頻圖中的矩形flash已不再垂直于時間軸,而是呈現(xiàn)傾斜甚至彎曲的形態(tài)。文獻[23]從電磁散射的角度對上述特殊形式的flash進行了初步解釋。進一步地,文獻[24]建立了風電機組葉片近場微動回波模型,分析了不同形狀葉片的微動回波特性,并通過風電機組模型進行實驗驗證。

    目前,已有較多近、遠場微動回波時頻分布特性結(jié)果的報道,為后續(xù)深入研究提供了有益的借鑒。不過,目前鮮有研究給出近、遠場微動回波時頻圖中各flash形狀的具體表達式,及其與扇葉、探測場景等各參數(shù)之間的對應(yīng)關(guān)系。為進一步實現(xiàn)微動回波時頻分布特性的精細化分析,以更好服務(wù)于后續(xù)深層次應(yīng)用,本文主要完成了以下工作:

    (1)從遠場微動回波信號模型中推導得到了微動回波瞬時頻率表達式,其各部分分別對應(yīng)時頻圖中的正弦型flash、零頻flash及矩形flash,并從積分運算性質(zhì)和電磁散射理論兩方面解釋了上述flash的形成機理。

    (2)推導了近場條件下,葉尖散射點、葉彀散射點及鏡面反射點的瞬時頻率表達式,證明了近場微動回波時頻圖中包含由上述3類散射點引入的類正弦型flash、零頻flash及部分余弦型flash。

    (3)給出了近、遠場微動回波中各種flash的位置分布、頻率拓展情況等特征與扇葉尺寸、轉(zhuǎn)速、數(shù)目及探測場景參數(shù)之間的關(guān)系表達式,并用仿真和實測數(shù)據(jù)驗證了上述分析結(jié)果。

    本文第2節(jié)建立了近場和遠場探測條件下扇葉轉(zhuǎn)動引入的微動回波信號模型。第3節(jié)對近、遠場探測條件下微動回波的時頻分布特性進行了詳細分析。第4節(jié)利用仿真和實測數(shù)據(jù)驗證了分析結(jié)果的正確性。最后對全文進行總結(jié)。

    2 旋轉(zhuǎn)葉片微動回波模型

    為簡潔直觀考慮,如圖1所示,以2維平面上單基地雷達為例進行分析。坐標原點Q表示扇葉旋轉(zhuǎn)中心,O表示雷達。φ(t)=ωt+φ0表示t時刻扇葉與X軸的夾角,φ0為初始位置扇葉與X軸的夾角。ω為扇葉轉(zhuǎn)速,本文中扇葉均按逆時針旋轉(zhuǎn)。|OQ|=R0表 示雷達與扇葉旋轉(zhuǎn)中心的距離,L為扇葉長度。當R0?L時,雷達處于遠場探測模式,此時扇葉上點P到雷達的距離可以近似為圖中紫色虛線,稱為遠場距離RF(t)。 當不滿足R0?L時,雷達處于近場探測模式,此時點P到雷達的距離稱為近場距離RN(t), 用圖中紅色虛線表示。設(shè)|PQ|=l,則RF(t)和RN(t)可以分別表示為

    圖1 近、遠場探測場景示意圖

    3 微動回波時頻分布特性

    3.1 遠場微動回波的時頻分布特性

    圖2 兩個信號乘積、加和的時頻分布結(jié)果

    f(t)=1/cos(t)

    此時I F[sF(t)]=0∪fmaxsin[φk(t)],時頻圖表現(xiàn)為零頻flash和K個正弦flash的組合,且相鄰正弦flash之間的相位差為2 π/Krad。

    綜上所述,與偶數(shù)扇葉微動回波時頻圖相比,奇數(shù)扇葉微動回波時頻圖中不僅包含矩形flash和正弦flash,還包含零頻flash。矩形flash和正弦flash在正負頻率面上呈現(xiàn)交錯分布形式。

    3.1.3 扇葉回波時頻分布特性的形成機制

    從上面的分析可知,當t→tm,k時第k個扇葉上所有散射點到雷達的距離相等,它們的回波可實現(xiàn)同相疊加,回波能量達到最大,在時頻圖中表現(xiàn)為垂直于零頻線的矩形flash。當不滿足t →tm,k時,各個散射點的回波互相抵消,整體回波由式(2)中的定積分運算確定。而定積分運算結(jié)果由積分表達式和積分限決定,因此最終扇葉回波的性質(zhì)將由積分上下限位置處(即葉尖散射點和葉轂散射點)的散射點決定。其中葉尖散射點回波對應(yīng)時頻圖中的正弦flash,葉轂散射點回波對應(yīng)時頻圖中的零頻flash。

    從電磁散射理論的角度分析,復(fù)雜目標的電磁散射來自局部效應(yīng)的散射中心,如鏡面反射中心、尖頂散射中心、蠕動波和行波效應(yīng)散射中心等[27]。當t→tm,k時,扇葉上所有散射點的入射波方向(雷達視線方向)均與扇葉垂直,進而產(chǎn)生鏡面反射,故此時扇葉上所有散射點均為鏡面反射中心,它們的回波一起在時頻圖中形成強的矩形flash。葉尖、葉彀散射點位于扇葉表面不連續(xù)處,構(gòu)成扇葉的尖頂散射中心,其中葉尖散射點回波在時頻圖中表現(xiàn)為正弦型flash,葉彀散射點回波在時頻圖中則表現(xiàn)為零頻flash。

    3.2 近場微動回波的時頻分布特性

    近場條件下,R0和L處于同一數(shù)量級,難以對式(2)中sN(t)進行近似化簡。從3.1節(jié)的分析可知,扇葉微動回波的時頻分布特性主要由鏡面反射點,葉尖、葉彀散射點等幾類局部散射中心決定。因此本節(jié)對這些散射中心進行分析以獲取近場微動回波的時頻分布特性。

    3.2.1 葉尖散射點引入的類正弦型flash

    由式(2)可知,葉尖散射點回波可以表示為

    圖4 近場探測條件下,葉尖散射點情況

    綜上所述,近場條件下偶數(shù)扇葉微動回波的時頻圖包含類正弦型flash和部分余弦型flash,且在正負頻率面上呈現(xiàn)斜對稱的分布形態(tài)。而奇數(shù)扇葉微動回波時頻圖則由正弦型flash、部分余弦型flash和零頻flash組成,且在在正負頻率面上呈現(xiàn)斜交錯分布的形態(tài)。

    4 仿真結(jié)果

    4.1 遠場探測情況

    圖5 近場探測條件下,鏡面反射點情況

    設(shè)信號中心頻率f0=722 MHz ,雷達與扇葉旋轉(zhuǎn)中心的距離R0=20000 m ,扇葉長度L=5 m,轉(zhuǎn) 速ω=120 r/min ,扇 葉 數(shù) 目 取K=2 和K=3兩種情況,第1個扇葉的初始相位φ0=0°。根據(jù)上述參數(shù)可得旋轉(zhuǎn)周期T=0.5 s, 信號波長λ=0.4155 m,最大多普勒頻率fmax≈302.43 Hz,仿真時長設(shè)置為0.5 s,對應(yīng)一個完整旋轉(zhuǎn)周期。得到圖6所示結(jié)果。

    根據(jù)3.1.1節(jié)分析結(jié)論可知在上述仿真條件下,2扇葉回波在一個旋轉(zhuǎn)周期內(nèi)的時頻圖中包含2個頻率范圍為–302.43 ~302.43 Hz的矩形flash,且各個矩形flash對應(yīng)的時間為kT/4,k ∈[1,2],2個相差為π,最大頻率為302.43 Hz的正弦flash。圖6(b)所示仿真時頻結(jié)果與圖6(a)所示理論時頻結(jié)果吻合,驗證了上述分析結(jié)論。

    圖6中矩形flash所在時間點回波能量處于sinc函數(shù)峰值位置,隨著時間偏移峰值時間點,回波能量急劇下降,因此圖中正弦flash的強度遠小于矩形flash,與理論分析結(jié)果一致。

    圖6 遠場條件下扇葉回波時頻分布特性仿真結(jié)果

    4.2 近場探測情況

    設(shè)信號中心頻率f0=722 MHz,雷達到扇葉旋轉(zhuǎn)中心的距離取R0=90 m 和R0=30 m兩種情況,扇葉長度L=60 m ,轉(zhuǎn)速ω=12 r/min,葉片數(shù)目取K=2和K=3兩種情況,第1個扇葉的初始相位φ0=0°。 根據(jù)上述參數(shù)可得旋轉(zhuǎn)周期T=5 s,信號波長λ=0.4155 m, 最大多普勒頻率fmax≈362.92 Hz,仿真時長設(shè)置為5 s,對應(yīng)一個完整旋轉(zhuǎn)周期。

    根據(jù)3.2節(jié)分析結(jié)果可知在一個周期內(nèi),每個扇葉均會引入1個正弦型flash和2個部分余弦型flash,且相鄰扇葉引入的同一類型flash之間的相位差為 (2π/K)rad 。當R0=90 m ,L=60 m時,正弦型flash和部分余弦型flash均能取得最大頻率絕對值fmax,且兩者取到最大頻率絕對值的時刻相同。此時每個扇葉引入的部分余弦flash的持續(xù)時間約為0.5807 s(約0.1161個周期),且當扇葉數(shù)目為偶數(shù)時,正負頻率面上相鄰的兩個部分余弦型flash在零頻處相接,組合成一個持續(xù)時間約為1.1614 s(約為0.2322個周期)的部分余弦型flash。當扇葉數(shù)目為奇數(shù)時,正負頻率面上的部分余弦型flash呈現(xiàn)交錯分布的狀態(tài),且此時存在零頻flash。R0>L條件下的理論和仿真時頻分布結(jié)果分別如圖7(a)~圖7(d)所示,與上述分析結(jié)論吻合。

    圖7 近場條件下扇葉回波時頻分布特性仿真結(jié)果

    當R0=30 m ,L=60 m時,能取到的最大頻率絕對值為fmax/2≈181.41 Hz,且正弦型flash和部分余弦型flash當頻率絕對值取得最大時對應(yīng)的時間點不同。此時,每個部分余弦型flash對應(yīng)的持續(xù)時間均為1/4個周期,即1.25 s。除上述特性外,R0<L情況下各種flash的數(shù)目及其分布與扇葉數(shù)目奇偶性之間的關(guān)系均與R0>L情 況下類似。R0<L條件下理論和仿真時頻分布結(jié)果分別如圖7(e)~圖7(h)所示,與上述分析結(jié)論吻合。

    5 實測結(jié)果

    本文實測數(shù)據(jù)均來自武漢大學自主研制的UHF頻段外輻射源雷達系統(tǒng)[28]。其中直升機和無人機旋翼的回波數(shù)據(jù)用于驗證3.1節(jié)中遠場探測條件下的時頻分布特性。風電機組扇葉回波數(shù)據(jù)則用于驗證3.2節(jié)中近場探測條件下的時頻分布特征。值得注意的是,外輻射源雷達是一種雙基地雷達,其微動回波受收發(fā)站及目標之間的位置關(guān)系影響(主要涉及雙基地角β、雙基角平分線與旋轉(zhuǎn)平面的夾角δ),扇葉上散射點的雷達視線方向由雙基角平分線決定[29]。

    5.1 遠場探測情況

    遠場探測條件下相關(guān)實驗配置和目標參數(shù)如表1所示。

    表1 實驗配置及目標參數(shù)

    由所示參數(shù),可以計算得到直升機旋翼的旋轉(zhuǎn)周期T=0.5 s ,最大多普勒頻率fmax≈850.86 Hz,相鄰矩形flash之間的時間間隔約為24.60 ms。圖8(a)給出了直升機2個完整旋轉(zhuǎn)周期內(nèi)的時頻分布結(jié)果,可以看到圖中包含12個沿正負多普勒平面交錯分布的矩形flash,且各矩形flash的相關(guān)參數(shù)與理論分析結(jié)果一致。圖中零頻線及其周圍分布著極強的信號成分,其中混合著直升機機身回波分量及微多普勒回波本身的零頻flash分量。值得注意的是,圖8(a)并未見到明顯的正弦flash,這是因為正弦型flash的能量遠小于矩形flash,在實測環(huán)境中可能被噪聲湮沒。

    由所示參數(shù),可以計算得到無人機旋翼的旋轉(zhuǎn)周期T=23.5 ms, 最大多普勒頻率fmax≈177.57 Hz,相鄰矩形flash之間的時間間隔約為11.76 ms。圖8(b)給出了無人機2個完整旋轉(zhuǎn)周期內(nèi)的時頻分布結(jié)果,可以看到圖中包含8個沿正負多普勒平面對稱分布的矩形flash,且各矩形flash的相關(guān)參數(shù)與理論分析結(jié)果一致。需要注意的是,圖中矩形flash的數(shù)目是8個,而不是理論計算的4個,這是因為實驗所用無人機為四旋翼無人機,包含4幅兩葉片的旋翼,其中兩幅葉片同步正轉(zhuǎn)、兩幅葉片同步反轉(zhuǎn),所以圖8(b)同時包含了正轉(zhuǎn)和反轉(zhuǎn)的旋翼微動回波。可以看到圖中不存在明顯的正弦flash,這是因為正弦flash的能量較微弱,被噪聲湮沒。

    圖8 遠場條件下實測時頻分布特性

    5.2 近場探測情況

    實驗所用信號頻率為722 MHz,風電機組扇葉長度為56.8 m,旋轉(zhuǎn)中心距地面高度約為80 m。機組與發(fā)射站之間的距離約為20.06 km,與接收站之間的距離約為24 km。數(shù)據(jù)采集期間,扇葉轉(zhuǎn)速約為12.70 r/min,扇面旋轉(zhuǎn)軸朝向約為南偏西10°(近似認為扇葉在垂直于水平面的平面上轉(zhuǎn)動)。根據(jù)上述參數(shù),可計算得到如圖9(a)所示的理論時頻分布特性??梢钥闯觯c單基雷達中,近場微動回波時頻分布規(guī)律一樣,圖9(a)包含多種不同形態(tài)的交錯分布的flash。由于雙基雷達中,目標回波的多普勒特性還受收發(fā)站與目標之間的位置關(guān)系的影響,因此圖9(a)呈現(xiàn)出一些更復(fù)雜的表現(xiàn),如:正負頻率面上多普勒拓展范圍不一樣(此時負頻率面上多普勒拓展范圍大于正頻率面上的多普勒拓展范圍);由鏡面反射點引入的flash在正負頻率面上的形態(tài)存在差異(此時正頻率面上鏡面反射flash彎曲程度大于負頻率面上的鏡面反射flash)。

    圖9(b)為實測風電機組扇葉微動回波時頻圖,將其與圖9(a)理論分析結(jié)果對比,兩者表現(xiàn)基本一致,驗證了理論分析結(jié)果。值得注意的是,理論分析模型中假設(shè)扇葉為理想線型、同時扇葉上所有散射點的散射強度保持一致。而實際情況中,扇葉形狀、各散射點的散射強度等因素均會對其回波特性造成影響,因此,扇葉實測回波數(shù)據(jù)的時頻分析結(jié)果將呈現(xiàn)更加復(fù)雜的細節(jié)表現(xiàn)。

    圖9 近場條件下實測風電機組回波的時頻分布特性

    6 結(jié)束語

    本文系統(tǒng)分析了近、遠場探測條件下,扇葉微動回波的時頻分布特性。結(jié)果表明,兩種條件下,時頻圖均表現(xiàn)為多種不同類型flash的組合。文中對這些flash的具體表現(xiàn)形式及其形成機理進行了詳細說明,并推導了上述flash的位置、分布、形態(tài)等特征與扇葉數(shù)目、尺寸、轉(zhuǎn)速等參數(shù)之間的關(guān)系。為目標微多普勒效應(yīng)的精細化建模,并服務(wù)于后續(xù)目標參數(shù)估計及分類識別等深層次應(yīng)用打下基礎(chǔ)。文中結(jié)論是在將扇葉簡化為直線型目標的前提下得到的,接下來將針對不同扇葉形狀對時頻特征的影響展開研究。

    猜你喜歡
    扇葉遠場微動
    基于RID序列的微動目標高分辨三維成像方法
    雷達學報(2018年5期)2018-12-05 03:13:08
    基于仿真與實測的列車遠場氣動噪聲分析
    基于稀疏時頻分解的空中目標微動特征分析
    互動小游戲
    大灰狼(2017年6期)2017-07-12 20:35:02
    風的力量
    某種陣列雷達發(fā)射通道遠場校準簡易方法
    微動橋橋足距離對微動裂紋萌生特性的影響
    沒有扇葉的風扇
    戰(zhàn)斗部遠場水下爆炸對艦船沖擊損傷評估
    遠場天線測試系統(tǒng)的研究與實現(xiàn)
    十八禁网站免费在线| 国产黄色免费在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品在线观看二区| 国产精品久久久久久精品电影小说| 午夜激情久久久久久久| 国产精品偷伦视频观看了| 亚洲天堂av无毛| 露出奶头的视频| www.999成人在线观看| 日韩欧美三级三区| 亚洲精品一二三| 亚洲精品一二三| 日韩视频在线欧美| 午夜精品国产一区二区电影| 一级毛片女人18水好多| 国产精品电影一区二区三区 | 国产黄色免费在线视频| 亚洲av第一区精品v没综合| 亚洲av第一区精品v没综合| 高清黄色对白视频在线免费看| 少妇裸体淫交视频免费看高清 | 精品亚洲成国产av| 亚洲专区中文字幕在线| 欧美+亚洲+日韩+国产| 一区福利在线观看| 搡老岳熟女国产| 亚洲第一av免费看| 欧美黑人精品巨大| 在线永久观看黄色视频| 精品国产超薄肉色丝袜足j| 久久婷婷成人综合色麻豆| 丝袜在线中文字幕| 十八禁高潮呻吟视频| 欧美av亚洲av综合av国产av| 日韩欧美国产一区二区入口| 天天躁狠狠躁夜夜躁狠狠躁| av福利片在线| 亚洲精品国产精品久久久不卡| 黄片小视频在线播放| 美女高潮到喷水免费观看| 亚洲人成77777在线视频| 国产福利在线免费观看视频| 大香蕉久久成人网| 最黄视频免费看| 亚洲国产av新网站| 日本撒尿小便嘘嘘汇集6| 久久久久国产一级毛片高清牌| 亚洲欧洲日产国产| 成人影院久久| 久久久久久久久免费视频了| √禁漫天堂资源中文www| 国产97色在线日韩免费| 亚洲精品自拍成人| 中文欧美无线码| 亚洲天堂av无毛| 国产成人一区二区三区免费视频网站| 精品久久久久久电影网| 国产成人欧美在线观看 | 两人在一起打扑克的视频| 国产精品 欧美亚洲| 久久99热这里只频精品6学生| 欧美 亚洲 国产 日韩一| 欧美人与性动交α欧美软件| 变态另类成人亚洲欧美熟女 | 日韩视频在线欧美| 亚洲精品久久成人aⅴ小说| 91麻豆精品激情在线观看国产 | 一区二区三区国产精品乱码| 成人精品一区二区免费| 一边摸一边抽搐一进一出视频| 日本黄色视频三级网站网址 | 久久久久久久久免费视频了| 国产视频一区二区在线看| 国产aⅴ精品一区二区三区波| 免费不卡黄色视频| 91麻豆精品激情在线观看国产 | 一级片'在线观看视频| 一级毛片电影观看| 日日爽夜夜爽网站| 在线十欧美十亚洲十日本专区| 午夜激情av网站| 国产成人欧美| 国产日韩欧美在线精品| 久久久欧美国产精品| 国产欧美日韩精品亚洲av| 日韩视频一区二区在线观看| 老鸭窝网址在线观看| 欧美乱妇无乱码| 亚洲伊人久久精品综合| www.精华液| 国产激情久久老熟女| 我的亚洲天堂| 一本大道久久a久久精品| 黑人巨大精品欧美一区二区mp4| 老司机午夜十八禁免费视频| 亚洲精品自拍成人| 国产av又大| 18禁黄网站禁片午夜丰满| 成年女人毛片免费观看观看9 | 成年人黄色毛片网站| 亚洲伊人久久精品综合| 欧美国产精品一级二级三级| 色视频在线一区二区三区| 国产在视频线精品| 在线观看66精品国产| 黑人巨大精品欧美一区二区蜜桃| 99re在线观看精品视频| 五月开心婷婷网| 久久天躁狠狠躁夜夜2o2o| 久久香蕉激情| 免费久久久久久久精品成人欧美视频| 欧美久久黑人一区二区| 精品国产国语对白av| 亚洲欧美日韩另类电影网站| 韩国精品一区二区三区| 成年动漫av网址| 久久久水蜜桃国产精品网| 黄频高清免费视频| 一边摸一边抽搐一进一小说 | 亚洲国产av新网站| 色尼玛亚洲综合影院| 怎么达到女性高潮| 久久精品亚洲熟妇少妇任你| 中文亚洲av片在线观看爽 | 国产一区有黄有色的免费视频| 国产免费av片在线观看野外av| av有码第一页| av超薄肉色丝袜交足视频| 精品国产乱码久久久久久小说| 午夜成年电影在线免费观看| 丁香六月欧美| 国产精品98久久久久久宅男小说| 日韩制服丝袜自拍偷拍| 午夜福利一区二区在线看| 日韩三级视频一区二区三区| 交换朋友夫妻互换小说| 国产一区二区在线观看av| 高清在线国产一区| 天天躁日日躁夜夜躁夜夜| 国产精品亚洲一级av第二区| 在线观看人妻少妇| 十八禁网站免费在线| 亚洲精华国产精华精| 久久久久国产一级毛片高清牌| 50天的宝宝边吃奶边哭怎么回事| 无限看片的www在线观看| 久久久久久人人人人人| 免费在线观看完整版高清| 久久中文字幕人妻熟女| 波多野结衣一区麻豆| 老司机福利观看| 99精品在免费线老司机午夜| 国产精品99久久99久久久不卡| videos熟女内射| 一本一本久久a久久精品综合妖精| 俄罗斯特黄特色一大片| 欧美日韩视频精品一区| 精品久久蜜臀av无| 99国产精品99久久久久| 久久久精品区二区三区| 精品国产一区二区三区四区第35| 五月开心婷婷网| 麻豆av在线久日| 精品久久蜜臀av无| 亚洲av成人一区二区三| 国产亚洲一区二区精品| 国产又色又爽无遮挡免费看| 久久人妻福利社区极品人妻图片| 国产国语露脸激情在线看| 最新美女视频免费是黄的| 欧美乱码精品一区二区三区| 亚洲国产欧美日韩在线播放| 午夜福利一区二区在线看| 日韩欧美免费精品| 99香蕉大伊视频| 成人影院久久| 欧美性长视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 99国产精品免费福利视频| 人成视频在线观看免费观看| 国产精品成人在线| 国产97色在线日韩免费| 国精品久久久久久国模美| 亚洲色图av天堂| 一进一出好大好爽视频| 日韩三级视频一区二区三区| 久久精品91无色码中文字幕| 欧美一级毛片孕妇| 搡老熟女国产l中国老女人| 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区三| 满18在线观看网站| 可以免费在线观看a视频的电影网站| 搡老熟女国产l中国老女人| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品美女久久av网站| 亚洲人成电影观看| 国产片内射在线| 在线观看66精品国产| 精品人妻在线不人妻| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成a人片在线一区二区| 搡老熟女国产l中国老女人| 女人高潮潮喷娇喘18禁视频| 国产精品自产拍在线观看55亚洲 | 午夜老司机福利片| 国产视频一区二区在线看| 免费高清在线观看日韩| 黄色a级毛片大全视频| 欧美精品高潮呻吟av久久| 欧美日韩成人在线一区二区| 色老头精品视频在线观看| 日韩三级视频一区二区三区| 欧美另类亚洲清纯唯美| 国产亚洲精品第一综合不卡| 黄片大片在线免费观看| 亚洲午夜精品一区,二区,三区| 岛国毛片在线播放| 久久热在线av| 精品人妻在线不人妻| 人人妻人人澡人人看| 国产精品99久久99久久久不卡| videosex国产| 精品第一国产精品| 精品视频人人做人人爽| 国产精品美女特级片免费视频播放器 | 国产淫语在线视频| 日韩一卡2卡3卡4卡2021年| 少妇 在线观看| 黄色视频,在线免费观看| 午夜久久久在线观看| 亚洲专区中文字幕在线| 欧美精品一区二区免费开放| 变态另类成人亚洲欧美熟女 | 色老头精品视频在线观看| 国产在线免费精品| 淫妇啪啪啪对白视频| 巨乳人妻的诱惑在线观看| 中文字幕人妻熟女乱码| 激情在线观看视频在线高清 | 国产国语露脸激情在线看| 水蜜桃什么品种好| 男女之事视频高清在线观看| 亚洲国产欧美网| 欧美黄色片欧美黄色片| 成人18禁在线播放| 国产xxxxx性猛交| 99香蕉大伊视频| 九色亚洲精品在线播放| 久久99一区二区三区| 国产免费视频播放在线视频| av欧美777| 日本vs欧美在线观看视频| 国产精品久久久久久精品电影小说| 亚洲中文av在线| 黑人猛操日本美女一级片| 十八禁高潮呻吟视频| 香蕉国产在线看| 成年版毛片免费区| 精品国产一区二区三区久久久樱花| 国产一区二区 视频在线| 日本欧美视频一区| e午夜精品久久久久久久| 热re99久久精品国产66热6| 精品免费久久久久久久清纯 | 亚洲国产欧美日韩在线播放| 亚洲精华国产精华精| 美国免费a级毛片| 一二三四社区在线视频社区8| 女人精品久久久久毛片| 亚洲av成人不卡在线观看播放网| 欧美日韩视频精品一区| 老司机亚洲免费影院| 国产一区二区激情短视频| 中文字幕最新亚洲高清| av视频免费观看在线观看| 在线av久久热| 国产日韩欧美亚洲二区| 老司机午夜十八禁免费视频| 国产一区二区三区视频了| 亚洲成人手机| 无遮挡黄片免费观看| 麻豆国产av国片精品| 九色亚洲精品在线播放| 国产精品久久久久久人妻精品电影 | 国产精品国产av在线观看| 少妇的丰满在线观看| 好男人电影高清在线观看| 免费人妻精品一区二区三区视频| 亚洲欧美精品综合一区二区三区| 亚洲av日韩在线播放| xxxhd国产人妻xxx| 国产日韩欧美亚洲二区| 国产男靠女视频免费网站| 国产精品久久久久久精品古装| 99re6热这里在线精品视频| 国产极品粉嫩免费观看在线| 久久精品国产综合久久久| 宅男免费午夜| 国产精品 欧美亚洲| 在线观看人妻少妇| e午夜精品久久久久久久| 女人被躁到高潮嗷嗷叫费观| 波多野结衣一区麻豆| 黑人巨大精品欧美一区二区mp4| 狠狠狠狠99中文字幕| 日日爽夜夜爽网站| 搡老岳熟女国产| 国产在线视频一区二区| 欧美午夜高清在线| 亚洲天堂av无毛| 亚洲va日本ⅴa欧美va伊人久久| 人妻 亚洲 视频| 99久久精品国产亚洲精品| 午夜两性在线视频| 久热这里只有精品99| 如日韩欧美国产精品一区二区三区| 搡老熟女国产l中国老女人| 国产av一区二区精品久久| 久久香蕉激情| 超碰成人久久| 欧美黑人欧美精品刺激| 国产成人免费观看mmmm| 一级黄色大片毛片| av有码第一页| 精品久久蜜臀av无| 美国免费a级毛片| 男女无遮挡免费网站观看| 久久久久国内视频| 国产免费视频播放在线视频| 精品亚洲成a人片在线观看| av超薄肉色丝袜交足视频| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 男女边摸边吃奶| 亚洲精品av麻豆狂野| 国产在线观看jvid| 男人操女人黄网站| 无限看片的www在线观看| a级片在线免费高清观看视频| 人人妻人人澡人人看| 老熟妇乱子伦视频在线观看| 91精品三级在线观看| e午夜精品久久久久久久| 成年版毛片免费区| 一二三四社区在线视频社区8| 亚洲伊人久久精品综合| 久久久精品区二区三区| 99国产综合亚洲精品| 国产片内射在线| 老熟妇仑乱视频hdxx| 久久国产精品大桥未久av| 一级毛片电影观看| 999久久久国产精品视频| 大陆偷拍与自拍| 国产欧美日韩综合在线一区二区| 一边摸一边抽搐一进一出视频| 国产欧美日韩精品亚洲av| 欧美日韩福利视频一区二区| 欧美性长视频在线观看| 18禁国产床啪视频网站| 女人精品久久久久毛片| 动漫黄色视频在线观看| 亚洲av第一区精品v没综合| 人人妻人人添人人爽欧美一区卜| 丁香六月天网| 久久久久久久精品吃奶| 久久青草综合色| 丁香六月天网| 在线观看免费日韩欧美大片| 天天躁日日躁夜夜躁夜夜| 另类亚洲欧美激情| 亚洲人成伊人成综合网2020| 国产精品久久久久久精品电影小说| 国产一区二区三区综合在线观看| av又黄又爽大尺度在线免费看| 亚洲精品中文字幕在线视频| 亚洲精品中文字幕一二三四区 | netflix在线观看网站| 久久久精品免费免费高清| 极品人妻少妇av视频| 国产精品1区2区在线观看. | 69av精品久久久久久 | 亚洲熟女精品中文字幕| 每晚都被弄得嗷嗷叫到高潮| 国产成人一区二区三区免费视频网站| 新久久久久国产一级毛片| 亚洲免费av在线视频| 久久久国产欧美日韩av| 精品久久蜜臀av无| 日韩欧美一区视频在线观看| 亚洲成国产人片在线观看| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 蜜桃在线观看..| 亚洲国产av新网站| 午夜精品国产一区二区电影| 高清毛片免费观看视频网站 | kizo精华| 亚洲综合色网址| 男人舔女人的私密视频| 男女边摸边吃奶| 丰满少妇做爰视频| 午夜福利一区二区在线看| 日本wwww免费看| 成人特级黄色片久久久久久久 | 免费不卡黄色视频| 亚洲精品一卡2卡三卡4卡5卡| 免费高清在线观看日韩| 人人澡人人妻人| h视频一区二区三区| 久久这里只有精品19| 麻豆乱淫一区二区| 少妇裸体淫交视频免费看高清 | 天堂俺去俺来也www色官网| 少妇裸体淫交视频免费看高清 | 国产欧美日韩一区二区精品| 日韩免费av在线播放| 色婷婷av一区二区三区视频| 丁香欧美五月| 亚洲精品国产色婷婷电影| 人妻一区二区av| 亚洲色图综合在线观看| 国产av精品麻豆| 欧美国产精品va在线观看不卡| av网站免费在线观看视频| 99国产精品一区二区三区| 亚洲av美国av| 亚洲,欧美精品.| 性少妇av在线| 97人妻天天添夜夜摸| 成人免费观看视频高清| 色综合欧美亚洲国产小说| 老熟女久久久| 桃红色精品国产亚洲av| 2018国产大陆天天弄谢| 亚洲 国产 在线| 高清毛片免费观看视频网站 | 免费看a级黄色片| 制服诱惑二区| av免费在线观看网站| 一边摸一边做爽爽视频免费| 国产免费av片在线观看野外av| 久久热在线av| 国产免费现黄频在线看| 免费日韩欧美在线观看| 超色免费av| 久久久久久人人人人人| 亚洲成国产人片在线观看| 午夜福利乱码中文字幕| 久久人人爽av亚洲精品天堂| 国产亚洲欧美在线一区二区| 我要看黄色一级片免费的| 国产精品 国内视频| 亚洲精品国产一区二区精华液| 午夜福利一区二区在线看| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品在线电影| 美女国产高潮福利片在线看| 在线观看免费日韩欧美大片| 老熟女久久久| 亚洲国产中文字幕在线视频| 亚洲久久久国产精品| 亚洲国产欧美网| 久久99一区二区三区| 日韩一区二区三区影片| 欧美人与性动交α欧美精品济南到| 乱人伦中国视频| 18禁裸乳无遮挡动漫免费视频| h视频一区二区三区| 黄片播放在线免费| 国产成人精品在线电影| 下体分泌物呈黄色| 麻豆国产av国片精品| 极品少妇高潮喷水抽搐| 男女高潮啪啪啪动态图| av有码第一页| 久久久久网色| 成人黄色视频免费在线看| 极品少妇高潮喷水抽搐| 高清在线国产一区| 性高湖久久久久久久久免费观看| 女性被躁到高潮视频| 久久精品亚洲av国产电影网| 亚洲少妇的诱惑av| 免费av中文字幕在线| 国产精品二区激情视频| 一区二区三区乱码不卡18| 啪啪无遮挡十八禁网站| 女人被躁到高潮嗷嗷叫费观| 大码成人一级视频| 男女床上黄色一级片免费看| 国产亚洲精品一区二区www | 大片电影免费在线观看免费| 十八禁网站网址无遮挡| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久免费高清国产稀缺| 久久久精品国产亚洲av高清涩受| av欧美777| 免费在线观看日本一区| 黄频高清免费视频| 精品一区二区三区视频在线观看免费 | 国产在线一区二区三区精| 亚洲精品国产一区二区精华液| 99国产精品免费福利视频| 国产亚洲精品一区二区www | 国产aⅴ精品一区二区三区波| 成人免费观看视频高清| 黄色视频在线播放观看不卡| 国产精品美女特级片免费视频播放器 | 91国产中文字幕| 亚洲精品在线观看二区| 男女高潮啪啪啪动态图| 久久精品人人爽人人爽视色| 午夜激情av网站| 久久香蕉激情| 国产又色又爽无遮挡免费看| av国产精品久久久久影院| 狠狠婷婷综合久久久久久88av| 在线观看免费视频日本深夜| 日韩一卡2卡3卡4卡2021年| 考比视频在线观看| 国产1区2区3区精品| 亚洲欧美一区二区三区黑人| 999精品在线视频| 狂野欧美激情性xxxx| 国产精品免费一区二区三区在线 | 亚洲一区中文字幕在线| 美女主播在线视频| av在线播放免费不卡| 日本av手机在线免费观看| 狠狠婷婷综合久久久久久88av| 建设人人有责人人尽责人人享有的| 99久久精品国产亚洲精品| 他把我摸到了高潮在线观看 | 国产成人精品无人区| 欧美亚洲日本最大视频资源| 亚洲第一欧美日韩一区二区三区 | 国产精品欧美亚洲77777| 欧美激情极品国产一区二区三区| 国产成+人综合+亚洲专区| 嫁个100分男人电影在线观看| 99国产极品粉嫩在线观看| www.熟女人妻精品国产| 国产成人精品在线电影| 国产亚洲av高清不卡| 操美女的视频在线观看| 国产伦理片在线播放av一区| 高清欧美精品videossex| 久久性视频一级片| 波多野结衣av一区二区av| 国产在线精品亚洲第一网站| 久久人妻福利社区极品人妻图片| 精品国产超薄肉色丝袜足j| 可以免费在线观看a视频的电影网站| 99精品欧美一区二区三区四区| 国产精品久久久久成人av| 91国产中文字幕| 国产在线观看jvid| 日日夜夜操网爽| 久久久久久久久久久久大奶| 丝瓜视频免费看黄片| 伦理电影免费视频| 国内毛片毛片毛片毛片毛片| 亚洲精品国产区一区二| 色综合欧美亚洲国产小说| 国产一区二区激情短视频| 国产精品一区二区精品视频观看| 菩萨蛮人人尽说江南好唐韦庄| 宅男免费午夜| 99国产精品一区二区三区| 最黄视频免费看| 日韩欧美一区视频在线观看| 亚洲午夜理论影院| 黄片大片在线免费观看| 日本欧美视频一区| 真人做人爱边吃奶动态| 最近最新免费中文字幕在线| 99久久人妻综合| 最近最新中文字幕大全电影3 | 麻豆av在线久日| videos熟女内射| 日本av手机在线免费观看| 国产精品久久久久久人妻精品电影 | 亚洲一区二区三区欧美精品| 男女午夜视频在线观看| 黄色a级毛片大全视频| 交换朋友夫妻互换小说| 精品亚洲乱码少妇综合久久| 亚洲avbb在线观看| 他把我摸到了高潮在线观看 | 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 成人免费观看视频高清| 中文字幕高清在线视频| 免费观看人在逋| 亚洲专区国产一区二区| 成人黄色视频免费在线看| 国产精品影院久久| 欧美成人免费av一区二区三区 | 国产欧美日韩一区二区三区在线| 黑人猛操日本美女一级片| 久久中文字幕一级| 成人av一区二区三区在线看| 多毛熟女@视频| 久久中文字幕一级| 成人av一区二区三区在线看| 在线播放国产精品三级| 国产伦人伦偷精品视频| 亚洲专区中文字幕在线| 久久人妻av系列| 亚洲一码二码三码区别大吗| av线在线观看网站| 人人妻人人添人人爽欧美一区卜|