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

    通道幅相不一致時的全向測角誤差分析*1

    2015-03-10 01:18:40田超,文樹梁
    現(xiàn)代防御技術(shù) 2015年6期

    ?

    通道幅相不一致時的全向測角誤差分析*1

    田超,文樹梁

    (北京無線電測量研究所,北京100854)

    摘要:在理想情況下,基于均勻圓陣的米波全向雷達(dá)可聯(lián)合-1階、0階和和1階相位模式實現(xiàn)全方位的無模糊測角。然而在實際工程中,各接收通道幅相特性往往不一致,這將導(dǎo)致激勵出的相位模式中包含誤差項從而引起明顯的測角誤差,而且該誤差無法補(bǔ)償,只能通過校正通道間的幅相誤差或選擇合適陣列參數(shù)的方式來減小。為此,通過理論推導(dǎo)得到測角誤差與各接收支路幅相誤差之間的解析關(guān)系式,明確了幅相誤差對測角誤差影響的同時也為合理選陣列參數(shù)以減小幅相誤差引起的測角誤差提供了理論依據(jù),仿真分析驗證了理論誤差分析的正確性并給出了最優(yōu)的陣列直徑取值。

    關(guān)鍵詞:米波雷達(dá);均勻圓陣;全向測角;相位模式;幅相不一致;測角誤差分析

    0引言

    米波全向雷達(dá)因具有反隱身、抗反輻射導(dǎo)彈的能力以及其尺寸小、機(jī)動性強(qiáng)和研制成本低的特點而逐漸受到重視[1-3]。米波全向雷達(dá)利用均勻圓陣激勵出的相位模式測角,但不同于一般基于均勻圓陣的超分辨測角算法利用了各階相位模式[4-7],其通過聯(lián)合利用-1階、0階和1階相位模式實現(xiàn)互耦條件下全方位無模糊測角。然而,實際工程中,各個接收通道的幅相特性不同,這使得激勵出的相位模式存在誤差項,從而導(dǎo)致測角誤差的產(chǎn)生?,F(xiàn)有關(guān)于均勻圓陣幅相誤差對測角誤差的分析一般只基于數(shù)值仿真且由于測角體制的不同,一些結(jié)論并不適用于米波全向雷達(dá),文獻(xiàn)[8]分析了幅相誤差對均勻圓陣系統(tǒng)性能的影響,但其針對的是有指向性合成方向圖,而米波全向雷達(dá)的天線波束則是方位全向的,文獻(xiàn)[9]將幅相誤差與互耦系數(shù)誤差進(jìn)行了等效,但對于米波全向雷達(dá),幅相誤差矩陣為對角矩陣,互耦矩陣為對稱Toeplitz矩陣,二者無法等效。而文獻(xiàn)[10-13]提出的各種幅相誤差自校正算法雖可有效緩解幅相不一致的影響,但其校正一般并不在相位模式空間中實現(xiàn),而米波全向雷達(dá)則是利用-1階,0階和1階相位模式實現(xiàn)目標(biāo)角度測量,而且幅相自校正的計算復(fù)雜度高,給信號處理能力提出了較高的要求。

    因此,本文針對幅相不一致對測角的影響,先通過理論推導(dǎo)得到幅相不一致與測角誤差間的解析關(guān)系,從而明確了幅相不一致引起的測角誤差不僅與幅相不一致的程度有關(guān),還與陣列的參數(shù)有關(guān),進(jìn)而考慮通過合理選擇陣列參數(shù)的方式,在雷達(dá)設(shè)計階段將幅相不一致的影響在統(tǒng)計意義上減至最小。

    1無模糊方位全向測角原理

    若發(fā)射信號為窄帶信號,則陣元k接收到的信號可表示為

    (1)

    式中:k=1~N;λ為波長;s0(t)為相位中心接收的信號。

    圖1 N陣元均勻圓陣示意圖Fig.1 Schematic map of uniform circular array   with N elements

    根據(jù)均勻圓陣相位模式的相關(guān)理論[14-15],各階相位模式可通過陣元接收信號的FFT來實現(xiàn),即有

    (2)

    若滿足πd/λ

    (3)

    式中:l=-N+1~N-1;β=πdcosφ/λ;Jl(x)為第一類l階貝塞爾函數(shù)。

    考慮到s0(t)的相位未知,僅利用F1無法得到準(zhǔn)確的目標(biāo)方位角,因此為消除s0(t)相位的影響,可利用F1和F0。

    (4)

    然而,米波全向雷達(dá)出于機(jī)動性的考慮,天線陣列的口徑往往較小,陣元間的互耦比較嚴(yán)重,在互耦系數(shù)未知的情況下,利用式(4)會產(chǎn)生較大的測角誤差。

    若陣列的互耦矩陣為

    (5)

    通過推導(dǎo)可以得到,考慮陣元間的互耦效應(yīng)后:

    (6)

    由式(6)知,陣元間存在互耦時,κ的相位即為利用F1和F0測角產(chǎn)生的誤差。

    文獻(xiàn)[1]證明,利用F1和F-1進(jìn)行測角可以消除互耦的影響,即在陣元間存在嚴(yán)重互耦時亦有

    (7)

    但由式(7)知,利用F1和F-1進(jìn)行測角雖消除了互耦的影響,但亦引入了π的測角模糊,因此需要聯(lián)合式(6)的測角結(jié)果來解模糊。

    具體的解模糊過程如下:記由式(7)求得的方位角在[0,π]范圍內(nèi)的值記為為θ1,由式(6)求得的方位角即為θ2,θ2∈[0,2π]。若|θ1-θ2|∈[0.5π,1.5π],則判斷方位角為θ1+π,否則,判斷方位角為θ1。

    綜上所述,要獲得精度高且無模糊的目標(biāo)方位角,一方面要求θ1的精度高,另一方面要求θ2的誤差不超過π/2。

    為敘述方便,稱利用F1和F0的測角方法為非對稱式測角方法,稱利用F1和F-1的測角方法為對稱式測角方法。

    2幅相不一致時的測角誤差分析

    實際工程中,各接收通道幅相特性并不完全一致,這將導(dǎo)致測角存在誤差。出于減小幅相不一致引起的測角誤差的需要,明確通道間幅相不一致與測角誤差的關(guān)系十分必要。

    若記第k個接收支路相對于參考通道的幅度為ak,相位為vk。假設(shè)ak均為[1-δ,1+δ]間的均勻分布,vk為[-σ,σ]間的均勻分布,且各ak,vk相互獨立。若令ak=1+?k,則?k為[-δ,δ]間的均勻分布。

    當(dāng)各通道間的幅相不一致較小,即δ和σ較小時,有

    akexp(jvk)≈1+?k+jvk.

    (8)

    考慮各接收支路的幅相不一致后,結(jié)合式(2)和式(8)可得新的對稱式測角模型為

    (9)

    式中:φk=2πk/N。

    根據(jù)貝塞爾函數(shù)的性質(zhì)可知

    (10)

    將式(10)代入式(9)并通過推導(dǎo)可得

    (11)

    式中:

    若記θΔ為測角誤差,則可由式(11)得到

    tan(2θΔ)=-2(A+B)/[NJ1(β)].

    (12)

    各通道間幅相的不一致較小時,|θΔ|亦較小,此時由式(12)近似可得

    θΔ≈-(A+B)/[NJ1(β)].

    (13)

    由于?k,vk的均值為0且相互獨立,因此容易得到A和B的均值亦均為0且相互獨立,于是有

    E(θΔ)=[E(A)+E(B)]/[NJ1(β)]=0,

    (14)

    Var(θΔ)=[Var(A)+Var(B)]/[NJ1(β)]2.

    (15)

    根據(jù)?k和vk分別為服從[-δ,δ]和[-σ,σ]間的均勻分布并利用各?k和vk間相互獨立的特性,通過推導(dǎo)可得

    (16)

    由式(16)不難看出,幅相不一致對對稱式測角誤差的影響可分為兩部分之和,一部分僅由幅度不一致決定,另一部分僅由相位不一致決定。

    對于米波全向雷達(dá),一般限定d/λ<1,此時,Me=1,Mo=2,式(16)可簡化為

    (17)

    由式(17)知,在δ和σ固定的前提下,幅相不一致引起的測角誤差可以通過N和β的合理取值盡量減小。

    雖然N越大,幅相不一致的影響越小,但N的選取受到陣列直徑、互耦的影響和測角方法的固有誤差等因素的限制,因此,在探討如何基于式(17),通過陣列參數(shù)的選擇來減小幅相不一致對測角精度的影響時,只重點考慮β。

    β與陣列直徑和波長的比值以及目標(biāo)的仰角有關(guān),然而,目標(biāo)仰角為不可控制的因素且目標(biāo)仰角的變化可等效為陣列直徑與波長比值(簡稱直徑波長比)的變化,因此,就對稱式測角方法而言,陣列直徑波長比的最優(yōu)選擇應(yīng)使得式(17)最小。

    盡管角度模糊可以解決的前提條件對非對稱式測角方法的測角誤差要求非常寬松,即不超過π/2,但考慮到陣元間的互耦嚴(yán)重,為了給天線設(shè)計提供盡量大的余量,往往希望包括幅相不一致在內(nèi)的因素對測角精度的影響盡量地小。因此,幅相不一致對非對稱測角精度的影響同樣需要分析。

    考慮各接收支路的幅相不一致后,結(jié)合式(2)和式(8)可得新的非對稱式測角模型為

    (18)

    由于κ為常數(shù),在考慮由幅相不一致引起的測角誤差時,可不失一般性地令κ=1。此時,可以得到

    (19)

    若記θΔ為幅相不一致引起的測角誤差,則可由式(19)得到

    tanθΔ=(C+D)/N.

    (20)

    當(dāng)各通道間幅相的不一致較小時,|θΔ|亦較小,此時由式(20)近似可得

    θΔ≈(C+D)/N.

    (21)

    由于?k,vk的均值為0且相互獨立,因此容易得到C和D的均值亦均為0且相互獨立,于是有

    E(θΔ)=[E(C)+E(D)]/N=0,

    (22)

    Var(θΔ)=[Var(C)+Var(D)]/N2.

    (23)

    根據(jù)?k和vk分別為服從[-δ,δ]和[-σ,σ]間的均勻分布并利用各?k和vk間相互獨立的特性,通過推導(dǎo)可得

    (24)

    由式(24)不難看出,幅相不一致對非對稱式測角誤差的影響同樣可分為2部分之和,一部分僅由幅度不一致決定,另一部分僅由相位不一致決定。若限定d/λ<1,則式(24)可簡化為

    (25)

    為了解決對稱式測角方法的測角模糊問題且給天線設(shè)計留有余量,同樣需要將幅相不一致引起的非對稱測角誤差盡量控制在最小。因此,有必要通過陣列參數(shù)的選擇,主要是陣列直徑波長比,使得式(25)最小。

    綜合式(17)和式(25)可以看到:①采用對稱式測角算法時,為減小幅度不一致的影響,陣列直徑波長比的選擇宜使U1(β)盡量大,為減小相位不一致的影響,陣列直徑波長比的選擇宜使U2(β)盡量大;②采用非對稱式測角算法時,為減小幅度不一致的影響,陣列直徑波長比的選擇宜使U3(β)盡量大,為減小相位不一致的影響,陣列直徑波長比的選擇宜使U4(β)盡量大。

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

    本節(jié)主要通過數(shù)值仿真為合理選擇陣列直徑波長比提供依據(jù)并驗證所得到的幅相不一致引起的測角誤差公式的正確性。

    仰角為0°時,U1(β)和U2(β)隨陣列直徑波長比的變化分別如圖2和圖3所示。

    圖2 U1(β)隨d/λ變化的曲線Fig.2 Curves of U1(β)changed with d/λ

    圖3 U2(β)隨d/λ變化的曲線Fig.3 Curves of U2(β)changed with d/λ

    對于對稱式測角算法,由圖2可知,d/λ越大,幅度不一致的影響越大,d/λ=0時幅度不一致影響最?。挥蓤D3可知,d/λ為0.7時相位不一致的影響最小。綜合圖2和圖3可知,不存在d/λ的取值使得幅度不一致和相位不一致的影響同時最小。當(dāng)需要綜合考慮幅度和相位不一致時,d/λ的取值宜使U1(β)U2(β)最小。U1(β)U2(β)隨陣列直徑波長比的變化如圖4所示。

    圖4 U1(β)U2(β)隨d/λ變化的曲線Fig.4 Curves of U1(β)U2(β)changed with d/λ

    由圖4可知,對于對稱式測角算法,若要求幅相不一致的影響綜合最小,d/λ宜取0.61。

    仰角為0°時,U3(β)和U4(β)隨陣列直徑波長比的變化分別如圖5和圖6所示。

    圖5 U3(β)隨d/λ變化的曲線Fig.5 Curves of U3(β)changed with d/λ

    圖6 U4(β)隨d/λ變化的曲線Fig.6 Curves of U4(β)changed with d/λ

    對于非對稱式測角算法,由圖5可知,d/λ為0.38時幅度不一致影響最小;由圖6可知,d/λ為0.4時相位不一致的影響最小。綜合圖5和圖6可知,對于非對稱式測角方法,同樣不存在d/λ的取值使得幅度不一致和相位不一致的影響同時最小。當(dāng)需要綜合考慮幅度和相位不一致時,d/λ的取值宜使U3(β)U4(β)最小。U3(β)U4(β)隨陣列直徑波長比的變化如圖7所示。

    圖7 U3(β)U4(β)隨d/λ變化的曲線Fig.7 Curves of U3(β)U4(β)changed with d/λ

    由圖7可知,對于非對稱式測角算法,若要求幅相不一致的影響綜合最小,d/λ宜取0.39。

    圖隨d/λ變化的曲線Fig.8 Curves ofUn(β)changed with d/λ

    由圖8可知,綜合考慮幅度和相位不一致分別對非對稱式和對稱式測角算法測角誤差的影響,陣列直徑波長比的最佳取值應(yīng)為0.44。在波長確定的情況下,陣列直徑即可確定。

    需要說明的是,當(dāng)要求對目標(biāo)有一定的仰角覆蓋時,陣列直徑波長比的選擇應(yīng)取得比0.44大。

    為驗證幅相不一致引起的測角誤差公式即式(17)和式(25)的正確性,將公式計算的測角誤差與數(shù)值統(tǒng)計的誤差結(jié)果進(jìn)行對比。

    假設(shè)仰角為0°,陣列直徑波長比為0.44,陣元數(shù)6,其中陣元數(shù)是按照固有誤差不超過1°時所需最少陣元數(shù)的準(zhǔn)則確定,具體可參考文獻(xiàn)[1]。僅存在幅度不一致和僅存在相位不一致時,兩種測角算法測角誤差的公式計算值和數(shù)據(jù)統(tǒng)計值分別隨最大幅度誤差δ和最大相位誤差σ變化的曲線對比如圖9和圖10所示。

    圖9 幅度不一致引起的測角誤差Fig.9 Angle error caused by amplitude inconsistency

    圖10 相位不一致引起的測角誤差Fig.10 Angle error caused by phase inconsistency

    由圖9和圖10可知,幅度和相位不一致分別引起的測角誤差,其公式計算值和誤差統(tǒng)計值基本相同,由此證明了本文所得到的誤差理論公式的正確性。另外,對比圖9和圖10可以發(fā)現(xiàn),對于對稱式測角,相位不一致比幅度不一致的影響更大;對于非對稱式測角,幅度不一致比相位不一致影響更大。

    本文出于為天線設(shè)計留有余量的考慮才對幅相不一致引起的非對稱式測角誤差提出了更高的要求。若實際工程中,陣元間互耦的影響可適當(dāng)控制在較小范圍,那么可以降低對幅相不一致引起的非對稱式測角誤差的要求,此時,陣列直徑波長比的選擇可主要以對稱式測角算法對各通道幅相誤差的要求為參考。

    4結(jié)束語

    聯(lián)合均勻圓陣激勵出的-1階、0階和1階相位模式可實現(xiàn)方位向的無模糊全向測角,但各接收通道存在幅相不一致時會引入測角誤差。本文通過分析幅相不一致與測角誤差的關(guān)系,推導(dǎo)得到了二者之間的解析表達(dá)式并通過數(shù)值仿真驗證其正確性,同時得到了幅相不一致影響下陣列參數(shù)的最優(yōu)取值,為系統(tǒng)陣列參數(shù)的最終確定提供了一定參考。

    參考文獻(xiàn):

    [1]田超, 文樹梁. 米波圓陣?yán)走_(dá)全向測角算法及其性能分析[J]. 系統(tǒng)工程與電子技術(shù), 2013, 35(12): 2507-2512.

    TIAN Chao, WEN Shu-Liang. Omni-Directional Angle Measurement Algorithm of Meter-Wave Circular Array Radar and Its Performance Analysis[J]. Systems Engineering and Electronics, 2013, 35(12): 2507-2512.

    [2]CHEN B X, ZHAO G H, ZHANG S H. Altitude Mea-surement Based on Beam Split and Frequency Diversity in VHF Radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(1): 3-13.

    [3]胡曉琴, 陳建文, 王永良,等. MIMO體制米波圓陣?yán)走_(dá)研究[J]. 國防科技大學(xué)學(xué)報, 2009, 54(1): 52-57.

    HU Xiao-qin, CHEN Jian-wen, WANG Yong-liang,et al. Research on the Meter-Wave Circular Array Radar for MIMO Systems[J]. Journal of National University of Defense Technology, 2009, 54(1): 52-57.

    [4]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.

    [5]PESAVENTOAND M, BOHME J F. Direction of Arrival Estimation in Uniform Circular Arrays Composed of Directional Elements[C]. Proc of Sensor Array and Multichannel Signal Processing Workshop Proceedings, 2002:503-507.

    [6]GOOSSENS R, ROGIER H. A Hybrid UCA-RARE/Root-MUSIC Approach for 2-D Direction of Arrival Estimation in Uniform Circular Arrays in the Presence of Mutual Coupling[J]. IEEE Transactions on Antennas and Propagation,2007, 55(3): 841-849.

    [7]WANG B H, HUI H T, LEONG M S. Decoupled 2D Direction of Arrival Estimation Using Compact Uniform Circular Arrays in the Presence of Elevation-Dependent Mutual Coupling[J]. IEEE Transactions on Antennas and Propagation, 2010, 58(3): 747-755.

    [8]李穎, 黃曉濤, 周智敏. 幅相誤差對均勻圓陣系統(tǒng)性能的影響[J]. 雷達(dá)科學(xué)與技術(shù), 2005, 3(2): 123-127.

    LI Ying, HUANG Xiao-tao, ZHOU Zhi-min. Effect of Amplitude and Phase Errors on the Performance of Uniform Circular Array[J]. Radar Science and Technology, 2005, 3(2): 123-127.

    [9]高火濤, 柯亨玉, 楊子杰,等. 互耦系數(shù)誤差與通道幅相誤差的等效關(guān)系[J]. 武漢大學(xué)學(xué)報:理學(xué)版, 2002, 48(1): 121-124.

    GAO Huo-tao, KE Heng-yu, YANG Zi-jie, et al.The Relationship Between Mutual Coefficient Errors and the Amplitude and Phase Errors of Channels[J]. Wuhan University Journal:Natural Science ed, 2002, 48(1): 121-124.

    [10]LU Zu-kun,JIANG Hong, GAO Ying. Amplitude and Phase Errors Self-Correcting Algorithm Based on the Uniform Circular Array[C]∥Proc of 2nd International Conference on Computer Science and Network Technology, Changchun, 2012: 136-140.

    [11]DU Yong-xing, XI Xiao-li, XI Wen-jing. The Correction of Mutual Coupling and the Amplitude and Phase error Based on Digital Beamforming Receiving Antenna[C]∥Proc of 10th International Symposium on Antennas,Propagation & EM Theory , 2012: 231-233.

    [12]魯祖坤, 高鷹, 肖劍,等. 基于均勻圓陣的幅相誤差自校正算法[J]. 火控雷達(dá)技術(shù), 2013,42(1): 1-5.

    LU Zu-kun, GAO Ying, XIAO Jian, et al.Amplitude and Phase Error Auto-Calibration Algorithm Based on Uniform Circular Array[J]. Fire Control Radar Technology, 2013, 42(1): 1-5.

    [13]陳德莉, 盧煥章, 張聰. 多徑條件下基于WSF的均勻圓陣幅相誤差自校正[J]. 系統(tǒng)仿真學(xué)報, 2008, 20(17): 4563-4565.

    CHEN De-li, LU Huang-zhang, ZHANG Cong.UCA Gain/Phase Self-Calibration in Presence of Multipath Based on Weighted Subspace Fitting[J]. Journal of System Simulation, 2008,20(17): 4563-4565.

    [14]DAVIES D E N. A Transformation Between the Phasing Techniques Required for Linear and Circular Aerial Arrays[C]∥Proc of the Institution of Electrical Engineers, 1965, 112(11): 2041-2045.

    [15]DU K L. Pattern Analysis of Uniform Circular Array[J]. IEEE Transactions on Antennas and Propagation, 2004, 52(4): 1125-1129.

    Error Analysis of Omni-Directional Angle Measurement with Amplitude and Phase Inconsistency between Channels

    TIAN Chao,WEN Shu-liang

    (Beijing Institute of Radio Measurement, Beijing 100854, China)

    Abstract:In ideal scenarios, unambiguous azimuth angle is available by using the phase mode -1, 0 and 1 based on the omni-directional VHF radar with uniform circular array. However, in engineering, the amplitude and phase character of each receiving channel is usually different. This inconsistency may cause some errors of the measured angle. Moreover, the errors can’t be compensated and can only be decreased by calibrating the inconsistency or choosing the array parameters. For this purpose, the specific relationship between angle error and amplitude and phase inconsistency is deduced theoretically so that a theoretical reference is provided for the optimal selection of array parameter in order to reduce the angle error. Simulation results demonstrate the correctness of the deduced formulas and the optimal diameter of the array is also given.

    Key words:VHF radar; uniform circular array; omni-directional angle measurement; phase mode; amplitude and phase inconsistency; angle error analysis

    中圖分類號:TN953+.5;TP391.9

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1009-086X(2015)-06-0104-08

    doi:10.3969/j.issn.1009-086x.2015.06.018

    通信地址:100854北京市142信箱203分箱1號E-mail:qctchao87@126.com

    作者簡介:田超(1987-),男,湖北黃岡人。博士生,主要研究方向為雷達(dá)總體技術(shù)。

    *收稿日期:2014-11-03;修回日期:2015-12-08

    欧美区成人在线视频| 亚洲精品粉嫩美女一区| 看免费成人av毛片| 精品久久久久久久久久免费视频| 国产精品美女特级片免费视频播放器| 日日撸夜夜添| 免费看光身美女| 老司机福利观看| 国产成人福利小说| 婷婷精品国产亚洲av在线| 不卡视频在线观看欧美| 国产91av在线免费观看| 俺也久久电影网| 成人特级av手机在线观看| 男女做爰动态图高潮gif福利片| 一区二区三区高清视频在线| 婷婷亚洲欧美| 国产高清三级在线| 国产亚洲欧美98| 99riav亚洲国产免费| 免费观看人在逋| 成年女人永久免费观看视频| 日韩精品有码人妻一区| 伊人久久精品亚洲午夜| 级片在线观看| 国产单亲对白刺激| 俺也久久电影网| 深爱激情五月婷婷| 免费高清视频大片| 九九久久精品国产亚洲av麻豆| 午夜福利18| 夜夜看夜夜爽夜夜摸| or卡值多少钱| 国产蜜桃级精品一区二区三区| 蜜臀久久99精品久久宅男| 国产精品乱码一区二三区的特点| 欧美潮喷喷水| 国产精品永久免费网站| 在线免费观看不下载黄p国产| 麻豆成人午夜福利视频| 亚洲av第一区精品v没综合| 久久久久久久午夜电影| 精品午夜福利在线看| 婷婷精品国产亚洲av在线| 国产熟女欧美一区二区| 大香蕉久久网| 成人高潮视频无遮挡免费网站| 亚洲丝袜综合中文字幕| 久久精品国产亚洲av天美| 狂野欧美白嫩少妇大欣赏| 99久国产av精品国产电影| 久久九九热精品免费| 女人十人毛片免费观看3o分钟| 午夜亚洲福利在线播放| 九九热线精品视视频播放| 国产蜜桃级精品一区二区三区| 蜜臀久久99精品久久宅男| 国产一级毛片七仙女欲春2| 免费av不卡在线播放| 久久久久国产精品人妻aⅴ院| 在线国产一区二区在线| 最近最新中文字幕大全电影3| 人妻久久中文字幕网| 欧美中文日本在线观看视频| 别揉我奶头~嗯~啊~动态视频| 俺也久久电影网| 精品免费久久久久久久清纯| 欧美性感艳星| 岛国在线免费视频观看| 91久久精品国产一区二区三区| 欧美日韩精品成人综合77777| 12—13女人毛片做爰片一| 白带黄色成豆腐渣| 国产精品亚洲美女久久久| 成人精品一区二区免费| 可以在线观看毛片的网站| 一个人免费在线观看电影| av福利片在线观看| 精品少妇黑人巨大在线播放 | 国产精品一区二区三区四区免费观看 | 搡老岳熟女国产| 欧美+亚洲+日韩+国产| 亚洲av中文av极速乱| 两个人的视频大全免费| 永久网站在线| 99视频精品全部免费 在线| 熟女电影av网| 日本黄大片高清| 天天躁夜夜躁狠狠久久av| 成年女人看的毛片在线观看| 精品熟女少妇av免费看| 麻豆久久精品国产亚洲av| 在线a可以看的网站| 午夜福利在线观看免费完整高清在 | 亚洲成人久久爱视频| 亚洲欧美日韩卡通动漫| av在线播放精品| 毛片一级片免费看久久久久| 嫩草影院入口| 美女cb高潮喷水在线观看| 欧美色视频一区免费| 在线免费十八禁| 国产视频内射| 极品教师在线视频| av黄色大香蕉| av在线老鸭窝| 精品久久久噜噜| 久久亚洲国产成人精品v| 亚洲欧美日韩东京热| 夜夜看夜夜爽夜夜摸| 22中文网久久字幕| 国产精品无大码| 91av网一区二区| 村上凉子中文字幕在线| 久久国内精品自在自线图片| 国产不卡一卡二| 精品人妻视频免费看| 久久热精品热| 婷婷精品国产亚洲av在线| 亚洲精品成人久久久久久| 久久婷婷人人爽人人干人人爱| 久久精品国产清高在天天线| 国产成人91sexporn| 黑人高潮一二区| 国产高清视频在线观看网站| 日韩欧美三级三区| 波多野结衣高清作品| 日韩三级伦理在线观看| 久久精品影院6| 久久人人爽人人爽人人片va| a级毛片免费高清观看在线播放| a级毛片免费高清观看在线播放| 乱人视频在线观看| 在线观看美女被高潮喷水网站| 九九热线精品视视频播放| 亚洲最大成人av| 乱人视频在线观看| 久久韩国三级中文字幕| 成人无遮挡网站| 午夜免费激情av| 别揉我奶头~嗯~啊~动态视频| 麻豆av噜噜一区二区三区| 干丝袜人妻中文字幕| 国产精品久久久久久久久免| 成人亚洲精品av一区二区| 免费看日本二区| 老司机午夜福利在线观看视频| 欧美中文日本在线观看视频| aaaaa片日本免费| 成人鲁丝片一二三区免费| 蜜桃久久精品国产亚洲av| 乱码一卡2卡4卡精品| 精品一区二区三区av网在线观看| 中国美白少妇内射xxxbb| 人人妻人人看人人澡| 人妻丰满熟妇av一区二区三区| 亚洲av中文字字幕乱码综合| 日韩成人伦理影院| 可以在线观看毛片的网站| 色5月婷婷丁香| 亚洲国产精品合色在线| 在线天堂最新版资源| 成人av一区二区三区在线看| 亚洲自拍偷在线| 日韩欧美一区二区三区在线观看| 成年女人永久免费观看视频| 少妇的逼好多水| 久久久久久大精品| 超碰av人人做人人爽久久| 久久天躁狠狠躁夜夜2o2o| 亚洲最大成人手机在线| 中文字幕精品亚洲无线码一区| 夜夜夜夜夜久久久久| 国产久久久一区二区三区| 国产乱人偷精品视频| 搡老熟女国产l中国老女人| 又黄又爽又免费观看的视频| 三级男女做爰猛烈吃奶摸视频| 99久久久亚洲精品蜜臀av| 最好的美女福利视频网| 久久精品国产亚洲av香蕉五月| 国产精品人妻久久久影院| 国产亚洲91精品色在线| 亚州av有码| av在线天堂中文字幕| 中国美女看黄片| 久久人人爽人人片av| 九九爱精品视频在线观看| 可以在线观看毛片的网站| 国产成人freesex在线 | 亚洲无线在线观看| 日韩强制内射视频| 亚洲av熟女| 非洲黑人性xxxx精品又粗又长| 日本在线视频免费播放| 国产伦精品一区二区三区四那| 黄色日韩在线| 波多野结衣高清无吗| 久久精品国产清高在天天线| 国产一区二区在线av高清观看| 最近中文字幕高清免费大全6| 老司机影院成人| 亚洲人成网站在线播| 看片在线看免费视频| 小说图片视频综合网站| 中文字幕av成人在线电影| 色综合亚洲欧美另类图片| 男女之事视频高清在线观看| 亚洲国产精品成人久久小说 | 精品人妻熟女av久视频| 成年av动漫网址| 中国美白少妇内射xxxbb| 精品一区二区三区人妻视频| eeuss影院久久| 丰满乱子伦码专区| 毛片女人毛片| 别揉我奶头 嗯啊视频| 在线播放无遮挡| 麻豆久久精品国产亚洲av| 国产成人91sexporn| 又黄又爽又刺激的免费视频.| av专区在线播放| 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 国产一区二区三区在线臀色熟女| 久久久久九九精品影院| 久久99热6这里只有精品| av专区在线播放| 国产老妇女一区| 深夜a级毛片| 久久久久久久午夜电影| 嫩草影院入口| 99久国产av精品| 日日摸夜夜添夜夜添小说| 波野结衣二区三区在线| 成人三级黄色视频| 成年女人永久免费观看视频| 国产亚洲精品久久久com| 亚洲美女搞黄在线观看 | 亚洲av成人av| 午夜老司机福利剧场| 如何舔出高潮| 亚洲,欧美,日韩| 此物有八面人人有两片| 国产精品人妻久久久影院| 亚洲婷婷狠狠爱综合网| 成年免费大片在线观看| 直男gayav资源| 国产在线男女| 精品一区二区三区人妻视频| 欧美又色又爽又黄视频| 美女黄网站色视频| 欧美日韩乱码在线| 日韩精品有码人妻一区| 亚洲av成人av| 色播亚洲综合网| 久久午夜福利片| 午夜福利视频1000在线观看| 欧美极品一区二区三区四区| 中文字幕人妻熟人妻熟丝袜美| 乱码一卡2卡4卡精品| 亚洲av二区三区四区| 我的老师免费观看完整版| 日日啪夜夜撸| 亚洲国产精品sss在线观看| 久久中文看片网| 国产亚洲91精品色在线| 搡女人真爽免费视频火全软件 | 99国产极品粉嫩在线观看| 精品人妻偷拍中文字幕| 22中文网久久字幕| 亚州av有码| 日韩强制内射视频| 久久久精品94久久精品| 嫩草影院入口| 免费搜索国产男女视频| 亚洲av免费高清在线观看| 直男gayav资源| 人人妻人人看人人澡| 国产老妇女一区| 日本 av在线| 免费电影在线观看免费观看| 深夜精品福利| 国产精品久久视频播放| 久久人人爽人人片av| 成人欧美大片| 亚洲无线观看免费| 日韩精品中文字幕看吧| 看十八女毛片水多多多| 村上凉子中文字幕在线| 亚洲欧美日韩高清在线视频| 日韩av不卡免费在线播放| 亚洲精品亚洲一区二区| 成人午夜高清在线视频| 日本色播在线视频| 午夜影院日韩av| 99热只有精品国产| 一级毛片我不卡| 乱系列少妇在线播放| 成人性生交大片免费视频hd| 久久午夜亚洲精品久久| 国产乱人偷精品视频| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区免费观看 | 国产精品人妻久久久久久| 国产av麻豆久久久久久久| 美女被艹到高潮喷水动态| 国产精品不卡视频一区二区| 国产av麻豆久久久久久久| 伊人久久精品亚洲午夜| 亚洲成人精品中文字幕电影| 日韩制服骚丝袜av| or卡值多少钱| 国产视频内射| 又黄又爽又刺激的免费视频.| 国产毛片a区久久久久| 男女下面进入的视频免费午夜| 深夜精品福利| 免费搜索国产男女视频| 免费黄网站久久成人精品| 老女人水多毛片| 亚洲精华国产精华液的使用体验 | 欧美高清性xxxxhd video| 国产av麻豆久久久久久久| 亚洲国产精品合色在线| 国产精品一区二区三区四区免费观看 | 男人的好看免费观看在线视频| 成人亚洲精品av一区二区| 欧美潮喷喷水| 国产精品久久久久久久电影| 亚洲自偷自拍三级| 狠狠狠狠99中文字幕| 观看美女的网站| 麻豆精品久久久久久蜜桃| 中文字幕免费在线视频6| 91狼人影院| 精华霜和精华液先用哪个| 精品久久久久久久久av| 欧美成人a在线观看| 校园春色视频在线观看| 啦啦啦观看免费观看视频高清| 国产精品一区二区三区四区免费观看 | 最近手机中文字幕大全| 亚洲一区高清亚洲精品| 精品国内亚洲2022精品成人| 国产爱豆传媒在线观看| avwww免费| 国产精品精品国产色婷婷| 成年免费大片在线观看| 亚洲不卡免费看| 久久亚洲国产成人精品v| 日本免费一区二区三区高清不卡| 久久久精品94久久精品| 国内精品宾馆在线| 内地一区二区视频在线| av天堂中文字幕网| 亚洲成人中文字幕在线播放| 毛片一级片免费看久久久久| 国产精品伦人一区二区| 午夜免费男女啪啪视频观看 | 人人妻人人澡欧美一区二区| 成人国产麻豆网| 免费搜索国产男女视频| 国产视频内射| 国产欧美日韩精品一区二区| 日韩精品青青久久久久久| 免费电影在线观看免费观看| 免费搜索国产男女视频| 亚洲图色成人| 日韩国内少妇激情av| 日本与韩国留学比较| 俄罗斯特黄特色一大片| 日韩亚洲欧美综合| 大又大粗又爽又黄少妇毛片口| 精品人妻偷拍中文字幕| 深夜精品福利| 久久久久九九精品影院| 最近最新中文字幕大全电影3| 久久久久国产网址| 亚洲欧美精品自产自拍| 国产精品免费一区二区三区在线| 久久午夜亚洲精品久久| 一级av片app| 久久精品国产亚洲av涩爱 | 久久久久免费精品人妻一区二区| 在线观看66精品国产| 国内揄拍国产精品人妻在线| 噜噜噜噜噜久久久久久91| 性欧美人与动物交配| 91在线观看av| 成年版毛片免费区| 国内揄拍国产精品人妻在线| 综合色丁香网| 亚洲欧美日韩高清专用| 国产午夜精品论理片| 午夜影院日韩av| 亚洲七黄色美女视频| 精品日产1卡2卡| 97超碰精品成人国产| 最新中文字幕久久久久| 久久6这里有精品| 中文字幕久久专区| 亚洲成人av在线免费| av免费在线看不卡| 一本一本综合久久| 国产久久久一区二区三区| 天天一区二区日本电影三级| 国产高清激情床上av| 国产男人的电影天堂91| 国产一级毛片七仙女欲春2| 搡老熟女国产l中国老女人| 在线a可以看的网站| 少妇的逼好多水| av女优亚洲男人天堂| 亚洲不卡免费看| 日韩国内少妇激情av| 午夜影院日韩av| 高清毛片免费观看视频网站| 国产91av在线免费观看| 亚洲人成网站高清观看| 日韩欧美 国产精品| 黑人高潮一二区| 日本一二三区视频观看| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 69av精品久久久久久| 欧美另类亚洲清纯唯美| 又黄又爽又免费观看的视频| 我的女老师完整版在线观看| 国产亚洲精品久久久com| 国产在视频线在精品| 嫩草影视91久久| 亚洲熟妇熟女久久| 午夜爱爱视频在线播放| 天天一区二区日本电影三级| 国产真实伦视频高清在线观看| a级毛片a级免费在线| 一卡2卡三卡四卡精品乱码亚洲| 黄色欧美视频在线观看| 老熟妇仑乱视频hdxx| 成人漫画全彩无遮挡| 国产精品嫩草影院av在线观看| 日本一本二区三区精品| 自拍偷自拍亚洲精品老妇| 国产一区二区在线av高清观看| 亚洲国产精品成人久久小说 | 欧美性猛交╳xxx乱大交人| 美女cb高潮喷水在线观看| 欧美高清成人免费视频www| 国产精品综合久久久久久久免费| 十八禁网站免费在线| 女人被狂操c到高潮| 国产午夜精品论理片| 亚洲精品色激情综合| 日本一二三区视频观看| 一卡2卡三卡四卡精品乱码亚洲| 国产一级毛片七仙女欲春2| 99久久精品热视频| 美女xxoo啪啪120秒动态图| 高清毛片免费看| 久久久久国产精品人妻aⅴ院| 久久中文看片网| 精品欧美国产一区二区三| 欧美高清性xxxxhd video| 成人高潮视频无遮挡免费网站| 久久久久国内视频| 日本欧美国产在线视频| 亚洲美女搞黄在线观看 | av国产免费在线观看| 日韩强制内射视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区三区四区免费观看 | 日韩欧美 国产精品| 免费大片18禁| 乱系列少妇在线播放| 国产伦精品一区二区三区四那| 最新中文字幕久久久久| av卡一久久| 国产乱人偷精品视频| 国产午夜精品久久久久久一区二区三区 | 日韩av不卡免费在线播放| 麻豆国产97在线/欧美| 国产久久久一区二区三区| 伊人久久精品亚洲午夜| av中文乱码字幕在线| 午夜福利视频1000在线观看| 国产av麻豆久久久久久久| 国产老妇女一区| 国产黄色小视频在线观看| 成人鲁丝片一二三区免费| 久久久久久久久久成人| 哪里可以看免费的av片| 午夜精品在线福利| av专区在线播放| 欧美bdsm另类| 六月丁香七月| 亚洲欧美日韩高清在线视频| 国内精品美女久久久久久| 一级毛片久久久久久久久女| 国产色爽女视频免费观看| 97碰自拍视频| videossex国产| 国内精品久久久久精免费| 日本免费一区二区三区高清不卡| 激情 狠狠 欧美| 精品午夜福利视频在线观看一区| 日韩人妻高清精品专区| 最近中文字幕高清免费大全6| 天堂av国产一区二区熟女人妻| 日韩人妻高清精品专区| 日韩av在线大香蕉| 97超碰精品成人国产| 一a级毛片在线观看| 18禁在线无遮挡免费观看视频 | 丰满人妻一区二区三区视频av| 亚洲精品乱码久久久v下载方式| 国产男靠女视频免费网站| 人妻少妇偷人精品九色| 亚洲第一区二区三区不卡| videossex国产| 熟妇人妻久久中文字幕3abv| 久久久久久伊人网av| 蜜桃亚洲精品一区二区三区| 亚洲精华国产精华液的使用体验 | 免费观看在线日韩| 一夜夜www| 成人亚洲精品av一区二区| 成年女人看的毛片在线观看| 99久久无色码亚洲精品果冻| 大又大粗又爽又黄少妇毛片口| 三级男女做爰猛烈吃奶摸视频| 日日干狠狠操夜夜爽| 99热这里只有是精品50| 午夜a级毛片| 国产美女午夜福利| 国产免费男女视频| a级毛片a级免费在线| 精品99又大又爽又粗少妇毛片| 国产精品久久久久久亚洲av鲁大| 日韩欧美免费精品| 淫妇啪啪啪对白视频| 两个人的视频大全免费| 偷拍熟女少妇极品色| av专区在线播放| 欧美日本亚洲视频在线播放| 亚洲熟妇中文字幕五十中出| 夜夜夜夜夜久久久久| 又爽又黄a免费视频| 一个人看视频在线观看www免费| www日本黄色视频网| 毛片女人毛片| 永久网站在线| 嫩草影院新地址| 在线播放无遮挡| 国产国拍精品亚洲av在线观看| 午夜久久久久精精品| av在线观看视频网站免费| 亚洲va在线va天堂va国产| 波野结衣二区三区在线| 嫩草影院精品99| 亚洲av中文字字幕乱码综合| 日日撸夜夜添| 高清毛片免费看| 麻豆久久精品国产亚洲av| 嫩草影院入口| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 中文字幕免费在线视频6| 亚洲精品在线观看二区| 成人av在线播放网站| h日本视频在线播放| 色哟哟·www| 99久久成人亚洲精品观看| 国产午夜精品久久久久久一区二区三区 | 俺也久久电影网| 婷婷精品国产亚洲av| 成人特级av手机在线观看| 免费观看精品视频网站| 日韩制服骚丝袜av| 高清毛片免费看| 日韩欧美精品v在线| 精品午夜福利视频在线观看一区| 国产乱人视频| 熟女人妻精品中文字幕| 久久精品国产亚洲av涩爱 | 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 99热精品在线国产| 有码 亚洲区| 午夜精品一区二区三区免费看| 美女黄网站色视频| 性欧美人与动物交配| 人人妻人人澡人人爽人人夜夜 | 免费在线观看成人毛片| 亚洲精品在线观看二区| 能在线免费观看的黄片| 亚洲精品亚洲一区二区| 亚洲国产精品成人综合色| 特级一级黄色大片| АⅤ资源中文在线天堂| 亚洲精品456在线播放app| 无遮挡黄片免费观看| 少妇熟女欧美另类| 日本与韩国留学比较| 男插女下体视频免费在线播放| av.在线天堂| 日本爱情动作片www.在线观看 | 欧美日韩在线观看h| 97热精品久久久久久| 亚洲中文字幕一区二区三区有码在线看| 免费大片18禁| 又黄又爽又免费观看的视频| 亚洲精品国产成人久久av| 欧美人与善性xxx| 露出奶头的视频|