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

    非零初始條件線性系統(tǒng)的Legendre多項式模型降階方法

    2016-12-05 00:50:57宋秋艷宋述剛
    長江大學學報(自科版) 2016年28期
    關(guān)鍵詞:降階狀態(tài)變量誤差

    宋秋艷,宋述剛

    (長江大學信息與數(shù)學學院,湖北 荊州 434023)

    ?

    非零初始條件線性系統(tǒng)的Legendre多項式模型降階方法

    宋秋艷,宋述剛

    (長江大學信息與數(shù)學學院,湖北 荊州 434023)

    討論了非零初始條件下線性系統(tǒng)的模型降階問題,并給出了一種基于Legendre正交多項式的時間域模型降階算法。該方法首先將系統(tǒng)的狀態(tài)變量在正交多項式空間中進行展開,然后由狀態(tài)方程得到展開系數(shù)的簡單遞推式,接著對其正交化,求得投影矩陣,通過正交投影變換得到降階系統(tǒng)。由該方法得到的降階系統(tǒng)可以匹配原始系統(tǒng)輸出變量一定數(shù)量的正交多項式展開系數(shù),從而保證了降階的精度。最后,通過2個數(shù)值算例驗證了該算法的有效性。

    模型降階;Legendre多項式;投影方法

    模型降階這一思想到現(xiàn)在已有40多年的歷史。簡而言之,模型降階就是在某種情況下將一個較大的復雜系統(tǒng)轉(zhuǎn)化為一個近似的較小系統(tǒng)的過程。模型降階是一種有效的降低動力系統(tǒng)復雜性的技術(shù)。該類技術(shù)能夠減少數(shù)據(jù)的存儲量和運算量,降低大型復雜系統(tǒng)的理論分析難度,加速系統(tǒng)的模擬計算,同時在一定的誤差范圍內(nèi)保持系統(tǒng)的某些重要屬性。模型降階方法已被成功地應用于許多工程應用領(lǐng)域和其他學科的分支中,如大規(guī)模集成電路模擬、自動化控制和機械工程等[1~3]。自模型降階方法被提出以來,已經(jīng)發(fā)展出多種方法,其中最主要的2類分別是Krylov子空間類方法和基于奇異值分解的平衡截斷模型降價方法。

    Krylov子空間方法是投影類模型降階方法,該類方法數(shù)值穩(wěn)定,算法實現(xiàn)簡單,并且計算量較小,受到很多科技工作者的青睞,但其沒有比較理想的誤差估計結(jié)果;平衡截斷模型降階方法可以直接得到降階模型的誤差估計結(jié)果,并能保持系統(tǒng)的穩(wěn)定性,但是該類方法在降階過程中需要求解大規(guī)模的Lyapunov方程,運算量較大,計算復雜度比較高,制約其應用的廣泛性。近些年來,由平衡截斷方法與投影類方法結(jié)合形成新型模型降階方法成為研究熱點之一。

    近年來,基于正交多項式(包括Chebyshev正交多項式、Legendre正交多項式和Laguerre正交多項式等)的模型降階方法也受到了廣泛的關(guān)注。這類算法的核心思想是首先將系統(tǒng)的狀態(tài)變量在以正交多項式為基底的空間中進行展開,然后由系統(tǒng)的狀態(tài)方程求得狀態(tài)變量的多項式展開系數(shù),最后通過該展開系數(shù)構(gòu)造標準列正交矩陣對原始系統(tǒng)進行降階。由此得到的降階系統(tǒng)一般都能夠匹配原始系統(tǒng)的輸出函數(shù)在正交多項式張成的空間中一定數(shù)量的展開系數(shù)。該類方法是一種時間域模型降階方法,已被成功地應用于線性系統(tǒng)、非線性系統(tǒng)以及一些特殊結(jié)構(gòu)系統(tǒng)的模型降階中[4~7]。

    傳統(tǒng)的模型降階方法,在降階過程中一般都只考慮系統(tǒng)的輸入輸出性態(tài),忽略初始條件的影響,或者為了簡便假設(shè)初始條件為零,這樣,使得原始系統(tǒng)的初始信息遭到破壞,給降階系統(tǒng)帶來不可預測的結(jié)果[8, 9]。因此,傳統(tǒng)的模型降階方法對非零初始條件的系統(tǒng)一般不太適用?;诖耍P者針對非零初始條件的線性系統(tǒng),提出一種基于Legendre正交多項式的時間域模型降階方法。

    1 Legendre正交多項式

    定義1 多項式:

    稱為Legendre多項式。

    性質(zhì)1 Legendre多項式在區(qū)間[-1,1]上滿足如下正交性:

    性質(zhì)2 對Legendre多項式,遞推公式(1)成立:

    (1)

    其中, P0(t)=1,P1(t)=t。

    Legendre多項式Pi(t)可以展開為如下的冪級數(shù):

    (2)

    式中, fij為冪級數(shù)tj的展開系數(shù)。

    將式(2)帶入式(1),可得:

    (3)

    比較式(3)兩邊關(guān)于t的各次冪的系數(shù),可得:

    任意一個可積函數(shù)x(t)均可以在Legendre正交多項式基底下近似展開:

    (4)

    由文獻[10] 可知,可測函數(shù)的正交多項式展開在Lebesgue意義下是一致收斂的,且正交多項式的近似展開在最小平方誤差意義下是最優(yōu)的,因此,相對低價的正交多項式近似可以達到較高的精度。

    可以將最小二乘法與Legendre正交多項式相結(jié)合來計算函數(shù)x(t)的冪級數(shù)展開系數(shù)xi。將式(2)代入到(4)中,比較兩邊t的各次冪的系數(shù),有:

    其中, αj的值可由Legendre多項式的正交性求得:

    在實際應用中,可以用文獻[11] 中復雜度為O(NlogN)的快速算法來計算函數(shù)x(t)的前N個Legendre正交多項式展開系數(shù)。

    2 基本算法

    下面,筆者將給出非零初始條件下線性系統(tǒng)的基于Legendre正交多項式的模型降階方法的具體過程。

    考慮如下單輸入單輸出線性系統(tǒng):

    (5)

    其中, A∈Rn×n;b,c∈Rn;x(t)∈Rn為系統(tǒng)的狀態(tài)變量; u(t),y(t)∈R分別是系統(tǒng)的輸入變量和輸出變量;n為系統(tǒng)的維數(shù)。

    為得到原始系統(tǒng)(5)的降階系統(tǒng),首先將系統(tǒng)的狀態(tài)變量x(t)和輸入變量u(t)近似展開為:

    (6)

    (7)

    其中,hi∈Rn,ui∈R分別為x(t)和u(t)的展開系數(shù)向量。

    將式(6)和(7)代入式(5),有:

    整理得:

    比較兩端常數(shù)項和ti(i=1,2,…,r-2)所對應的系數(shù),并忽略高階項tr-1,最終可得:

    (8)

    式中,h0為給定的初始條件x0。

    由遞推式(8)便可求得狀態(tài)變量x(t)的展開系數(shù)向量hi(i=1,2,…,r-1)。

    (9)

    3 數(shù)值算例

    下面,筆者通過2個數(shù)值算例來驗證上述模型降階方法的有效性。

    例1 考察一個實際的大氣風暴軌跡的地球大氣模型[3],由形如系統(tǒng)(5)的598階微分動力系統(tǒng)來描述,其初始條件為x0=[0,0,…,0,1]T。

    對該系統(tǒng)采用基于Legendre正交多項式的模型降階方法降至16階,原始系統(tǒng)與降階系統(tǒng)關(guān)于輸入函數(shù)u(t)=e-0.5tsin10t的瞬態(tài)響應及其相應的誤差如圖1、圖2所示。

    圖1 例1的瞬態(tài)響應

    圖2 例1中降階模型的相對誤差

    例2 考慮形如系統(tǒng)(5)的1006階微分動力系統(tǒng)[3],其中:

    A=diag{A1,A2,A3,A4}

    A3=diag{-1,-2,…,-1000}

    對該系統(tǒng)采用基于Legendre正交多項式的模型降階方法降至20階,原始系統(tǒng)與降階系統(tǒng)關(guān)于輸入函數(shù)u(t)=sint的瞬態(tài)響應及其相應的誤差如圖3、圖4所示。

    圖3 例2的瞬態(tài)響應

    圖4 例2中降階系統(tǒng)的絕對誤差

    由以上2個數(shù)值算例的模擬結(jié)果可以看出,基于Legendre正交多項式的模型降階方法得到的低階模型對原始模型有很好的近似效果。由于考慮了初始條件,該方法對于非零初始條件的線性系統(tǒng)是有效的。

    4 結(jié)語

    傳統(tǒng)的模型降階方法往往忽略初始條件,使得降階模型的精度無法保證。針對帶非零初始條件的線性系統(tǒng),提出了一種基于Legendre正交多項式的模型降階方法。該方法不僅考慮了初始條件,并且降階過程簡單高效,數(shù)值算例驗證了該方法的有效性。基于正交多項式的模型降階方法是一類時間域的模型降階方法,該類方法一般與系統(tǒng)的輸入函數(shù)有關(guān),如何削弱該類方法對輸入函數(shù)的依賴性,擴展其應用是值得進一步研究的問題。

    [1]蔣耀林. 模型降階方法[M]. 北京: 科學出版社, 2010.

    [2]Antoulas A C. Approximation of Large-Scale Dynamical Systems[M].Philadelphia:SIAM,2005.

    [3]Benner P, Mehrmann V, Sorensen D C. Dimension Reduction ofLarge-Scale Systems[M]. Springer-Verlag Berlin Heidelberg, 2005.

    [4]Wang J M, Chu C C, Yu Q J, et al. On projection-based algorithms formodel-order reduction of interconnects[J]. IEEE Trans Circuits Syst Regul Pap, 2002, 49(11): 1563~1585.

    [5]Jiang Y L, Chen H B. Time domain model order reduction of generalorthogonal polynomials for linear input-output systems[J]. IEEE Trans Autom Control, 2012, 57(2): 330~343.

    [6]Jiang YL, Chen H B. Application of general orthogonal polynomials tofast simulation of nonlinear descriptor systems throughpiecewise-linear approximation[J]. IEEE Trans Comput Aided Des Integr Circuits Syst, 2012, 31(5): 804~808.

    [7]Xiao Z H, Jiang Y L. Model order reduction of MIMO bilinear systems by multi-order Arnoldi method[J]. Syst Control Lett, 2016, 94: 1~10.

    [8] Heinkenschloss M, Reis T, Antoulas A C. Balanced truncation model reduction for systems with inhomogeneous initial conditions[J]. Automatica, 2011, 47(3): 559~564.

    [9] Shen J, Lam J. model reduction for discrete-time positive systems with inhomogeneous initial conditions[J]. Int J Robust Nonlin, 2015, 25(1): 88~102.

    [10] Szego G. Orthogonal polynomials[M]. New York city: American Mathematical Society, 1939.

    [11] Iserles A. A fast and simple algorithm for the computation ofLegendre coefficients[J]. Numer Math, 2011, 117(3): 529~553.

    [編輯] 洪云飛

    2016-06-26

    國家自然科學基金項目(11201039)。

    宋秋艷(1989-),女,碩士生,現(xiàn)主要從事應用數(shù)學方面的研究工作;通信作者:宋述剛,教授,2712281782@qq.com。

    O231

    A

    1673-1409(2016)28-0001-05

    [引著格式]宋秋艷,宋述剛.非零初始條件線性系統(tǒng)的Legendre多項式模型降階方法[J].長江大學學報(自科版),2016,13(28):1~5.

    猜你喜歡
    降階狀態(tài)變量誤差
    一階動態(tài)電路零狀態(tài)響應公式的通用拓展
    基于TwinCAT3控制系統(tǒng)的YB518型小盒透明紙包裝機運行速度的控制分析
    基于嵌套思路的飽和孔隙-裂隙介質(zhì)本構(gòu)理論
    單邊Lipschitz離散非線性系統(tǒng)的降階觀測器設(shè)計
    角接觸球軸承接觸角誤差控制
    哈爾濱軸承(2020年2期)2020-11-06 09:22:26
    Beidou, le système de navigation par satellite compatible et interopérable
    壓力容器制造誤差探究
    九十億分之一的“生死”誤差
    山東青年(2016年2期)2016-02-28 14:25:41
    降階原理在光伏NPC型逆變微網(wǎng)中的應用研究
    基于Krylov子空間法的柔性航天器降階研究
    √禁漫天堂资源中文www| 国产精品一区www在线观看| 纯流量卡能插随身wifi吗| 少妇被粗大的猛进出69影院 | 免费在线观看黄色视频的| 国国产精品蜜臀av免费| 久久久久网色| 国产日韩欧美在线精品| 国产免费一级a男人的天堂| 亚洲欧洲国产日韩| 王馨瑶露胸无遮挡在线观看| 男女午夜视频在线观看 | 欧美 亚洲 国产 日韩一| 狠狠精品人妻久久久久久综合| 色5月婷婷丁香| 亚洲精品视频女| 九色成人免费人妻av| 9191精品国产免费久久| av在线app专区| 久久久久久久久久久免费av| 日韩不卡一区二区三区视频在线| 看免费av毛片| 亚洲av综合色区一区| 国产精品不卡视频一区二区| 大香蕉久久成人网| 国产精品女同一区二区软件| 黄色 视频免费看| 久久久久久久久久成人| 哪个播放器可以免费观看大片| 男女国产视频网站| 在线观看www视频免费| 97超碰精品成人国产| 精品人妻偷拍中文字幕| 日韩av免费高清视频| 热99国产精品久久久久久7| 午夜福利视频在线观看免费| 满18在线观看网站| 夫妻午夜视频| 国产精品一二三区在线看| 另类亚洲欧美激情| freevideosex欧美| 免费大片18禁| 亚洲精品456在线播放app| 色5月婷婷丁香| 欧美日韩成人在线一区二区| 黄色怎么调成土黄色| 热re99久久精品国产66热6| 亚洲欧美成人精品一区二区| 哪个播放器可以免费观看大片| 色视频在线一区二区三区| 丰满饥渴人妻一区二区三| 久久亚洲国产成人精品v| 91午夜精品亚洲一区二区三区| av片东京热男人的天堂| 女性被躁到高潮视频| 午夜影院在线不卡| 亚洲国产精品一区二区三区在线| 一级毛片电影观看| 国产一级毛片在线| 90打野战视频偷拍视频| 精品久久久久久电影网| 日本欧美国产在线视频| 夜夜爽夜夜爽视频| 欧美xxⅹ黑人| 国产男人的电影天堂91| 五月天丁香电影| 精品一区二区三卡| 22中文网久久字幕| 人人澡人人妻人| 成人二区视频| 国产精品一区二区在线观看99| 亚洲欧洲国产日韩| 老司机亚洲免费影院| 日韩视频在线欧美| 久久99精品国语久久久| 少妇猛男粗大的猛烈进出视频| 午夜精品国产一区二区电影| 一级毛片黄色毛片免费观看视频| 国产精品欧美亚洲77777| 欧美xxxx性猛交bbbb| 精品一品国产午夜福利视频| 91国产中文字幕| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 满18在线观看网站| 9色porny在线观看| av国产久精品久网站免费入址| 男人舔女人的私密视频| 天堂中文最新版在线下载| 精品久久蜜臀av无| 九色亚洲精品在线播放| 捣出白浆h1v1| av免费在线看不卡| 在线免费观看不下载黄p国产| a级片在线免费高清观看视频| 亚洲精品,欧美精品| 日本爱情动作片www.在线观看| 欧美成人精品欧美一级黄| 黄色一级大片看看| 色婷婷av一区二区三区视频| 国产 精品1| 蜜桃在线观看..| 制服诱惑二区| 国产1区2区3区精品| 高清在线视频一区二区三区| 狠狠婷婷综合久久久久久88av| 大码成人一级视频| 99热这里只有是精品在线观看| 99re6热这里在线精品视频| 哪个播放器可以免费观看大片| 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 亚洲人成77777在线视频| 少妇人妻久久综合中文| 51国产日韩欧美| 欧美人与善性xxx| 欧美人与性动交α欧美精品济南到 | 中文欧美无线码| 亚洲欧美日韩卡通动漫| 丝袜喷水一区| 久久免费观看电影| 国产免费视频播放在线视频| 国产淫语在线视频| 美女主播在线视频| 美国免费a级毛片| 日韩视频在线欧美| 99香蕉大伊视频| videossex国产| 亚洲国产看品久久| 亚洲精品久久久久久婷婷小说| 黑人欧美特级aaaaaa片| 国产精品麻豆人妻色哟哟久久| 一级片免费观看大全| 免费大片18禁| 在线观看人妻少妇| 成人手机av| 亚洲色图综合在线观看| 一级爰片在线观看| 成人黄色视频免费在线看| 大片免费播放器 马上看| 国产又爽黄色视频| 精品国产乱码久久久久久小说| 午夜91福利影院| 亚洲国产看品久久| 国产极品粉嫩免费观看在线| 国产黄色视频一区二区在线观看| 免费av不卡在线播放| 18禁裸乳无遮挡动漫免费视频| 熟女人妻精品中文字幕| 精品酒店卫生间| 一本—道久久a久久精品蜜桃钙片| 人妻少妇偷人精品九色| 国产在线视频一区二区| 午夜激情久久久久久久| 狠狠精品人妻久久久久久综合| 精品熟女少妇av免费看| 天美传媒精品一区二区| 精品一品国产午夜福利视频| 亚洲国产av新网站| 国产精品国产三级国产专区5o| 免费大片黄手机在线观看| 成人国产av品久久久| 9色porny在线观看| 日本欧美视频一区| 成人国产麻豆网| 亚洲精品日本国产第一区| 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 丰满饥渴人妻一区二区三| 久久女婷五月综合色啪小说| 成人漫画全彩无遮挡| 高清毛片免费看| 国产成人精品无人区| 亚洲精品国产色婷婷电影| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 亚洲av电影在线观看一区二区三区| 中文字幕免费在线视频6| 久久99一区二区三区| 久久热在线av| 日韩视频在线欧美| 精品一区二区三区四区五区乱码 | 全区人妻精品视频| 日本av手机在线免费观看| 极品少妇高潮喷水抽搐| 亚洲精品久久久久久婷婷小说| 精品少妇内射三级| 中国美白少妇内射xxxbb| 视频区图区小说| 精品一区二区三区视频在线| 欧美激情极品国产一区二区三区 | 亚洲国产精品999| 一区二区av电影网| 激情视频va一区二区三区| 99精国产麻豆久久婷婷| 校园人妻丝袜中文字幕| av视频免费观看在线观看| 久久精品久久久久久噜噜老黄| 亚洲精品日本国产第一区| 捣出白浆h1v1| 欧美精品av麻豆av| 久久精品国产a三级三级三级| 人成视频在线观看免费观看| 9191精品国产免费久久| kizo精华| 九色成人免费人妻av| 三上悠亚av全集在线观看| 久久人妻熟女aⅴ| 视频中文字幕在线观看| 人妻 亚洲 视频| 26uuu在线亚洲综合色| 精品一区二区三区视频在线| 国产在线免费精品| 黑人欧美特级aaaaaa片| 秋霞伦理黄片| 国产精品人妻久久久影院| 亚洲成人一二三区av| 日本-黄色视频高清免费观看| 人人妻人人添人人爽欧美一区卜| 亚洲第一av免费看| 亚洲av成人精品一二三区| 曰老女人黄片| 日本免费在线观看一区| 国产综合精华液| 一级爰片在线观看| 中文字幕人妻熟女乱码| 亚洲国产色片| 精品熟女少妇av免费看| 好男人视频免费观看在线| 国产精品不卡视频一区二区| 最近最新中文字幕免费大全7| 青春草国产在线视频| 免费女性裸体啪啪无遮挡网站| av免费在线看不卡| 18禁动态无遮挡网站| 精品福利永久在线观看| 夜夜骑夜夜射夜夜干| 色5月婷婷丁香| 大码成人一级视频| 高清在线视频一区二区三区| 26uuu在线亚洲综合色| 侵犯人妻中文字幕一二三四区| 色网站视频免费| 9色porny在线观看| 亚洲人成网站在线观看播放| 美女国产视频在线观看| 成年av动漫网址| 亚洲av国产av综合av卡| 老司机亚洲免费影院| 亚洲中文av在线| 精品一区在线观看国产| 亚洲国产看品久久| 少妇的丰满在线观看| 国产精品蜜桃在线观看| av播播在线观看一区| 免费日韩欧美在线观看| 国产av精品麻豆| 免费大片18禁| 宅男免费午夜| 男女高潮啪啪啪动态图| 亚洲三级黄色毛片| 成人国产av品久久久| 午夜福利视频精品| 色5月婷婷丁香| 一边亲一边摸免费视频| 国产日韩欧美亚洲二区| 国产欧美日韩综合在线一区二区| 国产69精品久久久久777片| 最近最新中文字幕大全免费视频 | 欧美激情极品国产一区二区三区 | 丝瓜视频免费看黄片| 九草在线视频观看| 又大又黄又爽视频免费| 亚洲精品av麻豆狂野| 欧美成人精品欧美一级黄| 伦理电影大哥的女人| 天天躁夜夜躁狠狠久久av| 一级毛片电影观看| 亚洲精品第二区| 18+在线观看网站| 日本午夜av视频| 18禁在线无遮挡免费观看视频| 国产精品麻豆人妻色哟哟久久| 一级片免费观看大全| 国产综合精华液| 久久精品国产亚洲av天美| 国产成人精品久久久久久| xxx大片免费视频| 久久久久精品性色| 99re6热这里在线精品视频| 亚洲欧洲国产日韩| 夜夜骑夜夜射夜夜干| 亚洲精品久久成人aⅴ小说| 婷婷色麻豆天堂久久| 青青草视频在线视频观看| 国产一级毛片在线| 大香蕉久久网| 综合色丁香网| 欧美人与善性xxx| 男女国产视频网站| 国产av一区二区精品久久| 国产一区亚洲一区在线观看| 亚洲少妇的诱惑av| 少妇人妻久久综合中文| 午夜精品国产一区二区电影| 最黄视频免费看| 日韩一区二区三区影片| 日产精品乱码卡一卡2卡三| 欧美日韩av久久| 全区人妻精品视频| 中文字幕精品免费在线观看视频 | 曰老女人黄片| 亚洲一码二码三码区别大吗| 热99国产精品久久久久久7| 黑丝袜美女国产一区| 免费观看无遮挡的男女| 在线天堂中文资源库| 久久人人爽av亚洲精品天堂| 高清视频免费观看一区二区| 综合色丁香网| 国产一区二区三区综合在线观看 | 欧美最新免费一区二区三区| 男男h啪啪无遮挡| 国产在线视频一区二区| 久久青草综合色| 亚洲国产精品999| 一级毛片电影观看| 蜜桃在线观看..| 欧美另类一区| 97超碰精品成人国产| 欧美另类一区| 亚洲精品一二三| 免费av中文字幕在线| 亚洲第一av免费看| 高清av免费在线| 国产精品麻豆人妻色哟哟久久| 久久久久久伊人网av| 久久精品国产亚洲av天美| 99热国产这里只有精品6| 亚洲精品中文字幕在线视频| 成人黄色视频免费在线看| 激情视频va一区二区三区| 蜜臀久久99精品久久宅男| 51国产日韩欧美| 久久久久久久久久久久大奶| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 51国产日韩欧美| 国产精品 国内视频| 美女国产高潮福利片在线看| 成人毛片a级毛片在线播放| 老司机亚洲免费影院| 成年美女黄网站色视频大全免费| 国产亚洲精品久久久com| 九九在线视频观看精品| 高清欧美精品videossex| 成人亚洲欧美一区二区av| 久久午夜福利片| 成人亚洲精品一区在线观看| 9191精品国产免费久久| 日本免费在线观看一区| 国产片特级美女逼逼视频| 制服丝袜香蕉在线| 一本大道久久a久久精品| 国产一区二区三区综合在线观看 | av黄色大香蕉| av有码第一页| 欧美精品国产亚洲| 亚洲图色成人| 日韩一本色道免费dvd| 久久97久久精品| 中国美白少妇内射xxxbb| 男女边吃奶边做爰视频| 久久人人97超碰香蕉20202| 男女下面插进去视频免费观看 | 欧美激情极品国产一区二区三区 | 精品人妻在线不人妻| 日韩制服骚丝袜av| 成人二区视频| 91午夜精品亚洲一区二区三区| 我要看黄色一级片免费的| 国产精品.久久久| 免费av中文字幕在线| 老女人水多毛片| av免费观看日本| 中文乱码字字幕精品一区二区三区| 色94色欧美一区二区| 美女xxoo啪啪120秒动态图| 人妻 亚洲 视频| 我的女老师完整版在线观看| 亚洲av综合色区一区| 黄色一级大片看看| 久久人人爽人人爽人人片va| 亚洲五月色婷婷综合| 老司机影院毛片| 中文字幕最新亚洲高清| 少妇精品久久久久久久| 中国美白少妇内射xxxbb| 七月丁香在线播放| 热re99久久精品国产66热6| 欧美xxⅹ黑人| 美女福利国产在线| 午夜免费男女啪啪视频观看| 久久久久精品久久久久真实原创| 五月天丁香电影| 秋霞伦理黄片| 香蕉国产在线看| 最近手机中文字幕大全| 老女人水多毛片| 蜜臀久久99精品久久宅男| 我要看黄色一级片免费的| 国产成人精品福利久久| 亚洲精品国产色婷婷电影| 精品亚洲成国产av| 国产无遮挡羞羞视频在线观看| 成人亚洲欧美一区二区av| 高清av免费在线| 精品亚洲成国产av| 波多野结衣一区麻豆| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看 | 日本与韩国留学比较| 99久久精品国产国产毛片| 国产精品.久久久| 9色porny在线观看| 久久亚洲国产成人精品v| 国产精品一区二区在线不卡| 成年动漫av网址| 精品久久蜜臀av无| 国产又爽黄色视频| 亚洲av福利一区| 国产色爽女视频免费观看| 最近最新中文字幕免费大全7| 女的被弄到高潮叫床怎么办| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区国产| 男女免费视频国产| 中文乱码字字幕精品一区二区三区| 免费观看av网站的网址| 天天躁日日躁夜夜躁夜夜| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 久久狼人影院| 成人国语在线视频| 久久精品国产a三级三级三级| 久久这里只有精品19| 99国产极品粉嫩在线观看| bbb黄色大片| 美女视频免费永久观看网站| 超色免费av| 老司机在亚洲福利影院| 亚洲av第一区精品v没综合| 一本综合久久免费| 叶爱在线成人免费视频播放| 好看av亚洲va欧美ⅴa在| 美女高潮到喷水免费观看| 免费观看精品视频网站| 人妻一区二区av| 成人影院久久| 村上凉子中文字幕在线| 午夜老司机福利片| 亚洲五月色婷婷综合| 午夜亚洲福利在线播放| 人妻久久中文字幕网| 村上凉子中文字幕在线| 亚洲一区高清亚洲精品| 亚洲av成人不卡在线观看播放网| 老鸭窝网址在线观看| 精品久久久精品久久久| 亚洲av日韩精品久久久久久密| 高潮久久久久久久久久久不卡| 国产又爽黄色视频| 婷婷丁香在线五月| 久久久久久久国产电影| 在线观看免费日韩欧美大片| 欧美激情久久久久久爽电影 | 91精品三级在线观看| 精品久久蜜臀av无| 国产一区在线观看成人免费| 成人国语在线视频| 国产xxxxx性猛交| 一级片'在线观看视频| 亚洲欧美色中文字幕在线| 99香蕉大伊视频| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区免费开放| 超碰97精品在线观看| 国产av又大| 黄色 视频免费看| 国产成人av教育| 午夜免费成人在线视频| 国产单亲对白刺激| 天天影视国产精品| 色综合欧美亚洲国产小说| 精品一品国产午夜福利视频| 法律面前人人平等表现在哪些方面| 久久国产精品男人的天堂亚洲| 欧美一级毛片孕妇| 欧美av亚洲av综合av国产av| 在线观看免费午夜福利视频| 国产精品免费视频内射| 国产精品久久视频播放| 亚洲成人免费电影在线观看| 大香蕉久久成人网| 国产99白浆流出| 99香蕉大伊视频| 成人av一区二区三区在线看| 欧美日韩成人在线一区二区| 在线观看免费视频日本深夜| 日韩大码丰满熟妇| 变态另类成人亚洲欧美熟女 | 国产免费男女视频| 女人久久www免费人成看片| 亚洲在线自拍视频| 99国产精品免费福利视频| 国产aⅴ精品一区二区三区波| 男女床上黄色一级片免费看| 怎么达到女性高潮| 99在线人妻在线中文字幕 | 大香蕉久久网| 淫妇啪啪啪对白视频| 国产三级黄色录像| 王馨瑶露胸无遮挡在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 丰满饥渴人妻一区二区三| 青草久久国产| 婷婷成人精品国产| 99香蕉大伊视频| 少妇粗大呻吟视频| 热99re8久久精品国产| 99久久人妻综合| 极品人妻少妇av视频| 亚洲伊人色综图| 久久久久久亚洲精品国产蜜桃av| 欧美av亚洲av综合av国产av| 精品午夜福利视频在线观看一区| 日韩精品免费视频一区二区三区| 免费av中文字幕在线| 国产一卡二卡三卡精品| 国产成人av激情在线播放| 亚洲精品国产精品久久久不卡| 亚洲国产欧美网| 性色av乱码一区二区三区2| 国产精品1区2区在线观看. | 亚洲精品国产色婷婷电影| 精品国产美女av久久久久小说| 欧美精品av麻豆av| 久久天堂一区二区三区四区| av天堂久久9| 99热只有精品国产| 人人妻人人添人人爽欧美一区卜| 日韩熟女老妇一区二区性免费视频| 精品高清国产在线一区| 中文字幕人妻丝袜制服| 久久久久国产精品人妻aⅴ院 | 一二三四在线观看免费中文在| 大香蕉久久成人网| 午夜视频精品福利| 在线十欧美十亚洲十日本专区| 亚洲性夜色夜夜综合| 久久久久国产精品人妻aⅴ院 | 热99国产精品久久久久久7| 亚洲黑人精品在线| 女同久久另类99精品国产91| 久久久久久久久久久久大奶| 成年版毛片免费区| 久久性视频一级片| 亚洲色图av天堂| 人人妻人人澡人人爽人人夜夜| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩精品亚洲av| 自拍欧美九色日韩亚洲蝌蚪91| 男人舔女人的私密视频| 国产区一区二久久| 国产精品香港三级国产av潘金莲| √禁漫天堂资源中文www| 美女高潮喷水抽搐中文字幕| a级毛片在线看网站| 亚洲欧美日韩高清在线视频| 深夜精品福利| 国精品久久久久久国模美| 亚洲性夜色夜夜综合| 色老头精品视频在线观看| 电影成人av| 老司机在亚洲福利影院| av网站免费在线观看视频| 欧美日韩精品网址| 一边摸一边抽搐一进一小说 | 亚洲av日韩在线播放| 亚洲中文字幕日韩| 精品熟女少妇八av免费久了| 老熟女久久久| 国产单亲对白刺激| 亚洲视频免费观看视频| 在线国产一区二区在线| 亚洲av欧美aⅴ国产| 亚洲成人免费电影在线观看| 午夜精品在线福利| 一进一出抽搐gif免费好疼 | 欧美日韩成人在线一区二区| 欧美日韩黄片免| 老司机午夜十八禁免费视频| videos熟女内射| 久久久国产一区二区| 91麻豆精品激情在线观看国产 | 国产亚洲欧美精品永久| 黄色女人牲交| 欧美最黄视频在线播放免费 | 久久天堂一区二区三区四区| 我的亚洲天堂| 日韩视频一区二区在线观看| 女性被躁到高潮视频| 18禁裸乳无遮挡免费网站照片 | av线在线观看网站| 国产亚洲欧美98|