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

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

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

    郝天鐸, 崔 琛, 龔 陽, 孫從易

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

    0 引 言

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

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

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

    1 信號模型

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

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

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

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

    St+Sc+n=st+sc+n

    (1)

    (2)

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

    2 波形設(shè)計方法

    2.1 基于互信息方法的問題描述

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

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

    (3)

    式中,h(x|S)代表發(fā)射波形Toeplitz矩陣S已知時回波信號的熵;h(x|t,S)代表S和TIR均已知時回波信號的熵。當S已知時,TIR和回波信號服從聯(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)

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

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

    (9)

    (10)

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

    Rt=UΣtUH

    (11)

    Rc=UΣcUH

    (12)

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

    ?k,l∈(1,Nt)

    (13)

    同理,有

    Rn=VΣnVH

    (14)

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

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

    (15)

    (16)

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

    (17)

    至此,對時域上TIR和回波信號的互信息表達式進行了推導簡化,可以看出,式(17)與頻域上的互信息表達式較為相似。接下來將引入能量和PAR約束對時域發(fā)射波形進行求解。

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

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

    (18)

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

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

    (19)

    (20)

    (21)

    (22)

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

    2.3 引入SLP方法的問題求解

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

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

    (23)

    在該問題的約束條件下,通過對式(23)進行求解,可以得到此方程的最優(yōu)解x(1),接著將目標函數(shù)在x(1)處展開

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

    (24)

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

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

    (25)

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

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

    2.3.1 初始值求解

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

    (26)

    (27)

    (28)

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

    (29)

    (30)

    2.3.3 波形向量的逼近

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

    (31)

    d-sHg-gHs

    (32)

    (33)

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

    3 算法性能分析

    3.1 算法收斂性

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

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

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

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

    3.2 算法復雜度分析

    4 仿真實驗與分析

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

    4.1 實驗1:算法有效性驗證

    圖4 問題P和P2的收斂性比較Fig.4 Convergence comparison of P and P2

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

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

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

    圖6 不同算法的運算量對比Fig.6 Run-time of different algorithm

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

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

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

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

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

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

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

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

    5 結(jié) 論

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

    猜你喜歡
    優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    PEMFC流道的多目標優(yōu)化
    能源工程(2022年1期)2022-03-29 01:06:28
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
    圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
    事業(yè)單位中固定資產(chǎn)會計處理的優(yōu)化
    消費導刊(2018年8期)2018-05-25 13:20:08
    4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
    幾種常見的負載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    电影成人av| 久久午夜综合久久蜜桃| 亚洲人成电影观看| 久久精品91无色码中文字幕| 大型av网站在线播放| 国产日韩欧美在线精品| 人人澡人人妻人| 老司机午夜十八禁免费视频| 欧美成狂野欧美在线观看| 亚洲三区欧美一区| 蜜桃国产av成人99| av线在线观看网站| 亚洲美女黄片视频| 97在线人人人人妻| 午夜福利欧美成人| 国产伦理片在线播放av一区| 国产三级黄色录像| 成人免费观看视频高清| 香蕉丝袜av| 大陆偷拍与自拍| 国产亚洲欧美在线一区二区| 狠狠婷婷综合久久久久久88av| 久久99一区二区三区| 一级片免费观看大全| 欧美精品一区二区免费开放| 日本wwww免费看| av一本久久久久| 亚洲 欧美一区二区三区| 两个人看的免费小视频| 国产99久久九九免费精品| 天堂8中文在线网| 日日摸夜夜添夜夜添小说| 欧美性长视频在线观看| 日韩中文字幕欧美一区二区| 丝瓜视频免费看黄片| 丁香欧美五月| 欧美黑人精品巨大| 亚洲精品粉嫩美女一区| 欧美日韩亚洲综合一区二区三区_| 在线看a的网站| 黄频高清免费视频| 91av网站免费观看| 久久久久久久久久久久大奶| 欧美日韩成人在线一区二区| 午夜精品久久久久久毛片777| 午夜成年电影在线免费观看| 在线观看舔阴道视频| 国产日韩欧美视频二区| 国产精品久久久久久精品古装| 国产精品影院久久| 999久久久国产精品视频| 每晚都被弄得嗷嗷叫到高潮| 12—13女人毛片做爰片一| 一级黄色大片毛片| 国产高清视频在线播放一区| 久久人人97超碰香蕉20202| 无遮挡黄片免费观看| 丝袜喷水一区| 一区福利在线观看| 日韩欧美一区二区三区在线观看 | 一边摸一边做爽爽视频免费| 亚洲精品粉嫩美女一区| 在线av久久热| 久久精品亚洲熟妇少妇任你| 两个人免费观看高清视频| 亚洲国产欧美一区二区综合| 50天的宝宝边吃奶边哭怎么回事| 少妇粗大呻吟视频| 黄色视频不卡| 三级毛片av免费| 亚洲av成人一区二区三| 丁香欧美五月| 国产成人影院久久av| av天堂久久9| 麻豆乱淫一区二区| 久久久水蜜桃国产精品网| 精品国产乱码久久久久久小说| 日日夜夜操网爽| 激情视频va一区二区三区| 午夜激情av网站| 91国产中文字幕| 交换朋友夫妻互换小说| 国产精品 国内视频| 久久亚洲真实| 一夜夜www| 另类精品久久| 一边摸一边做爽爽视频免费| 亚洲成人免费电影在线观看| 一区二区三区国产精品乱码| 欧美乱码精品一区二区三区| 99riav亚洲国产免费| 男女床上黄色一级片免费看| 亚洲人成77777在线视频| 香蕉丝袜av| 国产免费福利视频在线观看| 亚洲自偷自拍图片 自拍| 咕卡用的链子| 黑人操中国人逼视频| 亚洲五月色婷婷综合| 丁香六月天网| 视频区欧美日本亚洲| 交换朋友夫妻互换小说| 久久影院123| 国产单亲对白刺激| 亚洲专区字幕在线| 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 热re99久久国产66热| 免费黄频网站在线观看国产| 纵有疾风起免费观看全集完整版| 国产在线精品亚洲第一网站| 国产精品影院久久| 一级,二级,三级黄色视频| 黄片播放在线免费| 色在线成人网| 高清毛片免费观看视频网站 | 一本一本久久a久久精品综合妖精| 中文欧美无线码| 啦啦啦视频在线资源免费观看| 韩国精品一区二区三区| 男女下面插进去视频免费观看| www.精华液| 亚洲熟女毛片儿| 久久久久精品人妻al黑| 成年女人毛片免费观看观看9 | 热re99久久精品国产66热6| 国产精品九九99| 两人在一起打扑克的视频| a级毛片黄视频| 国产精品免费一区二区三区在线 | 老司机福利观看| 九色亚洲精品在线播放| 亚洲精品国产精品久久久不卡| 色婷婷av一区二区三区视频| 亚洲中文字幕日韩| 女性生殖器流出的白浆| 色尼玛亚洲综合影院| 久久久久久免费高清国产稀缺| 欧美日韩中文字幕国产精品一区二区三区 | 免费久久久久久久精品成人欧美视频| 日韩中文字幕欧美一区二区| 免费观看人在逋| 久久久精品94久久精品| 欧美+亚洲+日韩+国产| 性色av乱码一区二区三区2| 午夜久久久在线观看| 另类精品久久| 国产成人精品无人区| 黄色丝袜av网址大全| 美女视频免费永久观看网站| 久久人妻福利社区极品人妻图片| 国产日韩一区二区三区精品不卡| 精品亚洲成a人片在线观看| cao死你这个sao货| 一级毛片女人18水好多| 欧美在线黄色| 69精品国产乱码久久久| 欧美乱码精品一区二区三区| 欧美黑人欧美精品刺激| 成年女人毛片免费观看观看9 | 老熟女久久久| 亚洲一区二区三区欧美精品| netflix在线观看网站| 99re在线观看精品视频| 国产极品粉嫩免费观看在线| 中文字幕精品免费在线观看视频| 亚洲久久久国产精品| 成人国产av品久久久| 欧美激情高清一区二区三区| 女警被强在线播放| 久久人妻av系列| 一本久久精品| 午夜两性在线视频| 少妇精品久久久久久久| 亚洲av第一区精品v没综合| 在线观看免费午夜福利视频| 在线av久久热| 我要看黄色一级片免费的| 天天躁狠狠躁夜夜躁狠狠躁| 91成人精品电影| 日本av手机在线免费观看| 淫妇啪啪啪对白视频| 老司机午夜福利在线观看视频 | 日韩中文字幕视频在线看片| 久久天堂一区二区三区四区| 国产精品.久久久| 久久九九热精品免费| 欧美成狂野欧美在线观看| 久久性视频一级片| 老鸭窝网址在线观看| 亚洲五月婷婷丁香| 午夜成年电影在线免费观看| 天天躁日日躁夜夜躁夜夜| 国产一区二区三区视频了| 一个人免费看片子| 99精品欧美一区二区三区四区| 在线看a的网站| 99热网站在线观看| 啦啦啦视频在线资源免费观看| 超碰97精品在线观看| 欧美黑人精品巨大| 丝袜人妻中文字幕| 亚洲成国产人片在线观看| 另类亚洲欧美激情| 午夜福利在线免费观看网站| 黑丝袜美女国产一区| 免费看十八禁软件| 久久这里只有精品19| 脱女人内裤的视频| 国产视频一区二区在线看| 成年动漫av网址| 久久中文字幕人妻熟女| 最黄视频免费看| 高潮久久久久久久久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 搡老岳熟女国产| 国产有黄有色有爽视频| 久久久国产成人免费| 精品少妇久久久久久888优播| 水蜜桃什么品种好| 成人免费观看视频高清| 两个人看的免费小视频| 久热这里只有精品99| 在线观看免费午夜福利视频| 国产精品亚洲一级av第二区| 日本a在线网址| 最近最新中文字幕大全电影3 | 亚洲欧美精品综合一区二区三区| 午夜福利免费观看在线| 久久午夜综合久久蜜桃| 欧美在线黄色| aaaaa片日本免费| 国产精品亚洲一级av第二区| 精品免费久久久久久久清纯 | 国产欧美日韩一区二区三区在线| 91成年电影在线观看| 俄罗斯特黄特色一大片| 99久久人妻综合| 久久精品熟女亚洲av麻豆精品| 欧美午夜高清在线| 亚洲av日韩在线播放| e午夜精品久久久久久久| 深夜精品福利| 国产成人精品久久二区二区91| 亚洲国产欧美在线一区| 欧美人与性动交α欧美软件| 又紧又爽又黄一区二区| 国产无遮挡羞羞视频在线观看| 国产精品99久久99久久久不卡| 一区福利在线观看| 久久中文字幕人妻熟女| 国产成人免费观看mmmm| av电影中文网址| 精品亚洲成a人片在线观看| 黄色视频,在线免费观看| 一区二区av电影网| 在线观看免费高清a一片| 伊人久久大香线蕉亚洲五| 一二三四社区在线视频社区8| 久9热在线精品视频| 国产成人av激情在线播放| 中文字幕制服av| 久久久久精品人妻al黑| 精品国产一区二区久久| 热re99久久国产66热| av天堂在线播放| 黄频高清免费视频| 成人18禁在线播放| 99久久99久久久精品蜜桃| 久久国产精品人妻蜜桃| 又紧又爽又黄一区二区| 蜜桃在线观看..| 波多野结衣一区麻豆| 脱女人内裤的视频| 精品人妻1区二区| 免费女性裸体啪啪无遮挡网站| 最新在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 色94色欧美一区二区| 欧美日韩福利视频一区二区| 99九九在线精品视频| 一个人免费看片子| 麻豆av在线久日| 亚洲精品国产一区二区精华液| 精品国产乱码久久久久久小说| 99久久99久久久精品蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av第一区精品v没综合| a在线观看视频网站| 久久久久久久精品吃奶| 777米奇影视久久| 搡老熟女国产l中国老女人| 日韩视频在线欧美| 国产成人精品久久二区二区91| 免费av中文字幕在线| 国产精品1区2区在线观看. | 黄色片一级片一级黄色片| 国产精品秋霞免费鲁丝片| 国产亚洲精品一区二区www | 国产xxxxx性猛交| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 一级黄色大片毛片| 久久久水蜜桃国产精品网| 99热网站在线观看| 色老头精品视频在线观看| 亚洲视频免费观看视频| 欧美日韩一级在线毛片| 午夜福利影视在线免费观看| 国产主播在线观看一区二区| 如日韩欧美国产精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 丝袜喷水一区| 男女边摸边吃奶| 久久99热这里只频精品6学生| 国产男靠女视频免费网站| 满18在线观看网站| 国产一区二区 视频在线| 国产男靠女视频免费网站| 国产成人精品无人区| 丰满迷人的少妇在线观看| 日韩制服丝袜自拍偷拍| 男女免费视频国产| 国产高清videossex| 丝袜在线中文字幕| 天天添夜夜摸| 亚洲欧洲精品一区二区精品久久久| 午夜视频精品福利| 亚洲av第一区精品v没综合| 久久国产精品人妻蜜桃| av天堂在线播放| 1024视频免费在线观看| 久久青草综合色| 人人妻,人人澡人人爽秒播| 久久免费观看电影| 国产免费视频播放在线视频| 国产精品麻豆人妻色哟哟久久| 亚洲精品自拍成人| 丁香六月欧美| 国产淫语在线视频| 国产极品粉嫩免费观看在线| 王馨瑶露胸无遮挡在线观看| 婷婷成人精品国产| 日韩免费高清中文字幕av| 女人高潮潮喷娇喘18禁视频| 国产一区二区三区综合在线观看| 午夜老司机福利片| 欧美性长视频在线观看| 午夜福利一区二区在线看| 香蕉国产在线看| 久久国产亚洲av麻豆专区| 成年动漫av网址| 在线观看免费日韩欧美大片| 高清在线国产一区| 999精品在线视频| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区在线不卡| 久久久精品区二区三区| 亚洲成人国产一区在线观看| 1024视频免费在线观看| 精品少妇久久久久久888优播| 国产男靠女视频免费网站| 桃红色精品国产亚洲av| 99国产精品99久久久久| 80岁老熟妇乱子伦牲交| 性少妇av在线| 激情视频va一区二区三区| 丝袜美腿诱惑在线| 久久国产精品大桥未久av| 1024视频免费在线观看| 丝袜美足系列| 在线永久观看黄色视频| 国产麻豆69| 国产精品一区二区精品视频观看| av有码第一页| 男女高潮啪啪啪动态图| 国产无遮挡羞羞视频在线观看| 男女之事视频高清在线观看| 久久人人爽av亚洲精品天堂| 丝袜喷水一区| 久久午夜亚洲精品久久| 不卡av一区二区三区| 一区福利在线观看| 国产精品亚洲av一区麻豆| 亚洲精品一卡2卡三卡4卡5卡| 91成年电影在线观看| 欧美大码av| 色在线成人网| av福利片在线| 欧美精品亚洲一区二区| 久久久水蜜桃国产精品网| 大片电影免费在线观看免费| 国产欧美日韩一区二区三区在线| 久久久久久久久久久久大奶| 一区二区三区国产精品乱码| 久热这里只有精品99| e午夜精品久久久久久久| 999精品在线视频| 亚洲精华国产精华精| 久久精品国产综合久久久| 人人妻人人爽人人添夜夜欢视频| 欧美成人免费av一区二区三区 | 日日爽夜夜爽网站| 日韩欧美一区二区三区在线观看 | 欧美人与性动交α欧美精品济南到| 下体分泌物呈黄色| 99精国产麻豆久久婷婷| 精品一区二区三区av网在线观看 | 桃红色精品国产亚洲av| 国产老妇伦熟女老妇高清| 国产成人精品在线电影| 精品一区二区三区av网在线观看 | 青草久久国产| 成年版毛片免费区| 精品国产乱码久久久久久男人| 99国产精品99久久久久| 亚洲久久久国产精品| 正在播放国产对白刺激| 亚洲第一av免费看| 亚洲国产av新网站| 精品一区二区三区视频在线观看免费 | 久久九九热精品免费| 水蜜桃什么品种好| 国产亚洲精品久久久久5区| 少妇猛男粗大的猛烈进出视频| 亚洲,欧美精品.| 成人三级做爰电影| 伊人久久大香线蕉亚洲五| 两性夫妻黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 男女边摸边吃奶| 一本大道久久a久久精品| 国产精品九九99| 精品一区二区三区av网在线观看 | 亚洲天堂av无毛| 少妇 在线观看| 男女下面插进去视频免费观看| 欧美日韩国产mv在线观看视频| 99在线人妻在线中文字幕 | 高清av免费在线| 无遮挡黄片免费观看| 午夜老司机福利片| 女人被躁到高潮嗷嗷叫费观| 99riav亚洲国产免费| 黄色视频在线播放观看不卡| 51午夜福利影视在线观看| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 夫妻午夜视频| 日韩视频在线欧美| 国产精品亚洲av一区麻豆| 99国产极品粉嫩在线观看| 成人黄色视频免费在线看| 一区二区三区激情视频| 日韩大码丰满熟妇| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看| 极品教师在线免费播放| 国产精品一区二区免费欧美| 一区二区日韩欧美中文字幕| av又黄又爽大尺度在线免费看| av天堂在线播放| 一进一出抽搐动态| 老汉色∧v一级毛片| 18在线观看网站| 777久久人妻少妇嫩草av网站| 免费女性裸体啪啪无遮挡网站| 亚洲av片天天在线观看| 成年人免费黄色播放视频| 久久人妻av系列| 久久天堂一区二区三区四区| 国产亚洲午夜精品一区二区久久| 国产日韩欧美在线精品| 欧美黄色淫秽网站| 精品人妻1区二区| 一进一出抽搐动态| 少妇粗大呻吟视频| 后天国语完整版免费观看| 国产视频一区二区在线看| 中亚洲国语对白在线视频| 亚洲色图av天堂| 精品免费久久久久久久清纯 | 少妇 在线观看| 母亲3免费完整高清在线观看| 热99国产精品久久久久久7| 久久中文字幕一级| 黑人操中国人逼视频| 久久精品人人爽人人爽视色| 黄色毛片三级朝国网站| 欧美激情极品国产一区二区三区| 日本av手机在线免费观看| 99精品欧美一区二区三区四区| 久久国产精品人妻蜜桃| 国产精品久久久久成人av| 久久国产精品大桥未久av| 大片免费播放器 马上看| videosex国产| 热re99久久国产66热| 午夜福利免费观看在线| 成人精品一区二区免费| 女性生殖器流出的白浆| 国产1区2区3区精品| 国产亚洲欧美在线一区二区| 精品少妇久久久久久888优播| 天堂动漫精品| 91大片在线观看| a级片在线免费高清观看视频| av一本久久久久| 在线永久观看黄色视频| 久久久久精品国产欧美久久久| 一区二区日韩欧美中文字幕| 欧美成人免费av一区二区三区 | 一本久久精品| 中文字幕人妻丝袜制服| 国产亚洲精品一区二区www | 国产成人精品在线电影| 中文字幕高清在线视频| 日韩视频一区二区在线观看| 天天躁日日躁夜夜躁夜夜| 中文字幕另类日韩欧美亚洲嫩草| 我要看黄色一级片免费的| 高清黄色对白视频在线免费看| 欧美 亚洲 国产 日韩一| videosex国产| 中文字幕高清在线视频| 国产在线免费精品| 久久国产亚洲av麻豆专区| 成人亚洲精品一区在线观看| 王馨瑶露胸无遮挡在线观看| 99精国产麻豆久久婷婷| 亚洲专区字幕在线| 亚洲一码二码三码区别大吗| 五月天丁香电影| 国产成人精品久久二区二区91| 欧美精品人与动牲交sv欧美| 一级毛片精品| 国产精品一区二区在线观看99| 亚洲成人手机| 欧美激情极品国产一区二区三区| 精品国产乱子伦一区二区三区| 国产免费av片在线观看野外av| av天堂久久9| 免费在线观看完整版高清| 美女扒开内裤让男人捅视频| 精品视频人人做人人爽| 久久久久精品国产欧美久久久| 亚洲男人天堂网一区| 青青草视频在线视频观看| 黄片小视频在线播放| 欧美精品av麻豆av| 宅男免费午夜| 精品国产乱子伦一区二区三区| 99在线人妻在线中文字幕 | 亚洲成国产人片在线观看| 日本欧美视频一区| 午夜久久久在线观看| 他把我摸到了高潮在线观看 | 国产精品二区激情视频| 99国产精品99久久久久| 丝袜美腿诱惑在线| 国产精品国产高清国产av | 老鸭窝网址在线观看| 亚洲熟女毛片儿| 欧美乱妇无乱码| 狠狠精品人妻久久久久久综合| 国产熟女午夜一区二区三区| 精品人妻1区二区| 久久久久视频综合| 午夜福利一区二区在线看| 99精品在免费线老司机午夜| 免费高清在线观看日韩| 亚洲欧美一区二区三区久久| 亚洲九九香蕉| 老司机在亚洲福利影院| 啪啪无遮挡十八禁网站| 精品高清国产在线一区| 黄色视频,在线免费观看| netflix在线观看网站| 成在线人永久免费视频| 国产亚洲精品一区二区www | 满18在线观看网站| 精品少妇一区二区三区视频日本电影| 制服诱惑二区| 国产精品影院久久| av网站在线播放免费| 蜜桃在线观看..| 国产av国产精品国产| 国产精品av久久久久免费| 黄色视频在线播放观看不卡| 大香蕉久久成人网| 纯流量卡能插随身wifi吗| 亚洲精品乱久久久久久| 在线观看免费视频网站a站| 俄罗斯特黄特色一大片| 中文字幕av电影在线播放| 国产男女超爽视频在线观看| 色播在线永久视频| 天天添夜夜摸| 亚洲午夜精品一区,二区,三区| a在线观看视频网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产av新网站| 天天影视国产精品| 亚洲精品国产色婷婷电影| 欧美激情极品国产一区二区三区| 99国产精品一区二区三区| a在线观看视频网站| 天天躁夜夜躁狠狠躁躁| 久久久久国内视频| 美女扒开内裤让男人捅视频| 成年动漫av网址| 国产高清国产精品国产三级| 人人妻人人爽人人添夜夜欢视频| 最近最新中文字幕大全电影3 | 啪啪无遮挡十八禁网站|