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

    船舶動(dòng)力定位反步逆最優(yōu)控制

    2017-02-09 05:23:37曌1
    關(guān)鍵詞:最優(yōu)控制無源坐標(biāo)系

    徐 海 祥, 瞿 洋, 余 文 曌1,

    ( 1.武漢理工大學(xué) 高性能船舶技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 湖北 武漢 430063;2.武漢理工大學(xué) 交通學(xué)院, 湖北 武漢 430063 )

    船舶動(dòng)力定位反步逆最優(yōu)控制

    徐 海 祥*1,2, 瞿 洋2, 余 文 曌1,2

    ( 1.武漢理工大學(xué) 高性能船舶技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 湖北 武漢 430063;2.武漢理工大學(xué) 交通學(xué)院, 湖北 武漢 430063 )

    針對(duì)動(dòng)力定位船舶定點(diǎn)定位過程中的高低頻運(yùn)動(dòng)特性,設(shè)計(jì)了一種基于無源濾波的反步逆最優(yōu)控制器.首先建立船舶的非線性方程,再利用線性化方程對(duì)應(yīng)的Riccati方程對(duì)船舶非線性方程進(jìn)行反步變換,最后基于Lyapunov函數(shù)設(shè)計(jì)了未建模擾動(dòng)環(huán)境力自適應(yīng)律和漸進(jìn)穩(wěn)定的控制律,同時(shí)滿足了局部性能指標(biāo)和全局性能指標(biāo).該方法從一個(gè)新的角度解決定點(diǎn)定位過程中船舶非線性系統(tǒng)最優(yōu)控制的問題,避免了傳統(tǒng)H∞魯棒控制需要求解Hamilton-Jacobi-Isaacs(HJI)方程的問題,保證了整個(gè)濾波和控制系統(tǒng)的穩(wěn)定性.同時(shí)還研究了誤差和控制權(quán)矩陣Q和R對(duì)系統(tǒng)性能的影響,確保了船舶低速定位過程中的能耗最優(yōu).動(dòng)力定位船舶定點(diǎn)定位的仿真結(jié)果驗(yàn)證了該方法的有效性.

    動(dòng)力定位;無源濾波;反步變換;最優(yōu)控制

    0 引 言

    動(dòng)力定位技術(shù)已廣泛用于深水鉆井、海底管線的檢測(cè)和調(diào)查、水下機(jī)器人作業(yè)、水下工程施工、油井增產(chǎn)措施和維修、平臺(tái)供應(yīng)、穿梭油船、浮式采油(有或無儲(chǔ)油設(shè)備)等作業(yè)中.對(duì)于海上長(zhǎng)時(shí)間定位的船舶,船舶的濾波和最優(yōu)控制問題對(duì)于節(jié)省燃料和減少污染具有重要意義.

    傳統(tǒng)動(dòng)力定位(DP)系統(tǒng)通常假定運(yùn)動(dòng)方程在36個(gè)艏向角每隔10°可線性化,并利用線性Kalman濾波濾除船舶的高頻運(yùn)動(dòng).對(duì)于每10°的線性化運(yùn)動(dòng)方程,最優(yōu)的Kalman濾波器相關(guān)參數(shù)和反饋控制器增益都需要隨著線性化過程不斷調(diào)整[1],并且相應(yīng)的噪聲協(xié)方差矩陣需要在線自適應(yīng)估計(jì)[2].相較于Kalman濾波,無源濾波參數(shù)易整定,并且實(shí)船試驗(yàn)也驗(yàn)證了其有效性[1].此后,F(xiàn)ossen論述了在低速和恒定速度時(shí),船舶系統(tǒng)具有良好的無源性,克服了無源理論在船舶控制系統(tǒng)中應(yīng)用的局限性[3].在無源理論的基礎(chǔ)上,卜德華等結(jié)合PID控制器驗(yàn)證了無源濾波器對(duì)濾除船舶位置和艏搖角高頻信息的有效性[4].針對(duì)船舶定點(diǎn)狀態(tài)環(huán)境最優(yōu)問題,F(xiàn)ossen等提出了環(huán)境最優(yōu)艏向控制,使船舶能根據(jù)外部緩變環(huán)境力自動(dòng)調(diào)節(jié)艏向,從而達(dá)到減小能耗的目的[5].對(duì)于固定艏向定位作業(yè),最優(yōu)控制問題通過設(shè)計(jì)H∞魯棒控制器來解決,但此方法需要求解復(fù)雜的Hamilton-Jacobi-Isaacs(HJI)方程[6-7].為了克服環(huán)境荷載對(duì)船舶系統(tǒng)的影響,Ngongi等利用H∞魯棒控制來抵抗外荷載,并用TS模糊模型來逼近DP系統(tǒng),最后通過求解LMI方程來保證該TS模糊控制器的穩(wěn)定性[8].相較于H∞魯棒控制,反步積分控制設(shè)計(jì)過程簡(jiǎn)單,被廣泛應(yīng)用于船舶動(dòng)力定位系統(tǒng)的設(shè)計(jì)當(dāng)中[9-10].為了滿足執(zhí)行機(jī)構(gòu)動(dòng)態(tài)特性的要求,Morishita等以及徐海祥等將執(zhí)行機(jī)構(gòu)的動(dòng)態(tài)特性考慮到反步積分控制器的設(shè)計(jì)過程中,使控制輸出更加平滑,減小了執(zhí)行機(jī)構(gòu)的磨損[11-12].由于反步積分控制器的設(shè)計(jì)基于船舶數(shù)學(xué)模型,模型的準(zhǔn)確性對(duì)控制效果有一定的影響.為了避免求解HJI方程,Ezal等針對(duì)單輸入單輸出系統(tǒng),首先提出了一種反步逆最優(yōu)控制器,將H∞最優(yōu)控制理論和反步積分控制理論相結(jié)合,同時(shí)滿足了局部最優(yōu)和全局最優(yōu)的性能指標(biāo)[13].在Ezal等的基礎(chǔ)上,Strand等將單輸入單輸出的反步逆最優(yōu)控制理論推廣到多輸入多輸出的錨泊輔助動(dòng)力定位中[14].隨后,Strand等又將反步逆最優(yōu)控制器和無源濾波相結(jié)合,進(jìn)一步推廣了反步逆最優(yōu)理論并保證了整個(gè)系統(tǒng)輸入狀態(tài)的穩(wěn)定(input-to-state stable)[15].由于多輸入多輸出的錨泊輔助動(dòng)力定位系統(tǒng)不具有一般性,Kim等在Ezal和Strand等的基礎(chǔ)上,推導(dǎo)了多輸入多輸出系統(tǒng)的一般理論形式[16].謝文博等在Strand等的基礎(chǔ)上,根據(jù)船舶動(dòng)力定位循跡控制中時(shí)變期望位置的特點(diǎn),推廣設(shè)計(jì)了時(shí)變反步逆最優(yōu)控制[17].

    以上文獻(xiàn)大都側(cè)重于濾波器或者控制器的設(shè)計(jì)穩(wěn)定,而較少考慮到濾波和控制器之間的相互影響,有時(shí)無法保證整個(gè)濾波和控制系統(tǒng)的穩(wěn)定性.在Strand等的研究基礎(chǔ)之上,本文結(jié)合動(dòng)力定位定點(diǎn)控制的具體要求,進(jìn)一步完善反步逆最優(yōu)控制理論在船舶動(dòng)力定位系統(tǒng)中的應(yīng)用,研究誤差和控制權(quán)矩陣Q和R對(duì)系統(tǒng)性能的影響,設(shè)計(jì)一種自適應(yīng)反步逆最優(yōu)控制器,以保證整個(gè)濾波和控制系統(tǒng)的穩(wěn)定性,降低船舶低速定位過程中的能耗.

    1 數(shù)學(xué)模型及其變換

    1.1 船舶動(dòng)力定位數(shù)學(xué)模型

    建立如圖1所示的3個(gè)坐標(biāo)系.在大地坐標(biāo)系XEYEZE下,船舶當(dāng)前低頻運(yùn)動(dòng)位置和船舶期望位置分別為η=(xyψ)T和ηd=(xdydψd)T.在船舶運(yùn)動(dòng)坐標(biāo)系XYZ下,船舶當(dāng)前低頻運(yùn)動(dòng)速度和期望速度分別為v=(uυr)T和vd=(udυdrd)T,其坐標(biāo)系原點(diǎn)位于船中縱剖面與重心G相距xg處.船舶固定坐標(biāo)系XDYDZD原點(diǎn)OD位于(xd,yd),相對(duì)于北東坐標(biāo)系旋轉(zhuǎn)了ψd.船舶固定坐標(biāo)系實(shí)際是為控制器的設(shè)計(jì)而引入的,可以認(rèn)為是一個(gè)“控制坐標(biāo)系”,其避免了船舶期望位置ηd對(duì)控制器的影響,使得船舶在期望位置ηd附近擁有良好的控制效果[14].

    大地坐標(biāo)系和船舶運(yùn)動(dòng)坐標(biāo)系下,船舶動(dòng)力定位數(shù)學(xué)模型為[18]

    (1)

    (2)

    (3)

    (4)

    y=η+Cωξ

    (5)

    式中:Aω為包含波浪主頻率ω0和相對(duì)阻尼ζ的定常矩陣;ωω和Eω分別為高斯白噪聲及其幅值矩陣;J(η)為坐標(biāo)系轉(zhuǎn)換矩陣;b∈R3×1,是作用于船體上的環(huán)境擾動(dòng)力,包括二階波浪慢漂力、風(fēng)和流的作用力;T為大時(shí)間常數(shù)矩陣;ωb為均值為零的高斯白噪聲矩陣;Eb為高斯白噪聲的幅值矩陣;M∈R3×3,為船舶慣性矩陣;D∈R3×3,為線性阻尼矩陣;τ∈R3×1,為作用于船體上的推力;y為船舶測(cè)量位置;Cωξ為船舶高頻運(yùn)動(dòng)位置,Cω=(0I).模型中相關(guān)矩陣的具體形式為

    M=m-Xu.000m-Yυ.mxg-Yr.0mxg-Nυ.Iz-Nr.?è?????÷÷÷,

    當(dāng)船舶前進(jìn)速度不為零時(shí),矩陣D一般是非線性的;但對(duì)于零速或者恒定速度,矩陣D可以假定為線性的[1].

    圖1 大地、船舶運(yùn)動(dòng)和船舶固定坐標(biāo)系

    1.2 基于數(shù)學(xué)模型的無源濾波器設(shè)計(jì)

    20世紀(jì)90年代,機(jī)器人等領(lǐng)域的無源控制已得到了較大的發(fā)展[19],而船舶無源控制直到90年代后期才得到應(yīng)用,限制無源理論在船舶中應(yīng)用的原因是船舶的附黏水質(zhì)量和阻尼項(xiàng)會(huì)隨著船舶在海浪中的遭遇頻率以及船速的變化而變化[3].在利用無源理論設(shè)計(jì)船舶的濾波器時(shí),需要考慮船舶是否具有無源性.Fossen已論證了當(dāng)船速較低時(shí)船舶具有良好的無源性[3].在一般工況下,環(huán)境擾動(dòng)力矩對(duì)船艏的擾動(dòng)角度在1°以內(nèi),而極限海況的擾動(dòng)角度在5°以內(nèi),因此可以做出如下的合理假設(shè):

    J(η)≈J(y)

    (6)

    根據(jù)文獻(xiàn)[1]和[18],利用式(6)并結(jié)合船舶動(dòng)力定位數(shù)學(xué)模型可設(shè)計(jì)如下船舶無源濾波器:

    ξ^.=Aωξ^+K1y

    (7)

    η^.=J(y)v^+K2y

    (8)

    b^.=-T-1b^+K3y

    (9)

    Mv^.+Dv^=τ+JT(y)b^+JT(y)K4y

    (10)

    (11)

    :y =y-y^

    式中;K1∈R6×3,K2,K3,K4∈R3×3,為濾波器增益矩陣.為使濾波器穩(wěn)定,結(jié)合Kalman-Yakubovich-Popov引理并利用頻域理論可以確定濾波器的相關(guān)增益矩陣[1].

    1.3 低頻運(yùn)動(dòng)數(shù)學(xué)模型的變換

    在大地坐標(biāo)系XEYEZE和船舶運(yùn)動(dòng)坐標(biāo)系XYZ下,無源濾波器設(shè)計(jì)用到的是低頻運(yùn)動(dòng)數(shù)學(xué)模型(8)和(10).為了避免船舶期望位置ηd對(duì)控制器設(shè)計(jì)的影響,控制器設(shè)計(jì)所用到的低頻運(yùn)動(dòng)模型為船舶固定坐標(biāo)系XDYDZD和船舶運(yùn)動(dòng)坐標(biāo)系XYZ下的變換模型.在XDYDZD下,偏差向量可以表示為[14]

    e=JT(ηd)(ηd-η)

    (12)

    記Jd=J(ηd),Je=JT(ηd)J(η),并考慮如下關(guān)系:

    (13)

    (14)

    對(duì)式(12)求導(dǎo)有

    (15)

    在定點(diǎn)定位和狀態(tài)保持控制中,期望位置偏差ed=(exdeydeψd)T和定位點(diǎn)處期望速度vd均為0,且濾波器的輸出是控制器的輸入.因此,式(15)可以簡(jiǎn)化為

    e^.=-Jev^-JTdK2y

    (16)

    同時(shí),式(10)可以寫為

    v^.=-M-1Dv^+M-1u+JTeJTdK4y

    (17)

    式中

    u=τ+JTeJTdb^

    (18)

    x.=f(x)+Bu+Hy

    (19)

    式中

    f(x)=-Jev^-M-1Dv^?è????÷÷,B=0B2?è???÷=0M-1?è???÷,

    x.=Ax+Bu+H0y

    (20)

    式中

    式(19)和(20)為船舶低頻運(yùn)動(dòng)方程(8)和(10)的變換形式,兩式均計(jì)及了船舶的期望位置ηd,避免了ηd對(duì)控制器設(shè)計(jì)的影響.

    2 控制器設(shè)計(jì)

    2.1 控制目標(biāo)

    (1)局部最優(yōu)目標(biāo)

    系統(tǒng)在所設(shè)定的期望位置xd=0附近時(shí),要實(shí)現(xiàn)的控制目標(biāo)是對(duì)于式(20)的線性系統(tǒng),設(shè)計(jì)一個(gè)穩(wěn)定的控制律,以滿足如下的性能指標(biāo):

    (21)

    式中:誤差權(quán)矩陣Q和控制權(quán)矩陣R均為對(duì)稱正定矩陣.對(duì)于最優(yōu)的衰減系數(shù)γ*,衰減系數(shù)γ應(yīng)滿足γ>γ*.

    (2)全局最優(yōu)目標(biāo)

    系統(tǒng)最終要實(shí)現(xiàn)的控制目標(biāo)是對(duì)于式(19)的非線性系統(tǒng),設(shè)計(jì)一個(gè)全局漸進(jìn)穩(wěn)定的控制律,以滿足如下的性能指標(biāo):

    (22)

    式中:q(x)和R*(x)為非線性系統(tǒng)中待確定的正定部分,其在反步逆最優(yōu)設(shè)計(jì)過程中滿足以下局部約束條件[15]:

    R*(0)=R

    (23)

    (24)

    因此,當(dāng)船舶靠近設(shè)定點(diǎn)xd=0時(shí),全局漸進(jìn)穩(wěn)定的控制律將退化為線性H∞(LQ)控制律,此時(shí)局部最優(yōu)目標(biāo)得到滿足.

    2.2 線性反步變換

    局部最優(yōu)性能指標(biāo)Jl所對(duì)應(yīng)的廣義Riccati代數(shù)方程為

    (25)

    正定矩陣P可以分解為

    P=LTΔL

    (26)

    式中:L為下三角矩陣,Δ為正定分塊對(duì)角矩陣,

    利用式(26),可將式(25)轉(zhuǎn)化為

    (27)

    式中

    設(shè)定新的變量z,通過反步變換,滿足全局目標(biāo)的非線性系統(tǒng)式(19)有如下映射關(guān)系[20]:

    (28)

    z.=Az+B0ul+H0y

    (29)

    考慮如下Lyapunov函數(shù):

    V=zTΔz

    (30)

    對(duì)式(30)求導(dǎo)并結(jié)合式(29)的線性系統(tǒng),局部最優(yōu)目標(biāo)所對(duì)應(yīng)的線性H∞控制律為

    (31)

    此時(shí),Lyapunov導(dǎo)函數(shù)滿足

    V.≤-zTQz-uTRu+γ2y Ty

    (32)

    在后文的控制器設(shè)計(jì)中,將會(huì)用到非線性反步逆最優(yōu)設(shè)計(jì)方法來構(gòu)建與式(29)相類似的非線性變換方程.在非線性反步逆最優(yōu)構(gòu)造過程中,都會(huì)確定每一步的虛擬控制量αi,最終存在全局最優(yōu)目標(biāo)所對(duì)應(yīng)的控制律:

    (33)

    滿足

    V.≤-q(z)-uT?R?(z)u?+γ2y Ty

    (34)

    其中q(z)和R*(z)的正定過程將會(huì)在后文中給出.

    2.3 非線性反步變換

    反步逆最優(yōu)控制器的設(shè)計(jì)首先通過構(gòu)建與式(29)相類似的非線性變換方程,再利用Lyapunov函數(shù)(30)推導(dǎo)反步逆最優(yōu)控制律.非線性變換方程的推導(dǎo)過程總共可分兩步進(jìn)行.

    步驟1 定義第一個(gè)偏差向量為

    (35)

    對(duì)式(35)求導(dǎo)得

    z.1=α1+(-Jev^-α1)+H1y

    (36)

    αh=-Π(z1)z1

    Ω(z1)=diag{3ρ1z1,3ρ2z2,3ρ3z3}

    式中:z1=(z1z2z3)T,ρ1、ρ2、ρ3為正常數(shù).

    定義第2個(gè)偏差向量為

    A12z2=-z2=-Jev^-α1

    (37)

    因此

    z.1=α1-z2+H1y

    (38)

    步驟2 對(duì)式(37)求導(dǎo)得

    (39)

    式中

    J.eJTe=(r-rd)S=(r^+δ-rd)S

    (40)

    (41)

    (42)

    式中

    G2(r^,z)=-JeM-1DJTe+r^S-A11-A22

    綜合式(38)和(42)可得與式(29)相類似的非線性變換方程:

    z.=Az+Bu+Gz+Hy +Σ

    (43)

    式中

    2.4 反步逆最優(yōu)控制器設(shè)計(jì)

    若用無源濾波器中的式(9)作為未建模環(huán)境力的估計(jì),則可選取和式(30)相同的Lyapunov函數(shù).雖然式(9)能夠?qū)崿F(xiàn)對(duì)未建模環(huán)境力的跟蹤,但是為了獲得較好的濾波效果,式(9)中的增益矩陣會(huì)受到限制,使得未建模環(huán)境力的跟蹤效率并不高.為了提高環(huán)境力的跟蹤效率,這里采用另外一種積分操作,可以取如下Lyapunov函數(shù):

    (44)

    (45)

    B2u+H2y +Σ2)+2zT1Δ1H1y -2b TΓ-1b^.

    (46)

    將式(27)代入上式,則

    (47)

    b^.=ΓJdJeBT2Δ2z2

    (48)

    (49)

    (50)

    式中

    反步逆最優(yōu)控制律可以取為

    (51)

    若令

    (52)

    則有

    (53)

    式中:q(z)應(yīng)該滿足正定條件,在假設(shè)船舶速度v和高頻運(yùn)動(dòng)艏向角速度δ有界的情況下,存在有界正實(shí)數(shù)β1和β2使得Υ滿足[15]

    (54)

    為了使整個(gè)系統(tǒng)穩(wěn)定,q(z)需滿足

    (55)

    結(jié)合式(54)和(55),存在0<θ<1,使式(53)滿足

    V.≤-(1-θ)q(z)-uT?R?u?+γ2y Ty

    (56)

    若q(z)滿足正定條件,則

    (57)

    方式一

    (58)

    方式二

    (59)

    式中χ為正定矩陣,可取為

    (60)

    方式三

    R-T?=u1(r^,z)I+u2(r^,z)R-1

    (61)

    (62)

    u2(r^,z)=1;λ1≥-λ2(1+λ1+λ2k)-1;k>0,λ1<-λ2{

    (63)

    (64)

    (65)

    式中:λmax(*)表示求取*的最大特征值.為了操作的簡(jiǎn)便,這里選用第3種選取方式.

    3 仿真結(jié)果與分析

    下面將通過計(jì)算機(jī)仿真來驗(yàn)證無源濾波和反步逆最優(yōu)算法的有效性.如圖2所示,本文仿真模型為縮尺比為1∶20的平臺(tái)供應(yīng)船模,船舶的慣性矩陣和阻尼矩陣通過CFD理論計(jì)算得到,船舶模型相關(guān)參數(shù)、環(huán)境力參數(shù)以及控制器等相關(guān)參數(shù)如表1所示,環(huán)境荷載的加載方式詳見文獻(xiàn)[18].

    為檢驗(yàn)未建模環(huán)境力自適應(yīng)律效果,使船舶狀態(tài)保持控制,初始時(shí)刻船舶位置為η0,期望位置也為η0,初始時(shí)刻船舶3個(gè)方向未建模環(huán)境力和力矩均設(shè)置為零.船舶在時(shí)刻t=200 s時(shí)受到表1中的環(huán)境力作用,在t=1 800 s環(huán)境力消失.船舶在東向、北向和艏向的位置如圖3所示.固定坐標(biāo)系下3個(gè)方向的自適應(yīng)未建模環(huán)境力如圖4所示.相應(yīng)的權(quán)矩陣選取為

    圖2 平臺(tái)供應(yīng)船模

    表1 相關(guān)參數(shù)設(shè)定值

    為檢驗(yàn)權(quán)矩陣Q和R對(duì)DP系統(tǒng)的能耗影響,使船舶作中近距離定點(diǎn)定位控制,初始時(shí)刻船舶位置為η0,期望位置為ηd=(5 m 5 m π/4)T,無環(huán)境力干擾.選取4組權(quán)矩陣.

    為了比較權(quán)矩陣對(duì)船舶能耗的影響,取如下總推力性能指標(biāo):

    (b) 船舶狀態(tài)保持北向位置

    (c) 船舶狀態(tài)保持艏向位置

    圖4 船舶狀態(tài)保持自適應(yīng)未建模環(huán)境力

    權(quán)矩陣為參數(shù)1、2、3和4時(shí),總推力性能指標(biāo)如圖5所示,船舶低頻東向位置、低頻北向位置和低頻艏向位置如圖6所示.

    針對(duì)本文中所涉及的無源濾波器,狀態(tài)保持仿真實(shí)驗(yàn)展示了相應(yīng)的濾波效果和未建模環(huán)境力

    圖5 總推力性能指標(biāo)

    (a) 船舶低頻東向位置

    (b) 船舶低頻北向位置

    (c) 船舶低頻艏向位置

    圖6 不同權(quán)矩陣下的控制效果

    Fig.6 The control results of different positive cost matrices

    自適應(yīng)律的效果.在突變未建模環(huán)境力的作用下,圖4展示了環(huán)境力自適應(yīng)律良好的跟蹤效果.圖3展示了在自適應(yīng)未建模環(huán)境力和實(shí)際環(huán)境力相近時(shí),船舶的位置和艏向?yàn)V波效果良好,在很大程度上濾除了高頻運(yùn)動(dòng).在t=200 s至t=1 800 s,自適應(yīng)未建模環(huán)境力和實(shí)際環(huán)境力存在較大偏差,導(dǎo)致船舶位置估計(jì)尤其是艏向角估計(jì)和實(shí)際低頻位置存在一定的偏差.造成未建模環(huán)境力和實(shí)際環(huán)境力有偏差的原因往往是突風(fēng)、大浪或者激流等.由于控制是響應(yīng)緩變環(huán)境力,變化較大且頻繁的外部環(huán)境力對(duì)船舶控制將會(huì)是一個(gè)巨大的挑戰(zhàn).

    圖5展示了不同權(quán)矩陣所對(duì)應(yīng)的總推力性能指標(biāo),船舶的低速運(yùn)動(dòng)意味著船舶較長(zhǎng)的響應(yīng)時(shí)間,與此同時(shí)相應(yīng)的能耗也較小.圖6給出了不同權(quán)矩陣對(duì)定位效果的影響.根據(jù)權(quán)矩陣的定義可知,參數(shù)Q的前3個(gè)元素對(duì)應(yīng)船舶的位置權(quán)重,而后3個(gè)元素對(duì)應(yīng)船舶的速度權(quán)重;參數(shù)R對(duì)應(yīng)控制器的輸出推力權(quán)重.參數(shù)2相較于參數(shù)1增大了位置權(quán)重,船舶的速度有所增大,相應(yīng)的超調(diào)量和響應(yīng)時(shí)間等系統(tǒng)響應(yīng)特性都得到了改善.參數(shù)3相較于參數(shù)1增大了推力權(quán)重,船舶位置的響應(yīng)時(shí)間變長(zhǎng)和超調(diào)量有所增加,這是由于增大推力權(quán)重會(huì)在一定程度上限制控制器的推力輸出.參數(shù)4相較于參數(shù)1增大了船舶速度權(quán)重,船舶位置的響應(yīng)時(shí)間明顯變長(zhǎng),且基本無超調(diào),這是因?yàn)樵龃笏俣葯?quán)重限制了船舶速度.值得注意的是,過大的船舶速度會(huì)破壞船舶系統(tǒng)的無源性,即船舶的附加質(zhì)量和阻尼明顯的非線性,將會(huì)對(duì)船舶系統(tǒng)的穩(wěn)定性有著明顯的影響.因此,為了使整個(gè)動(dòng)力定位控制系統(tǒng)穩(wěn)定,控制船舶的速度不宜過高是十分必要的.

    4 結(jié) 語

    本文針對(duì)動(dòng)力定位船舶定點(diǎn)定位問題設(shè)計(jì)了一種基于無源濾波的自適應(yīng)反步逆最優(yōu)控制器,避免了H∞魯棒控制器設(shè)計(jì)過程中需要求解HJI方程的問題,滿足了局部最優(yōu)和全局最優(yōu)的性能指標(biāo),對(duì)外界緩變未建模環(huán)境力進(jìn)行補(bǔ)償?shù)耐瑫r(shí),也保證了無源濾波和控制整個(gè)系統(tǒng)的穩(wěn)定性.不同權(quán)矩陣Q和R對(duì)動(dòng)力定位控制系統(tǒng)的響應(yīng)時(shí)間、超調(diào)量等動(dòng)態(tài)特性和控制推力有著直接的影響,在滿足系統(tǒng)穩(wěn)定性的條件下,可通過優(yōu)化性能權(quán)函數(shù)Q和R以達(dá)到提高控制精度與降低能耗的目的.

    [1] FOSSEN T I, STRAND J P. Passive nonlinear observer design for ships using Lyapunov methods:full-scale experiments with a supply vessel [J]. Automatica, 1999, 35(1):3-16.

    [2] FOSSEN T I, PEREZ T. Kalman filtering for positioning and heading control of ships and offshore rigs [J]. IEEE Control Systems, 2010, CST-29(6):32-46.

    [3] FOSSEN T I. Nonlinear passive control and observer design for ships [J]. Modeling, Identification and Control, 2000, 21(3):129-184.

    [4] 卜德華,徐海祥,李文娟,等. 基于非線性無源濾波器的船舶動(dòng)力定位仿真[J]. 武漢理工大學(xué)學(xué)報(bào), 2013, 35(10):69-73.

    BU Dehua, XU Haixiang, LI Wenjuan,etal. Simulation of ship dynamic positioning based on nonlinear passive filter [J]. Journal of Wuhan University of Technology, 2013, 35(10):69-73. (in Chinese)

    [5] FOSSEN T I, STRAND J P. Nonlinear passive weather optimal positioning control (WOPC) system for ships and rigs:experimental results [J]. Automatica, 2001, 37(5):701-715.

    [6] WU Huaining, LUO Biao. Neural network based online simultaneous policy update algorithm for solving the HJI equation in nonlinear control [J]. IEEE Transactions on Neural Networks and Learning Systems, 2012, 23(12):1884-1895.

    [7] CHEN W, ANDERSON B D O. On the convergence rate of the Leake-Liu algorithm for solving Hamilton-Jacobi-Bellman equation [J]. IFAC Proceedings Volumes, 2011, 44(1):8064-8069.

    [8] NGONGI W E, DU J, WANG R. Robust fuzzy controller design for dynamic positioning system of ships [J]. International Journal of Control Automation & Systems, 2015, 13(5):1294-1305.

    [9] Jr. ZAKARTCHOUK A, MORISHITA H M. A backstepping controller for dynamic positioning of ships:numerical and experimental results for a shuttle tanker model [J]. IFAC Proceedings Volumes, 2009, 42(18):394-399.

    [10] WITKOWSKA A. Control system design for dynamic positioning using vectorial backstepping [J]. Scientific Journals of the Maritime University of Szczecin, 2013, 108(36):182-187.

    [11] MORISHITA H M, SOUZA C E S. Modified observer backstepping controller for a dynamic positioning system [J]. Control Engineering Practice, 2014, 33:105-114.

    [12] 徐海祥,瞿 洋,余文曌,等. 基于動(dòng)態(tài)執(zhí)行機(jī)構(gòu)的欠驅(qū)船舶循跡控制[J]. 武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版), 2016, 40(1):6-10.

    XU Haixiang, QU Yang, YU Wenzhao,etal. Path following for underactuated marine craft with actuator dynamics [J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2016, 40(1):6-10. (in Chinese)

    [14] STRAND J P, EZAL K, FOSSEN T I,etal. Nonlinear control of ships:a locally optimal design [M] // Preprints of the IFAC NOLCOS 1998. Enschede:The Netherlands, 1998:732-738.

    [15] STRAND J P, FOSSEN T I. Nonlinear output feedback and locally optimal control of dynamically positioned ships:Experimental results [C] // Proceedings of IFAC Conference on Control Application in Marine Systems. Fukuoka:IFAC, 1998:89-95.

    [16] KIM H, BACK J, SHIM H,etal. Locally optimal and globally inverse optimal controller for multi-input nonlinear systems [C] // Proceedings of 2008 American Control Conference. Seattle: AACC, 2008:4486-4491.

    [17] 謝文博,付明玉,張 健,等. 動(dòng)力定位船舶自適應(yīng)反步逆最優(yōu)循跡控制[J]. 中國(guó)造船, 2013, 54(3):58-69.

    XIE Wenbo, FU Mingyu, ZHANG Jian,etal. Adaptive backstepping inverse optimal tracking control of dynamic positioning vessel [J]. Shipbuilding of China, 2013, 54(3):58-69. (in Chinese)

    [18] FOSSEN T I. Handbook of Marine Craft Hydrodynamics and Motion Control [M]. New York:John Wiley & Sons, 2011.

    [19] BERGHUIS H, NIJMEIJER H. Passivity approach to controller-observer design for robots [J]. IEEE Transactions on Robotics & Automation, 1993, 9(6):740-754.

    [20] FOSSEN T I, STRAND J P. Tutorial on nonlinear backstepping:Applications to ship control [J]. Modeling, Identification and Control, 1999, 20(2):83-135.

    Inverse optimal backstepping control of dynamic positioning ships

    XU Haixiang*1,2, QU Yang2, YU Wenzhao1,2

    ( 1.Key Laboratory of High Performance Ship Technology of Ministry of Education, Wuhan University of Technology, Wuhan 430063, China; 2.School of Transportation, Wuhan University of Technology, Wuhan 430063, China )

    In view of the characteristics of high and low frequency motion in ship dynamic locating and positioning, an inverse optimal backstepping controller is designed cascaded with a passive filter. Firstly, the nonlinear equation of the ship is established. Then, the Riccati equation corresponding to the linear equation is used to carry out the backstepping transformation of the nonlinear equation. Finally, the adaptive law and the asymptotic stable control law for disturbed environmental force without modeling are designed based on the Lyapunov function, which meets both the local cost function and global cost function. This new approach obtains a solution of the optimal control of nonlinear system in a new direction, avoids solving the Hamilton-Jacobi-Isaacs(HJI) equation in traditionalH∞controller design and simultaneously guarantees the stability of filter and controller. The influences of error and control weight matrixQandRon the system performance are discussed, and the optimal energy consumption during dynamic positioning of ship at slow speed is achieved. The simulation results of a dynamic locating and positioning ship show the effectiveness of this approach.

    dynamic positioning; passive filter; backstepping transformation; optimal control

    1000-8608(2017)01-0046-09

    2016-01-07;

    2016-09-20.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(61301279,51479158);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(163102006).

    徐海祥*(1975-),男,教授,博士生導(dǎo)師,E-mail:qukaiyang@163.com;瞿 洋(1991-),男,碩士生,E-mail:yangqu91@163.com;余文曌(1989-),男,講師,E-mail:yuwenzhao1989@gmail.com.

    U674.38

    A

    10.7511/dllgxb201701007

    猜你喜歡
    最優(yōu)控制無源坐標(biāo)系
    條件平均場(chǎng)隨機(jī)微分方程的最優(yōu)控制問題
    帶跳躍平均場(chǎng)倒向隨機(jī)微分方程的線性二次最優(yōu)控制
    一種三相無源逆變電源供電方案設(shè)計(jì)
    電子制作(2019年12期)2019-07-16 08:45:14
    Timoshenko梁的邊界最優(yōu)控制
    解密坐標(biāo)系中的平移變換
    坐標(biāo)系背后的故事
    基于PCH模型的航天器姿態(tài)無源控制
    基于重心坐標(biāo)系的平面幾何證明的探討
    無源互調(diào)干擾對(duì)TD-LTE系統(tǒng)的影響研究
    采用最優(yōu)控制無功STATCOM 功率流的解決方案
    久久青草综合色| av国产精品久久久久影院| av在线播放精品| 久久人人爽人人片av| 新久久久久国产一级毛片| 免费看不卡的av| 啦啦啦 在线观看视频| www.自偷自拍.com| 男人舔女人的私密视频| 国产真人三级小视频在线观看| 国产欧美日韩一区二区三 | 一本—道久久a久久精品蜜桃钙片| 免费在线观看黄色视频的| 又大又黄又爽视频免费| 亚洲国产av新网站| 99精品久久久久人妻精品| 波多野结衣一区麻豆| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲天堂av无毛| av福利片在线| 别揉我奶头~嗯~啊~动态视频 | 日韩av免费高清视频| 久久国产亚洲av麻豆专区| 国产精品免费视频内射| 人妻 亚洲 视频| 免费观看人在逋| 日本午夜av视频| 国产一区二区 视频在线| 欧美日韩av久久| 精品亚洲成国产av| 欧美激情高清一区二区三区| 七月丁香在线播放| 我要看黄色一级片免费的| 悠悠久久av| 亚洲视频免费观看视频| 国产精品偷伦视频观看了| 国产激情久久老熟女| 一本久久精品| 无遮挡黄片免费观看| 90打野战视频偷拍视频| 18禁黄网站禁片午夜丰满| 91老司机精品| 美女视频免费永久观看网站| www.av在线官网国产| 久久九九热精品免费| 男女国产视频网站| 水蜜桃什么品种好| 青春草视频在线免费观看| 欧美激情高清一区二区三区| 一本—道久久a久久精品蜜桃钙片| 国产成人免费无遮挡视频| 一区二区三区激情视频| 久久天堂一区二区三区四区| 免费高清在线观看视频在线观看| 制服人妻中文乱码| 亚洲精品成人av观看孕妇| avwww免费| 免费观看av网站的网址| xxx大片免费视频| 精品久久久精品久久久| 久久国产精品男人的天堂亚洲| 欧美激情高清一区二区三区| av网站在线播放免费| 色婷婷久久久亚洲欧美| 欧美老熟妇乱子伦牲交| 国产精品一国产av| 91精品伊人久久大香线蕉| 亚洲 国产 在线| 五月开心婷婷网| 别揉我奶头~嗯~啊~动态视频 | 午夜精品国产一区二区电影| 看免费成人av毛片| 中文字幕人妻熟女乱码| 一区二区三区乱码不卡18| 男女午夜视频在线观看| 欧美在线一区亚洲| 欧美中文综合在线视频| 麻豆乱淫一区二区| 日韩,欧美,国产一区二区三区| 日本午夜av视频| 亚洲精品av麻豆狂野| 纯流量卡能插随身wifi吗| 少妇猛男粗大的猛烈进出视频| 自线自在国产av| 成人亚洲欧美一区二区av| 操美女的视频在线观看| 老司机在亚洲福利影院| 国产不卡av网站在线观看| 丝袜人妻中文字幕| 99热全是精品| 国产在线免费精品| 中文字幕色久视频| 一级黄色大片毛片| 国产三级黄色录像| 18禁裸乳无遮挡动漫免费视频| 久久精品亚洲熟妇少妇任你| 1024香蕉在线观看| 精品熟女少妇八av免费久了| 亚洲自偷自拍图片 自拍| 国产在线观看jvid| 久久精品成人免费网站| 9色porny在线观看| 曰老女人黄片| 亚洲国产成人一精品久久久| 国产成人系列免费观看| 免费看不卡的av| 亚洲色图 男人天堂 中文字幕| 脱女人内裤的视频| 亚洲少妇的诱惑av| 日韩人妻精品一区2区三区| 国产男女内射视频| 亚洲精品成人av观看孕妇| 午夜福利,免费看| 久久性视频一级片| 国产99久久九九免费精品| 久久人人97超碰香蕉20202| 黑丝袜美女国产一区| 国产精品亚洲av一区麻豆| 亚洲自偷自拍图片 自拍| 久久亚洲精品不卡| 成人黄色视频免费在线看| 久久久久久久久免费视频了| 老司机在亚洲福利影院| 新久久久久国产一级毛片| 18禁国产床啪视频网站| 久久精品国产亚洲av高清一级| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩国产mv在线观看视频| 少妇精品久久久久久久| 青草久久国产| 一区二区三区乱码不卡18| 91字幕亚洲| 国产熟女欧美一区二区| 欧美中文综合在线视频| 男人舔女人的私密视频| 天堂中文最新版在线下载| 男女边吃奶边做爰视频| 国产精品久久久久久精品古装| 美女高潮到喷水免费观看| 99热网站在线观看| 午夜激情av网站| 亚洲国产精品999| 国产极品粉嫩免费观看在线| 一本综合久久免费| 亚洲第一青青草原| 最新在线观看一区二区三区 | 成在线人永久免费视频| 伦理电影免费视频| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区| 国产成人av激情在线播放| 亚洲 国产 在线| 国产欧美日韩一区二区三区在线| 亚洲黑人精品在线| 亚洲,一卡二卡三卡| 国产日韩欧美亚洲二区| 亚洲av国产av综合av卡| 亚洲国产最新在线播放| 宅男免费午夜| 国产97色在线日韩免费| 超碰97精品在线观看| 色94色欧美一区二区| 国产一区亚洲一区在线观看| 久久久精品区二区三区| 日本a在线网址| 久久ye,这里只有精品| 9热在线视频观看99| 婷婷成人精品国产| 久久国产亚洲av麻豆专区| 99国产精品99久久久久| 夫妻午夜视频| 国产成人a∨麻豆精品| 亚洲 国产 在线| 九草在线视频观看| 国产精品成人在线| 国产亚洲午夜精品一区二区久久| www日本在线高清视频| 91成人精品电影| 国产av精品麻豆| 亚洲伊人久久精品综合| 国产av国产精品国产| 亚洲精品国产区一区二| 午夜福利,免费看| 国产黄频视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产片特级美女逼逼视频| 亚洲精品一区蜜桃| 女人久久www免费人成看片| 一区在线观看完整版| 国产日韩一区二区三区精品不卡| 手机成人av网站| 国产免费又黄又爽又色| 久久久国产精品麻豆| 99香蕉大伊视频| 亚洲精品久久午夜乱码| 久久国产精品影院| 日韩欧美一区视频在线观看| 国产福利在线免费观看视频| 各种免费的搞黄视频| netflix在线观看网站| 自线自在国产av| 国产免费现黄频在线看| 中文字幕人妻熟女乱码| 国产免费视频播放在线视频| av在线老鸭窝| 又大又黄又爽视频免费| 国产精品免费视频内射| 亚洲九九香蕉| 国产1区2区3区精品| 青春草视频在线免费观看| 美女主播在线视频| 国产亚洲av高清不卡| 久久女婷五月综合色啪小说| 人人澡人人妻人| 国产在线一区二区三区精| 日韩av在线免费看完整版不卡| 免费在线观看日本一区| kizo精华| 欧美日韩亚洲国产一区二区在线观看 | 精品人妻1区二区| 男女边摸边吃奶| 99久久综合免费| 午夜福利一区二区在线看| 18禁黄网站禁片午夜丰满| 水蜜桃什么品种好| 男女之事视频高清在线观看 | 男女无遮挡免费网站观看| 99精品久久久久人妻精品| 丝袜脚勾引网站| 亚洲av日韩精品久久久久久密 | 亚洲第一av免费看| 日本wwww免费看| 欧美精品亚洲一区二区| 国产一卡二卡三卡精品| 欧美日韩成人在线一区二区| 男女之事视频高清在线观看 | 一区二区三区精品91| av一本久久久久| 欧美日韩国产mv在线观看视频| 亚洲专区中文字幕在线| 国产精品二区激情视频| 电影成人av| www.熟女人妻精品国产| 一二三四在线观看免费中文在| 久久精品熟女亚洲av麻豆精品| 交换朋友夫妻互换小说| 成年女人毛片免费观看观看9 | 秋霞在线观看毛片| 国产熟女欧美一区二区| 91精品国产国语对白视频| 美女午夜性视频免费| 亚洲人成电影观看| 中文字幕av电影在线播放| 女人被躁到高潮嗷嗷叫费观| 国产亚洲精品久久久久5区| 熟女少妇亚洲综合色aaa.| 亚洲国产精品成人久久小说| av在线老鸭窝| 久久99一区二区三区| 欧美另类一区| 欧美亚洲日本最大视频资源| 亚洲精品美女久久av网站| 在线av久久热| 激情视频va一区二区三区| 久久久国产精品麻豆| 亚洲国产中文字幕在线视频| 女人被躁到高潮嗷嗷叫费观| 一二三四社区在线视频社区8| 操美女的视频在线观看| 在线 av 中文字幕| 国产欧美日韩一区二区三 | 肉色欧美久久久久久久蜜桃| 免费久久久久久久精品成人欧美视频| 久久人人爽人人片av| 精品人妻熟女毛片av久久网站| 免费在线观看黄色视频的| 一本久久精品| 一区二区日韩欧美中文字幕| 超色免费av| 久久 成人 亚洲| av在线app专区| 国产成人精品久久二区二区91| 国产高清videossex| 欧美日韩亚洲高清精品| 美女脱内裤让男人舔精品视频| 91麻豆av在线| 亚洲欧美日韩另类电影网站| 老汉色av国产亚洲站长工具| 成人国产av品久久久| 天天躁夜夜躁狠狠久久av| 激情五月婷婷亚洲| 国产亚洲av高清不卡| 女人久久www免费人成看片| 婷婷丁香在线五月| 久久久久国产精品人妻一区二区| 国产av一区二区精品久久| 国产精品久久久久久人妻精品电影 | 热re99久久精品国产66热6| 久久久久精品人妻al黑| 女警被强在线播放| 精品久久久精品久久久| 国产成人系列免费观看| 久热这里只有精品99| 亚洲av美国av| 国产极品粉嫩免费观看在线| 丰满少妇做爰视频| 亚洲,欧美精品.| 免费不卡黄色视频| 午夜两性在线视频| 亚洲国产日韩一区二区| 婷婷色综合大香蕉| 国产精品久久久久成人av| 午夜福利免费观看在线| 国产1区2区3区精品| 高清视频免费观看一区二区| 两个人看的免费小视频| 欧美老熟妇乱子伦牲交| 国产视频首页在线观看| 久久毛片免费看一区二区三区| 亚洲精品乱久久久久久| 亚洲国产毛片av蜜桃av| 国产麻豆69| 黄色 视频免费看| 久久精品亚洲av国产电影网| 亚洲成人免费电影在线观看 | 美女国产高潮福利片在线看| 国产福利在线免费观看视频| 亚洲国产精品成人久久小说| 老鸭窝网址在线观看| 欧美 亚洲 国产 日韩一| 亚洲人成77777在线视频| av在线播放精品| 亚洲精品美女久久久久99蜜臀 | 99精国产麻豆久久婷婷| 欧美乱码精品一区二区三区| 中文字幕人妻丝袜制服| 欧美性长视频在线观看| 国产女主播在线喷水免费视频网站| 一本久久精品| 日本a在线网址| 欧美人与性动交α欧美软件| 精品免费久久久久久久清纯 | 男女之事视频高清在线观看 | av国产久精品久网站免费入址| 91成人精品电影| 国产女主播在线喷水免费视频网站| 青青草视频在线视频观看| 国产野战对白在线观看| 亚洲激情五月婷婷啪啪| 啦啦啦视频在线资源免费观看| 国产又爽黄色视频| 国产精品一国产av| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 18禁裸乳无遮挡动漫免费视频| 一区二区三区激情视频| 免费观看人在逋| 黄片播放在线免费| 波多野结衣av一区二区av| 嫁个100分男人电影在线观看 | 色94色欧美一区二区| 黄色毛片三级朝国网站| 久久久国产精品麻豆| 精品高清国产在线一区| 亚洲av电影在线进入| 真人做人爱边吃奶动态| 色网站视频免费| 欧美乱码精品一区二区三区| 成人亚洲精品一区在线观看| 激情五月婷婷亚洲| 欧美日韩黄片免| 免费日韩欧美在线观看| 免费看不卡的av| 亚洲欧洲国产日韩| 久久国产精品影院| 免费人妻精品一区二区三区视频| 精品国产超薄肉色丝袜足j| 亚洲人成电影观看| 老司机靠b影院| 99精品久久久久人妻精品| 国产高清视频在线播放一区 | 午夜久久久在线观看| 午夜福利乱码中文字幕| 久久午夜综合久久蜜桃| 国产一区二区三区综合在线观看| 国产成人a∨麻豆精品| 波野结衣二区三区在线| 婷婷色综合www| 国产熟女欧美一区二区| 国产精品久久久久久精品电影小说| 天天操日日干夜夜撸| 九色亚洲精品在线播放| 国产成人影院久久av| 久久久精品区二区三区| 操出白浆在线播放| 97在线人人人人妻| 欧美日韩亚洲国产一区二区在线观看 | 久久精品国产亚洲av高清一级| 在线 av 中文字幕| 丝袜人妻中文字幕| 晚上一个人看的免费电影| 又大又黄又爽视频免费| 国产亚洲av片在线观看秒播厂| 欧美黄色淫秽网站| 巨乳人妻的诱惑在线观看| 亚洲综合色网址| 成年动漫av网址| 欧美黄色淫秽网站| 老汉色∧v一级毛片| 香蕉丝袜av| 色视频在线一区二区三区| 岛国毛片在线播放| 青青草视频在线视频观看| 一区福利在线观看| 亚洲av男天堂| 精品亚洲成国产av| √禁漫天堂资源中文www| 1024香蕉在线观看| 日日爽夜夜爽网站| 欧美激情高清一区二区三区| 水蜜桃什么品种好| 国产熟女欧美一区二区| 热re99久久国产66热| 亚洲国产日韩一区二区| 天天躁日日躁夜夜躁夜夜| 欧美xxⅹ黑人| 久久国产精品影院| 99香蕉大伊视频| 亚洲国产精品国产精品| 午夜激情久久久久久久| 亚洲成人国产一区在线观看 | 99国产综合亚洲精品| av不卡在线播放| 亚洲黑人精品在线| 国产无遮挡羞羞视频在线观看| 精品熟女少妇八av免费久了| 久久久国产欧美日韩av| www.精华液| 欧美激情 高清一区二区三区| xxxhd国产人妻xxx| 波野结衣二区三区在线| 亚洲国产精品成人久久小说| 两性夫妻黄色片| 日韩制服丝袜自拍偷拍| 亚洲综合色网址| 久久天堂一区二区三区四区| 日日摸夜夜添夜夜爱| 国产高清不卡午夜福利| 黄色视频在线播放观看不卡| 成年人黄色毛片网站| 嫩草影视91久久| 18禁裸乳无遮挡动漫免费视频| 一区福利在线观看| 亚洲欧美日韩高清在线视频 | 黄色视频不卡| 欧美人与性动交α欧美软件| 丝袜人妻中文字幕| 亚洲 国产 在线| 精品久久蜜臀av无| 国产又色又爽无遮挡免| 久久国产亚洲av麻豆专区| 国产亚洲精品久久久久5区| 亚洲欧美一区二区三区国产| 别揉我奶头~嗯~啊~动态视频 | 51午夜福利影视在线观看| 九色亚洲精品在线播放| 中文字幕人妻熟女乱码| 亚洲第一av免费看| 久久热在线av| 亚洲人成77777在线视频| 美女国产高潮福利片在线看| 热99久久久久精品小说推荐| 9色porny在线观看| 色网站视频免费| 曰老女人黄片| 欧美人与性动交α欧美精品济南到| 亚洲伊人色综图| 国产男女内射视频| 在线精品无人区一区二区三| 咕卡用的链子| 1024视频免费在线观看| 午夜视频精品福利| 啦啦啦在线观看免费高清www| 亚洲视频免费观看视频| 午夜日韩欧美国产| 一级片免费观看大全| 免费看av在线观看网站| 成人国语在线视频| a级片在线免费高清观看视频| 久久 成人 亚洲| 亚洲欧美一区二区三区黑人| 亚洲精品在线美女| 欧美黄色淫秽网站| 亚洲第一青青草原| 女人爽到高潮嗷嗷叫在线视频| 伊人亚洲综合成人网| 国产一级毛片在线| 精品久久蜜臀av无| 亚洲成av片中文字幕在线观看| 国产成人精品久久二区二区免费| 久久久久久久久久久久大奶| 老司机靠b影院| 爱豆传媒免费全集在线观看| 亚洲男人天堂网一区| 欧美乱码精品一区二区三区| 18在线观看网站| 乱人伦中国视频| 亚洲国产中文字幕在线视频| 国产高清videossex| 性少妇av在线| 美女中出高潮动态图| 国产成人啪精品午夜网站| 一二三四社区在线视频社区8| 美国免费a级毛片| 日韩人妻精品一区2区三区| 黄色视频在线播放观看不卡| 亚洲黑人精品在线| 国产精品免费视频内射| av视频免费观看在线观看| 黄色视频不卡| 色视频在线一区二区三区| 国产精品国产av在线观看| 成人三级做爰电影| 午夜免费成人在线视频| 亚洲精品久久久久久婷婷小说| 国产亚洲av高清不卡| 1024视频免费在线观看| 亚洲专区国产一区二区| 久久免费观看电影| 少妇 在线观看| 中国美女看黄片| 午夜两性在线视频| 91国产中文字幕| 国产精品久久久久久精品古装| 别揉我奶头~嗯~啊~动态视频 | 欧美久久黑人一区二区| 在线亚洲精品国产二区图片欧美| 国产男人的电影天堂91| 九草在线视频观看| 黄色视频在线播放观看不卡| 人人妻人人澡人人看| 国产av一区二区精品久久| 亚洲av片天天在线观看| 国产精品国产三级专区第一集| 国产爽快片一区二区三区| e午夜精品久久久久久久| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 老熟女久久久| 国产精品亚洲av一区麻豆| 日本欧美视频一区| 欧美国产精品一级二级三级| 久久精品人人爽人人爽视色| 中文字幕人妻熟女乱码| av欧美777| 中文字幕高清在线视频| 18在线观看网站| 丁香六月天网| 日本av手机在线免费观看| 大型av网站在线播放| 久久久久久人人人人人| 多毛熟女@视频| 色94色欧美一区二区| 99热全是精品| 少妇的丰满在线观看| 国产精品二区激情视频| 成年人午夜在线观看视频| 久久精品久久久久久噜噜老黄| 日韩精品免费视频一区二区三区| 国产成人免费无遮挡视频| 麻豆av在线久日| 亚洲成人免费av在线播放| 日本av手机在线免费观看| 交换朋友夫妻互换小说| 久久精品国产a三级三级三级| 国产精品免费大片| 夫妻性生交免费视频一级片| 亚洲欧美激情在线| 国产精品人妻久久久影院| 岛国毛片在线播放| 侵犯人妻中文字幕一二三四区| 桃花免费在线播放| 亚洲成av片中文字幕在线观看| 男人爽女人下面视频在线观看| 女警被强在线播放| kizo精华| 国产精品熟女久久久久浪| 国产欧美日韩综合在线一区二区| 国产又色又爽无遮挡免| 老司机深夜福利视频在线观看 | 亚洲,一卡二卡三卡| 日韩一本色道免费dvd| 91成人精品电影| 啦啦啦中文免费视频观看日本| 欧美精品人与动牲交sv欧美| 国产亚洲av高清不卡| 观看av在线不卡| 在线观看人妻少妇| www.熟女人妻精品国产| 日本黄色日本黄色录像| 曰老女人黄片| 两人在一起打扑克的视频| 国产欧美日韩精品亚洲av| 欧美日韩av久久| 女人精品久久久久毛片| 成年人免费黄色播放视频| 亚洲欧美日韩另类电影网站| 日韩中文字幕视频在线看片| 久久精品人人爽人人爽视色| 热re99久久精品国产66热6| 国产成人啪精品午夜网站| 久久中文字幕一级| 飞空精品影院首页| 国产人伦9x9x在线观看| 满18在线观看网站|