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

    影響阻尼牛頓法收斂性的兩個重要參數(shù)

    2012-07-05 14:30:10張建軍李春泉張烈輝
    關(guān)鍵詞:線性方程組病態(tài)初值

    張建軍,李春泉,張烈輝

    (1.西南石油大學理學院,四川成都 610500;2.西南石油大學石油工程學院,四川成都 610500)

    影響阻尼牛頓法收斂性的兩個重要參數(shù)

    張建軍1,李春泉1,張烈輝2

    (1.西南石油大學理學院,四川成都 610500;2.西南石油大學石油工程學院,四川成都 610500)

    對阻尼牛頓算法作了適當?shù)母倪M,證明了新算法的收斂性.基于新算法,運用計算機代數(shù)系統(tǒng)Matlab,研究了迭代次數(shù)k,參數(shù)對(μ,λ)與初值x0三者間的依賴關(guān)系,研究了病態(tài)問題在新算法下趨于穩(wěn)定的漸變(瞬變)過程.數(shù)值結(jié)果表明:(1)阻尼牛頓迭代中,參數(shù)對(μ,λ)與迭代次數(shù)k間存在特有的非線性關(guān)系;(2)適當?shù)膮?shù)對(μ,λ)與阻尼因子α的共同作用能夠在迭代中大幅度地降低病態(tài)問題的Jacobi陣的條件數(shù),使病態(tài)問題逐漸趨于穩(wěn)定,從而改變原問題的收斂性與收斂速度.

    阻尼牛頓法;雅可比矩陣的條件數(shù);病態(tài)問題;參數(shù)對(μ,λ)

    1 引言

    目前,非線性方程組求解的數(shù)值方法有Newton法,同倫型法,單純形法與胞腔排除法[13].

    牛頓迭代法是一種非常實用的計算方法,迭代公式如下:選擇參數(shù)α的原則是使之滿足:其中μ∈(0,1)為某定數(shù),‖·‖表示向量或矩陣的任何一種范數(shù).

    按(3)和(4)式進行迭代修正的方法稱為阻尼牛頓法,參數(shù)α稱為阻尼因子,當α=1時,就是牛頓法.牛頓法的顯著缺點是,它是局部收斂的.它要求初始近似x0適當接近非線性方程組(1)的解x*,這個條件難以界定,阻尼牛頓法是基于下降算法,并具有大范圍收斂性的一種迭代法.文獻[4]給出了求解非線性方程組(1)的阻尼牛頓算法,并給出了特殊情形下阻尼牛頓法的大范圍收斂性結(jié)果.文獻[4]中涉及到兩個重要參數(shù)μ,λ.但在文獻所給算法的5個步驟中沒有明確參數(shù)λ的作用;文獻中關(guān)于算法的收斂性定理也與參數(shù)λ沒有直接的關(guān)聯(lián).

    本文結(jié)構(gòu)安排如下:

    第二節(jié)對文獻[4]的算法作了適當?shù)母倪M,證明了新算法的收斂性.新的算法與相應的收斂性定理均明確了參數(shù)μ與λ各自的作用與相互關(guān)系.第三節(jié)研究了阻尼牛頓迭代中迭代次數(shù)k,參數(shù)對(μ,λ)與初值x0三者間的關(guān)系;第三節(jié)還研究了病態(tài)問題在適當?shù)膮?shù)對(μ,λ)和阻尼因子α的共同作用下趨于穩(wěn)定的漸變(瞬變)過程.

    2 主要結(jié)果

    2.1 改進的阻尼牛頓算法

    設μ,λ為固定的數(shù),且滿足

    2.2 收斂性定理

    下面證明,上述算法的第4步到第5步循環(huán)可在有限步內(nèi)完成.

    從而由(9)式推得(4)式對這個擴充的參數(shù)域成立.

    這個定理指出,若F′(x)在D中非奇異,則阻尼牛頓法在每個迭代中,α僅改變有限多次.

    3 數(shù)值試驗

    3.1 迭代次數(shù)k、參數(shù)對(μ,λ)與初值x0三者間的關(guān)系

    為了考察新算法中參數(shù)μ與λ對迭代過程的影響,選擇文獻[4-5]中的非線性方程組為例加以說明.

    例1考慮非線性方程組

    兩組近似解分別為(0.4597,-0.9038,-0.5494)與(0.4404,-1.0948,-0.5546),精度要求:ε=10-5. (一)固定初值x0=(5,5,5)T,令參數(shù)對(μ,λ)取10種不同的線性關(guān)系(滿足(5)式)

    對每個參數(shù)對(μ,λ)計算取得近似解所需要的迭代次數(shù)k,參見圖1.圖中10條直線(與10條曲線相對應)反映的是參數(shù)μ與λ的不同線性關(guān)系,而10條曲線反映的是相應的直線上各點所對應的迭代次數(shù).不同線性關(guān)系下的最小,最大迭代次數(shù)分別為7,10,11,12,13,15,17,19,24,33與282,286,294,495,533,564,591,1000,1000,1000.

    圖1 迭代次數(shù)k、參數(shù)對(μ,λ)與初值x0的關(guān)系(1)

    計算結(jié)果表明:①在不同線性關(guān)系中,l越小,取得近似解所需的迭代次數(shù)越小.②在不同線性關(guān)系中,參數(shù)對(μ,λ)與迭代次數(shù)k間存在相似的非線性關(guān)系,即隨著直線的下降,迭代次數(shù)先減后增.在直線的前端或中部位置,迭代次數(shù)達到最小.在直線的末端,隨著直線的下降,迭代次數(shù)劇烈增加.

    (二)固定參數(shù)μ與λ的線性關(guān)系,令初值x0變化

    在(14)式中取l=1/10.初值(x0,y0,z0)分別取為(5,5,5),(-5,-5,-5),(20,20,20),(-20, -20,-20).四種初值下的最小迭代次數(shù)分別為7,9,15與12(參見圖2).計算結(jié)果表明:①初值越接近非線性方程組的解x*,取得近似解所需的迭代次數(shù)越小.②當初值偏離近似解較大時,參數(shù)對(μ,λ)與迭代次數(shù)k的依賴關(guān)系可能會出現(xiàn)較大的波動,如初值為(20,20,20),(-20,-20,-20)時,出現(xiàn)迭代次數(shù)超過1000的情形.

    圖2 迭代次數(shù)k、參數(shù)對(μ,λ)與初值x0的關(guān)系(2)

    值得注意的是,即使是對應同一初值,在某些參數(shù)對(μ,λ)下,迭代過程收斂于(0.4404, -1.0948,-0.5546);而在另一些參數(shù)對(μ,λ)下,迭代過程收斂于(0.4597,-0.9038,-0.5 494).

    為了說明新算法相對于牛頓法的優(yōu)越性,應用牛頓法求解例1給出的非線性方程組.在初值(5,5,5),(-5,-5,-5),(-20,-20,-20)下牛頓法收斂,迭代次數(shù)分別為9,9,11.與阻尼牛頓法相比,牛頓法的收斂速度更快.原因是上述初值下,牛頓迭代中方程組并沒有出現(xiàn)奇異或數(shù)值穩(wěn)定性問題.在初值(20,20,20)處,由于牛頓迭代中第一次近似解的雅可比矩陣的條件數(shù)為無窮大,計算機停止計算,牛頓法失敗.此時,阻尼牛頓法顯示出了優(yōu)越性,如圖2所示,在100對(μ,λ)中,97對收斂且迭代次數(shù)均小于400,最小迭代次數(shù)為15.

    3.2 病態(tài)問題趨于穩(wěn)定的漸變(瞬變)過程分析

    由于函數(shù)F在近似解處的Jacobi陣的條件數(shù)的大小能夠反映方程組的數(shù)值穩(wěn)定性.因此,本文通過計算條件數(shù)跟蹤研究了新算法下病態(tài)問題趨于穩(wěn)定的漸變(瞬變)過程.

    例2繼續(xù)研究例1中的非線性方程組.如前所述,在牛頓迭代中,一旦初值遠離所給非線性方程組的解x*,牛頓法就可能失敗.表1給出了4組初值下牛頓迭代的計算結(jié)果(其中:ck=cond(F′(xk)),k=1,2,···表示第k次近似解的Jacobi陣的條件數(shù)).所列4組初值下牛頓法不收斂.原因之一是計算的穩(wěn)定性,第一次迭代后,它們的近似解的Jacobi陣的條件數(shù)都很大(或由于Jacobi陣的某些元素為無窮大而使計算機無法確定其條件數(shù)),導致了計算過程的不穩(wěn)定.因此,它們均屬于病態(tài)問題.另一個原因是計算的效率問題,根據(jù)文獻[6-8],當函數(shù)F在近似解處的Jacobi陣的條件數(shù)很大時,可能導致迭代中產(chǎn)生很小的步長,從而降低算法的效率.

    表1 牛頓迭代中的病態(tài)問題舉例

    新算法下,上述病態(tài)問題都能在15步之內(nèi)取得滿足誤差精度的近似解.跟蹤計算了收斂過程中條件數(shù)ck的變化狀況.計算結(jié)果見表2.通過對該表的觀察與對比,發(fā)現(xiàn)條件數(shù)的變化呈現(xiàn)出一定的規(guī)律性:①通過選擇適當?shù)膮?shù)對(μ,λ)及對阻尼因子α的調(diào)整,使得第一次近似解的Jacobi陣的條件數(shù)都得到大幅度地下降.②從第二步迭代開始,各次近似解的Jacobi陣的條件數(shù)以兩種方式下降,直至收斂.對某些病態(tài)問題(如P01等),迭代中Jacobi陣的條件數(shù)逐漸變小(單調(diào)下降方式或小幅波動下降方式)而使迭代過程收斂.而對另一些病態(tài)問題(如P02、P03、P04等),Jacobi陣的條件數(shù)先逐漸變大,達到最大值后,陡然變小而使迭代過程收斂.

    表2 阻尼牛頓迭代中條件數(shù)ck的變化狀況

    參考文獻

    [1]Xu Z B,Zhang J S,Wang W.A cell exclusion algorithm for determining all the solutions of a nonlinear system of equation[J].Appl.Math.and Computation,1996,80:181-203.

    [2]Zhang J S,You Z Y,Xu Z B.A cell method for finding all solutions of a nonlinear equations and optimization problems[J].Math.Numerica Sinica,1994,16:195-203.

    [3]孟彪龍,邢志棟.非線性方程組求全部解的胞腔排除法[J].純粹數(shù)學與應用數(shù)學,1999,15(2):72-76.

    [4]黃象鼎,曾鐘鋼,馬亞南.非線性數(shù)值分析的理論與方法[M].武漢:武漢大學出版社,2004.

    [5]施妙根,顧麗珍.科學和工程計算基礎(chǔ)[M].北京:清華大學出版社,1999.

    [6]Deuflhard P.A modified Newton method for the solution of ill-conditioned systems of nonlinear equations with applications to multiple shooting[J].Numer.Math.,1974,22:289-315.

    [7]曾金平,劉星果.自然水平函數(shù)在求解病態(tài)非線性方程組的阻尼牛頓法中的應用[J].湖南大學學報:自然科學版,2003,30(1):8-10.

    [8]周碩,郭麗杰,吳柏生.Jacobi迭代預處理的條件數(shù)與迭代次數(shù)的關(guān)系[J].東北電力學院學報,2003,23(6):57-60.

    The effect of two important parameter upon the convergence in damped Newton method

    Zhang Jianjun1,Li Chunquan1,Zhang Liehui2

    (1.College of Science,Southwest Petroleum University,Chengdu610500,China; 2.College of Petroleum Engineering,Southwest Petroleum University,Chengdu,610500China)

    In this paper,the damped Newton method is improved suitably and the convergence for new method is proved.Based on the new algorithm,a program is proposed and fulfilled by numerical and symbolic computations in Matlab,we study the relation among the iteration degrees k,parameters(μ,λ)and initial value x0.We also study the gradual(transient)process of the ill-conditioned systems nonlinear equations tending stable.Numerical results show that there is a special nonlinear relation between the parameters(μ,λ)and the iteration degrees k for damped Newton method and that the suitable parameters(μ,λ)and damping coefficient α can greatly decrease condition number of Jacobi matrix of ill-conditioned systems nonlinear equations.The ill-conditioned problems can gradually become stable and thus the convergence and the convergence speed of the ill-conditioned systems nonlinear equations can be changed.

    damped Newton method,condition number of the Jacobi matrix, ill-conditioned systems equations,parameters(μ,λ)

    O241.7

    A

    201109018(2012)04-0433-07

    2011-09-24.

    四川省教育廳2011年重點科研項目(10ZA073).

    張建軍(1958-),博士,副教授,研究方向:石油工程信息分析與處理.

    2010 MSC:65k05,90C30

    猜你喜歡
    線性方程組病態(tài)初值
    具非定常數(shù)初值的全變差方程解的漸近性
    求解非線性方程組的Newton迭代與Newton-Kazcmarz迭代的吸引域
    一種適用于平動點周期軌道初值計算的簡化路徑搜索修正法
    病態(tài)肥胖對門診全關(guān)節(jié)置換術(shù)一夜留院和早期并發(fā)癥的影響
    病態(tài)肥胖對門診關(guān)節(jié)置換術(shù)留夜觀察和早期并發(fā)癥的影響
    三維擬線性波方程的小初值光滑解
    君子之道:能移而相天——王夫之《莊子解》對“社會病態(tài)”的氣論診療
    哲學評論(2018年1期)2018-09-14 02:34:46
    線性方程組解的判別
    保護私有信息的一般線性方程組計算協(xié)議
    基于Matlab實現(xiàn)線性方程組的迭代解法
    国产高清视频在线观看网站| 97超级碰碰碰精品色视频在线观看| 国产视频内射| 在线观看舔阴道视频| 动漫黄色视频在线观看| 美女高潮喷水抽搐中文字幕| 91久久精品国产一区二区成人| 最近在线观看免费完整版| 如何舔出高潮| 国产精品一区二区性色av| 一级a爱片免费观看的视频| 露出奶头的视频| 在现免费观看毛片| 日本一本二区三区精品| 国产精品一区二区性色av| 黄片小视频在线播放| 色播亚洲综合网| 久久草成人影院| 91九色精品人成在线观看| 日韩精品青青久久久久久| 有码 亚洲区| 俺也久久电影网| 亚洲avbb在线观看| 免费av不卡在线播放| 国产一区二区在线观看日韩| 欧美成人a在线观看| 成熟少妇高潮喷水视频| 国产精品乱码一区二三区的特点| 亚洲欧美日韩无卡精品| 极品教师在线免费播放| 国产精品乱码一区二三区的特点| 悠悠久久av| 深夜精品福利| 男人和女人高潮做爰伦理| 天堂影院成人在线观看| 日本黄大片高清| 一进一出好大好爽视频| 久久九九热精品免费| 国产 一区 欧美 日韩| 欧美日韩综合久久久久久 | 国产野战对白在线观看| 乱人视频在线观看| 亚洲av不卡在线观看| 免费看光身美女| 深夜a级毛片| 免费观看人在逋| 久久久久性生活片| 成人高潮视频无遮挡免费网站| 神马国产精品三级电影在线观看| 琪琪午夜伦伦电影理论片6080| 久久午夜亚洲精品久久| 国产在线精品亚洲第一网站| 12—13女人毛片做爰片一| 好男人电影高清在线观看| 国产私拍福利视频在线观看| 亚洲七黄色美女视频| 欧美xxxx性猛交bbbb| 国产色婷婷99| 1000部很黄的大片| 90打野战视频偷拍视频| 欧美3d第一页| 国语自产精品视频在线第100页| 成人三级黄色视频| 两个人视频免费观看高清| 最新在线观看一区二区三区| 黄色一级大片看看| 久久精品夜夜夜夜夜久久蜜豆| 男女之事视频高清在线观看| 精品日产1卡2卡| 久久精品久久久久久噜噜老黄 | 最近在线观看免费完整版| 怎么达到女性高潮| 又粗又爽又猛毛片免费看| 国产亚洲欧美在线一区二区| 国产成人欧美在线观看| 99国产精品一区二区三区| 久久国产乱子免费精品| 中亚洲国语对白在线视频| 国产久久久一区二区三区| 免费av不卡在线播放| 宅男免费午夜| 老司机午夜十八禁免费视频| 欧美激情久久久久久爽电影| 国产精品,欧美在线| 亚洲欧美清纯卡通| 午夜视频国产福利| 成人三级黄色视频| 毛片一级片免费看久久久久 | 久久久精品大字幕| 婷婷精品国产亚洲av在线| 久久久精品欧美日韩精品| 麻豆av噜噜一区二区三区| 99国产精品一区二区三区| 午夜福利欧美成人| 久久久久九九精品影院| 人妻丰满熟妇av一区二区三区| 精品一区二区三区视频在线观看免费| 长腿黑丝高跟| 热99re8久久精品国产| 国产精品人妻久久久久久| 99国产极品粉嫩在线观看| 99国产综合亚洲精品| 国产单亲对白刺激| 在线观看免费视频日本深夜| 色综合站精品国产| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 亚洲精品日韩av片在线观看| 免费搜索国产男女视频| 怎么达到女性高潮| 亚洲av成人不卡在线观看播放网| or卡值多少钱| 久久精品国产99精品国产亚洲性色| 久久久久免费精品人妻一区二区| 婷婷色综合大香蕉| 久久九九热精品免费| 日日摸夜夜添夜夜添小说| 婷婷精品国产亚洲av在线| 51午夜福利影视在线观看| 免费黄网站久久成人精品 | 国产在线男女| av欧美777| 亚洲精品久久国产高清桃花| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月| 一区二区三区四区激情视频 | 精品乱码久久久久久99久播| 少妇人妻精品综合一区二区 | 国产欧美日韩精品亚洲av| 18禁黄网站禁片午夜丰满| 91狼人影院| 色噜噜av男人的天堂激情| 国模一区二区三区四区视频| 看免费av毛片| 精华霜和精华液先用哪个| 嫩草影院入口| 麻豆一二三区av精品| 午夜影院日韩av| 亚洲黑人精品在线| 亚洲av中文字字幕乱码综合| 九九久久精品国产亚洲av麻豆| 欧美成人一区二区免费高清观看| 美女被艹到高潮喷水动态| 日韩欧美精品v在线| 国产精品久久久久久人妻精品电影| 国产久久久一区二区三区| 国产一区二区三区视频了| or卡值多少钱| 九色国产91popny在线| 国产高清视频在线播放一区| 啪啪无遮挡十八禁网站| 国产色婷婷99| 国内揄拍国产精品人妻在线| 18禁黄网站禁片午夜丰满| 99国产精品一区二区三区| 香蕉av资源在线| 一级作爱视频免费观看| 亚洲专区中文字幕在线| 久久6这里有精品| 久久国产精品人妻蜜桃| 村上凉子中文字幕在线| 99riav亚洲国产免费| 中文字幕久久专区| 国内揄拍国产精品人妻在线| 精品无人区乱码1区二区| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| 国产精品久久电影中文字幕| 午夜福利在线在线| 嫁个100分男人电影在线观看| 啦啦啦韩国在线观看视频| 亚洲真实伦在线观看| 毛片一级片免费看久久久久 | 在线播放无遮挡| 国内揄拍国产精品人妻在线| 亚洲欧美激情综合另类| 亚洲精品在线美女| 亚洲国产精品成人综合色| 亚洲在线自拍视频| 婷婷色综合大香蕉| 在线免费观看不下载黄p国产 | 五月伊人婷婷丁香| 久久久色成人| 伦理电影大哥的女人| 国产成年人精品一区二区| 99热这里只有是精品在线观看 | 久久久久九九精品影院| 国产精品乱码一区二三区的特点| 国产精品一区二区性色av| a级毛片a级免费在线| 精品日产1卡2卡| 身体一侧抽搐| 九九热线精品视视频播放| 婷婷丁香在线五月| 大型黄色视频在线免费观看| 成人无遮挡网站| 欧美最黄视频在线播放免费| 欧美bdsm另类| 免费在线观看成人毛片| 亚洲专区国产一区二区| 一个人免费在线观看电影| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕一区二区三区有码在线看| 女生性感内裤真人,穿戴方法视频| 成年免费大片在线观看| 999久久久精品免费观看国产| 丰满乱子伦码专区| 九色国产91popny在线| 亚洲美女黄片视频| 成年人黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看| 精品欧美国产一区二区三| 欧美潮喷喷水| 久久久久久久久大av| 美女 人体艺术 gogo| 精品福利观看| 国产一区二区亚洲精品在线观看| 日韩欧美在线二视频| 精品人妻偷拍中文字幕| 亚洲成av人片免费观看| 变态另类成人亚洲欧美熟女| 男女之事视频高清在线观看| 久久久久久久久中文| 久久精品国产亚洲av香蕉五月| 麻豆成人午夜福利视频| 欧美区成人在线视频| 免费大片18禁| 国产激情偷乱视频一区二区| 久久中文看片网| 成年免费大片在线观看| 大型黄色视频在线免费观看| 久久久色成人| 亚洲人成网站在线播放欧美日韩| 日本免费a在线| 高潮久久久久久久久久久不卡| 午夜影院日韩av| 亚洲成av人片免费观看| 少妇熟女aⅴ在线视频| 变态另类丝袜制服| 国产在线精品亚洲第一网站| 黄色女人牲交| 精品99又大又爽又粗少妇毛片 | 日本精品一区二区三区蜜桃| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 少妇的逼好多水| 亚洲美女搞黄在线观看 | 国产精品av视频在线免费观看| 久久午夜亚洲精品久久| 久久国产精品影院| 色播亚洲综合网| 亚洲午夜理论影院| 有码 亚洲区| 久久久久精品国产欧美久久久| 国产精品影院久久| 亚洲成av人片免费观看| 99久久九九国产精品国产免费| 亚洲内射少妇av| 亚洲精品乱码久久久v下载方式| av在线蜜桃| 国产精品久久久久久人妻精品电影| 国产精品电影一区二区三区| 男女床上黄色一级片免费看| av黄色大香蕉| 97人妻精品一区二区三区麻豆| 精品欧美国产一区二区三| 精品福利观看| 久久6这里有精品| 9191精品国产免费久久| 丰满乱子伦码专区| 国产精品电影一区二区三区| 波多野结衣高清作品| 国产成人啪精品午夜网站| 日本与韩国留学比较| 精品一区二区三区视频在线| 日本黄色视频三级网站网址| 精品久久久久久久末码| 婷婷亚洲欧美| 国产高清有码在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 欧美潮喷喷水| 精品久久久久久久久av| 国产成人a区在线观看| avwww免费| 亚洲欧美激情综合另类| 国产高清视频在线播放一区| 久9热在线精品视频| 欧美最新免费一区二区三区 | 精品国内亚洲2022精品成人| 亚洲18禁久久av| 日本在线视频免费播放| 亚洲片人在线观看| 99热这里只有精品一区| eeuss影院久久| 97超级碰碰碰精品色视频在线观看| 欧美三级亚洲精品| 免费搜索国产男女视频| 精品福利观看| 啪啪无遮挡十八禁网站| 亚洲av五月六月丁香网| 丝袜美腿在线中文| 黄色一级大片看看| 亚洲国产精品成人综合色| 国产精品99久久久久久久久| 欧美日韩乱码在线| 搞女人的毛片| 国产伦人伦偷精品视频| 伊人久久精品亚洲午夜| 亚洲七黄色美女视频| 久久久久免费精品人妻一区二区| 国产精品日韩av在线免费观看| 直男gayav资源| 少妇丰满av| 午夜视频国产福利| 九色成人免费人妻av| 国产 一区 欧美 日韩| 久久精品国产亚洲av香蕉五月| 成年免费大片在线观看| 神马国产精品三级电影在线观看| 能在线免费观看的黄片| 国产乱人视频| 精品人妻熟女av久视频| 午夜福利高清视频| 国产精品美女特级片免费视频播放器| 亚洲电影在线观看av| 97热精品久久久久久| 在线观看舔阴道视频| 国产蜜桃级精品一区二区三区| 男女下面进入的视频免费午夜| 男人舔奶头视频| www.www免费av| 国产精品不卡视频一区二区 | 97热精品久久久久久| 欧美+日韩+精品| 长腿黑丝高跟| 国产精品永久免费网站| 大型黄色视频在线免费观看| 一个人观看的视频www高清免费观看| 黄色丝袜av网址大全| 国产高潮美女av| 免费人成视频x8x8入口观看| 免费av观看视频| 国产私拍福利视频在线观看| 国产精品久久久久久久久免 | 伊人久久精品亚洲午夜| 免费在线观看影片大全网站| 天堂√8在线中文| 91麻豆av在线| 好男人电影高清在线观看| 免费观看的影片在线观看| 免费看光身美女| 亚洲精品粉嫩美女一区| 色av中文字幕| 婷婷精品国产亚洲av在线| 99国产精品一区二区蜜桃av| 欧美国产日韩亚洲一区| 亚洲精品在线美女| 最近中文字幕高清免费大全6 | 亚洲黑人精品在线| 日韩大尺度精品在线看网址| 91av网一区二区| 毛片女人毛片| 99国产精品一区二区三区| 国产欧美日韩精品亚洲av| 欧美成人免费av一区二区三区| 内地一区二区视频在线| 国产精品1区2区在线观看.| 少妇的逼好多水| 特大巨黑吊av在线直播| 国产高清视频在线播放一区| 色精品久久人妻99蜜桃| 在线播放无遮挡| 一进一出好大好爽视频| 美女高潮的动态| 永久网站在线| 99国产精品一区二区蜜桃av| 久久人妻av系列| 久久精品久久久久久噜噜老黄 | 非洲黑人性xxxx精品又粗又长| 在线十欧美十亚洲十日本专区| 99久久精品热视频| 亚洲av不卡在线观看| 久久久久久九九精品二区国产| 中亚洲国语对白在线视频| 久久久久国产精品人妻aⅴ院| 国产精品爽爽va在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 精品国产亚洲在线| 在线观看美女被高潮喷水网站 | 国产极品精品免费视频能看的| 嫩草影院精品99| 看片在线看免费视频| 精品午夜福利在线看| 国产综合懂色| 久久这里只有精品中国| 亚洲欧美精品综合久久99| 亚洲成人精品中文字幕电影| 熟女人妻精品中文字幕| 欧美三级亚洲精品| 毛片一级片免费看久久久久 | netflix在线观看网站| 一二三四社区在线视频社区8| 美女高潮的动态| 国产高潮美女av| 18美女黄网站色大片免费观看| 91狼人影院| 深爱激情五月婷婷| 真人做人爱边吃奶动态| 国产私拍福利视频在线观看| 极品教师在线视频| 大型黄色视频在线免费观看| 青草久久国产| 婷婷精品国产亚洲av| 免费大片18禁| 99久久精品一区二区三区| 久久精品国产清高在天天线| 久久精品国产99精品国产亚洲性色| 男人舔奶头视频| 欧美成人一区二区免费高清观看| 国产探花极品一区二区| 少妇裸体淫交视频免费看高清| 国产精品野战在线观看| 99热只有精品国产| 婷婷丁香在线五月| 波多野结衣高清无吗| 欧美乱妇无乱码| 国产精品98久久久久久宅男小说| 我要看日韩黄色一级片| 9191精品国产免费久久| 俄罗斯特黄特色一大片| 欧美高清成人免费视频www| 国产一区二区三区在线臀色熟女| 国产探花在线观看一区二区| 亚洲av美国av| 日韩欧美在线二视频| 精品久久久久久久久亚洲 | 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 国产久久久一区二区三区| av欧美777| 久久伊人香网站| 色综合亚洲欧美另类图片| 亚洲第一电影网av| 日本 av在线| 欧美另类亚洲清纯唯美| 热99在线观看视频| 日本三级黄在线观看| 午夜福利欧美成人| 亚洲男人的天堂狠狠| 高清毛片免费观看视频网站| 啪啪无遮挡十八禁网站| 99久久精品一区二区三区| 色吧在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲不卡免费看| 亚洲第一区二区三区不卡| 久久精品综合一区二区三区| 久久久久久久亚洲中文字幕 | 中国美女看黄片| 成人精品一区二区免费| 有码 亚洲区| 免费在线观看日本一区| 老司机午夜十八禁免费视频| 亚洲成av人片在线播放无| 久久精品国产亚洲av香蕉五月| 禁无遮挡网站| 99久久九九国产精品国产免费| 亚洲,欧美精品.| 日韩欧美国产在线观看| 深夜a级毛片| 美女高潮的动态| 亚洲专区中文字幕在线| 在线观看午夜福利视频| 伊人久久精品亚洲午夜| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 日日摸夜夜添夜夜添小说| 免费av毛片视频| 嫁个100分男人电影在线观看| 久久这里只有精品中国| 国产日本99.免费观看| 91久久精品电影网| 午夜精品在线福利| 亚洲中文字幕日韩| 免费人成视频x8x8入口观看| 黄色日韩在线| 嫩草影院入口| 亚洲成人免费电影在线观看| 在线观看午夜福利视频| 国产视频一区二区在线看| 免费搜索国产男女视频| 99国产精品一区二区三区| 日韩欧美免费精品| 中文字幕人妻熟人妻熟丝袜美| 嫩草影视91久久| 亚洲欧美清纯卡通| 国语自产精品视频在线第100页| 国产精品女同一区二区软件 | 欧美日韩福利视频一区二区| 欧美激情久久久久久爽电影| 亚洲av熟女| 色5月婷婷丁香| 亚洲内射少妇av| 国产大屁股一区二区在线视频| 自拍偷自拍亚洲精品老妇| 很黄的视频免费| 精品一区二区三区视频在线观看免费| av视频在线观看入口| 久久精品国产亚洲av涩爱 | 黄片小视频在线播放| 日韩欧美在线二视频| 国产主播在线观看一区二区| 欧美中文日本在线观看视频| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 亚洲中文字幕一区二区三区有码在线看| 国产高清有码在线观看视频| 成人高潮视频无遮挡免费网站| 成人特级av手机在线观看| 看片在线看免费视频| 90打野战视频偷拍视频| 亚洲国产日韩欧美精品在线观看| 天堂√8在线中文| 看免费av毛片| 99热只有精品国产| 国产色爽女视频免费观看| 精品一区二区免费观看| 免费在线观看成人毛片| 欧美一区二区精品小视频在线| 日韩欧美精品免费久久 | or卡值多少钱| 在线播放无遮挡| 午夜福利免费观看在线| 小蜜桃在线观看免费完整版高清| 天堂动漫精品| 脱女人内裤的视频| 久久久久久国产a免费观看| 丁香六月欧美| 露出奶头的视频| 欧美+日韩+精品| 69av精品久久久久久| 婷婷精品国产亚洲av| 俺也久久电影网| av在线蜜桃| 欧美xxxx性猛交bbbb| 久久人人精品亚洲av| 国产精品久久久久久久久免 | 一个人免费在线观看电影| 在线十欧美十亚洲十日本专区| 日韩欧美免费精品| 99国产精品一区二区三区| 夜夜爽天天搞| 亚洲中文字幕日韩| 国产淫片久久久久久久久 | 午夜福利免费观看在线| 亚洲一区二区三区不卡视频| 天美传媒精品一区二区| 黄色日韩在线| 欧美成狂野欧美在线观看| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 又粗又爽又猛毛片免费看| 国产一区二区三区在线臀色熟女| 首页视频小说图片口味搜索| 村上凉子中文字幕在线| 亚洲一区高清亚洲精品| 美女免费视频网站| 婷婷丁香在线五月| 午夜a级毛片| 日本三级黄在线观看| 亚洲激情在线av| 亚洲av成人精品一区久久| 国产精品三级大全| 欧美日韩国产亚洲二区| 久久6这里有精品| 国产精品亚洲一级av第二区| 成年人黄色毛片网站| 国产精品影院久久| 国产69精品久久久久777片| 丁香欧美五月| 国产精品永久免费网站| 香蕉av资源在线| 国产在线精品亚洲第一网站| 有码 亚洲区| 韩国av一区二区三区四区| 欧美成人性av电影在线观看| 久久国产精品影院| 精品欧美国产一区二区三| 波多野结衣高清无吗| 如何舔出高潮| 国产精品,欧美在线| 赤兔流量卡办理| 一个人免费在线观看电影| 欧美最新免费一区二区三区 | 蜜桃亚洲精品一区二区三区| 他把我摸到了高潮在线观看| 国产色婷婷99| 黄色视频,在线免费观看| 国产精品一区二区三区四区免费观看 | 天天躁日日操中文字幕| 久久精品91蜜桃| 51国产日韩欧美| 国产不卡一卡二| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品久久视频播放| 亚洲成人精品中文字幕电影| av专区在线播放| 桃红色精品国产亚洲av| 少妇被粗大猛烈的视频| 国产精品野战在线观看| 日本一二三区视频观看| 亚洲片人在线观看|