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

    機載相控陣PD雷達的MPRF設計與選擇

    2016-10-24 03:11:30連曉鋒湯子躍朱振波汪先超席秋實喬寧
    現(xiàn)代防御技術 2016年4期

    連曉鋒, 湯子躍, 朱振波,汪先超, 席秋實, 喬寧

    (1.空軍預警學院 三系,湖北 武漢 430019; 2.中國人民解放軍95028部隊,湖北 武漢 430079)

    ?

    機載相控陣PD雷達的MPRF設計與選擇

    連曉鋒1, 湯子躍1, 朱振波1,汪先超1, 席秋實2, 喬寧1

    (1.空軍預警學院 三系,湖北 武漢430019; 2.中國人民解放軍95028部隊,湖北 武漢430079)

    針對機載相控陣PD雷達波束寬度隨掃描角變化的特點,詳細闡述了中脈沖重復頻率(MPRF)設計與選擇的具體方法步驟。從M/N準則選擇、重頻范圍2方面確定深入研究了MPRF設計問題;從優(yōu)化模型、約束條件、優(yōu)化算法3方面深入探討了MPRF選擇問題,提出利用可調參數(shù)少、收斂速度快、全局尋優(yōu)能力強的粒子群算法選擇最優(yōu)MPRF。最后通過具體實例說明了機載相控陣PD雷達MPRF設計與選擇方法的有效性。

    機載預警雷達;相控陣;脈沖多普勒;MPRF設計;MPRF選擇;粒子群優(yōu)化算法

    0 引言

    機載相控陣脈沖多普勒(pulse-Doppler,PD)雷達多采用中脈沖重復頻率(medium pulse repetition frequency,MPRF)模式。機載相控陣PD雷達在原理上與機掃體制的機載PD雷達相同,但相掃體制的機載PD雷達不同掃描角分配的時間是不一樣的,而機掃體制的機載PD雷達所有方位分配的時間是均勻的,這就導致相掃體制與機掃體制的機載PD雷達的MPRF設計和選擇是不同的。目前國內外很多學者對機載PD雷達的MPRF設計與選擇進行了相關研究,大多針對機掃體制的機載PD雷達提出了很多PRF選擇方法。文獻[1]通過典型約束全面搜索最優(yōu)的MPRF組合,復雜、費時;文獻[2]提出了噪聲環(huán)境下的自適應PRF優(yōu)化選擇算法。文獻[3-5]運用遺傳算法對機載PD雷達的MPRF進行優(yōu)化選擇。文獻[6]提出運用模擬退火算法選擇最優(yōu)MPRF。但對采用相掃體制的機載PD雷達的MPRF設計和選擇的研究甚少,目前僅有文獻[7]對相控陣機載PD雷達PRF設計與選擇進行了相關研究,但該文獻中脈沖補償和重頻間隔的計算存在偏差,與實際情況不符,需要糾正。

    針對機載相控陣PD雷達波束寬度隨掃描角變化的特點,本文詳細闡述了MPRF設計與選擇的具體方法步驟。從M/N準則選擇、重頻(脈沖重復頻率)范圍確定深入研究了MPRF設計問題;從優(yōu)化模型、約束條件、優(yōu)化算法3方面深入探討了MPRF選擇問題,同時提出利用可調參數(shù)少、收斂速度快、全局尋優(yōu)能力強的粒子群算法選擇最優(yōu)MPRF。

    1 機載相控陣PD雷達的MPRF設計

    1.1時間資源分配

    對于機載相控陣PD雷達,波束寬度隨天線掃描角的增大而展寬,導致天線性能下降。為獲得全方位比較一致的探測性能,需要增加駐留時間來補償。工程上通常將相掃的-60°~60°區(qū)域分為7個扇區(qū),同一扇區(qū)內各波位的掃描間隔相同,波位駐留時間相同??紤]到不同波長在不同掃描角時波瓣展寬程度不同,所以采用的頻段不同,分區(qū)應有所不同[8]。對于頻率在3 GHz以下的機載 相控陣PD雷達,分區(qū)1為-60°~-45°,-45°~-30°,-30°~-15°,-15°~15°,15°~30°,30°~45°和45°~60°,對于頻率在3 GHz以上的機載相控陣PD雷達,分區(qū)2為-60°~-55°,-55°~-45°,-45°~-30°,-30°~30°,30°~45°,45°~55°和55°~60°。

    掃描區(qū)域內某一個波位上的駐留時間稱為景時間,每個波位上的景時間包括各個PRF對應的回波脈沖持續(xù)時間(稱為幀時間)[9],波位內各個重頻分配的時間是相等的。每個PRF的幀時間包含積累時間和填充時間,積累時間內的相干脈沖是用來進行目標檢測的。填充時間用以建立穩(wěn)定的AGC控制,會帶來一定的檢測損失,對于一般的相控陣雷達而言,各個波位上的填充時間一般是固定的,等于2倍的最大距離除以光速。填充時間內的填充脈沖主要是避免重頻切換時前一重頻的遠距離回波信號混入后一重頻的接收信號之中形成干擾,最遠處的回波返回之前的脈沖是不能拿來做處理的。

    為了方便闡述MPRF設計的過程,這里對機載相控陣PD雷達的相關參數(shù)作一個假設。假設某機載相控陣PD雷達采用“平衡木”天線外形,正側視放置,天線由長方陣組成,負責120°掃描空間,天線陣列為64×16,陣元間距半波長,波長0.15 m,載機速度>120 m/s,正常掃描模式下ts=3.333 3 s,最大作用距離400 km。根據(jù)相應的雷達參數(shù),建立分區(qū)1后總時間分配為

    Ntprf(pac1+pac2+…+pacK)+Ntfill,

    (1)

    (2)

    (3)

    從式(2),(3)可以看出,重頻f和重頻數(shù)N的確定即MPRF組的確定與機載預警雷達的時間資源的分配密切相關,這就決定了重頻f和重頻數(shù)N之間相互制約、相互影響。

    另外,根據(jù)相控陣雷達波束寬度計算公式可以確定各扇區(qū)掃描步長,進而確定各扇區(qū)的波位數(shù),同時提出了脈沖補償因子的概念用以分析相干脈沖。

    (1) 波位數(shù)

    各扇區(qū)內的波位數(shù)是根據(jù)掃描間隔確定的,可以利用相控陣雷達波束控制靈活的特點,根據(jù)需要控制掃描間隔,在大掃描角上適當?shù)卦龃髵呙栝g隔。相控陣方位半功率點波束寬度計算公式為[10]

    (4)

    式中:θd為波束掃描角;d為陣元間距。

    相掃的-60°~60°區(qū)域內方位向半功率點波束寬度的范圍在1.586 5°~3.172 8°之間,平均為2.00 7°,則平均波位數(shù)

    (5)

    (2) 脈沖補償

    相控陣雷達需要增加積累時間來補償因波束展寬引起的信噪比損失。各個方位的發(fā)射天線增益相比法線方向降低多少,則積累時間補償多少,其根據(jù)是發(fā)射天線增益和積累時間在雷達方程中對作用距離的影響是等效的[9]。增加積累時間通過增加相干脈沖數(shù)的方法(即脈沖補償)實現(xiàn)的,各扇區(qū)波位內脈沖補償程度如何確定,本文定義了脈沖補償因子。脈沖補償因子:對某一重頻,扇區(qū)內掃描步長與0°對應扇區(qū)內掃描步長的比值。根據(jù)計算得出的各扇區(qū)內掃描間隔與0°對應扇區(qū)內掃描間隔的比值,確定各扇區(qū)對應的脈沖數(shù)補償因子分別為2.0,1.4,1.2,1,1.2,1.4,2.0。則脈沖補償因子和某一重頻總的相干脈沖數(shù)P0之間關系為

    P0=α·p?K/2」,

    (6)

    式中:α為補償因子的組合;p?K/2」為0°掃描角時的相干脈沖數(shù),?·」表示向下取整,故

    P0=α·p?K/2」=(19+1.2×9×2+1.4×7×2+2×5×2)p?K/2」=80.2p?K/2」.

    1.2M/N準則的確立

    準則的選擇是重頻選擇算法的前提,M/N準則的確定涉及2方面:一是N的選擇,即總的重頻個數(shù)的確定,這主要從相控陣雷達時間分配角度考慮,確保滿足指標要求;二是M的選擇,即檢測所需重頻個數(shù),主要從虛警概率問題考慮。

    (1)N的確定

    由式(2),(3)可知,重頻f和N的確定與天線掃描周期、波位數(shù)、各個波位的相干脈沖數(shù)和填充時間均有關,天線掃描周期已經(jīng)給出,波位數(shù)、填充時間均可計算,則相干脈沖是N確定的關鍵因素。由分析可知,0°掃描角對應波位分配的時間最少,只要0°掃描角分配的駐留時間能夠滿足要求則其他掃描角均可以滿足,因各波位的填充時間相同且可計算,所以下面只考察0°掃描角的積累時間。

    B?K/2」=p?K/2」tprf.

    (7)

    主雜波寬度與相控陣列的波束寬度密切相關,但實際過程中決定主雜波譜寬還有一個重要參數(shù):相干脈沖數(shù)。圖1給出了相干脈沖與多普勒譜寬的曲線。

    圖1 脈沖數(shù)與歸一化譜寬的關系Fig.1 Relationship between pulse number and normalized spectrum width

    從圖1中可以看出,若要求主雜波譜寬小于20%,則所需相干脈沖數(shù)最少為28左右,所以p?K/2」≥cL=28,其中cL為0°掃描角所需的最小脈沖數(shù)。主雜波譜寬的確定就對最小脈沖數(shù)設置了一個下限,進而對總的重頻數(shù)設置了一個下限。M/N準則中的N的選擇應該滿足

    (8)

    當重頻范圍確定時就可以確定重頻數(shù)的范圍,具體重頻數(shù)的確定提供依據(jù)。另外,由于時間的限制,0°掃描角的脈沖數(shù)是有上限的,如果脈沖數(shù)取的太大會嚴重影響數(shù)據(jù)率,所以某一重頻0°掃描角的相干脈沖數(shù)應該滿足

    cL≤p?K/2」≤cH,cL

    (9)

    式中:cL為下限;cH為上限,一般考慮上限為平均脈沖數(shù)

    (10)

    正因為重頻f、重頻數(shù)N和相干脈沖數(shù)是相互制約、相互影響,由式(1)可確定另一相干脈沖數(shù)上限

    (11)

    (2)M的確定

    對MPRF而言,機載預警雷達一般采用二維CFAR進行恒虛警處理,可以提高目標的信干比,有利于目標檢測。這里對于M的選擇,只考慮一下M/N準則對虛警概率的影響。最終的虛警概率的計算公式如下:

    (12)

    式中:PF為最終的虛警概率;Pf是恒虛警處理一級所要求的虛警概率。圖2中為一些常用的M/N檢測準則與虛警率的關系。

    圖2 Pf與PF虛警率間的關系Fig.2 Relationship between Pf and PF

    從圖2中可以看出,在Pf一定的情況下,M相同,N減少會降低虛警率;N相同,M增大也會降低虛警率。還需要注意的是M的選擇一般在N/2左右的正整數(shù),如典型的3/8準則、3/5準則等。理論上利用兩重PRF就能對單個目標距離和多普勒頻率進行解模糊。但是,為了降低由于存在多目標所引起的虛警概率,至少需要用三重PRF進行檢測。

    1.3重頻范圍的選擇

    (1) 最小重頻的選擇

    最小重頻的選擇應從以下幾個方面入手:

    1) 機載PD雷達中使用MPRF的目的就是為了保證在副瓣雜波中具有良好的檢測性能,所以,選擇MPRF組時應盡量減小檢測單元內副瓣雜波的強度[11]。

    (13)

    式中:fmc為主瓣雜波抑制寬度(遮擋區(qū));τ為發(fā)射脈沖的寬度。

    2) 事實上,式(13)給出的是理想的最小重頻,但實際上,當M/N準則確定之后從式(2)也可以得出一個最小重頻

    (14)

    3) 另外,對中重頻而言,雷達參數(shù)確定之后一般要保證載機運動的多普勒頻率與重復頻率之間滿足[7]

    (15)

    式中:γ為正側面陣時空時二維雜波分布斜率;fd為載機巡航速度時的多普勒頻率,γ≥0.9。由式(15)可得

    fmin3≥2fdγ.

    (16)

    所以這里選取的最小重頻是

    (17)

    (2) 最大重頻的選擇

    最大重頻是重頻范圍的上限,一般考慮占空比的限制,兩者關系如下:

    (18)

    式中:du為占空比;τ為發(fā)射脈沖的寬度。

    事實上式(2)還給出了f1的一個限制,即一般情況下由于時間有限,式中脈沖數(shù)是有上限的,即由式(10)和(11)共同確定的脈沖數(shù)的上限。這樣,搜索重頻范圍中的最大值應為式(18),而f1重頻的范圍應該滿足

    (19)

    而第N個重頻則應該滿足

    (20)

    2 機載預警雷達MPRF選擇

    在機載預警雷達MPRF設計的基礎上,運用粒子群優(yōu)化算法選擇最優(yōu)MPRF組以滿足工程需要。重頻組選擇的目的旨在解距離模糊和速度模糊的同時,使得距離-速度二維盲區(qū)最小,即在M/N準則前提下在最小、最大重頻范圍內選擇最優(yōu)MPRF組使檢測的盲區(qū)覆蓋率最小。重頻組的選擇會受到雷達工作性能對其施加的諸多限制和約束,這實際上是一種解大規(guī)模組合優(yōu)化問題[12]。

    2.1優(yōu)化模型

    2.2約束條件

    機載預警雷達中MPRF的選擇還需定義約束條件,改變其中任意一個約束條件都會導致重頻選擇過程重新進行。本文從強制性約束和選擇性約束分別進行了討論。

    (1) 強制性約束

    1) 盲區(qū)約束

    a) 距離除不模糊的低速區(qū)外不存在盲區(qū)。

    b) 速度除不模糊的低速區(qū)外不存在盲區(qū)。

    實際上在距離-速度二維盲區(qū)最小的情況下,并不能保證距離除不模糊的低速區(qū)外、速度除不模糊的低速區(qū)外不存在盲區(qū),不可避免會出現(xiàn)一些盲區(qū)[12]。

    2) 重頻選擇范圍

    根據(jù)式(19),(20)可以確定最小頻率和最大頻率的范圍。

    3) 解模糊

    (21)

    (22)

    式中: Rmax最大作用距離;Dmax為最大多普勒頻率。式(21)和(22)保證解模糊的正確性。

    (2) 選擇性約束

    1) 重頻間隔

    重頻的間隔是為了保證在M/N準則下所有距離的特定速度上的盲區(qū)不存在。重頻間隔應滿足:

    (23)

    2) 數(shù)值約束

    數(shù)值約束是指對重復頻率值的約束,本文要求重復頻率值互為素數(shù)。

    2.3粒子群算法

    粒子群優(yōu)化(partical swarm optimization,PSO)算法是近年來發(fā)展起來的一種新的基于疊代的進化算法[14],種群中每個成員叫做粒子,代表著一個潛在的可行解[15],而食物的位置代表全局最優(yōu)解。每個粒子都有自己的位置和速度(決定飛行的方向和距離),以及一個由被優(yōu)化函數(shù)決定的適應度值。每次迭代中,粒子通過跟蹤2個最優(yōu)解動態(tài)調整自己的速度和位置,最終接近食物的位置,達到從全空間搜索最優(yōu)解的目的:一個是粒子自身找到的最優(yōu)解,叫做個體最優(yōu)解;另一個是種群找到的最優(yōu)解,叫做全局最優(yōu)解。粒子可以看作是D維空間中的一點,把第i個粒子的位置表示為xi(t)=(xi1,xi2,…,xiD),而把速度表示為vi(t)=(vi1,vi2,…,viD),個體最優(yōu)解表示為pi(t)=(pi1,pi2,…,piD),所有粒子中全局最優(yōu)解的下標用g表示。粒子速度和位置的更新方程如下[13]

    vid(t+1)=vid(t)+c1r1(pid(t)-xid(t))+

    c2r2(pgd(t)-xid(t)),

    (24)

    (25)

    xid(t+1)=xid(t)+vid(t+1),

    (26)

    圖3 PSO流程圖Fig.3 Flow chart of PSO

    3 仿真分析

    本文運用Matlab編程工具,對基于標準粒子群優(yōu)化算法的最優(yōu)PRF組的選擇方法進行仿真實驗,并與遺傳算法的求解結果進行對比。

    粒子群算法參數(shù):迭代次數(shù)M=200,種群規(guī)模N=30,Vmax=fmax-fmin,Vmin=fmin-fmax,c1=c2=1.494 45,以二維盲區(qū)遮蔽率作為個體適應值。

    則各種約束值為

    N最大值范圍3~8之間,從時間資源分配角度考慮N取5比較合適。選取3/5準則,不僅能滿足解模糊的需要,還能很好減小“重影”影響,則

    4 455.6 Hz,

    所以cH的范圍在36~48之間,這里選42(1/2頻率間隔)。即各重頻的范圍

    6 683.3 Hz≤f5≤fmax=7 700 Hz,

    fi-1≤fi≤fi+1,i=2,3,4,

    則149≤tprf1≤224,130≤tprf5≤150。

    運用標準粒子群優(yōu)化算法搜索得到的脈沖重復頻率為prf=[5 957,6 298,6 646,6 982,7 692] (Hz)。0°掃描角的相干脈沖數(shù)是38,40,42,44,48。運用遺傳算法搜索得到的脈沖重復頻率為prf=[5 210,6 329,6 771,7 127,7 439]。圖4給出了2種算法的最優(yōu)MPRF的迭代過程,圖5給出了2種算法的最優(yōu)MPRF對應的盲區(qū)圖。

    圖4 選擇最優(yōu)MPRF的迭代過程Fig.4 Iterative process of selecting the optimal MPRF

    圖5 最優(yōu)MPRF對應的盲區(qū)圖Fig.5 Blind spot corresponding to the optimal MPRF

    從圖4中可以看出,標準粒子群優(yōu)化算法相比遺傳算法收斂要快。從圖5中可以看出,標準粒子群優(yōu)化算法所得到的遮蔽率優(yōu)于基于遺傳算法得到的遮蔽率,所提算法能取得較好的效果。

    4 結束語

    針對機載相控陣PD雷達波束寬度隨掃描角變化的特點,本文詳細闡述了MPRF設計與選擇的具體方法步驟。實驗表明,本文所提方法能得到所需MPRF組,滿足機載相控陣PD雷達的工作需要,為機載相控陣PD雷達的MPRF波形設計提供了依據(jù)。

    [1]SIMPSON J. PRF Set Selection for Pulse Doppler Radars [C]∥The IEEE Region 5 Conference, 1988: 38-44.

    [2]MOORMAN R A, WESTERKAMP J J. Maximizing Noise-Limited Detection Performance in Medium PRF Radars by Optimizing PRF Visivility[C]∥Proc. IEEE National Aerospace and Electronics Conference, July 1993,Vol. l:288-293.

    [3]DAVIES P G, HUGHES E J. Medium PRF Set Selection Using Evolutionary Algorithms[J]. IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(3): 933-939.

    [4]ALABASTER C M, HUGHES E J, MATTHEW J H.Medium PRF Radar PRF Selection Using Evolutionary Algorithms[J]. IEEE Transactions on Aerospace and Electronic Systems, 2003,39(3): 990-1001.

    [5]ALABASTER C M, HUGHES E S. The Design of Medium PRF Radar Schedules for Optimum Delectability in Diverse Clutter Scenes[C]∥ Proc. IEEE Waveform Diversity and Design, 2006.

    [6]Soyeon Ahn, Heedeok Lee, Byungwook Jung. Medium PRF Set Selection for Pulsed Doppler Radars Using Simulated Annealing[C]∥IEEE National Radar Conference-Proceedings,2011,10(1):91-94.

    [7]侯耀清,陳輝,母其勇.相控陣機載PD雷達的波形設計問題[J]. 雷達科學與技術, 2007,5(1): 22-28.

    HOU Yao-qing, CHEN Hui, MU Qi-yong. Waveform Design for Phased-Array PD Radar[J]. Radar Science and Technology, 2007,5(1):22-28.

    [8]張曉波, 曹晨. 機載預警雷達最優(yōu)探測波長研究[J]. 中國電子科學研究院學報, 2011,6(3): 263-266.

    ZHANG Xiao-bo, CAO Chen. A Study on Airborne Early Warning Radar Optimal Detection Wavelength[J]. Journal of CAEIT, 2011, 6(3): 263-266.

    [9]曹晨. 機載相控陣PD雷達性能分析中的重要問題[J]. 現(xiàn)代雷達, 2008, 30(6):14-16.

    CAO Chen. Problems in Airborne Phased-Array PD Radar Performance Analysis[J]. Modern Radar, 2008, 30(6): 14-16.

    [10]張光義. 相控陣雷達原理[M].北京:國防工業(yè)出版社,2009:29-30.

    ZHANG Guang-yi.The Theory of Phase Array Radar[M].Beijing:National Defence Industry Press,2009:29-30.

    [11]陳國海. 機載脈沖多普勒雷達的中脈沖重復頻率的波形設計[J]. 現(xiàn)代雷達, 1999,21(1):14-18.

    CHEN Guo-hai. MPRF Wave-Form Design for Airborne PD Radar[J]. Modern Radar,1999,21(1):14-18.

    [12]韓偉. 機載預警雷達航跡綜合優(yōu)化技術研究[D]. 武漢: 空軍預警學院, 2009.

    HAN Wei. Study on the Synthesis Optimization of Airborne Early Warning Radar Track[D]. Wuhan: Airforce Early Warning Academy, 2009.

    [13]王鵬, 譚賢四, 王紅, 等. 基于遺傳算法的地面中重頻PD雷達PRF組選擇[J]. 現(xiàn)代雷達, 2008, 30(9):33-35.

    WANG Peng, TAN Xian-si, WANG Hong, et al. PRF Set Selection for Ground-Based MPRF PD Radars Based on Genetic Algorithm[J]. Modern Radar, 2008, 30(9): 33-35.

    [14]沈伋, 韓麗川, 沈益斌. 基于粒子群算法的飛機總體參數(shù)優(yōu)化[J]. 航空學報, 2008, 29(6):1538-1541.

    SHEN Ji, HAN Li-chuan, SHEN Yi-bin. Optimization of Airplane Primary Parameters Based on Particle Swarm Algorithm[J]. Acta Aeronautica et Astronautica Sinca, 2008, 29(6):1538-1541.

    [15]KENNEDY J, EBERHART R C. Particle Swarm Opti mization[C]∥Proc. IEEE International Conference on Neural Networks, 1995:1942-1948.

    MPRF Design and Selection of Airborne PD Radar

    LIAN Xiao-feng1, TANG Zi-yue1, ZHU Zhen-bo1, WANG Xian-chao1,XI Qiu-shi2, QIAO Ning1

    (1. Air Force Early Warning Academy,Department 3,Hubei Wuhan 430019, China;2. PLA,No.95028 Troop,Hubei Wuhan 430079, China)

    The detailed methods and steps of medium pulse repetition frequency(MPRF) design and selection are described based on the characteristics of phase array radar. The MPRF design problem is analyzed fromM/Nscheme selection and the pulse recurrence frequency(PRF) range, and The MPRF selection problem is discussed from three aspects of optimization model, constraints and optimization algorithm. Swarm optimization algorithm is developed to select optimum MPRF, for the algorithm has simple realization process, less adjustable parameters, higher convergence speed and stronger global optimization. Finally, an example is given to demonstrate the effectiveness of the proposed MPRF design and selection method.

    airborne early warning radar; phase array; pulse-Doppler(PD); medium pulse repetition frequency(MPRF) design; MPRF selection; particle swarm optimization

    2015-05-18;

    2015-10-15

    青年創(chuàng)新基金(2013QNCX0101);中國博士后科學基金(2015M572779)

    連曉鋒(1990-),男,河北石家莊人。碩士生,主要研究方向為目標檢測與識別。

    通信地址:430019武漢市江岸區(qū)黃浦大街288號空軍預警學院3系E-mail:1248656228@qq.com

    10.3969/j.issn.1009-086x.2016.04.021

    TN959.73;TN958.92

    A

    1009-086X(2016)-04-0129-07

    天堂影院成人在线观看| 国产精品香港三级国产av潘金莲| 成人三级黄色视频| 亚洲片人在线观看| 美女高潮喷水抽搐中文字幕| 亚洲中文av在线| 美女高潮喷水抽搐中文字幕| 亚洲中文av在线| 啦啦啦 在线观看视频| 男女床上黄色一级片免费看| 久久久久久大精品| 狠狠狠狠99中文字幕| 香蕉丝袜av| 亚洲性夜色夜夜综合| 欧美日韩亚洲综合一区二区三区_| 好男人在线观看高清免费视频 | 久热这里只有精品99| 叶爱在线成人免费视频播放| 欧美成狂野欧美在线观看| 满18在线观看网站| 中文字幕久久专区| 成熟少妇高潮喷水视频| 久久久国产成人精品二区| 非洲黑人性xxxx精品又粗又长| 中文字幕人妻熟女乱码| 久久九九热精品免费| 后天国语完整版免费观看| 一级作爱视频免费观看| 青草久久国产| 俺也久久电影网| 大香蕉久久成人网| 麻豆国产av国片精品| 十八禁网站免费在线| 久久精品国产99精品国产亚洲性色| 最近最新中文字幕大全免费视频| 国产黄a三级三级三级人| 国产精品九九99| av有码第一页| 亚洲av电影在线进入| or卡值多少钱| 少妇裸体淫交视频免费看高清 | 亚洲成人免费电影在线观看| 精品国产超薄肉色丝袜足j| 亚洲专区中文字幕在线| 国产欧美日韩精品亚洲av| 亚洲黑人精品在线| 亚洲精华国产精华精| 无人区码免费观看不卡| 在线国产一区二区在线| 精品电影一区二区在线| 国产一区二区激情短视频| 欧美乱妇无乱码| 亚洲av电影不卡..在线观看| 一级毛片高清免费大全| 嫩草影视91久久| 久久久久久久久久黄片| 久久午夜亚洲精品久久| 人人澡人人妻人| 一二三四在线观看免费中文在| 夜夜看夜夜爽夜夜摸| 18禁观看日本| 99在线人妻在线中文字幕| 国产精华一区二区三区| 脱女人内裤的视频| 亚洲国产欧美日韩在线播放| 亚洲精品色激情综合| 在线观看午夜福利视频| www.熟女人妻精品国产| 亚洲最大成人中文| 日本成人三级电影网站| 久久久久亚洲av毛片大全| 日本五十路高清| 搡老岳熟女国产| 一区福利在线观看| 亚洲精品国产精品久久久不卡| 国产亚洲欧美在线一区二区| 大型av网站在线播放| 99久久无色码亚洲精品果冻| 大香蕉久久成人网| 91麻豆av在线| 脱女人内裤的视频| 两性午夜刺激爽爽歪歪视频在线观看 | 美女高潮喷水抽搐中文字幕| 嫁个100分男人电影在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品99久久99久久久不卡| 国产精品香港三级国产av潘金莲| 亚洲aⅴ乱码一区二区在线播放 | 激情在线观看视频在线高清| 国产高清激情床上av| 久久亚洲真实| 亚洲一区高清亚洲精品| 色老头精品视频在线观看| 一a级毛片在线观看| 国产精品av久久久久免费| 中亚洲国语对白在线视频| 亚洲国产欧美日韩在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产日韩欧美精品在线观看 | 淫妇啪啪啪对白视频| 中文字幕高清在线视频| 亚洲av五月六月丁香网| 一进一出抽搐gif免费好疼| 午夜福利18| 他把我摸到了高潮在线观看| 国产精品九九99| 人成视频在线观看免费观看| 精品熟女少妇八av免费久了| 国产成人欧美在线观看| 91老司机精品| 两性午夜刺激爽爽歪歪视频在线观看 | 此物有八面人人有两片| 日本a在线网址| 欧美中文综合在线视频| 国产视频一区二区在线看| 国产精品久久久久久精品电影 | 男人舔女人的私密视频| 亚洲熟妇中文字幕五十中出| 国产精品,欧美在线| 国产精品亚洲美女久久久| 亚洲 欧美一区二区三区| 亚洲五月婷婷丁香| 大香蕉久久成人网| 首页视频小说图片口味搜索| 久久午夜综合久久蜜桃| 99久久99久久久精品蜜桃| 又紧又爽又黄一区二区| 久久精品人妻少妇| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美在线一区二区| 亚洲国产精品成人综合色| 精品熟女少妇八av免费久了| 黄色a级毛片大全视频| 亚洲中文日韩欧美视频| 男女那种视频在线观看| 极品教师在线免费播放| 免费高清在线观看日韩| 日本 欧美在线| 中文字幕精品亚洲无线码一区 | 久久人妻av系列| 国产伦一二天堂av在线观看| 色哟哟哟哟哟哟| 在线视频色国产色| 变态另类丝袜制服| 狠狠狠狠99中文字幕| 午夜精品久久久久久毛片777| 成人亚洲精品av一区二区| 一级毛片女人18水好多| 亚洲精品在线观看二区| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲欧美精品永久| 婷婷丁香在线五月| 中文字幕久久专区| 国产精品亚洲美女久久久| 亚洲国产精品合色在线| 悠悠久久av| 两个人免费观看高清视频| 久久中文字幕人妻熟女| 免费在线观看黄色视频的| 欧美一级毛片孕妇| ponron亚洲| 亚洲第一青青草原| 男女之事视频高清在线观看| 成年免费大片在线观看| 最新美女视频免费是黄的| 亚洲专区国产一区二区| 欧美日本亚洲视频在线播放| 中文字幕人妻丝袜一区二区| 亚洲五月天丁香| 久久久久久亚洲精品国产蜜桃av| 可以在线观看毛片的网站| 白带黄色成豆腐渣| 久久欧美精品欧美久久欧美| 亚洲国产欧洲综合997久久, | 十分钟在线观看高清视频www| 国内精品久久久久精免费| 久久久久国产精品人妻aⅴ院| 18禁观看日本| 色综合亚洲欧美另类图片| 国产精品精品国产色婷婷| 欧美成狂野欧美在线观看| 国产成人啪精品午夜网站| 91在线观看av| 老汉色∧v一级毛片| 久久精品国产清高在天天线| 亚洲 国产 在线| 变态另类丝袜制服| 国产成人啪精品午夜网站| 久久婷婷成人综合色麻豆| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 亚洲专区国产一区二区| 香蕉丝袜av| 伦理电影免费视频| 国产三级在线视频| 香蕉国产在线看| 搡老熟女国产l中国老女人| 午夜a级毛片| 香蕉国产在线看| 好看av亚洲va欧美ⅴa在| 国产蜜桃级精品一区二区三区| 午夜精品在线福利| 一本综合久久免费| 国产欧美日韩精品亚洲av| 国产区一区二久久| 国产精品免费视频内射| 少妇裸体淫交视频免费看高清 | 啦啦啦韩国在线观看视频| 色在线成人网| 大香蕉久久成人网| 成人国产综合亚洲| a级毛片a级免费在线| 亚洲avbb在线观看| 欧美日韩精品网址| 国产一区二区三区视频了| 午夜老司机福利片| svipshipincom国产片| 日日摸夜夜添夜夜添小说| 午夜激情av网站| 啦啦啦免费观看视频1| 嫩草影院精品99| 精品无人区乱码1区二区| 日日爽夜夜爽网站| 亚洲一区二区三区色噜噜| 精品国产一区二区三区四区第35| 51午夜福利影视在线观看| 亚洲免费av在线视频| 91在线观看av| 神马国产精品三级电影在线观看 | 久热这里只有精品99| 女警被强在线播放| 亚洲国产高清在线一区二区三 | 国产精品美女特级片免费视频播放器 | 久久久久久久午夜电影| 成人午夜高清在线视频 | 欧美+亚洲+日韩+国产| 国产欧美日韩精品亚洲av| 午夜精品在线福利| 成人18禁高潮啪啪吃奶动态图| 亚洲色图av天堂| 亚洲第一av免费看| av视频在线观看入口| 成人国语在线视频| 亚洲午夜理论影院| 午夜免费成人在线视频| 精品久久久久久久毛片微露脸| 亚洲人成电影免费在线| 男男h啪啪无遮挡| 2021天堂中文幕一二区在线观 | 精品卡一卡二卡四卡免费| 757午夜福利合集在线观看| 国产日本99.免费观看| 久久久久精品国产欧美久久久| av电影中文网址| 人人妻人人澡欧美一区二区| 美女高潮到喷水免费观看| 18美女黄网站色大片免费观看| 法律面前人人平等表现在哪些方面| 熟妇人妻久久中文字幕3abv| 日韩免费av在线播放| 哪里可以看免费的av片| 国产熟女午夜一区二区三区| 日韩精品免费视频一区二区三区| 日本免费a在线| 欧美成人免费av一区二区三区| 999精品在线视频| 亚洲av熟女| 啦啦啦观看免费观看视频高清| 久久久久久九九精品二区国产 | 一边摸一边做爽爽视频免费| 草草在线视频免费看| 一夜夜www| 在线av久久热| 国产免费男女视频| 午夜福利18| 满18在线观看网站| 久久伊人香网站| 国产亚洲精品久久久久久毛片| 亚洲激情在线av| 黄频高清免费视频| 91成年电影在线观看| 在线看三级毛片| 在线十欧美十亚洲十日本专区| 18禁国产床啪视频网站| 香蕉国产在线看| 精品国产国语对白av| 一区福利在线观看| 可以在线观看的亚洲视频| 757午夜福利合集在线观看| 亚洲,欧美精品.| 人人妻人人澡人人看| 欧美国产精品va在线观看不卡| 午夜免费观看网址| 中文字幕精品亚洲无线码一区 | www.自偷自拍.com| av视频在线观看入口| 成人一区二区视频在线观看| 亚洲欧美日韩无卡精品| 国产色视频综合| 亚洲自偷自拍图片 自拍| 草草在线视频免费看| 操出白浆在线播放| av电影中文网址| 国产又爽黄色视频| 国产精华一区二区三区| 免费一级毛片在线播放高清视频| 999久久久精品免费观看国产| e午夜精品久久久久久久| 成人欧美大片| 2021天堂中文幕一二区在线观 | 欧美国产日韩亚洲一区| 99久久久亚洲精品蜜臀av| 亚洲人成77777在线视频| 校园春色视频在线观看| 久久久久久免费高清国产稀缺| 嫩草影院精品99| 校园春色视频在线观看| 悠悠久久av| 亚洲性夜色夜夜综合| 久久久久国产一级毛片高清牌| 中文亚洲av片在线观看爽| 十八禁网站免费在线| 欧美色欧美亚洲另类二区| 精品不卡国产一区二区三区| 国产午夜精品久久久久久| 国产区一区二久久| 欧美日韩瑟瑟在线播放| 免费高清视频大片| 18禁黄网站禁片午夜丰满| 熟妇人妻久久中文字幕3abv| 日韩欧美一区二区三区在线观看| 欧美日韩一级在线毛片| 成人永久免费在线观看视频| 99re在线观看精品视频| 免费看十八禁软件| 一级毛片女人18水好多| 韩国精品一区二区三区| 女性生殖器流出的白浆| 校园春色视频在线观看| 男人舔女人的私密视频| 啪啪无遮挡十八禁网站| 亚洲五月婷婷丁香| 一夜夜www| 在线播放国产精品三级| 国产爱豆传媒在线观看 | 88av欧美| 久久午夜亚洲精品久久| 一夜夜www| 久久久精品欧美日韩精品| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 丁香欧美五月| 在线观看午夜福利视频| 久久热在线av| 午夜视频精品福利| 久久中文字幕人妻熟女| 久久久国产成人精品二区| 亚洲国产欧美网| 黄色 视频免费看| 天堂√8在线中文| 国产精品美女特级片免费视频播放器 | 午夜久久久久精精品| 一区二区三区激情视频| 久久久国产精品麻豆| 一区二区三区精品91| 免费看美女性在线毛片视频| 最新在线观看一区二区三区| 欧美日本视频| 久久久久久亚洲精品国产蜜桃av| 成人永久免费在线观看视频| 国产真人三级小视频在线观看| 亚洲精品美女久久久久99蜜臀| 日韩欧美一区二区三区在线观看| 美女 人体艺术 gogo| 女性被躁到高潮视频| 国产成年人精品一区二区| aaaaa片日本免费| 白带黄色成豆腐渣| 久久精品国产清高在天天线| 99在线视频只有这里精品首页| 日本一本二区三区精品| 在线观看免费视频日本深夜| 一二三四社区在线视频社区8| 久久久国产成人免费| 亚洲va日本ⅴa欧美va伊人久久| 在线永久观看黄色视频| 国产欧美日韩一区二区精品| 美国免费a级毛片| 午夜免费观看网址| 免费看美女性在线毛片视频| 久久久国产精品麻豆| 久久久久国产一级毛片高清牌| 亚洲人成网站在线播放欧美日韩| 99精品欧美一区二区三区四区| 国产精品亚洲av一区麻豆| 久久久久久免费高清国产稀缺| 国产成人精品无人区| 怎么达到女性高潮| 搡老岳熟女国产| 国产熟女午夜一区二区三区| 女性被躁到高潮视频| 成人亚洲精品一区在线观看| 午夜成年电影在线免费观看| 日本五十路高清| 在线观看午夜福利视频| 真人做人爱边吃奶动态| 村上凉子中文字幕在线| 狠狠狠狠99中文字幕| 午夜久久久久精精品| 好男人电影高清在线观看| cao死你这个sao货| 国产激情偷乱视频一区二区| 此物有八面人人有两片| a级毛片在线看网站| 成人精品一区二区免费| 男女下面进入的视频免费午夜 | 黄色片一级片一级黄色片| 久久久久久九九精品二区国产 | 免费一级毛片在线播放高清视频| 一个人免费在线观看的高清视频| 中文字幕精品免费在线观看视频| 大型av网站在线播放| 777久久人妻少妇嫩草av网站| a级毛片在线看网站| 国产亚洲av高清不卡| 大香蕉久久成人网| 久久久久久九九精品二区国产 | 一级毛片高清免费大全| 国产精品精品国产色婷婷| 一本综合久久免费| 啪啪无遮挡十八禁网站| 国产视频一区二区在线看| 俺也久久电影网| 亚洲熟女毛片儿| 午夜精品在线福利| 亚洲专区中文字幕在线| 男人操女人黄网站| 日本成人三级电影网站| 久久午夜亚洲精品久久| 97碰自拍视频| 99久久综合精品五月天人人| 国产精品久久久人人做人人爽| 欧美久久黑人一区二区| 国产一区二区在线av高清观看| 热re99久久国产66热| 免费观看人在逋| 久久久久国产精品人妻aⅴ院| 久久热在线av| 国产av一区二区精品久久| 欧美大码av| 亚洲熟妇熟女久久| 久久久久久久久中文| 免费看美女性在线毛片视频| 成年版毛片免费区| 欧美乱妇无乱码| 嫩草影视91久久| 在线av久久热| 少妇被粗大的猛进出69影院| 国产在线精品亚洲第一网站| 久久久久久国产a免费观看| 韩国av一区二区三区四区| 日日干狠狠操夜夜爽| 一级黄色大片毛片| av福利片在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产片内射在线| 国产精品免费一区二区三区在线| 精品久久久久久久人妻蜜臀av| 免费高清在线观看日韩| 精品福利观看| 男人操女人黄网站| 亚洲精品一区av在线观看| 欧美另类亚洲清纯唯美| 亚洲激情在线av| 窝窝影院91人妻| 极品教师在线免费播放| 婷婷精品国产亚洲av| 久久九九热精品免费| 此物有八面人人有两片| 亚洲精品国产一区二区精华液| 宅男免费午夜| 午夜日韩欧美国产| 满18在线观看网站| 亚洲中文字幕一区二区三区有码在线看 | 午夜福利一区二区在线看| 观看免费一级毛片| 免费高清在线观看日韩| 成年版毛片免费区| 久久久国产成人免费| 人人妻人人澡欧美一区二区| 色尼玛亚洲综合影院| 一二三四社区在线视频社区8| 欧美日本视频| 午夜福利在线观看吧| 悠悠久久av| 久久久国产成人精品二区| 两个人免费观看高清视频| 精品久久久久久久久久免费视频| 露出奶头的视频| 亚洲av五月六月丁香网| 麻豆国产av国片精品| 村上凉子中文字幕在线| 国产高清videossex| 天堂√8在线中文| 一本综合久久免费| a在线观看视频网站| 18禁观看日本| 熟女少妇亚洲综合色aaa.| 国产精品免费一区二区三区在线| 精品国产国语对白av| 亚洲黑人精品在线| 久久婷婷成人综合色麻豆| 亚洲成a人片在线一区二区| 久久精品aⅴ一区二区三区四区| 免费在线观看日本一区| 精品免费久久久久久久清纯| 亚洲色图av天堂| 1024视频免费在线观看| 又紧又爽又黄一区二区| 免费高清视频大片| 久久久久九九精品影院| 欧美日韩精品网址| 国产伦人伦偷精品视频| 日韩国内少妇激情av| 俄罗斯特黄特色一大片| 国产亚洲精品一区二区www| 久久人妻福利社区极品人妻图片| 午夜日韩欧美国产| 亚洲av电影在线进入| 免费观看精品视频网站| 亚洲熟妇中文字幕五十中出| 国产在线精品亚洲第一网站| 国产精品免费一区二区三区在线| 九色国产91popny在线| 黄片小视频在线播放| 亚洲人成网站高清观看| 深夜精品福利| 精品一区二区三区av网在线观看| 精品乱码久久久久久99久播| 久久草成人影院| 高潮久久久久久久久久久不卡| 国产av一区二区精品久久| 99热这里只有精品一区 | 国产av又大| 午夜福利高清视频| 久99久视频精品免费| 韩国精品一区二区三区| 亚洲免费av在线视频| 18禁黄网站禁片免费观看直播| 国产激情久久老熟女| 日韩精品免费视频一区二区三区| 国产久久久一区二区三区| 亚洲激情在线av| 久久热在线av| 亚洲第一青青草原| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三| 制服人妻中文乱码| 国产av不卡久久| 长腿黑丝高跟| 国产精品 欧美亚洲| 国产精品香港三级国产av潘金莲| 亚洲欧美精品综合久久99| 午夜福利18| 夜夜夜夜夜久久久久| 女人高潮潮喷娇喘18禁视频| 亚洲成a人片在线一区二区| 色在线成人网| 亚洲av电影不卡..在线观看| av欧美777| 国产aⅴ精品一区二区三区波| 国产精品,欧美在线| 欧美久久黑人一区二区| 天堂√8在线中文| 一级毛片高清免费大全| 丝袜美腿诱惑在线| 无限看片的www在线观看| cao死你这个sao货| 在线观看66精品国产| 久久午夜亚洲精品久久| 日本免费一区二区三区高清不卡| 午夜老司机福利片| 在线视频色国产色| 国产久久久一区二区三区| 美女高潮喷水抽搐中文字幕| 成人欧美大片| 久久中文字幕一级| 久久久久国产精品人妻aⅴ院| 18禁裸乳无遮挡免费网站照片 | 亚洲第一av免费看| 日韩欧美国产在线观看| 啦啦啦 在线观看视频| 亚洲av片天天在线观看| av在线天堂中文字幕| 欧美不卡视频在线免费观看 | 国产爱豆传媒在线观看 | 日本精品一区二区三区蜜桃| 亚洲av中文字字幕乱码综合 | www日本在线高清视频| 1024视频免费在线观看| 久久久久久九九精品二区国产 | 久久人妻av系列| 国产黄a三级三级三级人| 国产精品亚洲美女久久久| av超薄肉色丝袜交足视频| 不卡一级毛片| 日韩一卡2卡3卡4卡2021年| 日韩精品免费视频一区二区三区| 人妻丰满熟妇av一区二区三区| 国产精品av久久久久免费| 黑人巨大精品欧美一区二区mp4| 久99久视频精品免费| 一a级毛片在线观看| 看免费av毛片| 少妇裸体淫交视频免费看高清 |