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

    解析延拓課程教學(xué)改革淺探*
    ——兼談最大熵方法

    2017-02-15 08:57:04呂建平
    物理通報(bào) 2017年2期
    關(guān)鍵詞:松原關(guān)聯(lián)解析

    呂建平

    (安徽師范大學(xué)物理系 安徽 蕪湖 241000)

    解析延拓課程教學(xué)改革淺探*
    ——兼談最大熵方法

    呂建平

    (安徽師范大學(xué)物理系 安徽 蕪湖 241000)

    針對數(shù)學(xué)物理方法課程教學(xué)中學(xué)生普遍反映比較抽象的解析延拓部分,設(shè)計(jì)以數(shù)值模擬方法為技術(shù)路線的教學(xué)方案.通過物理模型的構(gòu)建,介紹線性響應(yīng)理論和最大熵解析延拓方法的基本原理,全面設(shè)計(jì)課堂教學(xué)環(huán)節(jié),展示對具體實(shí)例的計(jì)算結(jié)果.最后,對課程實(shí)踐過程中的有待提高之處進(jìn)行反思和探討.

    數(shù)學(xué)物理方法 解析延拓 最大熵方法

    數(shù)學(xué)物理方法課程是高等院校物理類本科專業(yè)必修的重要課程之一.然而,對中等程度的學(xué)生而言,學(xué)好這門課程并不容易.學(xué)生在學(xué)習(xí)過程中既需要高等數(shù)學(xué)等前期數(shù)學(xué)課程的良好基礎(chǔ),又需要對一些經(jīng)典物理模型有清晰的認(rèn)識[1,2].由于上述客觀原因,加之筆者水平所限,在承擔(dān)安徽師范大學(xué)卓越理科實(shí)驗(yàn)班數(shù)學(xué)物理方法課程之初,筆者常在課堂教學(xué)設(shè)計(jì)過程中感到困惑.通過課堂教學(xué)的反饋,常獲悉學(xué)生在學(xué)習(xí)過程中感到概念過于抽象,無法和具體的物理系統(tǒng)進(jìn)行有效聯(lián)系.

    和國內(nèi)多數(shù)高校一樣,我校數(shù)學(xué)物理方法課程在低年級開課,此時(shí)學(xué)生對一些重要物理問題的圖像和機(jī)理還沒有全面的認(rèn)識.這要求教師在教學(xué)過程中,在有限的教學(xué)時(shí)長內(nèi)補(bǔ)充介紹相關(guān)物理模型.而幾乎所有經(jīng)典教材在講述偏微分方程部分時(shí)都會引入包括電導(dǎo)、熱傳導(dǎo)在內(nèi)的輸運(yùn)過程作為主要物理模型.這引發(fā)筆者的思考:既然學(xué)生在數(shù)學(xué)物理方法課的學(xué)習(xí)過程中須對輸運(yùn)過程的物理圖像有基本概念,可否將輸運(yùn)過程的相關(guān)模型拓展到偏微分方程以外的教學(xué)內(nèi)容?此為筆者探索相關(guān)課程設(shè)計(jì)的背景.

    本文討論解析延拓的教學(xué)設(shè)計(jì),旨在借助于較為具體的物理圖像來避免教學(xué)內(nèi)容過于抽象化,以探索提升教學(xué)質(zhì)量的途徑.沿用吳崇試先生《數(shù)學(xué)物理方法》教材中關(guān)于解析延拓的概念[2]: 兩個函數(shù)f1(z)和f2(z) 分別在定義域g1和g2內(nèi)解析,且在g1和g2的公共區(qū)域等價(jià),則f2(z)可稱為f1(z)在g2區(qū)域的解析延拓,且f1(z)可稱為f2(z)在g1區(qū)域的解析延拓.又由同一性原理可知,對定義在同一個區(qū)域的兩個解析函數(shù),如果它們在該區(qū)域內(nèi)的無窮多個點(diǎn)上相等,且這些點(diǎn)中有聚點(diǎn)存在,這兩個函數(shù)必然在整個區(qū)域上等價(jià).在多體物理系統(tǒng)的研究中常須借助于多種關(guān)聯(lián)函數(shù).在頻域復(fù)平面上,延遲關(guān)聯(lián)函數(shù)和相應(yīng)的松原關(guān)聯(lián)函數(shù)在無窮多個點(diǎn)上相等,且這些點(diǎn)中存在聚點(diǎn)(此處聚點(diǎn)為虛軸上正方向的無窮遠(yuǎn)點(diǎn)i∞),因而必然在共同解析區(qū)域上等價(jià),且松原關(guān)聯(lián)函數(shù)可延拓到延遲關(guān)聯(lián)函數(shù)的解析區(qū)域.在實(shí)際計(jì)算中,由于一些數(shù)值算法在計(jì)算松原關(guān)聯(lián)函數(shù)時(shí)具有顯著的優(yōu)勢.因此,在通常情況下,我們總是先直接計(jì)算松原關(guān)聯(lián)函數(shù),再借助于解析延拓手段探索延遲關(guān)聯(lián)函數(shù).

    基于上述論證,同時(shí)借助于數(shù)值模擬方法,筆者將針對輸運(yùn)問題設(shè)計(jì)解析延拓教學(xué)方案,進(jìn)而在實(shí)際教學(xué)中進(jìn)行小范圍實(shí)踐.下文將對教學(xué)過程的設(shè)計(jì)原理和基本環(huán)節(jié)進(jìn)行介紹,并對教學(xué)實(shí)踐進(jìn)行反思和總結(jié).

    1 物理問題的提出

    輸運(yùn)過程有廣泛的范疇,在物理學(xué)的教學(xué)和研究中具有重要的地位.輸運(yùn)過程研究的常用理論框架之一為線性響應(yīng)理論.線性響應(yīng)理論對近平衡態(tài)系統(tǒng)的輸運(yùn)行為有良好的描述.根據(jù)線性響應(yīng)理論,我們只需要知道以時(shí)間為變量的電流-電流延遲關(guān)聯(lián)函數(shù),就可根據(jù)久保公式方便地計(jì)算出電導(dǎo)率[3].然而,如何求得電流-電流關(guān)聯(lián)函數(shù)?首先,借助于數(shù)值模擬可獲得電流-電流松原關(guān)聯(lián)函數(shù)[4].接下來,可根據(jù)解析延拓手段由松原關(guān)聯(lián)函數(shù)獲得相應(yīng)區(qū)域的延遲關(guān)聯(lián)函數(shù).關(guān)于電流-電流松原關(guān)聯(lián)函數(shù)的計(jì)算,筆者通過前期科研工作已有方法上的積累,既可將程序代碼以軟件的形式提供給學(xué)生使用[5~7],也可為學(xué)生直接提供松原關(guān)聯(lián)函數(shù)的數(shù)據(jù).無論哪種做法,都不會直接影響學(xué)生在解析延拓方面的訓(xùn)練.

    在對物理系統(tǒng)平衡態(tài)性質(zhì)的討論中,基于松原虛時(shí)的路徑積分蒙特卡洛方法已被廣泛采用.然而,本教學(xué)過程涉及輸運(yùn)過程的數(shù)值實(shí)驗(yàn),對蒙特卡洛方法有更高要求.原因是顯然的:此處蒙特卡洛數(shù)據(jù)只是解析延拓的輸入數(shù)據(jù),而不再是與實(shí)驗(yàn)可比的最終結(jié)果.一旦蒙特卡洛數(shù)據(jù)的精度不高,或抽樣數(shù)據(jù)之間的關(guān)聯(lián)性不能被明確把握,相關(guān)的解析延拓過程就難以獲得穩(wěn)定可靠的結(jié)果.因此,高效的蒙特卡洛算法是必要的技術(shù)手段.在課程教學(xué)中,我們將使用高效的蠕蟲型連續(xù)虛時(shí)蒙特卡洛算法,該算法通過新自由度的引入使得位形空間得以擴(kuò)展,進(jìn)而以馬爾科夫鏈的形式產(chǎn)生可供抽樣統(tǒng)計(jì)的位形.這種新位形的產(chǎn)生方式具有克服臨界慢化的功能,因而算法效率將獲顯著提高.在此過程中,學(xué)生須知道松原關(guān)聯(lián)函數(shù)的形式特點(diǎn):若考慮頻域復(fù)平面,在虛軸的松原頻率點(diǎn)上松原關(guān)聯(lián)函數(shù)可以直接計(jì)算,且等同于相應(yīng)的延遲關(guān)聯(lián)函數(shù).

    如前文所述,在物理上更重要的是實(shí)時(shí)延遲電流-電流關(guān)聯(lián)函數(shù).若已知松原關(guān)聯(lián)函數(shù)在上半平面內(nèi)一些點(diǎn)上的數(shù)值結(jié)果,如何通過解析延拓,獲得無限靠近實(shí)軸時(shí)的關(guān)聯(lián)函數(shù)數(shù)據(jù)?不失一般性,在已獲得松原關(guān)聯(lián)函數(shù)數(shù)據(jù)的情況下,可通過最大熵解析延拓辦法獲得以實(shí)頻為變量的延遲關(guān)聯(lián)函數(shù),進(jìn)而求得電導(dǎo)率.

    2 最大熵解析延拓方法

    通過多體物理理論中常用的萊曼(Lehmann)譜分解可知,在遠(yuǎn)離實(shí)軸的上半頻域復(fù)平面內(nèi),延遲關(guān)聯(lián)函數(shù)在形式上等價(jià)于相應(yīng)的松原關(guān)聯(lián)函數(shù).在松原頻率上,后者可由量子蒙特卡洛模擬直接測量.為探索實(shí)軸上的信息,可作解析延拓iωn→ω+i0+,其中ωn為松原頻率.我們擬采用的最大熵方法是最具代表性的解析延拓方法之一[8,9].最大熵方法的主要思路是運(yùn)用貝葉斯推理,把尋找最可能解析延拓結(jié)果的問題,轉(zhuǎn)化為求解“自由能”最小值的問題,它的穩(wěn)定性已被大量文獻(xiàn)結(jié)果所證實(shí).值得一提的是,最大熵方法對具有數(shù)值不確定性的數(shù)據(jù)具有強(qiáng)大的處理能力.在實(shí)際應(yīng)用中,我們借助于高效量子蒙特卡洛方法產(chǎn)生的松原關(guān)聯(lián)函數(shù)具有高度的精確性,完全可以通過最大熵方法獲得穩(wěn)定可靠的實(shí)頻輸出數(shù)據(jù).另外,我們擬定的課堂教學(xué)方案中也包含了可靠的誤差估計(jì)方法.他山之石,可以攻玉.在未來的教學(xué)過程中,我們也有可能借鑒其他的解析延拓方法,如隨機(jī)優(yōu)化方法[10]和隨機(jī)抽樣方法[11],將它們作為傳統(tǒng)最大熵方法的重要補(bǔ)充.

    最大熵方法最早發(fā)展于工程科學(xué)領(lǐng)域,后被物理學(xué)家借鑒,用來解決具有數(shù)值不確定性數(shù)據(jù)的解析延拓問題[8].首先,根據(jù)相關(guān)問題情境,引入“自由能”的概念.“自由能”的定義為F=χ2-AS,在形式上與熱力學(xué)中的自由能類似,其中χ2是描述擬合過程的協(xié)方差,S是相對于預(yù)設(shè)默認(rèn)模型的熵密度,A是常系數(shù),它的選取有不同的判定方法.而S類似于熱力學(xué)中的熵,此即為最大熵方法名稱的由來.我們不難想象,當(dāng)A=0時(shí),F(xiàn)中協(xié)方差χ2占主導(dǎo)地位,解析延拓過程就完全退化為數(shù)據(jù)擬合過程,擬合結(jié)果通常為噪聲結(jié)果,不具有預(yù)測能力;當(dāng)A很大時(shí),S占主導(dǎo)作用,擬合結(jié)果對應(yīng)于平庸的默認(rèn)模式.這兩種極限下,解析延拓結(jié)果都沒有意義.實(shí)際上,根據(jù)χ2隨熵權(quán)重變化的依賴關(guān)系可區(qū)分3個截然不同的區(qū)域:噪聲擬合區(qū)域、信息擬合區(qū)域和默認(rèn)模型區(qū)域.對這3個區(qū)域的探討和分析,將有助于我們正確使用并最終發(fā)展最大熵方法[9].

    在通過編程實(shí)現(xiàn)解析延拓的過程中,我們會盡可能地納入有效信息而完全避免噪聲信息.為此,我們把噪聲擬合區(qū)域和信息擬合區(qū)域的分界點(diǎn)作為熵權(quán)重的最佳選擇,并通過分析不同類型的熵權(quán)重選擇方法,探索最佳的解析延拓方案.從理論上講,無論輸入數(shù)據(jù)是否有數(shù)值上的不確定性,最大熵方法都應(yīng)該給出唯一穩(wěn)定的解析延拓結(jié)果.但值得注意的是,熵權(quán)重的選取會影響誤差估計(jì)的可靠性.而目前較為流行的熵權(quán)重選取方法是基于經(jīng)驗(yàn)判據(jù),即令協(xié)方差χ2等于樣本數(shù).該判據(jù)僅在樣本數(shù)很大且χ2精確可知的情況下較好地工作,因而不是我們的選擇.

    3 課堂教學(xué)的實(shí)例

    根據(jù)輸運(yùn)過程的基本特征,借助于相關(guān)的理論和方法,并考慮課堂教學(xué)的實(shí)際情況,我們擬定了如圖1所示的課堂教學(xué)基本流程.

    其中T為溫度,kB為玻爾茲曼常數(shù).若僅關(guān)注平衡態(tài)行為,通過附加的模擬可確定此時(shí)系統(tǒng)處于正常流體態(tài),其低溫下的超流特性被熱漲落破壞,但仍然保留普通流體的一些特征,包括擴(kuò)散性和可壓縮性.我們在前文中提到的蒙特卡洛方法對該系統(tǒng)具有較為強(qiáng)大的刻畫能力.

    圖1 課程教學(xué)的基本流程

    通過數(shù)值模擬,我們獲得了系統(tǒng)在平衡態(tài)的一系列數(shù)據(jù).在圖2中,我們展示了以虛頻為自變量的電流-電流松原關(guān)聯(lián)函數(shù).不難發(fā)現(xiàn),該函數(shù)以較快的速度衰減.又由線性響應(yīng)理論可知,以實(shí)頻為變量的電流-電流延遲關(guān)聯(lián)函數(shù)的虛部與電導(dǎo)率的規(guī)則項(xiàng)對應(yīng).為獲得該延遲關(guān)聯(lián)函數(shù)的虛部,我們采用上文中提到的、經(jīng)過優(yōu)化的最大熵方法進(jìn)行解析延拓.如圖3所示,我們獲得了電導(dǎo)率隨實(shí)頻率變化的關(guān)系.和我們的預(yù)期一致,相關(guān)的解析延拓手段展現(xiàn)出高度的精確性和穩(wěn)定性,這一點(diǎn)體現(xiàn)為解析延拓結(jié)果極窄的上下界范圍.最后,我們將直流(ω=0)電導(dǎo)率的規(guī)則項(xiàng)確定為0.085(2).

    圖2 高溫下正常流體態(tài)的松原流-流關(guān)聯(lián)強(qiáng)度(已設(shè)=1,

    圖3 高溫下正常流體態(tài)的電導(dǎo)率,所取參數(shù)與圖2相同

    4 總結(jié)和反思

    通過將數(shù)值模擬方法引入到解析延拓課程相關(guān)章節(jié)的教學(xué)過程中,學(xué)生可獲得對解析延拓理論及其實(shí)踐方法的直觀印象.通過與實(shí)際物理問題的結(jié)合,學(xué)生意識到相關(guān)教學(xué)內(nèi)容在現(xiàn)代物理研究中的重要性,因而在一定程度上增強(qiáng)了學(xué)習(xí)的積極性.然而,筆者也清楚地認(rèn)識到教學(xué)過程中的若干不足之處.例如,產(chǎn)生松原關(guān)聯(lián)函數(shù)的過程對學(xué)生而言類似于一個“黑匣子”,學(xué)生很難全面地理解數(shù)值模擬的細(xì)節(jié),這可能會在一定程度上削減學(xué)生參與相關(guān)數(shù)值實(shí)驗(yàn)的熱情.在數(shù)學(xué)物理方法課程的教學(xué)中,如何在有限的教學(xué)時(shí)長內(nèi),讓學(xué)生通過課本知識獲得完整解決一個實(shí)際物理問題的體驗(yàn)?以本文為起點(diǎn),在后續(xù)的教學(xué)實(shí)踐過程中,筆者將繼續(xù)修正相關(guān)的課程教學(xué)環(huán)節(jié),以期獲得教學(xué)效率的持續(xù)提高.

    1 四川大學(xué)數(shù)學(xué)學(xué)院高等數(shù)學(xué)、微分方程教研室.高等數(shù)學(xué)(第4冊)(第3版).北京: 高等教育出版社, 2010.76~98

    2 吳崇試.數(shù)學(xué)物理方法.北京:北京大學(xué)出版社, 1999.137~144

    3 列夫·達(dá)維多維奇·朗道,E·M·栗弗席茲.理論物理學(xué)教程:統(tǒng)計(jì)物理學(xué)(第5版).束仁貴,束莼,譯.北京:高等教育出版社,2011.337~347

    4K.Binder(edit.).TheMonteCarloMethodinConden-sedMatterPhysics.Germany:Springer, 2012

    5J.P.Lv,Q.H.Chen,Y.Deng.Two-specieshardc-orebosonsonthetriangularlattice:AquantumMonteCarlostudy.Phys.Rev.A, 2014,89:013628

    6J.P.Lv,G.Chen,Y.Deng.Z.Y.Meng,CoulombliquidphasesofbosonicclusterMottinsulatorsonapyrochlorelattice.Phys.Rev.Lett.115,2015:037202.

    7J.P.Lv.Kaleidoscopeofsupersolidphasesofinterac-tinghard-corebosonsonthedicelattice.TheEuropeanPhysicalJournalB,2015,88:13

    8M.Jarrell,J.E.Gubernatis.Bayesianinferenceandtheanalyticcontinuationofimaginary-timequantumMonteCarlodata.PhysicsReports1996, 269:133

    9D.Bergeron,A.M.S.Tremblay.Algorithmsforopti-mizedmaximumentropyanddiagnostictoolsforanalyticcontinuation.arXiv, 2015, 1507.01012

    10A.S.Mishchenko.StochasticOptimizationforAnalyt-icalContinuation:WhenaprioriKnowledgeisMissing, 14thchapterinCorrelatedElectrons:FromModelstoMaterials.ed.byE.Pavarini,W.Koch,F.AndersandM.Jarrell,VerlagdesForschungszentrumJülich, 2012: 1~24

    11A.W.Sandvik.Stochasticmethodforanalyticcontinu-ationofquantumMonteCarlodata.Phys.Rev.B,1998, 57:10287

    ExplorationontheTeachingReformofAnalyticContinuationCourse——AndTalkingabouttheMaximumEntropyMethod

    LvJianping

    (DepartmentofPhysics,AnhuiNormalUniversity,Wuhu,Anhui241000)

    For the analytic continuation part in Mathematical Physics course which is more or less abstract to students, we design the course teaching based on numerical simulation methods. In this paper, via the construction of physical model, we introduce the basic aspects for linear response theory and maximum entropy analytic continuation method. We design the specific procedure in course teaching process, and demonstrate the numerical results for a specific example. Finally, we discuss the possible directions to improve present course teaching process.

    Mathematical Physics; analytic continuation; maximum entropy method

    *國家自然科學(xué)基金,項(xiàng)目編號:11405003,11147013

    呂建平(1983- ),男,博士,主要從事數(shù)學(xué)物理方法、群論和大學(xué)物理等課程的教學(xué),以及凝聚態(tài)物理理論研究.

    2016-06-19)

    猜你喜歡
    松原關(guān)聯(lián)解析
    三角函數(shù)解析式中ω的幾種求法
    河南松原生物科技有限公司
    中國釀造(2019年9期)2019-10-08 05:44:02
    “一帶一路”遞進(jìn),關(guān)聯(lián)民生更緊
    愛的花旗松素 走進(jìn)吉林松原市場
    睡夢解析儀
    奇趣搭配
    電競初解析
    商周刊(2017年12期)2017-06-22 12:02:01
    智趣
    讀者(2017年5期)2017-02-15 18:04:18
    啊,松原
    黃河之聲(2016年13期)2016-11-09 01:42:16
    相機(jī)解析
    国产熟女午夜一区二区三区 | 全区人妻精品视频| 国产高清不卡午夜福利| 欧美日韩成人在线一区二区| 久久久久久久久大av| 欧美日韩视频高清一区二区三区二| 精品国产国语对白av| 国产精品久久久久久精品电影小说| 香蕉精品网在线| h视频一区二区三区| 欧美最新免费一区二区三区| 久久精品国产自在天天线| 午夜久久久在线观看| 精品久久久精品久久久| 国产乱人偷精品视频| 99re6热这里在线精品视频| 午夜福利视频在线观看免费| 黑丝袜美女国产一区| 丝袜喷水一区| 观看美女的网站| 一级片'在线观看视频| 美女国产视频在线观看| 国产永久视频网站| 人体艺术视频欧美日本| 又大又黄又爽视频免费| 少妇猛男粗大的猛烈进出视频| 国产精品一区www在线观看| 久久久a久久爽久久v久久| 性高湖久久久久久久久免费观看| 最新的欧美精品一区二区| 中文精品一卡2卡3卡4更新| 99九九在线精品视频| 蜜桃国产av成人99| 国产在线视频一区二区| 日韩中字成人| 国产成人免费无遮挡视频| 少妇人妻精品综合一区二区| 欧美老熟妇乱子伦牲交| 久热这里只有精品99| 天堂俺去俺来也www色官网| 黑人猛操日本美女一级片| 三级国产精品欧美在线观看| 一个人看视频在线观看www免费| 在线观看一区二区三区激情| 欧美 日韩 精品 国产| kizo精华| 国产精品成人在线| av国产精品久久久久影院| 青春草国产在线视频| 成年女人在线观看亚洲视频| 天天躁夜夜躁狠狠久久av| 五月玫瑰六月丁香| 女性被躁到高潮视频| 乱码一卡2卡4卡精品| 日韩熟女老妇一区二区性免费视频| 黄片无遮挡物在线观看| 在线观看美女被高潮喷水网站| 少妇 在线观看| 欧美97在线视频| av.在线天堂| av免费在线看不卡| 免费日韩欧美在线观看| 日韩一区二区视频免费看| 天堂中文最新版在线下载| 国产成人免费观看mmmm| 国产精品久久久久成人av| 在线亚洲精品国产二区图片欧美 | 亚洲无线观看免费| 国产精品一区二区在线观看99| 久久久久精品久久久久真实原创| 高清欧美精品videossex| 午夜福利视频在线观看免费| 免费观看av网站的网址| √禁漫天堂资源中文www| 成年人午夜在线观看视频| 夜夜爽夜夜爽视频| 18在线观看网站| 久久午夜综合久久蜜桃| 午夜精品国产一区二区电影| 全区人妻精品视频| av一本久久久久| 久久人人爽人人片av| 亚洲国产精品999| 亚洲精品久久久久久婷婷小说| 搡老乐熟女国产| 少妇熟女欧美另类| 国产在线免费精品| 国产精品国产三级国产av玫瑰| 夜夜骑夜夜射夜夜干| 欧美精品高潮呻吟av久久| 天堂中文最新版在线下载| 一边亲一边摸免费视频| 啦啦啦啦在线视频资源| 夫妻性生交免费视频一级片| 国产亚洲欧美精品永久| 啦啦啦啦在线视频资源| 欧美激情极品国产一区二区三区 | 久久韩国三级中文字幕| 亚洲第一av免费看| www.av在线官网国产| 精品视频人人做人人爽| 亚洲精品自拍成人| 精品久久久久久久久av| 超碰97精品在线观看| 国产一区二区三区综合在线观看 | 国产黄片视频在线免费观看| av在线app专区| 全区人妻精品视频| 久久久久人妻精品一区果冻| av线在线观看网站| 国产免费福利视频在线观看| 午夜av观看不卡| 亚洲综合色惰| 久久鲁丝午夜福利片| 大话2 男鬼变身卡| 哪个播放器可以免费观看大片| 国国产精品蜜臀av免费| 精品一区在线观看国产| 哪个播放器可以免费观看大片| 精品久久久久久久久亚洲| www.色视频.com| 亚洲婷婷狠狠爱综合网| a 毛片基地| 久久久久国产精品人妻一区二区| 国产日韩欧美视频二区| 新久久久久国产一级毛片| 麻豆成人av视频| 边亲边吃奶的免费视频| 另类亚洲欧美激情| 国产免费视频播放在线视频| 精品人妻一区二区三区麻豆| 天天操日日干夜夜撸| 久久久午夜欧美精品| 国产高清不卡午夜福利| av线在线观看网站| 久久久国产欧美日韩av| a级片在线免费高清观看视频| 两个人的视频大全免费| 日日撸夜夜添| 亚洲精品av麻豆狂野| 午夜91福利影院| 久久久久久久亚洲中文字幕| 成人手机av| 毛片一级片免费看久久久久| 两个人免费观看高清视频| 欧美精品一区二区免费开放| 亚洲精品av麻豆狂野| 在线观看人妻少妇| 大香蕉久久成人网| 99久久精品国产国产毛片| 久久久精品免费免费高清| 欧美亚洲 丝袜 人妻 在线| 国产av国产精品国产| 十八禁高潮呻吟视频| 人妻夜夜爽99麻豆av| 高清av免费在线| 啦啦啦啦在线视频资源| 人妻夜夜爽99麻豆av| 亚洲情色 制服丝袜| 国产深夜福利视频在线观看| 久久99精品国语久久久| 美女福利国产在线| 又黄又爽又刺激的免费视频.| 午夜日本视频在线| 我要看黄色一级片免费的| 亚洲国产欧美日韩在线播放| 国产av精品麻豆| 日韩中文字幕视频在线看片| 精品少妇内射三级| 99久久精品一区二区三区| 街头女战士在线观看网站| 国模一区二区三区四区视频| 欧美激情 高清一区二区三区| 国产精品一国产av| 国产精品女同一区二区软件| 午夜福利影视在线免费观看| 亚洲精品中文字幕在线视频| 国产不卡av网站在线观看| 午夜日本视频在线| 我的女老师完整版在线观看| 性色avwww在线观看| 国产在线一区二区三区精| 久久女婷五月综合色啪小说| 少妇人妻久久综合中文| 99久久综合免费| 欧美成人午夜免费资源| 中文字幕人妻丝袜制服| 精品一区在线观看国产| 亚洲精品中文字幕在线视频| 大香蕉久久网| 婷婷色综合大香蕉| 国产免费视频播放在线视频| 久久久久久人妻| 新久久久久国产一级毛片| 亚洲久久久国产精品| 国产欧美日韩综合在线一区二区| 免费大片18禁| 亚洲精品中文字幕在线视频| 午夜福利,免费看| 欧美日韩一区二区视频在线观看视频在线| 18在线观看网站| 久久ye,这里只有精品| 日韩欧美精品免费久久| 亚洲av不卡在线观看| 伦理电影大哥的女人| 精品国产露脸久久av麻豆| 欧美精品一区二区免费开放| 亚洲精品乱久久久久久| 久久精品国产亚洲av天美| 看十八女毛片水多多多| 亚洲av男天堂| 久久午夜福利片| 三上悠亚av全集在线观看| 人妻 亚洲 视频| 国产成人精品无人区| 啦啦啦在线观看免费高清www| 春色校园在线视频观看| 免费黄频网站在线观看国产| 久久99一区二区三区| 丰满迷人的少妇在线观看| 在线观看人妻少妇| 日韩精品有码人妻一区| 丰满饥渴人妻一区二区三| 插逼视频在线观看| 国产成人a∨麻豆精品| 欧美变态另类bdsm刘玥| 视频区图区小说| 午夜91福利影院| 亚洲图色成人| 熟女av电影| 国产有黄有色有爽视频| 十分钟在线观看高清视频www| 亚洲四区av| 秋霞伦理黄片| 91国产中文字幕| 啦啦啦视频在线资源免费观看| 男女无遮挡免费网站观看| 精品久久久久久久久亚洲| 人妻少妇偷人精品九色| 成人亚洲精品一区在线观看| 青青草视频在线视频观看| 丁香六月天网| 免费播放大片免费观看视频在线观看| 亚洲成人一二三区av| 制服丝袜香蕉在线| 国语对白做爰xxxⅹ性视频网站| a级片在线免费高清观看视频| 爱豆传媒免费全集在线观看| 国产乱人偷精品视频| 黄片无遮挡物在线观看| 熟女人妻精品中文字幕| 午夜福利在线观看免费完整高清在| 亚洲精品久久午夜乱码| 国产午夜精品一二区理论片| 大片电影免费在线观看免费| 少妇精品久久久久久久| 色婷婷av一区二区三区视频| 国产69精品久久久久777片| 成年美女黄网站色视频大全免费 | 国产在线一区二区三区精| 在线观看免费日韩欧美大片 | 免费人妻精品一区二区三区视频| 亚洲高清免费不卡视频| 乱码一卡2卡4卡精品| 免费高清在线观看日韩| 桃花免费在线播放| 一本大道久久a久久精品| 日本欧美视频一区| 在线看a的网站| 国国产精品蜜臀av免费| 欧美3d第一页| 老司机影院毛片| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 亚洲经典国产精华液单| 国产高清有码在线观看视频| 少妇被粗大的猛进出69影院 | 亚洲国产精品一区三区| kizo精华| 国产高清国产精品国产三级| 久热这里只有精品99| 爱豆传媒免费全集在线观看| 丁香六月天网| av免费在线看不卡| 高清在线视频一区二区三区| 天天操日日干夜夜撸| 97在线视频观看| 99热6这里只有精品| 91精品一卡2卡3卡4卡| 亚洲av中文av极速乱| av在线播放精品| 久久久久久久国产电影| 亚洲第一区二区三区不卡| 91在线精品国自产拍蜜月| 久久毛片免费看一区二区三区| xxx大片免费视频| 三上悠亚av全集在线观看| 久久99精品国语久久久| 欧美xxxx性猛交bbbb| 肉色欧美久久久久久久蜜桃| 夜夜骑夜夜射夜夜干| 亚洲经典国产精华液单| 亚洲婷婷狠狠爱综合网| .国产精品久久| 亚洲精品av麻豆狂野| 欧美+日韩+精品| 搡老乐熟女国产| 午夜免费观看性视频| 成人国产av品久久久| 亚洲精品日韩在线中文字幕| 狂野欧美白嫩少妇大欣赏| 日本免费在线观看一区| 我要看黄色一级片免费的| 夫妻性生交免费视频一级片| 一区二区三区免费毛片| 国产成人精品福利久久| 男男h啪啪无遮挡| 伊人久久精品亚洲午夜| av在线播放精品| 欧美激情国产日韩精品一区| 久久人人爽av亚洲精品天堂| 亚洲久久久国产精品| 中国国产av一级| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 赤兔流量卡办理| 日韩av免费高清视频| 国产精品麻豆人妻色哟哟久久| 精品一区二区免费观看| 亚洲少妇的诱惑av| 国产精品秋霞免费鲁丝片| 亚洲精品456在线播放app| 在线观看免费日韩欧美大片 | 人人妻人人澡人人看| 麻豆乱淫一区二区| 亚洲精品,欧美精品| 韩国av在线不卡| 日韩av不卡免费在线播放| 下体分泌物呈黄色| 精品人妻偷拍中文字幕| 国产黄色视频一区二区在线观看| 视频在线观看一区二区三区| 亚洲精品第二区| 男女边摸边吃奶| 黄色视频在线播放观看不卡| 久久热精品热| 久久精品国产亚洲网站| 亚洲美女黄色视频免费看| 国产亚洲一区二区精品| 亚洲天堂av无毛| 亚洲精品456在线播放app| 日本av免费视频播放| 国产精品一二三区在线看| 国产极品天堂在线| 少妇熟女欧美另类| 欧美精品一区二区免费开放| 91精品国产九色| av免费观看日本| 亚洲av综合色区一区| 免费人妻精品一区二区三区视频| 亚洲精品中文字幕在线视频| 26uuu在线亚洲综合色| 欧美成人精品欧美一级黄| 欧美日韩视频精品一区| 在线亚洲精品国产二区图片欧美 | 免费少妇av软件| 777米奇影视久久| 亚洲人与动物交配视频| 亚洲国产精品一区二区三区在线| 国产 精品1| 91久久精品电影网| 亚洲精品第二区| 这个男人来自地球电影免费观看 | 亚洲av.av天堂| 久久精品国产鲁丝片午夜精品| www.色视频.com| av在线app专区| 最近最新中文字幕免费大全7| 丁香六月天网| 大片电影免费在线观看免费| 考比视频在线观看| 免费大片18禁| 五月开心婷婷网| 精品少妇内射三级| 免费观看在线日韩| 欧美精品亚洲一区二区| 免费人妻精品一区二区三区视频| 香蕉精品网在线| 久久精品久久久久久噜噜老黄| 亚洲国产成人一精品久久久| 亚洲国产av新网站| 免费少妇av软件| 精品久久久久久久久av| 国产午夜精品一二区理论片| 综合色丁香网| 久久女婷五月综合色啪小说| 2021少妇久久久久久久久久久| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 制服诱惑二区| 亚洲不卡免费看| 极品人妻少妇av视频| 国产亚洲精品第一综合不卡 | 男人爽女人下面视频在线观看| 欧美精品人与动牲交sv欧美| 天堂俺去俺来也www色官网| 亚洲精品日韩av片在线观看| 成年av动漫网址| 亚洲怡红院男人天堂| 国产亚洲精品第一综合不卡 | 国产 精品1| 久久久午夜欧美精品| 另类精品久久| 成人国产麻豆网| 中文字幕人妻熟人妻熟丝袜美| 国产一区亚洲一区在线观看| 日韩欧美一区视频在线观看| 亚洲av电影在线观看一区二区三区| 亚洲,欧美,日韩| 91在线精品国自产拍蜜月| 久久精品人人爽人人爽视色| av不卡在线播放| 亚洲av欧美aⅴ国产| 国产片特级美女逼逼视频| 日韩人妻高清精品专区| 建设人人有责人人尽责人人享有的| 制服诱惑二区| 新久久久久国产一级毛片| 自拍欧美九色日韩亚洲蝌蚪91| 久热这里只有精品99| 9色porny在线观看| 国产精品一区二区在线观看99| 最近手机中文字幕大全| 日韩一本色道免费dvd| 久久精品久久久久久久性| 一区在线观看完整版| av线在线观看网站| 一边亲一边摸免费视频| 人人妻人人爽人人添夜夜欢视频| 大片免费播放器 马上看| 精品熟女少妇av免费看| 狂野欧美白嫩少妇大欣赏| 街头女战士在线观看网站| 久久午夜福利片| 久久久久久久久久人人人人人人| 老司机亚洲免费影院| 亚洲一区二区三区欧美精品| 中文字幕制服av| 我的女老师完整版在线观看| 考比视频在线观看| 亚洲精品av麻豆狂野| 国产熟女欧美一区二区| 九色亚洲精品在线播放| 美女主播在线视频| 乱码一卡2卡4卡精品| 有码 亚洲区| 欧美+日韩+精品| 99久久综合免费| 欧美xxⅹ黑人| 韩国高清视频一区二区三区| 精品久久久久久电影网| 日韩一区二区视频免费看| 欧美性感艳星| 亚洲人成77777在线视频| 国产老妇伦熟女老妇高清| 人成视频在线观看免费观看| 免费看av在线观看网站| 女人精品久久久久毛片| 久久久欧美国产精品| 免费观看无遮挡的男女| 在线看a的网站| 久久99热6这里只有精品| 少妇人妻 视频| 国产日韩欧美在线精品| 免费大片18禁| 国产有黄有色有爽视频| 麻豆乱淫一区二区| 亚洲国产精品专区欧美| 成人午夜精彩视频在线观看| 久久国产亚洲av麻豆专区| 大话2 男鬼变身卡| 成人二区视频| 免费不卡的大黄色大毛片视频在线观看| 欧美97在线视频| 在线观看美女被高潮喷水网站| 久久国产精品大桥未久av| 欧美日韩av久久| 成人黄色视频免费在线看| 国产乱人偷精品视频| 欧美国产精品一级二级三级| av有码第一页| 亚洲精品日本国产第一区| 国产精品无大码| 91成人精品电影| 中文字幕最新亚洲高清| 国产男女超爽视频在线观看| 精品人妻熟女av久视频| 久久久久精品久久久久真实原创| 一级毛片aaaaaa免费看小| 免费观看性生交大片5| 一级毛片黄色毛片免费观看视频| 麻豆精品久久久久久蜜桃| 成人国产麻豆网| 精品视频人人做人人爽| 国精品久久久久久国模美| 久久99一区二区三区| 欧美亚洲日本最大视频资源| av又黄又爽大尺度在线免费看| kizo精华| 国语对白做爰xxxⅹ性视频网站| 午夜日本视频在线| 久久精品国产鲁丝片午夜精品| 国产精品欧美亚洲77777| 亚洲中文av在线| 五月天丁香电影| 国产一区二区三区av在线| 午夜福利视频在线观看免费| 久久人人爽av亚洲精品天堂| av.在线天堂| 性色av一级| 18禁裸乳无遮挡动漫免费视频| 欧美激情 高清一区二区三区| 午夜91福利影院| 国产极品天堂在线| 日本91视频免费播放| 中国国产av一级| 国产av码专区亚洲av| 久久99热6这里只有精品| 色吧在线观看| 亚洲国产精品999| 夜夜骑夜夜射夜夜干| 国产精品成人在线| 在线观看一区二区三区激情| 在线观看美女被高潮喷水网站| 免费观看在线日韩| 美女视频免费永久观看网站| 精品久久蜜臀av无| 国产欧美亚洲国产| 美女国产高潮福利片在线看| 婷婷色综合大香蕉| 99久久精品一区二区三区| 精品少妇内射三级| 国产精品国产三级专区第一集| 亚洲精品一区蜜桃| 中文字幕最新亚洲高清| 高清午夜精品一区二区三区| 韩国高清视频一区二区三区| 国产一区亚洲一区在线观看| 久久久久国产网址| 简卡轻食公司| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一二三| 亚洲成人手机| 久久久午夜欧美精品| 在线观看人妻少妇| 日本爱情动作片www.在线观看| 在线看a的网站| 亚洲人成网站在线播| 国产亚洲午夜精品一区二区久久| 成人免费观看视频高清| 一级黄片播放器| 色吧在线观看| 伊人久久精品亚洲午夜| a 毛片基地| 少妇人妻久久综合中文| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费播放大片免费观看视频在线观看| 80岁老熟妇乱子伦牲交| 一级毛片我不卡| 国产精品成人在线| 男人添女人高潮全过程视频| 精品久久久精品久久久| 国产精品一国产av| 亚洲精品久久成人aⅴ小说 | 日本黄色日本黄色录像| 亚洲欧美中文字幕日韩二区| 亚洲国产av影院在线观看| 久久亚洲国产成人精品v| 99久久精品国产国产毛片| 3wmmmm亚洲av在线观看| 热99久久久久精品小说推荐| 我要看黄色一级片免费的| 欧美bdsm另类| 赤兔流量卡办理| 插阴视频在线观看视频| 国产一区二区三区av在线| 亚洲丝袜综合中文字幕| 最后的刺客免费高清国语| 黄片无遮挡物在线观看| h视频一区二区三区| 最后的刺客免费高清国语| 亚洲精品日韩av片在线观看| av专区在线播放| videossex国产| 成人黄色视频免费在线看| 久久久精品免费免费高清| 亚洲精品国产av成人精品| av免费观看日本| 久久久欧美国产精品| 日本欧美国产在线视频| 搡女人真爽免费视频火全软件| 久久av网站| 黄色视频在线播放观看不卡| 亚洲av.av天堂| 国产欧美亚洲国产| 国产老妇伦熟女老妇高清| 免费观看a级毛片全部| 国产精品麻豆人妻色哟哟久久| 制服人妻中文乱码| 99视频精品全部免费 在线| 日本av免费视频播放| 丝瓜视频免费看黄片| 在线观看三级黄色|