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

    基于數(shù)值仿真的X波段極化SAR海雜波統(tǒng)計建模與特性分析

    2021-11-11 06:14:24杜延磊
    關(guān)鍵詞:散射系數(shù)雜波海面

    杜延磊, 高 帆, 劉 濤, 楊 健,*

    (1. 中國科學(xué)院空天信息創(chuàng)新研究院, 北京 100101; 2. 清華大學(xué)電子工程系, 北京 100084;3. 海軍工程大學(xué)電子工程學(xué)院, 湖北 武漢 430033)

    0 引 言

    海雜波是制約極化合成孔徑雷達(dá)(synthetic aperture radar,SAR)海上目標(biāo)檢測和識別性能的主要因素之一,其物理機制復(fù)雜,且受雷達(dá)系統(tǒng)參數(shù)、觀測幾何參數(shù)及環(huán)境物理參數(shù)等諸多因素的影響,從而導(dǎo)致對其特征描述和抑制的難度大[1-2]。特別在中高海況的復(fù)雜海面條件下,對海雜波統(tǒng)計特性了解的欠缺成為在該背景下極化SAR海面目標(biāo)檢測、識別和跟蹤困難的主要問題。因此,深入了解極化SAR海雜波的統(tǒng)計分布特性有利于改善雜波抑制效果,提升對海面弱小目標(biāo)檢測性能,這對基于極化SAR的海面目標(biāo)探測預(yù)警、跟蹤識別、海上救援等軍事和民用都具有重要意義。

    極化SAR海雜波是隨機粗糙海面在各個極化通道對入射雷達(dá)電磁波的后向散射回波,因此,對極化SAR海雜波的統(tǒng)計建模和特性分析的前提是獲得大量不同極化狀態(tài)和參數(shù)條件下的海面回波數(shù)據(jù)。目前,獲得海雜波的手段主要分為兩類:一類是通過岸基或機載平臺對海雜波進(jìn)行實測獲得;另一類則主要基于海面形態(tài)建模和電磁散射求解的理論方法對海雜波進(jìn)行數(shù)值仿真[3]。前者的優(yōu)點在于所獲得的海雜波數(shù)據(jù)真實準(zhǔn)確,基于此所構(gòu)建的雜波模型對特定設(shè)備和特定應(yīng)用場景的針對性更強,實用性更好。因此,為了提升各型海用雷達(dá)裝備的海上目標(biāo)探測能力,國內(nèi)外開展了豐富多樣的雷達(dá)海雜波測量試驗,丁昊等對此進(jìn)行了較為系統(tǒng)全面的介紹[1,4]。然而,海雜波的測量試驗成本較高,且受限于試驗所用的設(shè)備、平臺性能以及所在的海域和時機,只能獲得特定場景和參數(shù)條件下的海雜波數(shù)據(jù)[1]。而海雜波的數(shù)值仿真則具有較好的靈活性和普適性,可以一定程度上彌補實測試驗的不足。這類方法首先利用海面形態(tài)建模方法生成隨機粗糙海面,進(jìn)而通過數(shù)值方法求解其散射方程(低頻類方法)[5]或通過引入近似假設(shè)獲得海面遠(yuǎn)場散射場的解析近似解(高頻類方法)[6-7]。目前,在實際的海雜波研究中,上述兩類海雜波獲取手段可以實現(xiàn)相互驗證、相互補充[1]。

    海雜波的統(tǒng)計分布模型是極化SAR海上目標(biāo)檢測的重要先驗知識[8-9]。對海雜波統(tǒng)計建模的研究經(jīng)歷了從高斯性到非高斯性、從單一極化到多極化、從純統(tǒng)計模型到物理/統(tǒng)計結(jié)合模型的發(fā)展過程[2,10]。較早提出的幾種適用于低分辨率單一極化海雜波的典型分布模型包括Rayleigh分布、Log-normal分布、Weibull分布、K分布等[2,11]。而隨著SAR技術(shù)的發(fā)展,早期的雜波分布模型難以準(zhǔn)確描述較高分辨率的極化SAR海雜波,尤其是其較長的拖尾表現(xiàn)出顯著的非高斯性[2]。目前,復(fù)合高斯模型(又稱乘積模型或紋理模型)是最為常用的一類極化SAR海雜波統(tǒng)計分布模型,其假設(shè)雷達(dá)雜波的散射向量可以分解為兩部分的乘積,一是描述目標(biāo)后向散射截面(radar cross section, RCS)空間變化的紋理特征,另一是服從復(fù)高斯分布的斑噪向量,并通過選擇不同的紋理分布模型,形成各種具體的極化SAR海雜波模型[10]。

    海雜波分布模型的參數(shù)估計也是其統(tǒng)計建模的關(guān)鍵問題之一。傳統(tǒng)的參數(shù)估計方法主要包括最大似然估計方法和矩估計方法[11-12]。其中,利用最大似然估計參數(shù)時往往需要用到數(shù)值方法進(jìn)行迭代求解,計算效率低[12]。而矩估計方法的特征函數(shù)基于Fourier變換獲得,在分布模型較為復(fù)雜的時候難以甚至無法獲得解析形式的特征函數(shù),此外,基于矩的估計方法對數(shù)據(jù)敏感性也常常較高,參數(shù)估計的精度較低[13]。針對SAR海雜波的乘積模型,Nicolas通過利用Mellin變換替換Fourier變換來構(gòu)建新的特征函數(shù)和累積量生成函數(shù),從而在此理論框架下構(gòu)建Mellin類統(tǒng)計量,并基于此進(jìn)行海雜波模型參數(shù)估計,具有計算簡單、精度高、通用性好的特點[14-15]。

    目前,對海雜波特性的研究主要利用實測數(shù)據(jù)的建模分析,而基于具有一定普適性的數(shù)值仿真建模研究相對較少,有限的相關(guān)研究存在海面仿真尺寸較小、極化信息不足的問題,且對所獲得雜波數(shù)據(jù)的統(tǒng)計建模和分析方法較為陳舊。鑒于此,本文利用成熟的Apel海浪譜模型,結(jié)合Monte Carlo方法仿真生成不同風(fēng)速下的大尺寸二維粗糙海面,進(jìn)而利用結(jié)合了物理光學(xué)(physical optics, PO)模型和彈跳射線法(shooting and bouncing ray,SBR)的電磁散射模型對多極化多角度的海雜波進(jìn)行仿真,最后基于Mellin類統(tǒng)計量對仿真結(jié)果進(jìn)行統(tǒng)計建模和特性分析。

    1 極化SAR海雜波數(shù)值仿真模型

    1.1 二維粗糙海面建模

    粗糙海面的空間形態(tài)建模是開展極化SAR海雜波數(shù)值仿真和分析的基礎(chǔ)。目前,海面形態(tài)建模最常用的方法是基于海浪譜模型的Monte Carlo方法,其假設(shè)隨機粗糙海面由無數(shù)的諧波疊加而成,從而可以通過對海浪譜先濾波再進(jìn)行逆Fourier變換獲得海面的高度起伏,因此,該方法也被稱為線性濾波方法[5,16]?;诖?二維粗糙海面高度f(x,y)可以表達(dá)為如下的Fourier級數(shù)展開:

    (1)

    式中:Lx和Ly分別為二維海面的長寬,bmn為高斯隨機變量。由于實際的海面仿真中,海面的仿真尺寸和剖分密度是有限的,其決定了對海浪譜濾波的上下限,并具有如下關(guān)系[17-18]:

    (2)

    (3)

    此外,根據(jù)式所示的海面高度Fourier級數(shù)展開特性,可以獲得其系數(shù)bmn的周期性和共軛對稱的特性,即

    b(m+M)(n+N)=bmn

    (4)

    (5)

    進(jìn)而結(jié)合粗糙面協(xié)方差的定義可以推導(dǎo)獲得bmn為[5]

    (6)

    圖1所示為利用Apel譜和Monte Carlo方法生成的風(fēng)速為6 m/s下的二維粗糙海面。

    圖1 利用Apel譜和Monte Carlo方法生成的6 m/s風(fēng)速下二維粗糙海面Fig.1 Two-dimensional rough ocean surface at 6 m/s generated using Apel spectrum and Monte Carlo method

    1.2 海面電磁散射計算

    基于Monte Carlo方法生成的二維粗糙海面,采用理論的電磁散射方法對海面雜波信號進(jìn)行仿真。現(xiàn)有常用的海面電磁散射計算方法大致可以分為低頻方法和高頻方法兩大類[11]。低頻方法主要是以矩量法[5]、時域有限差分法[22]、有限元法[23]等為代表的數(shù)值方法,其通過數(shù)值求解粗糙面問題的散射場,具有精度高、適用性廣的優(yōu)點,但由于其計算量巨大,實際中只能在較低的頻段使用[24]。而高頻方法在粗糙面散射問題的求解中通過引入一定的近似假設(shè),從而獲得散射場的解析形式表達(dá),可以在特定應(yīng)用范圍內(nèi)獲得較為準(zhǔn)確的結(jié)果,且求解速度快。

    由于本文主要關(guān)注3~5級海況下海面X波段的極化SAR海雜波特性,在此條件下,二維粗糙海面的尺寸和起伏都遠(yuǎn)遠(yuǎn)大于入射波波長,高精度的低頻數(shù)值方法須求解巨量的未知數(shù),對計算資源的消耗巨大,難以滿足三維問題下的海雜波仿真需求[24]。因此,本文采用高頻方法PO方法和SBR方法以分別考慮粗糙海面的單次散射和多次散射,從而開展極化海雜波的仿真。此外,由于海水具有較大的介電常數(shù),對海面模型采用理想電導(dǎo)體假設(shè)[25]。

    考慮如圖2所示的海面散射場景,入射電磁場Ei、Hi照射到海面一小面元s′,并形成散射場Es、Hs,并定義以海面模型質(zhì)心為原點的三維坐標(biāo)系O′X′Y′Z′。根據(jù)Maxwell方程組,散射場可由Stratton-Chu方程表達(dá)為[7]

    圖2 海面散射幾何示意圖Fig.2 Geometric configuration of ocean scattering

    (7)

    由于Stratton-Chu方程包含面積分項,無法獲得精確的解析解,因此需要引入一定的近似假設(shè),對方程進(jìn)行簡化求解。為求解Stratton-Chu方程,PO法引入3個近似假設(shè)[7]。① 高頻近似:假設(shè)入射波波長遠(yuǎn)小于目標(biāo)尺寸,從而近似認(rèn)為入射波只具備粒子性,忽略其在粗糙面棱邊及邊緣繞射等電磁波波動性,且認(rèn)為總的散射場僅由粗糙面上被入射波照射到的小面元散射場疊加而來,即,被遮擋的面元對總散射場無貢獻(xiàn)。② 遠(yuǎn)場近似:假設(shè)散射場觀測點遠(yuǎn)離目標(biāo)粗糙面,且二者距離也遠(yuǎn)大于粗糙面尺寸,即|r|?|r′|。③ 切平面近似:假設(shè)目標(biāo)表面任意一點及其附近表面曲率半徑比波長大得多,可以近似為平面。

    基于上述三個近似假設(shè),PO法可給出單站條件下(入射波方向i和散射波方向s相反,即i=-s)粗糙面元s′的后向散射場為

    (8)

    (9)

    式中:面積分項可以利用斯托克斯公式進(jìn)行計算,得到解析表達(dá)式為[26]

    (10)

    (11)

    式中:l1,l2,l3分別表示三角面元s′的3條邊矢量。

    由式(9)可以發(fā)現(xiàn),后向散射場的極化方向與入射場一致,因此,PO法無法給出交叉極化下的散射系數(shù)。鑒于此,本文進(jìn)一步采用SBR方法來計算粗糙面的多次散射。

    (12)

    類似的,根據(jù)幾何關(guān)系,可以獲得第n次反射的入射波方向為

    (13)

    式中:I為單位矩陣?;诖_定的反射路徑和PO法,可以獲得單站條件下,最終的SBR多次散射場為[27]

    (14)

    基于上述針對海面單次散射和多次散射計算獲得單站散射場,最終接收天線端的總散射場表示為

    (15)

    (16)

    式中:A表示仿真海面像元內(nèi)雷達(dá)波照射面積;θi為雷達(dá)入射角。為了仿真分析極化SAR幅度圖的海雜波統(tǒng)計特性[15,28],本文后續(xù)主要采用如下定義的后向散射系數(shù)的幅度進(jìn)行統(tǒng)計建模:

    (17)

    對于極化SAR系統(tǒng)而言,入射電場和散射電場也可由散射矩陣建立如下關(guān)系[8]:

    (18)

    式中:下標(biāo)v和h表示水平極化和垂直極化;S為散射矩陣,也可以表達(dá)為向量形式為

    s=vec(ST)=[Shh,Shv,Svh,Svv]T∈Cd

    (19)

    式中:vec(·)表示矩陣向量化操作;(·)T表示轉(zhuǎn)置;d為向量維度。

    2 基于Mellin類統(tǒng)計量的雜波建模

    2.1 極化SAR海雜波的典型乘積模型

    基于前文所介紹的海雜波數(shù)值仿真方法,可以計算獲得成像單元內(nèi)極化SAR散射向量s,其協(xié)方差矩陣定義為

    Σ=E{ssH}

    (20)

    式中:E{·}表示期望;(·)H表示Hermitian(轉(zhuǎn)置共軛)操作。對于視數(shù)為L的多視處理數(shù)據(jù)而言,多視極化協(xié)方差矩陣C為[15]

    (21)

    式中:Xf0表示矩陣正定;C表示復(fù)數(shù)集。

    而隨著雷達(dá)技術(shù)的發(fā)展,SAR圖像的空間分辨率得到了顯著的提升,傳統(tǒng)的Gaussian統(tǒng)計模型越來越難以滿足對高分辨率、非均勻海雜波的描述,尤其是大量數(shù)據(jù)表明海雜波分布的長拖尾表現(xiàn)出顯著的非高斯特性[2]。因此,諸多學(xué)者提出了一類復(fù)合Gaussian海雜波統(tǒng)計分布模型,又稱為乘積模型或紋理模型。這類模型假設(shè)雷達(dá)雜波的散射向量可以分解為兩部分的乘積,一是描述目標(biāo)后向散射截面空間變化的紋理特征,另一是服從復(fù)高斯分布的斑噪向量,并通過選擇不同的紋理分布模型,形成各種具體的極化SAR海雜波模型[10]。

    假設(shè)SAR海雜波各極化通道的紋理變化特性相同,則分辨單元內(nèi)的散射向量k可以寫作:

    (22)

    式中:τ為紋理參數(shù),其均值為1;s與前文定義一致,表示斑噪向量。對于多視數(shù)據(jù),其極化協(xié)方差矩陣表示為

    (23)

    對于給定的τ,散射向量k和極化協(xié)方差矩陣C的概率密度函數(shù)(probability density function, PDF)分別為[10]

    (24)

    (25)

    式中:tr(·)表示矩陣的跡;Γd(L)為歸一化因子,定義為

    (26)

    式中:Γ(·)表示gamma函數(shù)。

    當(dāng)紋理分布p(τ)已知時,海雜波分辨單元內(nèi)散射向量k和極化協(xié)方差矩陣C的PDF為

    (27)

    (28)

    目前,為了描述不同場景的極化SAR數(shù)據(jù)的雜波信號,學(xué)者們通過引入不同的紋理分布模型提出了多種極化SAR雜波的復(fù)合Gaussian統(tǒng)計模型。本文總結(jié)了其中較為常用的幾種典型雜波分布統(tǒng)計模型,包括K分布、G0分布、Kummer-U分布、W分布、M分布等,表1給出了這幾種分布模型紋理參量τ、散射向量k和極化協(xié)方差矩陣C的PDF。特別地,對于Wishart分布而言,可以認(rèn)為其紋理模型為常數(shù),即p(τ)=δ(τ-1)。

    表1中,式(32)中的Kν表示ν階第二類修正Bessel函數(shù),式(34)中U(a,b,z)為第二類合流超幾何函數(shù),式(35)中W為Whittaker W函數(shù),其為合流超幾何函數(shù)U(a,b,z)的函數(shù),式(36)中M為第一類合流超幾何函數(shù),亦即Kummer-M函數(shù)[10]。

    2.2 Mellin類統(tǒng)計量的定義與計算

    對于所獲得的一組雜波數(shù)據(jù),其統(tǒng)計建模一般來說包括兩個部分,一是確定其分布模型,二是對所確定分布模型的參數(shù)進(jìn)行估計。在極化SAR海雜波研究中,在經(jīng)過早期利用簡單的Rayleigh分布、Log-normal分布、Weibull分布等純統(tǒng)計模型的發(fā)展之后,越來越多的研究認(rèn)為中低分辨率條件下,極化SAR海雜波可以用K分布來描述,并基于此發(fā)展了大量K分布參數(shù)估計方法。其中,傳統(tǒng)的參數(shù)估計方法主要包括最大似然估計方法和矩估計方法[11-12]。而利用最大似然估計參數(shù)時往往需要使用數(shù)值方法進(jìn)行迭代求解,計算效率低。因此,基于矩的參數(shù)估計方法得到了廣泛地使用。

    矩估計方法的數(shù)學(xué)理論基礎(chǔ)是:對于給定隨機變量的分布函數(shù),其特征函數(shù)可以表示為關(guān)于矩的級數(shù)展開形式,當(dāng)該分布所有矩存在且有限,且其特征函數(shù)的冪級數(shù)展開在靠近原點處收斂,則該分布函數(shù)可由其特征函數(shù)唯一表示[15]。然而,矩估計方法的特征函數(shù)是基于Fourier變換獲得,在分布模型較為復(fù)雜的時候難以甚至無法獲得解析形式的特征函數(shù),此外,基于矩的估計方法對數(shù)據(jù)敏感性也常常較高,參數(shù)估計的精度較低[28]。因此,本文采用基于Mellin類統(tǒng)計量的方法來估計雜波分布模型的參數(shù)。

    不同于傳統(tǒng)參數(shù)估計方法的特征函數(shù)基于Fourier變換,Mellin類統(tǒng)計量的定義與計算采用Mellin變換。對于復(fù)矩陣元函數(shù)f(Z),其Mellin變換定義為[15]

    (29)

    基于Mellin變換,矩陣變量Z的Mellin類特征函數(shù)(characteristic function, CF)φZ和Mellin類累積量生成函數(shù)(cumulant generation function, CGF)φZ定義為

    φZ(χ)=E{Zχ-d}=M{pZ(Z)}(χ)=

    φZ(χ)=ln[φZ(χ)]

    (30)

    進(jìn)一步地,根據(jù)矩量和累積量的定義,Z的v階Mellin類矩μv(Z)和v階Mellin類累積量υv(Z)分別定義為

    (31)

    由(31)式可以發(fā)現(xiàn)μv(Z)=E{lnZ},因此Mellin類矩也被稱為對數(shù)矩,類似的,Mellin類累積量也被稱為對數(shù)累積量。當(dāng)Z的所有對數(shù)矩和對數(shù)累積量存在時,有

    (37)

    對數(shù)矩量和對數(shù)累積量之間可以由Faadi Bruno公式的組合形式表示為[29]

    (38)

    從而可獲得Z的前三階Mellin類累積量為

    (39)

    2.3 基于Mellin類統(tǒng)計量的雜波模型參數(shù)估計

    為了使用Mellin類統(tǒng)計量對海雜波模型參數(shù)進(jìn)行估計,我們需要推到獲得二者之間的嚴(yán)格關(guān)系。假設(shè)乘積模型的極化協(xié)方差矩陣寫作:

    (40)

    (41)

    (χ-d)ln|Σ|-(χ-d)dlnL

    (42)

    (43)

    進(jìn)一步地,對于紋理分量T,其Mellin類特征函數(shù)為

    (44)

    結(jié)合Mellin統(tǒng)計量定義和Faá di Bruno公式可得:

    (45)

    由于隨機矩陣變量X=UV的Mellin類特征函數(shù)和統(tǒng)計量有如下性質(zhì)[14]:

    (46)

    因此,乘積模型極化協(xié)方差矩陣C的Mellin類累積量為

    (47)

    代入式(43)最終有:

    (48)

    式(48)給出了極化SAR雜波的Mellin類累積量與乘積模型參數(shù)之間的關(guān)系,基于此,我們可以利用式(38)計算雜波數(shù)據(jù)的若干階累積量,從而可以對雜波分布模型參數(shù)進(jìn)行估算。對于式(48)中紋理模型的Mellin類統(tǒng)計量,表2給出了K分布、G0分布、Kummer-U分布、W分布、M分布等乘積模型的紋理分布對數(shù)累積量。

    表2 典型紋理分布模型的對數(shù)累積量Table 2 Log-cumulants of typical texture distribution models

    此外,在理論上,這5種典型乘積模型可以嚴(yán)格地劃分二階和三階Mellin類累積量所形成的平面空間,如圖3所示[9]。因此,利用樣本數(shù)據(jù)的Mellin類統(tǒng)計量還可以判斷雜波數(shù)據(jù)所服從的分布。

    圖3 典型雜波乘積模型在二階/三階Mellin類累積量空間的分布Fig.3 Typical product models in the space of 2nd/3rd Mellin-kind cumulants

    3 仿真結(jié)果與分析

    利用前文介紹的海面形態(tài)建模和雷達(dá)散射計算高頻方法,我們對海面X波段全極化海雜波進(jìn)行數(shù)值仿真,并基于Mellin類統(tǒng)計量方法開展海雜波統(tǒng)計建模和特性分析。

    由于本文旨在對3~5級海況下海面X波段(9.6 GHz)高分辨率極化海雜波進(jìn)行仿真建模和統(tǒng)計分析,因此,我們首先基于Apel海浪譜和Monte Carlo方法分別生成了6 m/s、9 m/s和12 m/s風(fēng)速下的二維粗糙海面,如圖4(a)所示。所生成海面的尺寸為128 m×128 m (4 096λ×4 096λ,λ為入射波波長),其中,為了仿真海面大尺度海浪,生成真實的中高海況海面,在海面仿真的過程中,本文使用了包含譜峰的海浪能量譜,這使得生成的海面包含了所有的大尺度風(fēng)浪,同時也使得仿真海面尺寸范圍較大,遠(yuǎn)遠(yuǎn)大于現(xiàn)有的低頻數(shù)值仿真方法可以實現(xiàn)的大小。尤其在仿真12 m/s風(fēng)速下海面時,為了包含海浪譜的譜峰,則生成的海面尺寸達(dá)到256 m×256 m。在后續(xù)的海面電磁散射計算中,為了與其他兩個風(fēng)速海面具有可比性,也為了減少運算量,我們對12 m/s風(fēng)速海面進(jìn)行裁剪,獲得尺寸為128 m×128 m的海面模型。由圖4(a)可以看出,隨著風(fēng)速增大,海面的粗糙起伏也逐漸變大,在本文的案例中,所生成的6 m/s、9 m/s和12 m/s風(fēng)速下的海面均方根高度分別為6.13λ、12.80λ和19.88λ。此外,圖4(a)中,風(fēng)向由y軸正坐標(biāo)指向負(fù)坐標(biāo)。為了保證對海面刻畫以及電磁散射仿真的精度,并考慮到海面小尺度粗糙度對雷達(dá)后向散射的影響,本文所仿真海面的剖分密度達(dá)到4個采樣點/λ。

    圖4 不同風(fēng)速下的二維粗糙海面和對應(yīng)的雷達(dá)后向散射系數(shù)Fig.4 Rough ocean surfaces at various wind speeds and corresponding radar backscattering coefficients

    基于所生成的不同風(fēng)速二維粗糙海面,本文利用所建立的海面電磁散射高頻仿真方法開展不同雷達(dá)入射角和方位角條件下X波段全極化海雜波數(shù)值仿真計算。在實際的計算中,根據(jù)SAR的成像空間分辨率概念,我們將海面剖分為64×64個大小為2 m×2 m的成像單元。圖4(b)所示為入射角為20°、雷達(dá)方位角為逆風(fēng)方向條件下,本文所計算的不同風(fēng)速下二維海面X波段的VV極化雷達(dá)后向散射系數(shù)。由圖可以發(fā)現(xiàn),隨著海面風(fēng)速的增大,海面的雷達(dá)后向散射系數(shù)也增大,這與海面粗糙度增大,海面雷達(dá)散射增強的基本常識相符合。此外,考慮到雷達(dá)觀測方向為逆風(fēng)方向(沿y軸負(fù)坐標(biāo)指向正坐標(biāo)),由圖4(a)和圖4(b)可以發(fā)現(xiàn),海面迎向入射波的面元后向散射系數(shù)較大,這與粗糙面雷達(dá)散射特性相一致。

    3.1 仿真海雜波的后向散射系數(shù)與分布驗證

    在開展海雜波的建模和分析之前,我們需要對仿真海雜波的后向散射系數(shù)和統(tǒng)計分布準(zhǔn)確性進(jìn)行驗證。本文選擇X波段廣泛使用的并由實測數(shù)據(jù)擬合獲得的地球物理模式函數(shù)XMOD-2對仿真海雜波的同極化后向散射系數(shù)進(jìn)行驗證[30]。對于交叉極化,受限于目前較少的X波段交叉極化實測數(shù)據(jù),我們采用Voronovich和Zavorotny在文獻(xiàn)中給出的二階小斜率近似模型(SSA-2)結(jié)果和測量數(shù)據(jù)對仿真海雜波的交叉極化后向散射系數(shù)進(jìn)行驗證[31-32]。圖5所示為仿真海雜波后向散射系數(shù)與XMOD-2、SSA-2模型結(jié)果及實測數(shù)據(jù)的對比。

    圖5 利用其他模型和實測數(shù)據(jù)對仿真海雜波后向散射系數(shù)精度進(jìn)行驗證Fig.5 Accuracy validation of simulated sea clutter using other models and measurement

    由圖5可以看出,在同極化條件下,仿真海雜波的后向散射系數(shù)與XMOD-2模型有較好的吻合,最大誤差不超過2 dB,這表明本文采用的海面雷達(dá)散射高頻方法在同極化條件下對X波段海雜波具有較高的仿真精度。而對于交叉極化,盡管在高風(fēng)速情況下仿真海雜波后向散射系數(shù)與SSA-2模型結(jié)果較為接近,但在低風(fēng)速情況下,仿真結(jié)果則與參考數(shù)據(jù)有相對較大的誤差,最大可達(dá)約10 dB??紤]到海面交叉極化后向散射能量遠(yuǎn)遠(yuǎn)低于同極化,且本文主要關(guān)注中高海況下海雜波的統(tǒng)計分布及定性的特性分析,這一誤差大小對于交叉極化仿真是可以接受的。

    為了驗證數(shù)值仿真所得到海雜波的分布特性準(zhǔn)確性,我們利用實測的海雜波分布與仿真海雜波分布進(jìn)行對比。本文采用由岸基X波段固態(tài)功放監(jiān)視/導(dǎo)航雷達(dá)于2021年1月6日15時6分利用凝視模式獲得的海面雜波數(shù)據(jù),該數(shù)據(jù)由海軍航空大學(xué)的“雷達(dá)對海探測數(shù)據(jù)共享計劃”發(fā)布,并包含了同步測量的氣象水文等輔助資料。獲得實測數(shù)據(jù)的雷達(dá)中心頻率為9.4 GHz,極化方式為HH極化,距離分辨率為6 m,入射角為80°,雷達(dá)觀測方位角為55.05°,其他的雷達(dá)系統(tǒng)參數(shù)在文獻(xiàn)[3,33]具體給出。此外,根據(jù)與同步測量的海洋氣象水文數(shù)據(jù)進(jìn)行時空匹配,可知所測海雜波數(shù)據(jù)對應(yīng)海面的風(fēng)速為9.57 m/s,風(fēng)向為61.73°,因此,考慮到雷達(dá)觀測方位角,雷達(dá)觀測方位相對于風(fēng)向約為6.68°。

    由于實測海雜波數(shù)據(jù)是未經(jīng)輻射定標(biāo)的散射強度結(jié)果,與仿真數(shù)據(jù)雷達(dá)后向散射系數(shù)無法進(jìn)行直接對比,因此,本文對仿真和實測海雜波數(shù)據(jù)的進(jìn)行了歸一化處理,并基于歸一化后的海雜波數(shù)據(jù)計算了其前三階對數(shù)累積量,圖6(a)所示為實測和仿真海雜波數(shù)據(jù)的二階/三階對數(shù)累積量對比。由圖中散點分布可以看出實測和仿真海雜波具有非常相近的分布特征,且較為服從Wishart分布,這表明在此條件下,所獲得的海雜波紋理特性并不顯著,考慮到本案例較大的雷達(dá)入射角,這一特性是較為合理的。

    圖6 仿真海雜波分布與實測海雜波分布的對比Fig.6 Comparison of distributions of simulated sea clutter and measured sea clutter

    基于圖6(a)結(jié)果,我們采用K分布模型對歸一化實測和仿真海雜波的雷達(dá)后向散射系數(shù)和散射場幅度的概率密度分布進(jìn)行擬合,并利用所計算的二階/三階對數(shù)累積量對模型參數(shù)進(jìn)行估計,最終仿真和實測海雜波分布的對比如圖6(b)所示??梢钥闯?仿真海雜波與實測海雜波的分布具有良好的匹配,證明了利用前文的海面雜波高頻數(shù)值仿真方法對海雜波分布仿真的準(zhǔn)確性,以及利用Mellin類統(tǒng)計量開展海雜波統(tǒng)計分布和特性分析的可行性。

    在本文的海雜波數(shù)值仿真中,盡管我們通過采用包含海浪譜譜峰考慮了所有大尺度風(fēng)浪的貢獻(xiàn),但是小尺度波浪的貢獻(xiàn)則由于剖分密度的限制進(jìn)行了一定的截斷。由式可知,4個采樣點/λ的剖分密度表明海面波數(shù)大于402 rad/m的小尺度波浪在仿真海面中被舍去。因此,對于本文采用的電磁散射方法和海面剖分密度是否足以描述海浪細(xì)節(jié)對海雜波的貢獻(xiàn),尤其在5級高海況條件下,需要進(jìn)一步進(jìn)行驗證。圖7給出了5級海況下,分別利用4個采樣點/λ和8個采樣點/λ仿真海面獲得的海雜波的二階/三階對數(shù)累積量及分布對比。特別地,這里海雜波的仿真像元大小為1 m×1 m。由圖7(a)可以看出,采用不同剖分密度海面獲得的極化SAR仿真海雜波具有相近的二階/三階對數(shù)累積量,且兩組海雜波的后向散射系數(shù)平均差異在同極化下可以忽略,在交叉極化下僅約0.72 dB。在圖7(b)中,由不同剖分密度海面獲得的海雜波在同極化和交叉極化都分別具有幾乎一致的分布。這表明海面較小尺度的波浪對海面雷達(dá)散射的貢獻(xiàn)較小,這也與我們前期的研究結(jié)論相符合[17],本文采用4個采樣點/λ的海面剖分密度可以滿足海雜波仿真中對海浪細(xì)節(jié)的描述。

    圖7 不同剖分密度海面獲得的仿真海雜波分布對比Fig.7 Comparison of distributions of simulated sea clutters from surfaces with different discretizing density

    3.2 X波段高分辨率極化SAR海雜波建模與分析

    在海雜波數(shù)值仿真的基礎(chǔ)上,我們基于Mellin類統(tǒng)計量構(gòu)建海雜波的分布模型,并分析了雷達(dá)入射角、海面風(fēng)速風(fēng)向?qū)ks波分布的影響。圖8所示為不同風(fēng)速下X波段同極化海雜波的分布,其中,仿真海雜波的雷達(dá)入射角為20°,風(fēng)向為逆風(fēng)方向。由圖8(a)首先可以看出,VV極化和HH極化海雜波的二階/三階對數(shù)累積量非常接近,這主要是因為雷達(dá)入射角較小,VV和HH極化的雷達(dá)后向散射系數(shù)較為接近,表明VV和HH極化海雜波有著相近的分布,尤其是二者的紋理分布接近,因此在圖8(b)只展示了VV極化結(jié)果。進(jìn)一步地,我們可以發(fā)現(xiàn)高分辨率的同極化X波段海雜波并不服從傳統(tǒng)的適用于中低分辨率海雜波的K分布模型,這說明在此條件下,海雜波的紋理特征不服從gamma分布。而在不同的風(fēng)速下,X波段同極化海雜波表現(xiàn)出不同的紋理分布,并最終使得其服從不同的分布模型。具體來說,6 m/s和9 m/s風(fēng)速下海雜波分別服從Kummer-U分布和M分布,但二者都比較接近G0分布,在12 m/s風(fēng)速下,X波段同極化海雜波服從M分布。基于海雜波所服從的分布模型,我們利用所計算的對數(shù)累積量對相應(yīng)的模型參數(shù)進(jìn)行了估計,并在圖8(b)中對比了不同風(fēng)速下海雜波幅度的概率密度分布。由圖8(b)可以看出,利用圖8(a)所指示的分布模型可以很好地擬合海雜波數(shù)據(jù)。隨著風(fēng)速的增大,海雜波分布中后向散射系數(shù)幅度較大部分占比逐漸升高,海雜波分布由陡逐漸變緩,這與物理事實相符合。

    圖8 不同風(fēng)速下X波段同極化海雜波的分布Fig.8 Distribution of X-band co-pol sea clutters at various wind speeds

    圖9展示了不同風(fēng)速下X波段交叉極化海雜波的分布,相關(guān)的海雜波仿真參數(shù)與圖8相同。類似于同極化,HV極化和VH極化海雜波的二階、三階對數(shù)累積量也非常接近,考慮到極化互易性原理,這是合理的,因此,在此處及后續(xù)的交叉極化海雜波分布模型中,我們僅展示HV極化。而與同極化不同的是,交叉極化海雜波的二階、三階對數(shù)累積量在不同風(fēng)速下都集中在Wishart分布附近,表明交叉極化條件下高分辨率X波段海雜波沒有顯著的紋理特征,因此,在高分辨率條件下,極化SAR海雜波在四個極化通道不適宜使用相同的紋理模型來描述,這與傳統(tǒng)針對中低分辨率海雜波在各極化通道采用相同紋理模型不一致。在圖9(b)中我們利用K分布擬合了不同風(fēng)速下交叉極化的海雜波分布,與同極化類似,隨著風(fēng)速的增大,海雜波分布由陡逐漸變緩。

    圖9 不同風(fēng)速下X波段交叉極化海雜波的分布Fig.9 Distribution of X-band cross-pol sea clutters at various wind speeds

    圖10所示為不同入射角下X波段同極化海雜波的分布,其中,仿真海雜波的海面風(fēng)速為9 m/s,風(fēng)向為逆風(fēng)方向。由圖10(a)可以看出,雷達(dá)入射角對高分辨率X波段同極化海雜波的分布有著顯著的影響,在入射角較小時,同極化海雜波有著顯著的紋理特征,且VV和HH極化的二階/三階對數(shù)累積量非常接近。而在雷達(dá)入射角較大時,海雜波的紋理特征顯著減弱,且二階/三階對數(shù)累積量在VV和HH極化間顯示出差異。在圖10(b)中我們分別利用M分布和K分布擬合雷達(dá)入射角為20°、40°、60°時的VV極化海雜波分布,可以看出所用模型可以很好地擬合海雜波分布,由于海面同極化后向散射能量隨著入射角增大而減小,圖中的分布模型也隨著入射角增大而逐漸變陡。

    圖10 不同入射角下X波段同極化海雜波的分布Fig.10 Distribution of X-band co-pol sea clutters at various incidence angles

    圖11所示為不同風(fēng)向下X波段同極化海雜波的分布,其中,仿真海雜波的雷達(dá)入射角為20°,海面風(fēng)速為9 m/s。由圖可以看出,該條件下海雜波的二階/三階對數(shù)累積量集中在G0分布附近,并在不同的風(fēng)向下稍有區(qū)別。在圖11(b)中,我們分別采用M分布、G0分布和Kummer-U分布對風(fēng)向為逆風(fēng)、側(cè)風(fēng)和45°時的海雜波分布進(jìn)行擬合,可以發(fā)現(xiàn)海雜波分布隨著風(fēng)向由側(cè)風(fēng)向逆風(fēng)變化逐漸變緩。

    圖11 不同風(fēng)向下X波段同極化海雜波的分布Fig.11 Distribution of X-band co-pol sea clutters in various wind directions

    圖12展示了雷達(dá)視數(shù)對X波段同極化海雜波分布的影響,其中,仿真海雜波的雷達(dá)入射角為20°,海面風(fēng)速為9 m/s,風(fēng)向為側(cè)風(fēng)方向。由圖12(a)可以看出,隨著雷達(dá)視數(shù)的增大,海雜波的紋理特征逐漸減弱,考慮到雷達(dá)多視處理的基本原理及其空間濾波的效果,海雜波的紋理特征被平滑是合理的。進(jìn)一步地,由圖12(b)也可以看出,多視處理使得海雜波的極大值和極小值減少,海雜波的分布逐漸向中間移動。

    圖12 不同雷達(dá)視數(shù)條件下X波段同極化海雜波的分布Fig.12 Distribution of X-band co-pol sea clutters under conditions with various radar looks

    4 結(jié) 論

    中高海況下極化SAR海雜波的統(tǒng)計建模是海面艦船、導(dǎo)彈及相關(guān)飛行器檢測預(yù)警的關(guān)鍵和難點問題。本文針對3~5級海況下X波段極化SAR海雜波開展數(shù)值仿真、統(tǒng)計建模與分析研究。首先利用成熟的Apel海浪譜模型和Monte Carlo方法建立了6 m/s、9 m/s和12 m/s風(fēng)速下二維粗糙海面模型,進(jìn)而采用高頻方法開展不同風(fēng)速、風(fēng)向、雷達(dá)入射角、極化方式下的X波段海雜波仿真計算。為了獲得準(zhǔn)確的海面極化散射系數(shù),我們采用PO方法和SBR方法以分別考慮粗糙海面的單次散射和多次散射?;诜抡婧ks波數(shù)據(jù),利用典型的乘積模型和Mellin類統(tǒng)計量開展高分辨率X波段極化SAR海雜波統(tǒng)計建模和特性分析。仿真結(jié)果表明:在2 m×2 m的高分辨率條件下,X波段同極化海雜波并不服從適用于中低分辨率海雜波的K分布,且在不同風(fēng)速風(fēng)向下呈現(xiàn)出不同的紋理特征;X波段交叉極化海雜波表現(xiàn)出較弱的紋理特征,服從Wishart分布。因此,X波段極化SAR海雜波在同極化和交叉極化條件下具有不同的紋理分布特征,這與傳統(tǒng)針對中低分辨率海雜波在各極化通道采用相同紋理模型不一致。此外,隨著入射角和雷達(dá)視數(shù)的增大,海雜波的紋理特征逐漸減弱。

    盡管本文基于高頻仿真方法和Mellin類統(tǒng)計量對X波段極化海雜波進(jìn)行仿真建模,給出了海面風(fēng)速、風(fēng)向、雷達(dá)入射角等雷達(dá)系統(tǒng)參數(shù)和環(huán)境參數(shù)對極化海雜波影響的初步分析結(jié)果,但受限于模型在交叉極化的仿真精度,相關(guān)的分析僅在定性層面展開,未給出定量化的結(jié)果。在后續(xù)工作中,可進(jìn)一步構(gòu)建精細(xì)化的海面模型,并考慮海浪尖峰的繞射等現(xiàn)象,提高模型在交叉極化的仿真精度,并在本文工作的基礎(chǔ)上增加仿真數(shù)據(jù),建立極化SAR海雜波統(tǒng)計模型與雷達(dá)系統(tǒng)參數(shù)和海洋環(huán)境參數(shù)之間的經(jīng)驗?zāi)J胶瘮?shù),服務(wù)海面極化SAR海面目標(biāo)檢測應(yīng)用。

    猜你喜歡
    散射系數(shù)雜波海面
    等離子體層嘶聲波對輻射帶電子投擲角散射系數(shù)的多維建模*
    STAR2000型空管一次雷達(dá)雜波抑制淺析
    海面床,輕輕搖
    北部灣后向散射系數(shù)的時空分布與變化分析
    第六章 邂逅“胖胖號”
    海面上的“一千座埃菲爾鐵塔”
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    相關(guān)廣義復(fù)合分布雷達(dá)海雜波仿真
    遙測遙控(2015年2期)2015-04-23 08:15:22
    一維帶限Weierstrass分形粗糙面電磁散射的微擾法研究
    久久国产精品影院| 国产精品综合久久久久久久免费| 日本成人三级电影网站| 亚洲专区国产一区二区| 777久久人妻少妇嫩草av网站| 免费一级毛片在线播放高清视频| 免费观看精品视频网站| 男男h啪啪无遮挡| 久久 成人 亚洲| 一区二区三区国产精品乱码| 9191精品国产免费久久| 搡老熟女国产l中国老女人| 亚洲午夜理论影院| 亚洲自偷自拍图片 自拍| 色老头精品视频在线观看| 一级a爱视频在线免费观看| 免费无遮挡裸体视频| 又紧又爽又黄一区二区| 精品一区二区三区视频在线观看免费| 老熟妇仑乱视频hdxx| 2021天堂中文幕一二区在线观 | 黄色视频不卡| 草草在线视频免费看| 欧美久久黑人一区二区| www.精华液| 亚洲国产中文字幕在线视频| 亚洲人成网站高清观看| 久久国产精品人妻蜜桃| 脱女人内裤的视频| 国产精品乱码一区二三区的特点| 亚洲五月色婷婷综合| 色在线成人网| 久久香蕉国产精品| 国产高清有码在线观看视频 | 色综合亚洲欧美另类图片| 午夜免费观看网址| 老汉色av国产亚洲站长工具| 亚洲人成网站高清观看| 欧洲精品卡2卡3卡4卡5卡区| 国产高清videossex| 怎么达到女性高潮| 日本熟妇午夜| 日韩精品青青久久久久久| 亚洲第一av免费看| 美女高潮喷水抽搐中文字幕| 欧美性猛交黑人性爽| 国产精品免费一区二区三区在线| 日本免费a在线| 好男人在线观看高清免费视频 | 最好的美女福利视频网| 久久久水蜜桃国产精品网| 一级片免费观看大全| 精品免费久久久久久久清纯| 禁无遮挡网站| 国产欧美日韩一区二区精品| 亚洲av电影不卡..在线观看| 欧美最黄视频在线播放免费| 麻豆成人av在线观看| avwww免费| 国产欧美日韩精品亚洲av| 欧美日韩亚洲综合一区二区三区_| 欧美人与性动交α欧美精品济南到| 国产激情久久老熟女| 黄片播放在线免费| 热re99久久国产66热| 黑人操中国人逼视频| 日本熟妇午夜| 女同久久另类99精品国产91| 国产爱豆传媒在线观看 | 国产蜜桃级精品一区二区三区| 国产精品 国内视频| 久久伊人香网站| 欧美精品亚洲一区二区| 黄色女人牲交| 欧美成人一区二区免费高清观看 | 免费在线观看亚洲国产| 精品卡一卡二卡四卡免费| 女人高潮潮喷娇喘18禁视频| 大型av网站在线播放| 法律面前人人平等表现在哪些方面| 欧美不卡视频在线免费观看 | 老鸭窝网址在线观看| 久久伊人香网站| 女同久久另类99精品国产91| 亚洲精品粉嫩美女一区| 首页视频小说图片口味搜索| 国产精品98久久久久久宅男小说| 亚洲久久久国产精品| 香蕉国产在线看| 50天的宝宝边吃奶边哭怎么回事| 女人高潮潮喷娇喘18禁视频| 国产精品久久久久久人妻精品电影| 熟女少妇亚洲综合色aaa.| 精品福利观看| 欧美激情 高清一区二区三区| 黄片大片在线免费观看| 成人精品一区二区免费| 亚洲精品久久成人aⅴ小说| 母亲3免费完整高清在线观看| 国产精品 欧美亚洲| 一级黄色大片毛片| 美女高潮到喷水免费观看| 99久久国产精品久久久| 国产伦一二天堂av在线观看| 色综合亚洲欧美另类图片| 成年人黄色毛片网站| 日韩大尺度精品在线看网址| 国产av在哪里看| 香蕉丝袜av| 午夜福利成人在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲欧美精品综合久久99| 女性被躁到高潮视频| av在线播放免费不卡| a级毛片在线看网站| 久久午夜综合久久蜜桃| 自线自在国产av| 丝袜美腿诱惑在线| 法律面前人人平等表现在哪些方面| 黑人巨大精品欧美一区二区mp4| 久久99热这里只有精品18| 一本久久中文字幕| 久久精品国产综合久久久| 天天躁狠狠躁夜夜躁狠狠躁| 视频在线观看一区二区三区| 亚洲国产精品久久男人天堂| 亚洲va日本ⅴa欧美va伊人久久| 日韩免费av在线播放| 欧美丝袜亚洲另类 | 久9热在线精品视频| 日本撒尿小便嘘嘘汇集6| 99久久综合精品五月天人人| 欧洲精品卡2卡3卡4卡5卡区| 久久久久免费精品人妻一区二区 | 国产午夜精品久久久久久| 青草久久国产| 欧美乱色亚洲激情| 一进一出好大好爽视频| 99在线视频只有这里精品首页| 色综合婷婷激情| 一级a爱片免费观看的视频| 亚洲精华国产精华精| 两个人看的免费小视频| 波多野结衣av一区二区av| 日本 欧美在线| 欧美久久黑人一区二区| 最好的美女福利视频网| 国产亚洲av嫩草精品影院| 99国产精品99久久久久| 久久精品国产亚洲av香蕉五月| 黄片大片在线免费观看| 午夜久久久久精精品| 在线观看66精品国产| 欧美黑人精品巨大| 99re在线观看精品视频| 在线观看免费午夜福利视频| 美女国产高潮福利片在线看| 天堂动漫精品| 婷婷丁香在线五月| 一边摸一边做爽爽视频免费| 黑丝袜美女国产一区| 亚洲男人的天堂狠狠| 动漫黄色视频在线观看| 51午夜福利影视在线观看| 国产片内射在线| 亚洲国产欧美日韩在线播放| 成人午夜高清在线视频 | 久久久国产精品麻豆| 欧美中文日本在线观看视频| 一级黄色大片毛片| 成在线人永久免费视频| 黄色丝袜av网址大全| 成人av一区二区三区在线看| 国产精品 国内视频| 精品国内亚洲2022精品成人| 老司机靠b影院| 精华霜和精华液先用哪个| 国产高清视频在线播放一区| 亚洲精品久久成人aⅴ小说| av在线天堂中文字幕| 18禁美女被吸乳视频| 亚洲成av人片免费观看| 成人18禁高潮啪啪吃奶动态图| 亚洲电影在线观看av| 美女大奶头视频| 欧美zozozo另类| 欧美成狂野欧美在线观看| 国产av一区在线观看免费| 精品一区二区三区av网在线观看| 免费无遮挡裸体视频| 无遮挡黄片免费观看| 久久久国产精品麻豆| 中文字幕人妻熟女乱码| 国产高清videossex| 一级a爱视频在线免费观看| 在线观看一区二区三区| 午夜亚洲福利在线播放| 最新在线观看一区二区三区| 好看av亚洲va欧美ⅴa在| 黄色女人牲交| 可以免费在线观看a视频的电影网站| 女同久久另类99精品国产91| 两个人视频免费观看高清| 国产成人啪精品午夜网站| 国产乱人伦免费视频| 91成人精品电影| 51午夜福利影视在线观看| 国产成人欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 成人亚洲精品av一区二区| 国产精品亚洲美女久久久| www日本黄色视频网| 久久香蕉国产精品| 国产av又大| 一级毛片女人18水好多| 在线观看www视频免费| 一进一出好大好爽视频| 一边摸一边抽搐一进一小说| 热re99久久国产66热| 黄色女人牲交| 女人被狂操c到高潮| 不卡一级毛片| 亚洲五月天丁香| 波多野结衣巨乳人妻| 久久亚洲精品不卡| 97超级碰碰碰精品色视频在线观看| 国产精品1区2区在线观看.| 国产午夜精品久久久久久| 91av网站免费观看| 精品少妇一区二区三区视频日本电影| 99国产综合亚洲精品| 一本精品99久久精品77| 亚洲国产精品999在线| 午夜福利在线观看吧| 黑人操中国人逼视频| 宅男免费午夜| 亚洲国产精品久久男人天堂| 丁香欧美五月| 国内精品久久久久久久电影| 国产精品综合久久久久久久免费| 嫁个100分男人电影在线观看| 中文资源天堂在线| 在线天堂中文资源库| 午夜成年电影在线免费观看| 国产精品综合久久久久久久免费| 午夜福利视频1000在线观看| 免费无遮挡裸体视频| 亚洲av电影在线进入| 正在播放国产对白刺激| 一区福利在线观看| 最新在线观看一区二区三区| 日韩欧美国产一区二区入口| 在线观看免费午夜福利视频| 欧美日韩黄片免| 欧美成人性av电影在线观看| 91在线观看av| 精华霜和精华液先用哪个| 亚洲精品美女久久久久99蜜臀| 一进一出抽搐gif免费好疼| 黄片小视频在线播放| 男女之事视频高清在线观看| www.精华液| 激情在线观看视频在线高清| 欧美日本亚洲视频在线播放| 欧美性猛交黑人性爽| 桃色一区二区三区在线观看| 欧美+亚洲+日韩+国产| 久久久久九九精品影院| 亚洲色图av天堂| 久久中文字幕人妻熟女| 男男h啪啪无遮挡| 黑丝袜美女国产一区| 怎么达到女性高潮| 男女床上黄色一级片免费看| 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 成人手机av| 亚洲av五月六月丁香网| 亚洲七黄色美女视频| 精品卡一卡二卡四卡免费| 亚洲国产欧洲综合997久久, | 麻豆成人av在线观看| 国产黄色小视频在线观看| 窝窝影院91人妻| 欧美激情 高清一区二区三区| 亚洲av电影不卡..在线观看| 美女免费视频网站| 狂野欧美激情性xxxx| 久久精品人妻少妇| 亚洲色图 男人天堂 中文字幕| 国产乱人伦免费视频| 制服丝袜大香蕉在线| 一区二区三区国产精品乱码| 欧美日韩中文字幕国产精品一区二区三区| 最新美女视频免费是黄的| 久久狼人影院| 亚洲第一欧美日韩一区二区三区| 人人妻人人澡欧美一区二区| 在线看三级毛片| 少妇 在线观看| 十八禁网站免费在线| 欧美中文综合在线视频| 无遮挡黄片免费观看| 非洲黑人性xxxx精品又粗又长| 精品一区二区三区视频在线观看免费| 在线观看免费日韩欧美大片| 黑人巨大精品欧美一区二区mp4| 满18在线观看网站| 日韩国内少妇激情av| 天天添夜夜摸| 香蕉久久夜色| 欧美性猛交╳xxx乱大交人| 国产99久久九九免费精品| 1024香蕉在线观看| 制服人妻中文乱码| 在线永久观看黄色视频| 亚洲国产中文字幕在线视频| 一区二区日韩欧美中文字幕| 色尼玛亚洲综合影院| 午夜福利欧美成人| 欧美一级毛片孕妇| 男人舔女人下体高潮全视频| 可以免费在线观看a视频的电影网站| 亚洲精品在线观看二区| x7x7x7水蜜桃| 最好的美女福利视频网| 女人被狂操c到高潮| 婷婷精品国产亚洲av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 麻豆久久精品国产亚洲av| 99riav亚洲国产免费| 黄色成人免费大全| 国产成人精品久久二区二区91| 好看av亚洲va欧美ⅴa在| 变态另类丝袜制服| 首页视频小说图片口味搜索| 亚洲国产精品合色在线| 日韩欧美国产在线观看| 侵犯人妻中文字幕一二三四区| 91麻豆精品激情在线观看国产| 无人区码免费观看不卡| 欧美国产精品va在线观看不卡| 欧美黑人精品巨大| 在线观看免费午夜福利视频| 欧美激情久久久久久爽电影| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区 | 精品久久久久久,| av片东京热男人的天堂| 2021天堂中文幕一二区在线观 | 国产精品久久久人人做人人爽| 亚洲男人的天堂狠狠| 久久狼人影院| 日本在线视频免费播放| 色综合亚洲欧美另类图片| 欧美日韩亚洲国产一区二区在线观看| 国产精品综合久久久久久久免费| 亚洲av美国av| 免费在线观看完整版高清| av福利片在线| 亚洲最大成人中文| 亚洲成av片中文字幕在线观看| 国产亚洲欧美98| 超碰成人久久| 色综合亚洲欧美另类图片| 99精品久久久久人妻精品| 国产免费男女视频| 男人舔女人下体高潮全视频| 别揉我奶头~嗯~啊~动态视频| 久久精品亚洲精品国产色婷小说| 岛国视频午夜一区免费看| 久久久久久亚洲精品国产蜜桃av| 一进一出抽搐gif免费好疼| 黄片大片在线免费观看| 色婷婷久久久亚洲欧美| 成人免费观看视频高清| 国产熟女xx| 黑丝袜美女国产一区| 中文资源天堂在线| 日日干狠狠操夜夜爽| 久久久久国产精品人妻aⅴ院| 成人午夜高清在线视频 | 午夜亚洲福利在线播放| 成年版毛片免费区| 十八禁人妻一区二区| 国产亚洲av高清不卡| 男人的好看免费观看在线视频 | 99热只有精品国产| 亚洲在线自拍视频| 久久久久久久午夜电影| 国产亚洲精品av在线| 黑丝袜美女国产一区| 成人永久免费在线观看视频| 两个人视频免费观看高清| 婷婷亚洲欧美| 久久精品成人免费网站| 在线观看www视频免费| 精品欧美国产一区二区三| 老司机靠b影院| 免费在线观看黄色视频的| 男人舔奶头视频| 两人在一起打扑克的视频| 无遮挡黄片免费观看| 特大巨黑吊av在线直播 | 脱女人内裤的视频| 日韩视频一区二区在线观看| 免费在线观看视频国产中文字幕亚洲| av超薄肉色丝袜交足视频| 国产亚洲av高清不卡| 欧美成狂野欧美在线观看| ponron亚洲| 亚洲人成伊人成综合网2020| 国产又爽黄色视频| 国产精品永久免费网站| 亚洲国产欧洲综合997久久, | 亚洲av熟女| 亚洲av五月六月丁香网| 国产精华一区二区三区| 好男人电影高清在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲国产精品999在线| 国产精品久久电影中文字幕| 91麻豆精品激情在线观看国产| 91老司机精品| 色综合站精品国产| 天天一区二区日本电影三级| 丝袜人妻中文字幕| 18禁观看日本| 三级毛片av免费| 91成人精品电影| 午夜福利成人在线免费观看| 亚洲,欧美精品.| 我的亚洲天堂| 久久精品国产清高在天天线| 亚洲欧洲精品一区二区精品久久久| 国产精品,欧美在线| 午夜免费鲁丝| 亚洲熟妇熟女久久| 国产成+人综合+亚洲专区| 日日干狠狠操夜夜爽| 欧美日韩精品网址| 十八禁人妻一区二区| 12—13女人毛片做爰片一| 热99re8久久精品国产| 在线视频色国产色| 老司机福利观看| 精品国产美女av久久久久小说| 美女扒开内裤让男人捅视频| 欧美日本亚洲视频在线播放| 久久伊人香网站| 天堂动漫精品| 好看av亚洲va欧美ⅴa在| 69av精品久久久久久| 熟女电影av网| 非洲黑人性xxxx精品又粗又长| 成人三级黄色视频| 精品久久久久久久久久久久久 | 色尼玛亚洲综合影院| 深夜精品福利| 色播在线永久视频| 大型av网站在线播放| 亚洲成人免费电影在线观看| 中文字幕精品亚洲无线码一区 | 99久久精品国产亚洲精品| 日韩 欧美 亚洲 中文字幕| 91大片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 少妇熟女aⅴ在线视频| 18禁裸乳无遮挡免费网站照片 | 十八禁人妻一区二区| 一进一出抽搐gif免费好疼| 亚洲免费av在线视频| 男女下面进入的视频免费午夜 | 日本 av在线| 黄色a级毛片大全视频| 日韩视频一区二区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 好男人在线观看高清免费视频 | 正在播放国产对白刺激| 久久久国产精品麻豆| 亚洲国产精品合色在线| 亚洲色图 男人天堂 中文字幕| 国产av在哪里看| 国产激情久久老熟女| 亚洲国产欧洲综合997久久, | 日韩大码丰满熟妇| 日本免费一区二区三区高清不卡| 精品久久蜜臀av无| 日韩高清综合在线| 免费在线观看影片大全网站| 可以免费在线观看a视频的电影网站| 国产欧美日韩一区二区三| 亚洲 欧美一区二区三区| 日本免费a在线| 一边摸一边做爽爽视频免费| 精品一区二区三区四区五区乱码| 国产又爽黄色视频| 国产成人av激情在线播放| 久久中文字幕人妻熟女| 一边摸一边做爽爽视频免费| 日韩一卡2卡3卡4卡2021年| 久热这里只有精品99| 窝窝影院91人妻| а√天堂www在线а√下载| 在线十欧美十亚洲十日本专区| 午夜亚洲福利在线播放| 欧美日韩亚洲国产一区二区在线观看| 人妻丰满熟妇av一区二区三区| 日本 av在线| 亚洲第一欧美日韩一区二区三区| 国产av在哪里看| 国产区一区二久久| 在线视频色国产色| 99久久综合精品五月天人人| www日本黄色视频网| 成人亚洲精品av一区二区| 在线观看一区二区三区| 亚洲国产中文字幕在线视频| 国产蜜桃级精品一区二区三区| 伦理电影免费视频| 宅男免费午夜| 国产成人av激情在线播放| 亚洲熟妇中文字幕五十中出| 脱女人内裤的视频| 国产三级黄色录像| 国产亚洲精品综合一区在线观看 | 在线天堂中文资源库| 黄色毛片三级朝国网站| 亚洲中文av在线| 看黄色毛片网站| 亚洲性夜色夜夜综合| 亚洲自偷自拍图片 自拍| 欧美zozozo另类| 亚洲五月色婷婷综合| 真人一进一出gif抽搐免费| 久久婷婷成人综合色麻豆| 久久久久九九精品影院| 午夜免费成人在线视频| 国产又黄又爽又无遮挡在线| 丝袜美腿诱惑在线| 久久久久久免费高清国产稀缺| 夜夜夜夜夜久久久久| 久久中文看片网| 久久精品亚洲精品国产色婷小说| 欧美乱妇无乱码| 亚洲自拍偷在线| 亚洲在线自拍视频| 最近最新中文字幕大全免费视频| 久久婷婷成人综合色麻豆| 老司机福利观看| 国产三级黄色录像| 久久国产精品影院| 亚洲免费av在线视频| 美女高潮喷水抽搐中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 最近在线观看免费完整版| 国产不卡一卡二| 成人欧美大片| 中国美女看黄片| 国产野战对白在线观看| 久久热在线av| 狠狠狠狠99中文字幕| 欧美一区二区精品小视频在线| 国产亚洲欧美精品永久| 视频区欧美日本亚洲| 麻豆成人av在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧洲综合997久久, | 1024手机看黄色片| 午夜a级毛片| 亚洲av中文字字幕乱码综合 | 久久欧美精品欧美久久欧美| 日本三级黄在线观看| 在线天堂中文资源库| 亚洲一区二区三区不卡视频| 搡老妇女老女人老熟妇| 色老头精品视频在线观看| 国产成年人精品一区二区| 久99久视频精品免费| 国产野战对白在线观看| www.熟女人妻精品国产| 国产午夜精品久久久久久| 人人妻人人看人人澡| 91av网站免费观看| 午夜激情av网站| 青草久久国产| 嫩草影视91久久| 人妻久久中文字幕网| 亚洲av美国av| 黑人操中国人逼视频| 亚洲一码二码三码区别大吗| 在线观看免费视频日本深夜| 中文字幕人妻熟女乱码| av天堂在线播放| 欧美一级毛片孕妇| 欧美午夜高清在线| 人人澡人人妻人| 亚洲精品美女久久av网站| 看黄色毛片网站| 精品免费久久久久久久清纯| 午夜激情av网站| a在线观看视频网站| 老熟妇仑乱视频hdxx| 亚洲一区高清亚洲精品| 久久99热这里只有精品18| 一个人观看的视频www高清免费观看 | 午夜成年电影在线免费观看| 欧美黑人欧美精品刺激| 国产激情偷乱视频一区二区| 听说在线观看完整版免费高清| 天天一区二区日本电影三级| 日韩高清综合在线| 亚洲全国av大片| 国产精品亚洲av一区麻豆| 午夜福利一区二区在线看|