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

    基于逆繞射拋物方程法的輻射源定位技術(shù)研究

    2011-08-08 03:11:00李德鑫楊日杰王鴻吉蔣志忠
    電波科學(xué)學(xué)報(bào) 2011年4期
    關(guān)鍵詞:電波拋物輻射源

    李德鑫 楊日杰 王鴻吉 蔣志忠

    (1.海軍航空工程學(xué)院電子信息工程系,山東 煙臺(tái)264001;2.空軍哈爾濱飛行仿真技術(shù)研究所,黑龍江 哈爾濱150001)

    1.引 言

    輻射源無(wú)源定位技術(shù)是電子戰(zhàn)的一個(gè)重要組成部分,無(wú)論是傳統(tǒng)的戰(zhàn)略戰(zhàn)術(shù)電子偵察,還是電子干擾,都離不開(kāi)利用目標(biāo)輻射源的信號(hào)或外輻射源的信號(hào)對(duì)目標(biāo)進(jìn)行定位。目前的無(wú)源定位系統(tǒng)采用的基本技術(shù)有到達(dá)角(AOA,Angle of Arrival)定位技術(shù)、到達(dá)時(shí)差(TDOA,Time Difference of Arrival)定位技術(shù)、到達(dá)頻差(FDOA,F(xiàn)requency Difference of Arrival)定位技術(shù),以及聯(lián)合定位技術(shù)等[1-2]。受環(huán)境因素影響,傳統(tǒng)的定位算法對(duì)于復(fù)雜環(huán)境非合作輻射源信號(hào)的來(lái)波方向、傳播時(shí)間、距離等參數(shù)的估計(jì)將會(huì)變得十分復(fù)雜。

    拋物方程法是一種前向全波法,能夠處理垂直與水平方向上介質(zhì)非均勻分布的環(huán)境,可方便地用于分析地形起伏變化、地表電磁特性變化及傳播空間媒質(zhì)折射率變化等條件下的電波傳播問(wèn)題[3-4]。澳大利亞Spencer等人在此基礎(chǔ)上提出了逆繞射拋物方程定位系統(tǒng)(IDPELS,Inverse Diffraction Parabolic Wave Equation Localization System),能夠?qū)崿F(xiàn)復(fù)雜環(huán)境下的輻射源定位,在實(shí)測(cè)中具有很好的收斂效果[5-6]。

    傳統(tǒng)的IDPELS沒(méi)有考慮遠(yuǎn)距離、大范圍的輻射源定位、誤差橢圓、電磁剖面取值方法等問(wèn)題,限制了該方法的實(shí)際應(yīng)用。本文將在理想大氣輻射源定位的基礎(chǔ)上,對(duì)一個(gè)和兩個(gè)輻射源在標(biāo)準(zhǔn)大氣、地形遮蔽條件下的定位問(wèn)題進(jìn)行分析,對(duì)定位誤差、數(shù)據(jù)點(diǎn)取值、誤差橢圓的繪制等問(wèn)題進(jìn)行研究。

    2.逆繞射拋物方程法

    設(shè)電磁場(chǎng)時(shí)諧因子為e-ikx,標(biāo)量ψ表示電磁場(chǎng)分量,k為真空中傳播常數(shù),n為折射指數(shù),直角坐標(biāo)系(x,y,z)下,假設(shè)ψ與y無(wú)關(guān),即大氣沿水平方向均勻分布,只隨距離(x軸)和高度(z軸)變化,則水平極化波中只有Ey為非零電場(chǎng)分量(ψ(x,z)=Ey(x,z)),垂直極化波中只有Hy為非零磁場(chǎng)分量(ψ(x,z)=Hy(x,z))。電波傳播過(guò)程中,ψ(x,z)滿(mǎn)足以下二維標(biāo)量波動(dòng)方程。

    定義沿x軸正向傳播的衰減函數(shù)為

    將式(2)代入式(1),并作適當(dāng)近似即可得到高度-距離空間的二維標(biāo)量拋物方程,采用Feit-Fleck近似法[7],得到 Feit-Fleck型拋物方程

    采用Leontovich邊界條件[8],視地海面為光滑理想邊界。將 Cosine-taper(Tukey)窗函數(shù)[9]作為計(jì)算域的上邊界,利用分步傅立葉變換(SSFT,Split-step Fourier transform)算法求解得

    式中:ae為地球半徑。定義反映路徑上障礙物的繞射效應(yīng),N= eik(n-2+z/ae)Δx反映傳播媒介的折射效應(yīng),對(duì)于不規(guī)則地形條件N=eik(n-2+z/a+z·?2T/?x2)Δx. e

    拋物方程法是已知發(fā)射源初始場(chǎng)求解電波的空間分布;逆繞射拋物方程法與之相反,是由遠(yuǎn)場(chǎng)某一觀測(cè)區(qū)域的場(chǎng)分布出發(fā),求解電波的逆繞射問(wèn)題,由于輻射源初始場(chǎng)是收斂的,從而可估計(jì)發(fā)射源的位置信息。電波在媒質(zhì)中傳播時(shí),部分能量向空中輻射或被媒質(zhì)和傳播邊界吸收導(dǎo)致?lián)p耗,且拋物方程是對(duì)橢圓型波動(dòng)方程忽略后向傳播的前向近似,因此無(wú)法對(duì)電波傳播完全求逆。但是當(dāng)邊界條件已知時(shí),該算法具有很好的收斂性。

    對(duì)式(4)做逆函數(shù)推導(dǎo)可得電波傳播的逆繞射拋物方程公式為

    式(5)中的初始場(chǎng)u(x0,z)與式(4)中的不同,為遠(yuǎn)場(chǎng)截獲的某一區(qū)域的場(chǎng)分布。

    用逆繞射拋物方程模型計(jì)算水平極化條件下,電磁波在光滑海面上視距傳播時(shí)的場(chǎng)分布特性。Leontovich邊界條件中,海水介電常數(shù)εr=80,海水導(dǎo)電率σ=4 S/m;其它仿真條件,輻射源1高度z0=300 m,電波頻率f=1 GHz,高斯方向圖,天線(xiàn)增益35 dB,天線(xiàn)3 dB帶寬為3.5°,觀測(cè)距離為30 km.輻射源2高度z0=100 m,電波頻率f=1.8 GHz,高斯方向圖,天線(xiàn)增益29.5 dB,天線(xiàn)3 dB帶寬為3°,觀測(cè)距離為20 km.海面上為理想均勻大氣n=1。

    圖1(a)為輻射源1的逆繞射場(chǎng)強(qiáng)分布圖,圖1(b)為場(chǎng)強(qiáng)等高線(xiàn)圖。圖中u為折射脊線(xiàn)(主要由折射場(chǎng)分量疊加而成),v為反射脊線(xiàn)(主要由海面反射場(chǎng)分量疊加而成),A(x=30000,z=300.79)為場(chǎng)強(qiáng)極大值點(diǎn)亦為兩條測(cè)向線(xiàn)的收斂點(diǎn),即算法的估值點(diǎn)。由于算法步長(zhǎng)(Δz≈1.6 m)的影響,對(duì)于高度z的估算產(chǎn)生了誤差。由圖1可以看出理想條件下算法具有很好的收斂性和估值精度。

    圖2為兩個(gè)輻射源的場(chǎng)強(qiáng)分布圖和等高線(xiàn)圖,共有兩個(gè)收斂點(diǎn)(20000,100.31),(29600,299.35);由于兩個(gè)輻射源相互干涉,對(duì)于輻射源1的距離估算產(chǎn)生了誤差。由圖2可以看出,對(duì)于兩個(gè)輻射源的定位問(wèn)題,算法仍具有很好的收斂性。

    3.數(shù)據(jù)點(diǎn)的選取

    選取觀測(cè)位置垂直剖面上的一組場(chǎng)分布數(shù)據(jù)作為觀測(cè)數(shù)據(jù),在實(shí)際情況下,大范圍數(shù)據(jù)點(diǎn)的獲取需要大型天線(xiàn)和平臺(tái)的大范圍機(jī)動(dòng)(如車(chē)載[5-6]或機(jī)載),并獲取重點(diǎn)區(qū)域的電磁數(shù)據(jù)。數(shù)據(jù)點(diǎn)的取值范圍及選取的位置會(huì)對(duì)定位造成誤差,以標(biāo)準(zhǔn)大氣條件下對(duì)輻射源1的定位為例,數(shù)據(jù)點(diǎn)的選取區(qū)域如圖3所示,橢圓圈定的垂直剖面為觀測(cè)數(shù)據(jù),仿真實(shí)驗(yàn)數(shù)據(jù)點(diǎn)數(shù)由上下兩邊界向中點(diǎn)遞減(圖中標(biāo)出了1024和512個(gè)數(shù)據(jù)點(diǎn)的取值區(qū)域),對(duì)于不同數(shù)據(jù)點(diǎn)情況下算法的估計(jì)值如表1所示。

    數(shù)據(jù)點(diǎn)數(shù)為1024時(shí),由于大氣折射的影響,折射定位線(xiàn)給估值帶來(lái)了誤差,同時(shí)由于海面反射的干涉和疊加,使得對(duì)于高度z的估計(jì)有較大的誤差。數(shù)據(jù)點(diǎn)數(shù)為512時(shí),由于截取數(shù)據(jù)點(diǎn)距海面較遠(yuǎn),受海面反射效應(yīng)的影響減弱,對(duì)于z估計(jì)的精度提高;但由于精度較高的電波直線(xiàn)傳播數(shù)據(jù)點(diǎn)數(shù)的減少,對(duì)于距離x的估計(jì)誤差增加。中部數(shù)據(jù)點(diǎn)所含位置信息相類(lèi)似,減少數(shù)據(jù)點(diǎn)數(shù)目對(duì)于估值并沒(méi)有太大的影響,體現(xiàn)出算法具有很好的穩(wěn)定性。算法的誤差在32個(gè)數(shù)據(jù)點(diǎn)時(shí)開(kāi)始急劇增加,小于16點(diǎn)時(shí),算法失效。數(shù)據(jù)點(diǎn)的選取與輻射源的高度和觀測(cè)距離有關(guān),隨著觀測(cè)距離的增加受海面反射影響逐步減弱,但同時(shí)數(shù)據(jù)點(diǎn)所含位置信息也在減少。

    表1 單個(gè)輻射源不同數(shù)據(jù)點(diǎn)數(shù)的估計(jì)值

    兩個(gè)輻射源定位的數(shù)據(jù)點(diǎn)選取情況與單個(gè)輻射源定位時(shí)相類(lèi)似,不同之處在于,兩個(gè)輻射源之間的干涉和疊加會(huì)對(duì)彼此的定位造成誤差。輻射場(chǎng)范圍較小的輻射源由于完全包含在大范圍輻射源的輻射場(chǎng)中,受干涉影響其誤差相對(duì)較大。對(duì)于輻射場(chǎng)范圍較大的輻射源,選取干涉臨界點(diǎn)以上的數(shù)據(jù)點(diǎn)會(huì)提高其定位精度,但同時(shí)會(huì)降低輻射場(chǎng)范圍較小的輻射源定位精度。兩個(gè)輻射源數(shù)據(jù)點(diǎn)選取與輻射源的高度、電波傳播仰角、觀測(cè)距離有關(guān)。需要指出的是,算法僅能解決不相干輻射源的多點(diǎn)定位問(wèn)題,對(duì)于相干輻射源,由于輻射場(chǎng)之間相互干擾,導(dǎo)致算法失效。

    對(duì)于大氣折射的特殊情況——大氣波導(dǎo),由于其形成具有區(qū)域性、不均勻和不確定性,對(duì)于非合作輻射源的定位,很難估算大范圍波導(dǎo)的真實(shí)折射和陷獲情況。常見(jiàn)的大氣波導(dǎo)陷獲角一般很小,電波異常傳播的區(qū)域一般局限在貼海面較低區(qū)域,為獲取較高精度的估值,應(yīng)盡量選取遠(yuǎn)離波導(dǎo)區(qū)域的數(shù)據(jù)點(diǎn)。

    因此,關(guān)于數(shù)據(jù)點(diǎn)的選取,應(yīng)選擇盡量遠(yuǎn)離海面反射的較高區(qū)域,同時(shí)盡可能多的獲得位置信息含量高的數(shù)據(jù)點(diǎn),以提高定位精度。

    4.誤差橢圓

    對(duì)于定位誤差,可利用誤差橢圓來(lái)說(shuō)明定位區(qū)域和收斂范圍。圖4中,定義u為折射脊線(xiàn),v為反射脊線(xiàn),A為輻射源真值點(diǎn),O為估算的收斂點(diǎn),θ為u、v的交角,γ為v、x的交角,E1、E2為測(cè)向線(xiàn)誤差。

    圖4 定位誤差分析

    式(6)為uOv坐標(biāo)系下的橢圓,與xOy坐標(biāo)系的變換公式為

    其中:

    將式(7)代入式(6)即可獲得直角坐標(biāo)系下的誤差橢圓公式。

    對(duì)于地球表面大氣條件下的輻射源定位問(wèn)題,按照國(guó)際電訊聯(lián)盟(ITU-R)[10]規(guī)定的對(duì)流層標(biāo)準(zhǔn)大氣條件?n/?h≈-40×10-6計(jì)算折射指數(shù)n,對(duì)單個(gè)輻射源和兩個(gè)輻射源的定位情況進(jìn)行仿真,其余仿真條件與理想大氣情況相同。繪制定位區(qū)域場(chǎng)強(qiáng)等高線(xiàn)圖,圖5為表1中512個(gè)數(shù)據(jù)點(diǎn)情況;圖6為兩個(gè)輻射源1024個(gè)數(shù)據(jù)點(diǎn)情況;圖7為兩個(gè)輻射源選取靠近上邊界的512個(gè)數(shù)據(jù)點(diǎn)情況(注意與圖5取值區(qū)域不同)。

    圖5~6中“*”為算法的估值點(diǎn),“Δ”為輻射源真實(shí)位置,由兩圖可以看出,誤差橢圓隨著估計(jì)誤差的增大面積隨之增大,能夠反映出收斂域的面積、方向等信息,很好地說(shuō)明了估值點(diǎn)與真值點(diǎn)的關(guān)系,直觀顯示出算法的收斂區(qū)域和定位精度。

    需要說(shuō)明的是,圖5~7繪制的是場(chǎng)強(qiáng)等高線(xiàn)圖,由于折射、反射疊加效應(yīng)的影響,折射脊線(xiàn)與反射脊線(xiàn)不再是圖1中理想情況下的等高線(xiàn),而是一條收斂帶中各區(qū)域場(chǎng)強(qiáng)極大值疊加而成的收斂線(xiàn)(由于各極值不相等,故不再是一條等高線(xiàn)),三維圖可以很好地顯示出這種收斂的趨勢(shì)。

    如前文所述,由于兩個(gè)輻射源間的相互干涉和疊加,對(duì)各自的定位估值帶來(lái)了較大的誤差。圖5與圖6中關(guān)于輻射源1的定位估值點(diǎn)分別為(29500,296.16)和(29400,278.65),相應(yīng)的誤差橢圓面積與方向有所不同,誤差橢圓面積隨誤差增加而增大,橢圓的傾角說(shuō)明定位值的誤差朝向;圖6輻射源2的定位值為(19900,85.98),該輻射源的輻射區(qū)域完全被輻射源1覆蓋,由反射折射疊加帶來(lái)的誤差較大,驗(yàn)證了前文的論述。

    圖7為兩個(gè)輻射源512個(gè)數(shù)據(jù)點(diǎn)時(shí)的定位情況,其余仿真條件與圖6相同,此處的數(shù)據(jù)點(diǎn)由于選取上部的512個(gè)點(diǎn),剔除了受海面反射影響較為強(qiáng)烈的下部數(shù)據(jù)點(diǎn)。因此,定位線(xiàn)收斂性強(qiáng),定位效果較好。兩個(gè)輻射源的定位值分別為(29400,283.43),(19600,101.91),與圖6的定位值相比,高度估值精度有較大提高;但是由于剔除掉了部分?jǐn)y有反射距離信息的數(shù)據(jù)點(diǎn),受折射效應(yīng)的影響,距離誤差有所增加。兩個(gè)輻射源相互間的干涉影響主要體現(xiàn)在海面反射效應(yīng)上,上述仿真實(shí)驗(yàn)中,輻射源2高度較低,觀測(cè)數(shù)據(jù)中由于海面反射效應(yīng)的混疊導(dǎo)致誤差增加,因而定位誤差增大。

    5.地形遮蔽條件下的輻射源定位

    對(duì)流層電波傳播中,不規(guī)則地形是影響電波傳播特性的另一主要因素,可利用類(lèi)似正弦包絡(luò)形狀的不規(guī)則地形分析地形遮蔽對(duì)定位算法的影響。地形模型為

    式中:h為包絡(luò)的最大高度;x1為包絡(luò)半高程所在的距離;w為包絡(luò)的半寬度。仿真設(shè)定值,地形1為h=600 m,x1=20000 m,w=1000 m;地形2為h=200 m,x1=10000 m,w=1000 m(此處設(shè)定值是相對(duì)于觀測(cè)區(qū)域而定)。

    標(biāo)準(zhǔn)大氣條件下,受地形遮蔽影響的兩個(gè)輻射源輻射場(chǎng)分布如圖8所示,此處為地形1,其它仿真條件與圖6相同。利用逆繞射拋物方程法的場(chǎng)分布圖如圖9所示。

    依前例,“*”為算法的估值點(diǎn),“Δ”為輻射源真實(shí)位置;由圖9可以看出,在地形遮蔽條件下,算法仍有很好的收斂性。圖10中,受地形遮蔽影響,下部分觀測(cè)值的信息損失很大,對(duì)于輻射源1的定位值為(29400,277.06),與圖6的精度相當(dāng),輻射源2的位置信息損失過(guò)大,無(wú)法定位。

    輻射源定位精度與地形最大高度有關(guān),一般情況下當(dāng)山峰的高度高于輻射源天線(xiàn)3 dB帶寬角所能達(dá)到的最大視距高度時(shí),定位精度將會(huì)急劇變差,直至算法失效。圖9中,山峰的高度高于3 d B帶寬角所能達(dá)到的最大視距高度,低于電波傳播的最大仰角視距高度,定位失效。

    圖10為圖9的等高線(xiàn)場(chǎng)分布圖,兩圖山峰位置處均有一個(gè)明顯的劈尖,利用該劈尖信息可以估算地形的最大高度。因此,當(dāng)輻射源位置已知情況下,該算法還可應(yīng)用于地形的估算。

    圖11為地形2條件下的局部等高線(xiàn)場(chǎng)分布圖,兩個(gè)輻射源的定位值分別為(29600,277.06)和(19900,82.80),由定位誤差橢圓可以看出,算法在較低山峰遮蓋下的定位精度與標(biāo)準(zhǔn)大氣條件下相當(dāng),體現(xiàn)了算法的穩(wěn)定性和實(shí)用性。

    6.結(jié) 論

    在理想大氣輻射源定位的基礎(chǔ)上,對(duì)一個(gè)和兩個(gè)輻射源在標(biāo)準(zhǔn)大氣、地形遮蔽條件下的定位問(wèn)題進(jìn)行了分析,研究了定位誤差、數(shù)據(jù)點(diǎn)取值、誤差橢圓繪制等問(wèn)題,提高了算法的定位精度,擴(kuò)大了算法的應(yīng)用范圍。

    逆繞射拋物方程法可對(duì)遠(yuǎn)距離、大范圍復(fù)雜條件下的電波傳播進(jìn)行預(yù)測(cè),對(duì)非合作輻射源進(jìn)行定位,對(duì)未知地形進(jìn)行估算,具有很好的應(yīng)用前景。

    逆繞射拋物方程法有著自身的缺陷,需要觀測(cè)大范圍的場(chǎng)強(qiáng)數(shù)據(jù)來(lái)保證估值的精度,進(jìn)而需要大型天線(xiàn)和平臺(tái)的大范圍機(jī)動(dòng)來(lái)保障。算法的收斂性與計(jì)算步長(zhǎng)有關(guān),步長(zhǎng)越小估值精度越高,但是算法的復(fù)雜度隨之增加,計(jì)算時(shí)間增大,實(shí)時(shí)性變差。

    后續(xù)的工作應(yīng)當(dāng)包括數(shù)據(jù)點(diǎn)選取建模、誤差概率統(tǒng)計(jì)及動(dòng)態(tài)海面、雜波、復(fù)雜地形條件下的建模,以期在最佳高度、距離觀測(cè)最少的數(shù)據(jù)點(diǎn)獲得較高的定位精度,提高算法的適用性和穩(wěn)定性。

    [1] 曹軼超,方建安,羅賢云.一個(gè)基于TDOA的無(wú)線(xiàn)定位新算法[J] .電波科學(xué)學(xué)報(bào),2008,23(5):841-846.CAO Yichao,F(xiàn)ANG Jianan,LUO Xianyun.Constrained least-squares and wireless location based on TDOA measurements[J] .Chinese Journal of Radio Science,2008,23(5):841-846.(in Chinese)

    [2] 王 強(qiáng),鐘丹星,鄧新蒲,等.基于圓錐向角和相位差變化率的運(yùn)動(dòng)單平臺(tái)無(wú)源定位算法和精度分析[J] .航空學(xué)報(bào),2009,30(8):1497-1502.WANG Qiang,ZHONG Danxing,DENG Xinpu,et al.Single moving observer passive localization algorithm using conical bearings and rate of phase difference change and its accuracy analysis[J] .Acta Aeronautica et Astronautica Sinica,2009,30(8):1497-1502.(in Chinese)

    [3] 胡繪斌,毛鈞杰,柴舜連.電波傳播中求解寬角拋物方程的誤差分析[J] .電波科學(xué)學(xué)報(bào),2006,21(2):199-203.HU Huibin,MAO Junjie,CHAI Shunlian.Error analysis of wide-angle parabolic equation for radio propagation[J] .Chinese Journal of Radio Science,2006,21(2):199-203.(in Chinese)

    [4] 姚景順,楊世興.拋物方程模型在海上電波傳播中的應(yīng)用[J] .電波科學(xué)學(xué)報(bào),2009,24(3):493-497.YAO Jingshun,YANG Shixing.A terrain parabolic equation model for propagation over the ocean[J] .Chinese Journal of Radio Science,2009,24(3):493-497.(in Chinese)

    [5] SPENCER T A,WALKER R A,HAWKES R M.Inverse diffraction parabolic wave equation localization system (IDPELS)[C] //International Symposium on GNSS/GPS.Sydney,Australia,2004:1-15.

    [6] SPENCER T A,WALKER R A,HAWKES R M.Inverse diffraction parabolic wave equation localization system (IDPELS)[J] .Journal of Global Positioning Systems,2005,4(1-2):245-257.

    [7] FEIT M D,F(xiàn)LECK J A.Light propagation in graded index fibers[J] .Appl Opt,1978,17(24):3990-3998.

    [8] KUTTER R,DOCKEY G.Theoretical description of parabolic approximation Fourier split-step method of representing electromagnetic propagation in the troposphere[J] .Radio Science,1991,26(2):381-393.

    [9] HARRIS F J.On the use of windows for harmonic analysis with the discrete Fourier transform[J] .IEEE Proc,1978,66(1):51-83.

    [10] 熊 皓.無(wú)線(xiàn)電波傳播[M] .北京:電子工業(yè)出版社,2000:361-375.

    猜你喜歡
    電波拋物輻射源
    高空拋物罪的實(shí)踐擴(kuò)張與目的限縮
    法律方法(2022年2期)2022-10-20 06:45:28
    永遠(yuǎn)的紅色電波
    The Speed of Light
    基于博弈論的GRA-TOPSIS輻射源威脅評(píng)估方法
    關(guān)于拋物-拋物Keller-Segel類(lèi)模型的全局解和漸近性
    瞌睡電波
    幽默大師(2019年11期)2019-11-23 08:47:58
    數(shù)字電視外輻射源雷達(dá)多旋翼無(wú)人機(jī)微多普勒效應(yīng)實(shí)驗(yàn)研究
    不要高空拋物!
    高空莫拋物
    外輻射源雷達(dá)直升機(jī)旋翼參數(shù)估計(jì)方法
    99香蕉大伊视频| 午夜福利免费观看在线| bbb黄色大片| 精品国产乱子伦一区二区三区 | 免费人妻精品一区二区三区视频| 欧美少妇被猛烈插入视频| 一本久久精品| av在线老鸭窝| 欧美在线一区亚洲| 亚洲avbb在线观看| 久久久久国产精品人妻一区二区| 宅男免费午夜| 制服诱惑二区| 国产免费视频播放在线视频| 亚洲av成人不卡在线观看播放网 | 又黄又粗又硬又大视频| 国产日韩欧美在线精品| 国产伦人伦偷精品视频| 精品欧美一区二区三区在线| 色老头精品视频在线观看| 亚洲欧美色中文字幕在线| 国产精品av久久久久免费| 热99re8久久精品国产| 一区二区三区四区激情视频| 丝袜在线中文字幕| 91成年电影在线观看| 99热国产这里只有精品6| 脱女人内裤的视频| 亚洲av成人不卡在线观看播放网 | 国产老妇伦熟女老妇高清| 天天躁狠狠躁夜夜躁狠狠躁| 久久人妻熟女aⅴ| 国产成人av教育| 中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放| 久久国产精品大桥未久av| 极品人妻少妇av视频| 午夜福利视频精品| 少妇精品久久久久久久| avwww免费| 久久精品国产亚洲av高清一级| 18禁裸乳无遮挡动漫免费视频| 人妻人人澡人人爽人人| 美女中出高潮动态图| 成年美女黄网站色视频大全免费| 自拍欧美九色日韩亚洲蝌蚪91| 不卡av一区二区三区| 欧美大码av| 又紧又爽又黄一区二区| 捣出白浆h1v1| 日本撒尿小便嘘嘘汇集6| 久久国产亚洲av麻豆专区| 天天躁夜夜躁狠狠躁躁| 日韩欧美免费精品| 午夜久久久在线观看| 精品福利观看| 久久久欧美国产精品| 一级,二级,三级黄色视频| 国产成人啪精品午夜网站| a级毛片黄视频| 美女午夜性视频免费| 国产男女超爽视频在线观看| 一本色道久久久久久精品综合| 久久精品国产综合久久久| 国产日韩欧美亚洲二区| 美女大奶头黄色视频| svipshipincom国产片| 国产精品免费视频内射| 精品少妇黑人巨大在线播放| 久久这里只有精品19| 丰满少妇做爰视频| 国产主播在线观看一区二区| 亚洲色图 男人天堂 中文字幕| 亚洲欧美日韩高清在线视频 | 国产欧美亚洲国产| 成人黄色视频免费在线看| 久久99热这里只频精品6学生| 每晚都被弄得嗷嗷叫到高潮| 中文字幕人妻丝袜制服| 一本色道久久久久久精品综合| 免费日韩欧美在线观看| 在线永久观看黄色视频| 99热国产这里只有精品6| 亚洲av电影在线观看一区二区三区| 热re99久久国产66热| 成年av动漫网址| 巨乳人妻的诱惑在线观看| 亚洲七黄色美女视频| 午夜久久久在线观看| 少妇 在线观看| 国产伦理片在线播放av一区| 亚洲一卡2卡3卡4卡5卡精品中文| 丰满少妇做爰视频| 80岁老熟妇乱子伦牲交| 99香蕉大伊视频| 热re99久久国产66热| 精品亚洲成a人片在线观看| 免费观看av网站的网址| 91大片在线观看| 国产成人欧美| 中文字幕另类日韩欧美亚洲嫩草| 精品久久蜜臀av无| 老司机亚洲免费影院| 日韩人妻精品一区2区三区| 欧美在线一区亚洲| 亚洲国产精品一区二区三区在线| 日本a在线网址| 男女边摸边吃奶| 国产一级毛片在线| 夜夜夜夜夜久久久久| 99久久综合免费| 欧美日韩亚洲高清精品| 80岁老熟妇乱子伦牲交| 精品免费久久久久久久清纯 | 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 99久久国产精品久久久| www.熟女人妻精品国产| 欧美乱码精品一区二区三区| 国产精品熟女久久久久浪| 涩涩av久久男人的天堂| 黄频高清免费视频| 女人精品久久久久毛片| 欧美日韩国产mv在线观看视频| 99国产极品粉嫩在线观看| 久久久久精品人妻al黑| 欧美成狂野欧美在线观看| 国产成人精品无人区| 老司机午夜福利在线观看视频 | 人人澡人人妻人| 手机成人av网站| 午夜免费观看性视频| 99精国产麻豆久久婷婷| 精品国产乱码久久久久久小说| 国产av国产精品国产| 丝瓜视频免费看黄片| 亚洲久久久国产精品| 这个男人来自地球电影免费观看| 老司机午夜福利在线观看视频 | 高清av免费在线| 欧美日韩av久久| 91成人精品电影| 99国产综合亚洲精品| 午夜老司机福利片| 窝窝影院91人妻| 午夜视频精品福利| 欧美黑人欧美精品刺激| 免费人妻精品一区二区三区视频| 亚洲精品国产一区二区精华液| 老司机在亚洲福利影院| 热99久久久久精品小说推荐| 交换朋友夫妻互换小说| 一进一出抽搐动态| 汤姆久久久久久久影院中文字幕| 欧美激情高清一区二区三区| 中文字幕人妻丝袜制服| 久久av网站| 国产免费av片在线观看野外av| 欧美激情 高清一区二区三区| 久久久国产欧美日韩av| 天天影视国产精品| 飞空精品影院首页| 久久人人爽人人片av| 欧美日韩黄片免| 性色av一级| 又紧又爽又黄一区二区| 精品亚洲成a人片在线观看| 少妇被粗大的猛进出69影院| 国产伦人伦偷精品视频| 日韩视频一区二区在线观看| 人成视频在线观看免费观看| 在线看a的网站| 18禁黄网站禁片午夜丰满| 多毛熟女@视频| av天堂久久9| 99九九在线精品视频| www.自偷自拍.com| 国产欧美日韩精品亚洲av| 女性生殖器流出的白浆| 一进一出抽搐动态| 欧美日韩亚洲国产一区二区在线观看 | av国产精品久久久久影院| 热99国产精品久久久久久7| 日本五十路高清| 亚洲精品中文字幕在线视频| 国产老妇伦熟女老妇高清| 亚洲精品国产av蜜桃| 亚洲久久久国产精品| 亚洲五月婷婷丁香| av片东京热男人的天堂| 91大片在线观看| 老司机亚洲免费影院| 女人高潮潮喷娇喘18禁视频| 亚洲精品粉嫩美女一区| 国产黄色免费在线视频| 国产激情久久老熟女| 一本大道久久a久久精品| 嫩草影视91久久| 18禁国产床啪视频网站| 一级片免费观看大全| 免费不卡黄色视频| 91老司机精品| 国产精品99久久99久久久不卡| 啦啦啦免费观看视频1| 精品免费久久久久久久清纯 | 欧美另类亚洲清纯唯美| 久久精品国产a三级三级三级| 日韩人妻精品一区2区三区| 丁香六月欧美| 国产精品香港三级国产av潘金莲| 国产成人一区二区三区免费视频网站| 黄网站色视频无遮挡免费观看| 俄罗斯特黄特色一大片| 一级毛片精品| 欧美黄色片欧美黄色片| www.精华液| 日韩制服丝袜自拍偷拍| 男人操女人黄网站| 黑人巨大精品欧美一区二区mp4| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 99国产精品一区二区三区| 久久久久精品人妻al黑| 正在播放国产对白刺激| 日韩电影二区| 各种免费的搞黄视频| 激情视频va一区二区三区| 黄频高清免费视频| 99九九在线精品视频| 久久av网站| 久久久久国产一级毛片高清牌| 飞空精品影院首页| 亚洲精品自拍成人| 女性被躁到高潮视频| av有码第一页| 精品一区在线观看国产| 秋霞在线观看毛片| 在线观看人妻少妇| 曰老女人黄片| 在线观看免费视频网站a站| 色精品久久人妻99蜜桃| 国产精品秋霞免费鲁丝片| 最黄视频免费看| av电影中文网址| 啦啦啦中文免费视频观看日本| 日韩免费高清中文字幕av| 精品亚洲乱码少妇综合久久| 美女脱内裤让男人舔精品视频| 最新在线观看一区二区三区| 不卡av一区二区三区| 午夜福利视频在线观看免费| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 999久久久精品免费观看国产| 满18在线观看网站| 一本久久精品| 国产精品自产拍在线观看55亚洲 | 两性夫妻黄色片| 在线 av 中文字幕| 肉色欧美久久久久久久蜜桃| 国产成人a∨麻豆精品| 欧美老熟妇乱子伦牲交| 久久久精品94久久精品| 精品一品国产午夜福利视频| 欧美黄色淫秽网站| av线在线观看网站| 成年人黄色毛片网站| 国产精品久久久久久精品电影小说| 国产三级黄色录像| 久久影院123| 免费少妇av软件| 99热国产这里只有精品6| 秋霞在线观看毛片| 大型av网站在线播放| 亚洲av成人一区二区三| 日本精品一区二区三区蜜桃| 久久国产精品人妻蜜桃| 国产不卡av网站在线观看| 久久久久视频综合| 国产亚洲精品第一综合不卡| 日本撒尿小便嘘嘘汇集6| 免费在线观看影片大全网站| 制服诱惑二区| 啦啦啦啦在线视频资源| 亚洲人成电影观看| 欧美精品亚洲一区二区| 欧美在线一区亚洲| 色老头精品视频在线观看| 国产精品九九99| 亚洲免费av在线视频| 国产黄色免费在线视频| 一区在线观看完整版| 啦啦啦在线免费观看视频4| 超碰97精品在线观看| 黄色视频,在线免费观看| 男女高潮啪啪啪动态图| 激情视频va一区二区三区| 国产免费现黄频在线看| 中文精品一卡2卡3卡4更新| 啦啦啦啦在线视频资源| 五月开心婷婷网| 欧美日韩福利视频一区二区| 美国免费a级毛片| 国产又爽黄色视频| 国产成人精品无人区| 国产精品一区二区在线不卡| 欧美国产精品va在线观看不卡| 十八禁高潮呻吟视频| 亚洲国产日韩一区二区| 操美女的视频在线观看| 不卡av一区二区三区| 在线永久观看黄色视频| 电影成人av| av国产精品久久久久影院| 老熟女久久久| 中文字幕人妻丝袜制服| 亚洲avbb在线观看| 久久久精品区二区三区| 久久久水蜜桃国产精品网| 亚洲伊人久久精品综合| 91九色精品人成在线观看| 美女中出高潮动态图| 亚洲精华国产精华精| 久久 成人 亚洲| 欧美在线黄色| 久久精品久久久久久噜噜老黄| 欧美精品啪啪一区二区三区 | 欧美一级毛片孕妇| 高潮久久久久久久久久久不卡| 激情视频va一区二区三区| 亚洲欧美精品自产自拍| av国产精品久久久久影院| av欧美777| 欧美人与性动交α欧美精品济南到| 在线观看免费高清a一片| 免费黄频网站在线观看国产| 色精品久久人妻99蜜桃| 亚洲色图综合在线观看| 一边摸一边做爽爽视频免费| 色94色欧美一区二区| 日本猛色少妇xxxxx猛交久久| 久久久久久免费高清国产稀缺| 国产成人精品久久二区二区91| 高清欧美精品videossex| 叶爱在线成人免费视频播放| 国产精品久久久久成人av| 国内毛片毛片毛片毛片毛片| av不卡在线播放| 中文字幕人妻丝袜制服| 超碰成人久久| 国产成人精品久久二区二区91| 国产精品熟女久久久久浪| 国产伦人伦偷精品视频| a级毛片黄视频| 亚洲欧美精品自产自拍| 美女国产高潮福利片在线看| 最近最新免费中文字幕在线| 亚洲三区欧美一区| 亚洲熟女精品中文字幕| 18禁观看日本| 操美女的视频在线观看| 天天影视国产精品| 亚洲欧洲日产国产| 国产精品久久久人人做人人爽| www.精华液| 久久精品aⅴ一区二区三区四区| 中文精品一卡2卡3卡4更新| 啦啦啦免费观看视频1| 自拍欧美九色日韩亚洲蝌蚪91| 啦啦啦免费观看视频1| 婷婷丁香在线五月| av视频免费观看在线观看| 久久久久精品国产欧美久久久 | 国产精品自产拍在线观看55亚洲 | 成人三级做爰电影| 热re99久久精品国产66热6| 99精国产麻豆久久婷婷| 久久国产精品人妻蜜桃| 日韩电影二区| 久久精品久久久久久噜噜老黄| 国产精品熟女久久久久浪| 日本黄色日本黄色录像| 亚洲精品乱久久久久久| 97人妻天天添夜夜摸| 日韩熟女老妇一区二区性免费视频| 国产精品免费大片| 99久久人妻综合| 久久免费观看电影| 国产1区2区3区精品| 热99久久久久精品小说推荐| 久久久精品国产亚洲av高清涩受| 国产精品99久久99久久久不卡| 大香蕉久久网| www.熟女人妻精品国产| 伊人久久大香线蕉亚洲五| 老鸭窝网址在线观看| 美国免费a级毛片| 99re6热这里在线精品视频| 中文字幕色久视频| 久久人人爽人人片av| 男女下面插进去视频免费观看| 高清欧美精品videossex| 男女免费视频国产| 亚洲五月色婷婷综合| 91国产中文字幕| 免费在线观看影片大全网站| 极品人妻少妇av视频| 在线观看人妻少妇| 亚洲欧洲日产国产| 99精品欧美一区二区三区四区| 精品国产乱子伦一区二区三区 | av电影中文网址| 男女之事视频高清在线观看| 水蜜桃什么品种好| tocl精华| 欧美+亚洲+日韩+国产| 桃花免费在线播放| 国产主播在线观看一区二区| 久久人人爽av亚洲精品天堂| www.av在线官网国产| 国产免费一区二区三区四区乱码| 亚洲av日韩在线播放| 丰满少妇做爰视频| 在线永久观看黄色视频| 国产伦理片在线播放av一区| 欧美黄色淫秽网站| 伊人久久大香线蕉亚洲五| 一区二区三区精品91| 精品高清国产在线一区| www.999成人在线观看| 成人18禁高潮啪啪吃奶动态图| 国产真人三级小视频在线观看| 国产欧美日韩一区二区精品| 美女福利国产在线| a级片在线免费高清观看视频| 人人妻人人澡人人看| 国产精品国产三级国产专区5o| 交换朋友夫妻互换小说| 亚洲欧美清纯卡通| 精品久久久久久久毛片微露脸 | 欧美激情久久久久久爽电影 | 欧美黑人精品巨大| 桃花免费在线播放| 久久人人97超碰香蕉20202| 亚洲色图 男人天堂 中文字幕| 精品少妇一区二区三区视频日本电影| 最近最新中文字幕大全免费视频| 精品久久久精品久久久| 日韩熟女老妇一区二区性免费视频| 夫妻午夜视频| 中文字幕人妻丝袜一区二区| 午夜91福利影院| 国产成人精品无人区| 亚洲欧美激情在线| 久久精品久久久久久噜噜老黄| 波多野结衣av一区二区av| 午夜成年电影在线免费观看| 捣出白浆h1v1| 欧美人与性动交α欧美软件| 亚洲精品久久久久久婷婷小说| 午夜视频精品福利| 精品久久蜜臀av无| 午夜福利,免费看| 中文字幕高清在线视频| 纵有疾风起免费观看全集完整版| 超碰97精品在线观看| 亚洲国产av影院在线观看| 啦啦啦免费观看视频1| 国产精品一区二区在线不卡| 人人妻人人澡人人爽人人夜夜| 久久久久网色| 日本av免费视频播放| 99国产极品粉嫩在线观看| 欧美另类亚洲清纯唯美| 精品人妻熟女毛片av久久网站| 在线观看舔阴道视频| 欧美精品人与动牲交sv欧美| 少妇粗大呻吟视频| 国产男人的电影天堂91| 日本精品一区二区三区蜜桃| 亚洲精品国产av成人精品| 午夜福利在线观看吧| 亚洲国产成人一精品久久久| 777米奇影视久久| 国产伦理片在线播放av一区| 91九色精品人成在线观看| 电影成人av| 国产高清视频在线播放一区 | 1024香蕉在线观看| 亚洲av成人一区二区三| 久久天躁狠狠躁夜夜2o2o| 国产欧美日韩一区二区三区在线| 亚洲综合色网址| 99精品欧美一区二区三区四区| 丰满少妇做爰视频| 又大又爽又粗| av天堂在线播放| 一级毛片电影观看| 色94色欧美一区二区| 下体分泌物呈黄色| 最黄视频免费看| 日本欧美视频一区| 一级毛片女人18水好多| 国产免费福利视频在线观看| 国产精品久久久久久精品古装| 最新在线观看一区二区三区| 黄片小视频在线播放| 午夜91福利影院| 中国美女看黄片| 日韩制服丝袜自拍偷拍| 免费av中文字幕在线| 国产一区二区激情短视频 | 国产欧美日韩一区二区精品| 天天影视国产精品| xxxhd国产人妻xxx| 51午夜福利影视在线观看| 啦啦啦啦在线视频资源| 成人国语在线视频| 免费在线观看完整版高清| 人妻人人澡人人爽人人| 国产精品偷伦视频观看了| 国产野战对白在线观看| 9色porny在线观看| 久久人妻熟女aⅴ| 热re99久久精品国产66热6| 久久久久久亚洲精品国产蜜桃av| 久久天堂一区二区三区四区| 丰满少妇做爰视频| 国产黄频视频在线观看| 国产精品麻豆人妻色哟哟久久| 中国国产av一级| 亚洲九九香蕉| 精品国产超薄肉色丝袜足j| 亚洲精品乱久久久久久| 国产欧美日韩精品亚洲av| 五月天丁香电影| 麻豆国产av国片精品| 青春草视频在线免费观看| 精品国产乱子伦一区二区三区 | 另类精品久久| 国产男人的电影天堂91| 亚洲精品国产精品久久久不卡| 五月天丁香电影| 在线观看免费午夜福利视频| 精品人妻在线不人妻| 丝袜美足系列| 欧美亚洲 丝袜 人妻 在线| 丁香六月天网| 欧美久久黑人一区二区| 亚洲欧美一区二区三区黑人| 久久久久久久大尺度免费视频| 午夜激情久久久久久久| 亚洲av日韩精品久久久久久密| 久热爱精品视频在线9| 性少妇av在线| 欧美激情 高清一区二区三区| 中文字幕精品免费在线观看视频| 波多野结衣一区麻豆| 精品视频人人做人人爽| 精品少妇内射三级| 欧美亚洲 丝袜 人妻 在线| 一个人免费看片子| av国产精品久久久久影院| 老鸭窝网址在线观看| 亚洲精品一区蜜桃| 欧美久久黑人一区二区| 狠狠狠狠99中文字幕| 狠狠精品人妻久久久久久综合| 久久久精品94久久精品| 精品国产一区二区三区久久久樱花| 久久久精品区二区三区| 亚洲综合色网址| 日韩熟女老妇一区二区性免费视频| 一级黄色大片毛片| 一本—道久久a久久精品蜜桃钙片| 99国产精品免费福利视频| 手机成人av网站| 精品久久久久久电影网| 深夜精品福利| 丝袜脚勾引网站| 亚洲国产av新网站| 欧美黑人欧美精品刺激| 下体分泌物呈黄色| 操美女的视频在线观看| 男女床上黄色一级片免费看| 岛国在线观看网站| 亚洲国产av新网站| 亚洲成av片中文字幕在线观看| 一边摸一边抽搐一进一出视频| 精品国产一区二区久久| 美女高潮喷水抽搐中文字幕| 老司机影院成人| 国产极品粉嫩免费观看在线| 欧美老熟妇乱子伦牲交| 女性被躁到高潮视频| 国产麻豆69| 中文字幕人妻丝袜一区二区| 97人妻天天添夜夜摸| 成人国产av品久久久| 亚洲精品中文字幕一二三四区 | 天天添夜夜摸| 美女视频免费永久观看网站| 免费av中文字幕在线| 中文字幕人妻丝袜制服| 999久久久国产精品视频| 成人18禁高潮啪啪吃奶动态图| 亚洲av日韩在线播放| 乱人伦中国视频| 精品少妇久久久久久888优播| 亚洲精华国产精华精| 欧美av亚洲av综合av国产av| 午夜福利视频精品| 亚洲美女黄色视频免费看|