• <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.

    亚洲欧美成人精品一区二区| 99在线视频只有这里精品首页| 欧美性猛交╳xxx乱大交人| 成人三级黄色视频| 欧美日本视频| 亚洲av成人精品一区久久| 中文字幕免费在线视频6| 性欧美人与动物交配| 美女cb高潮喷水在线观看| 亚洲av熟女| 91久久精品国产一区二区三区| 在线观看一区二区三区| 国产 一区精品| 久久国产乱子免费精品| 男人狂女人下面高潮的视频| 久久精品国产99精品国产亚洲性色| 国产探花在线观看一区二区| 国产毛片a区久久久久| 身体一侧抽搐| 在线观看一区二区三区| 中文字幕熟女人妻在线| 禁无遮挡网站| 久久草成人影院| 又黄又爽又刺激的免费视频.| 岛国毛片在线播放| 亚洲精品久久久久久婷婷小说 | 国产成人aa在线观看| 欧美极品一区二区三区四区| 欧美最黄视频在线播放免费| 国内精品一区二区在线观看| 天天躁夜夜躁狠狠久久av| 国产不卡一卡二| 丰满人妻一区二区三区视频av| 亚洲国产欧美人成| 亚洲婷婷狠狠爱综合网| 亚洲成a人片在线一区二区| 99热这里只有精品一区| 成人一区二区视频在线观看| 在线观看66精品国产| av天堂在线播放| 日日撸夜夜添| 美女xxoo啪啪120秒动态图| 性插视频无遮挡在线免费观看| 非洲黑人性xxxx精品又粗又长| 麻豆乱淫一区二区| 老女人水多毛片| 内射极品少妇av片p| 一级黄色大片毛片| 欧美三级亚洲精品| 亚洲成人久久性| 99国产极品粉嫩在线观看| eeuss影院久久| 国产高清视频在线观看网站| 日本爱情动作片www.在线观看| 日本一二三区视频观看| 99热这里只有是精品在线观看| 九九爱精品视频在线观看| 中文字幕久久专区| 免费看光身美女| 在线天堂最新版资源| 日韩人妻高清精品专区| 18禁在线播放成人免费| 国产91av在线免费观看| 国产午夜福利久久久久久| 国产精品国产三级国产av玫瑰| 18禁裸乳无遮挡免费网站照片| 高清在线视频一区二区三区 | 高清午夜精品一区二区三区 | 在线播放无遮挡| 亚洲精品日韩在线中文字幕 | 亚洲色图av天堂| 久久精品久久久久久噜噜老黄 | 亚洲精品乱码久久久v下载方式| 国产精品一区二区三区四区免费观看| 亚洲综合色惰| 久久久国产成人免费| 成人一区二区视频在线观看| 久久人人精品亚洲av| 免费看美女性在线毛片视频| 免费黄网站久久成人精品| 免费无遮挡裸体视频| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 99久久人妻综合| 国产高清不卡午夜福利| 欧美精品国产亚洲| 久99久视频精品免费| 国产伦精品一区二区三区四那| 色哟哟哟哟哟哟| 村上凉子中文字幕在线| 久久中文看片网| 一区二区三区四区激情视频 | 久久中文看片网| 亚洲最大成人手机在线| 成人美女网站在线观看视频| 久久欧美精品欧美久久欧美| 日韩av在线大香蕉| av在线老鸭窝| 日韩精品青青久久久久久| 噜噜噜噜噜久久久久久91| 色视频www国产| 国产伦一二天堂av在线观看| 国产伦精品一区二区三区视频9| 中国美女看黄片| 国产成人a区在线观看| 1000部很黄的大片| 久久久精品欧美日韩精品| 婷婷色av中文字幕| 国语自产精品视频在线第100页| 美女国产视频在线观看| 国产精品一区www在线观看| 国产精品人妻久久久久久| 亚洲乱码一区二区免费版| 日本免费a在线| 性欧美人与动物交配| 伦理电影大哥的女人| 国产伦精品一区二区三区视频9| 亚洲精品国产av成人精品| 欧美一级a爱片免费观看看| 1000部很黄的大片| а√天堂www在线а√下载| 免费黄网站久久成人精品| 亚洲中文字幕日韩| 日韩强制内射视频| 色尼玛亚洲综合影院| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久亚洲av鲁大| 日日摸夜夜添夜夜爱| 在线播放国产精品三级| 久久久久久九九精品二区国产| 噜噜噜噜噜久久久久久91| 男人狂女人下面高潮的视频| 成年女人看的毛片在线观看| 91精品一卡2卡3卡4卡| 午夜福利在线在线| av在线播放精品| 日韩欧美国产在线观看| 久久精品人妻少妇| 国产成人a∨麻豆精品| 成年av动漫网址| 国产亚洲精品久久久久久毛片| 天堂√8在线中文| 精品无人区乱码1区二区| 一个人看视频在线观看www免费| 日韩制服骚丝袜av| 国产美女午夜福利| 一个人看的www免费观看视频| av又黄又爽大尺度在线免费看 | 中出人妻视频一区二区| 精品熟女少妇av免费看| 18禁在线播放成人免费| 男人狂女人下面高潮的视频| 亚洲性久久影院| 久久99蜜桃精品久久| 天天躁日日操中文字幕| 中文字幕av在线有码专区| 亚洲精品成人久久久久久| 国产成人影院久久av| 国产激情偷乱视频一区二区| 热99在线观看视频| 三级毛片av免费| 天天一区二区日本电影三级| 免费电影在线观看免费观看| 99热6这里只有精品| 国产精品乱码一区二三区的特点| 久久人人爽人人爽人人片va| 亚洲av电影不卡..在线观看| 22中文网久久字幕| 能在线免费看毛片的网站| 日韩一区二区视频免费看| 久久国产乱子免费精品| 99在线视频只有这里精品首页| 国产伦一二天堂av在线观看| 亚洲一级一片aⅴ在线观看| 欧美xxxx性猛交bbbb| 亚洲欧美日韩东京热| av免费在线看不卡| 成人鲁丝片一二三区免费| 女同久久另类99精品国产91| 久久综合国产亚洲精品| 免费人成在线观看视频色| 18禁黄网站禁片免费观看直播| 日本一二三区视频观看| 亚洲国产精品合色在线| 综合色av麻豆| 观看免费一级毛片| 亚洲精品成人久久久久久| 国内精品宾馆在线| 婷婷亚洲欧美| 听说在线观看完整版免费高清| 国产91av在线免费观看| 日本五十路高清| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件| 精品一区二区三区视频在线| 在线观看66精品国产| 欧美性猛交黑人性爽| 一卡2卡三卡四卡精品乱码亚洲| 我的老师免费观看完整版| 综合色丁香网| 免费一级毛片在线播放高清视频| 18禁黄网站禁片免费观看直播| 久久九九热精品免费| 尾随美女入室| 国产成人精品婷婷| 蜜桃久久精品国产亚洲av| 在线观看av片永久免费下载| 日本一二三区视频观看| 亚洲,欧美,日韩| 亚洲欧美日韩无卡精品| 免费观看人在逋| 久久精品夜色国产| 国产精品野战在线观看| 亚洲成人久久性| 国产精品人妻久久久影院| 一本久久中文字幕| 干丝袜人妻中文字幕| 久久久久国产网址| 男女做爰动态图高潮gif福利片| 久久午夜福利片| 免费人成在线观看视频色| av在线观看视频网站免费| 日本爱情动作片www.在线观看| 日韩欧美国产在线观看| 日本-黄色视频高清免费观看| 九九久久精品国产亚洲av麻豆| 成人午夜高清在线视频| 日本黄色片子视频| 国产亚洲精品av在线| 悠悠久久av| 亚洲精品国产成人久久av| 校园春色视频在线观看| 精品熟女少妇av免费看| 激情 狠狠 欧美| 99热这里只有精品一区| 不卡一级毛片| 听说在线观看完整版免费高清| 少妇熟女欧美另类| 国产精品不卡视频一区二区| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 内地一区二区视频在线| 欧美最新免费一区二区三区| 99久久精品热视频| 美女内射精品一级片tv| 如何舔出高潮| 岛国毛片在线播放| 3wmmmm亚洲av在线观看| 美女内射精品一级片tv| 黄色欧美视频在线观看| 亚洲欧美精品自产自拍| 丰满乱子伦码专区| 人妻少妇偷人精品九色| 尤物成人国产欧美一区二区三区| av天堂中文字幕网| 日韩中字成人| 亚洲五月天丁香| 高清在线视频一区二区三区 | 日韩亚洲欧美综合| 亚洲国产精品久久男人天堂| 中文在线观看免费www的网站| 欧美最黄视频在线播放免费| 国产三级中文精品| 欧美日韩乱码在线| 啦啦啦韩国在线观看视频| 国产伦精品一区二区三区视频9| 欧美色欧美亚洲另类二区| 日本-黄色视频高清免费观看| 日日摸夜夜添夜夜爱| 国产精品人妻久久久久久| 日日撸夜夜添| 久久久久久久久久久丰满| 中文欧美无线码| 亚洲国产精品成人综合色| 精品久久久久久成人av| 欧美人与善性xxx| 国产中年淑女户外野战色| 赤兔流量卡办理| 国产女主播在线喷水免费视频网站 | 国产亚洲91精品色在线| 日韩欧美精品免费久久| 国产久久久一区二区三区| 精品熟女少妇av免费看| 99久久中文字幕三级久久日本| 少妇裸体淫交视频免费看高清| 乱系列少妇在线播放| 老师上课跳d突然被开到最大视频| 久久久久免费精品人妻一区二区| 亚洲欧美精品自产自拍| 婷婷六月久久综合丁香| 91在线精品国自产拍蜜月| 免费在线观看成人毛片| 国产综合懂色| 熟女电影av网| 精品午夜福利在线看| 亚洲成a人片在线一区二区| 小说图片视频综合网站| 国国产精品蜜臀av免费| 人妻久久中文字幕网| 网址你懂的国产日韩在线| 亚洲自偷自拍三级| 日本黄色片子视频| 日本熟妇午夜| 国产精品免费一区二区三区在线| 老司机影院成人| 人妻制服诱惑在线中文字幕| 亚洲国产欧洲综合997久久,| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 深夜精品福利| 一进一出抽搐动态| 中文字幕av在线有码专区| 只有这里有精品99| 久久午夜福利片| 可以在线观看毛片的网站| 国产 一区精品| 精品一区二区三区视频在线| 日韩,欧美,国产一区二区三区 | 97人妻精品一区二区三区麻豆| 卡戴珊不雅视频在线播放| 国内精品久久久久精免费| 久久99蜜桃精品久久| 春色校园在线视频观看| 成人午夜精彩视频在线观看| 国内揄拍国产精品人妻在线| 国产精品美女特级片免费视频播放器| 一级毛片aaaaaa免费看小| 在线天堂最新版资源| 少妇的逼好多水| 能在线免费观看的黄片| 久久亚洲国产成人精品v| 我的老师免费观看完整版| 18禁在线播放成人免费| 男人狂女人下面高潮的视频| 日韩,欧美,国产一区二区三区 | 亚洲欧美日韩无卡精品| 亚洲高清免费不卡视频| 大又大粗又爽又黄少妇毛片口| 最好的美女福利视频网| 日韩高清综合在线| 久久这里有精品视频免费| 久久99蜜桃精品久久| 亚洲aⅴ乱码一区二区在线播放| 欧美成人免费av一区二区三区| 小说图片视频综合网站| 亚州av有码| 国产精品久久久久久精品电影小说 | 午夜视频国产福利| 男的添女的下面高潮视频| 一边亲一边摸免费视频| 国产成人aa在线观看| 18+在线观看网站| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 老司机福利观看| 亚洲经典国产精华液单| 美女黄网站色视频| 1024手机看黄色片| 久久久久久久亚洲中文字幕| 亚洲中文字幕一区二区三区有码在线看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产三级在线视频| 国产高清有码在线观看视频| 九九爱精品视频在线观看| 99视频精品全部免费 在线| 爱豆传媒免费全集在线观看| 男人的好看免费观看在线视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人久久性| 激情 狠狠 欧美| 中文字幕免费在线视频6| 亚洲av二区三区四区| 国产伦精品一区二区三区视频9| 亚洲七黄色美女视频| 特大巨黑吊av在线直播| 三级经典国产精品| 国产精品,欧美在线| 久久九九热精品免费| 可以在线观看毛片的网站| 国产一区二区三区av在线 | 99精品在免费线老司机午夜| 亚洲精品久久久久久婷婷小说 | 久久久午夜欧美精品| 91久久精品国产一区二区三区| 国产精品久久久久久久久免| 亚洲成人av在线免费| 久久久久久久久久黄片| 亚洲欧美中文字幕日韩二区| 波多野结衣巨乳人妻| 婷婷色综合大香蕉| 18禁在线无遮挡免费观看视频| 中文字幕人妻熟人妻熟丝袜美| 欧美一区二区精品小视频在线| 淫秽高清视频在线观看| 久久99热6这里只有精品| 网址你懂的国产日韩在线| 淫秽高清视频在线观看| 嘟嘟电影网在线观看| 亚洲成人av在线免费| 亚洲av成人av| 麻豆乱淫一区二区| 亚洲av免费在线观看| 亚洲国产色片| 欧美日韩国产亚洲二区| 国产蜜桃级精品一区二区三区| 日韩高清综合在线| 亚洲在线自拍视频| 此物有八面人人有两片| 一个人免费在线观看电影| 亚洲天堂国产精品一区在线| 91午夜精品亚洲一区二区三区| 一夜夜www| 国内久久婷婷六月综合欲色啪| 日本黄大片高清| 听说在线观看完整版免费高清| 欧美极品一区二区三区四区| 你懂的网址亚洲精品在线观看 | 国产精品免费一区二区三区在线| 人妻少妇偷人精品九色| 国产精品不卡视频一区二区| 丰满人妻一区二区三区视频av| 亚洲欧美中文字幕日韩二区| 欧美+日韩+精品| 高清日韩中文字幕在线| 国产美女午夜福利| 白带黄色成豆腐渣| 丰满人妻一区二区三区视频av| 插阴视频在线观看视频| 欧美变态另类bdsm刘玥| 精华霜和精华液先用哪个| 国产亚洲av嫩草精品影院| 91精品一卡2卡3卡4卡| 99久久精品热视频| 99在线视频只有这里精品首页| 中文资源天堂在线| 国产精品久久久久久精品电影| 国产精品一区www在线观看| 亚洲欧美精品自产自拍| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久亚洲精品不卡| 日韩 亚洲 欧美在线| 婷婷精品国产亚洲av| 国产一区二区在线av高清观看| 亚洲图色成人| 亚洲精品日韩在线中文字幕 | 国产精品乱码一区二三区的特点| 精品久久国产蜜桃| 国产精品永久免费网站| 国内少妇人妻偷人精品xxx网站| 日日干狠狠操夜夜爽| 国产精品电影一区二区三区| 3wmmmm亚洲av在线观看| 国产真实乱freesex| 亚洲最大成人中文| 久久久久久大精品| 国产大屁股一区二区在线视频| 日本与韩国留学比较| 国产片特级美女逼逼视频| 看免费成人av毛片| 九九热线精品视视频播放| 国产av一区在线观看免费| 国产极品天堂在线| 好男人在线观看高清免费视频| h日本视频在线播放| 嫩草影院精品99| 日本一本二区三区精品| 婷婷色综合大香蕉| 亚洲国产色片| 国产高清不卡午夜福利| 日韩成人av中文字幕在线观看| 波多野结衣巨乳人妻| 国产成人91sexporn| 成人鲁丝片一二三区免费| 一区二区三区高清视频在线| av卡一久久| 国产成人一区二区在线| 久久久精品94久久精品| 色尼玛亚洲综合影院| 亚洲乱码一区二区免费版| 久99久视频精品免费| 网址你懂的国产日韩在线| 国产精品爽爽va在线观看网站| 国产 一区 欧美 日韩| 男女做爰动态图高潮gif福利片| 日韩一区二区三区影片| 亚洲中文字幕一区二区三区有码在线看| 亚洲婷婷狠狠爱综合网| 国产免费男女视频| 国产精品爽爽va在线观看网站| 国产色爽女视频免费观看| www日本黄色视频网| 久久人人爽人人片av| 亚洲中文字幕日韩| 欧美在线一区亚洲| 伦精品一区二区三区| 国产精品三级大全| 久久久精品大字幕| 精品一区二区免费观看| 成人亚洲精品av一区二区| 亚洲真实伦在线观看| 最近2019中文字幕mv第一页| 国产一级毛片七仙女欲春2| 嫩草影院入口| 国产成人freesex在线| 色综合站精品国产| 亚洲人成网站高清观看| 国产精品久久久久久精品电影| 国产单亲对白刺激| 日本五十路高清| 久久久国产成人免费| 日本成人三级电影网站| av卡一久久| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品国产成人久久av| 久久久精品大字幕| 色吧在线观看| 亚洲天堂国产精品一区在线| 婷婷色综合大香蕉| 国产精品伦人一区二区| 欧美又色又爽又黄视频| 麻豆成人午夜福利视频| 亚洲电影在线观看av| 少妇猛男粗大的猛烈进出视频 | 国产伦理片在线播放av一区 | 国产老妇伦熟女老妇高清| 黄色日韩在线| 欧美zozozo另类| 欧美日本视频| 男女边吃奶边做爰视频| 亚洲丝袜综合中文字幕| 97超视频在线观看视频| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| 亚洲精品456在线播放app| 日韩av不卡免费在线播放| 国产亚洲av嫩草精品影院| 波野结衣二区三区在线| 男女视频在线观看网站免费| 两性午夜刺激爽爽歪歪视频在线观看| 黄色视频,在线免费观看| 看黄色毛片网站| 国产国拍精品亚洲av在线观看| 亚洲av成人精品一区久久| 黑人高潮一二区| 国产成人精品久久久久久| 校园人妻丝袜中文字幕| 麻豆成人av视频| 久久久精品欧美日韩精品| 免费av观看视频| 97在线视频观看| 美女被艹到高潮喷水动态| 九色成人免费人妻av| 最近最新中文字幕大全电影3| 精品一区二区免费观看| 午夜精品国产一区二区电影 | 伦理电影大哥的女人| 日韩av在线大香蕉| 久久久久国产网址| 国产黄片视频在线免费观看| 一个人看视频在线观看www免费| 欧美性感艳星| 国产免费男女视频| 欧美日韩乱码在线| 久久久久国产网址| 嫩草影院入口| 亚洲av二区三区四区| 欧美xxxx黑人xx丫x性爽| 人人妻人人澡人人爽人人夜夜 | 免费黄网站久久成人精品| 久久精品国产自在天天线| 91精品国产九色| 色噜噜av男人的天堂激情| 成人亚洲精品av一区二区| 少妇丰满av| 精品久久久久久久久久久久久| 亚洲综合色惰| 亚洲人成网站高清观看| 亚洲精品色激情综合| 亚洲人成网站高清观看| 波多野结衣高清作品| 99精品在免费线老司机午夜| 精品久久国产蜜桃| 国产精品野战在线观看| 亚洲欧美日韩东京热| 欧美bdsm另类| 亚洲精品乱码久久久v下载方式| 丰满的人妻完整版| 亚洲天堂国产精品一区在线| 久久精品国产清高在天天线| 级片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本av手机在线免费观看| 日韩国内少妇激情av| 亚洲欧美日韩卡通动漫| 大香蕉久久网| 好男人视频免费观看在线| 国产熟女欧美一区二区| 给我免费播放毛片高清在线观看| 三级国产精品欧美在线观看| 亚洲人与动物交配视频| 自拍偷自拍亚洲精品老妇| 99久国产av精品| 日韩欧美一区二区三区在线观看| 国产亚洲av片在线观看秒播厂 | 久99久视频精品免费| 日韩一本色道免费dvd| 国产中年淑女户外野战色| 国产精品福利在线免费观看| 女的被弄到高潮叫床怎么办| 精品国产三级普通话版| 一本精品99久久精品77| АⅤ资源中文在线天堂| 看黄色毛片网站| 日本av手机在线免费观看|