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

    用于心臟電活動(dòng)成像的空間濾波器輸出噪聲抑制方法?

    2018-09-06 07:33:18周大方張樹(shù)林蔣式勤
    物理學(xué)報(bào) 2018年15期
    關(guān)鍵詞:譜估計(jì)增益濾波器

    周大方 張樹(shù)林 蔣式勤?

    1)(同濟(jì)大學(xué)電子與信息工程學(xué)院,上海 201804)2)(中國(guó)科學(xué)院上海微系統(tǒng)與信息技術(shù)研究所,信息功能材料國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200050)(2018年2月6日收到;2018年4月8日收到修改稿)

    1 引 言

    波束成形方法是一種空間濾波技術(shù).它可以通過(guò)構(gòu)造一個(gè)空間濾波器,提取該位置上感興趣的源強(qiáng)度信息,同時(shí)抑制來(lái)自其他位置上源的影響[1,2].這種方法采用分布源模型估計(jì)分布電流源偶極矩的強(qiáng)度及其空間位置信息,即分布源的空間譜.相比采用單電流偶極子源模型,波束成形方法解決了估計(jì)多源的計(jì)算問(wèn)題,可以細(xì)致地描述生物電活動(dòng),但是,面臨如何準(zhǔn)確重建多源的問(wèn)題.近年來(lái),最小方差波束成形(minimum variance beamforming,MVB)[1?7]被用于重建分布等效電流偶極子(equivalent current dipole,ECD)源的研究,例如,識(shí)別預(yù)激綜合征(wolf f-parkinson-white)病人心臟房室的旁路傳導(dǎo)[8,9],以及定位房顫(atrial fibrillation,AF)病人的AF病灶[10]和利用磁場(chǎng)圖(magnetocardiogram,MCG)仿真數(shù)據(jù)的心臟輪廓成像[7].在重建磁場(chǎng)分布電流源時(shí),MVB通過(guò)自適應(yīng)的空間濾波技術(shù),最小化空間濾波器輸出總功率及歸一化噪聲空間譜強(qiáng)度,降低了測(cè)量噪聲對(duì)電流源空間譜估計(jì)的影響.可以發(fā)現(xiàn),如果在MVB方法的基礎(chǔ)上,再針對(duì)性地約束空間濾波器輸出的噪聲功率增益,可以進(jìn)一步提高分布電流源空間譜估計(jì)的源分辨能力,從而增強(qiáng)心臟電活動(dòng)磁成像的分辨率.

    本文提出了一種可抑制空間濾波器輸出噪聲功率增益(suppressing spatial f i lter output noisepower gain,SONG)的波束成形方法.用一個(gè)對(duì)空間濾波器輸出功率有影響的低跡的半正定矩陣構(gòu)造一種新的濾波器權(quán)矩陣,可以約束空間濾波器的輸出噪聲功率的增益.該低跡的半正定矩陣滿(mǎn)足特征值不大于1,且矩陣的跡小于其階數(shù).這樣,可以提高分布電流源空間譜估計(jì)的分辨率.文中通過(guò)理論分析和仿真實(shí)驗(yàn)比較了SONG和MVB方法.給出了用兩個(gè)健康人的36通道MCG數(shù)據(jù)得到的心臟電活動(dòng)成像.結(jié)果表明,SONG方法分辨電流源的能力較強(qiáng),能夠觀察到Rpeak時(shí)刻健康人的心室內(nèi)有較強(qiáng)電活動(dòng)等明顯的電生理特征.

    2 波束成形方法

    2.1 問(wèn)題的提出

    假設(shè)第j個(gè)單位電流偶極子源的位置和方向分別是rj=(xj,yj,zj)(j=1,2,···,n)和[1,0,0]T,[0,1,0]T或[0,0,1]T.它產(chǎn)生的理想磁場(chǎng)列向量為lX,j,lY,j或 lZ,j,相應(yīng)的導(dǎo)聯(lián)場(chǎng)矩陣表示rj處電流源的測(cè)量靈敏度.已知心臟磁場(chǎng)的測(cè)量數(shù)據(jù)b(t)[11,12],求電流源偶極矩q(t,rj)的逆問(wèn)題[13,14],可用線性方程表示[15,16]:

    其中,b(t)=[b1(t),b2(t),···,bc(t)]T表示t時(shí)刻c個(gè)測(cè)量通道的磁場(chǎng)向量;v(t)是t時(shí)刻的測(cè)量噪聲向量;等效電流源的偶極矩q(t,rj)=q(t,rj)η(t,rj),j=1,2,···,n, 其中n是電流源的數(shù)目.源偶極矩的強(qiáng)度是標(biāo)量

    表示單位向量.rj=(xj,yj,zj)是第j個(gè)電流源位置的坐標(biāo).下文中,q(t,rj),η(t,rj),q(t,rj),v(t)和b(t)將簡(jiǎn)寫(xiě)為qj,ηj,qj,v和b.

    空間濾波是一種常用的分布源重建方法.將心臟磁場(chǎng)測(cè)量數(shù)據(jù)b(t)作為空間濾波器的輸入,估計(jì)的分布源的偶極矩作為輸出.空間濾波可用加權(quán)的線性運(yùn)算表示:

    其中,W(t,rj)是空間濾波的權(quán)矩陣.(2)式可簡(jiǎn)寫(xiě)為

    MVB方法的基本原理是先用空間濾波技術(shù)重建心臟的分布電流偶極子源.然后,根據(jù)可描述電流源偶極矩平均強(qiáng)度的分布電流源空間譜估計(jì),對(duì)心臟電流源成像.電流源偶極矩是源電流密度在鄰域的體積分,它與該位置上的源電流密度大小,也就是電流源強(qiáng)度有關(guān)[13,14].

    式中,前一項(xiàng)為空間濾波器輸出的所有電流源偶極矩qi的功率和;后一項(xiàng)為空間濾波器輸出的噪聲功率,其中為噪聲功率的增益.對(duì)開(kāi)方,可得估計(jì)的電流源偶極矩平均強(qiáng)度空間譜:

    2.2 SONG波束成形方法

    本文提出了一種SONG波束成形方法.令空間濾波權(quán)矩陣

    可以證明,其中的c階實(shí)對(duì)稱(chēng)矩陣V=E(bbT)是一個(gè)對(duì)空間濾波器輸出功率有影響的半正定矩陣,滿(mǎn)足“矩陣的跡低于其階數(shù),且其特征值不大于1”,(以下簡(jiǎn)稱(chēng)V是低跡半正定陣).因此,SONG方法中輸出噪聲功率增益可表示為

    令任意矩陣Γ =Wj,MVB,也可以證明,存在不等式

    當(dāng)V=E(bbT)時(shí),由(10)和(11)式可知

    也就是說(shuō),SONG方法可以或更好地降低空間濾波器的輸出噪聲功率增益.將E(bbT)=I代入(8)和(9)式可知,SONG的噪聲空間譜強(qiáng)度與MVB相同,均等于1.綜上,SONG方法不僅可以約束噪聲空間譜對(duì)分布源空間譜估計(jì)的影響,還可約束空間濾波器輸出噪聲功率的增益.相比MVB,SONG方法可以提高估計(jì)空間譜的分辨率.

    2.3 SONG方法的分析與比較

    采用波束成形方法,電流源的空間譜估計(jì)決定了電活動(dòng)成像的分辨能力.由于多電流源重建問(wèn)題比較復(fù)雜,本文比較了SONG和MVB單電流源重建的源分辨率.

    由(2)式可知,當(dāng)單電流源偶極矩方向已知時(shí),空間濾波器的源偶極矩估計(jì)可以退化為源偶極矩的強(qiáng)度估計(jì)是退化后空間濾波器的權(quán)向量.由(7)式可知,單源在任意位置上產(chǎn)生的估計(jì)空間譜強(qiáng)度為假設(shè)單源S位置rs上的估計(jì)空間譜強(qiáng)度為b ?s.為了分析估計(jì)空間譜的單源分辨率,定義任意位置rj(j=1,2,···,n)上的點(diǎn)擴(kuò)散函數(shù)(point spread function)為(rj)[4],可用歸一化后的估計(jì)空間譜強(qiáng)度簡(jiǎn)單表示為

    由(9)式可知,單源重建時(shí),SONG方法的權(quán)向量為

    式中,

    其中,lj是電流源產(chǎn)生的磁場(chǎng)列向量.

    由(7)和(13)式可得SONG方法估計(jì)的空間譜

    其中

    將(15)和(17)式代入(14)式,可得SONG方法在任意位置rj(j=1,2,···,n)的估計(jì)空間譜強(qiáng)度

    歸一化后的點(diǎn)擴(kuò)散函數(shù)為

    同理可得,MVB單源重建的估計(jì)空間譜強(qiáng)度和歸一化后的點(diǎn)擴(kuò)散函數(shù)

    比較(19)和(20)式,有

    通常,所有測(cè)量通道上的理想磁場(chǎng)信號(hào)平均功率大于噪聲平均功率[11,12].所以,可令α=根據(jù)施瓦茲不等式性質(zhì),可知由此,從(20)和(21)式可以得到

    可見(jiàn),當(dāng)rj=rs時(shí),單源位置上點(diǎn)擴(kuò)散函數(shù)和相等,為最大值1. 當(dāng)rj?=rs時(shí),其他空間位置上,SONG方法的點(diǎn)擴(kuò)散函數(shù)比MVB的小.點(diǎn)擴(kuò)散函數(shù)?j可以反映單源對(duì)空間其他位置的估計(jì)空間譜強(qiáng)度影響大小(rj?=rs)的取值小,說(shuō)明單源對(duì)鄰域的影響擴(kuò)散小.因此,相比MVB,SONG方法對(duì)單源的空間分辨率較高.

    3 仿真與實(shí)驗(yàn)

    3.1 仿真數(shù)據(jù)和結(jié)果

    相關(guān)電流源(correlated sources)是比較難分辨的,所以,文中利用仿真的磁場(chǎng)數(shù)據(jù),比較了SONG與MVB方法估計(jì)相關(guān)電流源的能力.

    假定軀干G0={(x,y,z)|x∈[?12.5,12.5],y∈[?12.5,12.5],z∈[3,20]}(cm).隨機(jī)給定兩個(gè)相關(guān)源的位置坐標(biāo)為(6.5,?2.5,11)和(?2.5,6.5,11)(cm),它們的相關(guān)系數(shù)為0.9824.G0被劃分為間距1 cm的10625個(gè)體素(voxel).并用ECD源模型產(chǎn)生兩組仿真的36通道Z軸測(cè)量磁場(chǎng)數(shù)據(jù),采樣頻率為1 kHz[7,12],如圖1(a)和圖1(b)所示.假設(shè)其中測(cè)量噪聲分別為根均方(root-mean-square,rms)信噪比(signal-to-noise ratio,SNR)20 dBrms和10 dBrms的高斯白噪聲[7,12].產(chǎn)生仿真的磁場(chǎng)數(shù)據(jù)和進(jìn)行電流源空間譜估計(jì),都要用到等效電流源的導(dǎo)聯(lián)場(chǎng)矩陣.文中采用水平分層導(dǎo)體作為軀干模型,并通過(guò)解磁場(chǎng)的正問(wèn)題求導(dǎo)聯(lián)場(chǎng)矩陣[7,13,14].

    圖1(c)和圖1(d)是SONG和MVB在XY平面(z=11 cm)上歸一化后的空間譜估計(jì)強(qiáng)度的等高線圖,等高線的步長(zhǎng)為1%.結(jié)果表明,SNR對(duì)相關(guān)電流源的空間譜估計(jì)有影響.與MVB相比,用SONG方法估計(jì)的相關(guān)源強(qiáng)度比鄰域的估計(jì)源強(qiáng)度有明顯的增強(qiáng).MVB方法中,表示濾波器的輸出噪聲功率增益,本文在此基礎(chǔ)上給出了SONG與MVB方法的輸出噪聲功率增益比

    如圖1(e)和圖1(f)所示.其中,橫坐標(biāo)是分布源位置的索引,縱坐標(biāo)是增益比β.將(12)式代入β可知β>0.當(dāng)β>0時(shí),說(shuō)明SONG方法抑制噪聲的效果比MVB好.圖1表明,測(cè)量信噪比SNR分別為20和10 dB時(shí),用SONG重建相關(guān)源時(shí),濾波器輸出噪聲增益比β分別為584—622和592—623.在SNR=10 dB時(shí)抑制噪聲增益的效果較SNR=20 dB時(shí)更明顯.因此,圖1(c)和圖1(d)表明,SONG方法有較強(qiáng)的相關(guān)源分辨能力,并且在SNR較低時(shí),其源分辨能力更加明顯.

    圖1 (a),(b)兩種SNR情況下,兩個(gè)相關(guān)電流源產(chǎn)生的仿真磁場(chǎng)數(shù)據(jù)(紅色虛線表示重建相關(guān)源的時(shí)刻);(c),(d)SONG和MVB在XY平面(z=11 cm)上估計(jì)空間譜強(qiáng)度的等高線圖(×號(hào)表示相關(guān)源的給定位置);(e),(f)用SONG方法的空間濾波器輸出噪聲功率增益比βFig.1.(a),(b)Simulated magnetic data generated by two correlated current sources with noise,respectively.The red dashed line denotes the time of source reconstruction.(c),(d)Estimated spatial spectrum intensity contour on XY plane(z=11 cm)using SONG and MVB.The black cross signs×denote the given locations of two correlated sources.(e),(f)Ratio β of noise-power gain of spatial f i lter output using SONG method.

    3.2 心臟電活動(dòng)成像

    分布電流偶極子源的偶極矩強(qiáng)度可以反映分布電流的強(qiáng)度,因此,可利用分布源空間譜估計(jì)的強(qiáng)度對(duì)心臟電活動(dòng)成像.文中還用兩個(gè)健康人的心臟磁場(chǎng)數(shù)據(jù)比較了SONG和MVB成像的效果.在求導(dǎo)聯(lián)場(chǎng)矩陣時(shí),采用水平分層導(dǎo)體作為軀干模型[7,13,14].用這種模型時(shí),導(dǎo)體邊界上的單位法向量均平行于心臟測(cè)量磁場(chǎng)的單位向量,所以,求解導(dǎo)聯(lián)場(chǎng)矩陣時(shí),軀干體電導(dǎo)的影響可以忽略[7,13,14].

    圖2是沿Z軸測(cè)量的兩個(gè)健康人的36通道單周期心臟磁場(chǎng)曲線.測(cè)量平面為25 cm×25 cm.36通道的SNR約為10—20 dB,帶寬0.01—100 Hz,采樣頻率1 kHz[12].

    圖2中心臟磁場(chǎng)信號(hào)的Rpeak時(shí)刻對(duì)應(yīng)心室的除極期.這時(shí)健康人心室電活動(dòng)較強(qiáng),心房的相對(duì)較弱[17],比較容易識(shí)別[11,12].因此,文中比較了SONG和MVB兩種方法的成像結(jié)果.將測(cè)量平面與健康人核磁共振影像(magnetic resonance imaging,MRI)的心臟冠狀位(coronal view)、水平位(transverse view)和矢狀位(sagittal view)圖的坐標(biāo)配準(zhǔn),然后,用MRI中心臟的位置作為房室位置的參考.圖3(a)和圖3(b)給出了歸一化后的分布源空間譜估計(jì)強(qiáng)度的等高線圖.其中,冠狀位視角的XY面(z=10 cm)、水平位視角的XZ面(y=2.5 cm)和矢狀位視角的Y Z面(x=0.5 cm)的交點(diǎn)(0.5,2.5,10)cm位于心室內(nèi),用藍(lán)色的小方框表示,譜強(qiáng)度等高線的步長(zhǎng)為1%.圖3(a)和圖3(b)表明,SONG方法的成像結(jié)果能夠反映健康人Rpeak時(shí)刻電活動(dòng)的特征.因?yàn)镽peak時(shí)刻,健康人的心房除極已結(jié)束,心室正處于除極期.根據(jù)圖中色標(biāo),可以觀察到心室內(nèi)黃色表示的電活動(dòng)強(qiáng)度明顯比心房?jī)?nèi)紅色表示的電活動(dòng)強(qiáng)度高,以及心室鄰域的電活動(dòng)相對(duì)較強(qiáng)[7,11,12,17,18].由于SONG增加了心臟內(nèi)外分布源空間譜估計(jì)強(qiáng)度的差異,心臟分布電流源的分辨率提高了.用MVB的成像結(jié)果相對(duì)模糊,特征不明顯.

    圖3(c)和圖3(d)的結(jié)果表明,Rpeak時(shí)刻用SONG方法,兩個(gè)健康人的濾波器輸出噪聲增益比β分別為206—228和215—232.也就是說(shuō),當(dāng)實(shí)測(cè)心磁數(shù)據(jù)SNR為10—20 dB時(shí),相比MVB,SONG能夠更好地抑制空間濾波器輸出噪聲.如圖3(a)和圖3(b)所示,SONG方法的源分辨能力以及成像效果相對(duì)較好.

    圖2 (a),(b)兩個(gè)健康人的36通道單周期心臟磁場(chǎng)數(shù)據(jù)及對(duì)應(yīng)的Rpeak時(shí)刻Fig.2.(a),(b)The 36-channel MCG data of single-cycle from two healthy subjects,as well as the time-points of Rpeakfor MCG imaging.

    4 討 論

    由(12)式可見(jiàn),SONG波束成形方法可以約束空間濾波器輸出噪聲的功率增益.用構(gòu)造濾波權(quán)矩陣Wj,SONG類(lèi)似的方法,還可以構(gòu)造其他的濾波權(quán)矩陣.雖然理論上其他空間濾波權(quán)矩陣對(duì)應(yīng)的噪聲空間譜強(qiáng)度也恒等于1,并可以約束空間濾波輸出噪聲功率增益,但是,仿真結(jié)果表明,其他濾波權(quán)矩陣會(huì)使兩個(gè)相關(guān)電流源定位到它們的中間位置.因此,必須利用測(cè)量信號(hào)的二階特征矩陣E(bbT)構(gòu)造空間濾波權(quán)矩陣[2,7].因?yàn)闇y(cè)量信號(hào)矩陣E(bbT)可以反映電流源的相關(guān)性.圖1中用SONG方法,兩個(gè)給定的相關(guān)源可以準(zhǔn)確定位.

    我們還利用MCG仿真數(shù)據(jù),研究了重建分布源的最小范數(shù)空間濾波(minimum norm spatial f i ltering,MN)方法[19,20],這種非自適應(yīng)的空間濾波方法,采用了Tikhonov正則化技術(shù).仿真結(jié)果表明,有可利用的心臟三維輪廓時(shí),MN的分布電流源成像結(jié)果較好[19].參考文獻(xiàn)[19]的圖3和圖4給出了重建的16個(gè)等效電流源,與本文圖3中SONG方法的成像結(jié)果類(lèi)似.

    圖3 (a),(b)兩個(gè)健康人的分布源空間譜估計(jì)強(qiáng)度的等高線圖;(c),(d)兩個(gè)健康人的空間濾波器輸出噪聲功率增益比βFig.3.(a),(b)Contour of the estimated intensity of distributed source spatial spectrum of two healthy subjects;(c),(d)spatial f i lter output noise-power gain ratio β of two healthy subjects.

    5 結(jié) 論

    本文在研究MVB方法的基礎(chǔ)上,提出了一種用于心臟電活動(dòng)成像的可抑制空間濾波器輸出噪聲功率增益的波束成形方法.該方法利用一種低跡半正定矩陣構(gòu)造了一個(gè)濾波器權(quán)矩陣,可以降低空間濾波器輸出的噪聲功率增益,提高重建分布電流源偶極矩強(qiáng)度分辨率即分布電流源空間譜估計(jì)的源分辨能力,從而增強(qiáng)心臟電活動(dòng)磁成像的分辨率.仿真實(shí)驗(yàn)和分析比較的結(jié)果表明,SONG方法優(yōu)于MVB方法.當(dāng)心磁信號(hào)的SNR不低于10 dB時(shí),采用該方法可以提高心臟電活動(dòng)成像的效果,將有助于相關(guān)的醫(yī)學(xué)研究和應(yīng)用.

    感謝中國(guó)科學(xué)院上海微系統(tǒng)與信息技術(shù)研究所的張懿教授、謝曉明教授和孔祥燕教授及其團(tuán)隊(duì)為本研究提供可用的心臟磁場(chǎng)數(shù)據(jù)與核磁共振影像數(shù)據(jù),以及有關(guān)的技術(shù)交流.

    猜你喜歡
    譜估計(jì)增益濾波器
    基于增益調(diào)度與光滑切換的傾轉(zhuǎn)旋翼機(jī)最優(yōu)控制
    基于單片機(jī)的程控增益放大器設(shè)計(jì)
    電子制作(2019年19期)2019-11-23 08:41:36
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    基于Multisim10和AD603的程控增益放大器仿真研究
    電子制作(2018年19期)2018-11-14 02:37:02
    開(kāi)關(guān)電源EMI濾波器的應(yīng)用方法探討
    電子制作(2018年16期)2018-09-26 03:26:50
    基于MATLAB實(shí)現(xiàn)的AR模型功率譜估計(jì)
    基于Canny振蕩抑制準(zhǔn)則的改進(jìn)匹配濾波器
    基于TMS320C6678的SAR方位向預(yù)濾波器的并行實(shí)現(xiàn)
    經(jīng)典功率譜估計(jì)方法的研究
    Welch譜估計(jì)的隨機(jī)誤差與置信度
    国产乱来视频区| 大陆偷拍与自拍| 国产色爽女视频免费观看| 国产精品久久久久久久久免| 色婷婷久久久亚洲欧美| 插阴视频在线观看视频| 一二三四中文在线观看免费高清| 舔av片在线| 亚洲国产精品一区三区| 日韩强制内射视频| 成人18禁高潮啪啪吃奶动态图 | 王馨瑶露胸无遮挡在线观看| 免费av中文字幕在线| 日日撸夜夜添| 一级毛片 在线播放| 国产真实伦视频高清在线观看| 免费高清在线观看视频在线观看| 久久久久久九九精品二区国产| 国产免费福利视频在线观看| 99re6热这里在线精品视频| 人妻一区二区av| 91狼人影院| 久久鲁丝午夜福利片| 欧美人与善性xxx| 免费播放大片免费观看视频在线观看| 精品国产露脸久久av麻豆| 欧美最新免费一区二区三区| 大陆偷拍与自拍| 高清视频免费观看一区二区| 久久av网站| 欧美激情国产日韩精品一区| 人妻制服诱惑在线中文字幕| 免费看不卡的av| 国产欧美亚洲国产| 插逼视频在线观看| 天堂中文最新版在线下载| 国产高潮美女av| 亚洲av中文av极速乱| 欧美日韩视频精品一区| 精华霜和精华液先用哪个| 精品人妻视频免费看| 26uuu在线亚洲综合色| 国产黄频视频在线观看| 国产精品av视频在线免费观看| 日韩制服骚丝袜av| 国产精品国产三级国产av玫瑰| 亚州av有码| 高清av免费在线| a级毛片免费高清观看在线播放| 国精品久久久久久国模美| 日本av手机在线免费观看| 久久韩国三级中文字幕| 99久国产av精品国产电影| 精品一区二区三区视频在线| 久久久色成人| 春色校园在线视频观看| 麻豆成人午夜福利视频| 免费在线观看成人毛片| 国产欧美日韩一区二区三区在线 | 国产精品国产三级专区第一集| 亚洲精品乱码久久久久久按摩| 亚洲精品乱码久久久久久按摩| 欧美亚洲 丝袜 人妻 在线| 精华霜和精华液先用哪个| 七月丁香在线播放| 国产国拍精品亚洲av在线观看| 精品久久国产蜜桃| 观看美女的网站| 欧美日韩视频高清一区二区三区二| 91久久精品国产一区二区成人| 一二三四中文在线观看免费高清| av国产精品久久久久影院| 一区二区三区免费毛片| 97热精品久久久久久| videossex国产| 大香蕉97超碰在线| 王馨瑶露胸无遮挡在线观看| 日韩亚洲欧美综合| 国产片特级美女逼逼视频| 少妇人妻精品综合一区二区| 美女视频免费永久观看网站| 免费黄网站久久成人精品| 99热这里只有精品一区| 亚洲精品中文字幕在线视频 | 欧美一区二区亚洲| 色视频在线一区二区三区| 视频中文字幕在线观看| 久久综合国产亚洲精品| 丰满迷人的少妇在线观看| 青青草视频在线视频观看| 亚洲欧美日韩卡通动漫| 青青草视频在线视频观看| 嫩草影院入口| 亚洲四区av| 在现免费观看毛片| 国产精品久久久久久精品电影小说 | 我要看日韩黄色一级片| 麻豆国产97在线/欧美| 久久ye,这里只有精品| 成人漫画全彩无遮挡| 久久 成人 亚洲| 午夜老司机福利剧场| 欧美xxxx黑人xx丫x性爽| 18禁裸乳无遮挡免费网站照片| 美女cb高潮喷水在线观看| 亚洲欧美一区二区三区国产| 一级av片app| 一级av片app| 免费黄网站久久成人精品| 亚洲av中文字字幕乱码综合| 看免费成人av毛片| 国产免费一级a男人的天堂| 日本色播在线视频| 国产精品av视频在线免费观看| 国产精品三级大全| 深爱激情五月婷婷| 亚洲色图综合在线观看| 一本一本综合久久| 免费在线观看成人毛片| 久久精品国产亚洲av天美| 一个人免费看片子| 国产午夜精品一二区理论片| 国产精品99久久久久久久久| 性高湖久久久久久久久免费观看| 欧美精品一区二区免费开放| 国产伦精品一区二区三区视频9| 少妇人妻久久综合中文| 99久久精品热视频| 一本—道久久a久久精品蜜桃钙片| 亚洲内射少妇av| 精品少妇久久久久久888优播| 王馨瑶露胸无遮挡在线观看| 亚洲欧美中文字幕日韩二区| 国产亚洲一区二区精品| 国产91av在线免费观看| 亚洲国产毛片av蜜桃av| 赤兔流量卡办理| 国产在线视频一区二区| 男女免费视频国产| 国产精品成人在线| 成人国产av品久久久| 97超碰精品成人国产| 黄色日韩在线| 91午夜精品亚洲一区二区三区| 亚洲欧美成人综合另类久久久| 亚洲自偷自拍三级| 久久精品久久精品一区二区三区| 国产精品99久久99久久久不卡 | 联通29元200g的流量卡| 99re6热这里在线精品视频| 国产有黄有色有爽视频| av在线观看视频网站免费| 亚洲av福利一区| 久久99热这里只频精品6学生| 97热精品久久久久久| 久久久色成人| 在线亚洲精品国产二区图片欧美 | 草草在线视频免费看| 高清av免费在线| 校园人妻丝袜中文字幕| 欧美日韩国产mv在线观看视频 | 一个人看的www免费观看视频| 国产爱豆传媒在线观看| 欧美日韩精品成人综合77777| 国产亚洲欧美精品永久| 综合色丁香网| 国产精品嫩草影院av在线观看| 蜜桃亚洲精品一区二区三区| 精品亚洲乱码少妇综合久久| 大又大粗又爽又黄少妇毛片口| 国产一区二区三区综合在线观看 | 男女边摸边吃奶| 26uuu在线亚洲综合色| 国产精品免费大片| 91在线精品国自产拍蜜月| 精品一区二区三区视频在线| 极品教师在线视频| 欧美一级a爱片免费观看看| 亚洲在久久综合| 91久久精品电影网| 又大又黄又爽视频免费| 一边亲一边摸免费视频| 国产黄片视频在线免费观看| 色婷婷久久久亚洲欧美| 女的被弄到高潮叫床怎么办| 日韩电影二区| 欧美日韩综合久久久久久| 国产极品天堂在线| 美女国产视频在线观看| 久久影院123| 国产亚洲精品久久久com| 国产爽快片一区二区三区| 精品一品国产午夜福利视频| 成人亚洲精品一区在线观看 | 夫妻午夜视频| 久久青草综合色| 国产成人aa在线观看| 日日啪夜夜爽| 亚洲欧美日韩东京热| 国产女主播在线喷水免费视频网站| 2022亚洲国产成人精品| 26uuu在线亚洲综合色| 欧美另类一区| 丰满少妇做爰视频| 成人二区视频| 一级毛片aaaaaa免费看小| 最近中文字幕2019免费版| 26uuu在线亚洲综合色| 国产高清国产精品国产三级 | 亚洲经典国产精华液单| 国产黄频视频在线观看| 国产乱人视频| 五月开心婷婷网| 亚洲av国产av综合av卡| 我的女老师完整版在线观看| 亚洲精品久久午夜乱码| 午夜激情福利司机影院| 国产成人a区在线观看| av福利片在线观看| 自拍偷自拍亚洲精品老妇| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美成人a在线观看| 国产精品国产三级国产专区5o| 99热这里只有是精品50| 女性生殖器流出的白浆| 在线观看免费日韩欧美大片 | 久久久精品免费免费高清| 国产黄片美女视频| 男的添女的下面高潮视频| 精品一区二区免费观看| 好男人视频免费观看在线| 国产精品蜜桃在线观看| 国产在线男女| 天天躁夜夜躁狠狠久久av| 十分钟在线观看高清视频www | 久久精品国产自在天天线| 麻豆国产97在线/欧美| 午夜老司机福利剧场| 18禁动态无遮挡网站| 秋霞伦理黄片| 欧美最新免费一区二区三区| 中文欧美无线码| 多毛熟女@视频| 大片电影免费在线观看免费| 久久精品国产亚洲av涩爱| 亚洲一级一片aⅴ在线观看| av专区在线播放| 亚洲精品日本国产第一区| 精品久久久久久久末码| 人人妻人人添人人爽欧美一区卜 | 国产免费视频播放在线视频| 久久97久久精品| 亚洲aⅴ乱码一区二区在线播放| 久久99热这里只频精品6学生| 久久久精品94久久精品| 国国产精品蜜臀av免费| 成人特级av手机在线观看| 日韩强制内射视频| 一级毛片久久久久久久久女| 亚洲国产最新在线播放| 精品一区二区三区视频在线| 日韩一区二区三区影片| 午夜精品国产一区二区电影| 国产爽快片一区二区三区| 伊人久久国产一区二区| 亚洲精品国产av蜜桃| 日日啪夜夜爽| 精品一区二区三卡| 日本黄色日本黄色录像| 精品熟女少妇av免费看| 这个男人来自地球电影免费观看 | 成人一区二区视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 晚上一个人看的免费电影| 男女无遮挡免费网站观看| 在线精品无人区一区二区三 | 男女国产视频网站| 人人妻人人添人人爽欧美一区卜 | 国产人妻一区二区三区在| 国产乱来视频区| 久久久久久久久久久丰满| 欧美一级a爱片免费观看看| 不卡视频在线观看欧美| 色网站视频免费| 天天躁夜夜躁狠狠久久av| 国产精品精品国产色婷婷| 高清av免费在线| 国产精品嫩草影院av在线观看| 啦啦啦啦在线视频资源| 国产深夜福利视频在线观看| 国产精品不卡视频一区二区| 亚洲成色77777| 中文字幕av成人在线电影| av国产免费在线观看| 男女边吃奶边做爰视频| 一区二区三区四区激情视频| 热99国产精品久久久久久7| 久久久久久人妻| 黄色一级大片看看| 久久精品人妻少妇| 一本—道久久a久久精品蜜桃钙片| 夫妻午夜视频| 在现免费观看毛片| 中文字幕免费在线视频6| 人妻一区二区av| 国产女主播在线喷水免费视频网站| 亚洲国产精品专区欧美| 99九九线精品视频在线观看视频| 春色校园在线视频观看| 欧美日韩亚洲高清精品| 80岁老熟妇乱子伦牲交| 99视频精品全部免费 在线| 久久久久久九九精品二区国产| 精品久久久久久电影网| 少妇裸体淫交视频免费看高清| 国产精品熟女久久久久浪| 久久久久久久国产电影| 午夜视频国产福利| 精品国产露脸久久av麻豆| 观看av在线不卡| 亚洲精品国产色婷婷电影| 一级毛片我不卡| 狠狠精品人妻久久久久久综合| 国产毛片在线视频| 欧美xxⅹ黑人| 国产伦精品一区二区三区视频9| 超碰av人人做人人爽久久| 99久国产av精品国产电影| 极品少妇高潮喷水抽搐| 国产精品99久久久久久久久| 亚洲精品亚洲一区二区| 免费久久久久久久精品成人欧美视频 | 特大巨黑吊av在线直播| 成人二区视频| 久久久久久人妻| 男女无遮挡免费网站观看| 欧美三级亚洲精品| 看非洲黑人一级黄片| 中国国产av一级| a 毛片基地| 99热这里只有是精品50| 日本免费在线观看一区| 国产一区二区在线观看日韩| 国产黄片视频在线免费观看| 国产视频内射| 亚洲丝袜综合中文字幕| 黄色日韩在线| 2022亚洲国产成人精品| 国产精品人妻久久久影院| 午夜福利视频精品| 干丝袜人妻中文字幕| 日韩成人伦理影院| 91aial.com中文字幕在线观看| 卡戴珊不雅视频在线播放| 亚洲精品国产av蜜桃| 插逼视频在线观看| 又粗又硬又长又爽又黄的视频| 97精品久久久久久久久久精品| 国产在线男女| 亚洲美女黄色视频免费看| 国产免费视频播放在线视频| 最近最新中文字幕大全电影3| 亚洲av日韩在线播放| 国产精品爽爽va在线观看网站| 超碰97精品在线观看| 女性生殖器流出的白浆| 久久久久精品久久久久真实原创| 午夜老司机福利剧场| kizo精华| 99热国产这里只有精品6| 一本一本综合久久| 国产在线男女| 久久精品久久久久久噜噜老黄| 日韩欧美精品免费久久| 国产av一区二区精品久久 | 久久99热这里只有精品18| 中文字幕精品免费在线观看视频 | 亚洲精品一二三| 午夜福利在线在线| av国产精品久久久久影院| 婷婷色av中文字幕| 国产探花极品一区二区| 亚洲成人中文字幕在线播放| 免费观看性生交大片5| 国产久久久一区二区三区| 99视频精品全部免费 在线| 少妇的逼好多水| 三级国产精品片| www.av在线官网国产| 99热网站在线观看| 国产亚洲精品久久久com| 国产精品不卡视频一区二区| 欧美日韩视频精品一区| 欧美xxxx黑人xx丫x性爽| 精品一品国产午夜福利视频| 国产一区二区三区av在线| 卡戴珊不雅视频在线播放| 成年女人在线观看亚洲视频| 日本色播在线视频| 国语对白做爰xxxⅹ性视频网站| 老女人水多毛片| 我的老师免费观看完整版| 天堂8中文在线网| 精品亚洲成国产av| 国产爱豆传媒在线观看| 综合色丁香网| 美女内射精品一级片tv| 国产成人a区在线观看| 国产精品偷伦视频观看了| 欧美成人一区二区免费高清观看| 国产色爽女视频免费观看| 亚洲真实伦在线观看| av在线播放精品| 精品亚洲成国产av| 亚洲精品久久久久久婷婷小说| 美女福利国产在线 | 在线天堂最新版资源| 搡老乐熟女国产| 欧美高清性xxxxhd video| 三级国产精品欧美在线观看| 精品少妇久久久久久888优播| 国产色爽女视频免费观看| 亚洲真实伦在线观看| 成人无遮挡网站| 成人午夜精彩视频在线观看| 国产精品偷伦视频观看了| 亚洲高清免费不卡视频| 日韩欧美一区视频在线观看 | freevideosex欧美| 在线观看免费高清a一片| 秋霞伦理黄片| 亚洲成人中文字幕在线播放| 亚洲av免费高清在线观看| 国产成人freesex在线| 99热网站在线观看| 日本黄色片子视频| 色网站视频免费| 我要看日韩黄色一级片| 国产在线一区二区三区精| 亚洲av国产av综合av卡| 女人十人毛片免费观看3o分钟| 99热这里只有是精品50| 日本与韩国留学比较| 插阴视频在线观看视频| 日韩在线高清观看一区二区三区| 成人国产麻豆网| 下体分泌物呈黄色| 成人免费观看视频高清| 日韩av免费高清视频| 国产欧美日韩一区二区三区在线 | 下体分泌物呈黄色| 日本免费在线观看一区| 少妇人妻 视频| 日本一二三区视频观看| 一级爰片在线观看| 3wmmmm亚洲av在线观看| 欧美xxxx性猛交bbbb| 久久女婷五月综合色啪小说| 亚洲精品国产色婷婷电影| 亚洲av成人精品一区久久| 免费播放大片免费观看视频在线观看| 在线免费观看不下载黄p国产| 色综合色国产| 久久久久人妻精品一区果冻| 乱码一卡2卡4卡精品| 国产 一区精品| 人妻少妇偷人精品九色| 黑丝袜美女国产一区| 日本爱情动作片www.在线观看| 国产成人午夜福利电影在线观看| 久久久国产一区二区| 少妇精品久久久久久久| 日本wwww免费看| 免费看日本二区| 汤姆久久久久久久影院中文字幕| 一本一本综合久久| 在线天堂最新版资源| 国产成人a∨麻豆精品| 亚洲,欧美,日韩| xxx大片免费视频| 又爽又黄a免费视频| 国内揄拍国产精品人妻在线| 日韩一区二区三区影片| 国产女主播在线喷水免费视频网站| 丰满少妇做爰视频| 成人高潮视频无遮挡免费网站| 我要看日韩黄色一级片| 久久婷婷青草| 一本—道久久a久久精品蜜桃钙片| 欧美极品一区二区三区四区| 亚洲av国产av综合av卡| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载| 丰满少妇做爰视频| 国产精品一区二区三区四区免费观看| 久久久久久久精品精品| 熟女av电影| 中国国产av一级| 国产精品人妻久久久久久| 国产亚洲欧美精品永久| 中文字幕人妻熟人妻熟丝袜美| 亚洲美女黄色视频免费看| 久久久久国产精品人妻一区二区| 91久久精品电影网| 女性被躁到高潮视频| 亚洲电影在线观看av| 午夜福利影视在线免费观看| 久久女婷五月综合色啪小说| av在线播放精品| 亚洲最大成人中文| 人人妻人人添人人爽欧美一区卜 | 亚洲精品日韩在线中文字幕| 简卡轻食公司| 交换朋友夫妻互换小说| 热99国产精品久久久久久7| 欧美三级亚洲精品| 婷婷色综合www| 亚洲电影在线观看av| 亚洲高清免费不卡视频| 国产精品一区二区三区四区免费观看| 久久久久久久精品精品| 国产成人freesex在线| 日本av手机在线免费观看| 午夜视频国产福利| 日本猛色少妇xxxxx猛交久久| 成年女人在线观看亚洲视频| 亚洲内射少妇av| 免费人成在线观看视频色| 久久久久视频综合| 涩涩av久久男人的天堂| 亚洲怡红院男人天堂| 亚洲不卡免费看| 少妇人妻一区二区三区视频| 91精品国产九色| 五月伊人婷婷丁香| 久久精品久久久久久久性| 91aial.com中文字幕在线观看| 美女高潮的动态| 自拍偷自拍亚洲精品老妇| 天堂俺去俺来也www色官网| 亚洲精品乱码久久久久久按摩| 水蜜桃什么品种好| 亚洲激情五月婷婷啪啪| 熟女人妻精品中文字幕| 我的老师免费观看完整版| 毛片女人毛片| 高清黄色对白视频在线免费看 | 国产在视频线精品| 人体艺术视频欧美日本| 婷婷色综合www| av在线蜜桃| 在线 av 中文字幕| av网站免费在线观看视频| 国产精品99久久99久久久不卡 | tube8黄色片| 免费不卡的大黄色大毛片视频在线观看| 老女人水多毛片| 国产精品人妻久久久影院| 国产 一区 欧美 日韩| 在线观看三级黄色| 成人综合一区亚洲| 国产亚洲欧美精品永久| 最近最新中文字幕大全电影3| 亚洲精品色激情综合| 国产亚洲av片在线观看秒播厂| 中国三级夫妇交换| 亚洲av成人精品一区久久| 又粗又硬又长又爽又黄的视频| 国产伦精品一区二区三区四那| 久久久国产一区二区| 久久毛片免费看一区二区三区| 久久久久人妻精品一区果冻| 久久精品国产亚洲av涩爱| 天天躁日日操中文字幕| 99视频精品全部免费 在线| 一本一本综合久久| 精品人妻偷拍中文字幕| 国产美女午夜福利| 久久久亚洲精品成人影院| 天堂俺去俺来也www色官网| 观看av在线不卡| 高清毛片免费看| 亚洲精品乱久久久久久| 日韩欧美一区视频在线观看 | 99九九线精品视频在线观看视频| 国产探花极品一区二区| 在线观看免费高清a一片| 最近中文字幕2019免费版| 男人爽女人下面视频在线观看| 精品国产露脸久久av麻豆| 色婷婷久久久亚洲欧美| 久久久久久九九精品二区国产| 日韩一区二区三区影片| 国产 一区精品| 精品亚洲乱码少妇综合久久| 久久久久久久久大av| 国产淫语在线视频| 成人18禁高潮啪啪吃奶动态图 | 国产高清国产精品国产三级 | 久久久久人妻精品一区果冻| 亚洲在久久综合| 亚洲色图综合在线观看| 人妻制服诱惑在线中文字幕| 欧美xxxx黑人xx丫x性爽| 制服丝袜香蕉在线| 51国产日韩欧美| 亚洲真实伦在线观看| 亚洲人成网站在线播| 日韩大片免费观看网站| 国产欧美日韩精品一区二区| 亚洲av电影在线观看一区二区三区| 亚洲性久久影院| 天天躁日日操中文字幕| 寂寞人妻少妇视频99o| 狂野欧美激情性bbbbbb|