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

    基于序列線性規(guī)劃的雷達(dá)低峰均比估計(jì)波形設(shè)計(jì)

    2018-10-15 09:42:34郝天鐸孫從易
    關(guān)鍵詞:優(yōu)化

    郝天鐸, 崔 琛, 龔 陽(yáng), 孫從易

    (1. 國(guó)防科技大學(xué)電子對(duì)抗學(xué)院, 安徽 合肥 230037; 2. 中國(guó)人民解放軍96630部隊(duì), 北京 102206)

    0 引 言

    波形設(shè)計(jì)是雷達(dá)的關(guān)鍵技術(shù)之一。在進(jìn)行面向參數(shù)估計(jì)的波形設(shè)計(jì)時(shí),現(xiàn)有大部分方法均在頻域[1-3]進(jìn)行研究,很難直接對(duì)波形包絡(luò)加以約束,從而無(wú)法兼顧放大器的非線性特性導(dǎo)致發(fā)射波形失真,無(wú)法發(fā)揮雷達(dá)發(fā)射機(jī)的最大效能,所以實(shí)際系統(tǒng)中還必須考慮波形的低峰均比(peak-to-average power ratio, PAR)或者恒模約束[4]。然而,恒模波形雖然是一種理想的發(fā)射波形,但其自由度較低,與恒模波形相比,能夠同時(shí)兼顧估計(jì)性能和雷達(dá)功率放大器的非線性特性[5-6]。因此,基于時(shí)域的低PAR估計(jì)波形更加符合實(shí)際工程需求。

    現(xiàn)階段,只有少量文獻(xiàn)對(duì)該問(wèn)題進(jìn)行了研究[4-6]。文獻(xiàn)[4]基于最小均方誤差(minimum mean square error,MMSE)準(zhǔn)則[7],在PAR約束下對(duì)信道參數(shù)進(jìn)行了估計(jì),但其并未考慮信號(hào)相關(guān)雜波的影響。文獻(xiàn)[8]在信號(hào)模型中加入了信號(hào)相關(guān)雜波,基于最小克拉美羅界(Cramer-Rao bound,CRB)準(zhǔn)則對(duì)目標(biāo)散射系數(shù)進(jìn)行估計(jì),但其假定發(fā)射波形和雜波協(xié)方差矩陣的卷積項(xiàng)是固定不變的,而實(shí)際中由于發(fā)射波形是不斷迭代更新的,所以卷積項(xiàng)也是不斷變化的。文獻(xiàn)[9]基于最大互信息(mutual information,MI)準(zhǔn)則[9],在頻域模型加入信號(hào)相關(guān)雜波,采用時(shí)頻域轉(zhuǎn)換的方法,將頻域所得能量約束下的最優(yōu)波形轉(zhuǎn)化到時(shí)域進(jìn)行PAR波形設(shè)計(jì),由于頻域波形不含相位信息,從頻域變換到時(shí)域后會(huì)帶來(lái)估計(jì)性能的下降。在信號(hào)相關(guān)雜波背景下,以上文獻(xiàn)均無(wú)法直接在時(shí)域合成低PAR估計(jì)波形,這主要是由于所求優(yōu)化問(wèn)題是一個(gè)較為復(fù)雜的非凸問(wèn)題,很難將其轉(zhuǎn)化為凸問(wèn)題進(jìn)行求解。

    針對(duì)上述問(wèn)題,本文引入序列線性規(guī)劃(sequence linear programming, SLP)的思想,用線性形式擬合原優(yōu)化問(wèn)題的目標(biāo)函數(shù)和約束條件,實(shí)現(xiàn)了非凸問(wèn)題向凸問(wèn)題的轉(zhuǎn)換,在信號(hào)相關(guān)雜波背景下直接在時(shí)域合成了低PAR估計(jì)波形。仿真實(shí)驗(yàn)驗(yàn)證了本文理論分析的正確性和算法的有效性。

    1 信號(hào)模型

    本文主要考慮相關(guān)雜波條件下針對(duì)擴(kuò)展目標(biāo)估計(jì)的波形設(shè)計(jì)問(wèn)題。目標(biāo)的散射特性用目標(biāo)沖激響應(yīng)(target impulse response, TIR)表示[10],相關(guān)雜波用雜波沖激響應(yīng)(clutter impulse response, CIR)表示[11],主要在時(shí)域?qū)﹄x散的時(shí)間信號(hào)進(jìn)行分析,信號(hào)模型如圖1所示。

    圖1 相關(guān)雜波下的信號(hào)模型Fig.1 Signal model

    圖1中,s∈CNs×1代表發(fā)射波形,其長(zhǎng)度為Ns;t∈CNt×1和c∈CNc×1分別代表TIR和CIR,根據(jù)文獻(xiàn)[12]的雷達(dá)信號(hào)模型,令TIR和CIR長(zhǎng)度相同,即Nt=Nc;n∈CNn×1代表噪聲和干擾的總和,Nn=Ns+Nc-1,x代表來(lái)自目標(biāo)和環(huán)境的回波,其長(zhǎng)度Nx=Nn。該模型可表示為

    x=t*s+c*s+n=Ts+Cs+n=

    St+Sc+n=st+sc+n

    (1)

    (2)

    本文假定噪聲向量n服從均值為0,協(xié)方差矩陣為單位陣的復(fù)高斯分布,同時(shí)假定TIR和CIR均為復(fù)高斯隨機(jī)向量,其中,t~CN(0Nt,Rt),c~CN(0Nc,Rc),Rt∈CNt×Nt,Rc∈CNc×Nc。

    2 波形設(shè)計(jì)方法

    2.1 基于互信息方法的問(wèn)題描述

    本文以回波和目標(biāo)之間最大化互信息為優(yōu)化準(zhǔn)則進(jìn)行波形設(shè)計(jì),互信息越大,估計(jì)越精確??紤]前述對(duì)噪聲、TIR和CIR的假設(shè),并假定三者相互獨(dú)立,則目標(biāo)t和回波信號(hào)x的互信息可表示如下[13]

    I(x;t|S)=h(x|S)-h(x|t,S)

    (3)

    式中,h(x|S)代表發(fā)射波形Toeplitz矩陣S已知時(shí)回波信號(hào)的熵;h(x|t,S)代表S和TIR均已知時(shí)回波信號(hào)的熵。當(dāng)S已知時(shí),TIR和回波信號(hào)服從聯(lián)合高斯分布,且滿足

    (4)

    令Rx=SRtSH+SRcSH+Rn,有

    (5)

    則可得

    ln det(SRtSH+SRcSH+Rn)+Nxln π+Nx

    (6)

    ln det(SRcSH+Rn)+Nxln π+Nx

    (7)

    把式(6)和式(7)代入式(3)中,可得

    I(x,t|s)=ln det[INx+SRtSH(SRcSH+Rn)-1]

    (8)

    應(yīng)用矩陣求逆引理[14],于是式(8)變?yōu)?/p>

    I(x;t|S)=ln det{INt-Rt[AB-1A-A]}

    (9)

    (10)

    假設(shè)噪聲、TIR和CIR均為廣義寬平穩(wěn)隨機(jī)過(guò)程,則Rt和Rc[15]可分別表示為

    Rt=UΣtUH

    (11)

    Rc=UΣcUH

    (12)

    式中,Σt=diag(λt,1,…,λt,Nt),Σc=diag(λc,1,…,λc,Nc),它們的對(duì)角元素代表目標(biāo)和雜波功率譜密度的采樣,酉矩陣U可表示為

    ?k,l∈(1,Nt)

    (13)

    同理,有

    Rn=VΣnVH

    (14)

    式中,Σn=diag(λn,1,…,λn,Nn)代表對(duì)噪聲功率譜密度的采樣;V也為酉矩陣。把式(11)、式(12)和式(14)代入式(10)中,可得

    ln det{INt+Σt[Σc+(ZHZ)-1]-1}

    (15)

    (16)

    式(16)也可通過(guò)文獻(xiàn)[17]的引理2得到,其中Λs=diag(λs,1,…,λs,Nt),則互信息表達(dá)式可變?yōu)?/p>

    (17)

    至此,對(duì)時(shí)域上TIR和回波信號(hào)的互信息表達(dá)式進(jìn)行了推導(dǎo)簡(jiǎn)化,可以看出,式(17)與頻域上的互信息表達(dá)式較為相似。接下來(lái)將引入能量和PAR約束對(duì)時(shí)域發(fā)射波形進(jìn)行求解。

    2.2 低PAR估計(jì)波形設(shè)計(jì)

    不失一般性,設(shè)發(fā)射波形的最大能量Ps=Ns,則PAR可定義為

    (18)

    式中,sk為波形向量s的第k個(gè)采樣。在加上PAR約束后,要求每一個(gè)發(fā)射波形元素sk的能量都小于某一固定的PAR值,即

    |sk|2≤μ,μ∈[1,Ns]

    (19)

    (20)

    (21)

    (22)

    問(wèn)題P′的不等式約束雖然變?yōu)橐粋€(gè)凸集,但目標(biāo)函數(shù)卻是非凸的,因此對(duì)該問(wèn)題求解較為困難。此時(shí),可引入SLP的思想對(duì)其進(jìn)行求解。

    2.3 引入SLP方法的問(wèn)題求解

    SLP是一種將非線性規(guī)劃問(wèn)題轉(zhuǎn)化為線性規(guī)劃問(wèn)題的研究方法。在實(shí)際優(yōu)化過(guò)程中,由于其理論簡(jiǎn)單、易于實(shí)現(xiàn)的特點(diǎn)被廣泛應(yīng)用,在凸規(guī)劃問(wèn)題中用線性形式擬合目標(biāo)函數(shù)和約束不等式方程方面均有很好的表現(xiàn)[18]。SLP的主要思想是通過(guò)泰勒公式,在初始值x(0)處對(duì)非線性目標(biāo)函數(shù)和約束條件進(jìn)行一階展開(kāi),略去二階和二階以上的高階項(xiàng),實(shí)現(xiàn)非線性問(wèn)題中目標(biāo)函數(shù)和約束條件的線性化處理,并通過(guò)求解該線性規(guī)劃問(wèn)題逼近原始問(wèn)題的解[19]。設(shè)非線性函數(shù)為f(x),在x(0)處的泰勒展開(kāi)式為

    f(x)≈f(x(0))+f′(x(0))(Δx)

    (23)

    在該問(wèn)題的約束條件下,通過(guò)對(duì)式(23)進(jìn)行求解,可以得到此方程的最優(yōu)解x(1),接著將目標(biāo)函數(shù)在x(1)處展開(kāi)

    f(x)≈f(x(1))+f′(x(1))(Δx)

    (24)

    以此類推,可得SLP的通式為

    f(x(n))≈f(x(n-1))+f′(x(n-1))(Δx)

    (25)

    對(duì)式(25)不斷求解更新,直到滿足收斂條件即可得到最優(yōu)解x(n)。

    然而,SLP有一定的適用范圍,即新的設(shè)計(jì)點(diǎn)與原設(shè)計(jì)點(diǎn)的距離Δx不可過(guò)遠(yuǎn),否則高階項(xiàng)不可忽略。當(dāng)Δx較小時(shí),近似精度高,但優(yōu)化速度慢;當(dāng)Δx較大時(shí),優(yōu)化速度快,但近似精度低,將會(huì)使近似函數(shù)偏離原始函數(shù)。因此,選擇合適的Δx尤為重要,在后面的方法設(shè)計(jì)中,對(duì)Δx的設(shè)定范圍要盡量合理,使得近似函數(shù)在靠近原函數(shù)的同時(shí)盡可能收斂到全局最優(yōu)值。綜上所述,由于序列線性規(guī)劃算法對(duì)于求解工程中的非線性問(wèn)題有著很好的效果,本文采用SLP的思想,對(duì)目標(biāo)函數(shù)進(jìn)行一階展開(kāi),將原始問(wèn)題轉(zhuǎn)化為一個(gè)凸問(wèn)題進(jìn)行求解。

    2.3.1 初始值求解

    首先,進(jìn)行初始值的求解。為了求得一個(gè)較好的初始極值,可以將約束條件中的不等式約束變?yōu)榈仁郊s束,采用拉格朗日法進(jìn)行求解。然而,第二個(gè)約束條件為PAR約束,若μ>1,則變?yōu)榈仁郊s束后會(huì)使發(fā)射波形的總能量超過(guò)限定的最大能量,所以在這里可以先略去第二個(gè)不等式約束,只將第一個(gè)不等式約束條件變?yōu)榈仁郊s束,優(yōu)化問(wèn)題變?yōu)?/p>

    (26)

    (27)

    (28)

    2.3.2 目標(biāo)函數(shù)優(yōu)化

    (29)

    (30)

    2.3.3 波形向量的逼近

    在問(wèn)題P2的求解過(guò)程中,并未對(duì)S進(jìn)行Toeplitz矩陣結(jié)構(gòu)(見(jiàn)式(2))的約束,因此優(yōu)化波形矩陣Sopt并不一定滿足信號(hào)模型中的Toeplitz矩陣結(jié)構(gòu),還需要對(duì)其進(jìn)行Toeplitz矩陣逼近后才可求得發(fā)射波形的向量解。這里采用最小F-范數(shù)的方法來(lái)逼近波形向量,有

    (31)

    d-sHg-gHs

    (32)

    (33)

    綜上,為了得到滿足低PAR要求的發(fā)射波形向量,需引入SLP的思想,將非線性優(yōu)化問(wèn)題轉(zhuǎn)化為線性問(wèn)題,快速求解出近似的波形向量,便于實(shí)際工程應(yīng)用。設(shè)門(mén)限值為τ,最大迭代次數(shù)為κ,優(yōu)化算法的步驟可總結(jié)如下。

    3 算法性能分析

    3.1 算法收斂性

    在步驟2迭代的過(guò)程中,記第k次迭代得到的低PAR波形矩陣為S(k),記問(wèn)題P2的目標(biāo)函數(shù)值為I2(S(k)),記原始問(wèn)題P的目標(biāo)函數(shù)值為I(S(k))。通常對(duì)于SLP問(wèn)題而言,I2(S)是單調(diào)收斂的。由于問(wèn)題P2是問(wèn)題P的近似,所以不能保證I(S)也是單調(diào)的。但是,隨著k的增加,S(k)會(huì)逐漸趨近于原始問(wèn)題的最優(yōu)解,這就保證了I(S)最終會(huì)收斂趨近于最優(yōu)值。因此,雖然不能確保I(S)每一步都是單調(diào)遞增的,但其總體趨勢(shì)可以認(rèn)為是收斂的。由于要證明I(S)的收斂性較為困難,因此,通過(guò)仿真實(shí)驗(yàn)對(duì)其收斂性進(jìn)行了驗(yàn)證。

    此外,在步驟3的逼近過(guò)程中,主要針對(duì)優(yōu)化結(jié)果進(jìn)行逼近,而不是針對(duì)目標(biāo)函數(shù)逼近。因此,對(duì)于非線性程度較高的目標(biāo)函數(shù),逼近結(jié)果不一定能夠保證目標(biāo)函數(shù)取到最優(yōu)。以圖2為例,圖中是一條隨機(jī)的非線性曲線,若sopt處于圖中目標(biāo)函數(shù)的凹區(qū)間,雖然其針對(duì)結(jié)果的逼近誤差比圖中實(shí)線框內(nèi)的波形向量要小,但由于約束條件和目標(biāo)函數(shù)非線性的原因,使得目標(biāo)函數(shù)值反而變小,導(dǎo)致F-范數(shù)逼近方法效果變差。但是在一般情況下,該逼近方法不失為一種波形優(yōu)化途徑,可認(rèn)為其是一種次優(yōu)化方法。

    另外,為得到較好的波形向量,在問(wèn)題P2中直接加入對(duì)S的Toeplitz約束條件,通過(guò)智能搜索算法也可以得到較為理想的優(yōu)化結(jié)果,但這種方法將使得優(yōu)化問(wèn)題的求解耗費(fèi)更長(zhǎng)的時(shí)間。

    圖2 非線性程度較高的目標(biāo)函數(shù)示例Fig.2 Example of objective function with high degree of non-linearity

    3.2 算法復(fù)雜度分析

    4 仿真實(shí)驗(yàn)與分析

    圖3 目標(biāo)和雜波功率譜密度采樣Fig.3 Power spectral density of target and clutter

    4.1 實(shí)驗(yàn)1:算法有效性驗(yàn)證

    圖4 問(wèn)題P和P2的收斂性比較Fig.4 Convergence comparison of P and P2

    圖5 不同波形估計(jì)性能對(duì)比Fig.5 Estimation performance comparison of different waveforms

    從圖5可以看出,采用逼近方法得到的滿足Toeplitz矩陣結(jié)構(gòu)的發(fā)射波形的估計(jì)性能比最優(yōu)波形略差,但其估計(jì)性能仍然優(yōu)于文獻(xiàn)[6]和常用波形,這是由于文獻(xiàn)[6]的最優(yōu)頻域波形不含有相位信息,從頻域變換到時(shí)域時(shí)會(huì)帶來(lái)估計(jì)性能的下降。此外,目標(biāo)與回波之間的MI隨著能量的增大而增大,隨著CNR的增大而減小,在低CNR時(shí)(CNR≤20)本文方法產(chǎn)生的優(yōu)化波形性能要優(yōu)于常用波形,并且當(dāng)雜波的能量達(dá)到一定門(mén)限時(shí)(CNR=30),MI會(huì)變?yōu)榱?接收端幾乎得不到關(guān)于目標(biāo)的信息。

    圖6將本文算法與文獻(xiàn)[6]算法的耗時(shí)進(jìn)行了比較,可以看出,在相同波形長(zhǎng)度下,經(jīng)過(guò)一定的迭代收斂后,本文方法與文獻(xiàn)[6]的運(yùn)算量較為接近,當(dāng)波形長(zhǎng)度小于103時(shí),兩種方法的耗時(shí)均較小;當(dāng)波形長(zhǎng)度大于103時(shí),兩種方法的耗時(shí)均比較大,不利于實(shí)時(shí)信號(hào)處理。因此,本文算法適合序列碼長(zhǎng)相對(duì)較短的波形。

    圖6 不同算法的運(yùn)算量對(duì)比Fig.6 Run-time of different algorithm

    4.2 實(shí)驗(yàn)2:算法性能與PAR之間的關(guān)系

    圖7給出了不同PAR約束下,最優(yōu)波形和滿足Toeplitz結(jié)構(gòu)的優(yōu)化波形的互信息取值??梢钥闯?隨著μ的增大,問(wèn)題P2的可行解區(qū)域逐漸變大,因此所求波形的估計(jì)性能也有所增強(qiáng)。當(dāng)μ≥3時(shí),所設(shè)計(jì)波形的互信息非常接近能量約束下的最大互信息值。此外,由于進(jìn)行了波形矩陣的逼近,最后得到優(yōu)化波形向量的互信息相比于最優(yōu)波形有輕微的減少。

    圖7 互信息與PAR的關(guān)系Fig.7 Mutual information versus PAR

    圖8為不同PAR下滿足Toeplitz結(jié)構(gòu)的優(yōu)化波形實(shí)部和虛部的表示。μ=Ns代表能量約束下的波形,此時(shí)圖中對(duì)應(yīng)點(diǎn)的分布半徑較大,表明波形幅度起伏較大,不能保證所設(shè)計(jì)波形的恒模特性,不利于實(shí)際應(yīng)用。當(dāng)μ=1時(shí),產(chǎn)生的點(diǎn)位于單位圓上,是恒模波形,當(dāng)μ=2時(shí),圖中點(diǎn)的分布半徑比恒模波形略大,估計(jì)性能要更強(qiáng)。

    圖8 不同PAR波形的實(shí)部和虛部Fig.8 Real and imaginary parts of different PAR

    4.3 實(shí)驗(yàn)3:算法性能與距離量Δλ之間的關(guān)系

    圖9 距離區(qū)間對(duì)SLP算法性能的影響Fig.9 Influence of distance intervals on performance of SLP

    圖9(a)給出了不同距離Δλ對(duì)SLP算法的影響,I代表原始函數(shù)的互信息取值,I2代表采用SLP算法后的一階近似函數(shù)的互信息取值,而最優(yōu)值通過(guò)注水法得到。從圖9中可以看出,當(dāng)Δλ取值區(qū)間較小時(shí)(η≤10-1),由于收斂速度較慢,每次迭代后函數(shù)值幾乎無(wú)變化,因此最終經(jīng)過(guò)κ次收斂后得到的I和I2的取值也并沒(méi)有收斂到最優(yōu)值;當(dāng)η=4時(shí),I可以收斂到最優(yōu)值并且I與I2的取值非常接近;當(dāng)Δλ取值區(qū)間較大時(shí)(η≥10),收斂速度較快,但容易收斂到局部最優(yōu)值,并且此時(shí)原始函數(shù)和一階近似函數(shù)相差較大,因此圖9(a)中I和I2的值也相差較大,同時(shí),可以注意到,此時(shí)的I和I2分別收斂到一個(gè)固定值,這是因?yàn)楫?dāng)η大于一定門(mén)限時(shí),Δλ的最優(yōu)值也將收斂到一個(gè)固定值Δλopt,圖9(b)給出了Δλopt隨的變化曲線圖,縱坐標(biāo)代表Δλopt的F范數(shù)取值,可以看出,當(dāng)η≥10時(shí),Δλopt收斂到一個(gè)固定值。

    5 結(jié) 論

    為了提高雷達(dá)發(fā)射機(jī)的效能,避免放大器的非線性使發(fā)射波形失真,提出了一種信號(hào)相關(guān)雜波條件下雷達(dá)低PAR估計(jì)波形設(shè)計(jì)方法,在能量與PAR的雙重約束下,追求互信息的最大化。從時(shí)域角度出發(fā)構(gòu)建了問(wèn)題模型,引入序列線性規(guī)劃的思想將非凸優(yōu)化問(wèn)題轉(zhuǎn)化為一個(gè)凸問(wèn)題進(jìn)行研究,并基于F-范數(shù)對(duì)最優(yōu)矩陣波形進(jìn)行逼近,從而求得優(yōu)化波形向量。本文方法在給定的低PAR區(qū)間內(nèi)可有效求解得到估計(jì)波形。在算法運(yùn)算量相同的情況下,與現(xiàn)有方法相比,所提方法產(chǎn)生的波形具有更好的估計(jì)性能。

    猜你喜歡
    優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    PEMFC流道的多目標(biāo)優(yōu)化
    能源工程(2022年1期)2022-03-29 01:06:28
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
    事業(yè)單位中固定資產(chǎn)會(huì)計(jì)處理的優(yōu)化
    4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
    幾種常見(jiàn)的負(fù)載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    九草在线视频观看| 亚洲欧美精品自产自拍| 亚洲伊人久久精品综合| 97超碰精品成人国产| 欧美xxxx性猛交bbbb| 国产精品人妻久久久久久| av在线观看视频网站免费| 国产毛片在线视频| 国产黄色视频一区二区在线观看| 亚洲综合精品二区| 美女福利国产在线| 交换朋友夫妻互换小说| 久久女婷五月综合色啪小说| 精品国产乱码久久久久久小说| 国产xxxxx性猛交| 男男h啪啪无遮挡| 国产国拍精品亚洲av在线观看| 高清视频免费观看一区二区| 一二三四在线观看免费中文在 | 五月伊人婷婷丁香| 少妇被粗大的猛进出69影院 | 高清毛片免费看| 99热全是精品| 国产日韩一区二区三区精品不卡| 国产成人午夜福利电影在线观看| 精品熟女少妇av免费看| 久久人人97超碰香蕉20202| 激情视频va一区二区三区| 国产精品久久久久久av不卡| 成人国产麻豆网| 美女大奶头黄色视频| 欧美精品高潮呻吟av久久| 久久 成人 亚洲| 大片电影免费在线观看免费| 久久这里只有精品19| 国产一区二区三区av在线| 免费黄频网站在线观看国产| 日日爽夜夜爽网站| 久久久精品94久久精品| 亚洲av欧美aⅴ国产| 18+在线观看网站| 国产精品一二三区在线看| 9191精品国产免费久久| 久久这里只有精品19| 国产永久视频网站| 亚洲精品一二三| 黄色怎么调成土黄色| 十八禁网站网址无遮挡| 久久久精品免费免费高清| 97超碰精品成人国产| 欧美+日韩+精品| 十八禁网站网址无遮挡| 免费看光身美女| 久久久久久久亚洲中文字幕| a级毛片黄视频| 老司机亚洲免费影院| 久久久精品94久久精品| 日本wwww免费看| 波多野结衣一区麻豆| 欧美日韩成人在线一区二区| 2022亚洲国产成人精品| 亚洲av国产av综合av卡| 欧美日韩综合久久久久久| 人人妻人人爽人人添夜夜欢视频| 天天躁夜夜躁狠狠躁躁| 成年人午夜在线观看视频| 黄色怎么调成土黄色| 自拍欧美九色日韩亚洲蝌蚪91| 一级,二级,三级黄色视频| 十八禁网站网址无遮挡| 国产精品麻豆人妻色哟哟久久| av女优亚洲男人天堂| 毛片一级片免费看久久久久| 狂野欧美激情性xxxx在线观看| 国产精品人妻久久久影院| 成人国产麻豆网| 91国产中文字幕| 国产成人精品无人区| 夫妻性生交免费视频一级片| 精品人妻熟女毛片av久久网站| av视频免费观看在线观看| 高清毛片免费看| 一级黄片播放器| 亚洲国产色片| 99九九在线精品视频| 国产av一区二区精品久久| 久久久久国产网址| 国产亚洲精品第一综合不卡 | 久久狼人影院| 久久人人爽人人片av| av在线播放精品| 最后的刺客免费高清国语| 国产精品久久久久久精品电影小说| www.色视频.com| av在线观看视频网站免费| 18在线观看网站| 一级毛片黄色毛片免费观看视频| 在线精品无人区一区二区三| av视频免费观看在线观看| 亚洲精品一区蜜桃| 一本久久精品| 亚洲婷婷狠狠爱综合网| 亚洲av电影在线观看一区二区三区| 亚洲综合色惰| 国产精品国产三级专区第一集| 午夜福利网站1000一区二区三区| 九九在线视频观看精品| 国产高清三级在线| 黄网站色视频无遮挡免费观看| 国产av国产精品国产| av国产精品久久久久影院| 亚洲av男天堂| 婷婷色麻豆天堂久久| 国产高清不卡午夜福利| 一区二区日韩欧美中文字幕 | 狠狠精品人妻久久久久久综合| 黄网站色视频无遮挡免费观看| 久久精品久久精品一区二区三区| 一区二区三区四区激情视频| 最近中文字幕高清免费大全6| 日韩av在线免费看完整版不卡| 免费av中文字幕在线| 女性生殖器流出的白浆| videosex国产| 捣出白浆h1v1| 国产成人精品福利久久| 国产一区二区激情短视频 | 成人影院久久| 午夜久久久在线观看| 男女国产视频网站| 中国国产av一级| 国产国拍精品亚洲av在线观看| 亚洲成国产人片在线观看| 亚洲av.av天堂| 少妇人妻久久综合中文| 高清av免费在线| 免费久久久久久久精品成人欧美视频 | 国产亚洲一区二区精品| 日韩 亚洲 欧美在线| 一级毛片我不卡| 中文字幕av电影在线播放| 精品第一国产精品| 国产高清国产精品国产三级| 精品久久蜜臀av无| 亚洲国产欧美在线一区| 草草在线视频免费看| 久久精品国产自在天天线| 日韩精品有码人妻一区| 国产一区二区在线观看av| 成人国语在线视频| 精品午夜福利在线看| 久久人人爽人人爽人人片va| 中国三级夫妇交换| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品色激情综合| 最近最新中文字幕免费大全7| 十分钟在线观看高清视频www| 桃花免费在线播放| 成人国语在线视频| 七月丁香在线播放| www.色视频.com| 日本色播在线视频| 大香蕉久久成人网| 欧美人与性动交α欧美精品济南到 | 美国免费a级毛片| 一级片'在线观看视频| 大片电影免费在线观看免费| 成人影院久久| 午夜91福利影院| 韩国精品一区二区三区 | 国产极品粉嫩免费观看在线| 亚洲一级一片aⅴ在线观看| 最近最新中文字幕免费大全7| 中文字幕人妻熟女乱码| 免费人妻精品一区二区三区视频| 一级片免费观看大全| freevideosex欧美| 香蕉丝袜av| 免费av不卡在线播放| 18禁国产床啪视频网站| 日本91视频免费播放| 亚洲国产毛片av蜜桃av| 欧美 日韩 精品 国产| 99热全是精品| 国产精品国产三级国产专区5o| 青春草亚洲视频在线观看| 欧美精品一区二区免费开放| 欧美精品一区二区免费开放| 日本与韩国留学比较| 最新中文字幕久久久久| 亚洲国产精品专区欧美| 又粗又硬又长又爽又黄的视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美一区二区三区黑人 | 久久人妻熟女aⅴ| 男女下面插进去视频免费观看 | 久久 成人 亚洲| 久久综合国产亚洲精品| 国产精品女同一区二区软件| 久久久久久久精品精品| 久久人人爽人人爽人人片va| 青春草亚洲视频在线观看| videossex国产| 久久人人爽人人片av| 女的被弄到高潮叫床怎么办| 亚洲欧美色中文字幕在线| 午夜福利,免费看| 久久久国产精品麻豆| 免费看不卡的av| 在线天堂最新版资源| 91午夜精品亚洲一区二区三区| 在线天堂中文资源库| 捣出白浆h1v1| 免费av不卡在线播放| 美女主播在线视频| 伦理电影免费视频| 免费日韩欧美在线观看| 亚洲成人一二三区av| 日本91视频免费播放| 中国美白少妇内射xxxbb| 全区人妻精品视频| 免费播放大片免费观看视频在线观看| 大片电影免费在线观看免费| 日本爱情动作片www.在线观看| 国产69精品久久久久777片| 久久精品久久久久久久性| 精品一区二区三区视频在线| av电影中文网址| 亚洲va日本ⅴa欧美va伊人久久| 久久久精品区二区三区| 精品视频人人做人人爽| 免费观看人在逋| 免费看十八禁软件| 制服诱惑二区| 免费不卡黄色视频| 久久久久国产一级毛片高清牌| 夫妻午夜视频| 淫妇啪啪啪对白视频| 脱女人内裤的视频| 欧美日韩国产mv在线观看视频| 一二三四在线观看免费中文在| 变态另类成人亚洲欧美熟女 | 美女午夜性视频免费| www.999成人在线观看| 亚洲av电影在线进入| 一级片免费观看大全| 亚洲欧美激情综合另类| 亚洲一区二区三区欧美精品| 在线观看免费午夜福利视频| 精品无人区乱码1区二区| 亚洲欧洲精品一区二区精品久久久| 久久精品亚洲精品国产色婷小说| 久久久久久久国产电影| 亚洲成人国产一区在线观看| 亚洲av欧美aⅴ国产| 大码成人一级视频| 亚洲欧美日韩高清在线视频| 成年人午夜在线观看视频| 少妇粗大呻吟视频| 国产一区二区激情短视频| 久久亚洲真实| 制服人妻中文乱码| ponron亚洲| 俄罗斯特黄特色一大片| 后天国语完整版免费观看| 人妻久久中文字幕网| 久久热在线av| 日日爽夜夜爽网站| 天堂俺去俺来也www色官网| 一级a爱片免费观看的视频| 日本欧美视频一区| 国产精品一区二区在线观看99| a在线观看视频网站| av网站在线播放免费| 精品久久久精品久久久| 久久久精品区二区三区| 巨乳人妻的诱惑在线观看| 大香蕉久久网| 久99久视频精品免费| 最新美女视频免费是黄的| 妹子高潮喷水视频| 精品一区二区三卡| 满18在线观看网站| 免费观看a级毛片全部| 亚洲精品国产一区二区精华液| 欧美最黄视频在线播放免费 | 午夜精品在线福利| 亚洲精品国产色婷婷电影| 欧美 亚洲 国产 日韩一| 久久精品熟女亚洲av麻豆精品| 狠狠婷婷综合久久久久久88av| 国产男女内射视频| 精品一区二区三区视频在线观看免费 | 最新在线观看一区二区三区| 人妻丰满熟妇av一区二区三区 | 亚洲成av片中文字幕在线观看| 一区在线观看完整版| 亚洲成国产人片在线观看| 亚洲欧美日韩高清在线视频| 精品久久久精品久久久| av电影中文网址| 不卡av一区二区三区| 成年人午夜在线观看视频| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美在线一区二区| 久久九九热精品免费| 日日摸夜夜添夜夜添小说| www.自偷自拍.com| 国产成人精品在线电影| 美女 人体艺术 gogo| 窝窝影院91人妻| 建设人人有责人人尽责人人享有的| 一进一出好大好爽视频| 最近最新中文字幕大全免费视频| xxx96com| 曰老女人黄片| 国产淫语在线视频| 国产精品国产高清国产av | 精品乱码久久久久久99久播| 大型av网站在线播放| 国产成人av教育| 久久 成人 亚洲| 国产欧美日韩一区二区三区在线| 亚洲成av片中文字幕在线观看| 一区在线观看完整版| 亚洲色图综合在线观看| 国产精品乱码一区二三区的特点 | 性色av乱码一区二区三区2| 操出白浆在线播放| 黑人欧美特级aaaaaa片| 亚洲 国产 在线| 亚洲精品粉嫩美女一区| 久久香蕉激情| 成人18禁高潮啪啪吃奶动态图| 国产精品99久久99久久久不卡| 免费看a级黄色片| 丝袜美腿诱惑在线| 女警被强在线播放| 人人澡人人妻人| 看黄色毛片网站| 高清欧美精品videossex| 每晚都被弄得嗷嗷叫到高潮| 在线观看一区二区三区激情| 国产亚洲精品久久久久久毛片 | 母亲3免费完整高清在线观看| 亚洲av熟女| 久久精品国产亚洲av高清一级| 男女床上黄色一级片免费看| 后天国语完整版免费观看| 9191精品国产免费久久| 老司机靠b影院| 十八禁网站免费在线| 波多野结衣一区麻豆| 淫妇啪啪啪对白视频| 婷婷成人精品国产| www.999成人在线观看| 后天国语完整版免费观看| 精品卡一卡二卡四卡免费| 久久精品国产a三级三级三级| 国产精品二区激情视频| 色精品久久人妻99蜜桃| 欧美中文综合在线视频| 水蜜桃什么品种好| 久久精品aⅴ一区二区三区四区| 久久久久久久精品吃奶| 久久影院123| 婷婷成人精品国产| 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 久热爱精品视频在线9| 久久久水蜜桃国产精品网| 日韩精品免费视频一区二区三区| 亚洲第一青青草原| 欧美日韩一级在线毛片| 看黄色毛片网站| 亚洲精品久久成人aⅴ小说| 99精品久久久久人妻精品| 日韩免费高清中文字幕av| 丁香欧美五月| 久久亚洲真实| 伦理电影免费视频| 久久香蕉精品热| 成人亚洲精品一区在线观看| 18禁黄网站禁片午夜丰满| 男男h啪啪无遮挡| 国产不卡一卡二| 又大又爽又粗| 免费在线观看日本一区| 搡老岳熟女国产| 少妇被粗大的猛进出69影院| 日韩人妻精品一区2区三区| 一级毛片女人18水好多| 日韩制服丝袜自拍偷拍| 天天添夜夜摸| 亚洲综合色网址| 国产乱人伦免费视频| 男女床上黄色一级片免费看| 国产成人一区二区三区免费视频网站| 亚洲专区字幕在线| 黄色女人牲交| a在线观看视频网站| 国产成人精品久久二区二区免费| 大型黄色视频在线免费观看| 久久精品人人爽人人爽视色| 亚洲欧美激情综合另类| 亚洲国产精品合色在线| 在线永久观看黄色视频| 一区在线观看完整版| 亚洲国产欧美网| 一级a爱片免费观看的视频| 亚洲国产欧美一区二区综合| 国产日韩一区二区三区精品不卡| 久久婷婷成人综合色麻豆| 大香蕉久久网| 十分钟在线观看高清视频www| 欧美不卡视频在线免费观看 | 国产成人av教育| 黄色成人免费大全| 精品久久久久久久毛片微露脸| 黑人猛操日本美女一级片| 国产精品久久久久久精品古装| 午夜亚洲福利在线播放| 国产精品 欧美亚洲| 热re99久久精品国产66热6| 叶爱在线成人免费视频播放| 大片电影免费在线观看免费| 热re99久久精品国产66热6| 黄色视频,在线免费观看| svipshipincom国产片| 亚洲 国产 在线| 久久精品91无色码中文字幕| 国产人伦9x9x在线观看| 丝瓜视频免费看黄片| 精品福利永久在线观看| 免费观看精品视频网站| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩成人在线一区二区| 精品高清国产在线一区| 精品国内亚洲2022精品成人 | 国产成人欧美在线观看 | 中国美女看黄片| 久久 成人 亚洲| av在线播放免费不卡| 视频区图区小说| 欧美国产精品va在线观看不卡| 国产精品二区激情视频| 中出人妻视频一区二区| 精品人妻1区二区| x7x7x7水蜜桃| 麻豆成人av在线观看| 叶爱在线成人免费视频播放| 9色porny在线观看| 久久国产乱子伦精品免费另类| 亚洲精品自拍成人| 一区二区三区国产精品乱码| av在线播放免费不卡| 国产精品秋霞免费鲁丝片| 精品福利观看| 宅男免费午夜| 亚洲免费av在线视频| 亚洲成人免费电影在线观看| 久久中文字幕人妻熟女| 在线播放国产精品三级| 黑人巨大精品欧美一区二区蜜桃| 国产精品 国内视频| 高潮久久久久久久久久久不卡| 99re在线观看精品视频| 亚洲色图 男人天堂 中文字幕| 亚洲午夜精品一区,二区,三区| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线免费观看视频4| 欧美激情 高清一区二区三区| 国产成人欧美在线观看 | 天天躁夜夜躁狠狠躁躁| av网站免费在线观看视频| 国产激情欧美一区二区| 精品国产一区二区三区四区第35| 久久香蕉激情| 亚洲色图av天堂| 极品教师在线免费播放| 国产99久久九九免费精品| 国产xxxxx性猛交| 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲| 热99久久久久精品小说推荐| 国产极品粉嫩免费观看在线| 天天躁日日躁夜夜躁夜夜| 亚洲美女黄片视频| 欧美日韩成人在线一区二区| 香蕉久久夜色| 亚洲国产精品sss在线观看 | 搡老岳熟女国产| 一级作爱视频免费观看| 下体分泌物呈黄色| av有码第一页| 高清在线国产一区| 精品亚洲成a人片在线观看| 国产成人影院久久av| 香蕉久久夜色| 老汉色av国产亚洲站长工具| 国产激情久久老熟女| 18禁观看日本| 最新的欧美精品一区二区| 国产免费男女视频| 亚洲熟妇熟女久久| 校园春色视频在线观看| 丝瓜视频免费看黄片| 黄片大片在线免费观看| www日本在线高清视频| 女人久久www免费人成看片| 丰满人妻熟妇乱又伦精品不卡| 乱人伦中国视频| 大香蕉久久成人网| 欧美不卡视频在线免费观看 | 久久香蕉精品热| 高清黄色对白视频在线免费看| 国产成人精品久久二区二区免费| 午夜福利,免费看| 国产成人啪精品午夜网站| 亚洲第一av免费看| 亚洲成a人片在线一区二区| 身体一侧抽搐| 69av精品久久久久久| 久久精品国产综合久久久| 国产一区二区三区综合在线观看| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 国产精品 国内视频| 一级毛片女人18水好多| 在线永久观看黄色视频| 成人三级做爰电影| 午夜免费成人在线视频| 成年人黄色毛片网站| 国产精品亚洲一级av第二区| 大型av网站在线播放| 色婷婷久久久亚洲欧美| 亚洲人成77777在线视频| 国精品久久久久久国模美| 久久狼人影院| 亚洲成a人片在线一区二区| 嫩草影视91久久| 亚洲色图 男人天堂 中文字幕| 丰满饥渴人妻一区二区三| 最新的欧美精品一区二区| 变态另类成人亚洲欧美熟女 | 777久久人妻少妇嫩草av网站| 大型av网站在线播放| 国产亚洲精品一区二区www | 国产片内射在线| 无遮挡黄片免费观看| 精品国产一区二区久久| 黄色 视频免费看| 国产熟女午夜一区二区三区| 免费在线观看影片大全网站| 青草久久国产| 女性被躁到高潮视频| 母亲3免费完整高清在线观看| 女性被躁到高潮视频| 中文字幕人妻丝袜制服| 国产在视频线精品| 叶爱在线成人免费视频播放| 色播在线永久视频| 国产亚洲欧美在线一区二区| 热re99久久国产66热| 美女高潮到喷水免费观看| 怎么达到女性高潮| 日本精品一区二区三区蜜桃| 久久热在线av| 欧美在线黄色| 一级片'在线观看视频| 免费少妇av软件| 国产亚洲欧美精品永久| 国产精品久久久av美女十八| 操美女的视频在线观看| 国产精品影院久久| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 99热国产这里只有精品6| 日韩免费av在线播放| 国产欧美日韩一区二区精品| 国产99白浆流出| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 91国产中文字幕| 精品国产亚洲在线| 国产精品成人在线| 久久久精品区二区三区| 久久九九热精品免费| 亚洲熟女毛片儿| 欧美亚洲 丝袜 人妻 在线| 欧美日韩精品网址| 高清黄色对白视频在线免费看| 热99国产精品久久久久久7| 性色av乱码一区二区三区2| 欧美 亚洲 国产 日韩一| 天堂中文最新版在线下载| 超碰成人久久| 黄频高清免费视频| 色婷婷久久久亚洲欧美| cao死你这个sao货| 国产视频一区二区在线看| 国产精品偷伦视频观看了| 国产精品综合久久久久久久免费 | 无人区码免费观看不卡| 午夜精品国产一区二区电影| 午夜老司机福利片| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美98| 搡老岳熟女国产| 亚洲avbb在线观看| 精品人妻在线不人妻| 国产亚洲欧美精品永久| 成人18禁高潮啪啪吃奶动态图| 欧美老熟妇乱子伦牲交| 久久国产精品人妻蜜桃|