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

    一種分布式衛(wèi)星系統(tǒng)基線優(yōu)化設(shè)計(jì)方法

    2023-08-09 13:35:42胡瑞賢
    關(guān)鍵詞:徑向速度蒙特卡羅定位精度

    胡瑞賢, 張 昭, 駱 成

    (中國(guó)電子科技集團(tuán)公司電子科學(xué)研究院, 北京 100041)

    0 引 言

    星載合成孔徑雷達(dá)(synthetic aperture radar, SAR)[1]不受氣象條件、光照條件和國(guó)界的限制,可對(duì)地面、海面進(jìn)行常態(tài)化觀測(cè),受到各國(guó)廣泛關(guān)注。開(kāi)展星載地面動(dòng)目標(biāo)指示(ground moving target indicator, GMTI)技術(shù),突破其中關(guān)鍵技術(shù)難題,在軍事和民用領(lǐng)域具有廣闊的應(yīng)用前景。以美國(guó)為例,美國(guó)空軍快速能力辦公室于2018年啟動(dòng)了GMTI雷達(dá)衛(wèi)星研究項(xiàng)目,旨在通過(guò)由小型GMTI雷達(dá)衛(wèi)星組成的星座偵察監(jiān)視地面動(dòng)目標(biāo)。2022年4月,美國(guó)太空軍公布了天基GMTI能力開(kāi)發(fā)項(xiàng)目的進(jìn)展情況,期望通過(guò)實(shí)施天基GMTI雷達(dá)衛(wèi)星項(xiàng)目,大幅提升美軍的實(shí)時(shí)情報(bào)、監(jiān)視與偵察能力。然而由于衛(wèi)星的高速運(yùn)動(dòng)使地面不同位置的散射體相對(duì)于雷達(dá)有不同的速度,在天線波束的照射范圍內(nèi),雜波譜嚴(yán)重展寬[1-3]。如果采用單顆衛(wèi)星[4-7],天線面積和沿航向長(zhǎng)度受到衛(wèi)星有效載荷限制,系統(tǒng)最小可檢測(cè)速度(minimum detectable velocity, MDV)很大,無(wú)法有效檢測(cè)低速目標(biāo)。分布式衛(wèi)星系統(tǒng)[8-10]利用多顆小衛(wèi)星編隊(duì)飛行、可實(shí)現(xiàn)大衛(wèi)星的功能,具有制造成本低、發(fā)射費(fèi)用低、抗毀性強(qiáng)、可靠性高的顯著優(yōu)點(diǎn)。其有效載荷雷達(dá)協(xié)同工作可以共同完成單項(xiàng)或多項(xiàng)任務(wù)。例如,采用分布式衛(wèi)星星座可以實(shí)現(xiàn)衛(wèi)星組網(wǎng)通信和增大SAR成像[11]面積覆蓋率。星鏈衛(wèi)星在2022年俄烏戰(zhàn)爭(zhēng)中為烏克蘭提供了高速的網(wǎng)絡(luò)服務(wù),標(biāo)志著衛(wèi)星組網(wǎng)通信技術(shù)開(kāi)始應(yīng)用于實(shí)戰(zhàn)。在遙感領(lǐng)域,美國(guó)的黑杰克項(xiàng)目在開(kāi)發(fā)低成本快速組網(wǎng)偵察監(jiān)視衛(wèi)星。目前投入實(shí)際運(yùn)營(yíng)的分布式衛(wèi)星星座數(shù)量不多,例如德國(guó)的TerraSAR-X和TanDEM-X星座。對(duì)于動(dòng)目標(biāo)檢測(cè),分布式雷達(dá)在同一時(shí)刻有多個(gè)相位中心協(xié)同工作,可以提供較大的橫向孔徑,結(jié)合空時(shí)自適應(yīng)處理[2-3,12-14](space time adaptive processing, STAP)和去模糊處理[15-18],將顯著降低雷達(dá)的MDV、提高目標(biāo)定位精度,使雷達(dá)系統(tǒng)獲得較高的GMTI性能[19]。

    系統(tǒng)GMTI性能由系統(tǒng)輸出信雜噪比(signal-to-clutter-plus-noise ratio, output SCNR)、MDV、盲速和目標(biāo)定位精度等指標(biāo)[20-22]描述。輸出SCNR表示雜波抑制后信號(hào)功率與雜波加噪聲功率的比值,它反映了系統(tǒng)抑制雜波并對(duì)目標(biāo)信號(hào)進(jìn)行積累的性能。MDV反映了雷達(dá)系統(tǒng)對(duì)進(jìn)入主瓣雜波或主瓣雜波附近的低速目標(biāo)的檢測(cè)能力,MDV值越小,系統(tǒng)對(duì)低速目標(biāo)的檢測(cè)能力越強(qiáng)。盲速問(wèn)題指回波的空時(shí)頻譜結(jié)構(gòu)與固定雜波相似的運(yùn)動(dòng)目標(biāo),經(jīng)過(guò)雜波抑制后會(huì)被消除或嚴(yán)重衰減;在檢測(cè)此類目標(biāo)時(shí),將會(huì)產(chǎn)生漏警。目標(biāo)定位精度是衡量運(yùn)動(dòng)目標(biāo)重新標(biāo)定到真實(shí)位置時(shí)是否準(zhǔn)確的指標(biāo)。定位精度差會(huì)導(dǎo)致信息誤判,不利于實(shí)現(xiàn)對(duì)運(yùn)動(dòng)目標(biāo)的跟蹤和精確打擊。SAR圖像中運(yùn)動(dòng)目標(biāo)的圖像會(huì)在方位上出現(xiàn)位置偏差,其偏移量和目標(biāo)徑向速度有關(guān)。運(yùn)動(dòng)目標(biāo)的位置不能從圖像中直接讀取,需重新進(jìn)行測(cè)量。

    雷達(dá)系統(tǒng)參數(shù)和天線構(gòu)型的選擇對(duì)于雷達(dá)GMTI性能起著舉足輕重的作用。通常需要將以上多項(xiàng)指標(biāo)進(jìn)行綜合分析,通過(guò)各種方案、參數(shù)情況下的性能比較得到性能較好或性價(jià)比較高的方案。文獻(xiàn)[23]以獲得最優(yōu)速度響應(yīng)為目的,提出一種基線優(yōu)化設(shè)計(jì)方法。該方法根據(jù)“短的沿航向基線能解決盲速問(wèn)題,而長(zhǎng)的沿航向基線對(duì)低速目標(biāo)敏感”的特點(diǎn),提出將單星上多個(gè)接收孔徑間的短基線和多星間長(zhǎng)基線相結(jié)合,以改善系統(tǒng)速度響應(yīng)。這種方法的不足之處在于如果星間基線過(guò)長(zhǎng),盲速問(wèn)題較嚴(yán)重,搜索結(jié)果可能不夠精確,另外該方法未考慮定位精度對(duì)基線選擇的影響。事實(shí)上,沿航向基線與目標(biāo)定位精度關(guān)系密切。Ender等[24]和Gierull等[25]指出,借助分布式衛(wèi)星結(jié)構(gòu)獲得長(zhǎng)的沿航向基線能極大地提高目標(biāo)定位精度。本文在對(duì)多顆衛(wèi)星接收天線進(jìn)行子陣分割[26-28]的基礎(chǔ)上,通過(guò)理論分析推導(dǎo)了目標(biāo)定位精度與雷達(dá)系統(tǒng)參數(shù)之間的關(guān)系,揭示了基線長(zhǎng)度和參數(shù)改變系統(tǒng)速度響應(yīng)和目標(biāo)定位精度的原理。針對(duì)目標(biāo)定位精度與盲速對(duì)系統(tǒng)基線的不同要求,本文提出了一種同時(shí)兼顧速度響應(yīng)和目標(biāo)定位精度的基線設(shè)計(jì)方法,并研究了兩種誤差因素對(duì)分布式衛(wèi)星星座GMTI性能的影響,一種是單顆衛(wèi)星內(nèi)部多個(gè)天線子孔徑之間的基線誤差,另一種是衛(wèi)星之間的基線誤差。

    本文第1節(jié)提出了分布式星載SAR系統(tǒng)的幾何構(gòu)型并分析了其輸出SCNR。第2節(jié)推導(dǎo)了運(yùn)動(dòng)目標(biāo)定位表達(dá)式。由第1節(jié)和第2節(jié)分析可知,輸出SCNR和定位誤差都是基線的函數(shù)。在第3節(jié)提出了一種星載SAR-GMTI星座最優(yōu)基線設(shè)計(jì)方法,該方法同時(shí)考慮運(yùn)動(dòng)目標(biāo)速度響應(yīng)和定位誤差。第4節(jié)對(duì)該方法的性能進(jìn)行了仿真分析,第5節(jié)詳細(xì)分析了基線誤差對(duì)本方法結(jié)果的影響。第6節(jié)對(duì)本文工作進(jìn)行了總結(jié)。

    1 多天線輸出SCNR

    圖1給出了單顆衛(wèi)星和由兩顆衛(wèi)星構(gòu)成的分布式衛(wèi)星系統(tǒng)的幾何示意圖。假定衛(wèi)星天線工作在正側(cè)視情況,L為衛(wèi)星大天線相位中心間距,其速度均為Vs。圖中大天線在方位向上均勻劃分得到多個(gè)子天線,相鄰子天線的相位中心間距為d。假設(shè)在地面上存在一個(gè)勻速運(yùn)動(dòng)的理想點(diǎn)目標(biāo),其徑向速度為Vr,且該目標(biāo)與天線之間的距離為Rc。為了最大限度利用星載SAR功率孔徑積,考慮全孔徑標(biāo)準(zhǔn)發(fā)射(standard transmission, STTX),以及如下兩種方式接收:① 單顆衛(wèi)星均勻分成5個(gè)子孔徑同時(shí)接收這種構(gòu)型用SPAN5[28]表示),相鄰子孔徑之間的基線為d;② 兩顆衛(wèi)星均勻分成10個(gè)子孔徑同時(shí)接收(這種構(gòu)型用SPAN10[28]表示),星間基線為L(zhǎng)。為便于討論,后面小節(jié)里兩種構(gòu)型縮寫(xiě)為STTX-SPAN5和STTX-SPAN10。

    圖1 幾何構(gòu)型示意圖Fig.1 Skematic diagram of geometric configuration

    (1)

    (2)

    當(dāng)β=1時(shí),此時(shí)對(duì)應(yīng)的徑向速度為盲速,式(1)可以表示為

    (3)

    2 目標(biāo)定位誤差分析

    不同子孔徑接收信號(hào)表達(dá)式為

    yi=Si+ωi=αexp(jφ)exp(jφi/2)+ωi

    (4)

    假設(shè)每個(gè)子孔徑接收的運(yùn)動(dòng)目標(biāo)信號(hào)yi相互獨(dú)立,且高斯白噪聲ωi滿足獨(dú)立同分布,則推導(dǎo)得到變量α、φ和Vr的條件概率函數(shù)[29-30]表達(dá)式如下:

    (5)

    式中:(·)*表示共軛,且

    Ω=yi-αexp(jφ)exp(jφi)

    (6)

    條件概率函數(shù)f(y1,y2,…,yN|α,φ,Vr)的對(duì)數(shù)似然函數(shù)表達(dá)式如下:

    (7)

    (8)

    式中:

    (9)

    (10)

    考慮到白噪聲服從獨(dú)立同分布,可以得到

    (11)

    任何無(wú)偏估計(jì)方差均不小于CRB的倒數(shù),即有

    (12)

    為進(jìn)一步討論,令

    (13)

    式中:σVr表示徑向速度的標(biāo)準(zhǔn)差。根據(jù)σVr和文獻(xiàn)[31],可以得到運(yùn)動(dòng)目標(biāo)定位誤差表達(dá)式:

    (14)

    式中:Rc表示衛(wèi)星和運(yùn)動(dòng)目標(biāo)之間的斜距。

    對(duì)式(14)進(jìn)行討論。第一,此處得到的目標(biāo)定位誤差為誤差邊界,揭示了基線和定位誤差之間的關(guān)系,表示可以得到的最優(yōu)目標(biāo)定位結(jié)果。第二,式(14)中定位誤差ΔX與σVr成正比??紤]到當(dāng)前國(guó)內(nèi)外在軌運(yùn)行的星載SAR系統(tǒng)的軌道高度范圍為350~1 200 km,即使微小的測(cè)速誤差在幾百千米處的定位誤差都可達(dá)幾百米。由于徑向速度標(biāo)準(zhǔn)差σVr與d成反比,增大天線孔徑間距d可以提高定位精度,單個(gè)天線由于子孔徑間距限制導(dǎo)致目標(biāo)定位精度較差。

    3 最優(yōu)基線搜索方法

    由第1節(jié)和第2節(jié)中的分析可知,輸出SCNR和目標(biāo)定位精度均為基線的函數(shù)。通過(guò)基線優(yōu)化設(shè)計(jì),可以影響系統(tǒng)GMTI性能。最優(yōu)基線搜索方法是指對(duì)于給定軌道高度、運(yùn)行速度、衛(wèi)星數(shù)量、功率孔徑積的分布式衛(wèi)星星座系統(tǒng),在感興趣的目標(biāo)徑向速度區(qū)間內(nèi),按照一定的優(yōu)化準(zhǔn)則,針對(duì)不同的基線配置得到的輸出SCNR、MDV和目標(biāo)定位誤差情況進(jìn)行綜合比較,確定多顆衛(wèi)星之間基線以及單顆衛(wèi)星多個(gè)子孔徑之間基線最優(yōu)配置的方法。根據(jù)上述定義,下面提出一種將輸出SCNR和目標(biāo)定位精度相結(jié)合的最優(yōu)基線搜索方法。為了獲取滿意的搜索性能,在搜索最優(yōu)基線時(shí)需滿足兩條準(zhǔn)則:① 在感興趣的速度范圍內(nèi),輸出SCNR小于某一檢測(cè)門(mén)限的速度為盲速;② 使感興趣的速度范圍內(nèi),盲速區(qū)間與定位誤差區(qū)間的速度范圍盡可能小。

    根據(jù)以上準(zhǔn)則,提出基線搜索方法:首先設(shè)定搜索的速度區(qū)間(-Vr_max,Vr_max)、基線長(zhǎng)度范圍、輸出SCNR門(mén)限SCNRth和定位精度門(mén)限ΔX_th。然后,針對(duì)不同的天線分割情況,在(-Vr_max,Vr_max)區(qū)間內(nèi)用式(1)和式(14)分別計(jì)算每種基線配置下系統(tǒng)的輸出SCNR和目標(biāo)定位精度。最后,對(duì)基線的性能進(jìn)行評(píng)估,并分別找出盲速區(qū)間(-Vr1,Vr1)和存在定位誤差的速度區(qū)間(-Vr2,Vr2),使這兩個(gè)區(qū)間盡可能小的基線即為最優(yōu)基線。

    關(guān)于盲速區(qū)間的說(shuō)明:首先,在(-Vr1,Vr1)以外,即區(qū)間(-Vr_max,-Vr1)和(Vr1,Vr_max)內(nèi),不存在盲速點(diǎn);在區(qū)間(-Vr1,Vr1)內(nèi)也可能存在著單個(gè)甚至多個(gè)速度模糊點(diǎn),但是由于無(wú)法準(zhǔn)確測(cè)量此區(qū)間內(nèi)的目標(biāo)速度,故而此區(qū)間也屬于盲速區(qū)間。

    圖2 確定盲速區(qū)域的方法示意圖Fig.2 Schematic diagram of the method of determining the region of blind velocity

    為了取得滿意的定位性能,首先要求目標(biāo)可以被檢測(cè)到;其次要求目標(biāo)定位精度滿足條件。為了獲得優(yōu)良的目標(biāo)檢測(cè)性能,則應(yīng)當(dāng)使輸出SCNR和目標(biāo)定位精度兩種情況的重合區(qū)間盡可能大。假設(shè)在缺乏先驗(yàn)信息條件下,目標(biāo)徑向速度在(-Vr_max,Vr_max)內(nèi)滿足均勻分布[32],則其概率密度函數(shù)為

    (15)

    令Vr_def=max(|Vr1|,|Vr2|),其中max(X,Y)表示兩個(gè)元素X和Y的最大值。定義區(qū)間A=(-Vr_def≤Vr≤Vr_def)表示(-Vr1,Vr1)和(-Vr2,Vr2)最大重疊區(qū)域,且P(A)的表達(dá)式如下:

    (16)

    式(16)表明,為求解最優(yōu)基線,需要對(duì)不同構(gòu)型的性能概率進(jìn)行比較,速度區(qū)間不重合概率最低時(shí)所對(duì)應(yīng)的基線長(zhǎng)度即為最優(yōu)基線。

    (17)

    4 仿真結(jié)果分析

    仿真參數(shù)如表1所示。為防止碰撞,衛(wèi)星間距L不能太短(比如小于50 m),也不能過(guò)長(zhǎng),否則各天線波束照射重合區(qū)域過(guò)小。本文仿真取1 000 m為分布式衛(wèi)星間距的上限。此外,假設(shè)最大徑向速度Vmax為35 m/s,輸出SCNR檢測(cè)門(mén)限為SCNRo=13 dB,定位門(mén)限ΔX_th=10 m。

    表1 仿真參數(shù)

    為了說(shuō)明星間基線對(duì)系統(tǒng)性能的影響,令L取固定值(L=200 m)。圖3表示的是單星五天線和雙星十天線構(gòu)型的輸出SCNR隨目標(biāo)徑向速度的變化曲線。由圖3可以看出,由于圖2(a)單顆衛(wèi)星的最長(zhǎng)基線有限,在感興趣的速度范圍內(nèi)不存在盲速問(wèn)題。但是當(dāng)目標(biāo)徑向速度-6.42

    圖3 STTX-SPAN5和STTX-SPAN10的輸出SCNRFig.3 Output SCNR for STTX-SPAN5 and STTX-SPAN10

    圖3(b)表示雙星十天線構(gòu)型輸出SCNR隨Vr的變化曲線。在圖3(b)上取徑向速度區(qū)間[-5 m/s,5 m/s]和輸出SCNR區(qū)間[-5 dB,25 dB]進(jìn)行局部放大,得到圖3(c)。綜合圖3(b)和圖3(c)可知,由于天線總面積增加,同時(shí)相鄰子天線間的短基線和星間長(zhǎng)基線相結(jié)合,使得這種構(gòu)型比單星構(gòu)型對(duì)于低速區(qū)目標(biāo)更加敏感,可以檢測(cè)徑向速度較小的運(yùn)動(dòng)目標(biāo)。需要注意的是,長(zhǎng)基線的存在,使得低速區(qū)域出現(xiàn)了高增益的響應(yīng),很多速度區(qū)間獲得了高的輸出SCNR,但是由于基線過(guò)長(zhǎng)(長(zhǎng)度為200 m),這種區(qū)域是不連續(xù)的,此時(shí)干涉相位對(duì)速度變化更為敏感。如當(dāng)徑向速度范圍為-3.446 m/s

    圖4 STTX-SPAN5和STTX-SPAN10的定位誤差Fig.4 Location error for STTX-SPAN5 and STTX-SPAN10

    圖4(a)為STTX-SPAN5構(gòu)型的目標(biāo)定位誤差圖。受系統(tǒng)參數(shù)限制,STTX-SPAN5的目標(biāo)定位精度較差,最高精度只能達(dá)到50.89 m。在實(shí)際情況中,在缺乏全球定位系統(tǒng)、場(chǎng)景信息等輔助知識(shí)條件下,單星SAR-GMTI系統(tǒng)實(shí)現(xiàn)精確定位非常困難。這其中的根本原因在于波長(zhǎng)與天線孔徑的比值受限,波束較寬,分辨率較差。要提高定位精度,就要增加天線孔徑長(zhǎng)度,但在目前的技術(shù)條件下, 此方法的成本太高,因此必須借助分布式衛(wèi)星的星間長(zhǎng)基線。

    圖4(b)為STTX-SPAN10的目標(biāo)定位誤差曲線圖。同樣,衛(wèi)星間距L為200 m。在圖4(b)上取徑向速度區(qū)間[-7 m/s,7 m/s]和定位誤差區(qū)間[0 m,40 m]進(jìn)行局部放大,得到圖4(c)。綜合圖4(b)和圖4(c)可知,目標(biāo)定位誤差為1.26~34.65 m。當(dāng)徑向速度為0.28 m/s時(shí),低速區(qū)目標(biāo)定位精度最差可以達(dá)到34.65 m,而這種構(gòu)型下的最優(yōu)定位精度甚至可以達(dá)到1.26 m。與單星情況相比,定位精度提高了1個(gè)數(shù)量級(jí),定位誤差顯著降低,表明采用星間長(zhǎng)基線可以明顯提升目標(biāo)定位精度。

    綜合圖3和圖4的結(jié)果,當(dāng)基線L=200 m時(shí),按照第3節(jié)所述的最優(yōu)基線搜索方法,綜合確定不可靠的速度區(qū)間為(-3.446 m/s,3.446 m/s)。星間基線長(zhǎng)度改變,相應(yīng)的輸出SCNR和定位誤差精度也隨之發(fā)生變化,因此需要通過(guò)搜索、比較不同方案的優(yōu)劣來(lái)確定最終的基線結(jié)果。由表1中的衛(wèi)星基線范圍,采用第3節(jié)所述方法搜索得到雙星基線配置結(jié)果如圖5所示。需要強(qiáng)調(diào)的是,下面結(jié)果基于理想假設(shè)(不考慮衛(wèi)星姿態(tài)誤差、衛(wèi)星位置誤差、軌道與隊(duì)形誤差、電波在低層大氣和電離層中的傳播擾動(dòng)的影響等)。從圖5(a)可以看出,當(dāng)星間基線長(zhǎng)度為101.4 m時(shí),計(jì)算所得的概率最低,即最優(yōu)基線長(zhǎng)度為101.4 m。

    圖5 最優(yōu)基線搜索結(jié)果Fig.5 Optimal baseline searching result

    在這種基線配置條件下,輸出SCNR和定位精度結(jié)果如圖5(b)和5(c)所示,其不可靠的速度區(qū)間為(-2.732 m/s,2.732 m/s)。通過(guò)比較可知,該結(jié)果檢測(cè)性能最優(yōu)。需要指出的是,搜索步長(zhǎng)越精細(xì),運(yùn)算量越大,所需要的時(shí)間越長(zhǎng),越需要權(quán)衡穩(wěn)健性和運(yùn)算量之間的關(guān)系。搜索步長(zhǎng)太小,則沒(méi)有必要。

    5 基線誤差影響分析

    第4節(jié)介紹了最優(yōu)基線的結(jié)果。事實(shí)上,衛(wèi)星系統(tǒng)的GMTI性能受多種因素的影響(例如天線單元安裝精度、天線方向圖一致性、不同天線單元之間的反饋與互耦現(xiàn)象等)。由此引申出的問(wèn)題是,需要進(jìn)一步研究第4節(jié)介紹的方法的最優(yōu)基線結(jié)果是否具有穩(wěn)健性。考慮到本文的關(guān)鍵問(wèn)題是最優(yōu)基線搜索,此處重點(diǎn)研究基線誤差的影響。其他一些因素的影響參見(jiàn)文獻(xiàn)[33-36],在此不做討論。一般情況下,基線誤差是由材料熱脹冷縮的現(xiàn)象以及測(cè)量精度的限制等原因造成的。對(duì)于基線誤差,本節(jié)考慮兩種情況,即單星多天線之間的基線誤差(the multiple apertures baseline error inside a single satellite, MABEISS)和星間基線誤差(inter-satellite baseline error, ISBE)。本節(jié)主要研究基線誤差對(duì)輸出SCNR及定位誤差的影響。最重要的是,將討論MABEISS和ISBE對(duì)最優(yōu)基線結(jié)果的影響。

    5.1 基線誤差對(duì)輸出SCNR的影響

    首先,分析基線誤差對(duì)輸出SCNR的影響。在上面分析的基礎(chǔ)上,存在基線誤差時(shí),輸出SCNR[2-3]的表達(dá)式可以修正為

    (18)

    式中:Sreal是運(yùn)動(dòng)目標(biāo)實(shí)際導(dǎo)向矢量;Suse是匹配濾波時(shí)的測(cè)量導(dǎo)向矢量。Suse可以表示為

    Suse=α·[exp(jφ1),exp(jφ2),…,exp(jφP)]H

    (19)

    (20)

    第三步建立方程(3),檢驗(yàn)加入中介變量后文化與旅游產(chǎn)業(yè)融合的直接效應(yīng)(γ1是否顯著,是否與β1γ7同號(hào))。

    將式(19)和式(20)代入式(18),可以得到

    (21)

    式中:

    (22)

    從式(21)可知,存在基線誤差時(shí),有Sreal≠Suse,則由于導(dǎo)向矢量不匹配導(dǎo)致輸出SCNR下降。

    下面通過(guò)仿真分析基線誤差對(duì)輸出SCNR的影響。仿真參數(shù)如表1所示。在仿真中,考慮到實(shí)際情況,假設(shè)基線誤差δdi服從修正的高斯分布,其上限和下限分別為Δd和-Δd。

    下面討論MABEISS和ISBE對(duì)輸出SCNR的影響。值得注意的是,本節(jié)為避免這兩個(gè)因素相互干擾,在研究MABEISS影響時(shí),令I(lǐng)SBE取值為零,反之亦然。在下面討論的所有情況下都是類似的,后面不再贅述。

    5.1.1 MABEISS的影響

    下面仿真了MABEISS的3種情況。第1種情況是Δd=3 mm(其中,Δd表示基線誤差的界限),其結(jié)果如圖6(a)和圖6(d)所示。如圖6(a)所示,只進(jìn)行了一次蒙特卡羅實(shí)驗(yàn),曲線表示輸出SCNR差值(測(cè)量SCNR與真實(shí)輸出SCNR之間的差值)隨徑向速度變化。在圖6(d)中,進(jìn)行了2 000次蒙特卡羅實(shí)驗(yàn),圖6中曲線表示輸出SCNR差值隨徑向速度的變化。值得注意的是,曲線代表了2 000個(gè)實(shí)驗(yàn)樣本的絕對(duì)值的平均值。第2種情況是Δd=5 mm,如圖6(b)和圖6(e)所示。第3種情況是Δd=10 mm,對(duì)應(yīng)圖6(c)和圖6(f)。

    圖6 MABEISS輸出SCNR誤差隨徑向速度變化Fig.6 Output SCNR error variation with radial velocity for MABEISS

    從圖6結(jié)果可知,輸出SCNR誤差在10-3dB量級(jí)。圖6(a)~圖6(c)給出了一次蒙特卡羅實(shí)驗(yàn)的結(jié)果。需要注意的是,該結(jié)果很大程度上取決于一次隨機(jī)基線誤差。例如,一次Δd=3 mm的蒙特卡羅實(shí)驗(yàn)可能比Δd=5 mm的蒙特卡羅實(shí)驗(yàn)導(dǎo)致更大的定位誤差。為平滑實(shí)驗(yàn)結(jié)果,應(yīng)開(kāi)展多次蒙特卡羅實(shí)驗(yàn)。圖6(d)~圖6(f)給出了2 000次蒙特卡羅實(shí)驗(yàn)結(jié)果,從而可以發(fā)現(xiàn)輸出SCNR誤差從圖6(d)到圖6(f)越來(lái)越大。與圖6(d)~圖6(f)相比,輸出SCNR誤差與基線誤差界限成正比。例如,圖6(d)~圖6(f)的最大值分別為0.002 0 dB、0.003 3 dB和0.006 6 dB,與基線誤差范圍3 mm、5 mm和10 mm基本成正比。需要指出,大量實(shí)驗(yàn)(例如遠(yuǎn)遠(yuǎn)超過(guò)2 000次的實(shí)驗(yàn))雖然有利于獲得更準(zhǔn)確的結(jié)果,但非常耗時(shí),需要在計(jì)算量和準(zhǔn)確性之間進(jìn)行權(quán)衡,2 000次基本滿足需要。

    綜上所述,可以從圖6中看出10-3m量級(jí)的MABEISS對(duì)輸出SCNR的影響相對(duì)較小。即使MABEISS大到10 mm,對(duì)輸出SCNR的影響幾乎可以忽略。因此,10 mm以下MABEISS對(duì)星載GMTI系統(tǒng)運(yùn)動(dòng)目標(biāo)的速度響應(yīng)影響較小。

    5.1.2 ISBE的影響

    由于現(xiàn)有測(cè)量方法的測(cè)量精度限制,ISBE通常大于MABEISS。本節(jié)考慮ISBE的3種情況。第1種情況,如圖7(a)和圖7(d)所示,為Δd=2 cm(Δd表示基線誤差的界限)對(duì)應(yīng)的輸出SCNR誤差。在圖7(a)中,只進(jìn)行了一次蒙特卡羅實(shí)驗(yàn),圖7中曲線表示輸出SCNR誤差隨徑向速度的變化情況。在圖7(d)中,重復(fù)了2 000次蒙特卡羅實(shí)驗(yàn),其中圖7(d)中的曲線表示2 000次輸出SCNR誤差的平均值的絕對(duì)值。第2種情況,如圖7(b)和圖7(e)所示,為Δd=10 cm對(duì)應(yīng)的輸出SCNR誤差。第3種情況,如圖7(c)和圖7(f)所示,為Δd=100 cm對(duì)應(yīng)的輸出SCNR誤差。

    圖7 ISBE輸出SCNR誤差隨徑向速度變化Fig.7 Output SCNR error variation with radial velocity for ISBE

    圖7中輸出SCNR誤差結(jié)果基本在10-3dB或102dB左右。圖7(a)~圖7(c)給出了一次蒙特卡羅實(shí)驗(yàn)的結(jié)果。圖7(d)~圖7(f)給出了2 000次蒙特卡羅實(shí)驗(yàn)結(jié)果,可以發(fā)現(xiàn)定位誤差從圖7(d)到圖7(f)越來(lái)越大。圖7(d)~圖7(f)給出了2 000次蒙特卡羅實(shí)驗(yàn)結(jié)果,可以發(fā)現(xiàn)輸出SCNR從圖7(d)到圖7(f)越來(lái)越大。與圖7(d)~圖7(f)相比,輸出SCNR與基線誤差界限成正比。例如,圖7(d)~圖7(f)的最大值分別為0.007 8 dB、0.039 5 dB和0.394 60 dB,它們與基線誤差范圍0.02 m、0.10 m和1.00 m基本成正比。

    從上述實(shí)驗(yàn)結(jié)果可以發(fā)現(xiàn),ISBE對(duì)輸出SCNR的影響大于MABEISS。此外,根據(jù)圖7(a)~圖7(f),可以得出如下結(jié)論:

    (1) ISBE為厘米量級(jí)時(shí),輸出SCNR誤差約為10-2dB量級(jí),此時(shí)ISBE對(duì)GMTI性能的影響可忽略不計(jì)。

    (2) ISBE為分米量級(jí)時(shí),輸出SCNR誤差約為0.1 dB,此時(shí)ISBE對(duì)GMTI性能影響較小。

    (3) ISBE達(dá)到米級(jí)時(shí),輸出SCNR誤差約為1 dB,此時(shí)ISBE對(duì)GMTI性能的影響不可忽略。需要開(kāi)展專門(mén)的設(shè)計(jì)工作降低星間基線誤差。

    5.2 基線誤差對(duì)定位誤差的影響

    下面討論基線誤差對(duì)定位精度的影響。根據(jù)沿航向干涉儀方法(along-track interferometry, ATI)[32,37-42],可以得到存在基線誤差時(shí)徑向速度表達(dá)式:

    (23)

    式中:

    (24)

    定位誤差表達(dá)式為

    (25)

    由基線誤差引起的徑向速度估計(jì)誤差表達(dá)式為

    (26)

    定位誤差表達(dá)式如下:

    (27)

    下面仿真分析基線誤差對(duì)定位誤差的影響,仿真參數(shù)如表1所示,此時(shí)衛(wèi)星間距取值仍為200 m。

    5.2.1 MABEISS的影響

    圖8中定位誤差的結(jié)果基本在10-3m量級(jí)。圖8(a)~圖8(c)給出了一次蒙特卡羅實(shí)驗(yàn)的結(jié)果。圖8(d)~圖8(f)給出了2 000次蒙特卡羅實(shí)驗(yàn)結(jié)果,可以發(fā)現(xiàn)定位誤差從圖8(d)到圖8(f)越來(lái)越大。與圖8(d)~圖8(f)相比,定位誤差與基線誤差界限成正比。例如,圖8(d)~圖8(f)定位誤差的最大值分別為0.004 2 m、0.007 3 m和0.014 0 m,它們與基線誤差范圍3 mm、5 mm和10 mm基本成正比。

    圖8 MABEISS定位誤差差值隨徑向速度變化Fig.8 Location error variation with radial velocity for MABEISS

    總之,可以得出結(jié)論,即MABEISS對(duì)圖8中定位誤差的影響很小。即使MABEISS達(dá)到10 mm,對(duì)定位誤差的影響也幾乎可以忽略不計(jì)。因此,可以得到如下結(jié)論:對(duì)于星載GMTI系統(tǒng),MABEISS對(duì)運(yùn)動(dòng)目標(biāo)的定位誤差差值影響不大。

    5.2.2 ISBE的影響

    圖9(a)~圖9(c)給出了一次蒙特卡羅實(shí)驗(yàn)的結(jié)果,圖9(d)~圖9(f)給出了2 000次蒙特卡羅實(shí)驗(yàn)的結(jié)果,圖9中定位誤差數(shù)值在10-2~10-1m量級(jí)。通過(guò)仿真分析可知,定位誤差與基線誤差邊界成正比。例如,圖9(d)~圖9(f)的最大值分別為0.044 1 m、0.220 7 m和2.250 5 m,它們與基線誤差邊界2 cm、10 cm和100 cm基本成正比。

    圖9 ISBE定位誤差差值隨徑向速度變化結(jié)果Fig.9 Location error verus radial velocity for ISBE

    總之,從圖9的仿真結(jié)果可知,ISBE對(duì)定位誤差差值的影響大于MABEISS。此外,根據(jù)圖9(a)~圖9(f),可以得出如下結(jié)論:

    (1) 當(dāng)ISBE為厘米量級(jí),定位誤差約為10-2m量級(jí),此時(shí)ISBE對(duì)GMTI性能的影響可忽略不計(jì)。

    (2) 當(dāng)ISBE為分米量級(jí),定位誤差約為0.1 m量級(jí),此時(shí)ISBE對(duì)GMTI性能的影響不大。

    (3) 當(dāng)ISBE為米量級(jí)時(shí),定位誤差約為1 m量級(jí),此時(shí)ISBE對(duì)GMTI性能的影響不可忽略。

    5.3 存在基線誤差時(shí)的最優(yōu)基線搜索結(jié)果

    本節(jié)分析存在基線誤差時(shí)本文方法的魯棒性并仿真了3種情況:① MABEISS=3 mm, ISBE=2 cm;② MABEISS=3 mm, ISBE=10 cm;③ MABEISS=3 mm, ISBE=100 cm。根據(jù)統(tǒng)計(jì)規(guī)律,采用500次蒙特卡羅實(shí)驗(yàn)搜索最優(yōu)基線。

    圖10為存在誤差時(shí)的最優(yōu)基線搜索結(jié)果。圖10(a)對(duì)應(yīng)基線誤差(MABEISS=3 mm,ISBE=2 cm)的最優(yōu)基線搜索結(jié)果,與圖5(a)中沒(méi)有基線誤差的結(jié)果幾乎相同。同樣,圖10(b)對(duì)應(yīng)基線誤差(MABEISS=3 mm,ISBE=10 cm)的結(jié)果,與圖5(a)的結(jié)果差異也較小。換言之,即使ISBE界限增大到10 cm,其影響仍然可以忽略不計(jì)。但是當(dāng)ISBE達(dá)到100 cm量級(jí)時(shí),可以觀察到第3種情況,即圖10(c)對(duì)應(yīng)基線誤差(MABEISS=3 mm,ISBE=100 cm)與圖5(a)存在明顯差異,此時(shí)得到的最佳基線為102 m。由此可知,當(dāng)ISBE達(dá)到1 m量級(jí)時(shí),最優(yōu)基線變化不超過(guò)1 m。

    圖10 存在誤差時(shí)的最優(yōu)基線搜索結(jié)果Fig.10 Optimal baseline searching results when baseline error exists

    綜上所述,從圖10中可以得出如下結(jié)論。

    (1) 目前MABEISS測(cè)量精度較高,對(duì)最優(yōu)基線搜索結(jié)果影響不大。

    (2) ISBE通常大于MABEISS,其影響需要根據(jù)不同情況進(jìn)行討論。

    (3) 對(duì)于雷達(dá)系統(tǒng)而言,細(xì)微的基線誤差不會(huì)造成嚴(yán)重影響,即一定程度的ISBE可以接受。

    由于影響較小,在幾米范圍內(nèi)調(diào)整沿軌道基線既浪費(fèi)能源,也不必要。衛(wèi)星間基線可以定期調(diào)整。因此,在進(jìn)行基線設(shè)計(jì)時(shí),基線長(zhǎng)度的選擇以符合實(shí)際情況要求為宜。

    6 結(jié) 論

    本文從輸出SCNR、MDV、盲速和目標(biāo)定位精度等描述系統(tǒng)GMTI性能的指標(biāo)入手,推導(dǎo)了基線長(zhǎng)度、系統(tǒng)參數(shù)與系統(tǒng)速度響應(yīng)和目標(biāo)定位精度之間的關(guān)系,在此基礎(chǔ)上提出了一種兼顧速度響應(yīng)和目標(biāo)定位精度的最優(yōu)基線設(shè)計(jì)方法,為分布式星載SAR系統(tǒng)GMTI性能的優(yōu)化提供了有益的探索。該方法可以在雷達(dá)系統(tǒng)投入使用前通過(guò)優(yōu)化設(shè)計(jì)得到最優(yōu)的GMTI性能,從而避免巨大投入后取得的收益無(wú)法滿足實(shí)際使用需求。另外,該方法在MABEISS和ISBE有限的情況下仍能獲得最優(yōu)基線,顯示出本文方法良好的魯棒性。此外,該方法可擴(kuò)展到3顆衛(wèi)星甚至多顆衛(wèi)星的基線優(yōu)化問(wèn)題,為解決分布式衛(wèi)星GMTI性能優(yōu)化提供重要參考。

    需要指出的是,盡管本文方法具有很強(qiáng)的理論意義,但在實(shí)際情況中還需要考慮地球自轉(zhuǎn)、大氣擾動(dòng)、平臺(tái)高速運(yùn)動(dòng)、衛(wèi)星軌道誤差和姿態(tài)誤差等各種因素[1, 33-36, 43-44],對(duì)以上方法進(jìn)行檢驗(yàn)和改進(jìn),還需要進(jìn)一步開(kāi)展研究工作。另外,用小衛(wèi)星群替代大衛(wèi)星仍存在著諸多問(wèn)題,如由基線去相干[45-47]、星間長(zhǎng)基線[48]造成的盲速、外界干擾[49]等,對(duì)此也有必要進(jìn)行深入研究。

    猜你喜歡
    徑向速度蒙特卡羅定位精度
    北斗定位精度可達(dá)兩三米
    軍事文摘(2023年4期)2023-04-05 13:57:35
    利用蒙特卡羅方法求解二重積分
    GPS定位精度研究
    組合導(dǎo)航的AGV定位精度的改善
    非圓形光纖研究進(jìn)展
    航空兵器(2017年6期)2018-01-24 15:00:10
    臺(tái)風(fēng)威馬遜造成云南文山州強(qiáng)降水天氣雷達(dá)回波分析
    探討蒙特卡羅方法在解微分方程邊值問(wèn)題中的應(yīng)用
    距離頻率ML方法無(wú)模糊估計(jì)動(dòng)目標(biāo)徑向速度
    復(fù)合型種子源125I-103Pd劑量場(chǎng)分布的蒙特卡羅模擬與實(shí)驗(yàn)測(cè)定
    同位素(2014年2期)2014-04-16 04:57:20
    基于蒙特卡羅仿真的CRC檢錯(cuò)能力驗(yàn)證
    黄色 视频免费看| 黑人猛操日本美女一级片| 国产免费福利视频在线观看| 亚洲成人av在线免费| 伦理电影大哥的女人| 不卡视频在线观看欧美| 黄色一级大片看看| 最近中文字幕2019免费版| 亚洲av电影在线进入| 国产伦人伦偷精品视频| 国产精品成人在线| 99热国产这里只有精品6| 成人国语在线视频| 婷婷色综合大香蕉| 最近最新中文字幕免费大全7| 夫妻午夜视频| 亚洲情色 制服丝袜| 国产片内射在线| 国产成人免费无遮挡视频| 国产成人a∨麻豆精品| 一本久久精品| 狂野欧美激情性bbbbbb| 蜜桃国产av成人99| 在线观看一区二区三区激情| 国产淫语在线视频| 欧美人与善性xxx| 国产精品av久久久久免费| 一级毛片黄色毛片免费观看视频| 国产精品 国内视频| 亚洲一区二区三区欧美精品| 国产精品一国产av| 日日撸夜夜添| 久久久精品94久久精品| 久久 成人 亚洲| 亚洲四区av| 久久久久久免费高清国产稀缺| 嫩草影视91久久| 日韩精品免费视频一区二区三区| 国产不卡av网站在线观看| 亚洲av日韩精品久久久久久密 | 精品国产一区二区三区四区第35| 亚洲成人国产一区在线观看 | 美女大奶头黄色视频| 又大又爽又粗| 少妇人妻久久综合中文| 桃花免费在线播放| 亚洲一级一片aⅴ在线观看| 桃花免费在线播放| 成人国语在线视频| 亚洲国产av新网站| svipshipincom国产片| 午夜福利视频精品| 一边摸一边做爽爽视频免费| 久久人人97超碰香蕉20202| 成人午夜精彩视频在线观看| 精品人妻熟女毛片av久久网站| 欧美激情 高清一区二区三区| 在线天堂中文资源库| 免费高清在线观看日韩| 日韩人妻精品一区2区三区| 亚洲人成电影观看| bbb黄色大片| 日本猛色少妇xxxxx猛交久久| 啦啦啦在线观看免费高清www| www.精华液| 中文字幕色久视频| 妹子高潮喷水视频| 老司机靠b影院| 国产一级毛片在线| tube8黄色片| 精品人妻熟女毛片av久久网站| 久久久亚洲精品成人影院| 欧美日韩成人在线一区二区| e午夜精品久久久久久久| av一本久久久久| 亚洲欧美一区二区三区久久| 9191精品国产免费久久| 中文字幕高清在线视频| 91老司机精品| 嫩草影院入口| 久久精品久久精品一区二区三区| 亚洲成人手机| av在线播放精品| 少妇猛男粗大的猛烈进出视频| av国产精品久久久久影院| 热re99久久国产66热| 日本欧美国产在线视频| 超碰97精品在线观看| 另类精品久久| 伊人久久国产一区二区| 午夜福利乱码中文字幕| 尾随美女入室| 欧美日韩综合久久久久久| 男女边摸边吃奶| 婷婷色麻豆天堂久久| 99热网站在线观看| 亚洲欧洲日产国产| 久久精品久久久久久久性| 亚洲国产精品999| 国产精品久久久久久人妻精品电影 | 久久天堂一区二区三区四区| 亚洲欧美激情在线| 国产老妇伦熟女老妇高清| 亚洲欧美一区二区三区国产| 亚洲精品在线美女| 国产在线一区二区三区精| 久久久久久人人人人人| h视频一区二区三区| 18禁国产床啪视频网站| 19禁男女啪啪无遮挡网站| 久久国产精品大桥未久av| 亚洲精品日韩在线中文字幕| 性少妇av在线| 久久人人爽av亚洲精品天堂| 国产老妇伦熟女老妇高清| 黑人巨大精品欧美一区二区蜜桃| 一区福利在线观看| 天天添夜夜摸| 国产精品偷伦视频观看了| 美女扒开内裤让男人捅视频| 久久97久久精品| 最近的中文字幕免费完整| 久久精品久久久久久久性| av在线app专区| 国产伦理片在线播放av一区| 亚洲av日韩精品久久久久久密 | 欧美激情 高清一区二区三区| 久久久国产精品麻豆| 大香蕉久久成人网| 国产精品亚洲av一区麻豆 | 不卡视频在线观看欧美| 大陆偷拍与自拍| 婷婷色综合大香蕉| 国产亚洲最大av| 国产在视频线精品| 免费黄色在线免费观看| 一级爰片在线观看| 色94色欧美一区二区| 国产精品偷伦视频观看了| 亚洲精品国产一区二区精华液| 国产成人a∨麻豆精品| 免费在线观看视频国产中文字幕亚洲 | 你懂的网址亚洲精品在线观看| 亚洲精品,欧美精品| 国产极品粉嫩免费观看在线| 最近最新中文字幕免费大全7| 中文精品一卡2卡3卡4更新| 亚洲欧洲精品一区二区精品久久久 | 老司机影院成人| 嫩草影视91久久| 欧美老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 日韩伦理黄色片| 人人妻人人添人人爽欧美一区卜| 水蜜桃什么品种好| 亚洲国产看品久久| 天堂中文最新版在线下载| 又大又爽又粗| 伦理电影免费视频| 国产成人a∨麻豆精品| 晚上一个人看的免费电影| 精品人妻熟女毛片av久久网站| 大香蕉久久成人网| 黄色怎么调成土黄色| 少妇人妻精品综合一区二区| 大陆偷拍与自拍| 欧美另类一区| 国产一级毛片在线| 国产成人精品无人区| 欧美久久黑人一区二区| www.自偷自拍.com| 亚洲人成电影观看| 亚洲国产最新在线播放| 国产av一区二区精品久久| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩精品网址| 精品国产一区二区三区久久久樱花| 美女扒开内裤让男人捅视频| 日日啪夜夜爽| 欧美日韩视频精品一区| 亚洲国产最新在线播放| 不卡视频在线观看欧美| 一级毛片电影观看| 日韩中文字幕视频在线看片| 18禁观看日本| 亚洲av日韩在线播放| 亚洲国产精品一区三区| 美女午夜性视频免费| 人体艺术视频欧美日本| √禁漫天堂资源中文www| 中文字幕制服av| 亚洲国产精品一区二区三区在线| 黄片播放在线免费| 精品国产一区二区三区四区第35| 中文字幕人妻熟女乱码| 日本欧美国产在线视频| 亚洲成av片中文字幕在线观看| 最新在线观看一区二区三区 | 韩国精品一区二区三区| 成人手机av| 国产在视频线精品| 欧美黄色片欧美黄色片| 国产男女内射视频| 国产在线一区二区三区精| a级毛片在线看网站| 国产免费视频播放在线视频| 亚洲视频免费观看视频| 亚洲国产欧美在线一区| 国产深夜福利视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲婷婷狠狠爱综合网| 老汉色av国产亚洲站长工具| 日韩av免费高清视频| 亚洲精品国产区一区二| 久久免费观看电影| 午夜91福利影院| 1024香蕉在线观看| 可以免费在线观看a视频的电影网站 | 香蕉丝袜av| 免费av中文字幕在线| 美女视频免费永久观看网站| 一级,二级,三级黄色视频| 黄色毛片三级朝国网站| 黄片小视频在线播放| 婷婷色麻豆天堂久久| 别揉我奶头~嗯~啊~动态视频 | 肉色欧美久久久久久久蜜桃| 巨乳人妻的诱惑在线观看| 国产av码专区亚洲av| 黑人欧美特级aaaaaa片| 久久久国产欧美日韩av| 如日韩欧美国产精品一区二区三区| √禁漫天堂资源中文www| 久久精品国产a三级三级三级| 久久久亚洲精品成人影院| 欧美另类一区| 极品少妇高潮喷水抽搐| 国产精品久久久久久人妻精品电影 | 亚洲第一区二区三区不卡| 日韩中文字幕视频在线看片| 蜜桃国产av成人99| 香蕉丝袜av| 久久午夜综合久久蜜桃| 久久精品国产a三级三级三级| 又大又爽又粗| 精品少妇内射三级| 亚洲图色成人| 中文欧美无线码| 人人澡人人妻人| 免费黄频网站在线观看国产| 婷婷色麻豆天堂久久| av又黄又爽大尺度在线免费看| 久久久久久免费高清国产稀缺| 黄片小视频在线播放| 黄色一级大片看看| 人人妻人人添人人爽欧美一区卜| 一区二区三区乱码不卡18| 国产精品二区激情视频| 国产麻豆69| 热re99久久精品国产66热6| 十八禁网站网址无遮挡| 午夜日本视频在线| 水蜜桃什么品种好| 国产有黄有色有爽视频| 国产精品欧美亚洲77777| 国产日韩欧美亚洲二区| a级毛片在线看网站| 99久久99久久久精品蜜桃| 亚洲国产精品国产精品| 欧美xxⅹ黑人| 亚洲欧美清纯卡通| 欧美在线黄色| 国产一区亚洲一区在线观看| 一边摸一边做爽爽视频免费| 亚洲国产精品一区二区三区在线| 综合色丁香网| 午夜福利一区二区在线看| 人成视频在线观看免费观看| 91精品伊人久久大香线蕉| 国产精品久久久久久久久免| 中文字幕另类日韩欧美亚洲嫩草| 国产黄色免费在线视频| 久久ye,这里只有精品| 婷婷色av中文字幕| 两个人免费观看高清视频| 国产成人精品久久二区二区91 | 91aial.com中文字幕在线观看| 亚洲情色 制服丝袜| 亚洲国产成人一精品久久久| 久久99一区二区三区| 国产精品欧美亚洲77777| 无遮挡黄片免费观看| 色精品久久人妻99蜜桃| 最新在线观看一区二区三区 | 亚洲欧洲精品一区二区精品久久久 | 一级,二级,三级黄色视频| 久久精品国产综合久久久| 精品人妻熟女毛片av久久网站| 精品一区在线观看国产| av视频免费观看在线观看| 色婷婷av一区二区三区视频| 亚洲欧美日韩另类电影网站| 交换朋友夫妻互换小说| 一本久久精品| 国产毛片在线视频| 丰满乱子伦码专区| 欧美av亚洲av综合av国产av | 亚洲精品国产av蜜桃| 大香蕉久久成人网| 久热爱精品视频在线9| 亚洲第一区二区三区不卡| 精品少妇内射三级| 国产成人精品在线电影| 欧美精品av麻豆av| 国产精品国产av在线观看| 日韩大片免费观看网站| 一个人免费看片子| 桃花免费在线播放| 91成人精品电影| 少妇被粗大猛烈的视频| 999久久久国产精品视频| 国产一区二区三区av在线| 亚洲欧美一区二区三区国产| 久久天躁狠狠躁夜夜2o2o | 你懂的网址亚洲精品在线观看| 中文字幕高清在线视频| 看十八女毛片水多多多| 午夜日韩欧美国产| 女人久久www免费人成看片| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 亚洲色图 男人天堂 中文字幕| 亚洲自偷自拍图片 自拍| 天天添夜夜摸| 国产在视频线精品| 成年美女黄网站色视频大全免费| 亚洲欧美一区二区三区久久| 精品少妇黑人巨大在线播放| 黄色怎么调成土黄色| 一本一本久久a久久精品综合妖精| 久久这里只有精品19| 男人舔女人的私密视频| 桃花免费在线播放| 精品亚洲成a人片在线观看| 亚洲欧洲国产日韩| 国产精品久久久久久精品古装| 日韩一本色道免费dvd| 中文字幕高清在线视频| 亚洲国产日韩一区二区| 欧美日韩一级在线毛片| 久久精品久久久久久久性| 欧美人与善性xxx| 亚洲一级一片aⅴ在线观看| 久久久久精品国产欧美久久久 | 黄片小视频在线播放| 免费黄频网站在线观看国产| 亚洲一卡2卡3卡4卡5卡精品中文| 天天躁夜夜躁狠狠躁躁| 色网站视频免费| 王馨瑶露胸无遮挡在线观看| 色94色欧美一区二区| 亚洲国产精品一区三区| 亚洲精品美女久久久久99蜜臀 | 日韩中文字幕视频在线看片| 侵犯人妻中文字幕一二三四区| 免费观看人在逋| www.熟女人妻精品国产| 午夜福利一区二区在线看| 欧美国产精品一级二级三级| 伊人久久国产一区二区| 可以免费在线观看a视频的电影网站 | 国产淫语在线视频| 欧美人与性动交α欧美精品济南到| 蜜桃国产av成人99| 观看av在线不卡| 女人爽到高潮嗷嗷叫在线视频| 女性生殖器流出的白浆| 久久热在线av| 国产精品国产av在线观看| 国产人伦9x9x在线观看| 在线观看一区二区三区激情| 亚洲美女视频黄频| 久久国产亚洲av麻豆专区| 19禁男女啪啪无遮挡网站| 免费不卡黄色视频| 一个人免费看片子| 亚洲精品国产一区二区精华液| 热re99久久精品国产66热6| 香蕉国产在线看| 国产精品一二三区在线看| 91精品三级在线观看| 亚洲色图综合在线观看| 少妇精品久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 日韩一区二区三区影片| 久久精品国产亚洲av高清一级| 一边摸一边抽搐一进一出视频| 女人爽到高潮嗷嗷叫在线视频| 欧美激情高清一区二区三区 | 男人添女人高潮全过程视频| 欧美国产精品一级二级三级| 亚洲国产精品成人久久小说| 久久人妻熟女aⅴ| 成人影院久久| 18禁国产床啪视频网站| 日本爱情动作片www.在线观看| 午夜福利乱码中文字幕| 两个人看的免费小视频| 老汉色av国产亚洲站长工具| 国产精品一区二区在线不卡| 国产探花极品一区二区| 中国国产av一级| 成年av动漫网址| 国产在线视频一区二区| 日本爱情动作片www.在线观看| 国产精品国产三级国产专区5o| 午夜av观看不卡| 欧美日韩亚洲高清精品| 成人手机av| 在线亚洲精品国产二区图片欧美| 菩萨蛮人人尽说江南好唐韦庄| 韩国高清视频一区二区三区| 国产国语露脸激情在线看| 国产老妇伦熟女老妇高清| 亚洲激情五月婷婷啪啪| 免费看av在线观看网站| 大话2 男鬼变身卡| 麻豆精品久久久久久蜜桃| 男女午夜视频在线观看| 美女扒开内裤让男人捅视频| 精品国产超薄肉色丝袜足j| 中文字幕人妻丝袜制服| 欧美国产精品va在线观看不卡| 熟妇人妻不卡中文字幕| 最近最新中文字幕免费大全7| 丰满饥渴人妻一区二区三| 制服诱惑二区| 夜夜骑夜夜射夜夜干| 色网站视频免费| 女性被躁到高潮视频| 青春草视频在线免费观看| 午夜久久久在线观看| 久久久久久免费高清国产稀缺| 国产欧美日韩综合在线一区二区| 久热爱精品视频在线9| 午夜av观看不卡| 久久午夜综合久久蜜桃| 久久久国产精品麻豆| 欧美 亚洲 国产 日韩一| 久久久久久久久久久久大奶| 亚洲精品久久午夜乱码| 女的被弄到高潮叫床怎么办| 少妇人妻 视频| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| 99久久精品国产亚洲精品| 日韩 欧美 亚洲 中文字幕| 欧美精品av麻豆av| 亚洲精品在线美女| 国产精品秋霞免费鲁丝片| 精品国产超薄肉色丝袜足j| 看免费成人av毛片| 夜夜骑夜夜射夜夜干| 免费看不卡的av| 成人漫画全彩无遮挡| 亚洲伊人色综图| 各种免费的搞黄视频| 久久久精品区二区三区| 天堂中文最新版在线下载| 最近的中文字幕免费完整| 香蕉国产在线看| 老司机深夜福利视频在线观看 | 午夜免费鲁丝| 在线观看免费高清a一片| 制服丝袜香蕉在线| 亚洲欧美一区二区三区国产| 哪个播放器可以免费观看大片| 日韩制服丝袜自拍偷拍| 国产成人欧美在线观看 | 色94色欧美一区二区| 免费日韩欧美在线观看| 亚洲av国产av综合av卡| 欧美 日韩 精品 国产| 视频区图区小说| 日韩制服丝袜自拍偷拍| 亚洲精品国产一区二区精华液| 中文字幕人妻丝袜一区二区 | 老鸭窝网址在线观看| 高清视频免费观看一区二区| √禁漫天堂资源中文www| 久久久久精品人妻al黑| 美女福利国产在线| 国产成人精品在线电影| 97人妻天天添夜夜摸| 亚洲在久久综合| 七月丁香在线播放| 午夜福利一区二区在线看| 亚洲国产中文字幕在线视频| 欧美人与善性xxx| 一本久久精品| 老熟女久久久| 欧美少妇被猛烈插入视频| 十分钟在线观看高清视频www| 国产亚洲精品第一综合不卡| 老鸭窝网址在线观看| 美女视频免费永久观看网站| 国产97色在线日韩免费| 97精品久久久久久久久久精品| av网站免费在线观看视频| 国产一区二区 视频在线| 国产毛片在线视频| 国产精品女同一区二区软件| 成人免费观看视频高清| 99国产综合亚洲精品| 精品国产一区二区久久| 啦啦啦啦在线视频资源| 青草久久国产| 秋霞在线观看毛片| 国产xxxxx性猛交| 国产精品免费大片| 国产精品.久久久| 欧美日韩视频高清一区二区三区二| 国产在线视频一区二区| 国产精品久久久久成人av| 久久99一区二区三区| 老司机靠b影院| 久久久久视频综合| 一级爰片在线观看| 亚洲精品国产色婷婷电影| 久久毛片免费看一区二区三区| 久久99热这里只频精品6学生| 久久久久精品人妻al黑| 精品福利永久在线观看| 蜜桃国产av成人99| 麻豆av在线久日| 欧美成人精品欧美一级黄| 一区在线观看完整版| 亚洲图色成人| 国产亚洲最大av| 精品亚洲成a人片在线观看| 亚洲国产精品一区三区| 黄频高清免费视频| 国精品久久久久久国模美| 亚洲欧美精品自产自拍| 一本大道久久a久久精品| 日本av免费视频播放| 成年av动漫网址| 精品福利永久在线观看| 美女大奶头黄色视频| 最近2019中文字幕mv第一页| 精品久久久久久电影网| 纯流量卡能插随身wifi吗| 国产日韩欧美在线精品| 菩萨蛮人人尽说江南好唐韦庄| 日韩一区二区三区影片| 亚洲成国产人片在线观看| 免费日韩欧美在线观看| 日韩伦理黄色片| 日韩电影二区| 亚洲美女视频黄频| 宅男免费午夜| 如日韩欧美国产精品一区二区三区| 天天影视国产精品| 黑人欧美特级aaaaaa片| 十八禁高潮呻吟视频| 成人影院久久| 欧美精品高潮呻吟av久久| 国产日韩一区二区三区精品不卡| 99国产综合亚洲精品| 国语对白做爰xxxⅹ性视频网站| 亚洲精品日韩在线中文字幕| 亚洲成人免费av在线播放| 女人精品久久久久毛片| 亚洲在久久综合| a 毛片基地| 男女下面插进去视频免费观看| 黄色 视频免费看| av网站免费在线观看视频| 亚洲伊人色综图| 久久久久久久久免费视频了| 秋霞伦理黄片| 一区二区三区四区激情视频| 在线看a的网站| 亚洲精品日本国产第一区| 亚洲av成人精品一二三区| 麻豆乱淫一区二区| 日本一区二区免费在线视频| 亚洲人成网站在线观看播放| 中文字幕人妻丝袜制服| 久热爱精品视频在线9| 捣出白浆h1v1| 国产成人免费无遮挡视频| 欧美成人精品欧美一级黄| e午夜精品久久久久久久| 国产亚洲精品第一综合不卡| 在线观看免费午夜福利视频| 高清在线视频一区二区三区| 99热国产这里只有精品6| 欧美av亚洲av综合av国产av | 亚洲精品国产色婷婷电影| 国产精品久久久久久精品古装| 日本wwww免费看| 九色亚洲精品在线播放| 街头女战士在线观看网站| 国产精品免费视频内射| 成人毛片60女人毛片免费| 成人18禁高潮啪啪吃奶动态图| 赤兔流量卡办理| 久久久久国产一级毛片高清牌| 亚洲精品久久成人aⅴ小说| 久久久国产精品麻豆| 卡戴珊不雅视频在线播放| 黄色 视频免费看| 80岁老熟妇乱子伦牲交| 国产欧美日韩综合在线一区二区|