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

    五階常微分方程的Petrov-Galerkin譜元法

    2017-06-05 14:17:28王金平莊清渠
    華僑大學學報(自然科學版) 2017年3期
    關(guān)鍵詞:華僑大學階數(shù)元法

    王金平, 莊清渠

    (華僑大學 數(shù)學科學學院, 福建 泉州 362021)

    五階常微分方程的Petrov-Galerkin譜元法

    王金平, 莊清渠

    (華僑大學 數(shù)學科學學院, 福建 泉州 362021)

    通過區(qū)間剖分,降低數(shù)值逼近多項式的階數(shù),構(gòu)造滿足試探函數(shù)空間和檢驗函數(shù)空間的基函數(shù),使得離散問題所對應的線性系統(tǒng)的系數(shù)矩陣是稀疏的,并可以進行有效地求解.數(shù)值算例驗證了五階常微分方程的Petrov-Galerkin譜元法的有效性和高精度.

    五階常微分方程; Petrov-Galerkin譜元法; 基函數(shù); 數(shù)值實驗

    譜方法作為數(shù)值求解微分方程的方法之一,已被廣泛應用于自然科學和工程技術(shù)問題的數(shù)值計算[1-7].近二十年來,奇數(shù)次高階微分方程基于譜方法的數(shù)值模擬也受到不少科研工作者的關(guān)注[8-13].由于奇次高階微分方程的高階項缺少對稱性,若使用配置點法求解奇次高階微分方程邊值問題,如果不能選取恰當?shù)呐渲命c或采用高階多項式進行逼近,容易導致條件數(shù)比較高,造成計算不穩(wěn)定[10].吳勝等[13]提出Legendre-Petrov-Galerkin譜元法對三階微分方程進行求解.由于將計算區(qū)間進行剖分,在每個區(qū)間上不需要使用太高的多項式階數(shù),將問題分解成一系列子問題后再用Schur補過程進行求解,計算穩(wěn)定,且具有高精度.本文研究五階常微分方程的Petrov-Galerkin譜元逼近,主要考慮方程的數(shù)值計算.

    1 Petrov-Galerkin譜元法

    記Λ=(-1,1),設(shè)α和β為正常數(shù).考慮五階常微分齊次邊值問題,即

    首先,將區(qū)間Λ剖分成K(K≥2)個子區(qū)間

    Λk=(ak-1,ak),k=1,2,…,K.

    上式中:-1=a0

    定義區(qū)間Λk到參考區(qū)間Λ的坐標變換為

    ?x∈Λk.

    定義分片多項式空間為

    PN;K(Λ)∶={u;u|Λk∈PN(Λk),k=1,2,…,K},

    上式中:PN為Λk上次數(shù)不超過N的全體多項式組成的空間.

    用N表示離散參數(shù)(N,K).定義試探函數(shù)空間和檢驗函數(shù)空間為

    VN={u;u|Λ∈PN,K(Λ);u∈C1(Λ),u(±1)=ux(±1)=ux,x(1)=0},

    WN={u;u|Λ∈PN+1,K(Λ);u∈C2(Λ),u(±1)=ux(±1)=ux,x(±1)=0}.

    則問題(1)的Petrov-Galerkin譜元逼近形式為找uN∈VN, 使得

    ?vN∈WN.

    (2)

    為了方便表達,將式(2)記為α(uN,vN)=(f,vN),記

    定義

    則可以驗證

    當k=1,2,…,K-1時,令

    而且有

    }.

    將文獻[13]求解三階微分方程的Petrov-Galerkin譜元方法的逼近過程推廣到式(2)的計算中,詳細計算過程包括如下4個步驟.

    ?

    ?

    步驟4 由式(7),(8)可得

    ?vN∈WN.

    (9)

    問題(2)被分解成兩套子區(qū)間問題(6),(7)及Schur補問題(8).因此,可以大大提高計算效率.

    (αAk-βBk+Ck)Uk=Fk,k=1,2,…,K.

    上式中:

    由Legendre多項式的正交性可知

    j=i+5,

    j=i+4,

    j=i+3,

    j=i+2,

    j=i+1,

    j=i,

    j=i-1,

    j=i-2,

    j=i-3,

    j=i-4,

    j=i-5,

    j=i-6,

    其他.

    j=i+1,

    j=i,

    j=i-1,

    j=i-2,

    j=i-3,

    其他.

    j=i-1,

    其他.

    2 數(shù)值實驗

    例1 問題(1)在x∈(-1,1)上的解析解為

    u(x)=(x-1)2sin2mπx.

    取計算參數(shù)α=β=m=1.在半log尺度下,h=1/2時,L2-誤差ε及H1-誤差隨多項式階數(shù)N的變化情況,如圖1(a)所示.由圖1(a)可知:誤差隨N呈指數(shù)收斂,這說明對于解析解,數(shù)值解具有譜收斂性質(zhì).在log-log尺度下,N=8時,L2-誤差及H1-誤差隨h的變化情況,如圖1(b)所示.由圖1(b)可知:誤差關(guān)于h呈代數(shù)衰減.

    (a) 誤差隨N的變化情況 (b) 誤差隨h的變化情況圖1 誤差分別隨N和h的變化圖Fig.1 Errors as a function of N and h respectively

    例2 考慮變系數(shù)五階齊次邊值問題,即

    此時,問題的Petrov-Galerkin譜元逼近形式為找uN∈VN, 使得

    ?vN∈WN.

    (10)

    該問題同樣可采用第1節(jié)提出的計算過程進行快速計算.但由于α(x),β(x)的存在,弱形式中前兩個積分項一般不可精確計算,因此,采用Gauss-Lobatto數(shù)值積分公式替代.

    取α(x)=x,β(x)=sin 10x,γ=1,u(x)=sin3πx進行計算.在半log尺度下,h=1/2時,誤差隨多項式階數(shù)N的變化情況,如圖2(a)所示.在log-log尺度下,N=8時,誤差隨h的變化情況,如圖2(b)所示.由圖2可知:誤差隨N呈指數(shù)收斂,誤差隨h呈代數(shù)收斂,收斂性質(zhì)跟常系數(shù)方程的計算結(jié)果類似.從而說明Petrov-Galerkin譜元逼近法對求解一些變系數(shù)五階常微分方程也是有效的.

    (a) 誤差隨N的變化情況 (b) 誤差隨h的變化情況圖2 依賴于N和h的誤差變化Fig.2 Errors as a function corresponding to N and h respectively

    3 結(jié)束語

    利用Petrov-Galerkin譜元法對五階常微分方程進行求解.首先,通過區(qū)間剖分,將變分問題轉(zhuǎn)化為一系列子問題;然后,構(gòu)造恰當?shù)脑囂胶瘮?shù)和檢驗函數(shù),由此得到稀疏的線性系統(tǒng)再進行求解.數(shù)值算例說明Petrov-Galerkin譜元逼近法對一些常系數(shù)、變系數(shù)的五階常微分方程都是有效的,且具有高精度.

    [1] GOTTLIEB D,ORSZAG S A.Numerical analysis of spectral methods: Theory and applications[M].Philadelphia:SIAM,1977:1-179.

    [2] BERNARDI C,MADAY Y.Approximations spectrales de probl`emes aux limites elliptiques[M].Berlin:Springer-Verlag,1992:1-248.

    [3] GUO Benyu.Spectral methods and their applications[M].Singapore:World Scientific,1998:100-358.

    [4] KARNIADAKIS G,SHERWIN S J.Spectral/hp element methods for computational fluid dynamics[M].London:Oxford University Press,2005:187-348.

    [5] SHEN Jie,TANG Tao.Spectral and high-order methods with applications[M].Beijing:Science Press of China,2006:183-298.

    [6] CANUTO C,HUSSAINI M Y,QUARTERONI A,etal.Spectral methods: Fundamentals in single domains[M].Berlin:Springer-Verlag,2006:401-470.

    [7] CANUTO C,HUSSAINI M Y,QUARTERONI A,etal.Spectral methods: Evolution to complex geometries and applications to fluid dynamics[M].Berlin:Springer Verlag,2007:237-357.

    [8] MA Heping,SUN Weiwei.A Legendre-Petrov-Galerkin and chebyshev collocation method for third-order differential equations[J].SIAM J NUMER ANAL,2000,38(5):1425-1438.

    [9] MA Heping,SUN Weiwei.Optimal error estiamtes of the Legendre-Petrov-Galerkin method for the Korteweg-De Vries equation[J].SIAM J NUMER ANAL,2001,39(4):1380-1394.

    [10] SHEN Jie.A new dual-Petrov-Galerkin method for third and higher odd-order differential equations: Application to the KDV equation[J].SIAM J NUMER ANAL,2003,41(5):1595-1619.

    [11] SHEN Jie,WANG Lilian.Laguerre and composite Legendre-Laguerre dual-Petrov-Galerkin methods for third-order equations[J].Discrete and Continuous Dynamical Systems Series B,2006,6(6):1381-1402.

    [12] ZHENG Chunxiong,WEN Xin,HAN Houde.Numerical solution to a linearized KDV equation on unbounded domain[J].Numerical Methods for Partial Differential Equations,2008,24(2):383-399.

    [13] 吳勝,莊清渠.三階微分方程的Legendre-Petrov-Galerkin譜元方法[J].華僑大學學報(自然科學版),2013,34(3):344-348.

    (責任編輯: 陳志賢 英文審校: 黃心中)

    Petrov-Galerkin Spectral-Element Method for Solving Fifth-Order Ordinary Differential Equations

    WANG Jinping, ZHUANG Qingqu

    (School of Mathematical Sciences, Huaqiao University, Quanzhou 362021, China)

    The polynomial order in the numerical approximation is reduced by partitioning the interval into several subintervals, and appropriate basis functions of the trial and test spaces are constructed. Which leads to a linear system with sparse coefficient matrix. Then, an efficient computational process is introduced to solve the linear system.Numerical experiment results demonstrate the high accuracy and effectiveness to the Petrov-Galerkin spectral-element method. Keywords:fifth-order ordinary differential equation; Petrov-Galerkin spectral-element method; basis functions; numerical experiments

    10.11830/ISSN.1000-5013.201703027

    2016-01-17

    莊清渠(1980-),男,副教授,博士,主要從事微分方程數(shù)值解法的研究.E-mail:qqzhuang@hqu.edu.cn.

    國家自然科學基金資助項目(11501224); 華僑大學中青年教師科研提升資助計劃(ZQN-PY201); 華僑大學研究生科研創(chuàng)新能力培育計劃項目(1400213008)

    O 241.8

    A

    1000-5013(2017)03-0435-06

    猜你喜歡
    華僑大學階數(shù)元法
    關(guān)于無窮小階數(shù)的幾點注記
    確定有限級數(shù)解的階數(shù)上界的一種n階展開方法
    換元法在解題中的運用
    基于離散元法的礦石對溜槽沖擊力的模擬研究
    重型機械(2019年3期)2019-08-27 00:58:46
    僑校雙子星:暨南大學VS華僑大學
    換元法在解題中的應用
    “微元法”在含電容器電路中的應用
    一種新的多址信道有效階數(shù)估計算法*
    關(guān)于動態(tài)電路階數(shù)的討論
    華僑大學香港校友會慶建國六十周年暨《祖國與我》聯(lián)歡晚會
    九九爱精品视频在线观看| 中文字幕熟女人妻在线| 国产三级在线视频| 国产女主播在线喷水免费视频网站 | 91久久精品电影网| 国产精品99久久久久久久久| 亚洲av成人精品一区久久| 亚洲图色成人| 51国产日韩欧美| 日韩高清综合在线| 最近中文字幕高清免费大全6| 舔av片在线| 国产国拍精品亚洲av在线观看| 久久精品国产自在天天线| 国产爱豆传媒在线观看| 欧美不卡视频在线免费观看| 久久久久免费精品人妻一区二区| 亚洲在线观看片| 亚洲图色成人| 蜜桃亚洲精品一区二区三区| 成人性生交大片免费视频hd| 亚洲精品国产av成人精品| 六月丁香七月| 一二三四中文在线观看免费高清| 舔av片在线| 久久精品综合一区二区三区| 亚洲国产欧美人成| 亚洲自偷自拍三级| 亚洲综合色惰| 热99在线观看视频| 亚洲精品日韩av片在线观看| 人人妻人人澡欧美一区二区| 精品熟女少妇av免费看| 精品国产露脸久久av麻豆 | 美女高潮的动态| 性色avwww在线观看| 69av精品久久久久久| 少妇裸体淫交视频免费看高清| 在线免费观看的www视频| 99久久中文字幕三级久久日本| 99久久精品热视频| 久久久久精品久久久久真实原创| av在线蜜桃| 人体艺术视频欧美日本| 久久久久久久久久久免费av| 亚洲欧美日韩东京热| 人妻制服诱惑在线中文字幕| 欧美激情国产日韩精品一区| 亚洲国产日韩欧美精品在线观看| 青春草视频在线免费观看| 大又大粗又爽又黄少妇毛片口| 国产高清三级在线| 99久久人妻综合| 高清毛片免费看| 亚洲美女搞黄在线观看| 啦啦啦啦在线视频资源| 免费看av在线观看网站| 欧美3d第一页| 99热这里只有是精品在线观看| 中文字幕亚洲精品专区| 日韩一区二区三区影片| 毛片女人毛片| 国产免费视频播放在线视频 | 亚洲欧洲日产国产| 亚洲三级黄色毛片| av在线老鸭窝| 久久人妻av系列| 亚洲在线观看片| 免费人成在线观看视频色| 久久精品国产亚洲av天美| 亚洲av不卡在线观看| 嘟嘟电影网在线观看| 午夜激情欧美在线| 麻豆乱淫一区二区| 成人二区视频| 极品教师在线视频| 91精品国产九色| 国产91av在线免费观看| 国产精品久久电影中文字幕| 久久精品久久久久久久性| 99久国产av精品| 国产精品蜜桃在线观看| 国产一级毛片在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成网站在线播| 天堂网av新在线| 一级爰片在线观看| 一级毛片aaaaaa免费看小| 国产一区有黄有色的免费视频 | 亚洲精品影视一区二区三区av| 亚洲av熟女| 18禁在线无遮挡免费观看视频| 一边摸一边抽搐一进一小说| 久久国内精品自在自线图片| 亚洲成人久久爱视频| 精品人妻熟女av久视频| 国产黄色视频一区二区在线观看 | 亚洲中文字幕日韩| 成年女人看的毛片在线观看| 日本免费在线观看一区| 老司机影院毛片| 国产老妇女一区| av.在线天堂| av在线蜜桃| 午夜日本视频在线| 国产三级在线视频| www日本黄色视频网| 97人妻精品一区二区三区麻豆| 日韩中字成人| 国产精品国产三级国产av玫瑰| 99在线人妻在线中文字幕| 久久久久久久午夜电影| 天美传媒精品一区二区| 欧美成人一区二区免费高清观看| 18禁裸乳无遮挡免费网站照片| 一本一本综合久久| 国产精品一二三区在线看| 婷婷色av中文字幕| 国产69精品久久久久777片| 亚洲欧美精品综合久久99| 国产精品久久电影中文字幕| 亚洲精品乱码久久久v下载方式| 日本三级黄在线观看| 久久久精品大字幕| 91精品一卡2卡3卡4卡| 国产免费男女视频| 欧美+日韩+精品| 色综合色国产| 国产熟女欧美一区二区| 久久精品人妻少妇| 久久国产乱子免费精品| 日本wwww免费看| 天堂√8在线中文| 亚洲av成人av| 在线观看美女被高潮喷水网站| or卡值多少钱| 欧美成人一区二区免费高清观看| 国产人妻一区二区三区在| 久久久久网色| 精品无人区乱码1区二区| 又爽又黄无遮挡网站| 高清在线视频一区二区三区 | 日本免费在线观看一区| 22中文网久久字幕| 国产午夜精品久久久久久一区二区三区| 热99在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品三级大全| 日韩成人av中文字幕在线观看| 午夜免费激情av| 国产69精品久久久久777片| videos熟女内射| 国产午夜精品一二区理论片| 免费看光身美女| 高清午夜精品一区二区三区| 又粗又爽又猛毛片免费看| av免费在线看不卡| 亚洲婷婷狠狠爱综合网| 午夜福利视频1000在线观看| 99热6这里只有精品| 亚洲国产成人一精品久久久| av在线老鸭窝| 伊人久久精品亚洲午夜| 一区二区三区高清视频在线| 国产免费男女视频| 国语自产精品视频在线第100页| 成人毛片60女人毛片免费| 亚洲欧美日韩无卡精品| 99热网站在线观看| 听说在线观看完整版免费高清| 免费观看性生交大片5| 日韩高清综合在线| 最近的中文字幕免费完整| 亚洲精品456在线播放app| 亚洲av成人精品一区久久| 看免费成人av毛片| 成人毛片a级毛片在线播放| 少妇被粗大猛烈的视频| 亚洲国产精品专区欧美| 国产亚洲午夜精品一区二区久久 | 色尼玛亚洲综合影院| 国产精品野战在线观看| 国产毛片a区久久久久| 精品国产露脸久久av麻豆 | 中国美白少妇内射xxxbb| 亚洲国产精品合色在线| 波野结衣二区三区在线| 亚洲精品aⅴ在线观看| 日韩精品青青久久久久久| 欧美bdsm另类| 国产黄a三级三级三级人| 免费大片18禁| 最近手机中文字幕大全| 热99re8久久精品国产| 亚洲不卡免费看| 国产免费又黄又爽又色| 黄色一级大片看看| 久久久久久九九精品二区国产| 欧美日本亚洲视频在线播放| 久久久久久伊人网av| 亚洲av电影在线观看一区二区三区 | 天美传媒精品一区二区| 国产精品熟女久久久久浪| 成人亚洲精品av一区二区| 亚洲av中文字字幕乱码综合| 黄色配什么色好看| 欧美一区二区精品小视频在线| 国产精品福利在线免费观看| 七月丁香在线播放| 日韩成人av中文字幕在线观看| 亚洲欧美日韩高清专用| 国产精品永久免费网站| a级一级毛片免费在线观看| 亚洲中文字幕日韩| 精品人妻视频免费看| 日韩欧美三级三区| 我要搜黄色片| 欧美色视频一区免费| 精品久久久久久久久久久久久| 色5月婷婷丁香| av在线蜜桃| 婷婷六月久久综合丁香| 激情 狠狠 欧美| 一级毛片我不卡| 一级av片app| 联通29元200g的流量卡| 黄色日韩在线| 欧美日本亚洲视频在线播放| АⅤ资源中文在线天堂| 超碰av人人做人人爽久久| 成人性生交大片免费视频hd| 国产91av在线免费观看| 久久久欧美国产精品| 三级男女做爰猛烈吃奶摸视频| 欧美又色又爽又黄视频| 在线免费观看的www视频| 99久久无色码亚洲精品果冻| 欧美丝袜亚洲另类| 国产乱人偷精品视频| 久久久久久久久久成人| 亚洲国产日韩欧美精品在线观看| 熟女人妻精品中文字幕| 成人毛片60女人毛片免费| 青春草国产在线视频| 熟妇人妻久久中文字幕3abv| 久久久久久久久大av| 亚洲欧美成人综合另类久久久 | 99久久精品国产国产毛片| 国产视频内射| 成人午夜精彩视频在线观看| 国产极品精品免费视频能看的| 色综合亚洲欧美另类图片| 99国产精品一区二区蜜桃av| 欧美丝袜亚洲另类| 国产伦精品一区二区三区四那| 白带黄色成豆腐渣| 欧美日韩精品成人综合77777| 深爱激情五月婷婷| 熟女人妻精品中文字幕| 一级爰片在线观看| 精品不卡国产一区二区三区| 久久久久久伊人网av| 亚洲av.av天堂| 有码 亚洲区| 国产成人aa在线观看| videos熟女内射| 亚洲最大成人中文| 波多野结衣高清无吗| 成人三级黄色视频| 久久精品91蜜桃| 亚洲美女视频黄频| 国产精品三级大全| 国产人妻一区二区三区在| 可以在线观看毛片的网站| 久久精品91蜜桃| 日日撸夜夜添| 国产精品久久久久久av不卡| 18禁动态无遮挡网站| 日日摸夜夜添夜夜添av毛片| 天天躁夜夜躁狠狠久久av| 亚洲国产欧美人成| 亚洲欧美中文字幕日韩二区| 国产综合懂色| 高清午夜精品一区二区三区| 精品人妻视频免费看| 国产精品女同一区二区软件| 免费在线观看成人毛片| 亚洲成av人片在线播放无| 一区二区三区高清视频在线| 日本猛色少妇xxxxx猛交久久| 国模一区二区三区四区视频| 色尼玛亚洲综合影院| 色5月婷婷丁香| 国产一区二区三区av在线| 日本爱情动作片www.在线观看| 久久精品久久久久久久性| 特级一级黄色大片| 黄片无遮挡物在线观看| 噜噜噜噜噜久久久久久91| 亚洲人与动物交配视频| 亚洲人成网站在线观看播放| 国产三级中文精品| 中文字幕熟女人妻在线| 一级黄片播放器| 人人妻人人看人人澡| 日韩精品青青久久久久久| 最近视频中文字幕2019在线8| 少妇熟女欧美另类| 人人妻人人看人人澡| 蜜臀久久99精品久久宅男| 国产亚洲最大av| 久99久视频精品免费| 免费一级毛片在线播放高清视频| 国产免费一级a男人的天堂| 麻豆国产97在线/欧美| 欧美xxxx性猛交bbbb| 免费人成在线观看视频色| 亚洲精品乱码久久久久久按摩| 久久久久久久久久黄片| 性插视频无遮挡在线免费观看| 欧美成人精品欧美一级黄| 97超碰精品成人国产| 亚洲精品日韩在线中文字幕| 欧美一区二区亚洲| 桃色一区二区三区在线观看| 国产精品乱码一区二三区的特点| 成人三级黄色视频| 91在线精品国自产拍蜜月| 99久国产av精品| 亚洲经典国产精华液单| 99热网站在线观看| 午夜免费激情av| 日韩欧美国产在线观看| 日本免费在线观看一区| 国国产精品蜜臀av免费| 最近2019中文字幕mv第一页| 秋霞伦理黄片| 国产极品天堂在线| 欧美一区二区亚洲| 欧美极品一区二区三区四区| .国产精品久久| 五月玫瑰六月丁香| 非洲黑人性xxxx精品又粗又长| 高清视频免费观看一区二区 | 久久欧美精品欧美久久欧美| 日韩成人av中文字幕在线观看| 欧美极品一区二区三区四区| 国产精品1区2区在线观看.| 爱豆传媒免费全集在线观看| 久久久色成人| 久久久欧美国产精品| 91av网一区二区| 一二三四中文在线观看免费高清| videos熟女内射| 边亲边吃奶的免费视频| 淫秽高清视频在线观看| 国产精品久久久久久久久免| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久久黄片| 水蜜桃什么品种好| 免费一级毛片在线播放高清视频| 国产精品无大码| av卡一久久| 国产亚洲av片在线观看秒播厂 | 男女那种视频在线观看| 亚洲第一区二区三区不卡| 色哟哟·www| 国产黄色视频一区二区在线观看 | 亚洲精品国产成人久久av| 干丝袜人妻中文字幕| 一个人看的www免费观看视频| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 波多野结衣巨乳人妻| 国产午夜精品久久久久久一区二区三区| 免费观看性生交大片5| 免费黄网站久久成人精品| 欧美三级亚洲精品| 国产高清国产精品国产三级 | 国产淫语在线视频| 国产极品精品免费视频能看的| 国产亚洲最大av| 久久草成人影院| 色网站视频免费| 欧美性猛交黑人性爽| 精品99又大又爽又粗少妇毛片| 熟女电影av网| 亚洲精品国产成人久久av| 午夜福利成人在线免费观看| 在线天堂最新版资源| 91av网一区二区| 国产在视频线精品| 亚洲丝袜综合中文字幕| 2022亚洲国产成人精品| 国模一区二区三区四区视频| 亚洲人成网站在线观看播放| eeuss影院久久| 毛片女人毛片| 中文字幕av在线有码专区| 国产午夜精品论理片| 亚洲精品成人久久久久久| 少妇猛男粗大的猛烈进出视频 | 亚洲av.av天堂| 亚洲成人久久爱视频| 大香蕉久久网| 日韩大片免费观看网站 | 美女被艹到高潮喷水动态| 黄片无遮挡物在线观看| 插阴视频在线观看视频| 免费一级毛片在线播放高清视频| 精品人妻视频免费看| 午夜福利视频1000在线观看| 久久久欧美国产精品| 草草在线视频免费看| 非洲黑人性xxxx精品又粗又长| 久久鲁丝午夜福利片| 亚洲综合色惰| 国产乱来视频区| 人妻制服诱惑在线中文字幕| 亚洲18禁久久av| 久久久久久久国产电影| 人妻夜夜爽99麻豆av| av在线天堂中文字幕| 久久久久久久久久久丰满| 亚洲激情五月婷婷啪啪| 一个人免费在线观看电影| 亚洲av成人精品一区久久| 2021天堂中文幕一二区在线观| 春色校园在线视频观看| 国产成人一区二区在线| 亚洲怡红院男人天堂| 国产美女午夜福利| 少妇高潮的动态图| 男女边吃奶边做爰视频| videos熟女内射| 少妇人妻一区二区三区视频| 在现免费观看毛片| 中国美白少妇内射xxxbb| 日韩欧美三级三区| 欧美丝袜亚洲另类| 精品国产露脸久久av麻豆 | 成人三级黄色视频| 久久久色成人| 精品少妇黑人巨大在线播放 | 国产高清三级在线| 精品久久久久久久久亚洲| 中文乱码字字幕精品一区二区三区 | 欧美xxxx性猛交bbbb| 黄色日韩在线| 噜噜噜噜噜久久久久久91| 最近最新中文字幕大全电影3| 国产精品福利在线免费观看| 岛国在线免费视频观看| 亚洲av电影不卡..在线观看| 日韩人妻高清精品专区| 最后的刺客免费高清国语| 免费av观看视频| 又粗又硬又长又爽又黄的视频| 亚洲人成网站在线播| 自拍偷自拍亚洲精品老妇| 最后的刺客免费高清国语| 亚洲最大成人中文| 在线免费观看的www视频| 国产一级毛片在线| 中文字幕av成人在线电影| 春色校园在线视频观看| 日韩国内少妇激情av| 亚洲av免费在线观看| 蜜桃亚洲精品一区二区三区| 身体一侧抽搐| av国产久精品久网站免费入址| 亚洲,欧美,日韩| 高清av免费在线| 亚洲国产精品国产精品| 亚洲乱码一区二区免费版| 好男人在线观看高清免费视频| 亚洲国产精品合色在线| 久久久久久国产a免费观看| 国产乱人偷精品视频| av免费观看日本| 99热这里只有是精品在线观看| 身体一侧抽搐| 97热精品久久久久久| 国产伦在线观看视频一区| 一边亲一边摸免费视频| av线在线观看网站| 精华霜和精华液先用哪个| 男人和女人高潮做爰伦理| www.av在线官网国产| 最近视频中文字幕2019在线8| 欧美97在线视频| 成人国产麻豆网| 亚洲欧美一区二区三区国产| 老司机影院毛片| 一个人免费在线观看电影| 久久久久久久久中文| 有码 亚洲区| 久久久精品大字幕| 久久精品夜夜夜夜夜久久蜜豆| 1024手机看黄色片| 黄色一级大片看看| 日本免费a在线| 国产又色又爽无遮挡免| 国产精品国产三级国产av玫瑰| 日本黄色片子视频| 午夜福利在线在线| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区三区四区免费观看| 亚洲欧洲国产日韩| 国产精品一二三区在线看| 精品不卡国产一区二区三区| 高清av免费在线| 又爽又黄a免费视频| 26uuu在线亚洲综合色| 女人被狂操c到高潮| 婷婷色麻豆天堂久久 | 99视频精品全部免费 在线| 六月丁香七月| 亚洲欧美日韩东京热| 成年女人永久免费观看视频| 啦啦啦观看免费观看视频高清| 精品国产三级普通话版| 老司机福利观看| 亚洲国产精品sss在线观看| 一级黄片播放器| 成人毛片60女人毛片免费| 男人的好看免费观看在线视频| 欧美最新免费一区二区三区| 久久韩国三级中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看性生交大片5| 黄色配什么色好看| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 长腿黑丝高跟| 又黄又爽又刺激的免费视频.| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产午夜精品论理片| 有码 亚洲区| 亚洲精品aⅴ在线观看| 卡戴珊不雅视频在线播放| 看免费成人av毛片| 久久久久久大精品| 看黄色毛片网站| 欧美精品国产亚洲| 啦啦啦观看免费观看视频高清| 美女黄网站色视频| 国产色爽女视频免费观看| 精品久久久久久电影网 | 久久99精品国语久久久| 少妇人妻精品综合一区二区| 亚洲精品乱码久久久久久按摩| 狂野欧美白嫩少妇大欣赏| 国产三级中文精品| 非洲黑人性xxxx精品又粗又长| 天天躁夜夜躁狠狠久久av| 久久久色成人| av在线蜜桃| 国产老妇伦熟女老妇高清| 亚洲精品aⅴ在线观看| 国产淫片久久久久久久久| 纵有疾风起免费观看全集完整版 | 国产av不卡久久| 免费观看精品视频网站| 国产黄a三级三级三级人| 久久这里有精品视频免费| 国产精品日韩av在线免费观看| 亚洲av中文字字幕乱码综合| 色综合站精品国产| 久久久久久久久中文| 成人漫画全彩无遮挡| 久久久久免费精品人妻一区二区| 舔av片在线| 久久精品国产亚洲av涩爱| 国产高清有码在线观看视频| www.av在线官网国产| 91久久精品国产一区二区成人| 三级经典国产精品| 成人无遮挡网站| 国产白丝娇喘喷水9色精品| 久久午夜福利片| 插逼视频在线观看| 大香蕉97超碰在线| 国内精品美女久久久久久| 国产淫片久久久久久久久| 日韩av在线免费看完整版不卡| 99热6这里只有精品| 亚洲精品,欧美精品| 欧美丝袜亚洲另类| 深夜a级毛片| 久久久久久久国产电影| 久久99热这里只有精品18| 直男gayav资源| 亚洲国产精品专区欧美| av在线老鸭窝| 长腿黑丝高跟| 国产精品人妻久久久久久| 国产精品女同一区二区软件| 日韩,欧美,国产一区二区三区 | 亚洲在久久综合| 岛国在线免费视频观看| av福利片在线观看| 欧美成人一区二区免费高清观看| 精品不卡国产一区二区三区| 午夜老司机福利剧场| av在线蜜桃| 国产亚洲一区二区精品| 少妇丰满av| 欧美成人午夜免费资源| 搡女人真爽免费视频火全软件| 天堂影院成人在线观看| 亚洲国产欧洲综合997久久,| 日韩欧美在线乱码| 亚洲欧美日韩卡通动漫| 成人亚洲欧美一区二区av|