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

    基于加速度反饋的單輸入時(shí)不變線性離散系統(tǒng)極點(diǎn)配置

    2015-10-28 11:33:55潘金文
    中國(guó)機(jī)械工程 2015年8期
    關(guān)鍵詞:離散系統(tǒng)標(biāo)準(zhǔn)型閉環(huán)

    梁 青 王 珍 潘金文 彭 程

    中國(guó)科學(xué)技術(shù)大學(xué),合肥,230026

    基于加速度反饋的單輸入時(shí)不變線性離散系統(tǒng)極點(diǎn)配置

    梁青王珍潘金文彭程

    中國(guó)科學(xué)技術(shù)大學(xué),合肥,230026

    討論了單輸入時(shí)不變線性離散系統(tǒng)用加速度信號(hào)進(jìn)行極點(diǎn)配置的可行性。首先, 給出了用狀態(tài)差分反饋進(jìn)行極點(diǎn)配置的定理和算法; 其次, 用3個(gè)仿真實(shí)例驗(yàn)證了該算法的有效性;最后, 在隔振器實(shí)驗(yàn)平臺(tái)上驗(yàn)證了用加速度信號(hào)進(jìn)行反饋控制的可行性。

    線性離散系統(tǒng); 加速度反饋; 狀態(tài)微分反饋; 狀態(tài)差分反饋; 極點(diǎn)配置

    0 引言

    極點(diǎn)配置是線性系統(tǒng)設(shè)計(jì)最重要的工具之一[1-4]。用狀態(tài)反饋實(shí)現(xiàn)極點(diǎn)配置是用狀態(tài)空間描述的線性系統(tǒng)(連續(xù)或離散)控制器設(shè)計(jì)的有效手段,學(xué)者們對(duì)此已經(jīng)做了大量研究工作[1-2,4-7],在這些研究中,假設(shè)系統(tǒng)的全部狀態(tài)可測(cè),即全狀態(tài)反饋。對(duì)于狀態(tài)不能直接可測(cè)的系統(tǒng),通過構(gòu)建全維或降維漸近狀態(tài)觀測(cè)器得到系統(tǒng)狀態(tài)的觀測(cè)值,并用觀測(cè)狀態(tài)實(shí)現(xiàn)狀態(tài)反饋,增加系統(tǒng)的階次。對(duì)于線性連續(xù)系統(tǒng),一種特殊的狀態(tài)反饋——狀態(tài)微分反饋近年來受到學(xué)者們的廣泛關(guān)注,其實(shí)質(zhì)是利用加速度信號(hào)實(shí)現(xiàn)反饋,這方面已有較多的研究成果[8-15]。文獻(xiàn)[12-13]分別討論了單輸入和多輸入線性系統(tǒng)極點(diǎn)配置問題的求解算法。文獻(xiàn)[14]給出了用狀態(tài)微分反饋以及輸出微分反饋設(shè)計(jì)LQR(linear quadratic regulator)的方法。而對(duì)于線性離散系統(tǒng),利用加速信號(hào)實(shí)現(xiàn)狀態(tài)反饋的研究還未見報(bào)道。

    在經(jīng)典控制理論中,微分可以增加系統(tǒng)阻尼,為實(shí)現(xiàn)預(yù)期控制目標(biāo),微分反饋有時(shí)是必不可少的[8],而研究狀態(tài)微分反饋的另一個(gè)動(dòng)機(jī)來源于減振機(jī)械系統(tǒng)[9-11],其主要的傳感器是加速度計(jì)。由加速度可以很容易地在可接受的誤差范圍內(nèi)獲得速度信號(hào),但由于加速度計(jì)的測(cè)量誤差及測(cè)量噪聲干擾,由加速度獲得位移信號(hào)會(huì)存在較大的誤差,因而可用于反饋的信號(hào)只有加速度和速度信號(hào),它們分別對(duì)應(yīng)于機(jī)械系統(tǒng)速度和位移的微分,即狀態(tài)微分。在離散系統(tǒng)框架下,加速度和速度對(duì)應(yīng)于狀態(tài)的差分。

    本文主要研究單輸入時(shí)不變線性離散系統(tǒng)的加速度反饋極點(diǎn)配置問題。首先對(duì)問題進(jìn)行數(shù)學(xué)描述,使用線性變換將一般線性離散系統(tǒng)變換為具有能控標(biāo)準(zhǔn)型的系統(tǒng);其次討論了能控標(biāo)準(zhǔn)型系統(tǒng)狀態(tài)差分反饋極點(diǎn)配置問題以及一般線性離散系統(tǒng)狀態(tài)差分反饋極點(diǎn)配置問題的解,給出了狀態(tài)差分反饋極點(diǎn)配置定理。再次,給出了減振器系統(tǒng)、球桿系統(tǒng)以及雙積分系統(tǒng)的狀態(tài)差分反饋控制仿真結(jié)果,驗(yàn)證了所提算法的正確性和有效性。最后,在實(shí)際隔振器實(shí)驗(yàn)平臺(tái)上驗(yàn)證了用加速度信號(hào)進(jìn)行反饋控制的可行性。

    1 加速度反饋極點(diǎn)配置問題的描述與能控標(biāo)準(zhǔn)型系統(tǒng)

    1.1問題描述

    考慮能控的單輸入時(shí)不變線性離散系統(tǒng)為

    x(k+1)=Gx(k)+hu(k)x(0)=x0

    (1)

    其中,x(k)是狀態(tài)變量,x(k)∈Rn×1;u(k)是控制量,u(k)∈R1×1;G是系統(tǒng)矩陣,G∈Rn×n;h是控制增益向量,h∈Rn×1。用加速度的采樣值實(shí)現(xiàn)反饋控制,即

    u(k)=-ksa(k)

    (2)

    其中,ks是加速度反饋增益向量,ks∈R1×n;a(k)為加速度的采樣值以及計(jì)算得到的速度采樣值,a(k)∈Rn×1。為研究其極點(diǎn)配置問題,有

    u(k)=-ks(x(k+1)-x(k))/Ts=

    -k(x(k+1)-x(k))

    (3)

    其中,Ts是采樣周期;k是差分反饋增益,k=(k0,k1,…,kn-1),并有ks=kTs。則極點(diǎn)配置問題變成狀態(tài)差分增益向量k的求解問題。系統(tǒng)矩陣G的特征多項(xiàng)式為

    det(zI-G)=zn+g0+g1z+g2z2+…+gn-1zn-1

    (4)

    其中,z為變換算子,I為單位矩陣,I∈Rn×n。定義g(g∈R1×n)為系統(tǒng)矩陣G的特征多項(xiàng)式系數(shù),g=[g0g1…gn-1]。設(shè)計(jì)的目標(biāo)是通過式(3)的狀態(tài)差分反饋將閉環(huán)系統(tǒng)的極點(diǎn)配置到期望的位置上,以穩(wěn)定系統(tǒng)并使系統(tǒng)具有期望的動(dòng)態(tài)特性。此時(shí),閉環(huán)系統(tǒng)為

    x(k+1)=Gx(k)-hk(x(k+1)-x(k))

    (5)

    通過選擇k使矩陣I+hk滿秩,則閉環(huán)系統(tǒng)可以寫成

    x(k+1)=(I+hk)-1(G+hk)x(k)

    (6)

    閉環(huán)系統(tǒng)的特征多項(xiàng)式為

    det(zI-Gc)=det(zI-(I+hk)-1(G+hk))

    (7)

    Gc=(I+hk)-1(G+hk)

    其中,Gc是閉環(huán)系統(tǒng)的系統(tǒng)矩陣。設(shè)計(jì)的問題就是通過尋找狀態(tài)差分反饋增益向量k,使得閉環(huán)系統(tǒng)(式(6))的極點(diǎn)在期望的位置上。

    1.2能控標(biāo)準(zhǔn)型系統(tǒng)

    采用文獻(xiàn)[12]的算法,對(duì)原系統(tǒng)(式(1))進(jìn)行非奇異狀態(tài)變換:

    (8)

    其中,z(k)是變換后的狀態(tài)向量,z(k)∈Rn×1;Q-1是變換矩陣,Q-1∈Rn×n,設(shè)變換后的系統(tǒng)矩陣為Gf,控制增益向量為hf,則有

    (9)

    變換矩陣的選取如下:

    Q-1=[qq Gq G2…q Gn-1]T

    (10)

    (11)

    (12)

    其中,向量q∈R1×n;en是單位向量;Mc是原系統(tǒng)(式(1))的能控性矩陣。

    變換后的系統(tǒng)是下友型能控標(biāo)準(zhǔn)型系統(tǒng):

    z(k+1)=Gfz(k)+hfu(k)

    (13)

    該下友型能控標(biāo)準(zhǔn)型系統(tǒng)能夠簡(jiǎn)化基于狀態(tài)差分反饋控制的極點(diǎn)配置問題的求解過程。

    2 極點(diǎn)配置問題的解

    2.1下友型能控標(biāo)準(zhǔn)型系統(tǒng)問題的解

    系統(tǒng)(式(13))的狀態(tài)差分反饋為

    u(k)=-l(z(k+1)-z(k))

    (14)

    l=[l0l1…ln]

    式中,l為狀態(tài)差分反饋向量;li為第i個(gè)狀態(tài)差分反饋系數(shù),i=0,1,…,n。

    其中,重寫控制律(式(14))為

    u(k)=-l0(z1(k+1)-z1(k))-l1(z2(k+1)-

    z2(k))-…-ln-1(zn(k+1)-zn(k))

    (15)

    由式(13)易得

    (16)

    將式(15)代入式(16)的最后一式并利用狀態(tài)之間的關(guān)系,可得

    zn(k+1)=-g0z1(k)-g1z2(k)-…-gn-1zn(k)+

    l0z1(k)-(l0-l1)z2(k)-(l1-l2)z3(k)-

    …-ln-1zn(k+1)=-(g0-l0)z1(k)-

    (g1+l0-l1)z2(k)-(g2+l1-l2)z3(k)-

    …-(gn-1+ln-2-ln-1)zn(k)-ln-1zn(k+1)

    (17)

    若1+ln-1≠0,則式(17)可以寫成標(biāo)準(zhǔn)狀態(tài)方程:

    zn(k+1)=-(1+ln-1)-1[(g0-l0)z1(k)+

    (g1+l0-l1)z2(k)+…+(gn-1+ln-2-ln-1)zn(k)]

    (18)

    則下友型標(biāo)準(zhǔn)型的狀態(tài)差分反饋閉環(huán)系統(tǒng)為

    z(k+1)=Gcfz(k)

    (19)

    Gcf=

    由閉環(huán)系統(tǒng)(式(6))和狀態(tài)變換方程(式(8)),可以得到下友型能控標(biāo)準(zhǔn)型的閉環(huán)系統(tǒng)矩陣與原閉環(huán)系統(tǒng)矩陣的關(guān)系:

    (20)

    下友型能控標(biāo)準(zhǔn)型狀態(tài)差分反饋閉環(huán)系統(tǒng)(式(19))的特征多項(xiàng)式為

    det(zI-Gcf)=(1+ln-1)-1[g0-l0+(g1+l0-l1)z+

    …+(gn-1+ln-2-ln-1)zn-1]+zn

    (21)

    期望極點(diǎn)可以由特征方程的根給出,也可以由特征多項(xiàng)式系數(shù)給出,即

    d(z)=(z-λ1)(z-λ2)…(z-λn)=

    zn+d0+d1z+d2z2+…+dn-1zn-1

    (22)

    其中,特征值λi∈C1×1,i=1,2,…,n,若有復(fù)根,則必須以共軛形式成對(duì)出現(xiàn)。期望特征多項(xiàng)式的系數(shù)d∈R1×n且d=[d0d1…dn-1]。由式(21)和式(22)關(guān)于z的多項(xiàng)式對(duì)應(yīng)項(xiàng)系數(shù)相等可以得到

    (23)

    因?yàn)?+ln-1≠0,整理可得

    (24)

    寫成矩陣形式如下:

    DlT=gT-dT

    (25)

    其中,系數(shù)矩陣D為

    矩陣D的行列式為

    (26)

    其中,e∈Rn×1且e=[111…1]T。

    式(25)的解分以下兩種情況討論:

    (1)det D≠0。即矩陣D是非奇異的,由式(25)可得下友型能控標(biāo)準(zhǔn)型極點(diǎn)配置問題的解:

    l=(g-d)D-T

    (27)

    事實(shí)上,不需要矩陣求逆,也可直接得到式(24)的解。將式(24)各式的左邊依次相加,右邊依次相加可得

    (1+de)ln-1=(g-d)e

    (28)

    解得

    ln-1=(g-d)e/(1+de)

    (29)

    求得ln-1后再倒序求解式(24),即求解順序?yàn)閘0,l1,…,ln-2,可得

    (30)

    γ=1+ln-1

    (2)det D=0。即矩陣D是奇異的,為使式(28)有解,必須滿足

    ge=de

    (31)

    這就表明,ge+1=0和de+1=0必須同時(shí)成立。由式(28)可知,此時(shí)ln-1是自由的,只需選擇ln-1≠-1即可。此時(shí)式(24)的解為

    (32)

    若1+ge≠0且1+de≠0,則

    (33)

    式(33)表明,原系統(tǒng)的所有特征根均不為1,這意味著對(duì)于任意的輸入原系統(tǒng)都有一個(gè)平衡點(diǎn)xeq。在式(1)中令x(k+1)=x(k)=xeq,解得

    xeq=(I-G)-1hu

    (34)

    對(duì)于期望的閉環(huán)系統(tǒng),也有相同的結(jié)論。這表明,狀態(tài)差分反饋無需考慮系統(tǒng)有無平衡點(diǎn)。

    2.2一般系統(tǒng)的解

    上文針對(duì)具有下友型能控標(biāo)準(zhǔn)型的單輸入時(shí)不變線性離散系統(tǒng),求得了基于狀態(tài)差分反饋的問題的解,下面求一般單輸入時(shí)不變線性離散系統(tǒng)的解。

    將式(8)代入式(14)可得

    u(k)=-l(z(k+1)-z(k))=

    -l Q-1(x(k+1)-x(k))

    (35)

    對(duì)比式(3)和式(35)可得原系統(tǒng)(式(1))的狀態(tài)差分反饋控制問題的解:

    k=l Q-1

    (36)

    推導(dǎo)狀態(tài)差分反饋增益向量k的具體表達(dá)式。有

    l=(g-γd)Mr

    (37)

    其中,Mr為上三角矩陣。得到狀態(tài)差分反饋增益向量:

    k=(g-γd)MrQ-1

    (38)

    加速度反饋增益向量為

    ks=Ts(g-γd)MrQ-1

    (39)

    由上述討論得出以下定理。

    定理1若系統(tǒng)(式(1))能控且不含極點(diǎn)z=1, 則可通過狀態(tài)差分反饋(式(3))在z平面任意配置極點(diǎn)到期望的位置上。

    定理2若系統(tǒng)(式(1))能控且含極點(diǎn)z=1,則可通過狀態(tài)差分反饋(式(3))在z平面任意配置極點(diǎn)到至少含有以z=1為其中一個(gè)極點(diǎn)的期望位置上。

    3 仿真實(shí)例

    例1考慮減振機(jī)械系統(tǒng),其動(dòng)態(tài)方程為

    其中,各參數(shù)的具體意義見文獻(xiàn)[12]。易知原系統(tǒng)是一個(gè)開環(huán)漸近穩(wěn)定的系統(tǒng), 但系統(tǒng)達(dá)到平衡狀態(tài)的時(shí)間較長(zhǎng)。設(shè)計(jì)狀態(tài)差分反饋使閉環(huán)系統(tǒng)穩(wěn)定且零輸入響應(yīng)的暫態(tài)過程時(shí)間縮短。取采樣周期Ts=0.5 ms,離散化后的系統(tǒng)矩陣和控制量增益分別為

    h=[00-5×10-65×10-5]T

    閉環(huán)系統(tǒng)期望極點(diǎn)

    λ1,2=exp(Ts(-5±65i))

    初始狀態(tài)取

    x0=[0.050.050.20.2]T

    解得加速度反饋增益

    ks=[260.3388294.1277-1.2023-0.2872]

    由圖1可以看出,閉環(huán)系統(tǒng)的瞬態(tài)響應(yīng)明顯優(yōu)于開環(huán)系統(tǒng);由圖2可以看出,加速度信號(hào)(a3(t),a4(t))的瞬態(tài)值較大,而速度信號(hào)(a1(t),a2(t))的瞬態(tài)值不至于過大;圖3所示為控制信號(hào)的變化曲線。

    1.閉環(huán)系統(tǒng) 2.原系統(tǒng)圖1 原系統(tǒng)和閉環(huán)系統(tǒng)零輸入響應(yīng)(例1)

    圖2 狀態(tài)差分的理論值(例1)

    圖3 控制信號(hào)(例1)

    例2考慮球桿系統(tǒng),動(dòng)態(tài)方程為

    其中,各u(t)參數(shù)的具體意義見文獻(xiàn)[12]。易知原系統(tǒng)是一個(gè)開環(huán)不穩(wěn)定的系統(tǒng),設(shè)計(jì)狀態(tài)差分反饋使得閉環(huán)系統(tǒng)漸近穩(wěn)定且暫態(tài)過程盡可能短。取采樣周期Ts=5 ms,離散化后的系統(tǒng)矩陣和控制量增益分別為

    h=[00.004 9500]T

    閉環(huán)系統(tǒng)期望極點(diǎn)為

    λ1,2=exp(Ts(-3±i))

    λ3=-2Tsλ4=-5Ts

    解得加速度反饋增益

    ks=[-0.1064-1.010512.85176.1787]

    初始狀態(tài)為

    x0=[0.1 0.1 0.01 0.01]T

    由圖4可以看出,閉環(huán)系統(tǒng)漸近穩(wěn)定且能較快到達(dá)平衡狀態(tài);圖5所示為狀態(tài)差分的理論值;圖6所示為控制信號(hào)的變化曲線。仿真結(jié)果表明,對(duì)于開環(huán)不穩(wěn)定的系統(tǒng),本文方法也是有效的。

    圖4 閉環(huán)系統(tǒng)零輸入響應(yīng)(例2)

    圖5 狀態(tài)差分的理論值(例2)

    圖6 控制信號(hào)(例2)

    例3考慮雙積分系統(tǒng),動(dòng)態(tài)方程為

    取采樣周期Ts=5 ms,離散化后的系統(tǒng)矩陣和控制量增益分別為

    圖7~圖9的意義與例2相同,從圖7可以看出,閉環(huán)系統(tǒng)仍然具有一階積分特性,即狀態(tài)差分反饋不能完全去除原系統(tǒng)在z=1處的極點(diǎn),從而驗(yàn)證了定理2。

    圖7 閉環(huán)系統(tǒng)零輸入響應(yīng)(例3)

    圖8 狀態(tài)差分的理論值(例3)

    圖9 控制信號(hào)(例3)

    4 實(shí)驗(yàn)結(jié)果

    本課題組一直致力于隔振器方面的研究。圖10所示是本課題組自主研發(fā)的磁懸浮隔振系統(tǒng)[16]。

    圖10 磁懸浮隔振器實(shí)驗(yàn)平臺(tái)原理圖

    由圖10可以看出,負(fù)載與基礎(chǔ)的加速度可以分別由加速度傳感器1和2直接測(cè)量獲得,可以直接用于反饋控制。為驗(yàn)證加速度反饋的可行性,振動(dòng)主動(dòng)控制的目標(biāo)是當(dāng)負(fù)載或基礎(chǔ)偏離自身平衡位置時(shí),電磁作動(dòng)器按照一定的規(guī)律產(chǎn)生電磁力,從而使負(fù)載和基礎(chǔ)分別回到自身的平衡位置。在對(duì)電磁作動(dòng)器輸出的電磁力與基礎(chǔ)及速度的關(guān)系進(jìn)行建模時(shí),選擇負(fù)載位移、負(fù)載速度、基礎(chǔ)位移、基礎(chǔ)速度作為狀態(tài),則由負(fù)載加速度得到的負(fù)載速度信號(hào)、負(fù)載加速度,由基礎(chǔ)加速度得到的基礎(chǔ)速度信號(hào)以及基礎(chǔ)加速度信號(hào)就對(duì)應(yīng)于狀態(tài)的微分,即可直接利用加速度信號(hào)實(shí)現(xiàn)反饋控制而無需構(gòu)建狀態(tài)觀測(cè)器。上述系統(tǒng)可以用圖11所示模型來描述。

    圖11 電磁作動(dòng)器作用下磁懸浮隔振系統(tǒng)原理圖

    系統(tǒng)的狀態(tài)空間表達(dá)式為

    其中,u為電磁力,s1為負(fù)載位移,s2為基礎(chǔ)位移。磁懸浮隔振器中各元件的標(biāo)稱值見表1。

    表1 磁懸浮隔振器元件標(biāo)稱值

    (a)負(fù)載的加速度測(cè)量值與速度重構(gòu)值

    (b)基礎(chǔ)的加速度測(cè)量值與速度重構(gòu)值圖12 系統(tǒng)無控制與施加控制時(shí)負(fù)載與基礎(chǔ)的加速度測(cè)量值與速度重構(gòu)值

    圖13 負(fù)載與基礎(chǔ)的相對(duì)位移和電磁力變化曲線

    初始時(shí),手動(dòng)使負(fù)載稍微偏離自身平衡位置,施加控制與不施加控制的變量對(duì)比見圖12和圖13??刂品椒ú捎帽疚乃岱椒?采樣周期取Ts=0.2 ms,閉環(huán)系統(tǒng)極點(diǎn)取λ1,2=exp(Ts(-5±15i)),λ3,4=exp(Ts(-7±10i))。由圖12可知,若不施加控制,負(fù)載和基礎(chǔ)都需要較長(zhǎng)時(shí)間才能回到平衡位置;而使用加速度信號(hào)作為反饋的控制,可以使系統(tǒng)盡快回到平衡位置,該結(jié)果驗(yàn)證了加速度反饋的有效性。

    5 結(jié)語

    本文從單輸入時(shí)不變線性系統(tǒng)的離散形式出發(fā),闡述了用加速度反饋配置系統(tǒng)極點(diǎn)的實(shí)際意義,通過轉(zhuǎn)化將其歸結(jié)為用狀態(tài)的差分作反饋的增益求取問題。詳細(xì)給出了下友型能控標(biāo)準(zhǔn)型系統(tǒng)基于狀態(tài)差分反饋極點(diǎn)配置問題的解的求取過程,并通過狀態(tài)變換得到一般時(shí)不變線性離散系統(tǒng)基于狀態(tài)差分反饋極點(diǎn)配置問題的解。由求解過程的討論得到基于狀態(tài)差分反饋的極點(diǎn)配置定理和算法。所舉的3個(gè)典型實(shí)例和隔振器實(shí)驗(yàn)都驗(yàn)證了本文方法的可行性和有效性。

    [1]DorfRC,BishopRH.ModernControlSystems[M].London:PrenticeHall,2011.

    [2]李嗣福.計(jì)算機(jī)控制基礎(chǔ)[M].合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社, 2006.

    [3]WonhamW.OnPoleAssignmentinMulti-inputControllableLinearSystems[J].IEEETransactionsonAutomaticControl, 1967,12(6): 660-665.

    [4]周兆敏.用極點(diǎn)配置方法設(shè)計(jì)滑??刂破鱗J].信息與控制,1988,17(4): 47-49.

    ZhouZhaomin.SlidingModeControllerDesignbyUsingPolePlacementmethod[J].InformationandControl,1988,17(4): 47-49.

    [5]段學(xué)超,袁俊,滿曰剛.兩輪自平衡機(jī)器人的動(dòng)力學(xué)建模與模糊進(jìn)化極點(diǎn)配置控制[J].信息與控制, 2013, 42(2):189-195.

    DuanXuechao,YuanJun,ManYuegang.DynamicalModelingandFuzzyEvolutionaryPolePlacementControlforaTwo-wheeledSelfbalancingRobot[J].InformationandControl, 2013, 42(2):189-195.

    [6]KautskyJ,NicholsNK,VanDP.RobustPoleAssignmentinLinearStateFeedback[J].InternationalJournalofControl, 1985,41(5):1129-1155.

    [7]GaoZunhai,YeXichen.Pole-placementforMulti-inputLinearSystemsbyStateFeedbackorState-derivativeFeedback[J].AdvancedMaterialsResearch, 2013,605/607:1765-1768.

    [8]LewisF,SyrmosV.AGeometricTheoryforDerivativeFeedback[J].IEEETransactionsonAutomaticControl, 1991, 36(9):1111-1116.

    [9]OlgacN,ElmaliH,HosekM,etal.ActiveVibrationControlofDistributedSystemsUsingDelayedResonatorwithAccelerationFeedback[J].JournalofDynamicSystems,MeasurementandControl, 1997, 119(3): 380-389.

    [10]KwakSK,WashingtonG,YedavalliRK.AccelerationFeedback-basedActiveandPassiveVibrationControlofLandingGearComponents[J].JournalofAerospaceEngineering, 2002,15(1): 1-9.

    [11]ChattwejeeS.VibrationControlbyRecursiveTime-delayedAccelerationFeedback[J].JournalofSoundandVibration, 2008,317(1): 67-90.

    [12]AbdelazizTHS,ValaekM.Pole-placementforSISOLinearSystemsbyState-derivativeFeedback[J].IEEProceedings:ControlTheoryandApplications, 2004, 151(4): 377-385.

    [13]AbdelazizTHS,ValaekM.DirectAlgorithmforPolePlacementbyState-derivativeFeedbackforMulti-inputLinearSystems-NonsingularCase[J].Kybernetika, 2005, 41(5): 637-660.

    [14]AbdelazizTHS,ValaekM.StateDerivativeFeedbackbyLQRforLinearTime-invariantSystems[C]//Proceedingsofthe16thTriennialWorldCongressofInternationalFederationofAutomaticControl.Prague, 2005:1-6.

    [15]AmriI,SoudaniD,BenrejebM.RobustState-derivativeFeedbackLMI-basedDesignsforTime-varyingDelaySystem[C]//Proceedingsofthe2011InternationalConferenceonCommunications,ComputingandControlApplications.Hammamet, 2011.

    [16]孟令雷.磁懸浮隔振器的建模與控制[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2009.

    (編輯陳勇)

    Pole Placement for Single Input Linear Time-invariant Discrete Systems by Acceleration Feedback

    Liang QingWang ZhenPan JinwenPeng Cheng

    University of Science and Technology of China,Hefei,230026

    The feasibility of pole placement by acceleration signal for single input time-invariant linear discrete systems was discussed. Firstly, the theorems and algorithms of pole placement by state-difference were provided. Secondly, the effectiveness of the proposed method was validated by three simulation examples. Finally, an experimental platform was setup to validate the feasibility of the proposed method.

    linear discrete system; acceleration feedback; state-derivative feedback; state-difference feedback; pole placement

    2013-09-16

    國(guó)家自然科學(xué)基金資助項(xiàng)目(61004017);中國(guó)科學(xué)技術(shù)大學(xué)青年創(chuàng)新基金資助項(xiàng)目(WK2100100016)

    TP273< class="emphasis_italic">DOI

    :10.3969/j.issn.1004-132X.2015.08.020

    梁青,男,1961年生。中國(guó)科學(xué)技術(shù)大學(xué)信息科學(xué)技術(shù)學(xué)院副教授。主要研究方向?yàn)轱w行器制導(dǎo)與控制、系統(tǒng)仿真、智能信息處理等。王珍,女,1989年生。中國(guó)科學(xué)技術(shù)大學(xué)信息科學(xué)技術(shù)學(xué)院碩士研究生。潘金文,男,1989年生。中國(guó)科學(xué)技術(shù)大學(xué)信息科學(xué)技術(shù)學(xué)院博士研究生。彭程,男,1978年生。中國(guó)科學(xué)技術(shù)大學(xué)信息科學(xué)技術(shù)學(xué)院副教授。

    猜你喜歡
    離散系統(tǒng)標(biāo)準(zhǔn)型閉環(huán)
    冪級(jí)數(shù)收斂半徑和收斂域的求解探討
    ——如何培養(yǎng)學(xué)生的創(chuàng)新思維
    關(guān)于非自治離散系統(tǒng)中敏感性的一些結(jié)論
    以代數(shù)思想為主線—線性代數(shù)和高等代數(shù)課程教學(xué)的相通與兼容
    單周期控制下雙輸入Buck變換器閉環(huán)系統(tǒng)設(shè)計(jì)
    黑龍江電力(2017年1期)2017-05-17 04:25:05
    “翻棋”
    雙閉環(huán)模糊控制在石化廢水處理中的研究
    標(biāo)準(zhǔn)型不高于五階若當(dāng)塊矩陣群的冪單性
    最優(yōu)價(jià)格與回收努力激勵(lì)的閉環(huán)供應(yīng)鏈協(xié)調(diào)
    一種基于全閉環(huán)實(shí)時(shí)數(shù)字物理仿真的次同步振蕩阻尼控制
    2-D連續(xù)-離散系統(tǒng)的迭代學(xué)習(xí)控制
    国产白丝娇喘喷水9色精品| 国产久久久一区二区三区| 色播亚洲综合网| 亚洲人成网站高清观看| 免费观看的影片在线观看| 97在线视频观看| 深夜精品福利| 日韩三级伦理在线观看| 国内揄拍国产精品人妻在线| 婷婷色综合大香蕉| 99热只有精品国产| 淫秽高清视频在线观看| 狠狠狠狠99中文字幕| 日韩高清综合在线| 麻豆国产av国片精品| 给我免费播放毛片高清在线观看| 精品免费久久久久久久清纯| 日韩强制内射视频| 日韩强制内射视频| 男女边吃奶边做爰视频| 99在线人妻在线中文字幕| 日本精品一区二区三区蜜桃| 国产精品一区二区三区四区久久| 欧美精品国产亚洲| 亚洲av成人av| 在线观看一区二区三区| 国产女主播在线喷水免费视频网站 | av福利片在线观看| 可以在线观看毛片的网站| 乱码一卡2卡4卡精品| 久久久久久久亚洲中文字幕| 美女被艹到高潮喷水动态| 国产欧美日韩一区二区精品| 色播亚洲综合网| 色播亚洲综合网| 亚洲四区av| 一区二区三区四区激情视频 | 超碰av人人做人人爽久久| av在线播放精品| 欧美日韩乱码在线| 国产黄色小视频在线观看| 国产黄色小视频在线观看| 亚洲高清免费不卡视频| 我要搜黄色片| 91久久精品国产一区二区三区| a级毛片a级免费在线| 搡老妇女老女人老熟妇| 日日摸夜夜添夜夜爱| 91久久精品国产一区二区成人| 男女视频在线观看网站免费| 精品久久久久久久末码| 久久久久九九精品影院| 女人十人毛片免费观看3o分钟| 日本爱情动作片www.在线观看 | 精品久久久久久久久久久久久| 精品久久久久久久久久久久久| 亚洲第一电影网av| 欧美国产日韩亚洲一区| 欧美高清性xxxxhd video| 久久精品国产亚洲av天美| 免费看av在线观看网站| 国产成人aa在线观看| 99热6这里只有精品| 欧美xxxx性猛交bbbb| 美女cb高潮喷水在线观看| 亚洲欧美清纯卡通| 精品欧美国产一区二区三| 欧美色视频一区免费| 精品人妻熟女av久视频| 简卡轻食公司| 亚洲欧美中文字幕日韩二区| 久久久久性生活片| 欧美中文日本在线观看视频| 亚洲欧美精品综合久久99| 欧美国产日韩亚洲一区| 精品久久久噜噜| 又爽又黄a免费视频| 淫妇啪啪啪对白视频| 国产 一区 欧美 日韩| 国产成人a区在线观看| 成人一区二区视频在线观看| 精品久久久久久久人妻蜜臀av| 日日摸夜夜添夜夜添小说| 九九爱精品视频在线观看| 国产精品亚洲美女久久久| 插逼视频在线观看| 一进一出抽搐gif免费好疼| av在线蜜桃| 成人三级黄色视频| 亚洲国产欧洲综合997久久,| a级毛片a级免费在线| 欧美成人a在线观看| 春色校园在线视频观看| 亚洲精品影视一区二区三区av| 成人性生交大片免费视频hd| 国产v大片淫在线免费观看| 国产 一区 欧美 日韩| 午夜精品一区二区三区免费看| av黄色大香蕉| 亚洲色图av天堂| 国产精品99久久久久久久久| 三级经典国产精品| 又黄又爽又免费观看的视频| 白带黄色成豆腐渣| 男女啪啪激烈高潮av片| 一级a爱片免费观看的视频| 国产av在哪里看| 伊人久久精品亚洲午夜| 蜜桃亚洲精品一区二区三区| 最近在线观看免费完整版| 成年女人毛片免费观看观看9| 黄色欧美视频在线观看| 人妻制服诱惑在线中文字幕| 国产精品乱码一区二三区的特点| 国产白丝娇喘喷水9色精品| 免费无遮挡裸体视频| 免费在线观看影片大全网站| 日产精品乱码卡一卡2卡三| 久久精品人妻少妇| 日日撸夜夜添| 亚洲性夜色夜夜综合| 91久久精品国产一区二区成人| 亚洲欧美清纯卡通| 中出人妻视频一区二区| 乱系列少妇在线播放| 天天躁日日操中文字幕| 美女内射精品一级片tv| 最新在线观看一区二区三区| 如何舔出高潮| 国产激情偷乱视频一区二区| 伊人久久精品亚洲午夜| 亚洲性久久影院| av.在线天堂| 午夜福利18| 国产精品,欧美在线| 中国国产av一级| 嫩草影院入口| 三级毛片av免费| 99久国产av精品| 内地一区二区视频在线| 成人美女网站在线观看视频| 午夜视频国产福利| 欧美最新免费一区二区三区| 村上凉子中文字幕在线| 亚州av有码| 日本免费a在线| 99在线视频只有这里精品首页| 老司机福利观看| 国产午夜精品论理片| 全区人妻精品视频| 亚洲真实伦在线观看| 一进一出好大好爽视频| 99热只有精品国产| 99九九线精品视频在线观看视频| 我要搜黄色片| 亚洲四区av| 99久国产av精品| 老女人水多毛片| av在线亚洲专区| 日韩精品青青久久久久久| 草草在线视频免费看| 国产一区亚洲一区在线观看| 69人妻影院| 99久国产av精品国产电影| 欧美丝袜亚洲另类| 精品午夜福利视频在线观看一区| 国产高清三级在线| 97超碰精品成人国产| 波野结衣二区三区在线| 亚洲精品日韩av片在线观看| 香蕉av资源在线| 麻豆国产97在线/欧美| 亚洲av熟女| 欧美xxxx性猛交bbbb| 欧美国产日韩亚洲一区| 中文字幕人妻熟人妻熟丝袜美| 老女人水多毛片| 国产aⅴ精品一区二区三区波| 乱码一卡2卡4卡精品| 国产aⅴ精品一区二区三区波| 蜜桃久久精品国产亚洲av| 青春草视频在线免费观看| 99久久精品一区二区三区| 18+在线观看网站| 不卡视频在线观看欧美| 欧美日韩综合久久久久久| 欧美最新免费一区二区三区| 少妇高潮的动态图| 亚洲中文字幕日韩| 国产精品久久久久久久电影| 国产免费一级a男人的天堂| 日韩人妻高清精品专区| 少妇被粗大猛烈的视频| 麻豆国产97在线/欧美| 久久综合国产亚洲精品| 日日干狠狠操夜夜爽| 高清毛片免费观看视频网站| 少妇人妻精品综合一区二区 | 久久欧美精品欧美久久欧美| 91午夜精品亚洲一区二区三区| 在线免费观看的www视频| 深爱激情五月婷婷| 高清毛片免费观看视频网站| 国产单亲对白刺激| 国产aⅴ精品一区二区三区波| 国产精品一区二区免费欧美| 亚洲在线观看片| 91精品国产九色| 人人妻人人看人人澡| 国产免费一级a男人的天堂| 一级a爱片免费观看的视频| 精品乱码久久久久久99久播| 最近视频中文字幕2019在线8| 久久韩国三级中文字幕| 日韩 亚洲 欧美在线| 又黄又爽又刺激的免费视频.| 免费av观看视频| 国产精品久久久久久久电影| 插逼视频在线观看| 亚洲乱码一区二区免费版| 我的女老师完整版在线观看| 久久久久久久久久黄片| 国产一区二区亚洲精品在线观看| 日韩av在线大香蕉| 男人和女人高潮做爰伦理| 中文字幕免费在线视频6| 赤兔流量卡办理| 亚洲国产高清在线一区二区三| 国产av在哪里看| 亚洲精品一区av在线观看| 日韩三级伦理在线观看| 日本免费a在线| 亚洲七黄色美女视频| 夜夜爽天天搞| 亚洲va在线va天堂va国产| 国产亚洲av嫩草精品影院| 久久久国产成人精品二区| 中文字幕人妻熟人妻熟丝袜美| 欧美绝顶高潮抽搐喷水| 22中文网久久字幕| 男人舔女人下体高潮全视频| 美女黄网站色视频| a级毛片免费高清观看在线播放| 露出奶头的视频| 综合色丁香网| 麻豆成人午夜福利视频| 精品一区二区三区视频在线| 99热6这里只有精品| 又黄又爽又刺激的免费视频.| 国产精品乱码一区二三区的特点| 国产蜜桃级精品一区二区三区| 午夜福利视频1000在线观看| 亚洲在线自拍视频| 男女啪啪激烈高潮av片| 午夜亚洲福利在线播放| 精品久久久久久久久久久久久| 成人特级黄色片久久久久久久| a级毛片a级免费在线| 晚上一个人看的免费电影| h日本视频在线播放| 国产私拍福利视频在线观看| 六月丁香七月| av在线老鸭窝| 1000部很黄的大片| 国产高清激情床上av| 一个人免费在线观看电影| 女的被弄到高潮叫床怎么办| 婷婷精品国产亚洲av| 亚洲av美国av| 午夜爱爱视频在线播放| 日韩成人伦理影院| 久久久久久久久大av| av女优亚洲男人天堂| 国产精品99久久久久久久久| 97碰自拍视频| 日本免费a在线| 国产大屁股一区二区在线视频| 精品一区二区三区视频在线观看免费| 亚洲成人中文字幕在线播放| 色播亚洲综合网| 国产精品一区二区免费欧美| 亚洲国产高清在线一区二区三| 亚洲真实伦在线观看| av在线播放精品| 国产综合懂色| 在线观看av片永久免费下载| 我要搜黄色片| 国产av一区在线观看免费| 亚洲一区二区三区色噜噜| 亚洲成人精品中文字幕电影| 舔av片在线| 春色校园在线视频观看| 国产精品久久久久久久久免| 啦啦啦观看免费观看视频高清| 久久久久国产精品人妻aⅴ院| 亚洲美女视频黄频| 精品免费久久久久久久清纯| 99在线视频只有这里精品首页| 中文资源天堂在线| 亚洲自拍偷在线| 18禁裸乳无遮挡免费网站照片| av专区在线播放| 日产精品乱码卡一卡2卡三| 男女视频在线观看网站免费| 中文字幕久久专区| 国产免费一级a男人的天堂| 国产黄a三级三级三级人| 亚洲av免费高清在线观看| av专区在线播放| 日日啪夜夜撸| 校园人妻丝袜中文字幕| 精品日产1卡2卡| 大香蕉久久网| 日日撸夜夜添| 成人av在线播放网站| 中文字幕av在线有码专区| 人人妻人人澡欧美一区二区| 我要搜黄色片| 久久精品国产亚洲av天美| 国内精品美女久久久久久| 免费看美女性在线毛片视频| 久久久a久久爽久久v久久| 久久亚洲精品不卡| 久久精品国产亚洲av天美| 亚洲成人中文字幕在线播放| 特大巨黑吊av在线直播| 看非洲黑人一级黄片| 99热这里只有是精品在线观看| 日韩av不卡免费在线播放| 日韩亚洲欧美综合| 99久久久亚洲精品蜜臀av| 一本久久中文字幕| 免费观看在线日韩| 免费观看的影片在线观看| 国产老妇女一区| 欧美成人a在线观看| 国内精品一区二区在线观看| 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看| 欧美性感艳星| 日韩一本色道免费dvd| 欧美成人a在线观看| 黄色配什么色好看| 一个人观看的视频www高清免费观看| 亚洲,欧美,日韩| 18禁在线播放成人免费| 日日摸夜夜添夜夜添av毛片| 久久久精品欧美日韩精品| 国产在线男女| 亚洲无线在线观看| 亚洲国产日韩欧美精品在线观看| 中国国产av一级| 在线国产一区二区在线| 99九九线精品视频在线观看视频| 成人永久免费在线观看视频| 欧美日韩乱码在线| 久久久久国产网址| 99久久无色码亚洲精品果冻| 在线观看66精品国产| av卡一久久| 国产黄片美女视频| 少妇猛男粗大的猛烈进出视频 | 少妇熟女aⅴ在线视频| av在线天堂中文字幕| 看黄色毛片网站| 丰满乱子伦码专区| 亚洲精华国产精华液的使用体验 | 好男人在线观看高清免费视频| 免费观看在线日韩| 亚洲国产精品成人久久小说 | АⅤ资源中文在线天堂| 免费人成在线观看视频色| 欧美精品国产亚洲| 国产在线男女| 欧美性感艳星| 久久久久国产网址| 午夜福利在线观看吧| 欧美不卡视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 日本成人三级电影网站| 久久九九热精品免费| 国产极品精品免费视频能看的| 毛片女人毛片| 午夜福利高清视频| 国产日本99.免费观看| 人人妻人人澡欧美一区二区| 国产精品久久久久久久久免| 久久这里只有精品中国| 国产又黄又爽又无遮挡在线| 天堂√8在线中文| 欧美一级a爱片免费观看看| 欧美精品国产亚洲| 赤兔流量卡办理| 噜噜噜噜噜久久久久久91| 亚洲欧美成人精品一区二区| 精品免费久久久久久久清纯| h日本视频在线播放| 国产精品久久视频播放| 国产午夜福利久久久久久| 国产精华一区二区三区| 成人毛片a级毛片在线播放| 午夜爱爱视频在线播放| 国产av在哪里看| 干丝袜人妻中文字幕| 国产一区二区三区av在线 | 人妻制服诱惑在线中文字幕| 亚洲婷婷狠狠爱综合网| 精品一区二区三区视频在线观看免费| 简卡轻食公司| 午夜影院日韩av| 干丝袜人妻中文字幕| 婷婷色综合大香蕉| 精品久久久久久久久久久久久| 国产精品,欧美在线| 精品久久久久久久久久免费视频| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 亚洲一级一片aⅴ在线观看| 亚洲美女搞黄在线观看 | 老女人水多毛片| 午夜福利在线观看免费完整高清在 | 日韩,欧美,国产一区二区三区 | 亚洲av二区三区四区| 亚洲无线在线观看| 国产精华一区二区三区| 日本熟妇午夜| 精品久久国产蜜桃| 校园春色视频在线观看| 欧美一区二区精品小视频在线| 国产成人福利小说| 尤物成人国产欧美一区二区三区| 国产三级在线视频| 老司机影院成人| 99久久九九国产精品国产免费| 久久久久久久午夜电影| 国产又黄又爽又无遮挡在线| 国产真实伦视频高清在线观看| 国产日本99.免费观看| 久久国产乱子免费精品| 久久久久精品国产欧美久久久| 亚洲精品色激情综合| 国产av麻豆久久久久久久| 舔av片在线| 成人高潮视频无遮挡免费网站| 性插视频无遮挡在线免费观看| 亚洲中文字幕日韩| 热99re8久久精品国产| 深夜a级毛片| 色av中文字幕| 黑人高潮一二区| 亚洲av中文字字幕乱码综合| 69人妻影院| 小说图片视频综合网站| 好男人在线观看高清免费视频| 神马国产精品三级电影在线观看| 中出人妻视频一区二区| 精品乱码久久久久久99久播| 国产精品一区二区免费欧美| 久久久久久大精品| 大香蕉久久网| 在现免费观看毛片| 日本与韩国留学比较| 成年av动漫网址| 国产精品,欧美在线| av专区在线播放| 日韩欧美免费精品| 精品不卡国产一区二区三区| 亚洲成人久久爱视频| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 国产成人福利小说| 日本熟妇午夜| 国产亚洲精品久久久com| 男女啪啪激烈高潮av片| 亚洲av电影不卡..在线观看| 天堂动漫精品| 伦精品一区二区三区| av卡一久久| av在线蜜桃| 成年免费大片在线观看| 国产成人91sexporn| 国产欧美日韩精品一区二区| 日本a在线网址| 中文亚洲av片在线观看爽| 51国产日韩欧美| 成人三级黄色视频| 男女那种视频在线观看| 欧美性猛交╳xxx乱大交人| 一个人看的www免费观看视频| 精品欧美国产一区二区三| 蜜臀久久99精品久久宅男| 一卡2卡三卡四卡精品乱码亚洲| 香蕉av资源在线| 精品久久久久久久人妻蜜臀av| 91在线精品国自产拍蜜月| 精品不卡国产一区二区三区| 少妇人妻精品综合一区二区 | 国产精品,欧美在线| 赤兔流量卡办理| 91久久精品国产一区二区成人| 亚洲中文字幕一区二区三区有码在线看| 热99re8久久精品国产| 色噜噜av男人的天堂激情| 欧美一级a爱片免费观看看| 18禁在线播放成人免费| 欧美区成人在线视频| 三级毛片av免费| 午夜a级毛片| 国产精品电影一区二区三区| 91在线精品国自产拍蜜月| 日本-黄色视频高清免费观看| 秋霞在线观看毛片| 丰满的人妻完整版| 国产精品人妻久久久久久| 99久久九九国产精品国产免费| 久久久色成人| 搡老岳熟女国产| 99热这里只有是精品50| 欧美一区二区亚洲| 欧美日韩综合久久久久久| 午夜激情欧美在线| 国产伦精品一区二区三区视频9| 青春草视频在线免费观看| 精品免费久久久久久久清纯| 亚洲中文字幕一区二区三区有码在线看| 夜夜夜夜夜久久久久| 国产伦精品一区二区三区视频9| 免费一级毛片在线播放高清视频| 亚洲国产精品成人综合色| 国产高清视频在线播放一区| 一级av片app| 在线免费观看的www视频| 在线国产一区二区在线| 国模一区二区三区四区视频| 99国产极品粉嫩在线观看| 国产av一区在线观看免费| 露出奶头的视频| 欧美最新免费一区二区三区| 非洲黑人性xxxx精品又粗又长| 久久婷婷人人爽人人干人人爱| 国产在线男女| 欧美成人一区二区免费高清观看| 变态另类成人亚洲欧美熟女| 国产精华一区二区三区| 久久久久久九九精品二区国产| 99热全是精品| 亚洲国产精品国产精品| 色在线成人网| 日日干狠狠操夜夜爽| 一个人看视频在线观看www免费| 一个人免费在线观看电影| 国产免费男女视频| 欧美激情在线99| 内射极品少妇av片p| 成年女人永久免费观看视频| 真实男女啪啪啪动态图| 国产男人的电影天堂91| 可以在线观看毛片的网站| 我的女老师完整版在线观看| 午夜久久久久精精品| 久久久久久伊人网av| 成人美女网站在线观看视频| 97在线视频观看| 免费av观看视频| 亚洲av第一区精品v没综合| 黄色一级大片看看| 欧美xxxx黑人xx丫x性爽| 亚洲第一区二区三区不卡| 最近2019中文字幕mv第一页| 国产三级在线视频| 老司机福利观看| 日日撸夜夜添| 国产探花在线观看一区二区| 欧美性猛交╳xxx乱大交人| 色综合色国产| 99久久精品国产国产毛片| 日韩av不卡免费在线播放| 欧美中文日本在线观看视频| 婷婷六月久久综合丁香| 精品久久久噜噜| 日韩强制内射视频| 国产毛片a区久久久久| 在线观看午夜福利视频| 在线天堂最新版资源| 国产精品一区www在线观看| 国产午夜精品久久久久久一区二区三区 | 男人的好看免费观看在线视频| 色哟哟哟哟哟哟| 欧美成人免费av一区二区三区| 天堂√8在线中文| 亚洲经典国产精华液单| 在线看三级毛片| 国产高清视频在线观看网站| 白带黄色成豆腐渣| 非洲黑人性xxxx精品又粗又长| 在线免费观看的www视频| 又黄又爽又免费观看的视频| 中文字幕精品亚洲无线码一区| 国产精品国产三级国产av玫瑰| 亚洲最大成人中文| 午夜爱爱视频在线播放| 久久久久久伊人网av| 一个人看视频在线观看www免费| 可以在线观看的亚洲视频| 国产精品一及| 看免费成人av毛片| 精品无人区乱码1区二区| 亚洲自偷自拍三级| 国产白丝娇喘喷水9色精品| 久久99热6这里只有精品| 久久精品影院6| 婷婷精品国产亚洲av| 性欧美人与动物交配| 久久久久久国产a免费观看| 亚洲中文字幕日韩|