• <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网站的网址| 久久精品91无色码中文字幕| 中国美女看黄片| av福利片在线| 在线观看免费午夜福利视频| bbb黄色大片| 欧美一级毛片孕妇| 黄色怎么调成土黄色| 美国免费a级毛片| 高清欧美精品videossex| 日韩欧美国产一区二区入口| 人人妻人人澡人人爽人人夜夜| 后天国语完整版免费观看| 最新美女视频免费是黄的| 手机成人av网站| 少妇精品久久久久久久| 69av精品久久久久久 | 国产高清视频在线播放一区| 久久久水蜜桃国产精品网| 亚洲av国产av综合av卡| 国产成人av激情在线播放| 一级黄色大片毛片| 日韩欧美三级三区| 这个男人来自地球电影免费观看| 国产在线精品亚洲第一网站| 美女视频免费永久观看网站| 在线天堂中文资源库| 成人手机av| 中文字幕色久视频| 国产成+人综合+亚洲专区| 99国产极品粉嫩在线观看| 在线天堂中文资源库| 亚洲天堂av无毛| 91大片在线观看| 黄色成人免费大全| 国产免费av片在线观看野外av| 男女免费视频国产| 日本av手机在线免费观看| av不卡在线播放| 日本五十路高清| 美女高潮到喷水免费观看| 亚洲天堂av无毛| 人人妻人人添人人爽欧美一区卜| 少妇粗大呻吟视频| 91精品三级在线观看| 亚洲人成电影观看| 国产aⅴ精品一区二区三区波| av免费在线观看网站| 欧美在线黄色| 亚洲av欧美aⅴ国产| 久久久久视频综合| 免费日韩欧美在线观看| 国产主播在线观看一区二区| 日韩人妻精品一区2区三区| 婷婷丁香在线五月| 老司机在亚洲福利影院| 欧美激情 高清一区二区三区| 久久天躁狠狠躁夜夜2o2o| 欧美日韩亚洲国产一区二区在线观看 | 久久人妻熟女aⅴ| 久久久久久亚洲精品国产蜜桃av| 久久国产亚洲av麻豆专区| 国产成人欧美| 国产伦理片在线播放av一区| 亚洲精华国产精华精| av在线播放免费不卡| 少妇粗大呻吟视频| 人成视频在线观看免费观看| 亚洲av日韩精品久久久久久密| 国产xxxxx性猛交| 亚洲国产看品久久| 精品国产亚洲在线| 成年人午夜在线观看视频| 王馨瑶露胸无遮挡在线观看| 亚洲中文日韩欧美视频| 日日夜夜操网爽| 国产成人啪精品午夜网站| 亚洲精品国产一区二区精华液| 高清av免费在线| 午夜激情久久久久久久| 久久午夜亚洲精品久久| 欧美日韩成人在线一区二区| 超色免费av| 嫁个100分男人电影在线观看| 午夜福利在线观看吧| 亚洲人成伊人成综合网2020| 久久国产精品人妻蜜桃| 色婷婷久久久亚洲欧美| 国产野战对白在线观看| 五月开心婷婷网| 久久天躁狠狠躁夜夜2o2o| 日本vs欧美在线观看视频| 多毛熟女@视频| 一区二区三区国产精品乱码| 另类精品久久| 国产不卡av网站在线观看| 免费高清在线观看日韩| 99热国产这里只有精品6| 人妻久久中文字幕网| 真人做人爱边吃奶动态| 又紧又爽又黄一区二区| 免费少妇av软件| 精品视频人人做人人爽| 亚洲av国产av综合av卡| 国产精品久久电影中文字幕 | 一区二区三区激情视频| 国产老妇伦熟女老妇高清| 成人国产一区最新在线观看| 欧美日本中文国产一区发布| 欧美精品一区二区免费开放| 国产男女内射视频| 亚洲自偷自拍图片 自拍| 欧美亚洲日本最大视频资源| 成年动漫av网址| 亚洲性夜色夜夜综合| 美女主播在线视频| 久久ye,这里只有精品| 高清在线国产一区| 在线 av 中文字幕| 国产精品麻豆人妻色哟哟久久| 午夜福利欧美成人| 午夜福利在线免费观看网站| 国产精品 国内视频| 99久久精品国产亚洲精品| 国产在线精品亚洲第一网站| 国产精品亚洲一级av第二区| 一区福利在线观看| 国产av精品麻豆| 欧美日韩福利视频一区二区| 国产一区二区 视频在线| 国产精品一区二区精品视频观看| 日本欧美视频一区| 黑丝袜美女国产一区| 成人三级做爰电影| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产毛片av蜜桃av| 又大又爽又粗| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久成人av| 亚洲全国av大片| 久久中文看片网| 男女之事视频高清在线观看| 精品国产一区二区久久| 午夜日韩欧美国产| www.自偷自拍.com| 亚洲欧美日韩另类电影网站| 亚洲成人国产一区在线观看| 人成视频在线观看免费观看| 午夜成年电影在线免费观看| 脱女人内裤的视频| 欧美亚洲日本最大视频资源| 久久久国产欧美日韩av| 大香蕉久久网| 亚洲国产av新网站| 精品卡一卡二卡四卡免费| 妹子高潮喷水视频| 国产不卡一卡二| 老司机影院毛片| 欧美日韩福利视频一区二区| 精品国内亚洲2022精品成人 | 19禁男女啪啪无遮挡网站| 大型av网站在线播放| 免费不卡黄色视频| 老司机亚洲免费影院| 久久亚洲精品不卡| cao死你这个sao货| 久久久久精品国产欧美久久久| 九色亚洲精品在线播放| 一区二区日韩欧美中文字幕| 一级,二级,三级黄色视频| 国产精品亚洲av一区麻豆| 国精品久久久久久国模美| 亚洲精品成人av观看孕妇| 大香蕉久久网| 国产成人一区二区三区免费视频网站| 午夜日韩欧美国产| 老熟女久久久| 免费观看人在逋| 99久久国产精品久久久| 一级毛片精品| 精品亚洲乱码少妇综合久久| 狠狠婷婷综合久久久久久88av| 别揉我奶头~嗯~啊~动态视频| 国产成人一区二区三区免费视频网站| 99香蕉大伊视频| 国产淫语在线视频| 国产xxxxx性猛交| 欧美黄色片欧美黄色片| 深夜精品福利| 亚洲,欧美精品.| 欧美变态另类bdsm刘玥| 少妇被粗大的猛进出69影院| 美女高潮喷水抽搐中文字幕| 十八禁网站网址无遮挡| 国产伦人伦偷精品视频| 欧美日韩中文字幕国产精品一区二区三区 | 十分钟在线观看高清视频www| 美女扒开内裤让男人捅视频| 黄色怎么调成土黄色| 在线观看66精品国产| 久久精品熟女亚洲av麻豆精品| 肉色欧美久久久久久久蜜桃| 久久久精品国产亚洲av高清涩受| 免费一级毛片在线播放高清视频 | 久久久精品免费免费高清| 免费观看人在逋| 最新在线观看一区二区三区| 老司机深夜福利视频在线观看| 乱人伦中国视频| 露出奶头的视频| 国产精品秋霞免费鲁丝片| 久久影院123| 亚洲欧洲日产国产| 19禁男女啪啪无遮挡网站| 少妇被粗大的猛进出69影院| 女性生殖器流出的白浆| 无人区码免费观看不卡 | 涩涩av久久男人的天堂| 午夜福利免费观看在线| 一区二区av电影网| 免费一级毛片在线播放高清视频 | 午夜精品国产一区二区电影| 国产在线免费精品| 最新美女视频免费是黄的| 999久久久精品免费观看国产| 91麻豆精品激情在线观看国产 | 国产精品影院久久| 亚洲一码二码三码区别大吗| 久久久久视频综合| 老汉色∧v一级毛片| 捣出白浆h1v1| av福利片在线| 在线观看舔阴道视频| 亚洲美女黄片视频| 久久亚洲真实| 国产精品一区二区在线不卡| 在线十欧美十亚洲十日本专区| 国产精品二区激情视频| 人人妻人人澡人人看| 男男h啪啪无遮挡| 男人舔女人的私密视频| 日韩熟女老妇一区二区性免费视频| 99久久精品国产亚洲精品| 丰满少妇做爰视频| 巨乳人妻的诱惑在线观看| 午夜福利在线免费观看网站| 国产成人欧美在线观看 | 色尼玛亚洲综合影院| 啦啦啦在线免费观看视频4| 国产精品一区二区精品视频观看| 在线永久观看黄色视频| 麻豆av在线久日| 午夜日韩欧美国产| 久久精品熟女亚洲av麻豆精品| 亚洲av成人不卡在线观看播放网| 麻豆av在线久日| 男女高潮啪啪啪动态图| 水蜜桃什么品种好| a级片在线免费高清观看视频| 国产亚洲精品久久久久5区| 国产成人免费观看mmmm| 亚洲欧美精品综合一区二区三区| 国产欧美日韩一区二区三| a在线观看视频网站| av一本久久久久| 大片免费播放器 马上看| 欧美+亚洲+日韩+国产| 久久精品人人爽人人爽视色| 少妇的丰满在线观看| 777久久人妻少妇嫩草av网站| 欧美变态另类bdsm刘玥| 在线观看免费日韩欧美大片| 大片电影免费在线观看免费| videos熟女内射| 99国产精品一区二区蜜桃av | 成人国产一区最新在线观看| 国产不卡一卡二| 欧美在线一区亚洲| 黄色怎么调成土黄色| 蜜桃国产av成人99| 啪啪无遮挡十八禁网站| 最近最新中文字幕大全电影3 | 久久久久久久久久久久大奶| 久久99热这里只频精品6学生| 国产男靠女视频免费网站| 精品国产一区二区三区久久久樱花| 亚洲精品乱久久久久久| 国产精品99久久99久久久不卡| 欧美午夜高清在线| 亚洲欧美一区二区三区黑人| 精品乱码久久久久久99久播| avwww免费| 国产一区二区三区在线臀色熟女 | 精品第一国产精品| 国产日韩欧美视频二区| 亚洲av电影在线进入| 黄片大片在线免费观看| 多毛熟女@视频| 不卡av一区二区三区| 大香蕉久久网| 欧美另类亚洲清纯唯美| 亚洲成人免费电影在线观看| 老汉色av国产亚洲站长工具| 99久久精品国产亚洲精品| 热99国产精品久久久久久7| 不卡av一区二区三区| 搡老岳熟女国产| av超薄肉色丝袜交足视频| 最新美女视频免费是黄的| 嫩草影视91久久| 日本a在线网址| 9色porny在线观看| 国产野战对白在线观看| av网站免费在线观看视频| 大片电影免费在线观看免费| 黄色 视频免费看| 亚洲久久久国产精品| 国产97色在线日韩免费| 成年版毛片免费区| 十八禁网站网址无遮挡| 国产单亲对白刺激| 一级毛片女人18水好多| 日韩一区二区三区影片| 热99re8久久精品国产| 亚洲国产欧美在线一区| 国产在线精品亚洲第一网站| 美女国产高潮福利片在线看| 91麻豆精品激情在线观看国产 | 制服人妻中文乱码| 另类精品久久| 国产伦人伦偷精品视频| 国产精品一区二区在线不卡| 久久久久久久大尺度免费视频| 国产成人精品在线电影| www.熟女人妻精品国产| 成人国产一区最新在线观看| 中亚洲国语对白在线视频| 亚洲精品国产色婷婷电影| 亚洲国产欧美日韩在线播放| 制服诱惑二区| 老司机亚洲免费影院| 美女高潮喷水抽搐中文字幕| av又黄又爽大尺度在线免费看| 亚洲欧美日韩另类电影网站| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产精品一级二级三级| 成人av一区二区三区在线看| 香蕉国产在线看| 久久精品亚洲精品国产色婷小说| 亚洲,欧美精品.| 狠狠狠狠99中文字幕| 一级毛片电影观看| 欧美+亚洲+日韩+国产| 亚洲av欧美aⅴ国产| 两个人看的免费小视频| 亚洲美女黄片视频| 成人黄色视频免费在线看| 黄色毛片三级朝国网站| 老汉色av国产亚洲站长工具| 亚洲精品国产一区二区精华液| 亚洲精品中文字幕在线视频| 国产熟女午夜一区二区三区| 国产亚洲精品一区二区www | 老司机影院毛片| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 亚洲国产精品一区二区三区在线| 免费在线观看日本一区| 久久久水蜜桃国产精品网| 色综合欧美亚洲国产小说| 欧美日韩黄片免| 超色免费av| 亚洲专区字幕在线| 亚洲精品粉嫩美女一区| 亚洲精品av麻豆狂野| 黄色a级毛片大全视频| 蜜桃在线观看..| www.熟女人妻精品国产| 嫩草影视91久久| av网站免费在线观看视频| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 精品一区二区三区视频在线观看免费 | 久久亚洲精品不卡| 日本欧美视频一区| 99在线人妻在线中文字幕 | 久久影院123| 亚洲熟女毛片儿| 十八禁高潮呻吟视频| 亚洲国产成人一精品久久久| 午夜福利乱码中文字幕| 99riav亚洲国产免费| 国产精品欧美亚洲77777| 午夜激情av网站| 另类亚洲欧美激情| 老鸭窝网址在线观看| 国产男女内射视频| 欧美激情久久久久久爽电影 | 涩涩av久久男人的天堂| 成在线人永久免费视频| 99热网站在线观看| 久久精品亚洲精品国产色婷小说| 精品亚洲成a人片在线观看| 日韩欧美一区视频在线观看| 欧美老熟妇乱子伦牲交| 捣出白浆h1v1| 精品国产乱码久久久久久小说| 99国产精品99久久久久| 欧美+亚洲+日韩+国产| 免费高清在线观看日韩| 亚洲av国产av综合av卡| 婷婷丁香在线五月| 国产精品电影一区二区三区 | 国产亚洲精品久久久久5区| 菩萨蛮人人尽说江南好唐韦庄| 午夜免费鲁丝| 精品少妇一区二区三区视频日本电影| 天天影视国产精品| 91av网站免费观看| 狠狠狠狠99中文字幕| 国产亚洲欧美精品永久| av不卡在线播放| av线在线观看网站| 人人妻人人澡人人爽人人夜夜| 日韩一卡2卡3卡4卡2021年| 女人爽到高潮嗷嗷叫在线视频| 纯流量卡能插随身wifi吗| 91成人精品电影| 久久人妻熟女aⅴ| 亚洲av日韩在线播放| 欧美 亚洲 国产 日韩一| 亚洲成人免费电影在线观看| 另类亚洲欧美激情| 久久久久久久精品吃奶| 亚洲黑人精品在线| av又黄又爽大尺度在线免费看| 高清欧美精品videossex| 日韩免费av在线播放| 99re在线观看精品视频| 亚洲人成77777在线视频| 黄频高清免费视频| 亚洲免费av在线视频| 黄色视频在线播放观看不卡| 男女之事视频高清在线观看| 亚洲综合色网址| av不卡在线播放| 成人亚洲精品一区在线观看| 国产成人av激情在线播放| 久久精品国产亚洲av香蕉五月 | 国产亚洲一区二区精品| 日韩成人在线观看一区二区三区| 91麻豆精品激情在线观看国产 | 亚洲熟妇熟女久久| 午夜老司机福利片| 精品国产乱子伦一区二区三区| 热re99久久精品国产66热6| 日韩欧美三级三区| 日本黄色日本黄色录像| 久久九九热精品免费| 日日爽夜夜爽网站| 色在线成人网| 黑丝袜美女国产一区| 在线观看免费日韩欧美大片| 夫妻午夜视频| 欧美国产精品va在线观看不卡| 少妇粗大呻吟视频| 一本一本久久a久久精品综合妖精| 久久精品91无色码中文字幕| 亚洲国产中文字幕在线视频| 中文字幕av电影在线播放| 下体分泌物呈黄色| 999久久久国产精品视频| kizo精华| 老鸭窝网址在线观看| 亚洲精品在线观看二区| 欧美av亚洲av综合av国产av| 欧美在线一区亚洲| 亚洲av日韩在线播放| 法律面前人人平等表现在哪些方面| 国产日韩一区二区三区精品不卡| 亚洲午夜精品一区,二区,三区| 丁香欧美五月| 久久亚洲真实| 久久久久视频综合| 99精品在免费线老司机午夜| 12—13女人毛片做爰片一| 精品人妻在线不人妻| 少妇粗大呻吟视频| 19禁男女啪啪无遮挡网站| 久9热在线精品视频| 亚洲专区国产一区二区| 久久中文字幕一级| 日韩免费高清中文字幕av| 啦啦啦免费观看视频1| 国产淫语在线视频| 日韩一卡2卡3卡4卡2021年| 老司机靠b影院| 超色免费av| 在线观看66精品国产| 欧美精品一区二区大全| 精品久久久精品久久久| 国产午夜精品久久久久久| 黑人欧美特级aaaaaa片| 丝瓜视频免费看黄片| 肉色欧美久久久久久久蜜桃| 色综合欧美亚洲国产小说| 亚洲av片天天在线观看| 少妇精品久久久久久久| 日韩大码丰满熟妇| 色尼玛亚洲综合影院| 成人av一区二区三区在线看| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区在线不卡| 国产精品国产高清国产av | 久热这里只有精品99| 丝袜美腿诱惑在线| bbb黄色大片| 伊人久久大香线蕉亚洲五| 国产av又大| 亚洲精品在线美女| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩一区二区三区精品不卡| 他把我摸到了高潮在线观看 | 欧美性长视频在线观看| 欧美黑人精品巨大| 精品午夜福利视频在线观看一区 | 欧美国产精品va在线观看不卡| 久久狼人影院| 大码成人一级视频| 亚洲第一av免费看| 老司机午夜福利在线观看视频 | 免费观看人在逋| av电影中文网址| 一个人免费看片子| 国产成人欧美在线观看 | 亚洲av片天天在线观看| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| 一边摸一边抽搐一进一出视频| 欧美在线黄色| 女警被强在线播放| www.精华液| 精品福利观看| 国产亚洲精品久久久久5区| 免费日韩欧美在线观看| 亚洲黑人精品在线| 成年版毛片免费区| 国产黄色免费在线视频| 色综合婷婷激情| 国产精品av久久久久免费| 又黄又粗又硬又大视频| 黄频高清免费视频| 国产精品免费视频内射| 18禁美女被吸乳视频| 婷婷丁香在线五月| 人人妻人人添人人爽欧美一区卜| 国产精品99久久99久久久不卡| 色在线成人网| 亚洲免费av在线视频| 午夜两性在线视频| 国产精品98久久久久久宅男小说| 免费高清在线观看日韩| 中文字幕色久视频| 精品少妇一区二区三区视频日本电影| 999久久久精品免费观看国产| 欧美激情 高清一区二区三区| 一区二区av电影网| av国产精品久久久久影院| 国产精品免费视频内射| 蜜桃在线观看..| 亚洲av欧美aⅴ国产| 久久 成人 亚洲| 国产在线一区二区三区精| 肉色欧美久久久久久久蜜桃| 狂野欧美激情性xxxx| 亚洲专区国产一区二区| 男女之事视频高清在线观看| 婷婷丁香在线五月| 欧美性长视频在线观看| 伊人久久大香线蕉亚洲五| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 亚洲国产欧美在线一区| av线在线观看网站| 丝袜人妻中文字幕| 亚洲精品久久午夜乱码| 久久亚洲精品不卡| 成人国产av品久久久| 天天添夜夜摸| 精品一区二区三区av网在线观看 | 超碰97精品在线观看| 国产日韩欧美亚洲二区| 91国产中文字幕| 亚洲精品在线美女| 国产欧美日韩综合在线一区二区| 新久久久久国产一级毛片| 日韩欧美免费精品|