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

    新型分裂步長時域有限差分法

    2015-08-14 21:39:23林智參班濤
    現(xiàn)代電子技術(shù) 2015年15期

    林智參+班濤

    摘 要: 提出一種新型的分裂步長時域有限差分(NSS?FDTD)法,并對其數(shù)值色散進(jìn)行分析。該方法基于Split?Step方案和Crank?Nicolson方案,采用新的矩陣分解形式,與傳統(tǒng)的FDTD算法、SS?FDTD算法相比,減少了計(jì)算復(fù)雜度。新型算法的推導(dǎo)程序簡單,且具有良好的數(shù)值色散特性,還加入了一階Mur吸收邊界條件,給出一階Mur吸收邊界差分方程。將數(shù)值實(shí)驗(yàn)的結(jié)果和傳統(tǒng)FDTD方法及理論值進(jìn)行比較,數(shù)值結(jié)果一致性較好。

    關(guān)鍵詞: 時域有限差分法; 分裂步長; Split?Step方案; 數(shù)值色散

    中圖分類號: TN802?34; O441 文獻(xiàn)標(biāo)識碼: A 文章編號: 1004?373X(2015)15?0117?03

    New method of finite difference time domain for split step

    LIN Zhican, BAN Tao

    (South China Normal University, Guangzhou 510006, China)

    Abstract: A new split step finite difference time domain (NSS?FDTD) algorithm is presented, and its numerical dispersion is analyzed. The method is based on the schemes of Split?Step and Crank?Nicolson, adopted new matrix decomposition form. Compared with traditional algorithms of FDTD and SS?FDTD, the proposed algorithm can reduce computational complexity, and has simple deduction procedure and better numerical dispersion characteristic. The first?order Mur absorbing boundary condition is added in this paper, and its difference equation is presented. The numerical experiment results were compared with traditional FDTD method and theoretical values. The consistence of numerical results is better.

    Keywords: FDTD method; split step; Split?Step scheme; numerical dispersion

    0 引 言

    時域有限差分法(Finite Difference Time Domain,F(xiàn)DTD)是一種簡單直觀的全波分析時域算法[1?3],該方法以Yee氏立體網(wǎng)格作為電磁場離散單元,將麥克斯韋方程轉(zhuǎn)化為差分方程,能夠方便有效地結(jié)合計(jì)算機(jī)技術(shù)處理復(fù)雜的電磁場問題,目前已經(jīng)在電磁學(xué)的各個領(lǐng)域中得到了廣泛應(yīng)用。然而,應(yīng)用傳統(tǒng)的FDTD方法也明顯體現(xiàn)出其不足之處,為減小差分近似帶來的數(shù)值色散[4],空間網(wǎng)格尺寸必須遠(yuǎn)小于波長, 這樣反而增加運(yùn)算負(fù)擔(dān),因而,出現(xiàn)了多鐘FDTD的改進(jìn)方法[5?7]。本文以TM波為例,提出一種基于Split?Step方案[8]和Crank?Nicolson方案[9]的時域有限差分法,所提出的算法采用新的矩陣分解方式,簡化計(jì)算復(fù)雜度、減少差分近似所帶來的數(shù)值色散。

    1 NSS?FDTD算法理論推導(dǎo)

    考慮空間一無源區(qū)域,在均勻無耗、各向同性介質(zhì)中,介電常數(shù)為[ε,]磁導(dǎo)率為[μ,]二維TM波Maxwell微分方程組的矩陣形式如下:

    [?u?t=Mu] (1)

    式中:[u=[Ez,Hx,Hy]T, M=0-1ε??y1ε??x-1μ??y001μ??x00。]

    將矩陣[M]分解成四個子矩陣,分別記為[A2,][B2,][C2,][D2,]矩陣形式分別如下:

    [A=001ε??x0001μ??x00 B=0-1ε??y0-1μ??y00000]

    [C=0-1ε??y00001μ??x00 D=001ε??x-1μ??y00000]

    公式(1)可以寫為:

    [?u?t=A2u+B2u+C2u+D2u] (2)

    利用Split?Step方案,將式(2)分解成四個子方程式來求解。同時,從[n~][n+1]時間步,將一個時間步長等間隔分成四個子時間分步,即[n→n+14,][n+14→][n+24,][n+24→n+34]和[n+34→n+1,]得到以下四個子矩陣方程式:

    分步1:

    [?u?t=4?A2u n→n+14] (3)

    分步2:

    [?u?t=4?B2u n+14→n+24] (4)

    分步3:

    [?u?t=4?C2u n+24→n+34] (5)

    分步4:

    [?u?t=4?D2u n+34→n+1] (6)

    利用Crank?Nicolson方案,對矩陣方程式(3)~式(6)的右端進(jìn)行近似,進(jìn)一步化簡得到以下形式:

    [I-Δt4Aun+14=I+Δt4Aun] (7)

    [I-Δt4Bun+24=I+Δt4Bun+14] (8)

    [I-Δt4Cun+34=I+Δt4Cun+24] (9)

    [I-Δt4Dun+1=I+Δt4Dun+34] (10)

    式中:[I]為[3×3]的單位矩陣。四個分步內(nèi)需要求解的方程式以第(1)步為例化簡如下:

    [Ezn+14i, j-Δt4εΔxHyn+14i+12, j-Hyn+14i-12, j= Ezni, j+Δt4εΔxHyni+12, j-Hyni-12, j] (11)

    [Hxn+14i, j+12=Hxni, j+12] (12)

    [Hyn+14i+12, j-Δt4μΔxEzn+14i+1, j-Ezn+14i, j= Hyni+12, j+Δt4μΔxEzni+1, j-Ezni, j] (13)

    在分步1內(nèi),只需要求解方程式(11)和(13),并且方程式(11)和(13)是互為耦合的方程式。將式(13)代入到式(11)中,消去[Hyn+14i+12, j,]得到關(guān)于[Ezn+14i, j]的三對角矩陣方程式,其具體形式如下:

    [1+Δt28μεΔx2Ezn+14i, j-Δt216μεΔx2Ezn+14i+1, j+Ezn+14i-1, j= 1-Δt28μεΔx2Ezni, j+Δt216μεΔx2Ezni+1, j+Ezni-1, j+ Δt2εΔxHyni+12, j-Hyni-12, j]

    由式(14)求出[Ezn+14i, j,]然后代入到式(13)中,此時式(13)為關(guān)于[Hyn+14i+12, j]的一個顯式方程,可直接進(jìn)行求解。因此,在分步1內(nèi)只需要求解一個隱式方程和一個顯式方程。分步2,3和4內(nèi)也采用與分步1內(nèi)相同的方式進(jìn)行處理,而分步3,4中得到的不是三對角矩陣方程式,是顯式方程,求解變得更簡單。

    2 數(shù)值色散分析

    利用Fourier方法[10]在第[n]個時間步內(nèi),場分量在空間區(qū)域內(nèi)的表達(dá)形式如下:

    [UnI, J=Une-j(kxIΔx+kyJΔy)] (15)

    將式(15)代入到式(7)~(10)中,并將式(7)~(10)進(jìn)一步整理為如下矩陣形式:

    [Un+14=A2dUn] (16)

    [Un+24=B2dUn+14] (17)

    [Un+34=C2dUn+24] (18)

    [Un+1=D2dUn+34] (19)

    將式(16)~(18)代入到式(19)中,得到一個完整時間步長內(nèi)的矩陣方程式:

    [Un+1=D2dC2dB2dA2dUn=ΩUn] (20)

    可以求得[Ω]的特征值,其結(jié)果如下:

    [λ1=1,λ2=λ3*=P+jQ2-P2Q] (21)

    其中:

    [P=-P4xP4yb4d4+16(P4x+P4y)b2d2+104P2xP2yb2d2-224bd(P2x+P2y)+128]

    [Q=8P2xP2yb2d2+32bd(P2x+P2y)+128Pα=-2ΔαsinkαΔα2,α=x,y;b=Δt2ε;d=Δt2μ ]

    利用von Neumann方法[11],假定一角頻率為[ω]的電磁波產(chǎn)生的電磁場滿足:

    [Enz=EzejωΔtn,Hnα=HαejωΔtn,α=x,y] (22)

    將式(22)代入到式(20)中,得:

    [ejωΔtI-ΩUn=0] (23)

    其中,[Un]與初始值[U0]相關(guān),其具體關(guān)系式如下:

    [Un=U0ejωΔtn] (24)

    為了使式(23)中 [Un]有非零解,[Un]的系數(shù)行列式的值應(yīng)為零,即:

    [det(ejωΔtI-Ω)=0] (25)

    由式(21)得[Ω]的特征值,可以得到 NSS?FDTD 算法的數(shù)值色散表達(dá)式,其形式如下:

    [tan2ωΔt2=P-QP+Q=bd(P4xP4yb3d3-16P4xbd-96P2xP2ybd-16P4ybd+256P2x+256P2y)P4xP4yb4d4-16(P4x+P4y)b2d2-112P2xP2yb2d2+192(P2x+P2y)bd-256]

    為了方便分析NSS?FDTD的數(shù)值色散特性,作如下定義:[S=cΔtΔx,Δx=Δy,N=λΔx,λ代表波長,][N]表示單位波長元胞數(shù),傳輸角度為θ。圖1為[S=0.5,][N=8]時歸一化數(shù)值相位速度隨傳輸角度θ的變化曲線。

    圖1 歸一化數(shù)值相位速度隨傳輸角度[θ]的變化曲線

    從圖1可以看出,NSS?FDTD算法的歸一化相位速度大于傳統(tǒng)的FDTD算法和傳統(tǒng)SS?FDTD算法的歸一化數(shù)值相位速度,并接近1,且曲線變化值在0.98~0.995范圍內(nèi),比較平緩,歸一化數(shù)值相位速度各向異性誤差也比較小。

    3 一階Mur吸收邊界條件

    在分步1內(nèi),電場分量[Ez]在[i=1,][i=Imax,][j=1]和[j=Jmax]上的一階Mur吸收邊界差分方程式為:

    [Ezn+141, j=Ezn2, j+vΔt-4ΔxvΔt+4ΔxEzn+142, j-Ezn1, j] (26)

    [Ezn+14Imax, j=EznImax-1, j+vΔt-4ΔxvΔt+4ΔxEzn+14Imax-1, j-EznImax, j](27)[Ezn+14i,1=Ezni,2+vΔt-4ΔxvΔt+4ΔxEzn+14i, 2-Ezni,1] (28)

    [Ezn+14i, Jmax=Ezni, Jmax-1+vΔt-4ΔxvΔt+4ΔxEzn+14i, Jmax-1-Ezni, Jmax] (29)

    在分步2,3和4內(nèi),電場分量[Ez]在[i=1,][i=Imax,][j=1]和[j=Jmax]上的一階Mur吸收邊界差分方程式與分步1內(nèi)的差分方程式類似,此處不再一一進(jìn)行展開說明。

    4 計(jì)算結(jié)果

    將NSS?FDTD算法運(yùn)算于尺寸101 cm×101 cm的自由空間,以一階Mur為邊界條件,以二維TM波為例,在中心區(qū)域[Ez]場分量上加正弦波激勵源[sin(2πft),]其中,[f=1.5] GHz。網(wǎng)格尺寸為[Δx=Δy=1 cm,]為激勵源最高頻率對應(yīng)波長的[120,]計(jì)算網(wǎng)格數(shù)為101×101。觀察點(diǎn)在中心區(qū)域和吸收邊界之間的中心位置。NSS?FDTD數(shù)值計(jì)算結(jié)果仿真如圖2~圖4所示。其中圖2,圖3兩圖為NSS?FDTD數(shù)值計(jì)算過程中[Ez]場分量的空間分布圖,圖4為傳統(tǒng)FDTD算法和NSS?FDTD算法比較圖。

    圖2 運(yùn)行50步時[Ez]的3D圖

    圖3 運(yùn)行200步時[Ez]的3D圖

    圖4 兩種FDTD算法在觀察點(diǎn)處的電場值[Ez]

    由圖可見,NSS?FDTD算法的計(jì)算結(jié)果與傳統(tǒng)的FDTD算法的計(jì)算結(jié)果吻合的很好,且符合電磁場理論,從而證實(shí)新型分裂步長時域有限差分法的可行性。

    5 結(jié) 論

    本文基于Split?Step和Crank?Nicolson方案提出了一種新型的二維FDTD算法,改進(jìn)算法采用新的分解形式,與傳統(tǒng)的SS?FDTD算法相比較,減少了計(jì)算復(fù)雜度,優(yōu)化了計(jì)算公式,使推導(dǎo)過程更簡單。結(jié)合算例使用Matlab對NSS?FDTD算法進(jìn)行編程分析,結(jié)果表明,該算法具有良好的預(yù)期效果。

    參考文獻(xiàn)

    [1] 葛德彪,閆玉波.電磁波時域有限差分方法[M].西安:西安電子科技大學(xué)出版社,2002.

    [2] 王秉中.計(jì)算電磁學(xué)[M].北京:科學(xué)出版社,2002.

    [3] YEE K S. Numerical solution of initial boundary value problems involving Maxwell′s equations in isotropic media [J]. IEEE Transactions on Antennas and Propagation, 1966, 14 (3): 302?307.

    [4] TAFLOVE A, HAGNESS S C. Computational electrodynamics: the finite?difference time?domain method [M]. 2nd ed. Boston: Artech House, 2000.

    [5] 徐利軍,衰乃昌.高階ADI?FDTD算法的數(shù)值色散分析[J].電子與信息學(xué)報(bào),2005,27(10):1662?1665.

    [6] 黨濤,鄭宏興.關(guān)于二維ADI?FDTD方法的數(shù)值色散分析[J].中國民航學(xué)院學(xué)報(bào),2004,22(2):42?46.

    [7] 夏冬,黨濤,鄭宏興.一維ADI?FDTD方法的數(shù)值色散分析[J].中國民航學(xué)院學(xué)報(bào),2005,23(2):38?41.

    [8] LEE J, FORNBERG B. A split step approach for the 3?D Maxwell′s equations [J]. Journal of Computational and Applied Mathematics, 2003, 158(2): 485?505.

    [9] SMITH G D. Numerical solution of partial differential equations: finite difference methods [M]. 3rd ed. Oxford: Oxford University Press, 1986.

    [10] SMITH G D. Numerical solution of partial differential equations [M]. Oxford: Oxford University Press, 1978.

    [11] PEREDA J A, VIELVA L A, VEGAS A, et al. Analyzing the stability of the FDTD technique by combining the Von Neumann method with the Routh?Hurwitz criterion [J]. IEEE Transactions on Microwave Theory and Techniques, 2001, 49(2): 377?381.

    国产乱来视频区| 欧美最新免费一区二区三区| av在线播放精品| 精品久久久久久电影网| 久久精品熟女亚洲av麻豆精品| 欧美变态另类bdsm刘玥| 色网站视频免费| 国产免费福利视频在线观看| 日韩av不卡免费在线播放| 人人妻人人看人人澡| 亚洲性久久影院| .国产精品久久| 中文天堂在线官网| 欧美日韩在线观看h| 一区二区av电影网| 国产成人精品婷婷| av在线天堂中文字幕| 亚洲aⅴ乱码一区二区在线播放| 嫩草影院入口| 国产精品一区www在线观看| 人妻系列 视频| 别揉我奶头 嗯啊视频| 国产免费又黄又爽又色| 日韩成人伦理影院| 国产高清不卡午夜福利| 日韩一本色道免费dvd| 日韩中字成人| 中文精品一卡2卡3卡4更新| 亚洲最大成人av| 性插视频无遮挡在线免费观看| 最近2019中文字幕mv第一页| 亚洲图色成人| 秋霞伦理黄片| 黄色配什么色好看| 在线观看av片永久免费下载| 又大又黄又爽视频免费| 97人妻精品一区二区三区麻豆| 亚洲精品自拍成人| 亚洲精品视频女| 午夜福利在线在线| 搡女人真爽免费视频火全软件| 各种免费的搞黄视频| 亚洲av电影在线观看一区二区三区 | 国产黄色视频一区二区在线观看| 丰满少妇做爰视频| 黄色视频在线播放观看不卡| 免费黄色在线免费观看| 99热这里只有是精品在线观看| 国产乱人偷精品视频| 国产成人精品一,二区| 波野结衣二区三区在线| 中文欧美无线码| 777米奇影视久久| 深爱激情五月婷婷| 18禁裸乳无遮挡动漫免费视频 | 亚洲欧美精品专区久久| 永久网站在线| 汤姆久久久久久久影院中文字幕| 一区二区三区精品91| 少妇人妻 视频| 少妇的逼好多水| 一区二区三区四区激情视频| 亚洲精品国产av成人精品| 午夜激情福利司机影院| 中文字幕制服av| 国产成人a区在线观看| 国产毛片在线视频| 欧美区成人在线视频| 国产精品国产三级专区第一集| 日韩伦理黄色片| 亚洲精品色激情综合| 欧美性猛交╳xxx乱大交人| 2021少妇久久久久久久久久久| 国产精品熟女久久久久浪| 99精国产麻豆久久婷婷| 午夜精品国产一区二区电影 | 高清毛片免费看| 日韩一区二区视频免费看| 天天一区二区日本电影三级| 免费av不卡在线播放| 亚洲av国产av综合av卡| 亚洲综合色惰| 舔av片在线| 日本一二三区视频观看| 成年免费大片在线观看| 久久精品国产亚洲av天美| 亚洲无线观看免费| av在线蜜桃| 日韩一本色道免费dvd| 婷婷色综合大香蕉| 久久久国产一区二区| 嫩草影院入口| 国产一区二区在线观看日韩| 亚洲国产av新网站| 精品一区在线观看国产| 亚洲精品影视一区二区三区av| 亚洲国产色片| 亚洲精品中文字幕在线视频 | 人妻制服诱惑在线中文字幕| av福利片在线观看| 国产伦在线观看视频一区| 国产伦在线观看视频一区| 观看免费一级毛片| 国产综合懂色| 99久久精品一区二区三区| 五月开心婷婷网| 交换朋友夫妻互换小说| av在线观看视频网站免费| av免费观看日本| 久久精品久久久久久噜噜老黄| 五月伊人婷婷丁香| 亚洲av福利一区| 久久韩国三级中文字幕| 亚洲精品第二区| 精品亚洲乱码少妇综合久久| 国产毛片在线视频| 国产极品天堂在线| 有码 亚洲区| 大片免费播放器 马上看| 色视频www国产| 国产色婷婷99| 99久久九九国产精品国产免费| 肉色欧美久久久久久久蜜桃 | 狂野欧美白嫩少妇大欣赏| 国产乱人偷精品视频| 国产探花极品一区二区| 日韩欧美 国产精品| 看非洲黑人一级黄片| 国产亚洲最大av| 热99国产精品久久久久久7| 成年av动漫网址| 91精品一卡2卡3卡4卡| 两个人的视频大全免费| 亚洲av.av天堂| 一二三四中文在线观看免费高清| 人妻系列 视频| 成人欧美大片| 看非洲黑人一级黄片| 午夜老司机福利剧场| 精品一区二区三卡| 欧美亚洲 丝袜 人妻 在线| 国精品久久久久久国模美| 又爽又黄a免费视频| 日韩欧美 国产精品| 国产精品精品国产色婷婷| 中文字幕制服av| av.在线天堂| 亚洲国产欧美人成| 别揉我奶头 嗯啊视频| 身体一侧抽搐| 亚洲欧美日韩卡通动漫| 成人一区二区视频在线观看| 日韩精品有码人妻一区| 久久久久国产精品人妻一区二区| 成人亚洲精品一区在线观看 | 少妇裸体淫交视频免费看高清| 女的被弄到高潮叫床怎么办| 三级国产精品欧美在线观看| 久久97久久精品| av在线老鸭窝| 亚洲激情五月婷婷啪啪| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 成人国产麻豆网| 男人和女人高潮做爰伦理| kizo精华| 高清日韩中文字幕在线| 国产成人aa在线观看| 久久久久久久久久成人| 亚洲av福利一区| 欧美丝袜亚洲另类| 国产男女超爽视频在线观看| 99九九线精品视频在线观看视频| 日本一本二区三区精品| 精华霜和精华液先用哪个| 一区二区三区乱码不卡18| 成人一区二区视频在线观看| 看非洲黑人一级黄片| 久久国内精品自在自线图片| 另类亚洲欧美激情| 少妇人妻 视频| 亚洲欧美中文字幕日韩二区| 亚洲经典国产精华液单| 日韩强制内射视频| 观看免费一级毛片| eeuss影院久久| 日韩免费高清中文字幕av| 国产精品不卡视频一区二区| 国产欧美另类精品又又久久亚洲欧美| 欧美 日韩 精品 国产| 精品酒店卫生间| 草草在线视频免费看| 日韩电影二区| .国产精品久久| 国内少妇人妻偷人精品xxx网站| 亚洲精品色激情综合| 2021少妇久久久久久久久久久| 国产成人福利小说| 色播亚洲综合网| 国产精品蜜桃在线观看| 亚洲性久久影院| 日韩欧美精品v在线| 午夜免费鲁丝| 18禁动态无遮挡网站| 精品少妇黑人巨大在线播放| 亚洲欧美精品专区久久| 观看美女的网站| 婷婷色av中文字幕| 久久精品国产亚洲av涩爱| 精品人妻视频免费看| 校园人妻丝袜中文字幕| 如何舔出高潮| 黄色一级大片看看| 日韩伦理黄色片| 一个人观看的视频www高清免费观看| av线在线观看网站| av福利片在线观看| 亚洲av.av天堂| 日本一二三区视频观看| 国产成人a区在线观看| 人妻夜夜爽99麻豆av| 亚洲欧洲日产国产| 最近2019中文字幕mv第一页| kizo精华| 欧美日韩精品成人综合77777| 国产精品女同一区二区软件| 欧美xxxx性猛交bbbb| 成人欧美大片| 纵有疾风起免费观看全集完整版| 亚州av有码| 国产在视频线精品| 少妇熟女欧美另类| 国产在线一区二区三区精| 中文资源天堂在线| 国产一区二区三区综合在线观看 | 久久久久国产精品人妻一区二区| 九九在线视频观看精品| 卡戴珊不雅视频在线播放| 可以在线观看毛片的网站| 又爽又黄a免费视频| 国产精品福利在线免费观看| 亚洲欧美一区二区三区国产| 在线看a的网站| 免费观看无遮挡的男女| 久久久久久久久久成人| 久久久精品欧美日韩精品| av在线老鸭窝| 99热网站在线观看| 欧美高清性xxxxhd video| 日韩av不卡免费在线播放| 日韩欧美精品免费久久| 大香蕉久久网| 听说在线观看完整版免费高清| videos熟女内射| 91aial.com中文字幕在线观看| 一本一本综合久久| 国产老妇女一区| 视频中文字幕在线观看| 日韩视频在线欧美| 在线观看av片永久免费下载| 噜噜噜噜噜久久久久久91| 亚洲欧美中文字幕日韩二区| 精品视频人人做人人爽| 男的添女的下面高潮视频| 久久久久久久精品精品| 熟女av电影| 久久久国产一区二区| 久久精品国产鲁丝片午夜精品| 国产免费一级a男人的天堂| 亚洲国产欧美人成| 99久久精品国产国产毛片| 亚洲精品一区蜜桃| 天天一区二区日本电影三级| 只有这里有精品99| 在线观看国产h片| 国产中年淑女户外野战色| 色5月婷婷丁香| 国产精品无大码| 91久久精品国产一区二区成人| 久久久久久久久久久丰满| 夜夜看夜夜爽夜夜摸| 久久人人爽人人爽人人片va| 99热6这里只有精品| 亚洲性久久影院| 国产伦精品一区二区三区视频9| 2021少妇久久久久久久久久久| 国产午夜精品一二区理论片| 亚洲欧美日韩无卡精品| 丰满少妇做爰视频| 80岁老熟妇乱子伦牲交| 男女下面进入的视频免费午夜| 日本欧美国产在线视频| 免费av观看视频| 又粗又硬又长又爽又黄的视频| 蜜桃亚洲精品一区二区三区| 黄片无遮挡物在线观看| 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 美女cb高潮喷水在线观看| 99久久九九国产精品国产免费| 超碰av人人做人人爽久久| 波多野结衣巨乳人妻| 亚洲自偷自拍三级| 国产成人精品一,二区| 少妇人妻久久综合中文| 中文字幕免费在线视频6| 国产精品福利在线免费观看| 韩国高清视频一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲人与动物交配视频| 寂寞人妻少妇视频99o| av国产久精品久网站免费入址| 久久久久国产网址| 啦啦啦中文免费视频观看日本| 尾随美女入室| 一级二级三级毛片免费看| 午夜视频国产福利| 国产精品三级大全| 亚洲丝袜综合中文字幕| 精品久久久久久久久av| 人妻 亚洲 视频| 各种免费的搞黄视频| 丝瓜视频免费看黄片| av国产免费在线观看| 天堂俺去俺来也www色官网| 人妻夜夜爽99麻豆av| 日本一二三区视频观看| 国产熟女欧美一区二区| 毛片一级片免费看久久久久| 国产成人精品福利久久| 啦啦啦中文免费视频观看日本| 国产成人一区二区在线| 一级片'在线观看视频| 又黄又爽又刺激的免费视频.| 久久久成人免费电影| 亚洲av不卡在线观看| 亚洲天堂av无毛| 国产视频内射| 日韩欧美精品v在线| 国产精品一二三区在线看| 在线观看av片永久免费下载| 久久久久精品性色| 国产av码专区亚洲av| 久久6这里有精品| 亚洲精华国产精华液的使用体验| 久久久久久国产a免费观看| 男人添女人高潮全过程视频| 日本猛色少妇xxxxx猛交久久| 国产免费又黄又爽又色| 色哟哟·www| 精品视频人人做人人爽| 久久久精品94久久精品| 午夜视频国产福利| 啦啦啦在线观看免费高清www| 国产大屁股一区二区在线视频| 国产一区二区在线观看日韩| 国产久久久一区二区三区| 视频区图区小说| 成年av动漫网址| 成人国产麻豆网| 蜜桃久久精品国产亚洲av| 欧美日韩精品成人综合77777| 久热久热在线精品观看| av免费观看日本| 极品教师在线视频| 久久久久精品久久久久真实原创| 久久精品国产鲁丝片午夜精品| 搡女人真爽免费视频火全软件| 蜜桃久久精品国产亚洲av| 成人国产麻豆网| 丝瓜视频免费看黄片| 国产伦精品一区二区三区视频9| 欧美成人午夜免费资源| 日韩大片免费观看网站| 国产免费又黄又爽又色| 久久久久久国产a免费观看| 午夜视频国产福利| 深爱激情五月婷婷| 黄色配什么色好看| 一区二区三区免费毛片| 你懂的网址亚洲精品在线观看| 久久久久九九精品影院| 制服丝袜香蕉在线| 亚洲国产精品成人久久小说| 国产黄色免费在线视频| 又爽又黄无遮挡网站| 深夜a级毛片| 国产v大片淫在线免费观看| 日韩中字成人| 国产在线一区二区三区精| 亚洲av男天堂| 国产日韩欧美在线精品| 久久久久久国产a免费观看| 午夜免费男女啪啪视频观看| 久久人人爽av亚洲精品天堂 | 神马国产精品三级电影在线观看| 一级黄片播放器| 汤姆久久久久久久影院中文字幕| 久热久热在线精品观看| 精品一区二区免费观看| 午夜精品一区二区三区免费看| 九色成人免费人妻av| 久久久精品免费免费高清| 国产亚洲午夜精品一区二区久久 | 国产成人免费无遮挡视频| 国产成人精品福利久久| 一本久久精品| 国产高潮美女av| 黄片无遮挡物在线观看| 看非洲黑人一级黄片| 在现免费观看毛片| 夜夜爽夜夜爽视频| 男人添女人高潮全过程视频| 亚洲av中文av极速乱| 久久97久久精品| 国产av国产精品国产| 国产精品一二三区在线看| 美女脱内裤让男人舔精品视频| 舔av片在线| 蜜桃久久精品国产亚洲av| 在线观看av片永久免费下载| 日本av手机在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 婷婷色麻豆天堂久久| 大片免费播放器 马上看| 人人妻人人爽人人添夜夜欢视频 | 亚洲,一卡二卡三卡| 波多野结衣巨乳人妻| 在线观看三级黄色| 国产精品伦人一区二区| 国产黄色免费在线视频| 99久国产av精品国产电影| 黄色日韩在线| 嫩草影院新地址| 免费av观看视频| 精品久久久久久久久av| 国产色爽女视频免费观看| 草草在线视频免费看| 晚上一个人看的免费电影| 亚洲精品,欧美精品| 久久久久国产精品人妻一区二区| 寂寞人妻少妇视频99o| 精品久久久精品久久久| 五月开心婷婷网| 久久99热这里只频精品6学生| 免费看不卡的av| 国产国拍精品亚洲av在线观看| av线在线观看网站| 亚洲激情五月婷婷啪啪| 国产av国产精品国产| 可以在线观看毛片的网站| 99久久精品热视频| 亚洲精品第二区| av福利片在线观看| 视频中文字幕在线观看| 久久久久久久久久成人| 亚洲高清免费不卡视频| av在线蜜桃| 麻豆精品久久久久久蜜桃| eeuss影院久久| 亚洲天堂国产精品一区在线| 亚洲欧洲日产国产| 国产一区二区三区综合在线观看 | 免费黄频网站在线观看国产| 久久久久久国产a免费观看| 高清视频免费观看一区二区| 久久午夜福利片| 成人毛片60女人毛片免费| 黄色视频在线播放观看不卡| 欧美精品一区二区大全| 亚洲欧美中文字幕日韩二区| 日本免费在线观看一区| 少妇猛男粗大的猛烈进出视频 | av在线天堂中文字幕| 中文字幕制服av| 欧美日韩视频精品一区| 超碰97精品在线观看| 2021天堂中文幕一二区在线观| 国产精品99久久99久久久不卡 | 日韩一区二区视频免费看| 少妇被粗大猛烈的视频| 最近中文字幕高清免费大全6| a级一级毛片免费在线观看| 黑人高潮一二区| av在线app专区| 久久久成人免费电影| 狠狠精品人妻久久久久久综合| 午夜日本视频在线| 夫妻午夜视频| 高清毛片免费看| 一区二区三区精品91| 国产在视频线精品| 午夜福利在线观看免费完整高清在| 欧美xxxx性猛交bbbb| 久久精品国产自在天天线| 在线观看一区二区三区激情| 一级二级三级毛片免费看| 波野结衣二区三区在线| 亚洲成人久久爱视频| 日本色播在线视频| 国产精品麻豆人妻色哟哟久久| 嫩草影院精品99| 亚洲成色77777| 国产亚洲av嫩草精品影院| 色网站视频免费| 一区二区三区免费毛片| 伊人久久精品亚洲午夜| 久久久国产一区二区| 欧美高清性xxxxhd video| 亚洲精品日韩av片在线观看| 丝袜脚勾引网站| 赤兔流量卡办理| 国产精品三级大全| 国产亚洲av片在线观看秒播厂| 亚洲国产精品专区欧美| 久久久久久久大尺度免费视频| 日本午夜av视频| 国产精品蜜桃在线观看| 久久久久精品久久久久真实原创| 成年av动漫网址| 99热网站在线观看| 国产精品蜜桃在线观看| 久久久久久九九精品二区国产| 亚洲国产精品成人综合色| 国产精品国产av在线观看| 在线观看三级黄色| 日本三级黄在线观看| 日韩不卡一区二区三区视频在线| 亚洲成人一二三区av| 亚洲精品久久久久久婷婷小说| 免费在线观看成人毛片| 国产免费一级a男人的天堂| 国产精品久久久久久精品电影小说 | 亚洲欧美一区二区三区黑人 | 久久精品国产a三级三级三级| 全区人妻精品视频| 欧美成人a在线观看| 日韩国内少妇激情av| 天天躁日日操中文字幕| 午夜日本视频在线| 精品亚洲乱码少妇综合久久| 国产高清国产精品国产三级 | 成人国产av品久久久| 男女下面进入的视频免费午夜| 国产精品一二三区在线看| 免费高清在线观看视频在线观看| 久久久久性生活片| 国产 一区精品| 免费黄网站久久成人精品| 久久精品国产自在天天线| 久久久久久久久久久丰满| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久鲁丝午夜福利片| 久久精品久久精品一区二区三区| 欧美高清性xxxxhd video| 日韩国内少妇激情av| 舔av片在线| 夫妻性生交免费视频一级片| videossex国产| 91久久精品电影网| 啦啦啦啦在线视频资源| 精品国产乱码久久久久久小说| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| 欧美老熟妇乱子伦牲交| 亚洲精品影视一区二区三区av| 国产成人精品久久久久久| 最近的中文字幕免费完整| 97人妻精品一区二区三区麻豆| 国产精品蜜桃在线观看| 高清日韩中文字幕在线| 天堂网av新在线| 国产高清国产精品国产三级 | 色5月婷婷丁香| 好男人视频免费观看在线| 久久热精品热| 国产乱来视频区| 亚洲国产精品专区欧美| 免费人成在线观看视频色| 午夜亚洲福利在线播放| 黄色配什么色好看| 一二三四中文在线观看免费高清| 一级毛片黄色毛片免费观看视频| 精品一区二区三区视频在线| 少妇人妻一区二区三区视频| 真实男女啪啪啪动态图| 在线观看国产h片| 不卡视频在线观看欧美| 最后的刺客免费高清国语| 22中文网久久字幕| 超碰97精品在线观看| 色吧在线观看| 午夜免费男女啪啪视频观看| 国产高清有码在线观看视频| 制服丝袜香蕉在线| 一级av片app| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 制服丝袜香蕉在线| 亚洲欧美日韩另类电影网站 | av在线观看视频网站免费| 超碰97精品在线观看| 亚洲综合色惰| 国产成人午夜福利电影在线观看| 青青草视频在线视频观看| 黄色日韩在线| 熟女人妻精品中文字幕| 精品国产乱码久久久久久小说| 男人添女人高潮全过程视频| 午夜亚洲福利在线播放| 熟女电影av网| 国产视频内射| 日本wwww免费看| 国产色婷婷99| 国产精品女同一区二区软件| 美女脱内裤让男人舔精品视频|