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

    基本解法求解反問題的正則化方法

    2015-04-18 07:18:30張耀明周愛華潘月君
    關(guān)鍵詞:源點位勢病態(tài)

    楊 振,張耀明, 周愛華,潘月君

    (山東理工大學 理學院,山東 淄博 255049)

    基本解法求解反問題的正則化方法

    楊 振,張耀明, 周愛華,潘月君

    (山東理工大學 理學院,山東 淄博 255049)

    針對基本解法在求解反問題時的病態(tài)特性,將截斷奇異值分解(TSVD)、Tikhonov正則化方法應用于所得病態(tài)系統(tǒng)方程的求解,采用L曲線法和GCV方法確定其正則參數(shù),并比較了4種組合方法求解的精確性和穩(wěn)定性.數(shù)值算例表明,TSVD方法、Tikhonov正則化方法結(jié)合L曲線法和GCV法可有效地處理反問題中的病態(tài)特性.

    基本解法;反問題;截斷奇異值法;Tikhonov法;L曲線;GCV

    數(shù)學物理反問題的研究興起于20世紀50年代前后,其研究的主要對象是與勘測、識別和設(shè)計等有關(guān)的應用問題.在實際工程應用中,有時由于條件所限,僅僅可以獲得部分邊界上的邊界條件,其余邊界的邊界條件均不可測量.從數(shù)學的角度來看,這屬于反問題的一種.這類問題由于具有不適定及病態(tài)特性,故其數(shù)值分析具有相當大的挑戰(zhàn)性[1-3].邊界元法(BEM)可有效地求解此類問題,但涉及奇異和幾乎奇異積分的處理[4-6],計算過程相當?shù)胤爆嵑秃臅r.基本解法(MFS)無需對區(qū)域及其邊界進行網(wǎng)格劃分,因此無需計算單元積分,具有計算精度高、收斂速度快、程序?qū)崿F(xiàn)簡單等優(yōu)良特性,如今已被廣泛應用于各種問題的求解,但都不可避免地涉及到病態(tài)線性方程組的處理[6-11].

    為了克服基本解法在求解反問題線性系統(tǒng)時的病態(tài)性,許多學者提出和發(fā)展了各種方法,其中最具影響的是正則化方法.Chen C S[8]分別用TSVD方法和高斯消去法來求解基本解法中的病態(tài)系統(tǒng),認為TSVD法比高斯消去法優(yōu)越.此后, Feng G R等[9]進一步考慮了一般的邊界條件和邊界幾何形狀的問題,得出TSVD法在處理病態(tài)問題時具有計算穩(wěn)定、精度高的特點.而其他的正則化方法如Tikhonov正則化方法、Landweber迭代法和共軛梯度法等[1-2,11]也有一些應用.考慮到目前仍沒有一種適合所有病態(tài)問題解算的最優(yōu)正則化方法,因此對不同的正則化方法進行比較研究是一項具有重要意義的工作.

    當我們選用一種直接正則化方法后,正則化參數(shù)的確定非常重要,它是改善病態(tài)問題解的精確性的關(guān)鍵.目前,正則化參數(shù)的選取準則有L曲線法、廣義交叉檢驗(GCV)準則等[1-2,11].本文以二維位勢反問題為例,將TSVD、Tikhonov兩種直接正則化方法與L曲線和GCV兩種正則參數(shù)選取方法組合而成的4種不同方法,應用于基本解法求解未知邊界條件信息.通過數(shù)值算例,分析比較不同組合方法的計算效率和精度,以期為實際應用提供合理的選擇依據(jù).

    1 二維位勢反問題的基本解法

    本文假定Ω是R2中的一個有界區(qū)域,Γ=?Ω是其邊界.n=(n1,n2)是區(qū)域Ω的邊界Γ在x點處的單位外法向量.

    1.1 二維位勢邊界條件識別反問題

    二維位勢問題的控制微分方程為

    邊界條件為

    式中:u為勢函數(shù);n為邊界外法向量;Γu, Γq分別是已知u和?u/?n的邊界.

    二維位勢問題控制方程的基本解為

    1.2 二維位勢反問題的基本解法

    基本解法(MFS)是一種無網(wǎng)格法.MFS的基本思想是在邊界之外一定距離處放置一定數(shù)量的虛擬源點,并設(shè)想在每個源點處存在一個虛擬的密度函數(shù)值,其在研究區(qū)域上某點產(chǎn)生的影響能夠表示成密度函數(shù)與原問題基本解的線性組合.由于源點分布在物理邊界之外的虛擬邊界上,從而避免了基本解的源點奇異性.目前該方法已被應用于各種領(lǐng)域,并取得了很好的效果.

    假設(shè)在虛擬邊界?!渖戏植糿個源點yj,j=1,2,…,n,若在邊界Γ1上選取m個配置點xi,i=(1,2,…,m).則邊界配置點處的位勢和法向梯度可表示為:

    (1)

    (2)

    這里Φ=(Φ1,Φ2,…,Φn)表示源點的密度函數(shù)值向量,u*(xi,yj)為問題的基本解.

    根據(jù)式(1)、式(2)以及已知的邊界條件可得如下一個線性方程組

    (3)

    由方程(3)可以計算出Φ=(Φ1,Φ2,…,Φn),則由此可通過以下方程計算出未知邊界上的位勢和通量.

    (4)

    (5)

    一般地,在求解邊界條件識別反問題時,方程組(3)往往是不適定的或者高度病態(tài)的,常規(guī)方法如高斯消去法很難求得準確解.因此在求解方程時,我們采用正則化方法求解線性方程組(3).

    1.3 正則化方法

    1.3.1 截斷奇異值分解(TSVD)

    假設(shè)所要求解的方程為

    Ax=b

    (6)

    其中A∈RM×N(M>N),x∈RN,b∈RM,則矩陣A的奇異值分解形式為

    (7)

    這個解法用K階矩陣Ak逼近N階矩陣A.下面給出它的形式

    (8)

    其中K是TSVD方法的截斷項數(shù),Σk=diag(σ1,σ2,…,σK,0,…,0),通過這種方式,方程(6)中的矩陣A被AK代替,得到下列方程:

    AKx=b

    (9)

    方程(9)的TSVD解可表示為xK,且

    (10)

    這里的濾波因子為

    1.3.2Tikhonov正則化方法(TR)

    Tikhonov正則化方法是一種非常有效而且普遍使用的求解病態(tài)問題的方法,其求解過程是找到一組適定問題的解,使得這組解與原不適定問題比較接近,然后用這組適定問題的解去逼近原不適定問題的解.Tikhonov正則化方法是把正則化泛函

    Jα(x)=‖Ax-b‖2+α2‖x‖2,α>0

    (11)

    的極小元xα作為方程(6)的正則化解.表示成如下形式:

    (12)

    1.4 正則化參數(shù)的確定

    1.4.1L曲線準則

    L曲線準則的主要思想是通過log-log尺度以反映‖x‖2與‖Ax-b‖2的曲線對比,然后根據(jù)對比得出的結(jié)果決定正則化參數(shù).因為曲線形狀和“L”的形狀通常比較接近, 所以它被稱為L曲線準則.

    在使用L曲線準則確定正則化參數(shù)時,若運用Tikhonov正則化方法,則L曲線是由α>0的所有點(‖Axα-b‖,‖xα‖)構(gòu)成的光滑曲線

    (13)

    若運用TSVD法,則L曲線是由一系列點(‖AxK-b‖,‖xK‖)組成的插值樣條曲線.

    (14)

    對于確定的α>0(或者K=1,2,…,N),我們解方程(10)(或者(12))就會得到解xα(或者xK),以余量?!珹xα-b‖(或者‖AxK-b‖)為橫坐標,解的?!瑇α‖(或者‖xK‖)為縱坐標,可畫出L曲線.找出拐點處的α或者K,即為要找的參數(shù)值.

    1.4.2 廣義交叉檢驗(GCV)準則

    GCV準則是用正則化參數(shù)α(或者K)作為參變量,來求解GCV函數(shù)的最小值,只要得到GCV函數(shù)的極小值,那么對應α(或者K)就是我們要確定的最優(yōu)正則化參數(shù).其計算公式為

    (15)

    或者

    (16)

    2 數(shù)值算例

    考慮兩個二維位勢反問題的數(shù)值算例,采用4種不同的正則化組合方法對算例進行求解.為了對比不同方法數(shù)值解的準確性,定義如下平均相對誤差

    (17)

    例1 如圖1所示方形閉域熱流問題,邊長為2.精確解為u=y2-x2,邊界條件為:邊AB、BC上的溫度和熱流已知,邊CD、DA上所有邊界條件未知.

    如圖1所示,在距離真實邊界距離為d的地方取一虛擬邊界Γ′,使其與熱流區(qū)域邊界相似.在虛擬邊界四邊各分布10個源點,共40個源點,如圖2所示.為了考察正則化方法對于位勢反問題求解的有效性,我們分別用4種組合方法對本例進行求解.我們?nèi)≡袋c分布在虛擬邊界?!渖?,選取不同的虛實邊界距離d分別為5和10.圖3~圖 6給出了d取不同值時,未知邊界上各點溫度和熱流的相對誤差,可以看出,本文給出的4種組合方法在未知邊界各點求得的數(shù)值解都獲得了較高精度.

    圖1 方形閉域熱流 圖2 節(jié)點編號

    圖3 未知邊界的溫度誤差(d=5)

    圖4 未知邊界的熱流誤差(d=5)

    圖5 未知邊界的溫度誤差(d=10)

    圖6 未知邊界的熱流誤差(d=10)

    例2 如圖7所示圓形閉域熱流問題,半徑為1.精確解為u=x2-y2,我們考慮的邊界條件為:1/4邊界上的溫度和熱流已知,另外3/4邊界上的邊界條件未知.此時我們所得方程是不適定的,因此常規(guī)方法很難求解.

    如圖7所示,在距離真實邊界距離為d的地方取一虛擬邊界?!?,使其與熱流區(qū)域邊界相似.并在虛擬邊界上均勻地布置100個源點.同樣我們分別用4種組合方法對本例進行求解.選取不同的虛實邊界距離d分別為5和10.表1、表2給出了d取不同值時,未知邊界上溫度和熱流的平均相對誤差.表1和表2更準確地反映出本文給出的4種組合方法計算結(jié)果差異,也不難看出,MFS結(jié)合4種組合方法都能夠有效地求解不同邊界幾何形狀的位勢邊界條件反問題.

    圖7 圓形閉域熱流

    表1 四種方法參數(shù)選取及平均相對誤差比較結(jié)果(d=5)

    正則化方法參數(shù)(α/K)平均相對誤差/%溫度熱流TSVD-LC173.400991×10-61.709062×10-6TSVD-GCV122.492231×10-53.790088×10-5TR-LC2.44×10-113.268636×10-61.602206×10-6TR-GCV2.01×10-68.815173×10-51.121432×10-4

    表2 四種方法參數(shù)選取及平均相對誤差比較結(jié)果(d=10)

    正則化方法參數(shù)(α/K)平均相對誤差/%溫度熱流TSVD-LC131.946931×10-79.724293×10-8TSVD-GCV122.417499×10-73.138218×10-7TR-LC3.37×10-111.577152×10-76.523213×10-8TR-GCV3.91×10-111.866861×10-79.110862×10-8

    3 結(jié)束語

    本文針對基本解法求解反問題時的病態(tài)特性,應用TSVD方法和Tikhonov正則化方法處理所得病態(tài)系統(tǒng)方程,L曲線法和GCV法確定其正則化參數(shù),并比較了不同正則化方法的計算精度.數(shù)值算例表明,針對不同邊界條件以及邊界幾何形狀,MFS結(jié)合TSVD方法和Tikhonov正則化方法都能夠有效求解二維位勢反問題,即使已知邊界信息非常有限,依然能夠獲得穩(wěn)定、高精度的數(shù)值解.

    [1]肖庭延. 反問題的數(shù)值解法[M]. 北京:科學出版社, 2003.

    [2]劉繼軍. 不適定問題的正則化方法及應用[M]. 北京:科學出版社, 2005.

    [3]HansenPC.Rank-deficientanddiscreteill-poseedproblems:numericalaspectsoflinearinversion[M].Philadelphia:SIAM,1998.

    [4]張耀明,谷巖,陳正宗. 位勢邊界元法中的邊界層效應與薄體結(jié)構(gòu)[J]. 力學學報,2010,42(2):219-227.

    [5]ZhangYM,SunCL.Ageneralalgorithmforthenumericalevaluationofnearlysingularboundaryintegralsintheequivalentnon-singularBIEswithindirectunknowns[J].JournaloftheChineseInstituteofEngineers, 2008, 31(3):437-447.

    [6]孫煥純. 無奇異邊界元法[M]. 大連:大連理工大學出版社, 1999.

    [7]KupradzeVD,AleksidzeMA.Themethodoffunctionalequationsfortheapproximatesolutionofcertainboundaryvalueproblems[J].USSRComput.MathPhy,1964,4:82-126.

    [8]ChenCS,HokwonA.Somecommentsontheill-conditioningofthemethodoffundamentalsolutions[J].EngineeringAnalysiswithBoundaryElements,2006,30:405-410

    [9]FengGR,LiM,ChenCS.Ontheill-conditioningoftheMFSforirregularboundarydatawithsufficientregularity[J].EngineeringAnalysiswithBoundaryElements,2014,41:98-102.

    [10]TikhonovAN,ArseninVY.SolutionsofILL-posedproblem[M].NewYork:JohnWileyandSons,1977.

    [11]WeiT,HonYC,LingLV.MethodoffundamentalsolutionswithregularizationtechniquesforCauchyproblemsofellipticoperators[J].EngineeringAnalysiswithBoundaryElements,2007,31:373-385.

    (編輯:郝秀清)

    The regularization method of the MFS for the inverse problem

    YANG Zhen, ZHANG Yao-ming, ZHOU Ai-hua, PAN Yue-jun

    (School of Science, Shandong University of Technology, Zibo 255049, China)

    In order to resolve ill-conditioned problems existing in the regularization MFS for the 2D boundary conditions identification potential problems, suitable regularization methods are needed. The truncated singular value decomposition(TSVD)method and Tikhonov regularization method are used to solve linear systems with a large number of conditions respectively. The L-curve and generalized cross validation (GCV) methods are employed to determine the optimal regularization parameters. Furthermore, the accuracy and robustness of regularization solution for four combined methods are investigated. Numerical results show that TSVD method and Tikhonov method can effectively solve the ill posed system caused by inverse problems. Through applying to L-curve method and GCV method, continuous regularization parameter for Tikhonov method can be confirmed reasonably.

    MFS; inverse problems; truncated singular value decomposition; Tikhonov method;L-curve; GCV

    2015-01-11

    山東省自然科學基金重點項目(ZR2010AZ003)

    楊振, 男, librayz@126.com; 通信作者:張耀明,男,zymfc@163.com

    1672-6197(2015)06-0020-05

    O342

    A

    猜你喜歡
    源點位勢病態(tài)
    含Hardy位勢的非線性Schr?dinger-Poisson方程正規(guī)化解的多重性
    一類帶強制位勢的p-Laplace特征值問題
    病態(tài)肥胖對門診全關(guān)節(jié)置換術(shù)一夜留院和早期并發(fā)癥的影響
    病態(tài)肥胖對門診關(guān)節(jié)置換術(shù)留夜觀察和早期并發(fā)癥的影響
    君子之道:能移而相天——王夫之《莊子解》對“社會病態(tài)”的氣論診療
    哲學評論(2018年1期)2018-09-14 02:34:46
    隱喻的語篇銜接模式
    外語學刊(2017年3期)2017-12-07 01:45:38
    含變號位勢的ρ-Kirchhoff型方程組無窮多個高能量解的存在性
    含位勢的非線性雙調(diào)和方程解的存在性
    首屆“絲路源點·青年學者研討會”主題論壇在我校成功舉辦
    淺析井控坐崗的源點
    国产午夜精品久久久久久| 久久香蕉激情| 欧美大码av| 欧美人与性动交α欧美软件| 丁香六月天网| 国产精品av久久久久免费| 欧美精品一区二区免费开放| 成年女人毛片免费观看观看9 | 大香蕉久久成人网| 多毛熟女@视频| 色播在线永久视频| 久久久久久人人人人人| 亚洲第一青青草原| 亚洲av日韩在线播放| 制服诱惑二区| 脱女人内裤的视频| 国产精品一区二区精品视频观看| 久久久国产欧美日韩av| 天天添夜夜摸| 首页视频小说图片口味搜索| 五月天丁香电影| 日韩制服丝袜自拍偷拍| 国产精品自产拍在线观看55亚洲 | 91老司机精品| 亚洲精品中文字幕在线视频| a级片在线免费高清观看视频| 日本a在线网址| 在线av久久热| 国产精品国产高清国产av | 欧美激情极品国产一区二区三区| 欧美成狂野欧美在线观看| 一级毛片电影观看| 久热爱精品视频在线9| 国产又色又爽无遮挡免费看| 久久午夜亚洲精品久久| 中文字幕最新亚洲高清| 少妇 在线观看| 99re在线观看精品视频| 麻豆av在线久日| 国产激情久久老熟女| 757午夜福利合集在线观看| tube8黄色片| 久久精品人人爽人人爽视色| 亚洲国产av新网站| 精品久久蜜臀av无| 中文字幕最新亚洲高清| 黑人欧美特级aaaaaa片| 99精品在免费线老司机午夜| 色播在线永久视频| 岛国在线观看网站| 成人免费观看视频高清| 亚洲欧美激情在线| 在线亚洲精品国产二区图片欧美| 久久精品国产a三级三级三级| 日韩成人在线观看一区二区三区| 日韩欧美一区视频在线观看| 免费在线观看完整版高清| 亚洲精品美女久久av网站| av有码第一页| 香蕉丝袜av| 久久亚洲真实| 青草久久国产| 后天国语完整版免费观看| 三级毛片av免费| 亚洲国产欧美在线一区| 人人妻人人澡人人看| 国产精品 欧美亚洲| av免费在线观看网站| 男女无遮挡免费网站观看| 欧美日韩中文字幕国产精品一区二区三区 | 啦啦啦中文免费视频观看日本| 国产精品久久久久成人av| 在线观看66精品国产| 欧美成人午夜精品| 国产av精品麻豆| 在线天堂中文资源库| 国产伦人伦偷精品视频| 久9热在线精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 青青草视频在线视频观看| 黄色视频在线播放观看不卡| 国产深夜福利视频在线观看| 国产伦理片在线播放av一区| 精品一区二区三区视频在线观看免费 | 国产深夜福利视频在线观看| 一区福利在线观看| 在线永久观看黄色视频| videosex国产| 欧美激情极品国产一区二区三区| 久久中文字幕人妻熟女| 亚洲精华国产精华精| 亚洲天堂av无毛| av片东京热男人的天堂| 多毛熟女@视频| 国产精品一区二区免费欧美| 黄色毛片三级朝国网站| 变态另类成人亚洲欧美熟女 | 一进一出抽搐动态| 国产精品免费大片| 高清欧美精品videossex| 久久久国产欧美日韩av| 天天操日日干夜夜撸| 亚洲国产中文字幕在线视频| 亚洲精品国产区一区二| 久久人妻福利社区极品人妻图片| 老司机福利观看| 9热在线视频观看99| h视频一区二区三区| 男女午夜视频在线观看| 国产成人av激情在线播放| 91九色精品人成在线观看| 一级黄色大片毛片| 国产成人系列免费观看| 夜夜爽天天搞| 久久久国产成人免费| 一本—道久久a久久精品蜜桃钙片| 青草久久国产| 精品午夜福利视频在线观看一区 | 国产日韩一区二区三区精品不卡| 国产精品一区二区在线不卡| 国产精品一区二区在线不卡| 久久 成人 亚洲| 水蜜桃什么品种好| xxxhd国产人妻xxx| 在线 av 中文字幕| 精品高清国产在线一区| 亚洲va日本ⅴa欧美va伊人久久| 一区二区三区激情视频| 亚洲九九香蕉| 麻豆成人av在线观看| 精品亚洲成国产av| 久久精品成人免费网站| 黑人猛操日本美女一级片| 国产精品一区二区精品视频观看| 99在线人妻在线中文字幕 | 在线播放国产精品三级| av线在线观看网站| 久久久久久久久久久久大奶| 一区二区三区国产精品乱码| 久久久国产欧美日韩av| 99riav亚洲国产免费| 久久99热这里只频精品6学生| 日韩免费av在线播放| 婷婷丁香在线五月| 香蕉久久夜色| 精品少妇黑人巨大在线播放| 无人区码免费观看不卡 | 一夜夜www| 久久久久精品国产欧美久久久| 国产精品亚洲一级av第二区| 亚洲精品美女久久av网站| 国产麻豆69| 欧美激情 高清一区二区三区| av超薄肉色丝袜交足视频| 国产一区二区 视频在线| 欧美+亚洲+日韩+国产| 亚洲九九香蕉| 怎么达到女性高潮| 国产精品麻豆人妻色哟哟久久| 99国产精品一区二区蜜桃av | 热99国产精品久久久久久7| 国产真人三级小视频在线观看| 最近最新中文字幕大全免费视频| 日韩有码中文字幕| 狂野欧美激情性xxxx| www.精华液| 91老司机精品| 久久精品成人免费网站| 久久久久久久精品吃奶| 国产高清国产精品国产三级| 女人被躁到高潮嗷嗷叫费观| 久久中文字幕人妻熟女| 日本av手机在线免费观看| 国产野战对白在线观看| 麻豆国产av国片精品| 免费在线观看视频国产中文字幕亚洲| 午夜精品国产一区二区电影| 亚洲精品中文字幕一二三四区 | 亚洲,欧美精品.| 99精品久久久久人妻精品| 国产精品久久久久成人av| 80岁老熟妇乱子伦牲交| 亚洲天堂av无毛| 9191精品国产免费久久| 一区二区日韩欧美中文字幕| 久久久国产成人免费| 夜夜骑夜夜射夜夜干| 国产真人三级小视频在线观看| 久久久久国内视频| 国产精品久久久久成人av| 久久九九热精品免费| 一级片'在线观看视频| 精品福利观看| 国产亚洲欧美在线一区二区| 久久热在线av| 国产xxxxx性猛交| 欧美精品一区二区免费开放| 免费日韩欧美在线观看| 亚洲成av片中文字幕在线观看| 99热国产这里只有精品6| 99九九在线精品视频| 日韩欧美一区视频在线观看| 性少妇av在线| 黄网站色视频无遮挡免费观看| 正在播放国产对白刺激| 久久久久久久大尺度免费视频| 老司机影院毛片| 女人精品久久久久毛片| 国产99久久九九免费精品| 99国产精品一区二区三区| 欧美日韩亚洲高清精品| 大陆偷拍与自拍| 啪啪无遮挡十八禁网站| 夜夜夜夜夜久久久久| 王馨瑶露胸无遮挡在线观看| 在线观看免费视频日本深夜| 啦啦啦免费观看视频1| av福利片在线| av福利片在线| 久久国产精品人妻蜜桃| 少妇 在线观看| 久久av网站| 久久久精品国产亚洲av高清涩受| 亚洲欧美日韩另类电影网站| 91av网站免费观看| 黑人欧美特级aaaaaa片| 国产主播在线观看一区二区| 无人区码免费观看不卡 | 中亚洲国语对白在线视频| 淫妇啪啪啪对白视频| 波多野结衣av一区二区av| 欧美乱妇无乱码| 757午夜福利合集在线观看| 在线天堂中文资源库| 国产激情久久老熟女| 成人国产一区最新在线观看| 亚洲一区二区三区欧美精品| 亚洲专区字幕在线| 91成年电影在线观看| 欧美国产精品一级二级三级| 亚洲成人免费av在线播放| 亚洲五月婷婷丁香| 99国产精品一区二区三区| 国产欧美亚洲国产| 热99re8久久精品国产| videos熟女内射| 国产一卡二卡三卡精品| 日本wwww免费看| 一边摸一边抽搐一进一出视频| 精品少妇黑人巨大在线播放| 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 三级毛片av免费| 中文亚洲av片在线观看爽 | 欧美日韩成人在线一区二区| 国产精品 欧美亚洲| 无限看片的www在线观看| 老汉色∧v一级毛片| 久久久久网色| 老司机午夜福利在线观看视频 | 亚洲欧美精品综合一区二区三区| 在线观看人妻少妇| 国产在线免费精品| 中文字幕人妻熟女乱码| 91老司机精品| 悠悠久久av| 人妻一区二区av| 亚洲欧美一区二区三区黑人| 国产精品久久电影中文字幕 | 90打野战视频偷拍视频| 日日爽夜夜爽网站| 亚洲精品美女久久久久99蜜臀| 亚洲成国产人片在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美日韩视频精品一区| 一边摸一边抽搐一进一出视频| 国产在线免费精品| 国产精品一区二区免费欧美| 亚洲成人国产一区在线观看| 久久久精品94久久精品| 一区福利在线观看| 飞空精品影院首页| 国产亚洲午夜精品一区二区久久| 欧美日韩国产mv在线观看视频| 一进一出好大好爽视频| 久久久久国产一级毛片高清牌| 不卡一级毛片| 欧美精品一区二区免费开放| 在线 av 中文字幕| 俄罗斯特黄特色一大片| 欧美日韩国产mv在线观看视频| 美女高潮喷水抽搐中文字幕| 国产男女内射视频| 搡老岳熟女国产| bbb黄色大片| 成人特级黄色片久久久久久久 | 色视频在线一区二区三区| 国产xxxxx性猛交| 国产欧美日韩综合在线一区二区| 色播在线永久视频| 少妇精品久久久久久久| 精品人妻熟女毛片av久久网站| av网站免费在线观看视频| 日韩欧美免费精品| 91字幕亚洲| 亚洲一码二码三码区别大吗| 免费久久久久久久精品成人欧美视频| 无限看片的www在线观看| 国产精品久久久人人做人人爽| 一区二区三区精品91| 久久人人爽av亚洲精品天堂| 在线亚洲精品国产二区图片欧美| 午夜福利免费观看在线| 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 久久久水蜜桃国产精品网| 伦理电影免费视频| 日韩欧美一区视频在线观看| 久久久久网色| 国产欧美日韩精品亚洲av| 国产高清国产精品国产三级| 久久毛片免费看一区二区三区| 亚洲伊人久久精品综合| 999久久久国产精品视频| 欧美日韩视频精品一区| 亚洲精品一卡2卡三卡4卡5卡| 丝袜美腿诱惑在线| 中文字幕精品免费在线观看视频| 午夜福利免费观看在线| 国精品久久久久久国模美| 久久久久国内视频| 在线观看免费视频日本深夜| 久久久久久久久久久久大奶| 亚洲欧美精品综合一区二区三区| 精品免费久久久久久久清纯 | 亚洲精品在线美女| 香蕉国产在线看| 亚洲中文字幕日韩| 久久午夜亚洲精品久久| 免费在线观看影片大全网站| 欧美日韩福利视频一区二区| 纯流量卡能插随身wifi吗| 欧美亚洲 丝袜 人妻 在线| 国产精品熟女久久久久浪| 午夜福利在线免费观看网站| 91精品国产国语对白视频| 人妻一区二区av| 久久久精品94久久精品| 伊人久久大香线蕉亚洲五| 涩涩av久久男人的天堂| 91精品三级在线观看| 午夜老司机福利片| 国产成+人综合+亚洲专区| 精品国产一区二区三区四区第35| 欧美日韩亚洲综合一区二区三区_| 精品午夜福利视频在线观看一区 | 久久精品熟女亚洲av麻豆精品| av片东京热男人的天堂| 午夜91福利影院| 欧美变态另类bdsm刘玥| 99热国产这里只有精品6| 在线 av 中文字幕| 久久国产精品大桥未久av| 亚洲av片天天在线观看| 欧美激情高清一区二区三区| 女性被躁到高潮视频| 黑丝袜美女国产一区| 少妇猛男粗大的猛烈进出视频| 日韩欧美三级三区| 精品国产乱子伦一区二区三区| av视频免费观看在线观看| 亚洲成国产人片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 大片电影免费在线观看免费| 国产不卡一卡二| 国产单亲对白刺激| 啦啦啦免费观看视频1| 成人三级做爰电影| 亚洲国产欧美一区二区综合| 亚洲专区字幕在线| 男女下面插进去视频免费观看| 国精品久久久久久国模美| 国产精品二区激情视频| 免费一级毛片在线播放高清视频 | 国产熟女午夜一区二区三区| 悠悠久久av| 日韩有码中文字幕| 汤姆久久久久久久影院中文字幕| 久久精品亚洲熟妇少妇任你| 国产精品 国内视频| 午夜激情av网站| 国产精品 国内视频| 一级黄色大片毛片| 成人黄色视频免费在线看| 亚洲精品在线美女| bbb黄色大片| 亚洲精品在线观看二区| 老司机午夜福利在线观看视频 | 国产精品美女特级片免费视频播放器 | 一本久久精品| 天天影视国产精品| 午夜福利在线观看吧| 黄片小视频在线播放| 另类亚洲欧美激情| 国产精品成人在线| 中文字幕精品免费在线观看视频| 久久av网站| 在线av久久热| 在线亚洲精品国产二区图片欧美| 国产成人精品无人区| 日日摸夜夜添夜夜添小说| 日韩人妻精品一区2区三区| 怎么达到女性高潮| 午夜福利一区二区在线看| 成人国语在线视频| 欧美日韩亚洲高清精品| av超薄肉色丝袜交足视频| 侵犯人妻中文字幕一二三四区| 最近最新免费中文字幕在线| 不卡一级毛片| 精品久久久久久电影网| 97在线人人人人妻| 最新在线观看一区二区三区| 男女无遮挡免费网站观看| 高清av免费在线| 成年动漫av网址| 久久ye,这里只有精品| 宅男免费午夜| 91精品国产国语对白视频| 国产免费现黄频在线看| 老司机深夜福利视频在线观看| 国产亚洲精品久久久久5区| 午夜福利欧美成人| 色婷婷av一区二区三区视频| 嫩草影视91久久| 国产精品亚洲av一区麻豆| 日韩精品免费视频一区二区三区| 国产成人啪精品午夜网站| 亚洲,欧美精品.| 亚洲精品av麻豆狂野| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美日韩一区二区精品| 亚洲va日本ⅴa欧美va伊人久久| 黄片小视频在线播放| 免费av中文字幕在线| 涩涩av久久男人的天堂| 久久午夜亚洲精品久久| 中文亚洲av片在线观看爽 | 免费在线观看完整版高清| 精品欧美一区二区三区在线| 女人高潮潮喷娇喘18禁视频| 免费人妻精品一区二区三区视频| 一本综合久久免费| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜添小说| 欧美激情 高清一区二区三区| 狂野欧美激情性xxxx| 午夜免费成人在线视频| 精品熟女少妇八av免费久了| 一区二区三区激情视频| 亚洲人成电影观看| 美女国产高潮福利片在线看| 成人亚洲精品一区在线观看| 人人妻人人添人人爽欧美一区卜| 免费少妇av软件| 男人操女人黄网站| 国产精品亚洲一级av第二区| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 亚洲成av片中文字幕在线观看| 美女扒开内裤让男人捅视频| 亚洲欧美色中文字幕在线| 日本撒尿小便嘘嘘汇集6| 日韩制服丝袜自拍偷拍| 久久午夜亚洲精品久久| 国产视频一区二区在线看| 久久久久久亚洲精品国产蜜桃av| 最新在线观看一区二区三区| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕大全免费视频| 91老司机精品| 亚洲第一青青草原| 色尼玛亚洲综合影院| avwww免费| 国产单亲对白刺激| tube8黄色片| 欧美人与性动交α欧美精品济南到| 97人妻天天添夜夜摸| 男女免费视频国产| 成人永久免费在线观看视频 | 亚洲人成电影观看| 啦啦啦中文免费视频观看日本| 久久精品熟女亚洲av麻豆精品| 亚洲免费av在线视频| 在线观看www视频免费| 国产亚洲精品久久久久5区| 色视频在线一区二区三区| 亚洲国产毛片av蜜桃av| 日本vs欧美在线观看视频| 极品少妇高潮喷水抽搐| 老司机午夜十八禁免费视频| 国产亚洲精品第一综合不卡| 老司机午夜福利在线观看视频 | 1024视频免费在线观看| 日韩一区二区三区影片| 一区二区三区精品91| 色婷婷久久久亚洲欧美| 99精品欧美一区二区三区四区| 最新在线观看一区二区三区| 不卡一级毛片| 夜夜爽天天搞| 亚洲午夜理论影院| 精品少妇黑人巨大在线播放| 欧美一级毛片孕妇| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 欧美在线黄色| 一级片'在线观看视频| 亚洲精品中文字幕在线视频| 精品国产一区二区三区久久久樱花| 性色av乱码一区二区三区2| 人人妻人人澡人人看| 国产成+人综合+亚洲专区| 人妻 亚洲 视频| 午夜福利在线观看吧| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| av一本久久久久| 精品少妇一区二区三区视频日本电影| 黄片小视频在线播放| 亚洲国产av新网站| 一本久久精品| 欧美老熟妇乱子伦牲交| 亚洲天堂av无毛| 亚洲国产欧美网| av在线播放免费不卡| 后天国语完整版免费观看| 国产在视频线精品| 天天影视国产精品| 成年人黄色毛片网站| 9191精品国产免费久久| 日韩欧美一区视频在线观看| 美女主播在线视频| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 欧美日韩福利视频一区二区| 久久精品亚洲熟妇少妇任你| 99国产精品免费福利视频| 精品国产一区二区久久| 国产精品98久久久久久宅男小说| 国产一卡二卡三卡精品| e午夜精品久久久久久久| 日韩一区二区三区影片| 亚洲九九香蕉| 精品一区二区三区av网在线观看 | 国产精品久久久av美女十八| 日韩大码丰满熟妇| 国产在视频线精品| 国产不卡av网站在线观看| 精品国产一区二区久久| 麻豆成人av在线观看| 俄罗斯特黄特色一大片| 国产在线精品亚洲第一网站| 国产亚洲精品第一综合不卡| 亚洲五月婷婷丁香| 国产精品熟女久久久久浪| 亚洲av美国av| 老司机福利观看| 国产深夜福利视频在线观看| 天天操日日干夜夜撸| 中文字幕高清在线视频| 他把我摸到了高潮在线观看 | 国产午夜精品久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 80岁老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| 欧美日韩av久久| 亚洲性夜色夜夜综合| 色视频在线一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人免费| 少妇 在线观看| 亚洲人成77777在线视频| 国产aⅴ精品一区二区三区波| 国产亚洲欧美在线一区二区| av片东京热男人的天堂| 亚洲精品中文字幕在线视频| 久久亚洲真实| 不卡av一区二区三区| 久久精品人人爽人人爽视色| 欧美久久黑人一区二区| 精品视频人人做人人爽| 亚洲色图 男人天堂 中文字幕| www日本在线高清视频| 国产男女超爽视频在线观看| 午夜两性在线视频| 男女免费视频国产| 久久 成人 亚洲| 蜜桃在线观看..| 成人国语在线视频| 丝瓜视频免费看黄片| 国产精品偷伦视频观看了| 久久精品亚洲熟妇少妇任你| 亚洲人成77777在线视频| 国产99久久九九免费精品| 欧美黄色片欧美黄色片| 欧美国产精品一级二级三级| 在线观看一区二区三区激情| 久久香蕉激情| 另类精品久久| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| 俄罗斯特黄特色一大片| 色综合婷婷激情|