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

    微分求積方法解時(shí)間分?jǐn)?shù)階擴(kuò)散方程

    2017-12-02 03:55:49張學(xué)瑩
    關(guān)鍵詞:微分插值步長(zhǎng)

    曹 煥, 張學(xué)瑩, 劉 薈

    (河海大學(xué) 理學(xué)院, 江蘇 南京 211100)

    微分求積方法解時(shí)間分?jǐn)?shù)階擴(kuò)散方程

    曹 煥, 張學(xué)瑩*, 劉 薈

    (河海大學(xué) 理學(xué)院, 江蘇 南京 211100)

    微分求積(DQ)法是一種基于徑向基函數(shù)(RBFs)插值的無網(wǎng)格方法.本文選MultiQ-uadrics(MQ)函數(shù)作為徑向基函數(shù),并采用微分求積方法解決時(shí)間分?jǐn)?shù)階擴(kuò)散方程。在離散過程中,采用有限差分法離散時(shí)間項(xiàng),采用微分求積方法離散空間項(xiàng).最后,結(jié)合數(shù)值求解的結(jié)果做出相應(yīng)的誤差分析.

    徑向基函數(shù); MQ函數(shù); 微分求積法; 時(shí)間分?jǐn)?shù)階擴(kuò)散方程

    0 引言

    近年來,反常擴(kuò)散現(xiàn)象在物理學(xué)和工程學(xué)等領(lǐng)域得到了廣泛的關(guān)注[1,2],比如:污染物的運(yùn)輸、滲流、損耗以及磁等離子體等問題.與正常的擴(kuò)散現(xiàn)象相比,反常擴(kuò)散現(xiàn)象具有較強(qiáng)的遠(yuǎn)程相互作用與歷史依賴性等特點(diǎn),而標(biāo)準(zhǔn)的整數(shù)階偏微分方程已經(jīng)無法準(zhǔn)確的描述這類反常擴(kuò)散行為[3].隨后,分

    數(shù)階導(dǎo)數(shù)被陸續(xù)的證實(shí)是一種可以準(zhǔn)確的描述這類反常擴(kuò)散現(xiàn)象的數(shù)學(xué)模型[4].所以,為了解決這些問題,引入了分?jǐn)?shù)階擴(kuò)散方程.

    由于分?jǐn)?shù)階偏微分方程中含有分?jǐn)?shù)階的導(dǎo)數(shù),不容易求解分?jǐn)?shù)階偏微分方程的解析解,因此很多人研究其數(shù)值解法.目前,有限差分方法[5]已被證實(shí)為對(duì)分?jǐn)?shù)階擴(kuò)散方程進(jìn)行時(shí)間和空間上離散的一種有效方法.對(duì)于分?jǐn)?shù)階擴(kuò)散方程的數(shù)值解法,也可以選擇其他合適的方法對(duì)其空間上的離散,例如:邊界元法[6],有限元法[7,8],傅里葉法[9],無網(wǎng)格徑向基函數(shù)法[10-13]等.與傳統(tǒng)的有限差分方法對(duì)空間的離散相比較,對(duì)于處理一些大規(guī)模問題,這些方法在一定程度上減少了計(jì)算量,同時(shí)又提高了計(jì)算效率.2003年,Shu C提出DQ法解整數(shù)階的偏微分方程,且取得較理想的結(jié)果.本文將這一方法運(yùn)用到解分?jǐn)?shù)階偏微分方程中,并給出誤差結(jié)果分析.

    1 基本理論

    1.1 時(shí)間分?jǐn)?shù)階擴(kuò)散方程

    (1)

    0lt;αlt;1,X∈Ω,t∈(0,T)

    邊界條件:

    u(X,t)=g(X,t),X∈?Ω,t∈(0,T)

    (2)

    初始條件:

    u(X,0)=u0(X),X∈Ω

    (3)

    式(1)~(3)中:Q(X,t),g(X,t),u0(x)是給定的函數(shù),T是總時(shí)間,?α/?tα是關(guān)于時(shí)間t的分?jǐn)?shù)階導(dǎo).

    本文采用Caputo時(shí)間分?jǐn)?shù)階導(dǎo)數(shù)的定義,定義如下:

    (4)

    這里0lt;αlt;1.

    1.2 FDM方法對(duì)時(shí)間的離散

    首先把區(qū)間[0,T]分割成K個(gè)小區(qū)間,定義tk=kτ,k=0,1,…,K,這里τ=T/K為時(shí)間步長(zhǎng).用u(X,tk)表示函數(shù)u(X,t)在tk處的精確解.本文利用有限差分方法在t=tk+1處對(duì)時(shí)間分?jǐn)?shù)階導(dǎo)數(shù)進(jìn)行離散,有

    (5)

    式(5)中:

    根據(jù)式(5),對(duì)方程式(1)進(jìn)行時(shí)間離散,得到

    Δuk+1+Qk+1=

    (6)

    1.3 微分求積方法對(duì)空間的離散

    設(shè)f(x)是光滑函數(shù),其中x=(x1,x2,…,xn)T.已知它在節(jié)點(diǎn)xi及其支撐域內(nèi)各支撐點(diǎn)xj,j=1,2,…,n上的函數(shù)值,則f(x)關(guān)于xk的m階導(dǎo)數(shù)可以表示為

    (7)

    (8)

    可將其寫成矩陣形式是

    (9)

    由文獻(xiàn)[14]可知式(9)中的矩陣是條件正定的,故系數(shù)矩陣是可逆的,從而可求得系數(shù).再將其代入式(7),就可以求出函數(shù)f(x)關(guān)于xk的m階導(dǎo)數(shù)的近似值.

    由上述過程,本文采用DQ方法對(duì)方程式(6)進(jìn)行空間離散,可知

    (10)

    再將式(10)代入式(6)得到對(duì)空間的離散,并對(duì)其整理可得

    (11)

    2 數(shù)值算例

    本文考慮在規(guī)則區(qū)域Ω=[0,0.5]×[0,0.5]內(nèi),檢驗(yàn)時(shí)間分?jǐn)?shù)階擴(kuò)散方程的有效性和準(zhǔn)確性.數(shù)值的準(zhǔn)確性由均方根誤差(RMSE),最大絕對(duì)誤差(MAE)和最大相對(duì)誤差(RAE)計(jì)算.其誤差定義如下:

    算例

    0lt;αlt;1,X∈Ω,t∈T

    邊界條件:

    u(X,t)=t2ex+y,

    X∈?Ω,t∈(0,T)

    初始條件:

    u(X,0)=0,X∈Ω.

    α=0.9,Ω=[0,0.5]2,

    Q(X,t)=(2t2-α/Γ(3-α)-2t2)ex+y.

    該方程的精確解為

    u(X,t)=t2ex+y,x∈Ω,t∈(0,T).

    圖1為選取網(wǎng)格點(diǎn)N=961,分?jǐn)?shù)階導(dǎo)數(shù)α=0.9,時(shí)間步長(zhǎng)dt=0.1,參數(shù)c=0.05時(shí)的絕對(duì)誤差和相對(duì)誤差圖.

    (a) N=961,α=0.9,dt=0.1,c=0.5時(shí)絕對(duì)誤差圖

    (b) N=961,α=0.9,dt=0.1,c=0.5時(shí)相對(duì)誤差分布圖

    圖2為選取分?jǐn)?shù)階導(dǎo)數(shù)α=0.9,參數(shù)c=0.02時(shí),隨著插值點(diǎn)個(gè)數(shù)增加,不同時(shí)間步長(zhǎng)下的均方差對(duì)比圖;圖3為選取分?jǐn)?shù)階導(dǎo)數(shù)α=0.9,時(shí)間步長(zhǎng)t=0.1時(shí),隨著插值點(diǎn)個(gè)數(shù)增加,不同參數(shù)下的均方差對(duì)比圖.

    圖2 不同時(shí)間步長(zhǎng)下的均方根誤差比較

    圖3 不同參數(shù)c下的均方根誤差比較

    從圖2和圖3可以看出,參數(shù)c和時(shí)間步長(zhǎng)的選取對(duì)數(shù)值結(jié)果有較大的影響,當(dāng)選取較大的時(shí)間步長(zhǎng)時(shí),精度有明顯的提高,當(dāng)插值點(diǎn)個(gè)數(shù)為256時(shí),精度達(dá)到最高,當(dāng)插值點(diǎn)個(gè)數(shù)的增加,精度變化幅度逐漸減小,最終趨于穩(wěn)定,但如果點(diǎn)再繼續(xù)增加,就可能會(huì)出現(xiàn)病態(tài);當(dāng)選取較小的參數(shù)c時(shí),隨著插值點(diǎn)個(gè)數(shù)的增加,精度開始有明顯的提高,但當(dāng)插值點(diǎn)個(gè)數(shù)再增加時(shí),精度變化幅度逐漸減小,最終趨于穩(wěn)定.因此,選取適當(dāng)?shù)臅r(shí)間步長(zhǎng)和參數(shù)c是提高數(shù)值結(jié)果精度的有效方法.

    表1為選取時(shí)間步長(zhǎng)t=0.1,參數(shù)c=0.02時(shí),不同分?jǐn)?shù)階導(dǎo)數(shù)下的絕對(duì)誤差和相對(duì)誤差的結(jié)果比較.

    表1 t=0.1,c=0.02時(shí)不同α下的絕對(duì)誤差和相對(duì)誤差比較

    從表1可以看出,當(dāng)取時(shí)間步長(zhǎng)t=0.1,參數(shù)c=0.02,分?jǐn)?shù)階的導(dǎo)數(shù)越接近1時(shí),隨著插值點(diǎn)個(gè)數(shù)的增加,誤差越來越精確.

    3 結(jié)論

    本文在結(jié)合有限差分方法的基礎(chǔ)上,應(yīng)用微分求積方法對(duì)規(guī)則區(qū)域上的時(shí)間分?jǐn)?shù)階擴(kuò)散方程進(jìn)行數(shù)值求解.根據(jù)Caputo分?jǐn)?shù)階導(dǎo)數(shù)的積分定義,利用有限差分方法對(duì)時(shí)間導(dǎo)數(shù)進(jìn)行離散,利用微分求積方法對(duì)空間導(dǎo)數(shù)進(jìn)行離散,并通過數(shù)值算例得出數(shù)值結(jié)果.

    數(shù)值結(jié)果表明,微分求積方法是解決時(shí)間分?jǐn)?shù)階擴(kuò)散方程比較可行的方法.在解決問題的過程中,選取了MQ函數(shù)作為徑向基函數(shù),從圖3可以看出參數(shù)c的選擇對(duì)數(shù)值結(jié)果具有很大的影響,當(dāng)選取適當(dāng)參數(shù)c時(shí),數(shù)值結(jié)果的精度更高.同時(shí),還可得出時(shí)間步長(zhǎng)和插值點(diǎn)個(gè)數(shù)也對(duì)其具有一定影響.

    后續(xù),還可進(jìn)一步討論在不規(guī)則區(qū)域上此方法的可行性以及怎樣選取參數(shù)c可使數(shù)值結(jié)果達(dá)到更好的精度.

    [1] R Gorenflo,F Mainardi,D Noretti,et al.Discrete random walk models for space time fractional diffusion[J].Chenmical Physics,2007,284(1-2):521-541.

    [2] L Baowen,W Jiao.Anomalous heat conduction and anomalous diffusion in onedimensional systems[J].Physical Review Letters,2003,91(4):318-322.

    [3] Zhang X Y,Zhu H Q,Kuo L H.A comparison study of the LMAPS method and the LDQ method for time-dependent problems[J].Engineering Analysis with Boundary Elements,2013,37:1 408-1 415.

    [4] F Liu,V Anh,I Turner,et al.Time fractional advection dispersion equation[J].Journal of Applied Mathematics and Computing,2003,13:233-245.

    [5] H Wang,K Wang,T Sircar.A direct O (Nlog2N) finite difference method for fractional diffusion equations[J].Journal of Computational Physics,2010,229(299):8 095-8 104.

    [6] Katsikadelis J T.The BEM for numerical solution of partial fractional differential equations[J].Computers amp; Mathematics with Applications,2011,62(3):891-901.

    [7] Li C,Zhao Z, Chen Y. Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion[J].Computers amp; Mathematics with Applications,2011,62(3):855-875.

    [8] Jiang Y,Ma J.High-order finite element methods for time-fractional partial differential equations[J].Journal of Computational amp; Applied Mathematics,2011,235(11):3 285-3 290.

    [9] Chen C M,Liu F,Turner I,et al.A Fourier method for the fractional diffusion equation describing sub-diffusion[J].Journal of Computational Physics,2007,227(2):886-897.

    [10] Chen W,Ye L,Sun H.Fractional diffusion equations by the Kansa method[J].Computers amp; Mathematics with Applications,2010,59(5):1 614-1 620.

    [11] Liu Q,Gu Y,Zhang P,et al.An implicit RBF meshless approach for time fractional diffusion equations[J].Computational Mechanics,2011,48(1):1-12.

    [12] Fu Z J,Wen Chen,Leevan Ling,el al.Method of approximate particular solutions for constantand variable-order fractional diffusion modes[J].Engineering Analysis with Boundary Elements,2015,57:37-46.

    [13] Chen W,Fu Z J.Recent advances on radial basis function collocation methods[M].Berlin:Springer Berlin Heidelberg,2014.

    [14] Micchelli C A.Interolation of scattered deta:Distance matrices and conditionally positive definite functions[J].Constr Approx,1986(2):11-22.

    【責(zé)任編輯:蔣亞儒】

    Thedifferentialquadraturemethodsolutionsforthetime-fractionaldiffusionequation

    CAO Huan, ZHANG Xue-ying*, LIU Hui

    (School of Science, Hohai University, Nanjing 211100, China)

    The differential quadrature method is an alternative radial basis functions (RBFs) meshless method.This article selects MQ function as the radial basis function,and applys the differential quadrature method to solve the time fractional diffusion equation.In the discretization formulation,a finite difference scheme and the DQ are used repectively to discretize time fractional derivative and spatial derivative terms.Finally we make the error analysis with the results of the numerical investigation example.

    RBFs; MQ function; DQ method; time fractional diffusion function

    2017-08-27

    教育部留學(xué)回國(guó)人員科研啟動(dòng)基金項(xiàng)目(20145003412); 江蘇省自然科學(xué)基金項(xiàng)目(BK20160853)

    曹 煥(1992-),女,山東菏澤人,在讀碩士研究生,研究方向:計(jì)算數(shù)學(xué)

    張學(xué)瑩(1973-),男,山東郯城人,副教授,碩士生導(dǎo)師,研究方向:計(jì)算數(shù)學(xué), zhangxy@hhu.edu.cn

    2096-398X(2017)06-0179-04

    O241.82

    A

    猜你喜歡
    微分插值步長(zhǎng)
    基于Armijo搜索步長(zhǎng)的BFGS與DFP擬牛頓法的比較研究
    擬微分算子在Hp(ω)上的有界性
    上下解反向的脈沖微分包含解的存在性
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    借助微分探求連續(xù)函數(shù)的極值點(diǎn)
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    對(duì)不定積分湊微分解法的再認(rèn)識(shí)
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥搜索算法
    Blackman-Harris窗的插值FFT諧波分析與應(yīng)用
    国产真人三级小视频在线观看| a在线观看视频网站| 两人在一起打扑克的视频| 中文字幕久久专区| 熟女电影av网| 男女做爰动态图高潮gif福利片| 国产精品自产拍在线观看55亚洲| 国产日本99.免费观看| 亚洲 欧美一区二区三区| 小说图片视频综合网站| 一a级毛片在线观看| videosex国产| 国产成人aa在线观看| 精品国产乱子伦一区二区三区| av福利片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久久久久九九精品影院| 最新美女视频免费是黄的| 免费无遮挡裸体视频| 少妇人妻一区二区三区视频| 国产高清有码在线观看视频 | 精品国产亚洲在线| 99热只有精品国产| 最近在线观看免费完整版| 正在播放国产对白刺激| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线| 男女视频在线观看网站免费 | 88av欧美| 日本一二三区视频观看| 亚洲国产精品久久男人天堂| 精品不卡国产一区二区三区| 精品乱码久久久久久99久播| 欧美人与性动交α欧美精品济南到| 欧美久久黑人一区二区| 999精品在线视频| 在线免费观看的www视频| 久久久久国内视频| 亚洲第一电影网av| 黑人巨大精品欧美一区二区mp4| 婷婷六月久久综合丁香| av视频在线观看入口| 欧美色视频一区免费| 亚洲av成人精品一区久久| 欧美精品亚洲一区二区| 国产精品免费一区二区三区在线| 在线观看午夜福利视频| 观看免费一级毛片| 精品电影一区二区在线| 91成年电影在线观看| 狂野欧美白嫩少妇大欣赏| 免费观看精品视频网站| 午夜福利在线在线| 两个人视频免费观看高清| 免费看日本二区| 国产免费男女视频| 亚洲精品久久成人aⅴ小说| 母亲3免费完整高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 不卡av一区二区三区| 成人三级黄色视频| 狠狠狠狠99中文字幕| 又爽又黄无遮挡网站| 在线观看舔阴道视频| 国产三级在线视频| 人成视频在线观看免费观看| 中文在线观看免费www的网站 | 国产精品美女特级片免费视频播放器 | 老司机午夜十八禁免费视频| 9191精品国产免费久久| 亚洲成人精品中文字幕电影| 69av精品久久久久久| 午夜福利成人在线免费观看| 熟女电影av网| 女生性感内裤真人,穿戴方法视频| 久99久视频精品免费| 亚洲精品av麻豆狂野| 日韩欧美 国产精品| 亚洲熟妇熟女久久| 亚洲免费av在线视频| 韩国av一区二区三区四区| 伊人久久大香线蕉亚洲五| 琪琪午夜伦伦电影理论片6080| 国产成人aa在线观看| 9191精品国产免费久久| 国产三级中文精品| 黄片小视频在线播放| 久久精品成人免费网站| 亚洲七黄色美女视频| 久99久视频精品免费| 人人妻人人看人人澡| 国产三级中文精品| 18禁美女被吸乳视频| 可以免费在线观看a视频的电影网站| 亚洲自偷自拍图片 自拍| av国产免费在线观看| 99国产精品一区二区蜜桃av| 国内久久婷婷六月综合欲色啪| 在线a可以看的网站| 男女做爰动态图高潮gif福利片| 国产精品久久视频播放| 一本精品99久久精品77| 日韩欧美 国产精品| 变态另类丝袜制服| www国产在线视频色| 亚洲最大成人中文| 丰满人妻一区二区三区视频av | 美女午夜性视频免费| 免费看美女性在线毛片视频| 亚洲国产中文字幕在线视频| 少妇粗大呻吟视频| 日韩欧美免费精品| 熟女少妇亚洲综合色aaa.| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩一区二区三| 给我免费播放毛片高清在线观看| 色综合婷婷激情| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 欧美另类亚洲清纯唯美| 老司机靠b影院| 亚洲av日韩精品久久久久久密| 亚洲免费av在线视频| 色老头精品视频在线观看| 黄色女人牲交| 两性夫妻黄色片| 成人欧美大片| 五月伊人婷婷丁香| 欧美成狂野欧美在线观看| 国产1区2区3区精品| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 成人18禁在线播放| www.自偷自拍.com| 国内毛片毛片毛片毛片毛片| 成人18禁高潮啪啪吃奶动态图| 日本黄大片高清| 久久久久久大精品| www.熟女人妻精品国产| 神马国产精品三级电影在线观看 | 全区人妻精品视频| 丝袜美腿诱惑在线| 午夜成年电影在线免费观看| 黄色 视频免费看| 精品久久久久久久毛片微露脸| 国产成人欧美在线观看| 久热爱精品视频在线9| 国产欧美日韩精品亚洲av| 很黄的视频免费| 欧美 亚洲 国产 日韩一| 国产主播在线观看一区二区| 一级作爱视频免费观看| 欧美激情久久久久久爽电影| 久久中文看片网| 99精品久久久久人妻精品| 一级毛片高清免费大全| 久久精品aⅴ一区二区三区四区| 18禁裸乳无遮挡免费网站照片| 国产99久久九九免费精品| 亚洲精品一区av在线观看| 91av网站免费观看| 高潮久久久久久久久久久不卡| 欧美+亚洲+日韩+国产| 狠狠狠狠99中文字幕| 观看免费一级毛片| 99国产极品粉嫩在线观看| 麻豆成人午夜福利视频| 一区二区三区国产精品乱码| 香蕉av资源在线| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品中文字幕看吧| 九色成人免费人妻av| 首页视频小说图片口味搜索| 日本黄大片高清| 免费观看人在逋| 正在播放国产对白刺激| 亚洲成人精品中文字幕电影| 99在线视频只有这里精品首页| 亚洲午夜理论影院| 亚洲自拍偷在线| 亚洲精品久久国产高清桃花| 国产成人精品久久二区二区免费| 香蕉av资源在线| 精品无人区乱码1区二区| 亚洲av五月六月丁香网| 亚洲无线在线观看| 国产aⅴ精品一区二区三区波| 少妇粗大呻吟视频| 欧美成人一区二区免费高清观看 | 国内少妇人妻偷人精品xxx网站 | 国产又黄又爽又无遮挡在线| 久久国产乱子伦精品免费另类| 免费看十八禁软件| 男男h啪啪无遮挡| 亚洲片人在线观看| x7x7x7水蜜桃| 欧美三级亚洲精品| 亚洲欧美一区二区三区黑人| 99热这里只有是精品50| 亚洲国产精品sss在线观看| 高清毛片免费观看视频网站| 欧美绝顶高潮抽搐喷水| 99在线人妻在线中文字幕| 欧美3d第一页| 我的老师免费观看完整版| 免费观看人在逋| 国产成人av激情在线播放| 久久久国产欧美日韩av| 久久香蕉激情| 一级黄色大片毛片| 制服丝袜大香蕉在线| 久久热在线av| 亚洲成人精品中文字幕电影| 精品熟女少妇八av免费久了| 麻豆成人av在线观看| 亚洲无线在线观看| 国内精品久久久久精免费| 日韩成人在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 超碰成人久久| 黄色视频不卡| 最新美女视频免费是黄的| 欧美日本亚洲视频在线播放| 亚洲欧美日韩高清专用| 黄色成人免费大全| 欧美日韩乱码在线| 99国产精品99久久久久| 亚洲国产欧洲综合997久久,| 亚洲国产中文字幕在线视频| 村上凉子中文字幕在线| 男插女下体视频免费在线播放| 久久精品国产综合久久久| 色综合亚洲欧美另类图片| 久久这里只有精品19| 三级男女做爰猛烈吃奶摸视频| 亚洲中文av在线| 久久伊人香网站| 久久久久久人人人人人| 国产精品久久视频播放| 免费在线观看视频国产中文字幕亚洲| 哪里可以看免费的av片| 观看免费一级毛片| 男女做爰动态图高潮gif福利片| 久久久精品国产亚洲av高清涩受| 人妻夜夜爽99麻豆av| 亚洲欧美一区二区三区黑人| 久久精品国产综合久久久| 亚洲精品色激情综合| 亚洲第一欧美日韩一区二区三区| 久久久久久九九精品二区国产 | 不卡一级毛片| 亚洲av第一区精品v没综合| 日韩大尺度精品在线看网址| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类 | 欧美黄色片欧美黄色片| 校园春色视频在线观看| 精品久久久久久,| 久久久国产成人免费| 757午夜福利合集在线观看| 美女高潮喷水抽搐中文字幕| 99在线视频只有这里精品首页| 久久久久久亚洲精品国产蜜桃av| 欧美zozozo另类| 国产成人欧美在线观看| 国产亚洲精品av在线| 午夜福利高清视频| 天天添夜夜摸| 欧美色欧美亚洲另类二区| 看黄色毛片网站| 午夜福利18| 亚洲九九香蕉| 午夜激情福利司机影院| 天堂av国产一区二区熟女人妻 | 在线免费观看的www视频| 亚洲av熟女| 亚洲第一欧美日韩一区二区三区| 国产激情偷乱视频一区二区| 少妇粗大呻吟视频| 亚洲国产欧洲综合997久久,| 国产精华一区二区三区| 国产成人啪精品午夜网站| 亚洲色图av天堂| 免费电影在线观看免费观看| 免费在线观看成人毛片| 精品电影一区二区在线| 别揉我奶头~嗯~啊~动态视频| 日本成人三级电影网站| 久久精品91无色码中文字幕| 国产乱人伦免费视频| 制服丝袜大香蕉在线| 麻豆一二三区av精品| 99热只有精品国产| 国内久久婷婷六月综合欲色啪| 美女大奶头视频| 美女 人体艺术 gogo| 久久精品91蜜桃| 日韩av在线大香蕉| 激情在线观看视频在线高清| 啦啦啦韩国在线观看视频| 一级a爱片免费观看的视频| 非洲黑人性xxxx精品又粗又长| 老汉色∧v一级毛片| 麻豆av在线久日| 成人国产综合亚洲| 亚洲aⅴ乱码一区二区在线播放 | 久久性视频一级片| 在线a可以看的网站| 大型av网站在线播放| 国产亚洲精品av在线| 亚洲精品国产精品久久久不卡| 制服人妻中文乱码| 他把我摸到了高潮在线观看| 国产片内射在线| 午夜福利欧美成人| 成人亚洲精品av一区二区| 村上凉子中文字幕在线| 成人国语在线视频| 少妇人妻一区二区三区视频| 狂野欧美激情性xxxx| av中文乱码字幕在线| 可以在线观看的亚洲视频| 久久精品91蜜桃| 久久久久国内视频| 全区人妻精品视频| 欧美乱妇无乱码| 熟女电影av网| 青草久久国产| 动漫黄色视频在线观看| 久久久久九九精品影院| 精品国产亚洲在线| 好男人在线观看高清免费视频| 全区人妻精品视频| 两性夫妻黄色片| 日韩中文字幕欧美一区二区| 男插女下体视频免费在线播放| 欧美日韩精品网址| 亚洲av熟女| 日本熟妇午夜| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 香蕉丝袜av| 黑人操中国人逼视频| 亚洲精品粉嫩美女一区| 国产1区2区3区精品| 19禁男女啪啪无遮挡网站| 国产成人av激情在线播放| 黄片小视频在线播放| 亚洲欧洲精品一区二区精品久久久| 国产主播在线观看一区二区| 黄色丝袜av网址大全| 草草在线视频免费看| 精品久久久久久成人av| 中文字幕人成人乱码亚洲影| 欧美极品一区二区三区四区| 亚洲片人在线观看| 亚洲欧洲精品一区二区精品久久久| 久久婷婷人人爽人人干人人爱| 国产激情久久老熟女| 18禁黄网站禁片免费观看直播| 女人被狂操c到高潮| 日韩欧美精品v在线| 亚洲人成伊人成综合网2020| 色在线成人网| 婷婷丁香在线五月| 99国产极品粉嫩在线观看| 1024视频免费在线观看| 精品日产1卡2卡| 99riav亚洲国产免费| 日本一本二区三区精品| 亚洲精品中文字幕在线视频| 国产精品,欧美在线| 在线观看日韩欧美| 国产亚洲精品综合一区在线观看 | a级毛片在线看网站| 国产久久久一区二区三区| 两个人的视频大全免费| 亚洲av美国av| 最近在线观看免费完整版| 亚洲国产日韩欧美精品在线观看 | 搡老熟女国产l中国老女人| 两个人免费观看高清视频| 久久午夜亚洲精品久久| 又黄又粗又硬又大视频| 99精品欧美一区二区三区四区| 久久午夜综合久久蜜桃| 波多野结衣高清作品| 免费在线观看日本一区| 亚洲无线在线观看| 亚洲一区二区三区不卡视频| 最近最新中文字幕大全免费视频| 白带黄色成豆腐渣| 久久久久九九精品影院| 老汉色∧v一级毛片| 成人av在线播放网站| 亚洲全国av大片| 国产成人一区二区三区免费视频网站| 午夜免费成人在线视频| 久久久国产成人精品二区| 丁香欧美五月| 国产精品 欧美亚洲| 在线永久观看黄色视频| 日本一本二区三区精品| 国产一区二区在线av高清观看| 亚洲精品av麻豆狂野| 久久精品综合一区二区三区| 黑人操中国人逼视频| 极品教师在线免费播放| 又黄又爽又免费观看的视频| 日韩大码丰满熟妇| 色哟哟哟哟哟哟| 一进一出抽搐动态| 无遮挡黄片免费观看| 成年女人毛片免费观看观看9| 18禁黄网站禁片免费观看直播| 成人亚洲精品av一区二区| 久久婷婷人人爽人人干人人爱| 91大片在线观看| 日本精品一区二区三区蜜桃| 窝窝影院91人妻| 亚洲av日韩精品久久久久久密| 亚洲av片天天在线观看| 日本a在线网址| 婷婷亚洲欧美| 高清毛片免费观看视频网站| 久久久久久九九精品二区国产 | avwww免费| 日本在线视频免费播放| 日本精品一区二区三区蜜桃| 久久久久久久久中文| 精品国产乱码久久久久久男人| 岛国在线观看网站| 高清毛片免费观看视频网站| 18禁观看日本| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片| 国产又黄又爽又无遮挡在线| 一二三四社区在线视频社区8| 亚洲中文日韩欧美视频| 欧美成人性av电影在线观看| 露出奶头的视频| 人人妻人人看人人澡| 精品国产乱码久久久久久男人| 婷婷亚洲欧美| 欧美乱妇无乱码| 天天躁夜夜躁狠狠躁躁| www国产在线视频色| 男女做爰动态图高潮gif福利片| 国产aⅴ精品一区二区三区波| 十八禁网站免费在线| 丁香六月欧美| 亚洲欧美一区二区三区黑人| 熟女电影av网| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| 色综合欧美亚洲国产小说| 久久精品国产清高在天天线| 一级毛片精品| av在线天堂中文字幕| 色老头精品视频在线观看| www.www免费av| ponron亚洲| 夜夜夜夜夜久久久久| 亚洲国产看品久久| 男女那种视频在线观看| 色综合站精品国产| 亚洲一区中文字幕在线| 给我免费播放毛片高清在线观看| 午夜福利在线在线| √禁漫天堂资源中文www| 亚洲国产欧洲综合997久久,| 中文资源天堂在线| 午夜福利在线观看吧| 老司机靠b影院| 国产97色在线日韩免费| 久久婷婷人人爽人人干人人爱| 高潮久久久久久久久久久不卡| 国产精品自产拍在线观看55亚洲| 成年版毛片免费区| 亚洲色图 男人天堂 中文字幕| 女人高潮潮喷娇喘18禁视频| 国产久久久一区二区三区| 国产欧美日韩精品亚洲av| 在线观看舔阴道视频| 亚洲国产欧美一区二区综合| 老熟妇仑乱视频hdxx| 国产高清videossex| 亚洲18禁久久av| 欧美日本视频| 一边摸一边抽搐一进一小说| 两个人的视频大全免费| 精品欧美一区二区三区在线| 国产私拍福利视频在线观看| 免费一级毛片在线播放高清视频| 精品免费久久久久久久清纯| 在线观看免费视频日本深夜| 精品无人区乱码1区二区| 日韩精品免费视频一区二区三区| 久久久精品大字幕| 一个人免费在线观看电影 | 在线免费观看的www视频| 中文字幕久久专区| 亚洲精品中文字幕一二三四区| 熟女电影av网| 日本撒尿小便嘘嘘汇集6| 国产蜜桃级精品一区二区三区| 国产激情欧美一区二区| 国产黄色小视频在线观看| 精品少妇一区二区三区视频日本电影| 成人av在线播放网站| av国产免费在线观看| 亚洲av成人精品一区久久| 嫩草影院精品99| 国产黄片美女视频| 国产成人aa在线观看| 国产精品香港三级国产av潘金莲| 成人一区二区视频在线观看| 国产高清有码在线观看视频 | 欧美成人一区二区免费高清观看 | 亚洲人成77777在线视频| 特级一级黄色大片| 男人的好看免费观看在线视频 | 色尼玛亚洲综合影院| 国产精品九九99| 欧美激情久久久久久爽电影| 99精品久久久久人妻精品| 婷婷精品国产亚洲av在线| svipshipincom国产片| 色综合站精品国产| 国产精品1区2区在线观看.| 黄色片一级片一级黄色片| 亚洲国产看品久久| 成人欧美大片| 母亲3免费完整高清在线观看| 国内少妇人妻偷人精品xxx网站 | 又紧又爽又黄一区二区| x7x7x7水蜜桃| 亚洲成人国产一区在线观看| 成人永久免费在线观看视频| 俺也久久电影网| 色在线成人网| 亚洲欧美日韩高清在线视频| 欧美日韩乱码在线| 我的老师免费观看完整版| 欧美一级毛片孕妇| 国产精品一区二区三区四区久久| 欧美绝顶高潮抽搐喷水| 精品欧美国产一区二区三| 亚洲国产精品999在线| 欧美久久黑人一区二区| 一个人免费在线观看电影 | 国产精品九九99| 手机成人av网站| 国产成人一区二区三区免费视频网站| 日韩欧美一区二区三区在线观看| 看黄色毛片网站| 99riav亚洲国产免费| 国产三级在线视频| 久久久精品大字幕| 亚洲av熟女| 亚洲全国av大片| 久久久久精品国产欧美久久久| 天天添夜夜摸| 国产一级毛片七仙女欲春2| 特大巨黑吊av在线直播| 国产黄a三级三级三级人| 手机成人av网站| 精品一区二区三区四区五区乱码| 他把我摸到了高潮在线观看| 国产精品亚洲美女久久久| 国产精华一区二区三区| 一级a爱片免费观看的视频| 亚洲精品av麻豆狂野| 久久天躁狠狠躁夜夜2o2o| 神马国产精品三级电影在线观看 | 国产精品精品国产色婷婷| 夜夜看夜夜爽夜夜摸| 少妇的丰满在线观看| 国产av不卡久久| 午夜激情av网站| e午夜精品久久久久久久| www.999成人在线观看| 岛国在线观看网站| 久久久精品国产亚洲av高清涩受| 丰满人妻一区二区三区视频av | 欧美国产日韩亚洲一区| 嫩草影院精品99| 不卡一级毛片| 国产精品久久电影中文字幕| 亚洲av电影在线进入| 国产亚洲精品久久久久久毛片| 成年免费大片在线观看| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 在线观看美女被高潮喷水网站 | 一进一出好大好爽视频| 女人爽到高潮嗷嗷叫在线视频| 毛片女人毛片| 久久精品夜夜夜夜夜久久蜜豆 | 国产三级黄色录像| 小说图片视频综合网站| 狂野欧美激情性xxxx| av有码第一页| 亚洲 欧美 日韩 在线 免费| 国产v大片淫在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 久久香蕉精品热| 特大巨黑吊av在线直播| 露出奶头的视频| 麻豆久久精品国产亚洲av| 熟女少妇亚洲综合色aaa.| 国产精品精品国产色婷婷| 悠悠久久av| 国产高清videossex| 狠狠狠狠99中文字幕|