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

    頻率步進信號距離-方位聯(lián)合超分辨成像方法

    2016-11-29 03:50:57呂明久李少東楊軍馬曉巖
    電波科學(xué)學(xué)報 2016年3期
    關(guān)鍵詞:方位重構(gòu)距離

    呂明久 李少東 楊軍 馬曉巖

    (1.空軍預(yù)警學(xué)院研究生隊,武漢 430019;2.空軍預(yù)警學(xué)院,武漢 430019)

    ?

    頻率步進信號距離-方位聯(lián)合超分辨成像方法

    呂明久1李少東1楊軍2馬曉巖2

    (1.空軍預(yù)警學(xué)院研究生隊,武漢 430019;2.空軍預(yù)警學(xué)院,武漢 430019)

    對距離向、方位向二維稀疏的頻率步進信號模型,基于壓縮感知理論,研究了一種適用于稀疏頻率步進回波的距離-方位聯(lián)合逆合成孔徑雷達(Inverse Synthetic Aperture Radar,ISAR)超分辨成像方法.首先對稀疏頻率步進回波進行建模,在此基礎(chǔ)上構(gòu)建了距離向和方位向二維稀疏時的聯(lián)合稀疏基,最后利用二維平滑0-范數(shù)法在矩陣域直接進行處理,得到最終的ISAR超分辨成像結(jié)果.并對算法復(fù)雜度、超分辨性能進行了分析,得出了相應(yīng)的結(jié)論.理論分析和仿真結(jié)果表明所提方法在不同稀疏方式、不同稀疏條件下具有更好成像質(zhì)量、更快處理速度的優(yōu)勢.

    稀疏頻率步進信號;二維平滑0-范數(shù);壓縮感知;ISAR超分辨成像

    DOI 10.13443/j.cjors.2015083103

    引 言

    合成孔徑雷達/逆合成孔徑雷達(Synthetic Aperture Radar/Inverse Synthetic Aperture Radar, SAR/ISAR)成像系統(tǒng)主要通過增大信號帶寬達到高的距離分辨率,方位高分辨則主要通過增加方位積累時間獲得.頻率步進信號(Frequency-SteppedSignal,FSS)作為一種大帶寬信號,主要通過發(fā)射多個載頻跳變的子脈沖實現(xiàn)大的信號帶寬,且具有系統(tǒng)復(fù)雜性低、成本少、易于工程實現(xiàn)等優(yōu)點,已成為高分辨率雷達技術(shù)的發(fā)展趨勢[1-5].

    由于FSS需要發(fā)射多組脈沖串才能實現(xiàn)成像,每組脈沖串又包含多個子脈沖,當(dāng)目標(biāo)在距離、方位向存在脈沖缺失或遭受干擾部分脈組無效等情況時,回波將變?yōu)槎S稀疏頻率步進信號(Sparse Frequency-Stepped Signal,SFSS),這種距離、方位同時稀疏的SFSS給平動補償以及后續(xù)的寬帶合成帶來困難,使得常規(guī)成像算法失效.近年來,隨著壓縮感知(Compressive Sensing,CS)理論[6]被用于SAR/ISAR成像系統(tǒng)中,為稀疏情況下的目標(biāo)高分辨成像提供了一種思路[7-8].

    目前,稀疏條件下的高分辨成像方法主要有兩種研究思路:一是通過統(tǒng)計外推的形式,估計空缺的脈沖后再進行成像,如文獻[9]針對距離向稀疏的SFSS信號,通過自回歸(Auto Regressive,AR)模型、柯西高斯模型對缺損的子脈沖進行估計后進行成像,但是這種成像方法往往精度不高.第二種是利用目標(biāo)在觀測空間是稀疏的這一特征,結(jié)合近年來的稀疏表示理論或CS直接進行成像.其中一類是將距離向、方位向成像分別處理的距離-方位解耦算法[10-12],這類方法處理速度相對較快,能夠?qū)崿F(xiàn)高分辨的成像,但是由于沒有充分利用回波數(shù)據(jù)距離向、方位向的耦合信息,成像質(zhì)量將會有損失.另一類方法是將距離向、方位向同時進行處理的距離-方位聯(lián)合算法[12-14].距離-方位聯(lián)合處理的好處在于充分利用了距離-方位耦合信息,成像質(zhì)量得到提升.但目前方法大多是通過行列堆疊方式將矩陣運算轉(zhuǎn)換為一維向量形式后再利用CS重構(gòu)算法進行重構(gòu),這樣對硬件的存儲提出了較高的要求,運算量較大,有時甚至難以實現(xiàn).

    針對上述問題,本文基于CS理論,研究了一種適用于稀疏頻率步進回波的距離-方位聯(lián)合ISAR超分辨成像算法.首先對SFSS進行建模,給出了SFSS回波模型.其次對SFSS回波進行平動補償,將其轉(zhuǎn)化為轉(zhuǎn)臺成像模型.然后結(jié)合目標(biāo)回波模型和壓縮感知理論之間的關(guān)系,得到了距離向和方位向聯(lián)合稀疏基,最后采用二維平滑0-范數(shù)法(Two Dimensional Smoothed L0-norm,2D-SL0)在矩陣域直接對回波數(shù)據(jù)矩陣進行距離-方位聯(lián)合超分辨成像,避免了常規(guī)距離-方位聯(lián)合成像時的行列堆疊過程.理論分析和仿真結(jié)果表明:該方法不僅適用于完整回波條件下的FSS,對于不同稀疏程度、不同稀疏方式下的SFSS同樣適用.理論分析和仿真結(jié)果驗證了所提方法具有更好的成像性能和更快的處理速度.

    1 SFSS回波建模

    SFSS可以看作為傳統(tǒng)FSS的稀疏采樣形式.因此,首先對FSS回波進行建模,在此基礎(chǔ)上得出SFSS回波模型.一組FSS可以表示為[1]

    (1)

    式中: fn=f0+nΔf,n=0,1,2,…,N-1,N為子脈沖數(shù); rect(·)為窗函數(shù); f0為載頻; t為全時間; T為子脈寬; Tr為重復(fù)周期; Δf為步進頻率間隔,此時等效的信號帶寬為B=NΔf.

    假設(shè)目標(biāo)共有I個散射點,在短觀測時間內(nèi),目標(biāo)沿雷達徑向的速度為ν,徑向加速度為a,那么t時刻第i(i=1,…,I)個散射點與雷達的距離可表示為

    xicos[θ(tm)]-yisin[θ(tm)].

    (2)

    式中: R0為初始時刻參考距離; xi,yi分別為第i個散射點在目標(biāo)坐標(biāo)系上的橫縱坐標(biāo); θ(tm)為與慢時間tm有關(guān)的目標(biāo)轉(zhuǎn)動角度.此時距離雷達Ri(t,tm)的各散射點回波可以表示為

    exp[j2πfn(t-2Ri(t,tm)/c)].

    (3)

    式中, σi(xi,yi)表示第i個散射點強度.

    混頻后的基帶信號為

    exp[j2πfn(-2Ri(t,tm)/c)].

    (4)

    假設(shè)目標(biāo)在觀察時間內(nèi)勻速轉(zhuǎn)動,則θ(tm)=mΔθ,m=0,1,…,M-1,Δθ為角轉(zhuǎn)動步長,M為脈組數(shù).對于FSS信號,只要目標(biāo)速度滿足式(5)即可忽略包絡(luò)走動帶來的影響[1].

    v≤cT/4NTr.

    (5)

    (6)

    式中,

    (7)

    是由目標(biāo)運動而引入的平動相位項,必須進行補償.當(dāng)在觀測時間內(nèi)目標(biāo)轉(zhuǎn)動的角度較小時,則cos(θ(tm))≈1,sin(θ(tm))≈θ(tm).因此,式(6)可寫成

    (8)

    將式(8)中的常數(shù)項并入幅度信息中,得到

    exp[j2π(f0+nΔf)2yimΔθ/c].

    (9)

    由于信號在距離向與方位向采樣點數(shù)分別為N、M,所以距離向分辨率為Δx=c/(2NΔf),方位向分辨率為Δy=λ0/(2MΔθ),其中λ0=c/f0.

    按照距離、方位分辨率將目標(biāo)進行離散,得到如圖1所示的目標(biāo)區(qū)域離散化形式σ=[σpq]N×M.由于λn≈λ0,所以式(9)可寫為

    (10)

    圖1 雷達與目標(biāo)關(guān)系示意圖

    將式(10)寫成矩陣形式為

    (11)

    本文研究的SFSS就是在FSS的基礎(chǔ)上距離向、方位向二維同時稀疏,圖2為其發(fā)射波形的示意圖.

    (a) 脈組內(nèi)稀疏 (b) 二維稀疏情況圖2 二維稀疏步進頻信號示意圖

    圖2(a)為一組脈沖內(nèi)部子脈沖稀疏情況示意圖,圖中虛線脈沖代表稀疏的子脈沖.圖2(b)為信號距離、方位同時稀疏時回波信號的示意圖,圖中黑色部分代表接收到的有用脈沖數(shù)據(jù),無色部分代表稀疏的脈沖數(shù)據(jù).當(dāng)距離向、方位向同時稀疏時,則SFSS回波模型可以寫成

    (12)

    式中: U為稀疏回波數(shù)據(jù)矩陣; Dr為距離向稀疏矩陣,如式(13)所示的對角矩陣形式,對角線上第n個元素為0代表距離向第n個子脈沖缺失,

    (13)

    Da為方位向稀疏矩陣,如式(14)所示的對角矩陣形式,同理對角線上第m個元素為0代表距離向第m組回波缺失,

    (14)

    顯然,式(12)是個病態(tài)方程,存在無窮多個解,可以利用壓縮感知重構(gòu)理論實現(xiàn)求解.

    2 SFSS距離-方位聯(lián)合超分辨成像

    ISAR成像的前提是對目標(biāo)平動分量進行補償,使其轉(zhuǎn)化為等效的轉(zhuǎn)臺模型.假設(shè)目標(biāo)已經(jīng)經(jīng)過平動補償(具體補償過程可參考文獻[8,15]),即式(12)可以轉(zhuǎn)化為式(15)所示的轉(zhuǎn)臺模型:

    (15)

    根據(jù)CS理論,如果目標(biāo)是稀疏或可壓縮的,則可以通過CS理論準(zhǔn)確重構(gòu)出目標(biāo)圖像且觀測值也可顯著減少[6].實際上,當(dāng)分辨率確定時,目標(biāo)僅包含有限個散射點,即目標(biāo)具有稀疏性,可以用CS理論對目標(biāo)進行重構(gòu)實現(xiàn)最終成像.

    當(dāng)劃分的網(wǎng)格間隔變小時,可進一步提高分辨率[11].假設(shè)距離、方位網(wǎng)格間隔各提高J倍,得到的目標(biāo)離散表示為σ=[σpq]P×Q,其中P=NJ,Q=MJ,此時回波信號模型為

    (16)

    此時距離向、方位向分辨率分別為Δx′=Δx/J,Δy′=Δy/J,比原始分辨率提高了J倍.將式(16)寫成式(15)的形式,得到

    (17)

    (18)

    對于式(18)的模型,可以對距離向、方位向進行解耦,得到:

    U′=ArF;

    (19)

    FT=Aaσ′T.

    (20)

    式中, F為距離像信息.利用式(19)、(20)先進行距離向成像后進行方位向成像稱之為距離-方位解耦成像方法[11].可以通過行列堆疊的方式將式(19)、(20)變?yōu)橐痪S向量的形式分別進行求解[10],也可以通過SL0多稀疏恢復(fù)[16-17](SL0 for Multiple Sparse Recovery,SML0)算法直接對式(19)、(20)進行操作,這樣可以減少運算量,提高運算速度.

    實際上,目標(biāo)的回波數(shù)據(jù)通常是耦合在一起的,將距離向、方位向分開進行重構(gòu)沒有充分利用距離、方位的耦合信息,成像質(zhì)量將會有所損失[11].因此,為進一步提高運算速度和成像質(zhì)量,可以直接對式(18)進行操作,即距離-方位聯(lián)合成像方法.利用CS解決此類模型時常用行列堆疊的方式將二維矩陣運算轉(zhuǎn)化為一維向量形式,再通過正交匹配追蹤(Orthogonal Matching Pursuit,OMP)、平滑0-范數(shù)[16](Smoothed L0-norm,SL0)法等稀疏重構(gòu)算法進行成像.SFSS波形行列堆疊形式可以表示為

    U=Θδ.

    (21)

    文獻[18-19]針對式(18)中矢量的0-范數(shù)最小化問題不是連續(xù)函數(shù),常規(guī)算法無法直接求解,提出利用高斯函數(shù)來近似0-范數(shù),將式(18)的NP問題轉(zhuǎn)化為求解目標(biāo)函數(shù)的最大值問題,即2D-SL0.該算法直接用式(18)的矩陣模型取代式(21)的行列堆疊模型進行求解,避免了行列堆疊后的大維度矩陣帶來的處理困難問題,減小了對處理器內(nèi)存和運算能力的要求,提高了成像效率.表1給出了2D-SL0算法主要處理步驟[18-19].

    表1 2D-SL0算法主要處理步驟

    3 算法性能分析

    3.1 內(nèi)存需求及運算量分析

    本文研究的距離-方位聯(lián)合超分辨成像方法主要基于2D-SL0算法在矩陣域直接進行成像.下面對本文算法的運算量進行分析.為便于比較性能優(yōu)劣,同時給出基于匹配追蹤類算法[14]進行距離-方位聯(lián)合重構(gòu)以及利用匹配追蹤類算法[10]、SL0算法[12]、MSL0[16]算法進行距離-方位解耦重構(gòu)的運算量.下文中匹配追蹤類算法均基于OMP算法進行分析.

    從內(nèi)存存儲需求來看:

    基于距離-方位聯(lián)合模型,本文方法直接對矩陣進行處理,存儲Ar與Aa所需的存儲空間大小為C1=O(NP)+O(MQ);文獻[14]通過行列堆疊后利用匹配追蹤類算法進行處理,存儲Θ所需要的存儲空間大小為C2=O(NPMQ).

    基于距離-方位解耦模型,文獻[10]通過行列堆疊后利用匹配追蹤類算法分別對距離向、方位向分開處理,所需的存儲空間大小為C3=O(NP)+O(MQ); 文獻[12]通過SL0算法進行成像處理,存儲空間同樣為C4=O(NP)+O(MQ);利用文獻[16]的MSL0算法對距離向、方位向進行處理,所需的存儲空間大小同樣為C5=O(NP)+O(MQ).

    下面再來分析各個算法的計算量:

    基于距離-方位聯(lián)合模型,本文方法主要計算量為矩陣求偽逆過程以及投影步驟.求偽逆步驟的計算量為O(NP)+O(MQ),投影步驟的計算量為O((N+M)PQ),總計算量為T1≈O((N+M)PQ);文獻[16]方法主要計算量集中在求內(nèi)積與矩陣求偽逆過程,且與算法的稀疏度設(shè)置有關(guān),總的計算量可以表示為T2=O(K1MNPQ),K1為設(shè)置的稀疏度大?。?/p>

    基于距離-方位解耦成像模型,文獻[10]方法總的計算量為T3=O(K2NPQ)+O(K3MPQ),K2、K3分別為距離向、方位向重構(gòu)時設(shè)置的稀疏度大?。墨I[12]通過SL0算法需要的總計算量約為T4=O(MP2)+O(PQ2).文獻[16]利用MSL0算法的總計算量約為T5=O(MP1.376)+O(PQ1.376).

    通過上述分析可知:

    在內(nèi)存需求上,C1=C3=C4=C5

    在計算量上,T1

    因此,本文利用2D-SL0算法進行SFSS-ISAR距離-方位聯(lián)合成像具有降低對硬件的要求以及減少重構(gòu)計算量的優(yōu)勢.當(dāng)信號子脈沖數(shù)、脈組數(shù)變大時,雖然可以提高雷達的分辨力,但是帶來的卻是算法的內(nèi)存需求、處理運算量明顯增加.在實際應(yīng)用中應(yīng)根據(jù)實際需求選擇合適的發(fā)射信號參數(shù),既可以滿足分辨率的需求,又可以減少處理時間,增強算法時效性.

    3.2 超分辨能力分析

    從式(16)可以看出,當(dāng)J>1時可以提高分辨率.因此CS的引入,可以進一步提高目標(biāo)成像的分辨率.研究表明:量測矩陣各列之間的相關(guān)性對最終重構(gòu)結(jié)果有較大影響.當(dāng)分辨率一定時,量測矩陣各列之間越不相關(guān),成像分辨率就越高[20].文獻[21]利用量測矩陣相關(guān)度

    (22)

    作為評價量測矩陣性能優(yōu)劣的指標(biāo),其中ρ為Θ的列向量.μ(Θ)值越小,說明量測矩陣性能越好.因此在進行高分辨成像時,不同的稀疏方式對成像性能影響很大.量測矩陣相關(guān)度越小,每次量測的信息交疊、干擾就越小,分辨效果便會越好.所以選擇相干性低的量測矩陣更有利于分辨能力的提高.另外,需要指出的是高分辨倍數(shù)不可能隨意增大,當(dāng)倍數(shù)增大后矩陣的維數(shù)和各列之間的相干性也會隨之增大,既增加了處理時間又會降低重構(gòu)概率,關(guān)于最優(yōu)的超分辨倍數(shù)研究將是下一步研究的重點.

    4 仿真分析與驗證

    一般ISAR成像雷達系統(tǒng)多工作于X波段[12,23],因此設(shè)置典型X波段雷達參數(shù)進行仿真研究: FSS共發(fā)射M=255組步進頻率信號,每組信號包含N=256個子脈沖,其載頻f0=10GHz,脈寬T=5μs,頻率步進量Δf=2MHz,每組信號重復(fù)頻率fPR=200Hz,參考距離R0=50.5km,J=3,K2=K3=50,信噪比RSN=10dB,添加噪聲方式同文獻[22].假設(shè)目標(biāo)模型如圖3所示,目標(biāo)初始飛行方向與X軸夾角為60°,飛行速度v=400m/s.

    根據(jù)上述參數(shù),文獻[14]研究的距離-方位聯(lián)合成像方法要求的內(nèi)存存儲容量已經(jīng)大大超出所用計算機硬件處理能力(所用計算機配置為:處理器為Intel酷睿E7500,主頻2.93GHz,內(nèi)存2GB).因此,下文分別對本文方法、文獻[10]方法(方法1)、文獻[12]方法(方法2)和文獻[16]方法(方法3)進行比較研究.

    圖3 仿真目標(biāo)模型

    仿真1 全數(shù)據(jù)條件下的成像效果仿真

    首先,當(dāng)距離向、方位向不稀疏時,即全部數(shù)據(jù)參與最終的成像處理,仿真結(jié)果如圖4所示.

    圖4(a)為距離-方位解耦OMP算法進行處理的結(jié)果,由于該方法對稀疏度大小較為敏感,因此,重構(gòu)的結(jié)果存在較多虛假散射點; 圖4(b)~(c)為距離-方位解耦SL0算法、MSL0算法進行成像的結(jié)果.但是由于沒有利用距離-方位耦合信息,因此成像結(jié)果有所損失.圖4(d)為本文方法進行距離-方位聯(lián)合處理的結(jié)果.由于本文方法利用了距離向、方位向的耦合信息,所以成像效果在四種方法中最好.

    (a) 方法1成像結(jié)果

    (b) 方法2成像結(jié)果

    (c) 方法3成像結(jié)果

    (d) 本文方法成像結(jié)果圖4 完整回波條件下運動目標(biāo)成像仿真

    仿真2 稀疏條件下的成像效果仿真

    為驗證本文方法對不同稀疏度條件下的SFSS成像效果,首先定義SFSS稀疏度為K,其計算方式為K=N′/N且K=M′/M,N′、M′分別為缺失的子脈沖數(shù)與脈組數(shù).本節(jié)主要研究以下兩種稀疏方式下的SFSS成像:

    方式1:距離向、方位向稀疏形式為各自隨機稀疏,即距離向、方位向各自按照稀疏度大小進行隨機稀疏.

    方式2:距離向、方位向稀疏形式為各自分塊稀疏,即在距離向、方位向的某個位置集中稀疏.

    圖5給出了K=0.25時兩種稀疏方式條件下的稀疏回波數(shù)據(jù)示意圖.

    針對圖5中的兩種稀疏方式,按照式(22)可量測矩陣稀疏度分別為μ1=0.056 5、μ2=0.225 1,可見方式1的量測矩陣相關(guān)性明顯小于方式2.因此稀疏方式1的重構(gòu)效果要好于稀疏方式2的重構(gòu)效果.下面以稀疏方式1條件下的SFSS為例,圖6為K=0.25時幾種不同算法的最終成像結(jié)果.

    (a) 稀疏方式1

    (b) 稀疏方式2 圖5 不同隨機方式下的回波數(shù)據(jù)示意圖

    (a) 方法1成像結(jié)果

    (b) 方法2成像結(jié)果

    (c) 方法3成像結(jié)果

    (d) 本文方法成像結(jié)果圖6 方式1稀疏時的成像結(jié)果

    從圖6可以看出:在方式1的隨機稀疏方式下且K=0.25時,幾種算法都能實現(xiàn)最終的成像.從成像效果來看,圖6(a)~(c)的成像結(jié)果具有較多的旁瓣,尤其是機頭兩個散射點分辨不清.另外,圖6(a)的成像結(jié)果存在大量的虛假重構(gòu)點,效果最差.而本文方法得到的成像結(jié)果凝聚性最好,旁瓣也最少,機頭的兩個散射點能夠清晰地分辨出來,因此性能也最優(yōu).

    仿真3SFSS成像性能效果仿真對比

    為驗證本文方法在不同稀疏方式、不同稀疏度條件下的成像性能,采用圖像熵(Entropy)以及目標(biāo)背景比(Target-to-BackgroundRatio,TBR)[23]作為成像質(zhì)量好壞的評價指標(biāo)對成像質(zhì)量進行仿真對比.成像質(zhì)量越好,熵值則越小,TBR值則越大.其中在計算TBR值時利用全數(shù)據(jù)條件下信噪比為20dB時傳統(tǒng)RD算法成像結(jié)果作為背景支撐區(qū).圖7為不同稀疏方式及稀疏度下的成像性能曲線.

    從圖7(a)可以看出:隨著稀疏度的增加,幾種方法的熵值逐漸增大,但是本文距離-方位聯(lián)合成像的熵值始終最小,說明成像質(zhì)量最好.對于不同的稀疏方式,稀疏方式1的熵值小于稀疏方式2的熵值,說明了在稀疏方式1的條件下,成像性能更好.對于圖7(b)中的TBR值,同樣可以看出,本文方法的成像結(jié)果TBR值始終大于其他算法,且稀疏方式1的TBR值要大于稀疏方式2的TBR值,說明方式1的稀疏方式重構(gòu)性能更好,這與量測矩陣相關(guān)度分析結(jié)果也相吻合.

    為進一步驗證本文方法在處理時間上的優(yōu)勢,表2給出了不同方法的成像時間對比.

    (a) 熵值Er大小對比圖

    (b) TBR值大小對比圖圖7 SFSS成像性能對比

    稀疏度稀疏方式1方法1方法2方法3本文方法0.1541972386978556600254184538389788746稀疏度稀疏方式2方法1方法2方法3本文方法01541187378397346990253818936868741659

    從表2的結(jié)果可以看出: 在不同的稀疏方式、稀疏度條件下,本文方法與方法3相比于其他兩種算法處理時間明顯減少; 相比于方法3,本文方法的處理時間進一步減少.這也驗證了本文方法在處理時間上的優(yōu)勢.

    5 結(jié) 論

    本文針對SFSS-ISAR超分辨成像問題提出了一種距離-方位聯(lián)合超分辨成像方法,該方法具有更好的成像質(zhì)量,成像處理速度也更快,并從理論和仿真兩方面進行了分析證明.當(dāng)信號稀疏時,傳統(tǒng)的運動補償方法補償難度將增大,這會對后續(xù)的成像造成不利影響,甚至無法成像,值得進一步進行研究.另外,文中研究了兩種隨機稀疏的方式,對于其他不同稀疏方式條件下的成像效果與成像方法也將作為下一步研究的一個重要內(nèi)容.

    [1] 毛二可, 龍騰, 韓月秋. 頻率步進雷達數(shù)字信號處理[J]. 航空學(xué)報,2001, 22(增刊): 16-25.

    MAO E K, LONG T, HAN Y Q. Digital signal processing of stepped frequency radar[J]. Acta aeronautica et astronautica sinica, 2001, 22(Sup): 16-25. (in Chinese)

    [2] SEYFRIED D, SCHOEBEL J. Stepped-frequency radar signal processing[J]. Journal of applied geophysics, 2015, 1(112): 42-51.

    [3] LEE. S H, KANG. M S, PARK. S H, et al. Inverse synthetic aperture radar imaging using stepped chirp waveform[J]. The journal of korean institute of electromagnetic engineering and science, 2014, 25(9): 930-937.

    [4] TANG L, ZHU Y F, FU Q. Stepped frequency SAR imaging based on compensation in Doppler domain[C]. Signal Processing, Communications and Computing: ICSPCC 2014: IEEE International Conference on Signal Processing, Communications and Computing, Guilin, August 5-8, 2014[C]. IEEE, 2014.

    [5] 黃瓊, 屈樂樂, 吳秉橫, 等. 壓縮感知在超寬帶雷達成像中的應(yīng)用[J]. 電波科學(xué)學(xué)報, 2010, 25(1): 77-82+202.

    HUANG Q, QU L L, WU B H, et al. Compressive sensing for ultra-wideband radar imaging[J]. Chinese journal of radio science, 2010, 25(1): 77-82+202. (in Chinese)

    [6] DONOHO D L. Compressed sensing[J]. IEEE transations on information theory, 2006, 52(4): 1289-1306.

    [7] 黃大榮, 郭新榮, 張磊, 等. 稀疏孔徑ISAR機動目標(biāo)成像與相位補償方法[J]. 航空學(xué)報, 2014, 35(7): 2019-2030.

    HUANG D R, GUO X R, ZHANG L, et al. ISAR phase compensation and imaging of maneuvering target with sparse apertures[J]. Acta aeronautica et astronautica sinica, 2014, 35(7): 2019-2030. (in Chinese)[8] 顧福飛, 張群, 梁穎, 等. 基于壓縮感知的地面運動目標(biāo)成像方法[J]. 電波科學(xué)學(xué)報, 2013, 28(4): 764-770.

    GU F F, ZHANG Q, LIANG Y, et al. Imaging method of ground moving targets based on compressed sensing[J]. Chinese journal of radio science, 2013 28(4): 764-770. (in Chinese)

    [9] LIU Q H, LIU H B, LONG T. Estimation of stepped frequency radar missing echo pulses[C]//Radar Conference, IET international Radar Conference. Guilin. April 20-22, 2009.

    [10] LIU J H, XU S K, GAO X Z, et al. Novel imaging methods of stepped frequency radar based on compressed sensing[J]. Journal of systems engineering and electronics, 2012, 23(1): 47-56.

    [11] 杜小勇,胡衛(wèi)東,郁文賢. 基于稀疏成份分析的逆合成孔徑雷達成像技術(shù)[J]. 電子學(xué)報, 2006, 34(3): 491-495.

    DU X Y, HU W D, YU W X. ISAR imaging techniques based on sparse component analysis[J]. ACTA electronica sinica, 2006, 34(3): 491-495. (in Chinese)

    [12]劉記紅. 基于壓縮感知的ISAR成像技術(shù)研究[D]. 湖南:國防科學(xué)技術(shù)大學(xué), 2012.

    LIU J H. Inverse synthetic aperture radar imaging technique based on compressed sensing[D]. Hunan: National University of Defense Technology, 2012. (in Chinese)

    [13]吳敏, 邢孟道, 張磊. 基于壓縮感知的二維聯(lián)合超分辨ISAR成像算法[J]. 電子與信息學(xué)報, 2014(1): 187-193.

    WU M, XING M D, ZHANG L. Two dimensional joint super-resolution ISAR imaging algorithm based on compressive sensing[J]. Journal of electronics & information technology, 2014(1): 187-193. (in Chinese)

    [14]YANG J G, THOMPSON J, HUANG X T, et al. Random-frequency SAR imaging based on compressed sensing[J]. IEEE transactions on geoscience and remote sensing, 2013, 51(2): 983-994.

    [15]KIM K T. Focusing of high range resolution profiles of moving targets using stepped frequency waveforms[J]. IET radar sonar navigation, 2010, 4(4): 564-575.

    [16]MOHIMANI H, BABAIE Z M, GORODNIITSY I, et al. Sparse recovery using smoothed L0 (SL0): convergence analysis[J/OL]. (2010-1-28)[2015-08-15]. http://120.52.73.76/arxiv.org/pdf/1001.5073v1.pdf.

    [17]李少東, 陳文峰, 楊軍, 等. 任意稀疏結(jié)構(gòu)的多量測向量快速稀疏重構(gòu)算法研究[J]. 電子學(xué)報, 2015, 43(4): 708-715.

    LI S D, CHEN W F, YANG J, et al. Study on the fast sparse recovery algorithm via multiple measurement vectors of arbitrary sparse structure[J]. ACTA electronica sinica, 2015, 43(4): 708-715. (in Chinese)

    [18]GHAFFARI A, BABAIE Z M, JUTTEN C. Sparse decomposition of two dimensional signals[C]//Proceeding IEEE International Conference on Acoustics, Speech and Signal Process, 2009: 3157-3160.

    [19]EFTEKHARI A, BABAIE M, MOGHADDAM H A. Two-dimensional random projection[J]. Signal processing, 2011, 91(7): 1589-1603.

    [20]李龍珍, 姚旭日, 劉雪峰, 等. 基于壓縮感知超分辨鬼成像[J]. 物理學(xué)報. 2014, 22(63): 224201(1-7).

    LI L Z, YAO X R, LIU X F, et al. Super-resolution ghost imaging via compressed sensing[J]. Acta physica sinica, 2014, 22(63): 224201(1-7). (in Chinese)

    [21]ALBERT F J, LIAO W J. Coherence pattern-guided compressive sensing with unresolved grids[J]. Siam journal on imaging sciences, 2012, 5177(1): 179-202.

    [22]YE W, YEO T S, BAO Z. Weighted least squares estimation of phase errors for SAR/ISAR autofocus[J]. IEEE transactions on geoscience & remote sensing, 1999, 37(5): 2487-2494.

    [23]張磊. 高分辨SAR/ISAR成像及誤差補償技術(shù)研究[D]. 西安:西安電子科技大學(xué),2012.

    ZHANG L. Study on high resolution SAR/ISAR imaging and error correction[D]. Xi’an: University of Electronics Science and Technology of China, 2012. (in Chinese)

    呂明久 (1985-),男,安徽人,空軍預(yù)警學(xué)院博士研究生,研究方向為信號檢測與識別.

    Range-azimuth joint ISAR super-resolution imaging method based on frequency-stepped signal

    Lü Mingjiu1LI Shaodong1YANG Jun2MA Xiaoyan2

    (1.Department of Graduate Management, Air Force Early Warning Academy, Wuhan 430019, China;2.AirForceEarlyWarningAcademy,Wuhan430019,China)

    For the sparse frequency-stepped (SFS) radar signal which is sparse both in range and azimuth, a novel range-azimuth joint ISAR Imaging method is proposed based on compressive sensing (CS) theory. Firstly, the echo model of the SFSS is established. Secondly, the range-azimuth joint sparse base is built, which is sparse both in the range and azimuth dimension. Finally, the super-resolution ISAR Imaging can be realized via the 2D-Smoothed L0 algorithm which can dispose this echo model in matrix domain directly. Some performances of the proposed algorithm such as complexity and the super resolution capability are also analyzed and some corresponding results are given. At last, theoretical analysis and simulation experiments verify that the proposed algorithm have a better imaging performance and a faster processing speed at different sparse conditions and sparse patterns.

    sparse frequency-stepped signal; two dimensional smoothed L0; compressive sensing; ISAR super-resolution imaging

    10.13443/j.cjors. 2015083103

    2015-08-31

    TN957

    A

    1005-0388(2016)03-0537-09

    呂明久, 李少東, 楊軍, 等. 頻率步進信號距離-方位聯(lián)合超分辨成像方法[J],電波科學(xué)學(xué)報,2016,31(3):537-545.

    Lü M J, LI S D, YANG J, et al. Range-azimuth joint ISAR super-resolution imaging method based on frequency-stepped signal[J]. Chinese journal of radio science,2016,31(3):537-545. (in Chinese). DOI: 10.13443/j.cjors.2015083103

    聯(lián)系人: 呂明久 E-mail: lv_mingjiu@163.com

    猜你喜歡
    方位重構(gòu)距離
    長城敘事的重構(gòu)
    攝影世界(2022年1期)2022-01-21 10:50:14
    認(rèn)方位
    幼兒園(2021年12期)2021-11-06 05:10:20
    北方大陸 重構(gòu)未來
    算距離
    北京的重構(gòu)與再造
    商周刊(2017年6期)2017-08-22 03:42:36
    論中止行為及其對中止犯的重構(gòu)
    借助方位法的拆字
    中國修辭(2016年0期)2016-03-20 05:54:32
    說方位
    幼兒100(2016年28期)2016-02-28 21:26:17
    每次失敗都會距離成功更近一步
    山東青年(2016年3期)2016-02-28 14:25:55
    基于TMS320C6678的SAR方位向預(yù)濾波器的并行實現(xiàn)
    我要看日韩黄色一级片| 欧美激情国产日韩精品一区| 国产精品久久久久成人av| 亚洲内射少妇av| 纵有疾风起免费观看全集完整版| 国产成人精品一,二区| 插阴视频在线观看视频| 高清av免费在线| 亚洲美女视频黄频| 97超碰精品成人国产| 成人二区视频| 我要看日韩黄色一级片| 婷婷色麻豆天堂久久| 看免费成人av毛片| 伊人久久国产一区二区| 97精品久久久久久久久久精品| 亚洲国产精品一区三区| 日韩人妻高清精品专区| 国产精品国产av在线观看| 免费看日本二区| 国国产精品蜜臀av免费| 国产乱人视频| 97超碰精品成人国产| 91精品国产九色| 在线观看免费高清a一片| 深夜a级毛片| 色综合色国产| 97在线人人人人妻| 大片免费播放器 马上看| 国产一区亚洲一区在线观看| 日韩欧美 国产精品| 日韩在线高清观看一区二区三区| 国产在线一区二区三区精| 王馨瑶露胸无遮挡在线观看| www.色视频.com| 日产精品乱码卡一卡2卡三| 日韩,欧美,国产一区二区三区| 热re99久久精品国产66热6| 日韩免费高清中文字幕av| 欧美一区二区亚洲| 精品人妻视频免费看| 日韩国内少妇激情av| 亚洲精品一二三| 亚洲精品一区蜜桃| 秋霞在线观看毛片| 免费观看无遮挡的男女| 美女国产视频在线观看| 久久久亚洲精品成人影院| 内射极品少妇av片p| 99热这里只有是精品在线观看| 视频中文字幕在线观看| 亚洲综合色惰| 蜜臀久久99精品久久宅男| 尾随美女入室| 精品视频人人做人人爽| 香蕉精品网在线| 下体分泌物呈黄色| av卡一久久| 久久av网站| 亚洲精品国产av成人精品| 七月丁香在线播放| 国产有黄有色有爽视频| 熟女av电影| 久久久久久人妻| 黄色视频在线播放观看不卡| 黄色怎么调成土黄色| 日韩精品有码人妻一区| 亚洲久久久国产精品| 高清毛片免费看| 成年免费大片在线观看| 午夜福利影视在线免费观看| 国产av国产精品国产| 91久久精品国产一区二区三区| 欧美日韩视频精品一区| 一级毛片黄色毛片免费观看视频| 国产一区二区三区综合在线观看 | 日韩不卡一区二区三区视频在线| 亚洲成色77777| 国产成人a区在线观看| 嘟嘟电影网在线观看| av不卡在线播放| 在线观看免费高清a一片| 久久人人爽人人爽人人片va| 日韩伦理黄色片| 中文字幕av成人在线电影| 国产av精品麻豆| 亚洲中文av在线| 亚洲国产日韩一区二区| 免费看光身美女| 日韩在线高清观看一区二区三区| 99久久精品国产国产毛片| 亚州av有码| 亚洲欧美日韩无卡精品| 国产高清不卡午夜福利| 人人妻人人看人人澡| 最近最新中文字幕大全电影3| 亚洲成人一二三区av| 日本一二三区视频观看| 91精品国产国语对白视频| 美女xxoo啪啪120秒动态图| 一边亲一边摸免费视频| 亚洲自偷自拍三级| 国产精品人妻久久久影院| 日韩大片免费观看网站| 国产精品福利在线免费观看| 男女啪啪激烈高潮av片| 日日啪夜夜撸| 永久网站在线| 伦理电影大哥的女人| 亚洲国产日韩一区二区| 黄色日韩在线| 永久网站在线| 91午夜精品亚洲一区二区三区| 交换朋友夫妻互换小说| 亚洲色图av天堂| 国产男女超爽视频在线观看| 美女脱内裤让男人舔精品视频| 欧美日韩综合久久久久久| 久久国产精品大桥未久av | 久久国产精品男人的天堂亚洲 | 建设人人有责人人尽责人人享有的 | 日本-黄色视频高清免费观看| 看十八女毛片水多多多| 色5月婷婷丁香| 看免费成人av毛片| 女人久久www免费人成看片| 国产精品一区二区三区四区免费观看| 色综合色国产| 国产男女超爽视频在线观看| 夫妻性生交免费视频一级片| 边亲边吃奶的免费视频| 男女边摸边吃奶| 日本色播在线视频| 99久久中文字幕三级久久日本| 哪个播放器可以免费观看大片| 亚洲激情五月婷婷啪啪| 久久久久久久久久人人人人人人| 成人高潮视频无遮挡免费网站| 五月开心婷婷网| 国产av码专区亚洲av| 成年美女黄网站色视频大全免费 | 国精品久久久久久国模美| 日本av手机在线免费观看| 久久97久久精品| 大香蕉97超碰在线| 最近2019中文字幕mv第一页| 18禁动态无遮挡网站| 欧美日韩亚洲高清精品| 涩涩av久久男人的天堂| 免费人成在线观看视频色| 午夜视频国产福利| 联通29元200g的流量卡| 伦精品一区二区三区| 国产男女超爽视频在线观看| 亚洲精品aⅴ在线观看| 3wmmmm亚洲av在线观看| 国产成人91sexporn| 精品久久国产蜜桃| 综合色丁香网| 午夜老司机福利剧场| 国产精品熟女久久久久浪| 国产亚洲一区二区精品| 国产亚洲一区二区精品| 国产深夜福利视频在线观看| 国产 一区 欧美 日韩| 久久精品熟女亚洲av麻豆精品| 人妻 亚洲 视频| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 这个男人来自地球电影免费观看 | 26uuu在线亚洲综合色| 亚洲精品乱久久久久久| 亚洲精品国产成人久久av| 欧美日韩视频高清一区二区三区二| 亚洲va在线va天堂va国产| 免费看av在线观看网站| 少妇精品久久久久久久| 国产精品不卡视频一区二区| 黄色日韩在线| 男女啪啪激烈高潮av片| 最黄视频免费看| 亚洲内射少妇av| 免费大片18禁| 成人毛片60女人毛片免费| 大香蕉久久网| 高清午夜精品一区二区三区| 亚洲真实伦在线观看| 免费不卡的大黄色大毛片视频在线观看| 国产成人精品婷婷| 女人十人毛片免费观看3o分钟| 在现免费观看毛片| 97精品久久久久久久久久精品| 久久久久久久国产电影| 尾随美女入室| 美女脱内裤让男人舔精品视频| 一个人看的www免费观看视频| 伊人久久国产一区二区| 一本一本综合久久| 国产黄频视频在线观看| 国产精品久久久久久精品电影小说 | 国产精品国产av在线观看| 精品一品国产午夜福利视频| 高清视频免费观看一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲美女视频黄频| 嫩草影院新地址| 狂野欧美白嫩少妇大欣赏| 国产大屁股一区二区在线视频| 久久久久国产网址| 国产永久视频网站| 色吧在线观看| 久久99热这里只频精品6学生| 大香蕉97超碰在线| 亚洲精品国产av蜜桃| 成人午夜精彩视频在线观看| 国产av精品麻豆| 精品人妻视频免费看| 国产精品一及| 七月丁香在线播放| 一级毛片久久久久久久久女| 亚洲国产色片| 极品少妇高潮喷水抽搐| 日日摸夜夜添夜夜添av毛片| 免费观看在线日韩| 精品国产乱码久久久久久小说| 久久ye,这里只有精品| 五月天丁香电影| 狠狠精品人妻久久久久久综合| 80岁老熟妇乱子伦牲交| 99视频精品全部免费 在线| 91aial.com中文字幕在线观看| 日韩精品有码人妻一区| 久久影院123| 最近2019中文字幕mv第一页| 高清在线视频一区二区三区| 久久久久久久久久成人| 婷婷色麻豆天堂久久| 少妇的逼好多水| 高清欧美精品videossex| 亚洲精品乱码久久久v下载方式| 国产一区二区三区综合在线观看 | 国产高潮美女av| 日本vs欧美在线观看视频 | 午夜福利视频精品| 晚上一个人看的免费电影| 欧美日本视频| 国产高潮美女av| 日韩中字成人| 欧美国产精品一级二级三级 | 国产精品熟女久久久久浪| 最近最新中文字幕大全电影3| 日本黄色日本黄色录像| 汤姆久久久久久久影院中文字幕| 啦啦啦在线观看免费高清www| 日韩欧美精品免费久久| 精品亚洲乱码少妇综合久久| 国产成人一区二区在线| 日韩欧美 国产精品| 美女国产视频在线观看| 日本av手机在线免费观看| 免费黄色在线免费观看| 免费观看性生交大片5| 欧美另类一区| 麻豆成人午夜福利视频| 国产一区有黄有色的免费视频| 人人妻人人澡人人爽人人夜夜| 国产午夜精品久久久久久一区二区三区| 高清av免费在线| 夜夜骑夜夜射夜夜干| 亚洲成人一二三区av| 免费人成在线观看视频色| 高清午夜精品一区二区三区| 久久99热这里只有精品18| 91在线精品国自产拍蜜月| 韩国av在线不卡| 少妇的逼水好多| 永久免费av网站大全| 久久99热这里只有精品18| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 国产精品99久久久久久久久| 精品人妻偷拍中文字幕| 久久久久精品久久久久真实原创| 在线观看国产h片| 亚洲国产色片| 又黄又爽又刺激的免费视频.| 精品国产露脸久久av麻豆| av免费在线看不卡| 亚洲欧美日韩卡通动漫| 久久精品久久久久久久性| 欧美高清成人免费视频www| 天堂8中文在线网| 婷婷色综合www| 尤物成人国产欧美一区二区三区| 国产成人freesex在线| 免费观看a级毛片全部| 久久99热这里只频精品6学生| 日本黄色片子视频| 一个人免费看片子| 99精国产麻豆久久婷婷| 亚洲精品视频女| 五月伊人婷婷丁香| 深爱激情五月婷婷| 久久久久久久精品精品| 少妇人妻精品综合一区二区| 少妇精品久久久久久久| 全区人妻精品视频| 国产精品一及| 老司机影院成人| 夜夜爽夜夜爽视频| 1000部很黄的大片| 中文精品一卡2卡3卡4更新| 一区二区三区乱码不卡18| 老熟女久久久| 女性生殖器流出的白浆| 黄片wwwwww| 午夜免费鲁丝| 亚洲精品中文字幕在线视频 | 女性被躁到高潮视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲av二区三区四区| 黄片wwwwww| av国产久精品久网站免费入址| 夜夜爽夜夜爽视频| 在线观看av片永久免费下载| 中文精品一卡2卡3卡4更新| 2021少妇久久久久久久久久久| 观看免费一级毛片| 大陆偷拍与自拍| 啦啦啦在线观看免费高清www| 男人爽女人下面视频在线观看| 精品一品国产午夜福利视频| 亚洲精品乱码久久久久久按摩| 我的女老师完整版在线观看| 在现免费观看毛片| 精品一区二区三卡| 永久免费av网站大全| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 91精品一卡2卡3卡4卡| 日本av免费视频播放| 国产精品久久久久成人av| 伊人久久国产一区二区| 国产精品久久久久久久电影| 国产男人的电影天堂91| 国产黄色视频一区二区在线观看| 人妻一区二区av| 亚洲怡红院男人天堂| 99久国产av精品国产电影| 成年女人在线观看亚洲视频| 亚洲一区二区三区欧美精品| 一级毛片电影观看| 国产视频内射| 精品酒店卫生间| 99热国产这里只有精品6| 国模一区二区三区四区视频| 99热这里只有是精品50| 欧美精品一区二区大全| 欧美国产精品一级二级三级 | 成人无遮挡网站| 一个人免费看片子| 亚洲欧美精品自产自拍| 有码 亚洲区| 一级av片app| 国产在线免费精品| 日本wwww免费看| 能在线免费看毛片的网站| 在线观看人妻少妇| 日韩人妻高清精品专区| 中国国产av一级| 亚洲精品乱码久久久v下载方式| 日本免费在线观看一区| 又粗又硬又长又爽又黄的视频| 色视频www国产| 一个人看的www免费观看视频| 2018国产大陆天天弄谢| 国产亚洲一区二区精品| 亚洲色图av天堂| 在线观看三级黄色| 青春草国产在线视频| 国产精品.久久久| av网站免费在线观看视频| 国产永久视频网站| 搡老乐熟女国产| 熟妇人妻不卡中文字幕| 亚洲欧美一区二区三区黑人 | 亚洲人成网站在线观看播放| 一区二区三区精品91| 王馨瑶露胸无遮挡在线观看| 国产淫语在线视频| 久久99精品国语久久久| 综合色丁香网| 又爽又黄a免费视频| 日韩人妻高清精品专区| 91精品国产国语对白视频| 少妇精品久久久久久久| 久久热精品热| 综合色丁香网| 中文字幕人妻熟人妻熟丝袜美| av在线观看视频网站免费| 久久6这里有精品| videos熟女内射| 亚洲精品aⅴ在线观看| 欧美一区二区亚洲| 99久久人妻综合| 中文天堂在线官网| 啦啦啦视频在线资源免费观看| 亚洲第一av免费看| 亚洲熟女精品中文字幕| 能在线免费看毛片的网站| 黄色欧美视频在线观看| 国产爱豆传媒在线观看| 久久99热这里只频精品6学生| 老女人水多毛片| 亚洲精品久久久久久婷婷小说| 王馨瑶露胸无遮挡在线观看| 久久久精品免费免费高清| 免费播放大片免费观看视频在线观看| 国产精品久久久久久精品电影小说 | 一级a做视频免费观看| 亚洲精品成人av观看孕妇| 黄色一级大片看看| 精品午夜福利在线看| 青春草国产在线视频| 久久久久久久久久人人人人人人| 亚洲国产欧美人成| 亚洲精品国产av蜜桃| 永久网站在线| 深爱激情五月婷婷| 国产成人a区在线观看| 街头女战士在线观看网站| 色婷婷久久久亚洲欧美| 成人高潮视频无遮挡免费网站| 自拍欧美九色日韩亚洲蝌蚪91 | 久久精品国产亚洲网站| 亚洲电影在线观看av| 观看美女的网站| 国产色婷婷99| 欧美日韩精品成人综合77777| 久久久色成人| 一个人看的www免费观看视频| 国产亚洲午夜精品一区二区久久| 中文字幕免费在线视频6| 日韩欧美 国产精品| 大片免费播放器 马上看| 91精品一卡2卡3卡4卡| 精品午夜福利在线看| 久久精品国产a三级三级三级| 国产乱人视频| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| 观看av在线不卡| 久久久久国产精品人妻一区二区| 国产高清不卡午夜福利| 在线观看美女被高潮喷水网站| 在现免费观看毛片| 一级毛片久久久久久久久女| 亚洲人成网站高清观看| 少妇丰满av| 三级国产精品片| 国产熟女欧美一区二区| 综合色丁香网| 亚洲丝袜综合中文字幕| 日本wwww免费看| 超碰97精品在线观看| 黑丝袜美女国产一区| 丰满迷人的少妇在线观看| 男人添女人高潮全过程视频| 看非洲黑人一级黄片| 欧美xxxx性猛交bbbb| 99精国产麻豆久久婷婷| 国产综合精华液| 精品视频人人做人人爽| 日本av免费视频播放| 亚洲人与动物交配视频| 日本-黄色视频高清免费观看| 久久国产亚洲av麻豆专区| 久久亚洲国产成人精品v| 久久人人爽人人爽人人片va| 99久久人妻综合| 国产一区二区三区av在线| 下体分泌物呈黄色| 亚洲av二区三区四区| 嫩草影院新地址| 熟女av电影| 亚洲国产欧美在线一区| 狠狠精品人妻久久久久久综合| 亚洲精品久久久久久婷婷小说| 色婷婷久久久亚洲欧美| 精品国产露脸久久av麻豆| 啦啦啦中文免费视频观看日本| 午夜视频国产福利| 国产男女内射视频| 涩涩av久久男人的天堂| 欧美日韩精品成人综合77777| 久久久久国产精品人妻一区二区| 九九在线视频观看精品| 在线观看免费高清a一片| 国产精品不卡视频一区二区| 一本色道久久久久久精品综合| 免费看av在线观看网站| 久久久亚洲精品成人影院| 一本久久精品| 亚洲中文av在线| 美女xxoo啪啪120秒动态图| 精品一品国产午夜福利视频| av又黄又爽大尺度在线免费看| 久久精品国产a三级三级三级| 国产一区二区三区av在线| 久久人人爽人人片av| 黑丝袜美女国产一区| 日韩强制内射视频| a级一级毛片免费在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美 日韩 精品 国产| av播播在线观看一区| 欧美精品亚洲一区二区| 免费观看在线日韩| 一区二区三区乱码不卡18| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 欧美人与善性xxx| 舔av片在线| 精品一品国产午夜福利视频| 777米奇影视久久| 国产色爽女视频免费观看| 国产成人91sexporn| 丝瓜视频免费看黄片| 亚洲欧洲国产日韩| 亚洲精品国产av蜜桃| av国产免费在线观看| 国产v大片淫在线免费观看| 国产午夜精品久久久久久一区二区三区| 免费黄色在线免费观看| 亚洲精品自拍成人| 国产精品99久久99久久久不卡 | 久久国产乱子免费精品| 国产视频内射| 日日啪夜夜撸| 少妇人妻精品综合一区二区| 纵有疾风起免费观看全集完整版| 国产在视频线精品| 人体艺术视频欧美日本| 国模一区二区三区四区视频| 免费av不卡在线播放| 男男h啪啪无遮挡| 精品酒店卫生间| 久久6这里有精品| 亚洲欧美日韩东京热| 免费观看性生交大片5| 亚洲无线观看免费| 狂野欧美激情性xxxx在线观看| 80岁老熟妇乱子伦牲交| 国产精品三级大全| 欧美97在线视频| 午夜激情久久久久久久| 在线观看一区二区三区激情| 亚洲国产av新网站| 国产乱人偷精品视频| 久久久久久久久大av| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 在线观看一区二区三区激情| 欧美激情国产日韩精品一区| 一本色道久久久久久精品综合| 久久精品国产a三级三级三级| 亚洲伊人久久精品综合| 尤物成人国产欧美一区二区三区| 久久久久久伊人网av| 中文在线观看免费www的网站| 免费不卡的大黄色大毛片视频在线观看| 中文字幕制服av| 99久久精品热视频| videos熟女内射| 亚洲av男天堂| 国产 一区 欧美 日韩| 久久久久视频综合| 精品视频人人做人人爽| 国产精品欧美亚洲77777| 国产在线男女| 搡老乐熟女国产| 中国美白少妇内射xxxbb| 日韩成人av中文字幕在线观看| 亚洲最大成人中文| 91精品伊人久久大香线蕉| 亚洲成人av在线免费| 我要看黄色一级片免费的| 身体一侧抽搐| 午夜激情福利司机影院| av在线观看视频网站免费| 亚洲成人一二三区av| 日本午夜av视频| 久久久久久久国产电影| 精品一区二区免费观看| 精品亚洲成国产av| 又粗又硬又长又爽又黄的视频| 欧美变态另类bdsm刘玥| 丰满迷人的少妇在线观看| 在线精品无人区一区二区三 | 欧美高清成人免费视频www| 青春草亚洲视频在线观看| 欧美激情国产日韩精品一区| 女人十人毛片免费观看3o分钟| 国产一区二区三区综合在线观看 | 一区二区三区免费毛片| 你懂的网址亚洲精品在线观看| 在现免费观看毛片| 精品人妻偷拍中文字幕| 韩国高清视频一区二区三区| 亚洲精品,欧美精品| 晚上一个人看的免费电影| 草草在线视频免费看| 又黄又爽又刺激的免费视频.|