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

    基于脈沖位置幅度調(diào)制的距離速度同時測量

    2021-07-30 10:34:16黃少偉黃婉琳雷閏龍毛雪松
    計(jì)算機(jī)應(yīng)用 2021年7期
    關(guān)鍵詞:碼元激光雷達(dá)信噪比

    黃少偉,黃婉琳,雷閏龍,毛雪松*

    (1.肇慶學(xué)院計(jì)算機(jī)科學(xué)與軟件學(xué)院/大數(shù)據(jù)學(xué)院,廣東肇慶 526061;2.武漢科技大學(xué)信息科學(xué)與工程學(xué)院/人工智能學(xué)院,武漢 430081)

    0 引言

    車輛行駛的智能化是汽車工業(yè)發(fā)展的重要方向,實(shí)現(xiàn)方式上包括自主式和網(wǎng)聯(lián)式兩個大類[1-2],其中自主式智能駕駛要求通過測量裝置獲取道路中障礙物的位置和速度,用于計(jì)算出一條安全的行駛路徑[3]。測量包括主動和被動兩種方式[4],主動測量方式由測量裝置發(fā)射電磁波,并通過接收反射電磁波來獲取障礙物的特征參數(shù);被動測量方式由測量裝置的接收器感知障礙物的自然特性,如反射自然光、輻射熱等來獲取障礙物的特征參數(shù)。被動測量方式通常能獲得很高的解像度,但受環(huán)境因素干擾較大;主動測量方式雖然不能獲得被動測量方式的解像度,但受外部因素干擾小,測量的準(zhǔn)確度更高。主動測量的典型裝置有毫米波雷達(dá)和激光雷達(dá),毫米波雷達(dá)采用發(fā)射頻率調(diào)制連續(xù)波的方式實(shí)現(xiàn)障礙物距離和速度的同時測量[5-6]。由于毫米波雷達(dá)的空間解像度較差,單位時間內(nèi)的測量次數(shù)不必很多,因此完成一次測量的發(fā)射信號長度可以較長[7-8]。與毫米波雷達(dá)相比,激光雷達(dá)的空間解像度可提高近兩個數(shù)量級,因而要求在更短的時間內(nèi)能夠完成一次測量[9]。市場上大多數(shù)激光雷達(dá)采用發(fā)射短脈沖的方式完成目標(biāo)距離的測量,但因脈沖長度過短,不能覆蓋移動目標(biāo)產(chǎn)生的多普勒信號一個周期,無法直接測量目標(biāo)的速度。增大脈沖長度可以實(shí)現(xiàn)激光雷達(dá)的速度測量,但這樣做一方面使距離分辨率劣化[10],另一方面向空間輻射大量無助于目標(biāo)檢測的電磁能量,因此需要研制一種在光電器件上可以實(shí)現(xiàn)的、長度介于頻率調(diào)制連續(xù)波和短脈沖之間的測量信號波形,使激光雷達(dá)能夠保持較高的空間解像度,并完成對道路場景中障礙物距離和速度的測量。

    測量信號波形分為連續(xù)波和脈沖序列兩大類,文獻(xiàn)[11-12]中研究了連續(xù)波方式實(shí)現(xiàn)激光雷達(dá)距離和速度同時測量;但這些研究沒有考慮智能駕駛要求完成一次測量的信號波形非常短,通常只有幾微秒,也沒有考慮連續(xù)波方式在駕駛環(huán)境中存在的干擾問題。相對來說,脈沖方式的干擾概率遠(yuǎn)低于連續(xù)波方式,因此本課題組提出了幾種基于脈沖序列調(diào)制的激光雷達(dá)距離和速度同時測量方法[13-15],并通過仿真實(shí)驗(yàn)驗(yàn)證了這些方法的正確性,從理論上討論了這些方法用于道路環(huán)境中目標(biāo)檢測的可行性與可靠性;但這些方法存在一個共性的問題,即計(jì)算外差信號多普勒頻率時,需要根據(jù)距離測量結(jié)果確定外差信號中目標(biāo)回波的頭位置。在實(shí)際應(yīng)用中,距離測量雖然使用數(shù)據(jù)累加方法增強(qiáng)目標(biāo)回波信號檢測能力,但對接收機(jī)輸出信號信噪比(Signal-to-Noise Ratio,SNR)的要求遠(yuǎn)遠(yuǎn)高于傅里葉解析。通常距離測量要求信噪比約為0 dB,而傅里葉解析在信噪比為-30 dB 的情況下仍能計(jì)算出淹沒在噪聲中信號的頻率。因此,設(shè)計(jì)一種方法使速度測量獨(dú)立于距離測量具有非常重要的現(xiàn)實(shí)意義,當(dāng)某個方向檢測到存在運(yùn)動目標(biāo),但未能檢測出距離時,可自動增加該方向測量信號的發(fā)射功率,提高激光雷達(dá)對環(huán)境的自適應(yīng)能力。

    針對上述問題和應(yīng)用需求,本文在前期工作的基礎(chǔ)上,提出了一種在發(fā)射端同時調(diào)制脈沖位置和幅度的方法,即根據(jù)偽隨機(jī)碼“1”碼的位置確定發(fā)射脈沖的位置,根據(jù)相鄰“1”碼的間隔確定發(fā)射脈沖的幅度,實(shí)現(xiàn)對接收端輸出外差信號做頻率計(jì)算時,無須根據(jù)距離測量結(jié)果確定外差信號的頭位置,也無須根據(jù)速度測量結(jié)果來修正目標(biāo)的距離,從而使距離測量和速度測量相互獨(dú)立。

    1 系統(tǒng)框架與問題描述

    激光雷達(dá)測量裝置的系統(tǒng)框圖如圖1 所示,包括激光器、發(fā)射端、接收端、調(diào)制部分和信號處理部分等模塊,圖中的PD表示光電轉(zhuǎn)換器。調(diào)制碼發(fā)生器產(chǎn)生發(fā)射波形所需的位置幅度調(diào)制碼型,并作用于光外部調(diào)制器上,使發(fā)射端發(fā)射位置幅度調(diào)制的光信號。通常光調(diào)制器輸出光信號幅度只有幾毫瓦,為提高激光雷達(dá)對環(huán)境的測量能力,需要在調(diào)制器后增加光放大器提高發(fā)射信號的峰值功率。接收端使用光分離器將接收到的光回波信號分為兩路:一路用于距離測量,另一路用于速度測量。選擇使用雙檢測器結(jié)構(gòu)是因?yàn)楣庑盘柕念l率很高,1 550 nm 波長激光的頻率高達(dá)1.9×1014Hz,因此必須使用相干檢測的方式輸出包含于回波信號中的多普勒頻率。另一方面,相干檢測輸出的多普勒頻率對回波的脈沖序列產(chǎn)生幅度調(diào)制效應(yīng),脈沖序列中靠近多普勒正弦信號零點(diǎn)的脈沖受調(diào)制影響幅度明顯降低,對距離測量性能產(chǎn)生影響。

    圖1 激光雷達(dá)原理框圖Fig.1 Laser radar priciple block diagram

    對于采用偽隨機(jī)碼調(diào)制的激光雷達(dá),文獻(xiàn)[13]提出了非等間隔采樣信號頻譜分析方法:

    可以計(jì)算出PD1 輸出外差信號的多普勒頻率,測量出目標(biāo)的速度,其中f(ti)表示偽隨機(jī)碼“1”碼位置對應(yīng)的信號幅值,ΔTi表示相鄰“1”碼之間的時間間隔。與傅里葉變換的計(jì)算公式對比,可以看出為了計(jì)算多普勒信號的頻率,需要將PD1 輸出外差信號“1”碼對應(yīng)幅度倍乘相鄰“1”碼的時間間隔。如果PD1 輸出的外差信號淹沒在接收器的噪聲中,則無法確定“1”碼的位置,故不能用式(1)計(jì)算外差信號的頻率。為確定外差信號頭位置(頭位置確定,則所有“1”碼位置就可以確定),要利用PD2輸出信號并結(jié)合PD2和PD1之間固定的時延等信息計(jì)算出目標(biāo)距離。如果目標(biāo)距離不能正確測量,則目標(biāo)速度就不能正確測量。

    2 距離速度測量方法

    2.1 速度測量

    PD1輸出回波信號與本地信號的混頻信號可表示為:

    其中:φ包含因時延引起的相位變化;c(·)表示調(diào)制碼;T表示碼元寬度;ΔT表示從發(fā)射到接收之間的時延;表示向下取整。在調(diào)制碼為“1”的碼元區(qū)間內(nèi),s2(t)的采樣值與多普勒余弦信號A2cos(2πfDt+φ)一致;但在“0”碼元區(qū)間內(nèi),s2(t)的采樣值為0,偏離多普勒余弦信號值。為了能夠在接收端利用快速傅里葉變換(Fast Fourier Transform,F(xiàn)FT)方法計(jì)算出多普勒頻率,需根據(jù)式(1)對“1”碼元區(qū)間內(nèi)采樣信號幅值做修正。文獻(xiàn)[14]即為利用該方法,根據(jù)距離計(jì)算結(jié)果確定回波信號的頭,從而確定每個“1”碼所在的位置。為了使速度測量不依賴于距離測量的結(jié)果,可以將幅值修正過程調(diào)整到發(fā)射端,構(gòu)成新的調(diào)制碼,如圖2 所示。發(fā)射脈沖的位置由偽隨機(jī)碼“1”碼的位置決定,脈沖的幅度由相鄰碼符號之間的關(guān)系決定:如果“1”碼的后續(xù)碼元為“1”碼,則保持原碼幅度不變;如果“1”碼的后續(xù)碼元為“0”碼,則統(tǒng)計(jì)“0”碼的個數(shù);如果“1”碼的后續(xù)碼元“0”碼個數(shù)為Ni,則將該“1”碼的幅度乘以Ni+1。如果連“0”碼的個數(shù)過大,則前面的“1”碼幅值會過大,導(dǎo)致光器件無法實(shí)現(xiàn)過大功率光信號,因此需要在連“0”碼過大的位置插入“1”碼,本文選擇最大連“0”碼個數(shù)為3。

    圖2 調(diào)制碼的格式Fig.2 Format of modulation codes

    電光調(diào)制器輸出光強(qiáng)度與調(diào)制電壓之間的關(guān)系[16]為:

    式中:a為插入損耗;Iq為輸入光強(qiáng)度;I(t)為輸出光強(qiáng)度;V(t)為射頻端口的調(diào)制電壓;Vπ為直流驅(qū)動電壓。電光調(diào)制器后的光纖放大器在輸入光功率變化較小的范圍內(nèi)可視為增益不變,因此放大后的光信號碼元之間的幅度倍數(shù)關(guān)系可以與放大前碼元之間倍數(shù)關(guān)系保持不變。為使調(diào)制器輸出光信號幅度按照圖2 的位置幅度關(guān)系變化,可以根據(jù)式(4)設(shè)置電光調(diào)制器的射頻端口調(diào)制電壓。

    2.2 距離測量

    激光雷達(dá)的發(fā)射端輸出信號為:

    其中:A1表示發(fā)射信號的幅度,ft表示發(fā)射光信號的頻率,φ1為初始相位。

    由于光頻率太高,接收端光電轉(zhuǎn)換器PD2 只對光強(qiáng)度敏感,因此輸出信號為:

    其中:A2為PD2 輸出電信號的幅度,n(t)表示接收電路的噪聲。接收端距離測量部分信號處理器的任務(wù)是從光電轉(zhuǎn)換器輸出信號中測量出時延,并根據(jù)時延計(jì)算出目標(biāo)的距離。由于信號s3淹沒在噪聲中,不能通過儀器直接給出時延的值??梢岳梦墨I(xiàn)[14]計(jì)算相關(guān)函數(shù)確定峰值位置的方法,但是相關(guān)函數(shù)法需要按碼元個數(shù)移位相乘再相加,運(yùn)算量大。本文采用類似文獻(xiàn)[15]處理Golomb 序列調(diào)制的數(shù)據(jù)累加方法來確定偽隨機(jī)碼調(diào)制的回波信號時延,可以降低運(yùn)算量。

    如圖3 所示,第一行為PD2 輸出信號,作為原理演示,忽略了噪聲。信號處理部分檢測調(diào)制碼所有“0”碼到“1”碼的跳變時間點(diǎn),記錄相鄰跳變點(diǎn)之間的時間間隔。然后,根據(jù)這些時間間隔,將PD2輸出信號向左分別移位相應(yīng)的時間,構(gòu)成第二行、第三行……直到最后一個跳變點(diǎn)對應(yīng)的行。最后將所有行相加,得到和信號,和信號中出現(xiàn)由最小值到最大值的跳變,該峰值跳變位置對應(yīng)于回波信號相對于發(fā)射信號的時延。相比相關(guān)函數(shù)的方法,這種數(shù)據(jù)累加的方法需要移位的次數(shù)與從“0”碼到“1”碼的跳變次數(shù)相等,移位次數(shù)大幅下降,而且對信號運(yùn)算使用加法,而不是乘法,運(yùn)算量也大幅下降。統(tǒng)計(jì)表明,偽隨機(jī)碼的“0”碼到“1”碼的跳變次數(shù)約為碼長的1/4。

    圖3 脈沖位置幅度調(diào)制的數(shù)據(jù)累加方法Fig.3 Data accumulation method for pulse position and amplitude modulation

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

    3.1 測量原理的仿真驗(yàn)證

    假設(shè)激光雷達(dá)1 s 內(nèi)完成測量20 萬次,最大測量范圍為150 m,則可取測量信號序列時長為4 μs,等待時長為1 μs。取調(diào)制碼的碼元寬度為2 ns,則完成單次測量的調(diào)制碼序列包含2 000 個碼元。另一方面,假設(shè)激光雷達(dá)的發(fā)射波長為1 550 nm,運(yùn)動目標(biāo)的相對速度在0~360 km/h 范圍內(nèi),則相對運(yùn)動產(chǎn)生的多普勒頻率在0~129 MHz范圍內(nèi)。對每個碼元采樣一次,則129 MHz 信號一個周期內(nèi)可采樣點(diǎn)數(shù)大于3,測量信號長度可覆蓋多普勒信號500 多個周期,因此可以完成高頻信號頻率的計(jì)算。多普勒頻率低于一定值的情況下,可以認(rèn)為被測對象相對靜止,這類目標(biāo)只需測出相對距離即可。智能駕駛中期望測量的最低速度為人的步行速度,取1 m/s,仿真中取測量信號長度為5 μs,因此可以測量出的最低速度為0.2 m/s,可以滿足智能駕駛對速度參數(shù)的要求。

    圖4 給出PD2 輸出信噪比為0 dB 情況下,利用圖3 所給的數(shù)據(jù)累加算法獲得的結(jié)果。在信噪比為0 dB 的情況下,信號淹沒在噪聲中。經(jīng)過數(shù)據(jù)累加后,則出現(xiàn)一個明顯的最小值躍遷到最大值,利用此特征來判斷測量信號的飛行時間對應(yīng)的信號延遲。由于激光光束很窄,非球面透鏡的發(fā)散角可以做到0.1 mrad 左右,這樣即使在100 m 遠(yuǎn)處,光束的寬度只有1 cm,因此在一個光束內(nèi)出現(xiàn)兩個目標(biāo)的概率可認(rèn)為是零。在不考慮回波信號二次反射的情況下,PD2 的輸出信號中不會出現(xiàn)兩個目標(biāo)反射信號的疊加。即使考慮二次反射,通常二次反射信號非常弱,其影響可忽略不計(jì),因此在一次測量中只需找出一個目標(biāo)即可。通過尋找數(shù)據(jù)累加得到的信號中最大值和最小值,并檢驗(yàn)最大值和最小值時間間隔的方式即可判斷出信號幅值最大值位置是否對應(yīng)目標(biāo)距離。

    圖4 信噪比為0 dB的情況下數(shù)據(jù)累加后的信號Fig.4 Signal obtained by data accumulation when SNR equals to 0 dB

    為完成速度測量,可分析PD1輸出信號的多普勒頻率,在無噪聲的情況下,PD1 輸出的信號為幅度被調(diào)制的余弦波。圖5 的下圖給出幅度被調(diào)制的、頻率為1.29 MHz 的多普勒信號樣本,上圖給出一個周期的余弦信號細(xì)節(jié)表示。對這種形式的信號采樣,雖然得到的不是余弦信號的等間隔采樣信號,但由于信號幅值被預(yù)調(diào)制,因此可以直接使用快速傅里葉變換算法計(jì)算出余弦信號的頻率。

    圖5 PD1輸出的外差信號Fig.5 PD1 output heterodyne signal

    從每個脈沖序列的發(fā)射時刻開始計(jì)時,對PD1 輸出信號采樣長度為5 μs 的數(shù)據(jù)樣本,采樣頻率為500 MHz。當(dāng)信號淹沒于噪聲中時,不同于文獻(xiàn)[14]利用距離信息定位PD1 輸出信號中回波的方式,本文方法對5 μs 信號直接做快速傅里葉變換求出多普勒頻率。類似于文獻(xiàn)[14-15],取多普勒頻率分別為1.29 MHz、5 MHz、10 MHz、20 MHz、50 MHz、129 MHz這6 個頻率做快速傅里葉變換,在信噪比為-13 dB 的情況下得到如圖6所示的頻率譜。

    圖6 多普勒信號頻譜分析結(jié)果Fig.6 Spectrum analysis results of Doppler signal

    3.2 可靠性分析

    運(yùn)用快速傅里葉變換和圖3 所示的數(shù)據(jù)累加方法,對回波信號在信噪比為-40 dB 到0 dB 范圍內(nèi)各取樣點(diǎn)處做10 萬次速度測量和距離測量驗(yàn)證,統(tǒng)計(jì)給出發(fā)生速度測量和距離測量錯誤的概率。距離測量方面,將脈沖序列預(yù)先設(shè)定一個延時疊加到噪聲中,然后通過改變脈沖序列的幅度改變信噪比,最后在不同信噪比下比較數(shù)據(jù)累加方法獲得的峰值位置與預(yù)先設(shè)定的延時,如果二者不一致,則統(tǒng)計(jì)為測量錯誤。

    速度測量方面,將傅里葉變換的正頻率范圍內(nèi)頻譜峰值對應(yīng)的頻率和負(fù)頻率范圍內(nèi)頻譜峰值對應(yīng)的頻率取絕對值并求平均,作為多普勒信號的頻率??紤]到傅里葉變換本身存在頻率計(jì)算誤差,如果傅里葉變換計(jì)算得到的頻率和實(shí)際多普勒頻率相差0.2MHz,則統(tǒng)計(jì)為頻率測量錯誤。對距離和速度的測量錯誤統(tǒng)計(jì)結(jié)果如圖7 所示,當(dāng)信噪比分別為0 dB和-13 dB 時,距離測量和速度測量的誤檢測概率降為0,表明速度測量所需的信號功率比距離測量所需的信號功率低一個數(shù)量級,圖1中PD2下方的光分離器可選分光比為90∶10。

    圖7 距離測量和速度測量在不同信噪比下的誤檢測概率Fig.7 Error detection rates of range and speed measurement under different SNR

    距離測量誤差取決于脈沖的寬度和檢測裝置的硬件電路兩個方面,硬件電路的討論超出本文的范圍。脈沖的寬度可以決定距離測量的分辨率,與前期工作[15]一致,本文選取的脈沖寬度為2 ns,因此在硬件電路相同的情況下,所得的距離分辨率相同。速度測量方面,本文仿真中選取的脈沖序列長度為4 μs,比文獻(xiàn)[15]中的1 024 ns 時間窗口更長,因此測量精度更高。圖8 給出利用FFT 分析余弦等間隔采樣數(shù)據(jù)和PD1輸出數(shù)據(jù)頻率在1 MHz 到129 MHz 頻率范圍內(nèi)的誤差對比,仿真中選取的頻率增量步長為0.1 MHz,因此即使進(jìn)一步降低頻率增量步長,PD1 輸出數(shù)據(jù)利用FFT 獲得的頻率誤差不會超過0.2 MHz,對應(yīng)速度誤差0.16 m/s,可以滿足智能駕駛車輛環(huán)境感知的需求。

    圖8 FFT分析PD1輸出數(shù)據(jù)與余弦等間隔采樣數(shù)據(jù)的頻率誤差對比Fig.8 Comparison of frequency error between PD1 output data and cosine uniform sampling data by FFT

    對PD1 輸出信號做傅里葉變換,得到的峰值可能來自多普勒信號,也可能來自噪聲。圖7 給出當(dāng)信噪比大于-13 dB時可以確保頻譜的峰值均來自多普勒信號,因此對含噪信號在大于信噪比-13dB 的情況下計(jì)算頻率誤差,并與無噪聲情況對比,結(jié)果表明噪聲不構(gòu)成頻率計(jì)算誤差的影響因素。

    4 結(jié)語

    激光雷達(dá)在道路環(huán)境信息感知的應(yīng)用中,目前只能測量目標(biāo)的距離,而不能測量目標(biāo)的速度。前期工作提出的偽隨機(jī)碼調(diào)制方法從理論上可以實(shí)現(xiàn)激光雷達(dá)距離和速度同步測量,但是速度的測量必須建立在距離正確測量的基礎(chǔ)上,如果距離測量出現(xiàn)較大隨機(jī)誤差或測量數(shù)據(jù)丟失,則不能完成速度的測量。本文提出了一種利用偽隨機(jī)碼調(diào)制激光雷達(dá)發(fā)射信號脈沖序列位置和幅度的方法,并針對PD2 輸出波形特征設(shè)計(jì)數(shù)據(jù)累加,獲得目標(biāo)反射回波信號的飛行時間。對PD1輸出采樣獲得的非等間隔多普勒信號直接使用快速傅里葉變換方法,可計(jì)算出多普勒頻率。用計(jì)算機(jī)仿真的方法在多個信噪比下對距離和速度測量的穩(wěn)定性進(jìn)行分析,并與前期基于偽隨機(jī)碼調(diào)制的方法對比,驗(yàn)證了方法的可靠性。與前期工作對比,本文的主要工作在于使速度測量獨(dú)立于距離測量,在距離測量失敗的情況下,也可獨(dú)立完成目標(biāo)速度測量。由于速度測量對信噪比要求較低,通常速度測量會更穩(wěn)定,當(dāng)檢測到某個方向存在運(yùn)動目標(biāo),但無法測出目標(biāo)距離時,可通過算法控制提高該方向發(fā)射信號功率,實(shí)現(xiàn)目標(biāo)距離的測量,增強(qiáng)激光雷達(dá)的智能化。

    猜你喜歡
    碼元激光雷達(dá)信噪比
    手持激光雷達(dá)應(yīng)用解決方案
    北京測繪(2022年5期)2022-11-22 06:57:43
    法雷奧第二代SCALA?激光雷達(dá)
    汽車觀察(2021年8期)2021-09-01 10:12:41
    LFM-BPSK復(fù)合調(diào)制參數(shù)快速估計(jì)及碼元恢復(fù)
    基于深度學(xué)習(xí)的無人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    基于激光雷達(dá)通信的地面特征識別技術(shù)
    基于激光雷達(dá)的多旋翼無人機(jī)室內(nèi)定位與避障研究
    電子制作(2018年16期)2018-09-26 03:27:00
    低信噪比下LFMCW信號調(diào)頻參數(shù)估計(jì)
    電子測試(2018年11期)2018-06-26 05:56:02
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    基于極大似然準(zhǔn)則的短猝發(fā)信號盲解調(diào)
    保持信噪比的相位分解反褶積方法研究
    欧美极品一区二区三区四区| 欧美一区二区国产精品久久精品| 欧美日韩亚洲国产一区二区在线观看| 嫩草影院入口| 欧美日韩中文字幕国产精品一区二区三区| 丁香欧美五月| 欧美不卡视频在线免费观看| 99riav亚洲国产免费| 男人舔女人下体高潮全视频| 每晚都被弄得嗷嗷叫到高潮| 日本熟妇午夜| 夜夜躁狠狠躁天天躁| 69av精品久久久久久| 亚洲av成人av| 2021天堂中文幕一二区在线观| 亚洲av五月六月丁香网| 国产激情偷乱视频一区二区| 夜夜夜夜夜久久久久| 亚洲aⅴ乱码一区二区在线播放| 舔av片在线| 日韩欧美精品v在线| ponron亚洲| 亚洲精品在线美女| 99国产精品一区二区三区| 高清在线国产一区| 99国产极品粉嫩在线观看| 久久精品国产亚洲av涩爱 | 国产高清视频在线播放一区| 久久久久久久久中文| 久久伊人香网站| 成人性生交大片免费视频hd| а√天堂www在线а√下载| 亚洲熟妇中文字幕五十中出| av在线观看视频网站免费| 国产老妇女一区| 窝窝影院91人妻| 成人性生交大片免费视频hd| 久久久久久久午夜电影| 无遮挡黄片免费观看| 亚洲va日本ⅴa欧美va伊人久久| 国产在视频线在精品| 日韩有码中文字幕| 国产高清有码在线观看视频| 国产精品一及| 午夜福利在线在线| 亚洲欧美日韩东京热| 色尼玛亚洲综合影院| 又爽又黄无遮挡网站| 久久久久九九精品影院| 亚洲av免费在线观看| 中文字幕av成人在线电影| 国产伦在线观看视频一区| 亚洲精品乱码久久久v下载方式| 国产美女午夜福利| 十八禁网站免费在线| 亚洲乱码一区二区免费版| 变态另类成人亚洲欧美熟女| 天堂影院成人在线观看| 免费观看的影片在线观看| 男人的好看免费观看在线视频| 国产三级黄色录像| 午夜a级毛片| 深夜a级毛片| 国产精品一区二区三区四区免费观看 | 国产淫片久久久久久久久 | 女生性感内裤真人,穿戴方法视频| 亚洲人成网站在线播| 一区福利在线观看| a级一级毛片免费在线观看| 精品欧美国产一区二区三| 黄色女人牲交| 9191精品国产免费久久| 久久精品综合一区二区三区| 老司机深夜福利视频在线观看| 成年免费大片在线观看| 国产精华一区二区三区| 一区二区三区免费毛片| 岛国在线免费视频观看| 丰满乱子伦码专区| 久久亚洲精品不卡| 久久久成人免费电影| 国产精品影院久久| 成人特级黄色片久久久久久久| 亚洲色图av天堂| 久久久久国产精品人妻aⅴ院| 亚洲成av人片在线播放无| 一级黄片播放器| 国产黄色小视频在线观看| 亚洲精品色激情综合| 免费电影在线观看免费观看| 午夜影院日韩av| 精品福利观看| 中文字幕av在线有码专区| 国产精品永久免费网站| 欧美黄色淫秽网站| 熟女人妻精品中文字幕| 欧美在线一区亚洲| 久久亚洲真实| 一夜夜www| 亚州av有码| 美女大奶头视频| 国内精品美女久久久久久| 国产午夜精品论理片| 男女那种视频在线观看| 免费看日本二区| 久9热在线精品视频| 久久久久久久精品吃奶| 琪琪午夜伦伦电影理论片6080| 色综合亚洲欧美另类图片| 2021天堂中文幕一二区在线观| 熟女电影av网| 两性午夜刺激爽爽歪歪视频在线观看| 国产综合懂色| 如何舔出高潮| 少妇人妻精品综合一区二区 | 国产乱人视频| 99视频精品全部免费 在线| 午夜精品一区二区三区免费看| 精品99又大又爽又粗少妇毛片 | 国产一区二区三区在线臀色熟女| 国产乱人视频| 婷婷亚洲欧美| 成人高潮视频无遮挡免费网站| 欧美精品啪啪一区二区三区| 国产精品乱码一区二三区的特点| 成人欧美大片| 九色国产91popny在线| 天堂网av新在线| 国产av麻豆久久久久久久| 国产精品野战在线观看| 国产熟女xx| 欧美高清成人免费视频www| 可以在线观看毛片的网站| 在线看三级毛片| 久久久国产成人免费| 久久久国产成人免费| 精品人妻熟女av久视频| bbb黄色大片| 51国产日韩欧美| 精品久久久久久久末码| 欧美zozozo另类| 久久精品国产亚洲av香蕉五月| 国产黄片美女视频| 夜夜看夜夜爽夜夜摸| 亚洲 国产 在线| 国产一区二区在线观看日韩| 亚洲成人中文字幕在线播放| 色噜噜av男人的天堂激情| 在线观看一区二区三区| 午夜老司机福利剧场| 九色成人免费人妻av| 极品教师在线免费播放| 日日干狠狠操夜夜爽| 日本一二三区视频观看| 久久午夜亚洲精品久久| 欧美黑人巨大hd| www.999成人在线观看| 九九久久精品国产亚洲av麻豆| 丰满人妻一区二区三区视频av| 麻豆国产97在线/欧美| 国产精品影院久久| 国产乱人视频| 日本免费一区二区三区高清不卡| 不卡一级毛片| 永久网站在线| 亚洲中文字幕日韩| 在线天堂最新版资源| 免费无遮挡裸体视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日本视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 啪啪无遮挡十八禁网站| 日本免费a在线| 一个人观看的视频www高清免费观看| 国产精品精品国产色婷婷| 日本a在线网址| 97热精品久久久久久| 美女 人体艺术 gogo| 亚洲美女搞黄在线观看 | 国产高潮美女av| 亚洲成人精品中文字幕电影| 日韩欧美三级三区| a级毛片免费高清观看在线播放| 一个人免费在线观看的高清视频| 久久中文看片网| 高清日韩中文字幕在线| 一进一出好大好爽视频| 国产黄a三级三级三级人| 我要看日韩黄色一级片| 一个人观看的视频www高清免费观看| 一边摸一边抽搐一进一小说| www.www免费av| 91久久精品国产一区二区成人| 性插视频无遮挡在线免费观看| 男插女下体视频免费在线播放| aaaaa片日本免费| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕久久专区| 小蜜桃在线观看免费完整版高清| 老司机深夜福利视频在线观看| 久久国产精品影院| 精品一区二区三区视频在线观看免费| 中文字幕av成人在线电影| 精品久久久久久久久亚洲 | 精品不卡国产一区二区三区| 内射极品少妇av片p| 成人永久免费在线观看视频| 免费看光身美女| 全区人妻精品视频| 亚洲美女黄片视频| 亚洲精品成人久久久久久| 少妇人妻一区二区三区视频| 国产一区二区在线观看日韩| 亚洲美女视频黄频| 香蕉av资源在线| 国产蜜桃级精品一区二区三区| 国产精品一及| 国产成人福利小说| 一二三四社区在线视频社区8| 少妇的逼好多水| 一个人观看的视频www高清免费观看| 中亚洲国语对白在线视频| 91在线精品国自产拍蜜月| 狠狠狠狠99中文字幕| 夜夜爽天天搞| 给我免费播放毛片高清在线观看| 成人无遮挡网站| 色精品久久人妻99蜜桃| 欧美黑人巨大hd| 中亚洲国语对白在线视频| 嫩草影院入口| 日韩高清综合在线| 最后的刺客免费高清国语| 男人和女人高潮做爰伦理| 精品人妻熟女av久视频| 夜夜夜夜夜久久久久| 亚洲av电影不卡..在线观看| 欧洲精品卡2卡3卡4卡5卡区| 成人精品一区二区免费| 久久香蕉精品热| 久久国产乱子免费精品| 97超视频在线观看视频| 亚洲国产欧美人成| 日本免费a在线| 十八禁人妻一区二区| 99久久精品国产亚洲精品| 午夜久久久久精精品| 亚洲最大成人中文| 此物有八面人人有两片| 久久这里只有精品中国| 色吧在线观看| 欧美激情在线99| 99久久精品国产亚洲精品| 国产色婷婷99| 日本熟妇午夜| 久久香蕉精品热| 欧美三级亚洲精品| 亚洲一区二区三区色噜噜| 日本五十路高清| 亚洲国产欧洲综合997久久,| 在线观看午夜福利视频| 99久久精品国产亚洲精品| 欧美3d第一页| 国产精品一区二区免费欧美| 欧美乱色亚洲激情| 又紧又爽又黄一区二区| 身体一侧抽搐| 91麻豆精品激情在线观看国产| 国产高清有码在线观看视频| 午夜激情欧美在线| 日日夜夜操网爽| 三级男女做爰猛烈吃奶摸视频| 少妇人妻一区二区三区视频| 国产一区二区三区在线臀色熟女| 三级毛片av免费| 国产欧美日韩精品一区二区| h日本视频在线播放| 欧美在线黄色| 18禁在线播放成人免费| 成人鲁丝片一二三区免费| 欧美高清成人免费视频www| 搞女人的毛片| 麻豆久久精品国产亚洲av| 精品人妻偷拍中文字幕| av福利片在线观看| 中文字幕av成人在线电影| 国产成人福利小说| 久久久久国内视频| 欧美中文日本在线观看视频| 特大巨黑吊av在线直播| 亚洲欧美日韩高清在线视频| 18+在线观看网站| 在线观看免费视频日本深夜| 亚洲成人精品中文字幕电影| 亚洲熟妇中文字幕五十中出| 亚洲av熟女| 国产成人av教育| 成人亚洲精品av一区二区| 在线观看舔阴道视频| 精品人妻偷拍中文字幕| 国产国拍精品亚洲av在线观看| 嫩草影院入口| 欧美xxxx性猛交bbbb| 国产久久久一区二区三区| 可以在线观看毛片的网站| 白带黄色成豆腐渣| 在线观看av片永久免费下载| 黄色女人牲交| 听说在线观看完整版免费高清| 十八禁人妻一区二区| 首页视频小说图片口味搜索| 欧美一级a爱片免费观看看| 国产探花极品一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲精品亚洲一区二区| 亚洲熟妇熟女久久| 中国美女看黄片| aaaaa片日本免费| 夜夜躁狠狠躁天天躁| 成人欧美大片| 一进一出抽搐gif免费好疼| 亚洲片人在线观看| 性插视频无遮挡在线免费观看| 国产成人欧美在线观看| 久久久久久久久中文| 身体一侧抽搐| 深夜精品福利| 美女cb高潮喷水在线观看| 国产色爽女视频免费观看| 老司机深夜福利视频在线观看| 久久久久国产精品人妻aⅴ院| 亚洲精品亚洲一区二区| 长腿黑丝高跟| 国产美女午夜福利| 欧美高清成人免费视频www| av在线观看视频网站免费| 欧美性感艳星| 色av中文字幕| 9191精品国产免费久久| 日韩欧美 国产精品| 日韩欧美三级三区| 一边摸一边抽搐一进一小说| 国产aⅴ精品一区二区三区波| 国产熟女xx| 亚洲人成网站在线播放欧美日韩| 亚洲精品日韩av片在线观看| 成年女人永久免费观看视频| 精品一区二区三区视频在线| 欧美一级a爱片免费观看看| 日韩中文字幕欧美一区二区| 精品国产三级普通话版| 精品久久久久久久末码| 男女那种视频在线观看| 欧美日韩黄片免| 美女黄网站色视频| 一个人观看的视频www高清免费观看| 成熟少妇高潮喷水视频| 一个人看的www免费观看视频| 国产精品久久久久久精品电影| 国产精品久久视频播放| 亚洲人成网站在线播| 欧美激情国产日韩精品一区| av天堂在线播放| 免费大片18禁| 国产精品久久电影中文字幕| 小说图片视频综合网站| 五月玫瑰六月丁香| 久9热在线精品视频| 久久精品国产亚洲av天美| 国产精品国产高清国产av| 午夜福利欧美成人| 亚洲av成人av| 欧美中文日本在线观看视频| 亚洲欧美精品综合久久99| 日本三级黄在线观看| 欧美一级a爱片免费观看看| 国产亚洲精品综合一区在线观看| 欧美激情久久久久久爽电影| h日本视频在线播放| 亚洲专区国产一区二区| 亚洲 欧美 日韩 在线 免费| 三级国产精品欧美在线观看| 日韩免费av在线播放| 波多野结衣高清无吗| 日韩欧美精品v在线| 久久婷婷人人爽人人干人人爱| 黄片小视频在线播放| 丰满乱子伦码专区| 精品免费久久久久久久清纯| 中国美女看黄片| 亚洲国产色片| 性欧美人与动物交配| 高清在线国产一区| 欧美最黄视频在线播放免费| 亚洲国产日韩欧美精品在线观看| 欧美一区二区精品小视频在线| 精品久久国产蜜桃| 热99re8久久精品国产| 亚洲成av人片在线播放无| 欧美日本亚洲视频在线播放| 在线观看美女被高潮喷水网站 | 精品日产1卡2卡| 国产白丝娇喘喷水9色精品| 午夜a级毛片| 日韩欧美精品免费久久 | 国产一区二区在线观看日韩| 亚洲aⅴ乱码一区二区在线播放| 免费黄网站久久成人精品 | 亚洲欧美日韩卡通动漫| 亚洲人与动物交配视频| 一本综合久久免费| 深夜a级毛片| 亚洲欧美清纯卡通| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 一个人免费在线观看电影| 最后的刺客免费高清国语| 怎么达到女性高潮| 午夜免费成人在线视频| 欧美成人a在线观看| 亚洲在线自拍视频| 久久精品人妻少妇| 最后的刺客免费高清国语| 国产亚洲精品综合一区在线观看| 亚洲欧美日韩无卡精品| 欧美xxxx黑人xx丫x性爽| 99精品在免费线老司机午夜| 国产野战对白在线观看| 麻豆成人午夜福利视频| 国产黄片美女视频| 少妇的逼好多水| 国产成人av教育| 欧美日本视频| 免费在线观看影片大全网站| 亚洲午夜理论影院| 精品久久久久久成人av| av在线蜜桃| 一个人免费在线观看电影| 午夜亚洲福利在线播放| 欧美极品一区二区三区四区| 在线观看66精品国产| 久久亚洲精品不卡| 最近在线观看免费完整版| 欧美绝顶高潮抽搐喷水| 亚洲精品在线观看二区| 国模一区二区三区四区视频| 午夜精品一区二区三区免费看| 国产一级毛片七仙女欲春2| 国产一区二区三区在线臀色熟女| 国产69精品久久久久777片| 亚洲色图av天堂| 精品久久久久久久久亚洲 | 亚洲 欧美 日韩 在线 免费| 久9热在线精品视频| 亚洲激情在线av| 免费观看的影片在线观看| 中文在线观看免费www的网站| 三级国产精品欧美在线观看| 亚洲自偷自拍三级| 性色av乱码一区二区三区2| 亚洲成a人片在线一区二区| 老女人水多毛片| 成人特级黄色片久久久久久久| 久久久久久久午夜电影| 午夜激情欧美在线| 国产亚洲精品av在线| 国产免费av片在线观看野外av| 99国产精品一区二区蜜桃av| 午夜福利成人在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩中文字幕国产精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文字幕日韩| 在线免费观看的www视频| 一区二区三区免费毛片| 国产精品电影一区二区三区| 夜夜看夜夜爽夜夜摸| 久久99热这里只有精品18| 成人av在线播放网站| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 日本撒尿小便嘘嘘汇集6| 国内精品久久久久精免费| 少妇人妻一区二区三区视频| av在线天堂中文字幕| 老司机午夜十八禁免费视频| 国语自产精品视频在线第100页| 亚洲美女视频黄频| av女优亚洲男人天堂| 国产成人a区在线观看| 特级一级黄色大片| 宅男免费午夜| 久久午夜亚洲精品久久| 亚洲一区二区三区不卡视频| 自拍偷自拍亚洲精品老妇| 18美女黄网站色大片免费观看| 亚洲人成网站在线播| 成人av一区二区三区在线看| 成年女人永久免费观看视频| 久久草成人影院| 精品一区二区三区人妻视频| 久久精品综合一区二区三区| 高清在线国产一区| 国产免费一级a男人的天堂| 欧美激情在线99| 久久久精品大字幕| 男人舔女人下体高潮全视频| a在线观看视频网站| 欧美绝顶高潮抽搐喷水| av国产免费在线观看| 久9热在线精品视频| 18+在线观看网站| 欧美一区二区精品小视频在线| 亚洲最大成人av| 欧美丝袜亚洲另类 | 亚洲人成电影免费在线| 久久久久免费精品人妻一区二区| 国产91精品成人一区二区三区| 国产欧美日韩一区二区精品| 成年版毛片免费区| 国产一区二区亚洲精品在线观看| 天堂影院成人在线观看| 三级国产精品欧美在线观看| 国产欧美日韩精品亚洲av| 老司机午夜福利在线观看视频| 在线a可以看的网站| 亚洲性夜色夜夜综合| 欧美高清性xxxxhd video| 午夜福利18| 亚洲欧美清纯卡通| 欧美激情久久久久久爽电影| 亚洲激情在线av| 男人的好看免费观看在线视频| 欧美激情久久久久久爽电影| 国产精品免费一区二区三区在线| 老司机福利观看| 三级国产精品欧美在线观看| 欧美zozozo另类| 成人美女网站在线观看视频| 免费av毛片视频| 舔av片在线| 久久伊人香网站| 在线观看av片永久免费下载| 久久午夜亚洲精品久久| 91在线观看av| 九色成人免费人妻av| 首页视频小说图片口味搜索| 一个人看的www免费观看视频| 亚洲七黄色美女视频| 国产精品一区二区免费欧美| 精品国产三级普通话版| 男人的好看免费观看在线视频| 99热精品在线国产| 十八禁网站免费在线| 国产 一区 欧美 日韩| 欧美又色又爽又黄视频| 亚洲aⅴ乱码一区二区在线播放| 搡老熟女国产l中国老女人| 3wmmmm亚洲av在线观看| 精品无人区乱码1区二区| 韩国av一区二区三区四区| 亚洲美女视频黄频| 成人亚洲精品av一区二区| 99久久久亚洲精品蜜臀av| 亚洲最大成人手机在线| 在线播放无遮挡| 18禁黄网站禁片免费观看直播| 非洲黑人性xxxx精品又粗又长| h日本视频在线播放| 两个人视频免费观看高清| 亚洲av成人精品一区久久| 国产男靠女视频免费网站| 无人区码免费观看不卡| 激情在线观看视频在线高清| 俄罗斯特黄特色一大片| 国产在视频线在精品| 成人av一区二区三区在线看| 日本一二三区视频观看| 国产午夜精品论理片| 床上黄色一级片| 色尼玛亚洲综合影院| 日本 av在线| 精品一区二区三区视频在线| 午夜免费男女啪啪视频观看 | 午夜两性在线视频| 亚洲精华国产精华精| 国产精品电影一区二区三区| 亚洲激情在线av| 免费在线观看成人毛片| 国产欧美日韩精品一区二区| 五月玫瑰六月丁香| 蜜桃久久精品国产亚洲av| 3wmmmm亚洲av在线观看| 久99久视频精品免费| 久久久久久九九精品二区国产| 亚洲精品久久国产高清桃花| 国产精品亚洲一级av第二区| 麻豆国产97在线/欧美| 亚洲av成人精品一区久久| 国产亚洲精品久久久久久毛片| 中国美女看黄片| 亚洲,欧美,日韩| 少妇裸体淫交视频免费看高清| 国产精品一区二区免费欧美| 国产av在哪里看| 国产精华一区二区三区| 国产熟女xx| 在线观看美女被高潮喷水网站 | 制服丝袜大香蕉在线| 国产成+人综合+亚洲专区| 一区二区三区四区激情视频 | 亚洲国产高清在线一区二区三| 深爱激情五月婷婷| 久久久久久久久久成人|