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

    基于多迭代變分模態(tài)分解與復(fù)合多尺度散布熵的生物組織變性識(shí)別方法*

    2020-01-02 06:33:44胡偉鵬趙新民錢盛友
    傳感技術(shù)學(xué)報(bào) 2019年12期
    關(guān)鍵詞:背散射變性信噪比

    胡偉鵬,鄒 孝,劉 備,趙新民,錢盛友

    (湖南師范大學(xué)物理與電子科學(xué)學(xué)院,長(zhǎng)沙 410081)

    高強(qiáng)度聚焦超聲治療就是將超聲波的能量聚焦于靶區(qū),使靶區(qū)內(nèi)產(chǎn)生高溫,利用蛋白質(zhì)在一定的溫度下會(huì)產(chǎn)生變性使組織細(xì)胞壞死的原理,達(dá)到治療或組織切除目的[1]。在治療過(guò)程中,隨著超聲波聚焦于靶區(qū),靶區(qū)的細(xì)胞逐漸失去活性[2],位于超聲波傳播路徑上的正常組織也會(huì)吸收超聲能量,且越接近靶區(qū)的組織吸收的能量越多,因此使用超聲治療中需要實(shí)時(shí)監(jiān)測(cè)治療區(qū)域情況。迄今為止,超聲領(lǐng)域的研究人員已經(jīng)從超聲衰減系數(shù),聲速和熵等參數(shù)對(duì)采集的超聲信號(hào)進(jìn)行研究,希望找到能準(zhǔn)確反映生物組織特性的參數(shù)[3-6]。但采集的超聲背散射回波信號(hào)含有大量的噪聲,對(duì)判斷生物組織是否變性造成了很大的困難。減少超聲背散射信號(hào)中的噪聲,對(duì)準(zhǔn)確判斷生物組織是否發(fā)生變性非常重要。

    經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)與變分模態(tài)分解(Variational Mode Decomposition,VMD)憑借其強(qiáng)大的分析和去噪能力廣泛應(yīng)用于各個(gè)領(lǐng)域[7-11]。但是EMD得到的本征模函數(shù)(Intrinsic Mode Function,IMF)存在模態(tài)混疊,微弱信號(hào)難以提取,容易產(chǎn)生虛假信號(hào)[12-13],對(duì)于頻率相近的分量難以正確分離。與EMD相比VMD的魯棒性較好[14],在分析頻率相近的分量時(shí)更加容易分離,但是VMD分解的模態(tài)數(shù)是一個(gè)定值,在分解過(guò)程中部分噪聲信號(hào)會(huì)與有用信號(hào)分在同一個(gè)IMF中。本文基于VMD提出了多迭代變分模態(tài)分解(Multi-Iterative Variational Mode Decomposition,MIVMD),MIVMD在計(jì)算過(guò)程中,計(jì)算了各個(gè)IMF的能量比,將低能量比的IMF濾除,并根據(jù)重構(gòu)信號(hào)中有用頻率的能量比進(jìn)行多次分解重構(gòu)。

    由于超聲波在生物組織內(nèi)的非線性傳播特性,實(shí)驗(yàn)采集的超聲背散射信號(hào)是一種非線性信號(hào)[15]。信號(hào)的復(fù)雜程度會(huì)根據(jù)生物組織變性程度的改變而改變。目前為止,我們通常通過(guò)計(jì)算信號(hào)的熵值來(lái)評(píng)估信號(hào)的復(fù)雜程度[16-18]。散布熵(Dispersion Entropy,DE)是一種新的評(píng)估時(shí)間序列動(dòng)態(tài)特征的參數(shù),DE不會(huì)產(chǎn)生未定義的熵值,且時(shí)間序列的微小改變不會(huì)對(duì)DE造成很大影響[19]。復(fù)合多尺度散布熵(Composite Multi-scale Dispersion Entropy,CMDE)是DE的改進(jìn),進(jìn)行多尺度處理,能從更多尺度下對(duì)時(shí)間序列的復(fù)雜性進(jìn)行分析[20]。復(fù)合化能克服多尺度散布熵的缺點(diǎn),能結(jié)合多個(gè)粗?;蛄械男畔?降低熵值的標(biāo)準(zhǔn)偏差使熵值更可靠。針對(duì)生物組織在一定溫度下逐漸產(chǎn)生變性的特性,GK模糊聚類可以根據(jù)輸入特征產(chǎn)生一個(gè)聚類中心,根據(jù)特征對(duì)聚類中心的隸屬度判斷輸入特征屬于哪一類,因此可以判斷生物組織是否變性。

    通過(guò)上述分析,本文提出了MIVMD-CMDE-GK模糊聚類的生物組織變性識(shí)別方法。將采集的超聲背散射信號(hào)進(jìn)行MIVMD重構(gòu),計(jì)算重構(gòu)信號(hào)的CMDE,并將得到的所有CMDE作為特征參數(shù)通過(guò)GK模糊聚類得到兩個(gè)聚類中心,計(jì)算熵值與聚類中心的貼近度,根據(jù)貼近度的大小判斷生物組織是否變性。實(shí)驗(yàn)結(jié)果表明,本文所提的MIVMD-CMDE-GK方法有效提高了超聲背散射信號(hào)對(duì)生物組織變性識(shí)別的能力,為超聲治療診斷提供了參考。

    1 原理與方法

    1.1 模態(tài)分解

    1.1.1 變分模態(tài)分解

    VMD算法中,原始信號(hào)分解得到的IMF為調(diào)頻-調(diào)幅信號(hào),表達(dá)式為:

    uk(t)=Ak(t)cos[φk(t)]

    (1)

    (2)

    式中:uk(t)為VMD分解得到的第k個(gè)IMF,Ak(t)為瞬時(shí)幅值,φk(t)為瞬時(shí)相位,ω′k(t)為瞬時(shí)頻率。

    VMD計(jì)算過(guò)程中,不斷更新IMF的中心頻率和帶寬,在各個(gè)IMF的帶寬和最小時(shí)停止計(jì)算,且各IMF之和為原始輸入信號(hào)。約束變分問題函數(shù)如下:

    (3)

    式中:{ωk}={ω1,ω2,…,ωk}為IMF的中心頻率。通過(guò)使用增廣拉格朗日函數(shù)對(duì)約束變分問題進(jìn)行求解,如式(4)所示:

    (4)

    式中:α為二次懲罰參數(shù),λ(t)為拉格朗日懲罰算子。使用交替方向乘數(shù)法求解式(4)的鞍點(diǎn),同時(shí)更新uk,ωk,λ三種參數(shù)。將原始信號(hào)x(t),通過(guò)VMD分解成k個(gè)IMF的方法如下[21]:

    ②更新uk,ωk,λ:

    (5)

    (6)

    (7)

    ③設(shè)置精度ω,滿足當(dāng)式(8)時(shí),輸出結(jié)果uk。否則重復(fù)②步。

    (8)

    1.1.2 多迭代變分模態(tài)分解

    由于噪聲頻率分布范圍廣,而VMD分解重構(gòu)得到的模態(tài)數(shù)有限,因此得到每個(gè)IMF中都含有噪聲。為了更好的將信號(hào)中的噪聲濾除,提出了MIVMD,通過(guò)多次選擇有用信號(hào)能量占比高的分量重構(gòu),達(dá)到濾除噪聲的目的。MIVMD計(jì)算方法入下:

    ①對(duì)輸入信號(hào)進(jìn)行VMD分解,IMF數(shù)為k。

    ②計(jì)算各IMF中的主頻能量占此IMF能量的比,能量比越大說(shuō)明噪聲信號(hào)越少。

    ③將能量比小的IMF分量濾除,使用剩余分量重構(gòu)信號(hào)。

    ④計(jì)算重構(gòu)信號(hào)的頻譜能量,計(jì)算最大的前k-1個(gè)頻譜能量占重構(gòu)信號(hào)總能量的比。

    ⑤設(shè)定占比閾值,當(dāng)能量占比小于閾值時(shí),將重構(gòu)信號(hào)視為輸入信號(hào)并重復(fù)步驟(1)~步驟(4)。當(dāng)能量占比大于閾值時(shí),停止迭代,輸出重構(gòu)信號(hào)。

    1.2 復(fù)合多尺度散布熵

    對(duì)長(zhǎng)度為N采集信號(hào)x(t),計(jì)算信號(hào)的DE如下[22]:

    ①通過(guò)正太分布函數(shù),將采集信號(hào)x(t)映射成y(t)

    (9)

    式中:σ為標(biāo)準(zhǔn)差,mean為采集信號(hào)的平均值。

    ②將y(t)映通過(guò)式(10)映射成到1至c的集合z(t),c為類別數(shù)。

    zc(t)=round[cy(t)+0.5]

    (10)

    (11)

    ④每種分散模式的概率為:

    (12)

    ⑤計(jì)算DE:

    (13)

    CMDE是對(duì)DE的改進(jìn),先對(duì)原始輸入信號(hào)x(t)進(jìn)行粗?;叨纫蜃訛棣拥奶幚?得到第K個(gè)粗粒化序列:

    (14)

    然后對(duì)得到的τ個(gè)粗?;蛄猩⒉检厍缶?得到該粗粒化尺度因子下的CMDE為:

    (15)

    1.3 GK模糊聚類

    GK模糊聚類是一種距離協(xié)方差矩陣能自適應(yīng)動(dòng)態(tài)度量的模糊聚類算法,假設(shè)輸入的特征為B={b1,b2,…,bN},GK模糊聚類的目標(biāo)函數(shù)為[23]:

    (16)

    可以得到e個(gè)聚類中心Oi(i=1,2,…,e);f為模糊指數(shù),模糊指數(shù)會(huì)影響聚類效果,f太大會(huì)導(dǎo)致各類之間相互重疊;U為隸屬度矩陣U=[μij]e×N;μij表示第j個(gè)元素隸屬于第i類的隸屬度且滿足:

    (17)

    Dij表示第j個(gè)元素與聚類中心Oi的距離泛數(shù):

    (18)

    (19)

    常用拉格朗日乘法對(duì)GK模糊聚類的目標(biāo)函數(shù)進(jìn)行優(yōu)化,得到最小值點(diǎn)(U,V),其必要條件為:

    (20)

    (21)

    1.4 評(píng)價(jià)參數(shù)

    為了便于比較不同方法的去噪效果,使用信噪比和均方根誤差來(lái)進(jìn)行評(píng)價(jià)。

    信噪比:

    圖1 仿真信號(hào)與重構(gòu)信號(hào)的頻譜對(duì)比

    (22)

    Ps,Pn分別代表有用信號(hào)和噪聲的有效功率。信噪比越大說(shuō)明去噪效果越好。

    均方根誤差:

    (23)

    signal′,signal,l分別為重構(gòu)信號(hào)、未加噪信號(hào)、信號(hào)長(zhǎng)度。均方根誤差越小說(shuō)明去噪效果越好。

    本文采用劃分系數(shù)PC、劃分熵系數(shù)PE和Xie-beni指數(shù)評(píng)估聚類效果[24]。PC,PE均與隸屬度矩陣有關(guān),當(dāng)PC越接近1與PE越大接近0時(shí),聚類的效果越好。XB是用來(lái)評(píng)估類間距離的一種指數(shù),值越小,聚類類間距離越大,聚類效果越好。

    劃分系數(shù)PC:

    (24)

    劃分熵系數(shù)PE:

    (25)

    Xie-beni指數(shù)表達(dá)式為:

    (26)

    式中:θ為類的平均方差,Jmin是類間最短模糊距離。

    2 實(shí)驗(yàn)分析

    2.1 信號(hào)去噪

    使用VMD、MIVMD、EMD,分別對(duì)加噪信號(hào)進(jìn)行重構(gòu)分解。未加噪仿真信號(hào)為x(t)=3sin(40πt)+2sin(200πt),加噪仿真信號(hào)為Q(t)=x(t)+η,η為高斯白噪聲。對(duì)加噪仿真信號(hào)分別進(jìn)行VMD分解重構(gòu)、MIVMD分解重構(gòu)并與EMD分解重構(gòu)進(jìn)行比較。圖1(a)為信噪比為2.31 dB時(shí)加噪仿真信號(hào)頻譜圖,圖1(b)~1(d)分別為VMD,MIVMD,EMD三種方法重構(gòu)信號(hào)的頻譜圖。從圖1可以發(fā)現(xiàn),EMD和VMD均能抑制中高頻噪聲,但是對(duì)兩個(gè)有用頻率之間的噪聲抑制效果較差。而本文提出的MIVMD不僅能抑制中高頻噪聲,還能有效地抑制兩個(gè)有用頻率之間的噪聲。

    為了驗(yàn)證MIVMD對(duì)低信噪比信號(hào)有更好的濾波能力,使用EMD、VMD和MIVMD對(duì)信噪比不同的仿真信號(hào)進(jìn)行重構(gòu),并計(jì)算重構(gòu)信號(hào)的信噪比和均方根誤差如圖2所示。由圖2可知,MIVMD重構(gòu)信號(hào)的信噪比高于EMD與VMD重構(gòu)信號(hào)的信噪比,MIVMD重構(gòu)信號(hào)的均方根誤差低于EMD與VMD重構(gòu)信號(hào)的均方根誤差。同時(shí)可以發(fā)現(xiàn),VMD與EMD重構(gòu)信號(hào)的信噪比隨輸入信號(hào)信噪比增加而增加,而MIVMD在計(jì)算過(guò)程中,對(duì)各個(gè)IMF分量的能量進(jìn)行了分析,因此MIVMD重構(gòu)信號(hào)的信噪比在某一定值附近浮動(dòng)。

    圖3(a)、3(b)分別為通過(guò)光纖水聽器采集的正常與變性組織的超聲背散射信號(hào),使用MIVMD對(duì)采集的超聲背散射信號(hào)進(jìn)行分解重構(gòu),模態(tài)數(shù)設(shè)定為5,結(jié)果分別如圖3(c)、3(d)所示??梢杂^察到通過(guò)MIVMD重構(gòu)獲得的回波信號(hào)有明顯的脈沖波形,且波形震蕩衰減。

    圖2 不同輸入信噪比下不同方法的去噪效果比較

    圖3 超聲背散射回波信號(hào)

    2.2 變性識(shí)別

    為了驗(yàn)證本文提出的生物組織變性識(shí)別方法的有效性,本文將該方法的判定結(jié)果與實(shí)際切片判斷結(jié)果進(jìn)行比較。實(shí)驗(yàn)使用HIFU對(duì)新鮮離體豬肉組織輻照來(lái)改變生物組織特性,使用熱敏電阻測(cè)量聲焦域處的生物組織溫度,最大溫度在90 ℃以內(nèi);考慮到HIFU環(huán)境中強(qiáng)超聲的影響,使用光纖水聽器(FOPH2000;RP acoustics)獲取監(jiān)控超聲的回波信號(hào),并經(jīng)數(shù)字示波器(Model MDO3032;Tektronix)轉(zhuǎn)化為數(shù)字信號(hào)后進(jìn)行保存。采集了16組樣本共288例實(shí)驗(yàn)數(shù)據(jù)。為找到CMDE的最佳尺度,研究了3到20尺度下,第一組樣本的超聲背散射信號(hào)。為驗(yàn)證MIVMD處理能使提取的特征參數(shù)更加有效,使用CMDE對(duì)經(jīng)VMD和MIVMD處理的重構(gòu)信號(hào)進(jìn)行分析。在使用CMDE進(jìn)行分析時(shí),為了避免將所有數(shù)據(jù)視為一種散布模式,需要將類別數(shù)設(shè)定大于1,本文將變性與未變性組織的信號(hào)視為兩種類別,因此設(shè)定類別數(shù)為2。計(jì)算CMDE時(shí)散布模式數(shù)應(yīng)小于信號(hào)長(zhǎng)度,而信號(hào)進(jìn)行粗粒化處理將會(huì)減少信號(hào)長(zhǎng)度,為了在分析不同尺度下CMDE時(shí)不受維數(shù)影響,將維數(shù)設(shè)定為2。設(shè)定延時(shí)是為了從不同時(shí)間分辨率下對(duì)信號(hào)進(jìn)行分析,但延時(shí)過(guò)大時(shí),信號(hào)的有用信息會(huì)被丟失且產(chǎn)生混疊,通過(guò)比較將CMDE的延時(shí)設(shè)定為2。

    由圖4(a)可以發(fā)現(xiàn),使用VMD-CMDE方法得到的變性與未變性組織熵值之間存在差值,但是在尺度為6、16、19時(shí)變性與未變性組織熵值之間的差值小。通過(guò)圖4(b)可以發(fā)現(xiàn),由MIVMD-CMDE得到的變性與未變性組織熵值差值比VMD-CMDE方法得到的熵值差值大,在尺度8到11時(shí),變性組織與未變性組織的熵值差較大,其中在尺度為9時(shí)差值最大。經(jīng)過(guò)上述分析,經(jīng)過(guò)MIVMD重構(gòu)的實(shí)際超聲背散射信號(hào)的CMDE值能用于識(shí)別生物組織變性識(shí)別。

    為了證明尺度為9時(shí)CMDE能較好的區(qū)分變性組織與未變性組織,分別計(jì)算了尺度為9時(shí),經(jīng)MIVMD處理的16組實(shí)驗(yàn)的變性組織與未變性組織的超聲回波信號(hào)的CMDE均值。同時(shí)與文獻(xiàn)[16-18]中的近似熵(Approximate entropy,ApEn)、樣本熵(Sample Entropy,SE)、模糊熵(FuzzyEntropy,FE)進(jìn)行比較。如圖5所示,可以發(fā)現(xiàn),部分樣本組的SE與ApEn中出現(xiàn)變性時(shí)的熵值比未變性時(shí)的熵值小。所有樣本組的變性與未變性的FE不存在交疊,但是部分樣本組的變性與未變性的熵值差較小。所有樣本組的變性與未變性的CMDE值不存在交疊;相較于FE,變性與未變性的CEMD差值較大。變性組織的超聲回波信號(hào)的CMDE均值比未變性組織的超聲回波信號(hào)的CMDE均值高出0.110 8,約為8.12%。

    圖4 不同信號(hào)處理方式時(shí)尺度對(duì)CMDE值的影響

    圖5 變性與未變性組織的回波信號(hào)的不同熵

    將所有特征參量視為未知參量進(jìn)行GK聚類得到未變性與變性兩個(gè)聚類中心,計(jì)算CMDE與兩個(gè)聚類中心的貼近度,根據(jù)貼近度判斷樣本是否變性,并分別計(jì)算MIVMD-SE-GK、MIVMD-ApEn-GK、MIVMD-FE-GK和MIVMD-CMDE-GK四種方法的PC、PE、XB及變性識(shí)別率,結(jié)果如表1所示。

    表1 特征參數(shù)對(duì)聚類效果的影響

    通過(guò)對(duì)MIVMD-SE-GK、MIVMD-ApEn-GK、MIVMD-FE-GK和MIVMD-CMDE-GK方法分析可以發(fā)現(xiàn),MIVMD-ApEn-GK方法的PC指數(shù)最小且PE指數(shù)最大,聚類劃分在四種方法中最模糊,XB指數(shù)最大說(shuō)明MIVMD-ApEn-GK方法的類間分離最小,聚類離散程度最大。與MIVMD-ApEn-GK方法相比,MIVMD-SE-GK和MIVMD-FE-GK方法的PC、PE、XB評(píng)價(jià)指數(shù)都有所改善。在四種方式中,MIVMD-CMDE-GK方法的PC指數(shù)最大且PE指數(shù)最小,聚類劃分在四種方法中最清晰,XB指數(shù)最小說(shuō)明MIVMD-CMDE-GK方法的類間分離最大,聚類離散程度最小。在四種方法中,MIVMD-CMDE-GK方法的識(shí)別率最高,說(shuō)明CMDE對(duì)變性生物組織的識(shí)別能力比ApEn、SE、FE要強(qiáng)。

    3 結(jié)論

    超聲治療中采集的超聲背散射信號(hào)含有大量的噪聲,且噪聲頻率范圍廣,對(duì)判斷生物組織是否變性造成了影響。針對(duì)上述問題提出了MIVMD,同時(shí)采用了CMDE算法對(duì)超聲背散射回波信號(hào)進(jìn)行特征提取并結(jié)合GK模糊聚類進(jìn)行生物組織變性識(shí)別。通過(guò)仿真和實(shí)例驗(yàn)證可以得到一下結(jié)論:①M(fèi)IVMD分解重構(gòu)的去噪能力比EMD和VMD更好。②利用CMDE對(duì)處理過(guò)后的超聲背散射信號(hào)進(jìn)行了特征提取,結(jié)果表明:當(dāng)延時(shí)、類別、嵌入維數(shù)均設(shè)定為2時(shí),變性組織超聲背散射信號(hào)的熵值要高于未變性組織超聲背散射信號(hào)的熵值。經(jīng)MIVMD處理的超聲背散射信號(hào)的CMDE值在尺度為9時(shí)兩者差別最大。③對(duì)MIVMD分解重構(gòu)處理的超聲背散射信號(hào)進(jìn)行SE、ApEn、FE和CMDE特征提取并進(jìn)行GK聚類發(fā)現(xiàn),CMDE有更好地表征生物組織是否變性的能力。

    猜你喜歡
    背散射變性信噪比
    晉州市大成變性淀粉有限公司
    基于深度學(xué)習(xí)的無(wú)人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    征兵“驚艷”
    低信噪比下LFMCW信號(hào)調(diào)頻參數(shù)估計(jì)
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    當(dāng)變性女遇見變性男 一種奇妙的感覺產(chǎn)生了
    基于PSO-GRG的背散射模式掃描電鏡的數(shù)字處理及應(yīng)用
    變性淀粉在酸奶中的應(yīng)用
    西藏科技(2015年10期)2015-09-26 12:10:16
    保持信噪比的相位分解反褶積方法研究
    小型移動(dòng)背散射X射線安全檢查設(shè)備簡(jiǎn)介
    日本wwww免费看| 午夜福利在线免费观看网站| 亚洲第一av免费看| 天天躁狠狠躁夜夜躁狠狠躁| 成人三级黄色视频| 天天影视国产精品| 国产精品国产av在线观看| 十八禁网站免费在线| 欧美日韩乱码在线| 十八禁人妻一区二区| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器 | 国产精品一区二区在线不卡| 最近最新中文字幕大全电影3 | 丰满的人妻完整版| 1024香蕉在线观看| 精品国产美女av久久久久小说| av有码第一页| 亚洲欧美日韩无卡精品| 欧美老熟妇乱子伦牲交| 视频区欧美日本亚洲| 51午夜福利影视在线观看| 嫁个100分男人电影在线观看| 久久久精品国产亚洲av高清涩受| 国产精品永久免费网站| av国产精品久久久久影院| 欧美乱色亚洲激情| 久久精品91蜜桃| 99久久人妻综合| 国产欧美日韩一区二区三区在线| 热re99久久精品国产66热6| 国产亚洲精品第一综合不卡| a级片在线免费高清观看视频| 淫秽高清视频在线观看| 亚洲精华国产精华精| 成人18禁在线播放| 人妻丰满熟妇av一区二区三区| 精品卡一卡二卡四卡免费| 久久久久国产精品人妻aⅴ院| 久久久久国产精品人妻aⅴ院| 妹子高潮喷水视频| 亚洲少妇的诱惑av| 久久国产乱子伦精品免费另类| 国产免费av片在线观看野外av| 91成年电影在线观看| 久久99一区二区三区| 亚洲人成电影免费在线| 精品久久久久久成人av| 在线播放国产精品三级| 91精品国产国语对白视频| 国产精品乱码一区二三区的特点 | 自线自在国产av| 日本 av在线| 老汉色av国产亚洲站长工具| 中文字幕另类日韩欧美亚洲嫩草| 中国美女看黄片| 精品久久久久久电影网| 国内久久婷婷六月综合欲色啪| 亚洲精品久久午夜乱码| 91在线观看av| 亚洲成人免费电影在线观看| 91成人精品电影| 午夜福利,免费看| 久久中文字幕人妻熟女| 亚洲人成电影观看| 少妇的丰满在线观看| 精品国产乱子伦一区二区三区| 欧美乱色亚洲激情| 老司机靠b影院| 国产不卡一卡二| 天堂√8在线中文| 18禁美女被吸乳视频| 欧美精品一区二区免费开放| 婷婷丁香在线五月| 窝窝影院91人妻| 日韩精品中文字幕看吧| 亚洲成人免费电影在线观看| 欧美日本亚洲视频在线播放| 亚洲激情在线av| 国产精品av久久久久免费| 亚洲一区二区三区欧美精品| 日本撒尿小便嘘嘘汇集6| 成人影院久久| 无限看片的www在线观看| 12—13女人毛片做爰片一| 久久精品国产亚洲av香蕉五月| 俄罗斯特黄特色一大片| 青草久久国产| 久久久国产成人精品二区 | 88av欧美| 日韩人妻精品一区2区三区| 一级作爱视频免费观看| 亚洲欧洲精品一区二区精品久久久| 国产伦一二天堂av在线观看| 久久精品91蜜桃| 巨乳人妻的诱惑在线观看| 亚洲中文av在线| 午夜福利在线观看吧| aaaaa片日本免费| 国产aⅴ精品一区二区三区波| 中文字幕精品免费在线观看视频| 亚洲精品在线美女| 亚洲国产精品sss在线观看 | 青草久久国产| 熟女少妇亚洲综合色aaa.| 国产成+人综合+亚洲专区| 女人被躁到高潮嗷嗷叫费观| 男人舔女人下体高潮全视频| 女性生殖器流出的白浆| 最近最新免费中文字幕在线| 老司机午夜十八禁免费视频| 久久人人爽av亚洲精品天堂| 亚洲欧美精品综合一区二区三区| 久久欧美精品欧美久久欧美| 亚洲男人天堂网一区| 亚洲av第一区精品v没综合| 在线观看日韩欧美| 黑人猛操日本美女一级片| 国产精品乱码一区二三区的特点 | 欧美av亚洲av综合av国产av| 黄片小视频在线播放| 久久精品国产亚洲av高清一级| 制服诱惑二区| 色尼玛亚洲综合影院| 久久精品亚洲精品国产色婷小说| 一级a爱视频在线免费观看| 韩国精品一区二区三区| 精品乱码久久久久久99久播| 日日干狠狠操夜夜爽| 亚洲一卡2卡3卡4卡5卡精品中文| 美女扒开内裤让男人捅视频| 午夜福利影视在线免费观看| 好男人电影高清在线观看| 他把我摸到了高潮在线观看| 欧美在线一区亚洲| 免费在线观看黄色视频的| 欧美成人免费av一区二区三区| 亚洲av片天天在线观看| 丝袜在线中文字幕| avwww免费| 啦啦啦 在线观看视频| av超薄肉色丝袜交足视频| www国产在线视频色| 亚洲 欧美一区二区三区| 一级a爱视频在线免费观看| 国产又爽黄色视频| 俄罗斯特黄特色一大片| 成人三级做爰电影| 啦啦啦在线免费观看视频4| 国产精品98久久久久久宅男小说| 最好的美女福利视频网| 亚洲中文av在线| 国产熟女午夜一区二区三区| 婷婷丁香在线五月| 久久精品亚洲精品国产色婷小说| 久久欧美精品欧美久久欧美| 19禁男女啪啪无遮挡网站| 国产av又大| 精品国产乱子伦一区二区三区| 亚洲五月色婷婷综合| 日韩欧美一区二区三区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 精品国产乱子伦一区二区三区| 亚洲五月色婷婷综合| 久久精品国产亚洲av香蕉五月| 免费人成视频x8x8入口观看| 精品欧美一区二区三区在线| 99国产精品99久久久久| 亚洲 国产 在线| 啦啦啦 在线观看视频| 久9热在线精品视频| 久久国产精品人妻蜜桃| 国产视频一区二区在线看| avwww免费| av网站免费在线观看视频| 色综合欧美亚洲国产小说| 搡老熟女国产l中国老女人| 日韩精品青青久久久久久| 亚洲av日韩精品久久久久久密| 欧美人与性动交α欧美精品济南到| 高清在线国产一区| 国产精品国产高清国产av| 日韩 欧美 亚洲 中文字幕| 国产成人精品无人区| 亚洲精品国产一区二区精华液| 欧美国产精品va在线观看不卡| 国产亚洲精品第一综合不卡| 女性生殖器流出的白浆| 99热国产这里只有精品6| 久久青草综合色| 麻豆久久精品国产亚洲av | 欧美最黄视频在线播放免费 | 久久中文字幕一级| 操美女的视频在线观看| 亚洲色图综合在线观看| 欧美日韩亚洲高清精品| 露出奶头的视频| 老司机午夜十八禁免费视频| 91九色精品人成在线观看| 亚洲自拍偷在线| 日本免费a在线| 极品人妻少妇av视频| 欧美乱码精品一区二区三区| av天堂在线播放| 久久中文字幕人妻熟女| 十分钟在线观看高清视频www| 在线十欧美十亚洲十日本专区| 一级a爱片免费观看的视频| 一夜夜www| svipshipincom国产片| 精品国产超薄肉色丝袜足j| 深夜精品福利| 91国产中文字幕| 久久国产精品人妻蜜桃| 久久久久久免费高清国产稀缺| 一二三四社区在线视频社区8| 一级作爱视频免费观看| 一区二区三区国产精品乱码| 国产男靠女视频免费网站| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 一边摸一边抽搐一进一出视频| 在线观看一区二区三区| 国产精品久久久久久人妻精品电影| svipshipincom国产片| 欧美成人午夜精品| 欧美老熟妇乱子伦牲交| 如日韩欧美国产精品一区二区三区| 日日夜夜操网爽| 国产男靠女视频免费网站| 日日干狠狠操夜夜爽| 黑人猛操日本美女一级片| 亚洲 欧美 日韩 在线 免费| 亚洲精品美女久久久久99蜜臀| av福利片在线| 中文字幕av电影在线播放| 亚洲少妇的诱惑av| 日日干狠狠操夜夜爽| 国产激情久久老熟女| 久久国产精品影院| 亚洲片人在线观看| 1024视频免费在线观看| 久久午夜综合久久蜜桃| 亚洲成人免费电影在线观看| 十八禁网站免费在线| 精品一区二区三区四区五区乱码| 久久精品国产综合久久久| 亚洲少妇的诱惑av| 一夜夜www| 伦理电影免费视频| 色哟哟哟哟哟哟| 免费高清视频大片| 黄色片一级片一级黄色片| 午夜福利影视在线免费观看| 美女午夜性视频免费| 1024视频免费在线观看| 免费观看精品视频网站| 婷婷精品国产亚洲av在线| 国产成人精品无人区| 搡老乐熟女国产| 国产欧美日韩一区二区三区在线| 叶爱在线成人免费视频播放| 亚洲国产中文字幕在线视频| 亚洲欧美日韩另类电影网站| 长腿黑丝高跟| 精品熟女少妇八av免费久了| 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 亚洲aⅴ乱码一区二区在线播放 | 一区二区三区国产精品乱码| 国产在线观看jvid| 亚洲一区二区三区欧美精品| 亚洲精品久久午夜乱码| 国产精品综合久久久久久久免费 | 亚洲久久久国产精品| 国产高清国产精品国产三级| 午夜激情av网站| 18禁国产床啪视频网站| 真人一进一出gif抽搐免费| 午夜日韩欧美国产| 成人三级黄色视频| 久久久久国内视频| 日日夜夜操网爽| 在线观看一区二区三区激情| 又大又爽又粗| 免费在线观看日本一区| 久久人人97超碰香蕉20202| 久久中文字幕人妻熟女| 美女扒开内裤让男人捅视频| 丰满饥渴人妻一区二区三| 视频区欧美日本亚洲| 欧美日韩亚洲国产一区二区在线观看| 中文亚洲av片在线观看爽| 90打野战视频偷拍视频| 99久久综合精品五月天人人| 国产视频一区二区在线看| 亚洲av日韩精品久久久久久密| 校园春色视频在线观看| 欧美日韩黄片免| 男女之事视频高清在线观看| 丰满迷人的少妇在线观看| 最好的美女福利视频网| 久久精品91无色码中文字幕| 精品国产乱子伦一区二区三区| 村上凉子中文字幕在线| 亚洲aⅴ乱码一区二区在线播放 | 国产精品影院久久| 日韩人妻精品一区2区三区| 99精品在免费线老司机午夜| 天天影视国产精品| 一区二区三区精品91| 欧美乱码精品一区二区三区| 丰满的人妻完整版| 女警被强在线播放| 水蜜桃什么品种好| 在线天堂中文资源库| 大码成人一级视频| 成人特级黄色片久久久久久久| 欧美激情极品国产一区二区三区| 老司机福利观看| 成人av一区二区三区在线看| 男女下面插进去视频免费观看| 淫妇啪啪啪对白视频| 欧美日本亚洲视频在线播放| 嫁个100分男人电影在线观看| 亚洲欧美一区二区三区久久| 免费在线观看完整版高清| av片东京热男人的天堂| 侵犯人妻中文字幕一二三四区| 久久久久久大精品| 国产av一区在线观看免费| 日韩欧美国产一区二区入口| 麻豆av在线久日| 搡老熟女国产l中国老女人| 宅男免费午夜| 免费不卡黄色视频| 亚洲专区国产一区二区| 亚洲av日韩精品久久久久久密| 国产片内射在线| 日本五十路高清| 两性夫妻黄色片| 老司机亚洲免费影院| 国产欧美日韩精品亚洲av| 色综合欧美亚洲国产小说| 悠悠久久av| 性色av乱码一区二区三区2| 淫妇啪啪啪对白视频| 国产成人系列免费观看| 两性夫妻黄色片| 欧美激情久久久久久爽电影 | 在线看a的网站| 黄色毛片三级朝国网站| 日本 av在线| 99久久精品国产亚洲精品| 热re99久久精品国产66热6| 大型av网站在线播放| 国产1区2区3区精品| 精品国产美女av久久久久小说| 国产区一区二久久| 亚洲国产欧美日韩在线播放| 久久久久久久精品吃奶| 免费观看精品视频网站| 美女午夜性视频免费| 国产国语露脸激情在线看| 国产精品98久久久久久宅男小说| 欧美激情极品国产一区二区三区| 亚洲男人的天堂狠狠| 神马国产精品三级电影在线观看 | 亚洲色图综合在线观看| 精品一品国产午夜福利视频| 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸 | 日韩欧美在线二视频| av天堂在线播放| 桃红色精品国产亚洲av| 国产免费男女视频| 免费高清视频大片| 亚洲精品国产一区二区精华液| 99久久人妻综合| 欧美精品一区二区免费开放| 欧美国产精品va在线观看不卡| 亚洲av成人av| 高清在线国产一区| 亚洲欧美日韩无卡精品| 日韩有码中文字幕| 黑人巨大精品欧美一区二区mp4| 国产伦人伦偷精品视频| 可以在线观看毛片的网站| 欧美黑人精品巨大| 美女 人体艺术 gogo| 精品国产亚洲在线| cao死你这个sao货| 热re99久久国产66热| a级片在线免费高清观看视频| 99久久精品国产亚洲精品| 国产亚洲av高清不卡| 中文字幕另类日韩欧美亚洲嫩草| 精品久久久久久久久久免费视频 | 在线看a的网站| 在线观看免费高清a一片| 日韩欧美三级三区| 大码成人一级视频| 国产男靠女视频免费网站| 午夜精品久久久久久毛片777| 亚洲精品成人av观看孕妇| 黄频高清免费视频| 一级毛片精品| 18禁国产床啪视频网站| 中文字幕色久视频| 欧美一区二区精品小视频在线| 18美女黄网站色大片免费观看| 韩国av一区二区三区四区| 亚洲av成人av| 午夜成年电影在线免费观看| 亚洲少妇的诱惑av| 久久香蕉激情| 成人特级黄色片久久久久久久| 一个人观看的视频www高清免费观看 | 91九色精品人成在线观看| 在线永久观看黄色视频| 女人精品久久久久毛片| 纯流量卡能插随身wifi吗| 国产黄色免费在线视频| 高清欧美精品videossex| 国产熟女xx| 岛国视频午夜一区免费看| 黑人猛操日本美女一级片| 久久国产精品男人的天堂亚洲| av福利片在线| 精品福利观看| 一夜夜www| 香蕉久久夜色| 啦啦啦 在线观看视频| 天堂影院成人在线观看| 欧美在线黄色| 十分钟在线观看高清视频www| 亚洲 国产 在线| 中文欧美无线码| 国产精品99久久99久久久不卡| 日日干狠狠操夜夜爽| 久久精品亚洲精品国产色婷小说| 五月开心婷婷网| 欧美黑人精品巨大| 精品一区二区三区四区五区乱码| 热re99久久精品国产66热6| a级毛片在线看网站| 久久久久国产精品人妻aⅴ院| 欧美日韩亚洲综合一区二区三区_| 丁香六月欧美| 亚洲男人天堂网一区| 级片在线观看| 男女下面进入的视频免费午夜 | 国产av在哪里看| 精品无人区乱码1区二区| 久久久久久人人人人人| 中文欧美无线码| 亚洲人成网站在线播放欧美日韩| 免费女性裸体啪啪无遮挡网站| netflix在线观看网站| 99在线人妻在线中文字幕| 淫妇啪啪啪对白视频| 91麻豆精品激情在线观看国产 | 男女做爰动态图高潮gif福利片 | 久久久久九九精品影院| 又紧又爽又黄一区二区| 交换朋友夫妻互换小说| 国产伦一二天堂av在线观看| 一个人免费在线观看的高清视频| 窝窝影院91人妻| 亚洲久久久国产精品| 久久国产精品人妻蜜桃| 丰满的人妻完整版| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费视频网站a站| 国产精品亚洲av一区麻豆| 好男人电影高清在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线美女| 男人的好看免费观看在线视频 | 宅男免费午夜| 国产人伦9x9x在线观看| 午夜福利一区二区在线看| 亚洲avbb在线观看| 波多野结衣av一区二区av| 国产亚洲欧美98| 免费看a级黄色片| x7x7x7水蜜桃| 日本五十路高清| 免费在线观看亚洲国产| 欧美亚洲日本最大视频资源| 中文欧美无线码| 妹子高潮喷水视频| 欧美乱妇无乱码| 日本黄色视频三级网站网址| 99riav亚洲国产免费| 97人妻天天添夜夜摸| 在线看a的网站| 又黄又爽又免费观看的视频| 高清黄色对白视频在线免费看| 亚洲专区字幕在线| 黄频高清免费视频| 国产精品乱码一区二三区的特点 | 高清毛片免费观看视频网站 | 亚洲av日韩精品久久久久久密| 级片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 久热爱精品视频在线9| 91成人精品电影| 男女之事视频高清在线观看| 中文字幕色久视频| 好看av亚洲va欧美ⅴa在| 天天躁狠狠躁夜夜躁狠狠躁| 香蕉丝袜av| 亚洲午夜理论影院| 深夜精品福利| 国产成人av激情在线播放| 91精品三级在线观看| 亚洲精品av麻豆狂野| 成人免费观看视频高清| 国产成人精品久久二区二区免费| 国产免费现黄频在线看| 99国产精品99久久久久| 国产成人精品在线电影| 99久久99久久久精品蜜桃| 99国产综合亚洲精品| 国产成+人综合+亚洲专区| 亚洲人成77777在线视频| 搡老熟女国产l中国老女人| 日本免费a在线| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品久久久久5区| 亚洲一区二区三区不卡视频| 极品教师在线免费播放| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 久久久久国产精品人妻aⅴ院| 电影成人av| 久热爱精品视频在线9| 黄片大片在线免费观看| 老汉色∧v一级毛片| 99精品久久久久人妻精品| 国产深夜福利视频在线观看| e午夜精品久久久久久久| 亚洲国产毛片av蜜桃av| 老司机福利观看| 男女做爰动态图高潮gif福利片 | 欧美人与性动交α欧美精品济南到| 麻豆久久精品国产亚洲av | 成人国产一区最新在线观看| 美女 人体艺术 gogo| 无遮挡黄片免费观看| 高清av免费在线| 一进一出好大好爽视频| 久久午夜综合久久蜜桃| 亚洲成av片中文字幕在线观看| 日日夜夜操网爽| 国产精品av久久久久免费| a级毛片黄视频| 十八禁人妻一区二区| 久久婷婷成人综合色麻豆| 国产精品美女特级片免费视频播放器 | 中文字幕色久视频| 国产主播在线观看一区二区| av天堂久久9| 亚洲久久久国产精品| 亚洲av成人一区二区三| 最近最新中文字幕大全免费视频| 欧美日韩精品网址| 色精品久久人妻99蜜桃| 成人特级黄色片久久久久久久| av福利片在线| 亚洲美女黄片视频| 亚洲国产中文字幕在线视频| 日本欧美视频一区| 黄色丝袜av网址大全| 日韩三级视频一区二区三区| 韩国av一区二区三区四区| 69av精品久久久久久| 日韩中文字幕欧美一区二区| 免费在线观看黄色视频的| 久久亚洲真实| 午夜两性在线视频| av福利片在线| 精品久久久久久,| 国产精品 欧美亚洲| www.999成人在线观看| 日日夜夜操网爽| 91av网站免费观看| 久热爱精品视频在线9| 麻豆国产av国片精品| 日本五十路高清| 法律面前人人平等表现在哪些方面| 美女高潮喷水抽搐中文字幕| 精品高清国产在线一区| 日日干狠狠操夜夜爽| 天堂俺去俺来也www色官网| 午夜老司机福利片| av在线播放免费不卡| 91大片在线观看| 久久久久久久午夜电影 | 日韩欧美免费精品| 国产乱人伦免费视频| 午夜免费激情av| 波多野结衣高清无吗| 校园春色视频在线观看| 99热只有精品国产| 中文字幕av电影在线播放| 精品免费久久久久久久清纯| 国产亚洲精品久久久久久毛片| bbb黄色大片| 99re在线观看精品视频| 淫妇啪啪啪对白视频| 丝袜美腿诱惑在线|