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

    無源時間反轉(zhuǎn)聚焦方法

    2016-12-07 09:35:53劉偉付永慶許達(dá)
    關(guān)鍵詞:信源觀測點方根

    劉偉,付永慶,許達(dá)

    (1. 哈爾濱工程大學(xué) 信息與通信工程學(xué)院,黑龍江 哈爾濱,150001;2. 哈爾濱工程大學(xué) 水聲工程學(xué)院,黑龍江 哈爾濱,150001)

    無源時間反轉(zhuǎn)聚焦方法

    劉偉1,付永慶1,許達(dá)2

    (1. 哈爾濱工程大學(xué) 信息與通信工程學(xué)院,黑龍江 哈爾濱,150001;2. 哈爾濱工程大學(xué) 水聲工程學(xué)院,黑龍江 哈爾濱,150001)

    為了提高時間反轉(zhuǎn)聚焦方法在電子戰(zhàn)中的適用性,提出一種基于電磁波波動方程的無源時間反轉(zhuǎn)聚焦方法。該方法以平面電磁波波動方程的時間反轉(zhuǎn)不變性為基礎(chǔ),利用時間反轉(zhuǎn)陣列接收并記錄一段時間的信源發(fā)射信號,然后對該信號進(jìn)行時間反轉(zhuǎn)處理,利用計算已知目標(biāo)域內(nèi)各觀測點處的虛擬接收和信號及其能量代替時間反轉(zhuǎn)后信號的物理傳輸過程,實現(xiàn)時間反轉(zhuǎn)后信號在信源處的重聚焦,其中能量最大點所對應(yīng)的觀測點即為信源點。在理論上驗證該方法的正確性,然后以均勻圓形時間反轉(zhuǎn)陣列為例進(jìn)行仿真實驗。研究結(jié)果表明:該無源時間反轉(zhuǎn)聚焦算法具有有效性和準(zhǔn)確性,在低信噪比下具有良好聚焦性能。

    無源時間反轉(zhuǎn);波動方程;聚焦;觀測點

    時間反轉(zhuǎn)由光學(xué)中的相位共軛引申而來,F(xiàn)INK等[1]提出并驗證了時間反轉(zhuǎn)技術(shù)在超聲實驗中的有效性,隨后SONG等[2]對該技術(shù)在海洋實驗中的應(yīng)用進(jìn)行了研究,直到2004年,時間反轉(zhuǎn)開始被引入到電磁領(lǐng)域[3]。到目前為止,時間反轉(zhuǎn)技術(shù)已經(jīng)在超聲碎石[4]、無損探傷[5]、目標(biāo)探測[6?8]、無線通信[9?10]、成像[11?12]和雷達(dá)[13?14]等領(lǐng)域得到廣泛研究。其基本原理為:收發(fā)合置的傳感器陣列接收輻射源發(fā)射的信號并對該信號進(jìn)行時間反轉(zhuǎn),然后將時間反轉(zhuǎn)后的信號物理發(fā)射到空間中,則經(jīng)過反向物理傳輸,時間反轉(zhuǎn)后信號會在信源處實現(xiàn)重聚焦[15]。從以上原理可以看出:經(jīng)典時間反轉(zhuǎn)理論需要將時間反轉(zhuǎn)后信號進(jìn)行物理傳輸,并在信源處設(shè)置觀察陣列來實現(xiàn)信號重聚焦,不能直接用于實際電子戰(zhàn)中。水聲領(lǐng)域提出了一種虛擬時間反轉(zhuǎn)鏡來實現(xiàn)對輻射源目標(biāo)的被動檢測[16],但該方法需要預(yù)先對信道的沖擊響應(yīng)函數(shù)進(jìn)行估計,不適合用于復(fù)雜多變的電子戰(zhàn)環(huán)境。為了解決以上問題,本文作者提出了一種無源時間反轉(zhuǎn)聚焦方法,該方法只利用接收陣列接收一段信源的發(fā)射信號,接著將該信號進(jìn)行時間反轉(zhuǎn)處理,然后利用時間反轉(zhuǎn)后信號計算已知目標(biāo)域中各觀測點的虛擬接收和信號及其能量,即可實現(xiàn)在信源處的信號重聚焦。以平面電磁波波動方程為基礎(chǔ),從理論上分析了無源時間反轉(zhuǎn)聚焦方法的正確性,并給出了算法步驟,最后以均勻圓形時間反轉(zhuǎn)陣列為例,驗證了該算法的有效性。

    1 波動方程的時間反轉(zhuǎn)不變性

    時間反轉(zhuǎn)在時域中是對信號的一種逆序操作,即將信號在時序上進(jìn)行反轉(zhuǎn),在頻域內(nèi)可等效為相位共軛。以波動方程的時間反轉(zhuǎn)不變性為基礎(chǔ),把接收天線陣列接收到的目標(biāo)信號進(jìn)行時間反轉(zhuǎn)處理并重新發(fā)射到空間中,則會實現(xiàn)在目標(biāo)處的信號重聚焦。

    假設(shè)電磁場中的標(biāo)量位函數(shù)為?(r,t),則在無源場區(qū)域滿足的標(biāo)量波動方程為[17]:

    式中:μ和ε分別為介質(zhì)的介電常數(shù)和磁導(dǎo)率;r為電磁波傳播的位移矢量。

    考慮遠(yuǎn)場平面電磁波,則式(1)所示標(biāo)量波動方程的通解為:

    式中:f1和f2為任意函數(shù);,為電磁波傳播速度。

    可以看出,式(3)和(4)分別表示傳播方向相反的兩類波形,并且它們都是波動方程的解。

    對式(3)和(4)進(jìn)行時間反轉(zhuǎn)處理,可得:

    可見,?1(r,?t)與?2(r,?t)也是波動方程的解,因此,波動方程具有時間反轉(zhuǎn)不變性。這一特性說明由目標(biāo)激勵的每一個電磁波,經(jīng)過介質(zhì)傳播之后的場,必定會存在一個時間上與之相反的場沿著相同的路徑返回到目標(biāo)的位置,實現(xiàn)電磁波信號在目標(biāo)處的聚焦。

    2 無源時間反轉(zhuǎn)聚焦方法

    2.1時間反轉(zhuǎn)聚焦原理

    時間反轉(zhuǎn)鏡(time reversal mirror, TRM)陣列將接收到的信號進(jìn)行時間反轉(zhuǎn),然后重新發(fā)射到空間中,反向傳輸?shù)臅r間反轉(zhuǎn)信號就會實現(xiàn)在信源處的重聚焦,其原理如圖1所示。

    圖1 時間反轉(zhuǎn)聚焦模型Fig. 1 Model of time reversal focusing

    令用()Sω為信號頻譜,()Hω為信道的傳輸函數(shù),則信號經(jīng)過信道傳輸,到達(dá)TRM接收陣列時的頻譜可以表示為

    時域內(nèi)的時間反轉(zhuǎn)在頻域內(nèi)表示為復(fù)共軛,令RTR(ω)是R(ω)的復(fù)共軛,即

    將時間反轉(zhuǎn)后信號反向發(fā)射,并按原信道進(jìn)行反向傳輸可得到:

    式中:H?(ω)H(ω)為實、偶、正函數(shù),它在時間零點的傅里葉變換同相相加,得到最大值,實現(xiàn)信號重聚焦[18]。

    2.2無源時間反轉(zhuǎn)聚焦算法

    假設(shè)圖1中TRM陣列由N個天線陣元組成,且信源的發(fā)射信號為s(t)。定義信源與TRM陣列中第n個陣元天線的沖激響應(yīng)函數(shù)為hn( rn, t)(其中:1≤n≤N,rn為信源與第n個陣元之間的傳輸距離)。因此,第n個陣元天線接收到的信號為[19]

    寫成頻域形式為

    對各陣元的接收信號在時域上進(jìn)行時間反轉(zhuǎn),相當(dāng)于對其頻域信號進(jìn)行相位共軛,于是第n個陣元天線的時間反轉(zhuǎn)后信號可以表示為

    根據(jù)圖1所示的時間反轉(zhuǎn)聚焦原理,若想實現(xiàn)信號在信源處的重聚焦,需要將間反轉(zhuǎn)后信號進(jìn)行反向物理傳輸至目標(biāo)域,并在目標(biāo)域設(shè)置接收陣列以觀察時間反轉(zhuǎn)聚焦現(xiàn)象。但在非合作信源情況下(如敵方雷達(dá)或通信電臺等),實現(xiàn)時間反轉(zhuǎn)后信號的反向物理傳輸是沒有意義的。因此,為了解決非合作信源情況下的時間反轉(zhuǎn)聚焦問題,可將已知目標(biāo)域平均劃分成M個觀測點,并利用觀測點與各天線陣元之間的相對位置來計算各觀測點處的虛擬接收和信號來實現(xiàn)時間反轉(zhuǎn)聚焦。

    令目標(biāo)域中第k (1≤k≤M)個觀測點與第n個陣元天線之間的傳輸函數(shù)為Hkn(rkn,ω),其中,rkn為第k個觀測點與第n個陣元之間的傳輸距離,則可通過計算得到目標(biāo)域中第k個觀測點處虛擬接收到的第n個陣元反向傳輸?shù)男盘枮?/p>

    由式(21)可知:第n個接收天線的時間反轉(zhuǎn)信號達(dá)到最大值的時間僅與信源的發(fā)射信號有關(guān),且各天線陣元的時間反轉(zhuǎn)信號達(dá)到最大值的時間相同,都是在觀測點與信源點重合時達(dá)到最大值,即在信源處獲得重聚焦。因此可通過計算目標(biāo)域中各觀測點處虛擬接收到的N個天線陣元時間反轉(zhuǎn)信號的和信號的能量值來實現(xiàn)在信源處的重聚焦,其中最大能量點所對應(yīng)的觀測點即為信源點。

    由電磁平面波的波動方程的解可將目標(biāo)與TRM陣列中第n個陣元天線的沖激響應(yīng)函數(shù)寫為以下形式:

    計算目標(biāo)域中所有觀測點的信號能量,其中能量最大值所對應(yīng)的觀測點即為目標(biāo)信源點。因此,通過計算目標(biāo)域中各觀測點的信號能量值,即可實現(xiàn)時間反轉(zhuǎn)信號在目標(biāo)處的重聚焦。

    通過以上分析可將無源時間反轉(zhuǎn)聚焦算法分為以下步驟:

    1) 利用TRM陣列接收并記錄一段時間為T的信源發(fā)射信號;

    2) 將接收信號進(jìn)行時間反轉(zhuǎn)處理;

    3) 確定已知目標(biāo)搜索域中各目標(biāo)觀測點位置;

    4) 通過計算得到各目標(biāo)觀測點的虛擬接收和信號;

    5) 計算各觀測點的信號能量,實現(xiàn)在信源處的信號重聚焦。

    3 仿真分析

    考慮圖2所示均勻圓形時間反轉(zhuǎn)陣列,忽略其相干性,通道不一致性等,其中N個天線陣元均勻分布在半徑為r的圓周上,陣列的幾何中心與坐標(biāo)原點重合。為實現(xiàn)無源時間反轉(zhuǎn)聚焦,選取圖3所示為目標(biāo)搜索域,其中M個觀測點均勻分布在半徑為r′的圓周上,且每個觀測點的坐標(biāo)為Pk( r′,θk)。假設(shè)目標(biāo)位于目標(biāo)搜索域內(nèi)且方位角為82°,目標(biāo)信源發(fā)射頻率為10 MHz的正弦電磁波。仿真中取N=10,r=20 m,圖4所示為波束指向90°時的時間反轉(zhuǎn)陣列方向圖。目標(biāo)搜索域取r′=20 km,觀測點數(shù)M=360,即搜索步進(jìn)為1°,時間反轉(zhuǎn)陣列的信號記錄時間為1 ms,則目標(biāo)域內(nèi)各觀測點的虛擬接收和信號歸一化能量值分布如圖5所示。從圖5可以看出:目標(biāo)域內(nèi)的歸一化能量值在方位角為82°的觀測點處取得最大值,此觀測點即為信源點,驗證了文中所提無源時間反轉(zhuǎn)聚焦算法的有效性。

    圖2 均勻圓形時間反轉(zhuǎn)陣列Fig. 2 Uniform circular time reversal array

    圖3 目標(biāo)搜索域Fig. 3 Targets searching area

    圖4 時間反轉(zhuǎn)陣列方向圖Fig. 4 Directivity diagram of time reversal array

    圖5 目標(biāo)域內(nèi)歸一化能量值分布Fig. 5 Distribution of normalized energy values in target area

    改變時間反轉(zhuǎn)陣列陣元數(shù)N,分別取8,10和15,其他參數(shù)保持不變。將時間反轉(zhuǎn)陣列接收到的信號進(jìn)行時間反轉(zhuǎn),然后計算目標(biāo)域內(nèi)各觀測點的虛擬接收和信號能量值,圖6所示為目標(biāo)域內(nèi)各觀測點的和信號能量分布情況。從圖6可以看出:不同陣元數(shù)情況下,無源時間反轉(zhuǎn)聚焦算法都可以有效實現(xiàn)在信源處的信號重聚焦,聚焦方位為82°,而且陣元數(shù)越大,信源處的聚焦峰值也越大,這說明陣元數(shù)增加可以增強(qiáng)無源時間反轉(zhuǎn)算法在目標(biāo)處的聚焦效果。

    為觀察噪聲對無源時間反轉(zhuǎn)聚焦算法的影響,在仿真中加入信噪比不同的高斯白噪聲,其中時間反轉(zhuǎn)陣列的陣元數(shù)N分別取8,10和15,其他參數(shù)保持不變,算法時間反轉(zhuǎn)聚焦的均方根誤差eRMSE計算公式為

    圖6 不同陣元數(shù)下的無源時間反轉(zhuǎn)聚焦Fig. 6 Passive time reversal focusing in different numbers of element

    式中:L為樣本數(shù);Xi為樣本點;1≤i≤L;為樣本平均值。

    圖7所示為執(zhí)行200次蒙特卡洛實驗時,不同信噪比情況下無源時間反轉(zhuǎn)聚焦算法的均方根誤差。從圖7可以看出:隨著信噪比和陣元數(shù)的增大,該算法聚焦的均方根誤差都有逐漸減小的趨勢,而且在信噪比為0 dB時,3種情況下的聚焦均方根誤差都在0.8°以下,顯示了該方法在較低信噪比下的良好聚焦性能。圖8所示為陣元數(shù)為8的情況下,本文方法與傳統(tǒng)MUSIC和RB-MUSIC[20]在不同信噪比下的均方根誤差對比曲線。從圖8可以看出:文中所給出的無源時間反轉(zhuǎn)算法具有更小的均方根誤差,可以獲得更好的目標(biāo)檢測性能。

    為研究不同搜索步進(jìn),即不同觀測點數(shù)對算法聚焦性能的影響,在信噪比為0 dB的情況下取搜索步進(jìn)分別為1°,0.5°和0.1°,然后通過計算各觀測點處的虛擬接收和信號能量值實現(xiàn)在信源處的信號重聚焦。執(zhí)行200次蒙特卡洛實驗可得不同搜索步進(jìn)情況下的聚焦均方根誤差如表1所示。從表1可以看出:隨著搜索步進(jìn)的減小,即觀測點數(shù)的增加,3種陣元數(shù)情況下的聚焦誤差也隨之減小,同時,表1也顯示出了該無源時間反轉(zhuǎn)聚焦算法在信噪比為0 dB情況下的有效準(zhǔn)確聚焦性能。

    圖7 不同信噪比下的聚焦均方根誤差Fig. 7 RMSE of focusing in different SNR

    圖8 3種方法在不同信噪比下的均方根誤差Fig. 8 RMSE of three methods in different SNR

    表1 不同搜索步進(jìn)下的聚焦均方根誤差Table 1 RMSE of focusing in different searching steppings

    4 結(jié)論

    1) 提出了一種無源時間反轉(zhuǎn)聚焦方法,該方法不需要進(jìn)行時間反轉(zhuǎn)后信號的物理傳輸,也不需要對信道的沖擊響應(yīng)函數(shù)進(jìn)行預(yù)先估計,而是通過計算目標(biāo)域中各觀測點處虛擬接收到的N個天線陣元時間反轉(zhuǎn)信號的和信號的能量來實現(xiàn)在信源處的重聚焦,其中最大能量所對應(yīng)的觀測點即為信源點。

    2) 本文方法具有有效性,在低信噪比情況下具有良好聚焦性能。

    3) 無源時間反轉(zhuǎn)聚焦方法只需要對接收信號進(jìn)行時間反轉(zhuǎn)處理和虛擬能量計算即可實現(xiàn)聚焦,因此更適合應(yīng)用于電子戰(zhàn)中對敵方輻射源的偵測,也可用于任何需要對目標(biāo)進(jìn)行被動檢測的應(yīng)用環(huán)境。

    [1] FINK M, PRADA C, WU F, et al. Self focusing in inhomogeneous media with time reversal acoustic mirrors[C]// IEEE Ultrasonics Symposium Proceedings. Montreal: IEEE Press, 1989: 681?686.

    [2] SONG H C, KUPERMAN W A, HOGKISS W S, et al. Iterative time reversal in the ocean[J]. Journal of the Acoustical Society of America, 1999, 105(6): 3176?3184.

    [3] LEROSEY G, de ROSNY J, TOURIN A, et al. Time reversal of electromagnetic waves[J]. Physical Review Letters, 2004, 92 (19): 193904-1?193904-3.

    [4] THOMAS J L, FINK M. Self focusing on extended objects with time reversal mirror, application to lithotripsy[C]// IEEE Ultrasonics Symposium Proceedings. Cannes, France: IEEE Press, 1994: 1809?1814.

    [5] 趙乃至, 閻石, 齊霽. 基于時間反轉(zhuǎn)法的管道周向裂縫損傷檢測[J]. 沈陽建筑大學(xué)學(xué)報, 2013, 29(1): 44?49. ZHAO Naizhi, YAN Shi, QI Ji. The pipeline circumferential cracks damage detection based on time reversal method[J]. Journal of Shenyang Jianzhu University, 2013, 29(1): 44?49.

    [6] 時潔, 楊德森, 時勝國. 基于多途模型匹配的淺海噪聲源聚焦定位方法[J]. 電子學(xué)報, 2013, 41(3): 575?581. SHI Jie, YANG Desen, SHI Shengguo. Research on noise sound source location method in shallow water based on the multi-path model match[J]. ACTA Electronica Sinica, 2013, 41(3):575?581.

    [7] FU Yongqing, JIANG Yulei, LIU Zhanya. Near-field source localization method and application using the time reversal mirror technique[J]. Journal of Electronics (China), 2011, 28(4): 531?538.

    [8] 李壯, 喬鋼, 王健培, 等. 基于虛擬時間反轉(zhuǎn)鏡的短基線定位研究[J]. 應(yīng)用聲學(xué), 2012, 31(4): 256?261. LI Zhuang, QIAO Gang, WANG Jianpei, et al. Short baseline positioning based on virtual time reversal mirror[J]. Applied Acoustics, 2012, 31(4): 256?261.

    [9] KHAN W N, ZUBAIR M, WYNE S, et al. Performance evaluation of time-reversal on measured 60 GHz wireless channels[J]. Wireless Personal Communications, 2013, 71(1): 707?701.

    [10] CHEN Yan, YANG Yuhan, HAN Feng, et al. Time-reversal wideband communications[J]. IEEE Signal Processing Letters, 2013, 20(12): 1219?1222.

    [11] CHOI H, OGAWA Y T, OHGANE T. Time-reversal MUSIC imaging with time-domain gating technique[J]. IEICE Transactions on Communications, 2012, E95-B(7): 2377?2385.

    [12] LIU Xiaofei, WANG Bingzhou. Transmitting-mode time reversal imaging using MUSIC algorithm for surveillance in wireless sensor network[J]. IEEE Transactions on Antennas and Propagation, 2012, 60(1): 220?230.

    [13] JIN Y W, JOSE M F M, NICHOLAS O D. Time reversal in multiple-input multiple-output radar[J]. IEEE Journal of Selected Topics in Signal Processing, 2010, 4(1): 210?225.

    [14] FOROOZAN F, ASIF A, JIN Y W, et al. Direction finding algorithms for time reversal MIMO radars[C]// IEEE Workshop on Statistical Signal Processing. Nice: IEEE Press, 2011: 433?436.

    [15] FINK M. Time-reversed acoustics[J]. Reports on Progress in Physics, 2000, 63(5): 91?97.

    [16] 殷敬偉, 惠俊英. 虛擬時間反轉(zhuǎn)鏡技術(shù) [C]// 全國水聲學(xué)學(xué)術(shù)會議論文集. 福建: 中國聲學(xué)學(xué)會, 2005: 98?100. YIN Jingwei, HUI Junying. The technique of virtual time reversal mirror[C]// Underwater Acoustics Conference Proceedings. Fujian: China Acoustics Institute, 2005: 98?100.

    [17] 譚思浩. 基于時間反轉(zhuǎn)電磁場逆的超寬帶雷達(dá)成像算法研究[D]. 南昌: 南昌大學(xué)物理系, 2013: 38?41. TAN Sihao. Research on ultrawideband radar imaging algorithm using time-reversal electromagnetic inversion technique[D]. Nanchang: Department of Physics, Nanchang University, 2013: 38?41.

    [18] 生雪莉, 惠俊英, 梁國龍. 時間反轉(zhuǎn)鏡用于被動檢測技術(shù)的研究[J]. 應(yīng)用聲學(xué), 2005, 24(6): 351?358. SHENG Xueli, HUI Junying, LIANG Guolong. Study on passive detection technology based on an acoustic time-reversal mirror[J]. Applied Acoustics, 2005, 24(6): 351?358.

    [19] 楊伏洲, 王海燕, 申曉紅, 等. 基于時間反轉(zhuǎn)的非均勻線陣列超指向性陣元分布模型[J]. 上海交通大學(xué)學(xué)報, 2013, 47(12): 1907?1917. YANG Fuzhou, WANG Haiyan, SHEN Xiaohong, et al. Super-direction element distribution model of NLA based on TR[J]. Journal of Shanghai Jiao Tong University, 2013, 47(12): 1907?1917.

    [20] MATHEWS C P, ZOLTOWSKI M D. Eigenstructure techniques for 2-D angle estimation with uniform circular arrays[J]. IEEE Transactions on Signal Processing. 1994, 42(9): 2395?2407.

    (編輯 趙俊)

    Passive time reversal focusing method

    LIU Wei1, FU Yongqing1, XU Da2
    (1. College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China; 2. College of Underwater Acoustic Engineering, Harbin Engineering University, Harbin 150001, China)

    In order to improve the applicability for time reversal focusing method in electronic warfare, a passive time reversal focusing method based on wave equation of electromagnetic wave was proposed. Based on the time reversal invariance of wave equation in electromagnetic plane wave, a part of signals emitted from source were

    and recorded by the time reversal array, and then the recorded signals were time reversed. After that, by calculating the virtually

    sum signals and energy values of every observed point in the given targets area instead of the physical transmission process of time-reversed signals, are refocused at source and the observed point whose energy value of virtual

    sum signals is maximal is the source point. The correctness of the proposed method was verified theoretically, and then, the simulations with the time reversal array which is a uniform circular array were performed. The results show that the passive time reversal focusing method has effectiveness and accuracy, and the proposed method has good focusing performance at low SNR.

    passive time reversal; wave equation; focusing; observed points

    TN911.7

    A

    1672?7207(2016)03?0800?07

    10.11817/j.issn.1672-7207.2016.03.013

    2015?03?18;

    2015?06?10

    國家自然科學(xué)基金資助項目(61172038) (Project(61172038) supported by the National Natural Science Foundation of China)

    付永慶,教授,博士生導(dǎo)師,從事陣列信號處理的研究;E-mail: fuyongqing@hrbeu.edu.cn

    猜你喜歡
    信源觀測點方根
    方根拓展探究
    基于極化碼的分布式多信源信道聯(lián)合編碼
    無線電工程(2022年4期)2022-04-21 07:19:44
    高速公路網(wǎng)連續(xù)式交通量調(diào)查觀測點布設(shè)方法研究
    智能城市(2021年3期)2021-04-12 04:40:50
    洛陽市老城區(qū)西大街空間形態(tài)與熱環(huán)境耦合關(guān)系實測研究
    綠色科技(2019年12期)2019-07-15 11:13:02
    信源控制電路在功率容量測試系統(tǒng)中的應(yīng)用
    電子世界(2017年16期)2017-09-03 10:57:36
    均方根嵌入式容積粒子PHD 多目標(biāo)跟蹤方法
    揭開心算方根之謎
    張掖市甘州區(qū)代表性觀測點地下水位變化特征分析
    信源自動切換裝置的設(shè)計及控制原理
    基于升降溫全曲線的鋼筋混凝土梁溫度場分析
    男人狂女人下面高潮的视频| 99久久成人亚洲精品观看| 麻豆精品久久久久久蜜桃| 热99re8久久精品国产| 在线免费十八禁| 嫩草影院入口| 一区二区三区免费毛片| 国产黄片视频在线免费观看| 日日干狠狠操夜夜爽| 嫩草影院入口| 毛片女人毛片| 亚洲最大成人手机在线| 成人性生交大片免费视频hd| 在线免费观看不下载黄p国产| 少妇熟女欧美另类| 一本一本综合久久| 国产人妻一区二区三区在| 三级毛片av免费| 国产视频内射| 内地一区二区视频在线| 中文在线观看免费www的网站| 精品久久久久久久久久久久久| 高清在线视频一区二区三区 | 免费看美女性在线毛片视频| 亚洲精品国产成人久久av| 免费观看人在逋| 男女视频在线观看网站免费| 日韩高清综合在线| 久久久色成人| 国产真实乱freesex| 久久久国产成人免费| 久久精品国产亚洲网站| 一级av片app| 国产精品永久免费网站| 久久久久久久亚洲中文字幕| 国产激情偷乱视频一区二区| 97超碰精品成人国产| 精品少妇黑人巨大在线播放 | 午夜福利在线观看免费完整高清在 | 日韩,欧美,国产一区二区三区 | 永久网站在线| 自拍偷自拍亚洲精品老妇| 欧美日韩精品成人综合77777| 男女做爰动态图高潮gif福利片| av免费在线看不卡| 国产精品一区二区三区四区久久| 国产一区二区亚洲精品在线观看| 国产精品精品国产色婷婷| 亚洲欧美日韩卡通动漫| 嫩草影院入口| 国产片特级美女逼逼视频| 国产黄色视频一区二区在线观看 | 亚洲一区高清亚洲精品| 精品日产1卡2卡| 老师上课跳d突然被开到最大视频| 一级毛片aaaaaa免费看小| 亚洲av成人精品一区久久| 亚洲av中文字字幕乱码综合| 国产av麻豆久久久久久久| 日本三级黄在线观看| 国产久久久一区二区三区| 丝袜喷水一区| 不卡视频在线观看欧美| 在线观看66精品国产| 国产精品久久久久久久电影| 天堂中文最新版在线下载 | av天堂在线播放| 在线a可以看的网站| 成人毛片a级毛片在线播放| 两个人的视频大全免费| 又粗又爽又猛毛片免费看| 久久国产乱子免费精品| 97人妻精品一区二区三区麻豆| 欧美成人免费av一区二区三区| 校园人妻丝袜中文字幕| 黄色视频,在线免费观看| 变态另类丝袜制服| 91久久精品国产一区二区三区| 国产av不卡久久| 在线观看美女被高潮喷水网站| 99九九线精品视频在线观看视频| 听说在线观看完整版免费高清| 久久久精品欧美日韩精品| 男女视频在线观看网站免费| 精品欧美国产一区二区三| 午夜免费激情av| 最近2019中文字幕mv第一页| 日韩 亚洲 欧美在线| 成人鲁丝片一二三区免费| 精品一区二区三区视频在线| 桃色一区二区三区在线观看| 别揉我奶头 嗯啊视频| 久久精品国产99精品国产亚洲性色| 亚洲精品456在线播放app| 午夜精品国产一区二区电影 | 欧美xxxx黑人xx丫x性爽| 久久久久九九精品影院| 能在线免费观看的黄片| 欧美另类亚洲清纯唯美| 啦啦啦韩国在线观看视频| 国产精品乱码一区二三区的特点| 成人毛片60女人毛片免费| 美女大奶头视频| 深夜精品福利| 成人午夜高清在线视频| 国产精品久久久久久精品电影小说 | 男女啪啪激烈高潮av片| 97人妻精品一区二区三区麻豆| 国产男人的电影天堂91| 深爱激情五月婷婷| 午夜免费激情av| 亚洲一区高清亚洲精品| 天堂影院成人在线观看| 亚洲欧美日韩高清在线视频| 97超碰精品成人国产| 白带黄色成豆腐渣| 欧美激情国产日韩精品一区| h日本视频在线播放| 大型黄色视频在线免费观看| 久久久久久伊人网av| 噜噜噜噜噜久久久久久91| 性插视频无遮挡在线免费观看| 成人性生交大片免费视频hd| 精品99又大又爽又粗少妇毛片| 午夜视频国产福利| 久久鲁丝午夜福利片| 在线观看av片永久免费下载| 日韩欧美一区二区三区在线观看| 久久久久久国产a免费观看| 亚洲精品久久国产高清桃花| 麻豆成人午夜福利视频| 精品人妻熟女av久视频| 国产在视频线在精品| 国产午夜精品论理片| 国内少妇人妻偷人精品xxx网站| 亚洲av男天堂| 偷拍熟女少妇极品色| 在线免费十八禁| 中文字幕人妻熟人妻熟丝袜美| 国产成人freesex在线| 国产爱豆传媒在线观看| 久久婷婷人人爽人人干人人爱| 日本av手机在线免费观看| 可以在线观看的亚洲视频| 国产成人a区在线观看| 久久鲁丝午夜福利片| 蜜臀久久99精品久久宅男| 国产av不卡久久| 精品久久久久久久末码| 男人狂女人下面高潮的视频| 美女被艹到高潮喷水动态| 小说图片视频综合网站| 国产一区二区在线观看日韩| 最近视频中文字幕2019在线8| 久久久久网色| 国产在视频线在精品| 老熟妇乱子伦视频在线观看| 国产真实伦视频高清在线观看| 国产精品三级大全| 日本一二三区视频观看| 人妻久久中文字幕网| 亚洲成av人片在线播放无| 国产精品一区二区性色av| 高清日韩中文字幕在线| 91久久精品国产一区二区成人| a级毛色黄片| 色尼玛亚洲综合影院| 99热精品在线国产| 人人妻人人澡人人爽人人夜夜 | 国产成人aa在线观看| 国产精品麻豆人妻色哟哟久久 | 国产日本99.免费观看| 成人综合一区亚洲| 国产亚洲精品久久久久久毛片| 精华霜和精华液先用哪个| 中国国产av一级| 国产麻豆成人av免费视频| 精品久久久噜噜| 一级黄色大片毛片| 男女做爰动态图高潮gif福利片| 69人妻影院| 国产视频内射| 久久久欧美国产精品| 色哟哟哟哟哟哟| 久久久久久伊人网av| 嫩草影院入口| 免费黄网站久久成人精品| 成人永久免费在线观看视频| 久久热精品热| 有码 亚洲区| 亚洲国产色片| 国产色爽女视频免费观看| 免费av不卡在线播放| avwww免费| 免费人成视频x8x8入口观看| 日韩一区二区三区影片| 可以在线观看毛片的网站| 永久网站在线| 爱豆传媒免费全集在线观看| 美女cb高潮喷水在线观看| 亚洲第一电影网av| 亚洲国产精品合色在线| 午夜福利在线在线| 青春草亚洲视频在线观看| 国产片特级美女逼逼视频| 51国产日韩欧美| 国产精品久久电影中文字幕| 国产伦在线观看视频一区| 黄色一级大片看看| 在线观看av片永久免费下载| 久久99精品国语久久久| av国产免费在线观看| 午夜亚洲福利在线播放| 久久人人爽人人片av| 国内精品久久久久精免费| 久久久久国产网址| 国产三级中文精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩强制内射视频| 好男人在线观看高清免费视频| 午夜久久久久精精品| 日本熟妇午夜| 大型黄色视频在线免费观看| 亚洲国产欧美人成| 日本一二三区视频观看| 国产又黄又爽又无遮挡在线| 亚洲国产精品合色在线| www日本黄色视频网| 丰满乱子伦码专区| 麻豆国产av国片精品| 高清在线视频一区二区三区 | 青春草视频在线免费观看| 夫妻性生交免费视频一级片| 免费观看a级毛片全部| 一级毛片久久久久久久久女| 精品无人区乱码1区二区| 中国美白少妇内射xxxbb| 国产午夜福利久久久久久| 99热只有精品国产| 男插女下体视频免费在线播放| 久久久久久久久久久丰满| 国产老妇女一区| 久久久久久久久大av| 男人的好看免费观看在线视频| 亚洲图色成人| 一区二区三区高清视频在线| 亚洲自偷自拍三级| 亚洲av免费高清在线观看| 国产成人午夜福利电影在线观看| 91狼人影院| 一个人看的www免费观看视频| 午夜a级毛片| av视频在线观看入口| 欧美丝袜亚洲另类| 久久精品国产99精品国产亚洲性色| 少妇熟女欧美另类| 美女被艹到高潮喷水动态| 精品一区二区三区人妻视频| 亚洲精品粉嫩美女一区| av女优亚洲男人天堂| 精品久久久久久久人妻蜜臀av| 日韩人妻高清精品专区| 亚洲无线在线观看| 国产视频内射| 黄色配什么色好看| 好男人在线观看高清免费视频| 欧美xxxx黑人xx丫x性爽| 中文字幕人妻熟人妻熟丝袜美| 搞女人的毛片| 2022亚洲国产成人精品| 男女下面进入的视频免费午夜| 一个人看视频在线观看www免费| 日韩av不卡免费在线播放| 综合色av麻豆| 在现免费观看毛片| 内地一区二区视频在线| 99热全是精品| 国产视频内射| 黄色配什么色好看| 99热网站在线观看| 赤兔流量卡办理| 亚洲精品乱码久久久久久按摩| 成人特级av手机在线观看| 国产成人aa在线观看| 日韩成人伦理影院| 日韩中字成人| 99riav亚洲国产免费| 亚洲人成网站高清观看| 久久精品人妻少妇| 欧美性感艳星| 欧美高清性xxxxhd video| 欧美日韩一区二区视频在线观看视频在线 | 午夜激情欧美在线| av在线亚洲专区| 少妇被粗大猛烈的视频| 丝袜美腿在线中文| 国内久久婷婷六月综合欲色啪| 高清午夜精品一区二区三区 | 久久这里有精品视频免费| 此物有八面人人有两片| 久久综合国产亚洲精品| 亚洲av不卡在线观看| 全区人妻精品视频| 噜噜噜噜噜久久久久久91| 国产精品蜜桃在线观看 | 亚洲自偷自拍三级| 一区福利在线观看| 久久热精品热| 亚洲国产日韩欧美精品在线观看| 天堂网av新在线| 尾随美女入室| 国产三级在线视频| 日日干狠狠操夜夜爽| 国产又黄又爽又无遮挡在线| 欧美一区二区精品小视频在线| 日韩强制内射视频| 国产精品伦人一区二区| 老熟妇乱子伦视频在线观看| 日韩欧美精品v在线| 久久久久久国产a免费观看| 韩国av在线不卡| 亚洲在线观看片| 91久久精品电影网| 少妇裸体淫交视频免费看高清| 亚洲av不卡在线观看| 在线观看66精品国产| 久久精品国产清高在天天线| 亚洲国产精品成人久久小说 | 久久99蜜桃精品久久| 久久久久久九九精品二区国产| 偷拍熟女少妇极品色| 亚洲五月天丁香| 日本在线视频免费播放| 国产午夜精品论理片| 亚洲精品成人久久久久久| videossex国产| 亚洲国产精品国产精品| 少妇熟女欧美另类| 观看美女的网站| 高清毛片免费看| 亚洲欧美精品自产自拍| 亚洲国产欧美人成| 最近视频中文字幕2019在线8| 精品久久久久久成人av| 日韩大尺度精品在线看网址| 成人高潮视频无遮挡免费网站| 免费不卡的大黄色大毛片视频在线观看 | 国产精品蜜桃在线观看 | 国内精品美女久久久久久| 久久精品影院6| 女人被狂操c到高潮| 在线观看免费视频日本深夜| eeuss影院久久| 18禁在线无遮挡免费观看视频| av专区在线播放| 1000部很黄的大片| 少妇的逼好多水| 午夜激情欧美在线| 国产精品.久久久| 亚洲美女视频黄频| 婷婷亚洲欧美| 日韩成人伦理影院| 自拍偷自拍亚洲精品老妇| 一区二区三区免费毛片| 欧美潮喷喷水| 黄片wwwwww| av女优亚洲男人天堂| 国产精品国产高清国产av| av免费在线看不卡| 嫩草影院入口| 国产探花在线观看一区二区| 在线观看av片永久免费下载| 成人漫画全彩无遮挡| 欧美xxxx性猛交bbbb| 国产伦在线观看视频一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人免费观看视频高清| 国产国拍精品亚洲av在线观看| 亚洲国产最新在线播放| 久久亚洲国产成人精品v| 国产精品偷伦视频观看了| 黑丝袜美女国产一区| 青春草亚洲视频在线观看| videos熟女内射| 大香蕉久久成人网| 亚洲怡红院男人天堂| 欧美xxxx性猛交bbbb| 夫妻性生交免费视频一级片| 亚洲激情五月婷婷啪啪| 卡戴珊不雅视频在线播放| 日韩电影二区| 日韩人妻高清精品专区| 99九九在线精品视频| 国产亚洲精品第一综合不卡 | 亚洲av福利一区| 久久精品夜色国产| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线 | 天天影视国产精品| 日韩电影二区| 菩萨蛮人人尽说江南好唐韦庄| 国模一区二区三区四区视频| 欧美精品亚洲一区二区| 欧美激情极品国产一区二区三区 | 精品久久久久久久久av| av在线播放精品| 两个人免费观看高清视频| 高清午夜精品一区二区三区| 免费播放大片免费观看视频在线观看| 最近2019中文字幕mv第一页| av有码第一页| 蜜臀久久99精品久久宅男| 一边亲一边摸免费视频| 亚洲av不卡在线观看| 男女无遮挡免费网站观看| av免费观看日本| 国产精品免费大片| 久久精品国产亚洲av天美| 天天操日日干夜夜撸| av国产精品久久久久影院| 国产不卡av网站在线观看| 国产精品偷伦视频观看了| 久久97久久精品| av视频免费观看在线观看| 大香蕉久久成人网| 国产片特级美女逼逼视频| 高清在线视频一区二区三区| 街头女战士在线观看网站| 国产在视频线精品| 欧美激情国产日韩精品一区| 亚洲国产欧美在线一区| 国产视频内射| 亚洲精品乱码久久久久久按摩| 亚洲国产欧美在线一区| 九九爱精品视频在线观看| 国产一区有黄有色的免费视频| 午夜福利在线观看免费完整高清在| 另类亚洲欧美激情| 精品国产一区二区久久| 国产精品熟女久久久久浪| 久久久久久久久久成人| 午夜久久久在线观看| 亚洲精品乱码久久久v下载方式| av视频免费观看在线观看| 爱豆传媒免费全集在线观看| 夜夜骑夜夜射夜夜干| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| av一本久久久久| 一级毛片aaaaaa免费看小| 2022亚洲国产成人精品| a 毛片基地| 国产精品久久久久久精品古装| 我的老师免费观看完整版| 午夜久久久在线观看| 99久久人妻综合| 国产精品国产三级专区第一集| 我要看黄色一级片免费的| 亚洲国产精品999| 国产视频首页在线观看| 午夜福利视频在线观看免费| 搡女人真爽免费视频火全软件| 中文字幕精品免费在线观看视频 | 国产有黄有色有爽视频| 久久99一区二区三区| 亚洲欧美清纯卡通| 欧美成人午夜免费资源| 少妇被粗大猛烈的视频| 日韩 亚洲 欧美在线| 边亲边吃奶的免费视频| 久久久久久久精品精品| 亚洲精品一二三| 97在线视频观看| 男女无遮挡免费网站观看| 欧美精品一区二区大全| 在现免费观看毛片| 国产成人精品福利久久| 久久99蜜桃精品久久| 亚洲精品成人av观看孕妇| 亚洲欧美精品自产自拍| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站在线播| 欧美日韩成人在线一区二区| 99久久精品国产国产毛片| 在线精品无人区一区二区三| 国产一区二区在线观看av| 亚洲人成网站在线播| 妹子高潮喷水视频| 91国产中文字幕| 精品人妻熟女av久视频| av.在线天堂| 性色avwww在线观看| 黑丝袜美女国产一区| 视频中文字幕在线观看| 最近的中文字幕免费完整| 亚洲精品乱码久久久久久按摩| 免费大片18禁| 亚洲欧美成人精品一区二区| 亚洲av.av天堂| 十八禁高潮呻吟视频| 九九在线视频观看精品| 国产免费一区二区三区四区乱码| 在线播放无遮挡| 黑丝袜美女国产一区| 欧美日韩国产mv在线观看视频| 大陆偷拍与自拍| 午夜免费鲁丝| 欧美一级a爱片免费观看看| 日韩伦理黄色片| 如何舔出高潮| 久久久国产欧美日韩av| 成人国语在线视频| 免费黄色在线免费观看| 午夜福利在线观看免费完整高清在| 亚洲情色 制服丝袜| kizo精华| 成人漫画全彩无遮挡| 亚洲av综合色区一区| 成人手机av| 婷婷成人精品国产| 3wmmmm亚洲av在线观看| 欧美日韩国产mv在线观看视频| 国产av精品麻豆| 久久97久久精品| 成人亚洲精品一区在线观看| 亚洲精品aⅴ在线观看| 国产免费一区二区三区四区乱码| 老熟女久久久| 桃花免费在线播放| 亚洲五月色婷婷综合| 2022亚洲国产成人精品| 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 国产高清国产精品国产三级| 日韩成人伦理影院| 成人影院久久| 国产不卡av网站在线观看| 黄片播放在线免费| 国产爽快片一区二区三区| 最新中文字幕久久久久| 91精品国产九色| 亚洲欧美成人综合另类久久久| 少妇人妻精品综合一区二区| 2022亚洲国产成人精品| 水蜜桃什么品种好| 国产精品偷伦视频观看了| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 免费观看无遮挡的男女| 尾随美女入室| xxxhd国产人妻xxx| 日韩成人伦理影院| 国产综合精华液| 亚洲精品乱码久久久v下载方式| 美女国产视频在线观看| 午夜福利视频在线观看免费| 久久久久精品久久久久真实原创| 18禁动态无遮挡网站| 国产一区二区在线观看av| 亚洲精品乱码久久久v下载方式| 桃花免费在线播放| 狂野欧美激情性bbbbbb| 国产精品成人在线| 久久久久久久国产电影| av免费观看日本| 高清在线视频一区二区三区| 欧美xxⅹ黑人| 中文乱码字字幕精品一区二区三区| 久久免费观看电影| 久久国内精品自在自线图片| 欧美激情极品国产一区二区三区 | 一本色道久久久久久精品综合| 内地一区二区视频在线| 一级毛片 在线播放| 国产精品麻豆人妻色哟哟久久| 欧美国产精品一级二级三级| 97精品久久久久久久久久精品| 国产 一区精品| 免费黄网站久久成人精品| 人人妻人人澡人人看| 久久精品国产a三级三级三级| 亚洲av.av天堂| 精品视频人人做人人爽| 免费观看在线日韩| 曰老女人黄片| 黄片无遮挡物在线观看| 特大巨黑吊av在线直播| 成人毛片a级毛片在线播放| 日韩亚洲欧美综合| 国产精品一区二区在线观看99| 国产精品秋霞免费鲁丝片| 国产日韩欧美亚洲二区| 视频在线观看一区二区三区| 特大巨黑吊av在线直播| 欧美日韩综合久久久久久| 又粗又硬又长又爽又黄的视频| 两个人免费观看高清视频| 日韩制服骚丝袜av| 自线自在国产av| 搡女人真爽免费视频火全软件| 999精品在线视频| 七月丁香在线播放| 亚洲三级黄色毛片| 人人澡人人妻人| 午夜福利,免费看| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| 亚洲四区av| 国产成人精品久久久久久| 丝袜脚勾引网站| 国产 一区精品| 91成人精品电影| av免费在线看不卡| 久久久欧美国产精品| 色视频在线一区二区三区|