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

    基于時(shí)間調(diào)制的超聲波自主實(shí)時(shí)定位系統(tǒng)

    2016-06-06 01:06:37林舟杰陸澤櫞衷衛(wèi)聲

    熊 劍,林舟杰,陸澤櫞,郭 杭,衷衛(wèi)聲

    (1. 南昌大學(xué)信息工程學(xué)院測(cè)控教研室,江西 南昌 330031;2. 中國(guó)科學(xué)院合肥物質(zhì)科學(xué)研究院醫(yī)學(xué)物理中心,安徽 合肥 230031;3. 南昌大學(xué)空間科學(xué)與技術(shù)研究院,江西 南昌 330031)

    ?

    基于時(shí)間調(diào)制的超聲波自主實(shí)時(shí)定位系統(tǒng)

    熊劍1,林舟杰1,陸澤櫞2,郭杭3,衷衛(wèi)聲1

    (1. 南昌大學(xué)信息工程學(xué)院測(cè)控教研室,江西 南昌330031;2. 中國(guó)科學(xué)院合肥物質(zhì)科學(xué)研究院醫(yī)學(xué)物理中心,安徽 合肥230031;3. 南昌大學(xué)空間科學(xué)與技術(shù)研究院,江西 南昌330031)

    摘要:針對(duì)超聲波-射頻室內(nèi)定位系統(tǒng)信道分配問題復(fù)雜、實(shí)時(shí)性低下的問題,提出基于時(shí)間調(diào)制的超聲波自主實(shí)時(shí)定位系統(tǒng)。該系統(tǒng)引入直接序列擴(kuò)頻技術(shù),用時(shí)間和終端ID組成的二進(jìn)制序列調(diào)制超聲波實(shí)現(xiàn)信息攜帶,接收端在無射頻模塊的情況下通過恢復(fù)信息來估計(jì)超聲波渡越時(shí)間,識(shí)別發(fā)射源并計(jì)算位置坐標(biāo)。通過simulink對(duì)超聲波信息攜帶方法進(jìn)行了仿真驗(yàn)證,實(shí)驗(yàn)結(jié)果表明:系統(tǒng)通過擴(kuò)頻調(diào)制解調(diào)能夠準(zhǔn)確恢復(fù)時(shí)間信息,有效識(shí)別終端ID,區(qū)分不同終端信號(hào)源位置,并且對(duì)多址效應(yīng)和多路徑效應(yīng)有較好的抗干擾效果,從而降低系統(tǒng)復(fù)雜性,提高響應(yīng)時(shí)間并實(shí)現(xiàn)自主定位。

    關(guān)鍵詞:超聲波;時(shí)間調(diào)制;擴(kuò)頻;自主定位

    0引言

    在現(xiàn)代大型商場(chǎng)、高度自動(dòng)化工廠或不適合人員直接操作的室內(nèi)工作環(huán)境,常常需要借助室內(nèi)定位技術(shù)獲得目標(biāo)的精確位置?;诔暡ǖ奶綔y(cè)技術(shù)安裝使用方便,成本低,抗電磁干擾能力強(qiáng),在測(cè)距中應(yīng)用廣泛,能獲得較高的精度[1-2]。在室內(nèi)定位技術(shù)領(lǐng)域中,由于超聲波以上優(yōu)點(diǎn),也具有較大的研究意義。

    目前的超聲波室內(nèi)定位系統(tǒng)大多數(shù)都是以超聲-射頻組合的形式進(jìn)行定位[3-4]。需要借助射頻模塊控制超聲波的發(fā)射,命令接收模塊計(jì)時(shí)來估計(jì)超聲信號(hào)的渡越時(shí)間。為了避免信號(hào)的互相串?dāng)_,射頻模塊還要控制信道的分配,不僅給系統(tǒng)帶來了較高的復(fù)雜性,其定位實(shí)時(shí)性也大打折扣[5]。針對(duì)此問題,本文提出了基于時(shí)間調(diào)制的超聲波自主實(shí)時(shí)定位系統(tǒng),在不借助無線射頻模塊的情況下實(shí)現(xiàn)定位,降低系統(tǒng)復(fù)雜性,提高定位響應(yīng)時(shí)間。

    1常規(guī)超聲波定位原理

    1.1Cricket定位系統(tǒng)

    Cricket是經(jīng)典的基于超聲波和無線射頻的室內(nèi)定位系統(tǒng)[6]。目前室內(nèi)超聲波定位系統(tǒng)結(jié)構(gòu)大多與此相同。其模塊節(jié)點(diǎn)分為L(zhǎng)istener和Beacon。Beacon被安裝在固定位置,同時(shí)發(fā)送射頻信號(hào)和超聲波脈沖;Listener安裝在移動(dòng)目標(biāo)上,接收到Beacon的射頻信號(hào)開始計(jì)時(shí),等待超聲波信號(hào)到達(dá)則停止計(jì)時(shí)。設(shè)這兩種信號(hào)的時(shí)間差為Δt,則

    s=v·Δt

    (1)

    獲得三個(gè)以上Beacon信標(biāo)時(shí)計(jì)算出Listener和信標(biāo)的距離s1,s2,s3,隨后根據(jù)三邊定位等定位算法即可計(jì)算出移動(dòng)目標(biāo)位置[7]。

    1.2信道分配和實(shí)時(shí)性問題

    由于所有Beacon均發(fā)射同頻超聲信號(hào),Listener無法區(qū)分不同Beacon信號(hào)來源,從而可能計(jì)算出錯(cuò)誤的位置。因此,Circket采取了隨機(jī)延遲技術(shù),即在信號(hào)發(fā)射之前先監(jiān)聽一段隨機(jī)時(shí)間T,如果沒有收到其他Beacon射頻信號(hào)則開始發(fā)射超聲波信號(hào)。這種監(jiān)聽機(jī)制降低了超聲信號(hào)發(fā)射的頻率,雖然減小了因超聲信號(hào)串?dāng)_導(dǎo)致目標(biāo)位置計(jì)算出錯(cuò)的可能性,但是難以保證響應(yīng)時(shí)間,此外,由于超聲回波信號(hào)的反射等引起的多路徑效應(yīng),仍然會(huì)造成距離測(cè)量錯(cuò)誤,得出錯(cuò)誤坐標(biāo)。

    1.3基于超聲射頻定位系統(tǒng)的實(shí)時(shí)性改進(jìn)方法

    為了改進(jìn)超聲-射頻定位系統(tǒng)的實(shí)時(shí)性,本文利用調(diào)制與解調(diào)和直接序列擴(kuò)頻技術(shù)[8],令發(fā)射的超聲波具備信息攜帶能力,將時(shí)間信息調(diào)制到超聲信號(hào)中,接收裝置通過解調(diào)恢復(fù)發(fā)射時(shí)間,與到達(dá)時(shí)間作差,無須借助射頻模塊即可完成渡越時(shí)間估計(jì),實(shí)現(xiàn)距離測(cè)量和目標(biāo)位置的計(jì)算。此外,信號(hào)還攜帶自身ID號(hào)碼,通過ID識(shí)別可區(qū)分信號(hào)來源,與Cricket原理類似的定位系統(tǒng)則可以免去信道分配問題,有利于改善系統(tǒng)定位的響應(yīng)時(shí)間,提高定位動(dòng)態(tài)效率。

    1.4m序列擴(kuò)頻原理

    研究表明,偽隨機(jī)二進(jìn)制序列最適于超聲發(fā)射的脈沖壓縮信號(hào),其中使用最多的就是m序列[9-10]。本系統(tǒng)采用m序列實(shí)現(xiàn)超聲信號(hào)的擴(kuò)頻過程。

    m序列全稱是最長(zhǎng)線性反饋移位寄存器序列,是偽隨機(jī)序列中最重要的一種,具有偽噪聲碼的一般特性。用x表示每一級(jí)移位寄存器,系數(shù)Ci表示寄存器的反饋狀態(tài),則其相應(yīng)的反饋線連接狀態(tài)可用多項(xiàng)式表示為:

    f(x)=C0x0+C1x1+…+Cnxn

    (2)

    式中C0=1,Cn=1。當(dāng)上式滿足條件[11]

    1)f(x)為不可約的;

    2)f(x)可整除(1+xp),p = 2n-1;

    3)f(x)除不盡(1+xn),q < p。

    f(x)稱為m序列的本原多項(xiàng)式。本原多項(xiàng)式具體求法不作贅述,參考相關(guān)文獻(xiàn)可查詢常用多項(xiàng)式系數(shù)[12]。

    查詢的系數(shù)多為八進(jìn)制,求出對(duì)應(yīng)二進(jìn)制系數(shù)即可確定本原多項(xiàng)式。

    根據(jù)m序列的平移等價(jià)特性和平衡特性,其歸一化自相關(guān)函數(shù)為:

    (3)

    式(3)中p 為m序列的長(zhǎng)度,j為相關(guān)函數(shù)R(j)的時(shí)延,j取0,±1,±2,…等離散值[13]。

    由式(3)可知,接收端相關(guān)器使用本地m序列與接收到的超聲波信號(hào)作自相關(guān)運(yùn)算,當(dāng)本地m序列相位相同并且同步時(shí),在發(fā)射端的鍵控作用便被相互抵消,輸出一個(gè)解擴(kuò)后的窄帶信號(hào)。

    2超聲波自主識(shí)別定位系統(tǒng)

    2.1系統(tǒng)結(jié)構(gòu)

    超聲波自主識(shí)別定位系統(tǒng)結(jié)構(gòu)如圖1所示,系統(tǒng)僅由超聲波發(fā)射終端、接收終端組成,在系統(tǒng)運(yùn)行期間,兩者都使用高精度時(shí)鐘模塊進(jìn)行計(jì)時(shí)。在系統(tǒng)啟動(dòng)前,首先對(duì)發(fā)射和接收終端進(jìn)行時(shí)鐘同步校準(zhǔn)。

    發(fā)射終端:固定在天花板上,構(gòu)成一個(gè)坐標(biāo)系,并對(duì)信號(hào)發(fā)射時(shí)間和終端ID號(hào)進(jìn)行調(diào)制,并持續(xù)發(fā)射超聲波信號(hào)。

    接收終端:固定在定位物體之上,接收來自發(fā)射終端的超聲波信號(hào)并進(jìn)行數(shù)據(jù)處理和定位解算。

    圖1 超聲波自主識(shí)別定位系統(tǒng)Fig.1 Autonomously recognizable ultrasonic position system

    2.2信號(hào)傳輸原理

    系統(tǒng)工作流程如圖2所示。

    圖2 基于時(shí)間調(diào)制的超聲室內(nèi)定位方法流程Fig.2 The process of ultrasonic indoor position based on time modulation

    發(fā)射終端準(zhǔn)備發(fā)射超聲波,同時(shí)記錄此時(shí)發(fā)射時(shí)刻t1,使用偽隨機(jī)碼拓展信號(hào)頻譜,把時(shí)間信息t1和終端ID號(hào)n(用d(t)表示)調(diào)制在偽隨機(jī)序列當(dāng)中,原始數(shù)據(jù)d(t)與m序列c(t)相乘得到擴(kuò)頻后的復(fù)合信號(hào)

    g(t)=d(t)c(t)

    (4)

    使用復(fù)合信號(hào)g(t)調(diào)制超聲載波信號(hào),即當(dāng)原始數(shù)據(jù)二進(jìn)制序列為0時(shí)發(fā)送與m序列反相的序列,為1時(shí)發(fā)射與m序列同相的序列,得到調(diào)制信號(hào)s(t),其時(shí)域表達(dá)式可表示為:

    (5)

    式(5)中,(n-1)TB

    (6)

    完成信號(hào)擴(kuò)頻調(diào)制后,通過發(fā)射終端的超聲波驅(qū)動(dòng)電路發(fā)射經(jīng)擴(kuò)頻調(diào)制后的超聲波信號(hào)s(t)。

    圖3 發(fā)射端擴(kuò)頻調(diào)制框圖Fig.3 Spread spectrum modulation simulation of sender

    超聲波接收終端接收到任意一個(gè)發(fā)射終端的超聲波信號(hào),立刻記錄本地接收時(shí)刻t2;與此同時(shí),依照?qǐng)D4,該接收終端使用本地m序列與收到的調(diào)制信號(hào)s(t)作相關(guān)解擴(kuò),當(dāng)本地與發(fā)射端m序列結(jié)構(gòu)完全相同,即c(t)與c'(t)具有相同的周期,碼元同步,并且相位重合時(shí),由于m序列良好的自相關(guān)特性,根據(jù)式(5)可得

    c(t)·c′(t)=1

    (7)

    圖4 超聲波接收端解擴(kuò)解調(diào)框圖Fig.4 Dispreading modulation simulation of receiver

    接收終端根據(jù)超聲波信號(hào)發(fā)射時(shí)刻t1以及接收s(t)的時(shí)刻t2,以及公式

    (8)

    計(jì)算出該接收終端到n號(hào)發(fā)射終端的距離s。

    發(fā)射終端理論上可連續(xù)不斷地發(fā)射超聲擴(kuò)頻信號(hào),目標(biāo)接收終端可以隨時(shí)接收不同發(fā)射終端發(fā)出的信號(hào),只要同時(shí)收到3個(gè)有效信號(hào),經(jīng)過信號(hào)解擴(kuò)和解調(diào),即可根據(jù)三維定位算法[4],解算出目標(biāo)接收終端的位置。同時(shí)識(shí)別不同發(fā)射終端源,解決了多個(gè)發(fā)射終端帶來的多址干擾和超聲波信號(hào)反射帶來的多徑干擾問題。

    3系統(tǒng)實(shí)驗(yàn)與結(jié)果分析

    3.1超聲時(shí)間信息及ID數(shù)據(jù)傳輸系統(tǒng)的仿真

    以3路超聲波定位系統(tǒng)為例,采用Simulink對(duì)超聲時(shí)間信息及ID數(shù)據(jù)傳輸系統(tǒng)進(jìn)行了系統(tǒng)仿真和分析。該系統(tǒng)仿真如圖5所示。

    首先對(duì)每一路的超聲波發(fā)射時(shí)刻跟發(fā)射終端ID進(jìn)行編碼,編碼數(shù)據(jù)仿真以隨機(jī)二進(jìn)制序列生成器表示。

    系統(tǒng)仿真主要由擴(kuò)頻調(diào)制、信息調(diào)制、相關(guān)解擴(kuò)、信息解調(diào)4部分組成,其中信息調(diào)制與解調(diào)采用BPSK二進(jìn)制相移鍵控方式。仿真中針對(duì)工作頻率為40 kHz的超聲波換能器,使用頻率為40 kHz,振幅5 V的正弦波信號(hào)。

    圖5 超聲時(shí)間及ID傳輸-識(shí)別系統(tǒng)仿真Fig.5 Time & ID transfer-recognition system simulation

    3路超聲波信號(hào)經(jīng)過擴(kuò)頻調(diào)制后通過同一個(gè)高斯噪聲信道AWNG Channel,模擬多個(gè)超聲波發(fā)射終端帶來的多址干擾,同時(shí)對(duì)第一路信號(hào)進(jìn)行二次疊加模擬超聲波信號(hào)經(jīng)反射引起的多路徑效應(yīng)。

    圖5中擴(kuò)頻模塊和解擴(kuò)模塊對(duì)應(yīng)的仿真見圖6。

    圖6 擴(kuò)頻與解擴(kuò)子模塊Fig.6 Sub-module of spread and dispreading

    由于采用3路信號(hào),偽隨機(jī)碼選擇5級(jí)m序列即可滿足系統(tǒng)要求。根據(jù)表1所列,求出三個(gè)本原多項(xiàng)式二進(jìn)制系數(shù)分別為[1 0 1 0 0 1],[1 1 1 0 1 1],[1 0 1 1 1 1],設(shè)置每一路擴(kuò)頻子模塊中PN Sequence Generator模塊的參數(shù)Generator polynomial,解擴(kuò)時(shí)本地PN Sequence Generator使用參數(shù)設(shè)置一致的Generator polynomial。

    3.2仿真結(jié)果

    為了比較效果,實(shí)驗(yàn)對(duì)未引入DSSS技術(shù)的單路傳輸系統(tǒng)進(jìn)行驗(yàn)證,見圖7,其波形如圖8所示。由圖8中波形2對(duì)比波形1原始數(shù)據(jù)可見,單路超聲傳輸系統(tǒng)尚能準(zhǔn)確地恢復(fù)超聲波發(fā)射端的時(shí)間信息和ID數(shù)據(jù),結(jié)果較為準(zhǔn)確。加入多路徑效應(yīng)后,波形3與波形1對(duì)比可見,無DSSS的傳輸系統(tǒng)受到較大影響,接收端已無法準(zhǔn)確恢復(fù)發(fā)射端的數(shù)據(jù)。

    在三路超聲波系統(tǒng)共同工作時(shí),超聲波信號(hào)互相干擾,其中第一路收發(fā)波形如圖8波形4所示,接收端數(shù)據(jù)已經(jīng)無法完全恢復(fù),系統(tǒng)受到嚴(yán)重的多址干擾和多路徑干擾,在多路超聲信號(hào)串?dāng)_下無法準(zhǔn)確恢復(fù)數(shù)據(jù),從而會(huì)得出錯(cuò)誤目標(biāo)位置。

    圖9為引入DSSS技術(shù)的三路超聲時(shí)間信息及終端ID數(shù)據(jù)傳輸系統(tǒng)的仿真結(jié)果圖。三路發(fā)射終端發(fā)射不同的二進(jìn)制序列,分別對(duì)比各路發(fā)射終端可發(fā)現(xiàn),三路發(fā)射端時(shí)間信息和ID號(hào)構(gòu)成的二進(jìn)制序列經(jīng)過擴(kuò)頻調(diào)制以后,在接收端均成功地得到恢復(fù),解決了同頻信號(hào)串?dāng)_,信號(hào)反射帶來的多路徑干擾,順利實(shí)現(xiàn)了超聲波信息攜帶功能,從而為超聲波定位系統(tǒng)無間斷連續(xù)定位提供了理論基礎(chǔ),有利于提高系統(tǒng)的定位響應(yīng)時(shí)間。

    圖7 無DSSS的超聲時(shí)間及ID傳輸-識(shí)別系統(tǒng)仿真Fig.7 Time & ID transfer-recognition system simulation without DSSS

    圖8 無DSSS收發(fā)波形圖Fig.8 Waveform between sender and receiver without DSSS

    圖9 引入DSSS的收發(fā)波形圖Fig.9 Waveform between sender and receiver with DSSS

    4結(jié)論

    本文提出了基于時(shí)間調(diào)制的超聲波自主實(shí)時(shí)定位系統(tǒng)。該系統(tǒng)使用m序列對(duì)超聲波發(fā)射端數(shù)據(jù)進(jìn)行擴(kuò)頻調(diào)制,通過傳輸超聲波發(fā)射時(shí)間信息和終端ID號(hào)以實(shí)現(xiàn)自主定位。仿真結(jié)果證實(shí)了超聲波攜帶信息的可行性,能準(zhǔn)確恢復(fù)時(shí)間信息并識(shí)別終端ID,區(qū)分不同終端信號(hào)源位置,而且具有良好的抗多址干擾和抗多路徑干擾效果。系統(tǒng)能夠不借助無線模塊實(shí)現(xiàn)超聲波自主定位,減小定位系統(tǒng)的復(fù)雜性,有利于改善超聲波定位系統(tǒng)的動(dòng)態(tài)響應(yīng),對(duì)于提高超聲波定位系統(tǒng)實(shí)時(shí)性具備一定的研究和應(yīng)用價(jià)值。

    隨著時(shí)間的累積,超聲波發(fā)射和接收終端的時(shí)鐘誤差會(huì)逐漸增大,這部分誤差可以借鑒GPS鐘差,通過超聲波傳輸導(dǎo)航電文來消除,這也將是后續(xù)研究工作的重要內(nèi)容,有待進(jìn)一步研究。

    參考文獻(xiàn):

    [1]Liu H, Darabi H, Banerjee P, et al. Survey of wireless indoor positioning techniques and systems[J]. Systems, Man, and Cybernetics, Part C: Applications and Reviews, IEEE Transactions, 2007, 37(6): 1067-1080

    [2]趙軍, 李鴻斌, 王智. 無線網(wǎng)絡(luò)室內(nèi)定位系統(tǒng)研究[J]. 信息與控制, 2008, 37(4):465-471.

    [3]Koyuncu H, Yang S H. A survey of indoor positioning and object locating systems[J]. IJCSNS International Journal of Computer Science and Network Security, 2010, 10(5): 121-128.

    [4]韓霜,羅海勇,陳穎,等. 基于TDOA的超聲波室內(nèi)定位系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J]. 傳感技術(shù)學(xué)報(bào),2010(3):347-353.

    [5]韓剛. 基于無線電和超聲波的室內(nèi)定位技術(shù)響應(yīng)時(shí)間的研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2009.

    [6]Priyantha N B. The cricket indoor location system[D]. Massachusetts Institute of Technology, 2005.

    [7]王鴻鵬,王耀寬. 基于Cricket傳感器網(wǎng)絡(luò)室內(nèi)定位系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J]. 計(jì)算機(jī)工程與應(yīng)用,2008(2):211-214.

    [8]Cho H, Kim S W. Mobile robot localization using biased chirp-spread-spectrum ranging[J]. Industrial Electronics, IEEE Transactions on, 2010, 57(8): 2826-2835.

    [9]程曉暢, 蘇紹景, 王躍科,等. 偽隨機(jī)碼超聲擴(kuò)頻測(cè)距系統(tǒng)設(shè)計(jì)與算法[J]. 測(cè)試技術(shù)學(xué)報(bào), 2007, 21(1): 79-83.

    [10]蘇紹璟,伍文君,黃芝平,等.含錯(cuò)m序列本原多項(xiàng)式的高階統(tǒng)計(jì)測(cè)定算法[J]. 兵工學(xué)報(bào),2010,12:1593-1598.

    [11]Li J, Zeng X, Tang X, et al. A family of quadriphase sequences of period 4 (2 n-1) with low correlation and large linear span[J]. Designs, codes and cryptography, 2013, 67(1): 19-35.

    [12]Tesoriero R, Tebar R, Gallud J A, et al. Improving location awareness in indoor spaces using RFID technology[J]. Expert Systems with Applications, 2010, 37(1): 894-898.

    [13]Enrique González, Laura Prados, Antonio J. Rubio, et al. ATLINTIDA: A Robust Indoor Ultrasound Location System: Design and Evaluation[J]. Advances in Soft Computing, 2009: 180-190.

    Autonomous Real-time Ultrasonic Position System Based on Time-modulation

    XIONG Jian1, LIN Zhoujie1, LU Zeyuan2,GUO Hang3, ZHONG Weisheng1

    (1. Department of Measuring & Control Engineering, Academy of Information Engineering,Nanchang University,Nanchang 330031, China; 2. Center of Medical Physics and Technology,Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China;3. Institute of Space Science and Technology, Nanchang University, Nanchang 330031, China)

    Abstract:For the problem of complexity of channel allocation and low real-time based on ultrasonic-RF indoor localization system, an independent real-time ultrasonic positioning system based on time modulation was put forward. The system applied direct sequence spread spectrum technology, binary sequence consisted of time and terminal ID modulate ultrasonic to implement information-carrying. The receiver estimated ultrasonic transit time, identified and calculated the position of the transmission source coordinates by restoring the information without RF module. Simulation experiment results showed that the system could accurately restored time information by spreading modulation and demodulation, valid identification terminal ID, to distinguish between different terminals source position and it had good effect on multi-access and multipath interference, thereby reducing system complexity, improving response time and achieve self-positioning.

    Key words:ultrasound; time modulation; spread spectrum; autonomous positioning.

    中圖分類號(hào):TP96

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

    文章編號(hào):1008-1194(2016)02-0043-05

    作者簡(jiǎn)介:熊劍(1977—),男,江西省撫州人,博士,講師,研究方向:多源信息融合、非線性濾波。E-mail:xiongjian@ncu.edu.cn。

    基金項(xiàng)目:江西省教育廳科技項(xiàng)目資助(GJJ14125;GJJ13056);國(guó)家自然科學(xué)基金項(xiàng)目資助(41164001;41374039)

    *收稿日期:2015-10-26

    如何舔出高潮| 成人漫画全彩无遮挡| 国产av码专区亚洲av| 欧美精品高潮呻吟av久久| 国产精品免费视频内射| 久久人人爽av亚洲精品天堂| 成人国语在线视频| 性高湖久久久久久久久免费观看| 久久久久久久亚洲中文字幕| 国产精品久久久久久精品电影小说| 亚洲精品国产色婷婷电影| 免费播放大片免费观看视频在线观看| 精品福利永久在线观看| 国产精品国产av在线观看| 久久人人97超碰香蕉20202| 亚洲精品日本国产第一区| 人人澡人人妻人| 99久久中文字幕三级久久日本| 亚洲精品日韩在线中文字幕| videosex国产| 美女中出高潮动态图| 自线自在国产av| a 毛片基地| 欧美成人午夜精品| 国产精品人妻久久久影院| 国产福利在线免费观看视频| 在线观看一区二区三区激情| 亚洲国产最新在线播放| 日韩不卡一区二区三区视频在线| 国产探花极品一区二区| 蜜桃国产av成人99| 国产亚洲午夜精品一区二区久久| 夫妻性生交免费视频一级片| 自线自在国产av| 国产成人a∨麻豆精品| 免费观看a级毛片全部| 18+在线观看网站| 99热全是精品| 人人妻人人添人人爽欧美一区卜| 丝袜人妻中文字幕| 18禁国产床啪视频网站| 亚洲人成网站在线观看播放| 国产亚洲av片在线观看秒播厂| 精品人妻一区二区三区麻豆| 国产xxxxx性猛交| 女人被躁到高潮嗷嗷叫费观| 香蕉精品网在线| 少妇 在线观看| 天天躁日日躁夜夜躁夜夜| 久久久久久久精品精品| 一级毛片电影观看| 亚洲内射少妇av| 日本vs欧美在线观看视频| 视频区图区小说| 69精品国产乱码久久久| 欧美精品一区二区免费开放| 国产精品 国内视频| 交换朋友夫妻互换小说| 欧美xxⅹ黑人| 波野结衣二区三区在线| 亚洲中文av在线| 1024香蕉在线观看| 日韩人妻精品一区2区三区| 男女无遮挡免费网站观看| 国产野战对白在线观看| 18禁动态无遮挡网站| 丰满饥渴人妻一区二区三| 一边亲一边摸免费视频| 亚洲av成人精品一二三区| 日韩av不卡免费在线播放| 免费高清在线观看视频在线观看| 欧美日韩精品成人综合77777| 观看av在线不卡| 国产极品粉嫩免费观看在线| 少妇的丰满在线观看| 国产成人91sexporn| 国产精品无大码| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区黑人 | 在线看a的网站| 最近中文字幕2019免费版| av有码第一页| 国产一区二区三区综合在线观看| 欧美变态另类bdsm刘玥| 亚洲第一av免费看| 亚洲少妇的诱惑av| 国产成人欧美| 最近的中文字幕免费完整| 欧美日韩精品网址| av不卡在线播放| 国语对白做爰xxxⅹ性视频网站| 日日爽夜夜爽网站| 搡老乐熟女国产| 日韩大片免费观看网站| 国产国语露脸激情在线看| 日产精品乱码卡一卡2卡三| 亚洲激情五月婷婷啪啪| 欧美日韩综合久久久久久| 女人久久www免费人成看片| 午夜福利,免费看| 欧美激情 高清一区二区三区| 水蜜桃什么品种好| 夜夜骑夜夜射夜夜干| 亚洲婷婷狠狠爱综合网| 少妇猛男粗大的猛烈进出视频| 精品第一国产精品| 男女边吃奶边做爰视频| 色吧在线观看| 亚洲伊人久久精品综合| 久久久久精品久久久久真实原创| 亚洲伊人久久精品综合| 香蕉国产在线看| 寂寞人妻少妇视频99o| 久久亚洲国产成人精品v| 国产xxxxx性猛交| 熟妇人妻不卡中文字幕| 天堂俺去俺来也www色官网| 久久综合国产亚洲精品| 青草久久国产| 老司机影院毛片| 啦啦啦啦在线视频资源| 久久午夜综合久久蜜桃| 一边亲一边摸免费视频| 黄片播放在线免费| 久久精品久久精品一区二区三区| 午夜免费鲁丝| 午夜福利视频精品| 最近手机中文字幕大全| 一级a爱视频在线免费观看| 韩国av在线不卡| 久久人人97超碰香蕉20202| 免费少妇av软件| 国产淫语在线视频| 成年动漫av网址| 老汉色av国产亚洲站长工具| 美女脱内裤让男人舔精品视频| 国产精品一国产av| 18在线观看网站| 欧美精品一区二区免费开放| 亚洲美女视频黄频| 成人国产av品久久久| 亚洲精品国产av成人精品| 婷婷色综合www| 亚洲欧美日韩另类电影网站| 女人被躁到高潮嗷嗷叫费观| 亚洲图色成人| 午夜福利乱码中文字幕| 欧美人与善性xxx| 中文字幕人妻丝袜一区二区 | 伦精品一区二区三区| 老司机影院毛片| 亚洲精品自拍成人| 波多野结衣一区麻豆| 26uuu在线亚洲综合色| 精品一区二区免费观看| 在线精品无人区一区二区三| www.av在线官网国产| 欧美激情高清一区二区三区 | 亚洲精品久久久久久婷婷小说| 成人国语在线视频| 少妇猛男粗大的猛烈进出视频| 黄色视频在线播放观看不卡| 老司机亚洲免费影院| 国产精品嫩草影院av在线观看| 在现免费观看毛片| 一二三四在线观看免费中文在| 最近的中文字幕免费完整| 啦啦啦啦在线视频资源| 亚洲国产成人一精品久久久| 自线自在国产av| 在线观看人妻少妇| 亚洲精品久久成人aⅴ小说| 80岁老熟妇乱子伦牲交| 国产精品二区激情视频| 国产免费现黄频在线看| 免费观看av网站的网址| 日韩人妻精品一区2区三区| 日本av免费视频播放| 亚洲av成人精品一二三区| 一级毛片 在线播放| 中文精品一卡2卡3卡4更新| www.精华液| 国产精品 欧美亚洲| 久久午夜综合久久蜜桃| 黄色毛片三级朝国网站| 超碰成人久久| 免费黄色在线免费观看| 国产精品一二三区在线看| 女人被躁到高潮嗷嗷叫费观| 久久精品国产亚洲av天美| 涩涩av久久男人的天堂| 欧美日韩精品网址| 免费黄色在线免费观看| 飞空精品影院首页| 久久久a久久爽久久v久久| 午夜福利一区二区在线看| 深夜精品福利| 久久久久久久久久久免费av| 肉色欧美久久久久久久蜜桃| 男的添女的下面高潮视频| 大香蕉久久成人网| 中国三级夫妇交换| 各种免费的搞黄视频| 国产成人91sexporn| 日本wwww免费看| 少妇人妻久久综合中文| 亚洲av.av天堂| 日日撸夜夜添| 搡女人真爽免费视频火全软件| 国产精品秋霞免费鲁丝片| 看非洲黑人一级黄片| 国产极品粉嫩免费观看在线| 又大又黄又爽视频免费| 免费播放大片免费观看视频在线观看| 精品午夜福利在线看| 大片电影免费在线观看免费| 少妇的丰满在线观看| 久久97久久精品| 丝袜脚勾引网站| 99久久人妻综合| 亚洲欧美成人精品一区二区| 免费高清在线观看日韩| 久久精品亚洲av国产电影网| 少妇的丰满在线观看| 天堂8中文在线网| 九色亚洲精品在线播放| 国产爽快片一区二区三区| 人妻 亚洲 视频| 人人妻人人澡人人爽人人夜夜| 亚洲精品久久久久久婷婷小说| 亚洲国产欧美在线一区| 王馨瑶露胸无遮挡在线观看| 999精品在线视频| 少妇的逼水好多| 黑丝袜美女国产一区| 国产精品成人在线| 欧美日韩亚洲国产一区二区在线观看 | 国精品久久久久久国模美| 一级毛片我不卡| 欧美日韩亚洲高清精品| 蜜桃在线观看..| 久久人妻熟女aⅴ| 免费播放大片免费观看视频在线观看| 成人18禁高潮啪啪吃奶动态图| 天天躁日日躁夜夜躁夜夜| 久久久久久久大尺度免费视频| 另类亚洲欧美激情| 国产成人精品无人区| 又大又黄又爽视频免费| 国产成人av激情在线播放| 人妻少妇偷人精品九色| 国产亚洲精品第一综合不卡| 国产精品一区二区在线观看99| 成人黄色视频免费在线看| 色哟哟·www| 美女国产高潮福利片在线看| 久久精品aⅴ一区二区三区四区 | 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 韩国精品一区二区三区| 十八禁高潮呻吟视频| 少妇被粗大的猛进出69影院| 波野结衣二区三区在线| 青草久久国产| 亚洲内射少妇av| 午夜91福利影院| 亚洲精品在线美女| 久久婷婷青草| 午夜福利在线免费观看网站| 久久99精品国语久久久| 国产成人精品婷婷| av国产精品久久久久影院| av女优亚洲男人天堂| 高清在线视频一区二区三区| 99re6热这里在线精品视频| 欧美精品人与动牲交sv欧美| 精品国产一区二区久久| 精品一区在线观看国产| 国产精品久久久久成人av| 香蕉精品网在线| kizo精华| 亚洲精品aⅴ在线观看| 亚洲美女视频黄频| 天堂8中文在线网| 一级a爱视频在线免费观看| 欧美av亚洲av综合av国产av | 免费在线观看完整版高清| 亚洲欧美精品自产自拍| 亚洲美女视频黄频| 久久久久国产一级毛片高清牌| 中文字幕av电影在线播放| 一级黄片播放器| 在线观看一区二区三区激情| 伦理电影大哥的女人| 两性夫妻黄色片| 伦理电影免费视频| 亚洲国产色片| 久热久热在线精品观看| 国产av一区二区精品久久| 成人毛片60女人毛片免费| 最黄视频免费看| 亚洲精品乱久久久久久| 啦啦啦在线观看免费高清www| 国精品久久久久久国模美| 9热在线视频观看99| 晚上一个人看的免费电影| 国产日韩欧美亚洲二区| 人人妻人人爽人人添夜夜欢视频| 精品国产一区二区三区久久久樱花| 免费人妻精品一区二区三区视频| 一本—道久久a久久精品蜜桃钙片| √禁漫天堂资源中文www| 嫩草影院入口| 欧美成人午夜免费资源| 黄色配什么色好看| 亚洲精品,欧美精品| 观看美女的网站| 欧美成人午夜精品| 亚洲,欧美精品.| 亚洲精品视频女| 亚洲av电影在线观看一区二区三区| 欧美成人午夜精品| 人人妻人人爽人人添夜夜欢视频| 两性夫妻黄色片| 亚洲精品久久午夜乱码| 国产精品久久久久久久久免| 黄色配什么色好看| 亚洲美女视频黄频| 最近手机中文字幕大全| 亚洲第一av免费看| 天堂中文最新版在线下载| 亚洲国产av影院在线观看| 国产男女内射视频| 少妇人妻精品综合一区二区| 国产精品国产av在线观看| 久久久久精品人妻al黑| 一区二区三区精品91| 色婷婷av一区二区三区视频| 精品第一国产精品| 蜜桃国产av成人99| 国产一区二区激情短视频 | 日产精品乱码卡一卡2卡三| 亚洲av电影在线观看一区二区三区| 亚洲精品成人av观看孕妇| 欧美老熟妇乱子伦牲交| 丝袜人妻中文字幕| 中国国产av一级| 亚洲成人一二三区av| 久久久久精品人妻al黑| 69精品国产乱码久久久| 亚洲精品av麻豆狂野| 国产成人欧美| 叶爱在线成人免费视频播放| 丰满乱子伦码专区| 成人二区视频| 九草在线视频观看| 菩萨蛮人人尽说江南好唐韦庄| 黄色怎么调成土黄色| 中文字幕制服av| 黑丝袜美女国产一区| 国产精品三级大全| 啦啦啦啦在线视频资源| 免费高清在线观看日韩| 精品午夜福利在线看| 美国免费a级毛片| av在线老鸭窝| 亚洲av日韩在线播放| 久久久久视频综合| 久久久久网色| 80岁老熟妇乱子伦牲交| 国产精品国产三级专区第一集| 久久久久人妻精品一区果冻| 久久久精品区二区三区| 捣出白浆h1v1| 亚洲欧美成人综合另类久久久| 久久国产亚洲av麻豆专区| 99久久人妻综合| 香蕉丝袜av| 各种免费的搞黄视频| 在线 av 中文字幕| 97在线人人人人妻| 国产成人精品久久久久久| 电影成人av| 国产成人精品婷婷| 日韩 亚洲 欧美在线| 熟女av电影| 欧美日韩精品成人综合77777| 高清不卡的av网站| 新久久久久国产一级毛片| 亚洲欧美精品自产自拍| 色哟哟·www| 性色av一级| 看免费成人av毛片| 日韩伦理黄色片| 丝袜美足系列| 亚洲第一av免费看| 日韩av不卡免费在线播放| 美女xxoo啪啪120秒动态图| 晚上一个人看的免费电影| 亚洲av电影在线进入| 精品午夜福利在线看| 少妇人妻精品综合一区二区| 午夜福利,免费看| 久久免费观看电影| 久久ye,这里只有精品| 黄色 视频免费看| 色网站视频免费| 亚洲精品国产一区二区精华液| 一边摸一边做爽爽视频免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩中字成人| 国产免费一区二区三区四区乱码| 午夜影院在线不卡| 亚洲精品日韩在线中文字幕| 久久99热这里只频精品6学生| 欧美成人午夜精品| 国产成人一区二区在线| 波多野结衣一区麻豆| 国产成人免费无遮挡视频| 18禁观看日本| 色94色欧美一区二区| 女人精品久久久久毛片| 超碰成人久久| 久久这里有精品视频免费| 视频在线观看一区二区三区| 久久热在线av| 久热久热在线精品观看| 色94色欧美一区二区| 色吧在线观看| 女性被躁到高潮视频| 在线观看免费视频网站a站| 日韩制服骚丝袜av| 三级国产精品片| 美女高潮到喷水免费观看| 久久久精品94久久精品| 欧美成人精品欧美一级黄| 中文字幕人妻丝袜制服| 中文字幕最新亚洲高清| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久二区二区91 | 久久99一区二区三区| 国产高清国产精品国产三级| 韩国精品一区二区三区| av网站免费在线观看视频| av国产精品久久久久影院| 国产日韩欧美在线精品| 黄片无遮挡物在线观看| 国产在视频线精品| 久久毛片免费看一区二区三区| 欧美日韩精品成人综合77777| 2021少妇久久久久久久久久久| 国产不卡av网站在线观看| 亚洲国产av影院在线观看| 亚洲欧美精品综合一区二区三区 | 永久免费av网站大全| 免费观看性生交大片5| av在线老鸭窝| 免费高清在线观看视频在线观看| 99九九在线精品视频| 亚洲视频免费观看视频| 日本wwww免费看| 最近的中文字幕免费完整| 人人妻人人澡人人爽人人夜夜| 熟女少妇亚洲综合色aaa.| 美女视频免费永久观看网站| 久久久a久久爽久久v久久| 伦理电影大哥的女人| 成人漫画全彩无遮挡| 日韩电影二区| 日韩一卡2卡3卡4卡2021年| 人妻系列 视频| 亚洲第一区二区三区不卡| 美女福利国产在线| 午夜福利,免费看| 国产亚洲午夜精品一区二区久久| 精品人妻一区二区三区麻豆| 又大又黄又爽视频免费| 久久精品久久久久久久性| 边亲边吃奶的免费视频| 免费在线观看视频国产中文字幕亚洲 | 在线观看人妻少妇| 哪个播放器可以免费观看大片| 最近2019中文字幕mv第一页| 免费观看在线日韩| 国产精品女同一区二区软件| 亚洲国产精品国产精品| 久久久亚洲精品成人影院| videossex国产| 深夜精品福利| 一边亲一边摸免费视频| 久久精品国产鲁丝片午夜精品| 亚洲一区二区三区欧美精品| 国产在视频线精品| 日韩,欧美,国产一区二区三区| 啦啦啦在线观看免费高清www| 久久久久久久国产电影| 久久久久久久久久久久大奶| 久久午夜福利片| 亚洲第一青青草原| 在线观看三级黄色| 亚洲美女搞黄在线观看| 亚洲精品视频女| 老司机亚洲免费影院| 午夜福利在线免费观看网站| 人人澡人人妻人| 两个人免费观看高清视频| 七月丁香在线播放| 999精品在线视频| 日韩av不卡免费在线播放| 最近手机中文字幕大全| 五月开心婷婷网| 久久综合国产亚洲精品| 久久精品aⅴ一区二区三区四区 | 久久人人爽人人片av| 制服丝袜香蕉在线| 亚洲国产欧美网| 国产探花极品一区二区| 久久ye,这里只有精品| 精品国产一区二区三区四区第35| 亚洲美女黄色视频免费看| 亚洲欧美一区二区三区国产| 免费观看在线日韩| 9191精品国产免费久久| 中文字幕人妻丝袜一区二区 | 国产男女内射视频| 久久久国产欧美日韩av| 搡女人真爽免费视频火全软件| 亚洲av综合色区一区| freevideosex欧美| 丝袜脚勾引网站| 亚洲欧美精品自产自拍| 丰满少妇做爰视频| 久久精品人人爽人人爽视色| 国产人伦9x9x在线观看 | 久久久久久伊人网av| 天天躁夜夜躁狠狠躁躁| 成人黄色视频免费在线看| 国产精品国产三级国产专区5o| 成人毛片a级毛片在线播放| 免费在线观看视频国产中文字幕亚洲 | 麻豆精品久久久久久蜜桃| 伊人亚洲综合成人网| 欧美 亚洲 国产 日韩一| 国产精品av久久久久免费| 免费观看无遮挡的男女| 捣出白浆h1v1| 我的亚洲天堂| 亚洲国产精品999| 免费av中文字幕在线| 日韩精品免费视频一区二区三区| 日韩中文字幕视频在线看片| 国产亚洲av片在线观看秒播厂| 日韩制服骚丝袜av| 国产深夜福利视频在线观看| 日韩伦理黄色片| 久久久久久久大尺度免费视频| 两性夫妻黄色片| 激情视频va一区二区三区| 在线观看免费日韩欧美大片| 国产精品国产三级国产专区5o| 2022亚洲国产成人精品| 亚洲五月色婷婷综合| 亚洲激情五月婷婷啪啪| 97人妻天天添夜夜摸| 男的添女的下面高潮视频| 男人爽女人下面视频在线观看| 在线天堂最新版资源| 大片免费播放器 马上看| 国产白丝娇喘喷水9色精品| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| av不卡在线播放| 汤姆久久久久久久影院中文字幕| 欧美成人午夜精品| 美女主播在线视频| 亚洲少妇的诱惑av| 国产精品免费大片| 1024香蕉在线观看| 亚洲情色 制服丝袜| a级片在线免费高清观看视频| 国产精品免费大片| 午夜免费观看性视频| √禁漫天堂资源中文www| 人妻 亚洲 视频| www.熟女人妻精品国产| a级片在线免费高清观看视频| 人人澡人人妻人| 捣出白浆h1v1| 91精品国产国语对白视频| 亚洲经典国产精华液单| av国产久精品久网站免费入址| 人体艺术视频欧美日本| 热99国产精品久久久久久7| 国产成人精品久久二区二区91 | 国产激情久久老熟女| 精品人妻在线不人妻| 国产精品一区二区在线不卡| 午夜福利,免费看| 免费观看无遮挡的男女| 亚洲成色77777| 国产在视频线精品| 在线看a的网站| 秋霞在线观看毛片| 青草久久国产| 女的被弄到高潮叫床怎么办| 搡女人真爽免费视频火全软件| 极品少妇高潮喷水抽搐| 免费不卡的大黄色大毛片视频在线观看| 美女午夜性视频免费| av一本久久久久| 99国产精品免费福利视频| 久久 成人 亚洲| 叶爱在线成人免费视频播放| 国产精品熟女久久久久浪| 电影成人av|