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

    基于FPGA 的鐵磁物質(zhì)探測(cè)系統(tǒng)設(shè)計(jì)與研究

    2021-02-04 06:53:40
    軟件導(dǎo)刊 2021年1期
    關(guān)鍵詞:鐵磁磁通信號(hào)處理

    (上海理工大學(xué)醫(yī)療器械與食品學(xué)院,上海 200093)

    0 引言

    磁共振成像(Magnetic Resonance Imaging,MRI)[1]與目前臨床上應(yīng)用普遍的醫(yī)學(xué)影像檢查手段X 線CT 檢查相比,雖無(wú)電離輻射危害但設(shè)備所在物理空間有強(qiáng)磁場(chǎng)存在(0.2~3.0T),第一大安全隱患就是鐵磁性物質(zhì)誤入造成的設(shè)備損壞和人員傷亡事故。由于MRI 區(qū)域內(nèi)強(qiáng)磁環(huán)境影響,對(duì)于隱藏、植入、附著于患者/工作人員/訪客身體的危險(xiǎn)細(xì)小鐵磁金屬物體檢測(cè)更進(jìn)一步增加了難度[2]。

    基于磁通門(mén)傳感器原理的鐵磁探測(cè)系統(tǒng)[3]被廣泛應(yīng)用于航空探潛、海洋監(jiān)測(cè)、地下和水下鐵磁物體的探測(cè)、地震預(yù)測(cè)以及醫(yī)學(xué)核磁共振等各領(lǐng)域。

    磁通門(mén)信號(hào)常用檢測(cè)方法[4]有脈沖幅值法、脈沖間隔法等。對(duì)于脈沖幅值法[5],由于磁芯微觀磁性能、應(yīng)變應(yīng)力、熱騷動(dòng)等外界因素影響,存在一些很難計(jì)值的尖峰,會(huì)直接威脅脈沖幅值法的零輸入穩(wěn)定性和精確性。而對(duì)于脈沖間隔法[6],被測(cè)磁場(chǎng)的橫向分量與激勵(lì)磁場(chǎng)、被測(cè)磁場(chǎng)在探頭軸向分量的疊加作用使磁芯加深飽和度,干擾了磁通門(mén)信號(hào)脈沖間隔測(cè)量的準(zhǔn)確性。

    經(jīng)理論分析,任意偶次諧波都可以反映被測(cè)磁場(chǎng)的大小,而二次諧波[7]在各偶次諧波中幅值最大,故選取其作為被測(cè)磁場(chǎng)的度量。如果激磁信號(hào)是非正弦信號(hào),信號(hào)中就會(huì)存在偶次諧波成分,偶次諧波成分會(huì)給磁通門(mén)帶來(lái)系統(tǒng)噪聲。因此,采用正弦信號(hào)作為激磁信號(hào),具有較強(qiáng)的抗畸變能力,且不含偶次諧波成分,對(duì)提高磁通門(mén)靈敏度的激磁電路設(shè)計(jì)具有重要意義。

    目前,國(guó)外醫(yī)用鐵磁探測(cè)主流產(chǎn)品為英國(guó)Metrasens公司的ferroguard 系列MRI 檢測(cè)系統(tǒng),其性能優(yōu)異彌補(bǔ)了傳統(tǒng)金屬檢測(cè)儀器在MRI 檢測(cè)中的缺點(diǎn)和不足,但是其價(jià)格高昂,且技術(shù)保密性嚴(yán)格。而國(guó)內(nèi)廣泛采用的模擬式信號(hào)檢測(cè)處理方法易受電子器件和外部強(qiáng)磁環(huán)境影響,系統(tǒng)的穩(wěn)定性及檢測(cè)準(zhǔn)確性和靈敏度不高。本文在上述理論分析基礎(chǔ)上,設(shè)計(jì)了一種提高M(jìn)RI 檢測(cè)區(qū)域內(nèi)強(qiáng)磁場(chǎng)環(huán)境下鐵磁探測(cè)系統(tǒng)靈敏度和整體性能的方法:根據(jù)磁通門(mén)傳感器測(cè)磁原理,以FPGA 作為系統(tǒng)主控制器,改變傳統(tǒng)模擬信號(hào)的檢測(cè)方式,采用基于二次諧波原理的數(shù)字信號(hào)處理并嘗試增加濾波算法以提高系統(tǒng)鐵磁性物質(zhì)探測(cè)靈敏度,可將此處理方法推廣到其它弱磁探測(cè)領(lǐng)域,應(yīng)用前景廣泛。

    1 系統(tǒng)整體方案設(shè)計(jì)

    基于磁通門(mén)傳感器的鐵磁探測(cè)系統(tǒng)整體方案如圖1所示,主要包括磁通門(mén)傳感器模塊、激勵(lì)信號(hào)模塊、信號(hào)調(diào)理模塊、數(shù)據(jù)采集模塊、FPGA 內(nèi)部數(shù)字信號(hào)處理和算法模塊及報(bào)警模塊??删幊踢壿嬈骷﨔PGA 作為主控器驅(qū)動(dòng)外部DDS 模塊產(chǎn)生正弦激勵(lì)信號(hào)并傳給磁通門(mén)探頭1,在磁通門(mén)探頭鐵芯到達(dá)過(guò)飽和后,將外部磁場(chǎng)調(diào)制到激勵(lì)的交變磁場(chǎng)中,經(jīng)過(guò)探頭轉(zhuǎn)化為電壓信號(hào)。轉(zhuǎn)化完成后的電信號(hào)經(jīng)過(guò)信號(hào)調(diào)理后,再由ADC 采集到FPGA 主控中進(jìn)行數(shù)字信號(hào)處理:①首先對(duì)由探頭1 采集到的電壓信號(hào)做基于數(shù)字正交鎖相的二次諧波式檢波,探頭2 的信號(hào)處理方式與探頭1 完全相同;②將探頭1 處理好的信號(hào)作為輸入信號(hào),探頭2 的信號(hào)作為參考信號(hào),同時(shí)進(jìn)入濾波模塊進(jìn)行算法處理,最后將數(shù)據(jù)輸出作為評(píng)判標(biāo)準(zhǔn)驅(qū)動(dòng)語(yǔ)音燈光報(bào)警。

    Fig.1 Overall scheme design of ferromagnetic detection system圖1 鐵磁探測(cè)系統(tǒng)總體方案設(shè)計(jì)

    2 硬件系統(tǒng)方案設(shè)計(jì)

    系統(tǒng)硬件從信號(hào)傳感器端到信號(hào)處理端,主要包括:磁通門(mén)傳感器模塊、激勵(lì)信號(hào)模塊、信號(hào)調(diào)理模塊以及數(shù)據(jù)采集模塊。系統(tǒng)主控制模塊使用可編程邏輯器件FP?GA,型號(hào)為Xilinx_ZYNQ7000。

    2.1 磁通門(mén)傳感器模塊

    磁通門(mén)現(xiàn)象[8]實(shí)際上屬于變壓器效應(yīng)的衍生現(xiàn)象,因此磁通門(mén)傳感器與變壓器結(jié)構(gòu)類(lèi)似。如圖2 所示,取高磁導(dǎo)率、易飽和材料作為磁芯,在其上分別繞兩組線圈:一組作為激勵(lì)線圈,另一組作為感應(yīng)線圈。在交流信號(hào)的激勵(lì)下,使磁芯處于周期性飽和狀態(tài),當(dāng)磁芯處于非飽和磁場(chǎng)中,其磁導(dǎo)率變化緩慢,而當(dāng)磁芯達(dá)到飽和時(shí),其磁導(dǎo)率變化明顯,此時(shí)被測(cè)磁場(chǎng)被調(diào)制進(jìn)感應(yīng)電勢(shì)中,可以通過(guò)測(cè)量磁通門(mén)傳感器感應(yīng)電勢(shì)中能夠反映被測(cè)磁場(chǎng)的量來(lái)度量磁場(chǎng)。磁通門(mén)傳感器的工作過(guò)程中,磁芯的飽和點(diǎn)[9]貌似一道“門(mén)”,通過(guò)這道“門(mén)”,被測(cè)磁場(chǎng)被調(diào)制成偶次諧波的感應(yīng)電動(dòng)勢(shì),從而反映待測(cè)磁場(chǎng)的大小。

    Fig.2 Fluxgate sensor structure圖2 磁通門(mén)傳感器結(jié)構(gòu)

    磁通門(mén)傳感器測(cè)量磁場(chǎng)的結(jié)構(gòu)原理如圖3 所示。運(yùn)動(dòng)的鐵磁性物質(zhì)在交變激勵(lì)磁場(chǎng)下磁化產(chǎn)生磁性,即為待測(cè)磁場(chǎng)信號(hào)[10]。當(dāng)傳感器磁芯達(dá)到飽和時(shí),待測(cè)磁場(chǎng)信號(hào)被調(diào)制,從而產(chǎn)生反映其磁場(chǎng)強(qiáng)度大小的感應(yīng)動(dòng)勢(shì)(偶次諧波)[11]。信號(hào)處理部分完成對(duì)偶次諧波的檢測(cè)和提取。

    Fig.3 Structure principle of magnetic field measured by fluxgate圖3 磁通門(mén)測(cè)量磁場(chǎng)結(jié)構(gòu)原理

    2.2 激勵(lì)信號(hào)模塊

    信號(hào)激勵(lì)模塊采用主控芯片F(xiàn)PGA 驅(qū)動(dòng)外部直接頻率合成芯片AD9833[12]以產(chǎn)生20KHz 的正弦激勵(lì)信號(hào)。AD9833 為ADI 公司生產(chǎn)的一款直接數(shù)字頻率合成器,其價(jià)格功耗低且外圍電路簡(jiǎn)單、頻率和相位可編程。FPGA與AD9833 接口配置如圖所示,由于AD9833 為三線SPI接口,因此兩者之間通過(guò)SPI 協(xié)議進(jìn)行通信。在系統(tǒng)時(shí)鐘和使能信號(hào)正常工作后,F(xiàn)PGA 將配置好的16 位命令控制字通過(guò)spido 接口寫(xiě)入AD9833 相關(guān)的寄存器中,從而在VOUT 端輸入所需要的正弦激勵(lì)信號(hào)。由于AD9833 內(nèi)部集成了10bits 的數(shù)模轉(zhuǎn)化器(DAC),直接輸出為可用的模擬信號(hào),無(wú)需再外置DAC。

    2.3 信號(hào)調(diào)理模塊

    通常被測(cè)的物體體積較小,經(jīng)磁通門(mén)傳感器感應(yīng)出來(lái)的信號(hào)比較微弱,且檢測(cè)環(huán)境對(duì)于被測(cè)信號(hào)有一定干擾,因此要對(duì)被測(cè)信號(hào)進(jìn)行信號(hào)調(diào)理后再由AD 采集到FPGA中作進(jìn)一步處理。理論上講,信號(hào)調(diào)理模塊需要對(duì)信號(hào)進(jìn)行選頻放大[13],但是增加一級(jí)濾波會(huì)引入一些隨機(jī)噪聲,且在FPGA 中會(huì)對(duì)信號(hào)作正交鎖相放大處理,因此此模塊只對(duì)信號(hào)進(jìn)行放大處理。如圖5 所示,傳感器信號(hào)放大模塊放大部分,分別為一級(jí)放大和兩級(jí)放大,放大倍數(shù)都為10 倍,因此整體放大100 倍。

    Fig.4 Interface configuration between FPGA and AD9833圖4 FPGA 與AD9833 接口配置

    Fig.5 Signal amplification module圖5 放大模塊

    2.4 數(shù)據(jù)采集模塊

    數(shù)據(jù)采集模塊使用24bits 的AD7760 芯片,具有寬輸入帶寬、高信噪比(SNR)、集成度高以及適合高速采集等特點(diǎn)。AD7760[14]與FPGA 的接口配置如圖6 所示。系統(tǒng)上電復(fù)位后,首先等待至少2 個(gè)MCLK 周期;然后通過(guò)PARD[15:0]接口按照相應(yīng)的地址寫(xiě)入4 個(gè)寄存器的值,先寫(xiě)寄存器地址再寫(xiě)寄存器值以設(shè)置正確的時(shí)鐘頻率和數(shù)據(jù)輸出速率;最后進(jìn)行數(shù)據(jù)采集,先輸出低16 位有效位,再輸出高16 位有效位,并且都是低8 位在前,高8 位在后,其中低16 位有效位的低8 位為狀態(tài)位,因此要從AD 讀取一個(gè)轉(zhuǎn)換結(jié)果,需進(jìn)行兩次16 位讀操作。而且每次有新的轉(zhuǎn)換結(jié)果,DRDY 引腳則輸出一個(gè)低電平有效的脈沖,隨后嚴(yán)格按給定時(shí)序先將nRD/WR 置低,然后將nCS 置低,在其均變?yōu)榈碗娖讲痪?,?shù)據(jù)總線將變?yōu)橛行顟B(tài),然后輸出高16 位有效位,保持一定時(shí)間的低電平后先將nCS 電平置高,再將nRD/WR 電平置高,則完成一次16 讀操作,高電平需至少保持一個(gè)完整的ICLK 周期才能進(jìn)行第二次低16 位有效位的讀取。將讀出的數(shù)據(jù)寄存在FIFO 中作為緩存,以待后續(xù)數(shù)據(jù)處理。

    3 信號(hào)處理及算法方案設(shè)計(jì)

    從ADC 中采集到的數(shù)據(jù)先暫存在FIFO 中,然后對(duì)其進(jìn)行數(shù)字信號(hào)處理和濾波算法處理。此模塊對(duì)于感應(yīng)信號(hào)的提取和處理采用了基于數(shù)字正交鎖相的二次諧波法和自適應(yīng)濾波算法。

    Fig.6 Interface configuration between AD7760 and FPGA圖6 AD7760 與FPGA 的接口配置

    3.1 基于數(shù)字正交鎖相的二次諧波法

    經(jīng)由交變激勵(lì)磁場(chǎng)調(diào)制的待測(cè)磁場(chǎng)經(jīng)過(guò)磁通門(mén)傳感器感應(yīng)出電信號(hào),感應(yīng)信號(hào)的偶次諧波可反映磁場(chǎng)大小,但其中二次諧波分量最大,因而一般采用二次諧波提取法表征磁場(chǎng)大小。一般而言,提取特定頻率的信號(hào)有選頻放大和鎖相放大[15-16]兩種方法。鎖相放大的結(jié)構(gòu)和考慮因素較選頻放大復(fù)雜,但是實(shí)際效果要比選頻放大好,故本文采用數(shù)字式的正交鎖相放大。信號(hào)來(lái)自于ADC 采集的數(shù)字信號(hào),全部在FPGA 內(nèi)進(jìn)行處理,與模擬式的信號(hào)處理相比,受器件性能及外部溫度環(huán)境影響小,抗干擾能力強(qiáng)。

    基于數(shù)字正交鎖相的二次諧波法如圖7 所示。

    Fig.7 Second harmonic method based on digital orthogonal phase-locked圖7 基于數(shù)字正交鎖相的二次諧波法

    已經(jīng)由ADC 采集好的兩路探頭轉(zhuǎn)化信號(hào):輸入信號(hào)和參考信號(hào)。參考信號(hào)通過(guò)相移電路之后得到一組正交信號(hào),再將其與輸入信號(hào)相乘,得到兩路信號(hào)VI和VQ如式(1)和式(2)所示,將所測(cè)頻率的頻譜搬移到基帶上以方便進(jìn)行預(yù)濾波處理。

    然后,經(jīng)過(guò)一個(gè)窄帶的低通濾波器濾掉基帶以外的頻率分量,得到所需信號(hào);最后將濾波器輸出進(jìn)行平方、求和,再開(kāi)方,所求頻點(diǎn)信號(hào)幅度V如式(3)所示。

    3.2 自適應(yīng)濾波算法

    在微弱信號(hào)檢測(cè)和信號(hào)去噪領(lǐng)域[17],自適應(yīng)濾波器[18]在某種程度上可以最優(yōu)化方式消除包含在有用信號(hào)中的未知干擾。有用信號(hào)用作自適應(yīng)濾波器的期望響應(yīng),參考信號(hào)用作濾波器的輸入。其中,參考信號(hào)可以來(lái)自于定位的某一傳感器或一組傳感器[19],并以承載的新信號(hào)是微弱的或基本不可預(yù)測(cè)的方式,供給基本信號(hào)。

    自適應(yīng)濾波算法如圖8 所示,在本文中,基本輸入d(n)中含正交鎖相解調(diào)后反映磁場(chǎng)大小的信號(hào)和背景干擾噪聲,參考輸入u(n)只含背景干擾噪聲。首先計(jì)算線性濾波器對(duì)輸出及輸入信號(hào)的響應(yīng)y(n),然后通過(guò)比較產(chǎn)生誤差e(n),根據(jù)誤差自動(dòng)調(diào)整濾波器參數(shù)[20],在收斂后,參考輸入也即背景噪聲便從輸入信號(hào)中消除。本文調(diào)整濾波器參數(shù)的校正因子采用最小均方(LMS)更新算法[21]。

    Fig.8 Adaptive filtering algorithm module圖8 自適應(yīng)濾波算法模塊

    4 系統(tǒng)測(cè)試與結(jié)果分析

    實(shí)驗(yàn)檢測(cè)系統(tǒng)如圖9 所示,將搭建好的簡(jiǎn)易檢測(cè)系統(tǒng)置于MRI 室外進(jìn)行實(shí)驗(yàn)檢測(cè),將探頭垂直固定,兩個(gè)探頭經(jīng)由電路系統(tǒng)處理后將數(shù)據(jù)通過(guò)串口連接至PC 端,在上位機(jī)波形顯示窗口實(shí)時(shí)觀察信號(hào)波形,從而判斷是否有鐵磁性物質(zhì)通過(guò)。在本實(shí)驗(yàn)中,分別將手機(jī)、小刀(代替手術(shù)刀)、鑰匙和硬幣作為被測(cè)物檢測(cè)系統(tǒng)性能,如圖10 所示。

    Fig.9 Construction of detection system圖9 檢測(cè)系統(tǒng)搭建

    首先排除檢測(cè)距離對(duì)實(shí)驗(yàn)結(jié)果的影響,統(tǒng)一在離檢測(cè)系統(tǒng)60cm 處進(jìn)行標(biāo)記,實(shí)驗(yàn)人員手持被測(cè)物從標(biāo)記處走過(guò),得到不同被測(cè)物的信號(hào)波形,如圖11 所示。為方便觀察波形,將ADC 采集得到的信號(hào)數(shù)字碼進(jìn)行數(shù)字信號(hào)處理和自適應(yīng)濾波算法處理,用MATLAB 繪制出波形圖。如果沒(méi)有鐵磁物質(zhì)經(jīng)過(guò)傳感器,信號(hào)波形比較平穩(wěn),近乎為直線輸出;當(dāng)有鐵磁物質(zhì)靠近時(shí),平穩(wěn)波形會(huì)發(fā)生波動(dòng),從而引發(fā)報(bào)警系統(tǒng)示警。由圖11 可知,在被測(cè)物離探測(cè)系統(tǒng)60cm 處,手機(jī)對(duì)傳感器的擾動(dòng)最為明顯,硬幣最小。而沒(méi)有進(jìn)行數(shù)字信號(hào)和自適應(yīng)濾波算法處理的系統(tǒng)無(wú)法檢測(cè)到50cm 處的硬幣,因此證明此方案對(duì)于靈敏度檢測(cè)有一定提高。

    Fig.10 Measured objects圖10 被測(cè)物

    Fig.11 Signal waveforms of different measured objects圖11 不同被測(cè)物信號(hào)波形

    本文設(shè)計(jì)實(shí)驗(yàn)檢測(cè)系統(tǒng)靈敏度,在頻率為100~400Hz的激勵(lì)磁場(chǎng)下,在與傳感器平行位置放置數(shù)字特斯拉計(jì),用來(lái)測(cè)量檢測(cè)物通過(guò)時(shí)傳感器輸出的磁場(chǎng)大小,并用FUL?KE 5 位半數(shù)字電壓表測(cè)量傳感器的輸出電壓,靈敏度單位為uV/nT。記錄如表1 所示。

    Table 1 System sensitivity parameters表1 系統(tǒng)靈敏度參數(shù)

    在100~400Hz 的激勵(lì)范圍內(nèi),得到系統(tǒng)的平均靈敏度為1.66uV/nT,相比于原系統(tǒng)1.49uV/nT,靈敏度提高了10.3%。

    接下來(lái)檢驗(yàn)檢測(cè)距離對(duì)結(jié)果的影響,在離檢測(cè)系統(tǒng)40cm、50cm、60cm、70cm、80cm 以及90cm 處做標(biāo)記,按照上述實(shí)驗(yàn)分別做同樣檢測(cè)動(dòng)作,觀察不同被測(cè)物在不同檢測(cè)距離下的波形情況如表2 所示。由測(cè)試結(jié)果可知,此系統(tǒng)對(duì)于磁共振檢測(cè)中常見(jiàn)的鐵磁性物質(zhì)的探知識(shí)別距離可以達(dá)到1m 左右,原系統(tǒng)的檢測(cè)距離只有0.5m。

    Table 2 Waveforms of different measured objects at different detection distances表2 不同被測(cè)物在不同檢測(cè)距離下的波形情況

    5 結(jié)語(yǔ)

    本文設(shè)計(jì)了一種基于數(shù)字信號(hào)處理方法和濾波算法的可提高核磁共振室外鐵磁探測(cè)精度和準(zhǔn)確度的系統(tǒng)。在眾多模擬式信號(hào)處理方法原理基礎(chǔ)上,采用數(shù)字式信號(hào)處理方式并通過(guò)算法優(yōu)化極大提高系統(tǒng)精度和可靠性。通過(guò)實(shí)驗(yàn)結(jié)果可知,對(duì)于核磁共振檢測(cè)中常見(jiàn)的鐵磁性物質(zhì)的檢測(cè)靈敏度提高效果顯著,檢測(cè)距離可達(dá)1m??蓪⒃摲N方法可推廣至其它微弱信號(hào)檢測(cè)領(lǐng)域,例如食品金屬檢測(cè)、航空探潛、海洋監(jiān)測(cè)、地下和水下鐵磁物體探測(cè)等,應(yīng)用前景廣闊。

    猜你喜歡
    鐵磁磁通信號(hào)處理
    關(guān)于兩類(lèi)多分量海森堡鐵磁鏈模型的研究
    軸向磁通電勵(lì)磁雙凸極電機(jī)及容錯(cuò)運(yùn)行控制策略
    基于LabVIEW的微型磁通門(mén)磁強(qiáng)計(jì)測(cè)試系統(tǒng)搭建
    《信號(hào)處理》征稿簡(jiǎn)則
    《信號(hào)處理》第九屆編委會(huì)
    《信號(hào)處理》征稿簡(jiǎn)則
    《信號(hào)處理》第九屆編委會(huì)
    基于磁通門(mén)原理的零磁通交直流電流傳感器
    你好,鐵磁
    你好,鐵磁
    麻豆国产av国片精品| 成人国产综合亚洲| 国产成人系列免费观看| 嫩草影院精品99| 日韩高清综合在线| av黄色大香蕉| 亚洲色图av天堂| av中文乱码字幕在线| 亚洲av免费在线观看| 欧美成人一区二区免费高清观看 | 日韩欧美精品v在线| 色综合亚洲欧美另类图片| av国产免费在线观看| 色尼玛亚洲综合影院| 又爽又黄无遮挡网站| 国产成人啪精品午夜网站| 国产伦在线观看视频一区| 午夜视频精品福利| 国产欧美日韩精品亚洲av| a在线观看视频网站| 久久伊人香网站| 欧美精品啪啪一区二区三区| 村上凉子中文字幕在线| 91麻豆精品激情在线观看国产| 黑人欧美特级aaaaaa片| 99精品在免费线老司机午夜| 欧美绝顶高潮抽搐喷水| 老熟妇乱子伦视频在线观看| 国产黄片美女视频| 国产亚洲精品久久久久久毛片| 国产成+人综合+亚洲专区| 久久天躁狠狠躁夜夜2o2o| 久久九九热精品免费| 99久久国产精品久久久| 别揉我奶头~嗯~啊~动态视频| 亚洲人与动物交配视频| 欧美色欧美亚洲另类二区| 午夜激情欧美在线| 非洲黑人性xxxx精品又粗又长| 嫩草影院入口| 国产精华一区二区三区| 日本一二三区视频观看| 欧美成人一区二区免费高清观看 | 亚洲国产日韩欧美精品在线观看 | 亚洲色图av天堂| 亚洲中文字幕日韩| 日韩免费av在线播放| 美女免费视频网站| 欧美zozozo另类| 免费看日本二区| 在线观看舔阴道视频| 一级毛片女人18水好多| 草草在线视频免费看| 国产高清激情床上av| 男女做爰动态图高潮gif福利片| 久久久成人免费电影| 欧美黑人欧美精品刺激| 露出奶头的视频| 岛国在线免费视频观看| 亚洲欧美日韩高清专用| 国产精品亚洲av一区麻豆| 国产精品 欧美亚洲| 国产激情偷乱视频一区二区| 99久久99久久久精品蜜桃| 老司机在亚洲福利影院| 国产一级毛片七仙女欲春2| 久久久久亚洲av毛片大全| 一二三四社区在线视频社区8| 麻豆国产av国片精品| 啦啦啦韩国在线观看视频| 日韩欧美 国产精品| 99在线人妻在线中文字幕| 麻豆国产av国片精品| 美女被艹到高潮喷水动态| 成人av一区二区三区在线看| 欧美性猛交╳xxx乱大交人| aaaaa片日本免费| 国产午夜精品论理片| 少妇人妻一区二区三区视频| 成人无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 一级毛片高清免费大全| 亚洲精品久久国产高清桃花| 久久久久久久久久黄片| 麻豆国产97在线/欧美| 久久精品91无色码中文字幕| 麻豆成人午夜福利视频| 中亚洲国语对白在线视频| 亚洲欧美激情综合另类| 国内精品久久久久精免费| 在线观看66精品国产| 好男人电影高清在线观看| 亚洲国产看品久久| 国产一区在线观看成人免费| 亚洲精品色激情综合| 日日干狠狠操夜夜爽| 欧美日韩一级在线毛片| 欧美中文综合在线视频| 波多野结衣巨乳人妻| 国产视频一区二区在线看| 国产一区在线观看成人免费| 美女免费视频网站| 亚洲国产欧洲综合997久久,| 91在线精品国自产拍蜜月 | 丝袜人妻中文字幕| 久久久色成人| 国产激情久久老熟女| 一本精品99久久精品77| www日本在线高清视频| 97碰自拍视频| 丰满的人妻完整版| 露出奶头的视频| 久久久久久国产a免费观看| 又大又爽又粗| 国产淫片久久久久久久久 | 亚洲 国产 在线| 亚洲片人在线观看| 他把我摸到了高潮在线观看| 成在线人永久免费视频| 99国产综合亚洲精品| 嫁个100分男人电影在线观看| 老司机午夜十八禁免费视频| 欧美黄色淫秽网站| 高潮久久久久久久久久久不卡| 亚洲最大成人中文| 精品电影一区二区在线| av欧美777| 女生性感内裤真人,穿戴方法视频| 中国美女看黄片| 国产视频内射| 国产成人啪精品午夜网站| 怎么达到女性高潮| 香蕉久久夜色| 亚洲 欧美 日韩 在线 免费| 亚洲乱码一区二区免费版| 亚洲18禁久久av| 国产伦人伦偷精品视频| 无人区码免费观看不卡| 日韩欧美国产在线观看| www日本在线高清视频| 人人妻,人人澡人人爽秒播| 国产成年人精品一区二区| 国产成人精品无人区| 在线观看免费午夜福利视频| 日本撒尿小便嘘嘘汇集6| 亚洲在线自拍视频| 亚洲熟女毛片儿| 久久婷婷人人爽人人干人人爱| 天堂动漫精品| 久久久久九九精品影院| 成人三级做爰电影| 一本综合久久免费| 国产三级在线视频| 国产精品日韩av在线免费观看| 欧美成人一区二区免费高清观看 | 五月伊人婷婷丁香| 淫秽高清视频在线观看| 中文字幕精品亚洲无线码一区| 午夜久久久久精精品| 小蜜桃在线观看免费完整版高清| 亚洲一区二区三区色噜噜| 日韩三级视频一区二区三区| 国产精品 国内视频| 国产v大片淫在线免费观看| 一级毛片高清免费大全| 动漫黄色视频在线观看| 丰满的人妻完整版| 亚洲性夜色夜夜综合| 精品一区二区三区视频在线 | 午夜福利成人在线免费观看| 91麻豆av在线| 成人高潮视频无遮挡免费网站| 91av网一区二区| 免费看美女性在线毛片视频| 午夜福利视频1000在线观看| 免费看日本二区| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 国产91精品成人一区二区三区| 免费人成视频x8x8入口观看| av片东京热男人的天堂| 99热这里只有是精品50| 又紧又爽又黄一区二区| 国产又色又爽无遮挡免费看| 男女之事视频高清在线观看| 久久久色成人| 日韩大尺度精品在线看网址| 99视频精品全部免费 在线 | 狂野欧美白嫩少妇大欣赏| 国产精品 欧美亚洲| 成人高潮视频无遮挡免费网站| 色在线成人网| 舔av片在线| 后天国语完整版免费观看| 免费观看精品视频网站| 全区人妻精品视频| 一级a爱片免费观看的视频| 国产精品爽爽va在线观看网站| 国产成人aa在线观看| a在线观看视频网站| 久久午夜综合久久蜜桃| 一个人看视频在线观看www免费 | 国产乱人视频| 国产欧美日韩一区二区三| 一个人免费在线观看电影 | 国产精品影院久久| 极品教师在线免费播放| 免费在线观看亚洲国产| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 丰满人妻一区二区三区视频av | 久久99热这里只有精品18| 一个人免费在线观看的高清视频| 国产淫片久久久久久久久 | 欧美日韩亚洲国产一区二区在线观看| 18禁观看日本| 十八禁人妻一区二区| 久久久久免费精品人妻一区二区| 国产成年人精品一区二区| 成人鲁丝片一二三区免费| 亚洲精品在线美女| 中亚洲国语对白在线视频| 中文字幕av在线有码专区| 丰满的人妻完整版| 老司机午夜福利在线观看视频| 全区人妻精品视频| 母亲3免费完整高清在线观看| 真人做人爱边吃奶动态| 热99在线观看视频| 欧美黄色淫秽网站| 日韩欧美国产一区二区入口| 最近在线观看免费完整版| 欧美精品啪啪一区二区三区| 黄色视频,在线免费观看| 亚洲天堂国产精品一区在线| 国产精品久久久av美女十八| 99在线视频只有这里精品首页| 曰老女人黄片| 男女之事视频高清在线观看| 9191精品国产免费久久| 欧美黄色片欧美黄色片| 亚洲第一电影网av| av国产免费在线观看| 久99久视频精品免费| 精品福利观看| 亚洲av第一区精品v没综合| 成人一区二区视频在线观看| 国产免费av片在线观看野外av| 国内精品久久久久久久电影| 成人av一区二区三区在线看| 青草久久国产| 欧美激情久久久久久爽电影| 亚洲九九香蕉| 波多野结衣巨乳人妻| 国产91精品成人一区二区三区| 久久久国产成人精品二区| 999精品在线视频| 免费人成视频x8x8入口观看| 一个人看视频在线观看www免费 | 99久久精品国产亚洲精品| 亚洲电影在线观看av| 欧美激情在线99| 日日摸夜夜添夜夜添小说| av欧美777| 精品一区二区三区四区五区乱码| 全区人妻精品视频| 成人无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 国产欧美日韩精品亚洲av| 久久伊人香网站| 亚洲成人精品中文字幕电影| 99热6这里只有精品| 熟女电影av网| 啪啪无遮挡十八禁网站| 日本熟妇午夜| 国内毛片毛片毛片毛片毛片| 波多野结衣巨乳人妻| 久久久久久久久久黄片| 美女扒开内裤让男人捅视频| 一级a爱片免费观看的视频| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 亚洲五月婷婷丁香| 一级黄色大片毛片| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区黑人| 一区二区三区国产精品乱码| 国产午夜精品久久久久久| 国产高清videossex| 成人三级黄色视频| 国产单亲对白刺激| 亚洲国产精品sss在线观看| 狂野欧美白嫩少妇大欣赏| 精品国产亚洲在线| 久久性视频一级片| 老熟妇乱子伦视频在线观看| 久久久久久久精品吃奶| 又黄又爽又免费观看的视频| 国产高清有码在线观看视频| 免费在线观看视频国产中文字幕亚洲| 一个人免费在线观看的高清视频| 啦啦啦韩国在线观看视频| 宅男免费午夜| 亚洲av电影不卡..在线观看| 巨乳人妻的诱惑在线观看| 国产成人一区二区三区免费视频网站| 成人av一区二区三区在线看| 久久久国产成人免费| 欧美日韩综合久久久久久 | 欧美乱码精品一区二区三区| 黄色视频,在线免费观看| 国产1区2区3区精品| 国产精品一区二区三区四区久久| 最近最新中文字幕大全电影3| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 欧美日本亚洲视频在线播放| 日本黄色片子视频| 国产男靠女视频免费网站| 国产真实乱freesex| 国产视频一区二区在线看| 夜夜爽天天搞| 国产精品99久久99久久久不卡| 搡老妇女老女人老熟妇| 久久精品国产综合久久久| 国产午夜精品久久久久久| 两个人看的免费小视频| 国产精品一及| 观看免费一级毛片| 99热只有精品国产| 黑人欧美特级aaaaaa片| 日本精品一区二区三区蜜桃| 黄色丝袜av网址大全| 在线观看一区二区三区| 美女扒开内裤让男人捅视频| 国产精品av视频在线免费观看| 亚洲欧美激情综合另类| 男人舔女人的私密视频| 女警被强在线播放| 神马国产精品三级电影在线观看| 成人av在线播放网站| 成年女人看的毛片在线观看| 国产精华一区二区三区| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 熟女电影av网| 久久午夜亚洲精品久久| 成人特级黄色片久久久久久久| a级毛片在线看网站| 国产毛片a区久久久久| 久久中文字幕人妻熟女| 不卡av一区二区三区| 免费看美女性在线毛片视频| 超碰成人久久| 91久久精品国产一区二区成人 | 国产极品精品免费视频能看的| 久久久久九九精品影院| 不卡一级毛片| 精品一区二区三区视频在线 | xxx96com| 久久久成人免费电影| 国内久久婷婷六月综合欲色啪| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 久久国产乱子伦精品免费另类| 亚洲成a人片在线一区二区| aaaaa片日本免费| 午夜免费成人在线视频| 亚洲精品456在线播放app | 精品午夜福利视频在线观看一区| 久久精品91蜜桃| 亚洲欧美日韩卡通动漫| 国产精品久久电影中文字幕| 日本黄色视频三级网站网址| 亚洲无线在线观看| 天堂影院成人在线观看| 日本与韩国留学比较| 香蕉国产在线看| 亚洲午夜精品一区,二区,三区| 在线免费观看的www视频| 久久这里只有精品19| 黄色成人免费大全| 露出奶头的视频| 一卡2卡三卡四卡精品乱码亚洲| 黄色成人免费大全| 中文字幕熟女人妻在线| 亚洲第一欧美日韩一区二区三区| 日韩中文字幕欧美一区二区| 亚洲精品色激情综合| 久久精品国产99精品国产亚洲性色| 亚洲精品在线美女| 麻豆成人午夜福利视频| 男人舔女人的私密视频| 欧美黄色淫秽网站| 日本 av在线| 日本三级黄在线观看| 久久这里只有精品中国| 亚洲欧洲精品一区二区精品久久久| 午夜免费观看网址| 久久久精品大字幕| 欧美性猛交黑人性爽| 黑人操中国人逼视频| 国产精品精品国产色婷婷| 无限看片的www在线观看| а√天堂www在线а√下载| 在线十欧美十亚洲十日本专区| 日日干狠狠操夜夜爽| h日本视频在线播放| 又大又爽又粗| 一进一出抽搐gif免费好疼| 日韩有码中文字幕| 亚洲真实伦在线观看| 桃色一区二区三区在线观看| 国产成人系列免费观看| 色尼玛亚洲综合影院| 美女黄网站色视频| 脱女人内裤的视频| 俄罗斯特黄特色一大片| 色老头精品视频在线观看| 免费搜索国产男女视频| 亚洲精品在线美女| 在线免费观看不下载黄p国产 | 亚洲国产中文字幕在线视频| 欧美黄色片欧美黄色片| 亚洲av中文字字幕乱码综合| 制服丝袜大香蕉在线| 激情在线观看视频在线高清| 又大又爽又粗| 国内少妇人妻偷人精品xxx网站 | 国产精品免费一区二区三区在线| 一级毛片女人18水好多| 黄色视频,在线免费观看| 免费看光身美女| 亚洲国产欧美一区二区综合| 不卡一级毛片| 97超级碰碰碰精品色视频在线观看| 国产精品久久久av美女十八| 黄色丝袜av网址大全| 淫妇啪啪啪对白视频| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 国产成人av教育| 亚洲熟女毛片儿| 99riav亚洲国产免费| 亚洲熟女毛片儿| 欧美成狂野欧美在线观看| 精品国产超薄肉色丝袜足j| 日本一二三区视频观看| 中文亚洲av片在线观看爽| 亚洲精品久久国产高清桃花| 国产精品98久久久久久宅男小说| 国产成人精品无人区| 国产毛片a区久久久久| 51午夜福利影视在线观看| 国产午夜精品论理片| 99在线人妻在线中文字幕| 级片在线观看| 国内揄拍国产精品人妻在线| 亚洲国产精品久久男人天堂| 日韩欧美精品v在线| 最新在线观看一区二区三区| 午夜影院日韩av| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 亚洲国产欧美一区二区综合| 国产毛片a区久久久久| 亚洲七黄色美女视频| 9191精品国产免费久久| 亚洲中文av在线| 可以在线观看毛片的网站| 99精品欧美一区二区三区四区| 日韩欧美国产在线观看| 国产精品久久久人人做人人爽| 欧洲精品卡2卡3卡4卡5卡区| 久久国产乱子伦精品免费另类| 国产日本99.免费观看| 久久久久久九九精品二区国产| 麻豆国产av国片精品| 国产乱人伦免费视频| 麻豆国产97在线/欧美| 亚洲av成人精品一区久久| www日本黄色视频网| 国产日本99.免费观看| 亚洲自拍偷在线| 在线看三级毛片| 99国产极品粉嫩在线观看| 欧美色视频一区免费| 亚洲 国产 在线| 久久这里只有精品19| 婷婷精品国产亚洲av在线| 欧美在线一区亚洲| 国产精品国产高清国产av| 日本a在线网址| 午夜精品一区二区三区免费看| 日本精品一区二区三区蜜桃| 国产高清三级在线| 久久天堂一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 九九久久精品国产亚洲av麻豆 | 国产熟女xx| 最近视频中文字幕2019在线8| 啦啦啦韩国在线观看视频| 琪琪午夜伦伦电影理论片6080| 天堂av国产一区二区熟女人妻| 一个人免费在线观看的高清视频| 黑人操中国人逼视频| 性色av乱码一区二区三区2| 亚洲专区字幕在线| 久久久久久大精品| 搡老熟女国产l中国老女人| 免费大片18禁| 日韩国内少妇激情av| 91在线精品国自产拍蜜月 | 日韩免费av在线播放| 欧美乱妇无乱码| 国产极品精品免费视频能看的| 两性夫妻黄色片| 熟女电影av网| 免费搜索国产男女视频| 成人三级黄色视频| aaaaa片日本免费| 国产三级在线视频| 国产亚洲欧美98| 国产野战对白在线观看| 国产精品亚洲av一区麻豆| 亚洲真实伦在线观看| 97超级碰碰碰精品色视频在线观看| 人人妻,人人澡人人爽秒播| 免费在线观看亚洲国产| 国产av在哪里看| 亚洲一区高清亚洲精品| 美女免费视频网站| 成人亚洲精品av一区二区| av欧美777| 丰满人妻熟妇乱又伦精品不卡| 99久久精品热视频| 欧美国产日韩亚洲一区| 宅男免费午夜| 亚洲人成网站高清观看| 97超视频在线观看视频| 天堂网av新在线| 九色成人免费人妻av| netflix在线观看网站| 欧美成人一区二区免费高清观看 | 欧美成人免费av一区二区三区| 久久久久久久久久黄片| 成人永久免费在线观看视频| 久久性视频一级片| 亚洲自拍偷在线| 国产综合懂色| 黄频高清免费视频| 动漫黄色视频在线观看| 婷婷丁香在线五月| 国产午夜精品久久久久久| 网址你懂的国产日韩在线| 99精品久久久久人妻精品| 日本 av在线| 国产av麻豆久久久久久久| 亚洲国产欧美一区二区综合| 日韩免费av在线播放| 日韩欧美在线二视频| 19禁男女啪啪无遮挡网站| 国内精品一区二区在线观看| 成人性生交大片免费视频hd| 国产主播在线观看一区二区| 午夜日韩欧美国产| 亚洲国产欧美一区二区综合| 日本五十路高清| 黄色丝袜av网址大全| 久久中文字幕一级| 亚洲 欧美 日韩 在线 免费| 精品午夜福利视频在线观看一区| 19禁男女啪啪无遮挡网站| 黄色成人免费大全| 99热6这里只有精品| 在线免费观看的www视频| 国产伦在线观看视频一区| 亚洲成av人片免费观看| 国产真人三级小视频在线观看| 成人欧美大片| 国模一区二区三区四区视频 | 中文在线观看免费www的网站| 国产精品av久久久久免费| 亚洲电影在线观看av| 婷婷丁香在线五月| 久久久国产精品麻豆| 国产精品98久久久久久宅男小说| 桃色一区二区三区在线观看| 亚洲,欧美精品.| 日韩精品中文字幕看吧| 麻豆成人午夜福利视频| 久久久成人免费电影| 黑人操中国人逼视频| 一本一本综合久久| 精华霜和精华液先用哪个| 一级毛片高清免费大全| 国产精品久久视频播放| 2021天堂中文幕一二区在线观| 国产aⅴ精品一区二区三区波| 亚洲av第一区精品v没综合| av在线蜜桃| 亚洲第一电影网av| 亚洲性夜色夜夜综合| 我的老师免费观看完整版| 久久伊人香网站| 叶爱在线成人免费视频播放| 1000部很黄的大片| 国产高清videossex| 俄罗斯特黄特色一大片| 成人av一区二区三区在线看| 欧美另类亚洲清纯唯美| 俄罗斯特黄特色一大片| 成年女人看的毛片在线观看| 蜜桃久久精品国产亚洲av|