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

    結(jié)合前饋調(diào)參與迭代學(xué)習(xí)的數(shù)據(jù)驅(qū)動(dòng)控制方法

    2020-07-15 02:24:58薄雨蒙曹明生高慧斌
    控制理論與應(yīng)用 2020年6期
    關(guān)鍵詞:控制精度濾波器軌跡

    薄雨蒙,曹明生,高慧斌

    (1.中國(guó)科學(xué)院長(zhǎng)春光學(xué)精密機(jī)械與物理研究所,吉林長(zhǎng)春 130000;2.中國(guó)科學(xué)院大學(xué),北京 100049)

    1 引言

    在運(yùn)動(dòng)控制系統(tǒng)中,前饋控制得到了廣泛的應(yīng)用,特別是在對(duì)控制精度要求很高的軌跡跟蹤任務(wù)中,前饋控制可以取得良好的控制效果.前饋控制器基于先驗(yàn)知識(shí)獲得前饋控制信號(hào)對(duì)誤差進(jìn)行補(bǔ)償,因此可以有效提高控制精度.傳統(tǒng)的前饋控制方法包括迭代學(xué)習(xí)控制(iterative learning control,ILC),基于模型的前饋控制和迭代前饋調(diào)參(iterative feedforward tuning,IFT)等.

    迭代學(xué)習(xí)控制適合于具有重復(fù)性的控制任務(wù),它由Arimoto等人[1]于1984年提出,之后幾十年發(fā)展出了多種擴(kuò)展的方法[2-3],并可以與其他控制方法,如自適應(yīng)控制[4-6]、魯棒控制[7-9]和模糊控制[10-11]等進(jìn)行結(jié)合,根據(jù)系統(tǒng)特性和控制需求,可以達(dá)到更好的控制效果.

    迭代學(xué)習(xí)控制利用以前迭代周期獲得的信息,通過(guò)學(xué)習(xí)方法估計(jì)并補(bǔ)償系統(tǒng)中存在的不隨迭代次數(shù)變化的擾動(dòng).但在控制過(guò)程中,非重復(fù)性擾動(dòng)存在積累現(xiàn)象,從而影響控制誤差的收斂性[12].傳統(tǒng)的學(xué)習(xí)律學(xué)習(xí)的是系統(tǒng)的控制輸入,這要求期望軌跡在迭代域上嚴(yán)格重復(fù),一旦軌跡發(fā)生變化,迭代學(xué)習(xí)控制將會(huì)失效[13].針對(duì)期望軌跡變化的情況,文獻(xiàn)[14]提出了分段ILC,在文獻(xiàn)[15]中以分段ILC為基礎(chǔ),將任務(wù)分解為多個(gè)子任務(wù)進(jìn)行單獨(dú)學(xué)習(xí),但其控制性能受分段方式影響,并且增加了算法實(shí)現(xiàn)的復(fù)雜度.

    相反,基于模型的前饋控制也可以達(dá)到較高的控制精度,但期望軌跡的變化不會(huì)引起控制效果的惡化[16-17].基于模型的前饋控制使用系統(tǒng)逆模型或系統(tǒng)敏感度函數(shù)建立學(xué)習(xí)濾波器,利用系統(tǒng)信息進(jìn)行前饋補(bǔ)償,獲得的控制效果依賴于系統(tǒng)模型的質(zhì)量[18].但在工程應(yīng)用中,對(duì)復(fù)雜系統(tǒng)建立精確模型難度大,成本高,更傾向于使用數(shù)據(jù)驅(qū)動(dòng)控制方法,即直接使用輸入輸出數(shù)據(jù)進(jìn)行控制器設(shè)計(jì)的方法[19-20].

    采用基函數(shù)可以將迭代學(xué)習(xí)控制與基于模型的前饋控制的優(yōu)點(diǎn)相結(jié)合,獲得適用于變軌跡任務(wù),高控制精度且不依賴于系統(tǒng)模型的前饋控制方法[21].通過(guò)基函數(shù)的引入,信號(hào)被投影到基函數(shù)張成的空間中,組成的前饋控制器可用于估計(jì)系統(tǒng)的逆模型[22].

    使用基函數(shù)將前饋控制器參數(shù)化,可以將系統(tǒng)模型辨識(shí)問(wèn)題轉(zhuǎn)換為參數(shù)化前饋控制器的迭代前饋調(diào)參問(wèn)題,并通過(guò)最優(yōu)化方法得到最優(yōu)前饋控制器[23-24].文獻(xiàn)[25]使用數(shù)據(jù)驅(qū)動(dòng)的梯度下降法進(jìn)行調(diào)參,但文獻(xiàn)[25-26]證明此方法在存在噪聲時(shí)將得到參數(shù)的有偏估計(jì),因此引入輔助變量來(lái)獲得參數(shù)的無(wú)偏估計(jì),文獻(xiàn)[27]對(duì)參數(shù)估計(jì)的精確度進(jìn)行了分析.文獻(xiàn)[28]提出了高階的優(yōu)化方法,使用過(guò)去所有迭代周期獲得的誤差信號(hào)進(jìn)行參數(shù)估計(jì),算法具有高精度和高噪聲容忍度.文獻(xiàn)[29-30]使用有理基函數(shù)以獲得更高的精度.

    將ILC與IFT進(jìn)行比較,發(fā)現(xiàn)ILC對(duì)于重復(fù)性任務(wù)可以得到很高的控制精度,但期望軌跡的變化會(huì)引起極大的效果惡化[30];反之,IFT則適用于變期望軌跡任務(wù),不依賴于系統(tǒng)模型,但控制精度低于ILC.

    為了在保證高控制精度的同時(shí),改善ILC在期望軌跡變化時(shí)的效果惡化問(wèn)題,本文針對(duì)執(zhí)行變軌跡步進(jìn)掃描運(yùn)動(dòng)的系統(tǒng),提出了一種將IFT與ILC進(jìn)行結(jié)合的數(shù)據(jù)驅(qū)動(dòng)的控制結(jié)構(gòu).通過(guò)引入IFT,得到最優(yōu)系統(tǒng)前饋控制器,消除期望軌跡引入的擾動(dòng),提高軌跡變化時(shí)的控制精度;通過(guò)ILC,消除系統(tǒng)重復(fù)性擾動(dòng),進(jìn)一步提高控制精度.在IFT的計(jì)算過(guò)程中,引入輸入整形濾波器,在保證梯度下降法計(jì)算精度的同時(shí)簡(jiǎn)化了計(jì)算過(guò)程.本文提出的算法在軌跡變化時(shí)仍能保證高控制精度,同時(shí)是一種數(shù)據(jù)驅(qū)動(dòng)控制方法,可以直接應(yīng)用于不同工業(yè)系統(tǒng)而不依賴于系統(tǒng)的先驗(yàn)信息.

    2 問(wèn)題提出

    2.1 控制結(jié)構(gòu)

    控制系統(tǒng)框圖如圖1所示,其中:C是反饋控制器,Kff是參數(shù)化前饋控制器,Ky是輸入整形濾波器,KILC是迭代學(xué)習(xí)控制器,KIFT是系統(tǒng)前饋控制器.實(shí)際系統(tǒng)的系統(tǒng)模型P可以看作離散單輸入單輸出線性時(shí)不變系統(tǒng).系統(tǒng)傳遞函數(shù)可以寫作

    其中B0(z),A0(z)∈R(F).此外,已知參考信號(hào)r是對(duì)預(yù)先規(guī)劃好的n階連續(xù)可導(dǎo)的軌跡以采樣時(shí)間Ts采樣得到的N個(gè)點(diǎn)的序列,ry是經(jīng)過(guò)輸入整形濾波器得到的整形參考信號(hào),y是系統(tǒng)輸出,ey是整形參考信號(hào)ry與輸出y之間的誤差,uff為Kff輸出的前饋控制信號(hào),ufb為反饋控制器C輸出的控制信號(hào),uILC為迭代學(xué)習(xí)控制器離線輸出的前饋控制信號(hào),uIFT為KIFT輸出的系統(tǒng)前饋控制信號(hào).所有信號(hào)都為以采樣時(shí)間Ts進(jìn)行采樣的n維離散序列.

    圖1 控制結(jié)構(gòu)示意圖Fig.1 Control configuration

    本文提出的控制方法在參數(shù)估計(jì)階段和迭代學(xué)習(xí)控制階段使用不同的控制結(jié)構(gòu).在系統(tǒng)控制框圖中,對(duì)虛線框表示的迭代調(diào)參前饋控制器Kff和輸入整形濾波器Ky通過(guò)基函數(shù)進(jìn)行參數(shù)化,并通過(guò)最優(yōu)化迭代調(diào)參算法迭代估計(jì)收斂的最優(yōu)參數(shù)值.當(dāng)?shù){(diào)參過(guò)程結(jié)束后,移除控制框圖中的虛線部分,并將最優(yōu)參數(shù)值輸入到系統(tǒng)前饋控制器KIFT中,此時(shí)的控制結(jié)構(gòu)是由前饋控制器和迭代學(xué)習(xí)控制器組合成的迭代學(xué)習(xí)控制結(jié)構(gòu).

    2.2 參數(shù)化的前饋控制器和輸入整形濾波器

    如圖2所示,前饋控制信號(hào)uff可以由期望軌跡r和參數(shù)化前饋控制器Kff(θf(wàn)f)表示為

    同時(shí),整形參考信號(hào)ry可以由r和參數(shù)化輸入整形濾波器Ky(θy)表示為

    其中Kff(θf(wàn)f)和Ky(θy)由定義1給出.

    圖2 加入輸入整形濾波器的迭代前饋調(diào)參框圖Fig.2 Configuration of IFT with input shaping filter

    定義1通過(guò)引入基函數(shù),參數(shù)化前饋控制器Kff(θf(wàn)f)和輸入整形濾波器Ky(θy)可表示為

    其中:

    多項(xiàng)式基函數(shù)表示為

    各項(xiàng)基函數(shù)為

    其中nff和ny為構(gòu)成Kff和Ky的基函數(shù)的個(gè)數(shù).通過(guò)基函數(shù),可以將輸入信號(hào)分解為各階導(dǎo)數(shù).因?yàn)樗械摩譱(z)都是關(guān)于z的函數(shù),為了簡(jiǎn)化表示,在后文中省略z.此外,為了簡(jiǎn)化,還將Ψy和Ψff組合為一個(gè)向量Ψ,表示為

    類似地,將參數(shù)向量θy和θf(wàn)f組合為一個(gè)向量θ,同上式表示為

    通過(guò)引入輸入整形濾波器Ky,輸入整形后的參考信號(hào)ry可以表示為

    根據(jù)式(2),可得到如下條件:

    由式(3)可推得,ry與r在平穩(wěn)段相等,而在移動(dòng)段有一段延遲,延遲長(zhǎng)度等于Ky中基函數(shù)的個(gè)數(shù).因此,在平穩(wěn)段,e=r-y=ry-y=ey.若只需要在平穩(wěn)段的跟蹤誤差最小,控制性能的評(píng)判標(biāo)準(zhǔn)可表示為

    最優(yōu)化的θf(wàn)f和θy需要滿足如下定義:

    定義2前饋控制器和輸入整形濾波器參數(shù)向量是最優(yōu)的,如果在控制系統(tǒng)中引入使用相應(yīng)參數(shù)向量建立的前饋控制器和輸入整形濾波器測(cè)得的誤差信號(hào)滿足其中S=(I +PC)-1是系統(tǒng)敏感度函數(shù).

    建立參數(shù)化前饋控制器Kff(θf(wàn)f)和輸入整形濾波器Ky(θy),通過(guò)最優(yōu)化迭代調(diào)參算法,算法收斂時(shí)可以得到滿足定義2的最優(yōu)參數(shù)向量θ*和相應(yīng)的的和此時(shí)Vy(θk)達(dá)到最小值.詳細(xì)算法將在下節(jié)中給出.

    3 算法設(shè)計(jì)

    3.1 最優(yōu)化迭代前饋調(diào)參算法

    根據(jù)第2.2節(jié)建立的參數(shù)化前饋控制器和輸入整形濾波器,在圖2中,第k次迭代測(cè)得的信號(hào)和uk可以使用第k次迭代的參數(shù)向量θk擴(kuò)展為參數(shù)化形式,表示如下:

    式中 K(θk)=Kff(θk)+CKy(θk).

    根據(jù)式(6)-(7),通過(guò)每次迭代的uk和yk,可以得到S和(SP)的估計(jì)值Sest和(SP)est:

    假設(shè)1K-1(θk)是穩(wěn)定的.

    注1如果K-1(θk)不穩(wěn)定,可以通過(guò)使用文獻(xiàn)[25]中的穩(wěn)定求逆算法計(jì)算得到K-1(θk).

    將式(9)-(12)代入式(8)中,得到梯度向量的估計(jì)值:

    其中ρ為學(xué)習(xí)增益.設(shè)置學(xué)習(xí)律為

    根據(jù)以上計(jì)算,通過(guò)梯度下降法迭代求解最優(yōu)化問(wèn)題(13)的步驟可以概括如下:

    程序1迭代估計(jì)θk.

    初始化:

    1) 確定合適的學(xué)習(xí)增益ρ,0 <ρ <1,使算法在收斂速度和噪聲容忍度中取得平衡;

    2) 選擇合適的基函數(shù)Ψ,用于建立參數(shù)化前饋控制器Kff,輸入整形濾波器Ky和參考信號(hào)r;

    3) 確定合適的參數(shù)向量初始值θ0;

    主程序:

    1) 使用參數(shù)向量θk運(yùn)行控制程序,獲得uk和yk;

    3) 計(jì)算學(xué)習(xí)矩陣

    4) 計(jì)算下一次迭代的θk+1=θk+Lkek;

    5) k →k+1.返回步驟1.

    注2通過(guò)程序1,根據(jù)梯度下降法的收斂性,可以得到一個(gè)收斂的參數(shù)向量θ,且這個(gè)參數(shù)向量滿足定義1,收斂性證明詳見文獻(xiàn)[24,26].將這個(gè)收斂的參數(shù)向量記作θ*.

    3.2 通過(guò)迭代前饋調(diào)參得到前饋控制器

    注意到,誤差ey是整形參考信號(hào)ry與輸出信號(hào)y之間的差值,并不是跟蹤期望軌跡r得到的跟蹤誤差.因此,最小化的Vy(θk)并不意味著對(duì)整個(gè)期望軌跡的V(θk)的最小化.

    由式(5),可以得出,誤差ey等于0的條件為

    由式(14),可以得到如下結(jié)論:建立Kff和Ky是為了得到系統(tǒng)逆模型的估計(jì)值,使ey為0.如果建立的Kff和Ky使ey為0,此時(shí)就可以作為系統(tǒng)逆模型的估計(jì)值.這與文獻(xiàn)[22]中的迭代學(xué)習(xí)系統(tǒng)辨識(shí)概念類似.由梯度下降法的最優(yōu)性,當(dāng)使用梯度下降法的程序1收斂時(shí),可認(rèn)為結(jié)果滿足定義2,即=0.此時(shí)由收斂時(shí)的組合成的可以作為系統(tǒng)逆模型的估計(jì)值.因此在控制結(jié)構(gòu)中引入KIFT=作為前饋控制器.

    在Kff和Ky收斂為后,移除原本的前饋控制器Kff和輸入整形濾波器Ky,將KIFT=作為前饋控制器加入到系統(tǒng)中,補(bǔ)償期望軌跡引入的擾動(dòng).KIFT輸出的前饋控制量uIFT為

    注3如果不穩(wěn)定,可以通過(guò)零相位跟蹤技術(shù)[31]得到穩(wěn)定的KIFT.

    使用組合的KIFT=而不是直接使用作為前饋控制器,是因?yàn)閱为?dú)的Kff只能估計(jì)系統(tǒng)傳遞函數(shù)的分母部分,無(wú)法準(zhǔn)確估計(jì)系統(tǒng)逆模型[24].而不直接將Kff參數(shù)化為分子和分母都為多項(xiàng)式的形式,即文獻(xiàn)[29]中的有理基函數(shù):

    是因?yàn)閷?duì)有理基函數(shù)參數(shù)化的Krat進(jìn)行最優(yōu)化迭代運(yùn)算時(shí),算法復(fù)雜度高,且需要選擇合適的使算法收斂的權(quán)重矩陣,而這個(gè)矩陣是N×N維[29]或N×n維的[30],維度高導(dǎo)致待定參數(shù)過(guò)多,難以選擇使算法收斂的參數(shù)矩陣.而將Krat分解為本文的Kff和Ky,不需要設(shè)置高維權(quán)重參數(shù)矩陣進(jìn)行,顯著降低了算法復(fù)雜度,同時(shí)保證算法效果.

    3.3 迭代前饋調(diào)參與迭代學(xué)習(xí)控制結(jié)合算法

    對(duì)比ILC與IFT,發(fā)現(xiàn)IFT魯棒性很強(qiáng),期望軌跡變化時(shí)控制效果依然穩(wěn)定,而ILC在變軌跡時(shí)存在嚴(yán)重的效果惡化現(xiàn)象;然而,ILC經(jīng)過(guò)幾次學(xué)習(xí)就可以達(dá)到比IFT更高的控制精度,這是由于ILC還去除了KIFT無(wú)法消除的其他重復(fù)性擾動(dòng).

    僅使用第3.2節(jié)中獲得的最優(yōu)前饋控制器KIFT得到的控制精度較低.為了提高控制精度,引入迭代學(xué)習(xí)控制器KILC,采集上一個(gè)迭代周期的誤差ek-1,離線更新前饋控制量消除其他重復(fù)性擾動(dòng).在本文中,選用PD型迭代學(xué)習(xí)控制器,第k次迭代時(shí)KILC輸出的前饋控制量為

    其中kp和kd為學(xué)習(xí)律的學(xué)習(xí)參數(shù).

    根據(jù)上文所述的最優(yōu)化梯度下降迭代前饋調(diào)參方法,得到前饋控制器KIFT,并與PD型迭代學(xué)習(xí)控制結(jié)合,可以將完整算法流程圖表示如圖3.算法的具體步驟表示如程序2.

    圖3 本文提出算法的流程圖Fig.3 Flow chart of proposed algorithm

    程序2

    1) 運(yùn)行程序1,獲得收斂的最優(yōu)參數(shù)向量θ*;

    2) 移除Kff和Ky;

    3) 建立KIFT=

    4) 在控制結(jié)構(gòu)中加入KIFT和KILC;

    5) 根據(jù)新的控制結(jié)構(gòu)進(jìn)行實(shí)驗(yàn),直到滿足停止條件.

    注4由于算法只是預(yù)先消除了由期望軌跡引入的誤差,對(duì)迭代過(guò)程沒(méi)有影響,因此穩(wěn)定性與PD型迭代學(xué)習(xí)律的穩(wěn)定性相同.證明詳見文獻(xiàn)[14].

    4 迭代前饋調(diào)參仿真

    在本段中,對(duì)結(jié)合前饋調(diào)參與迭代學(xué)習(xí)的算法進(jìn)行仿真,驗(yàn)證算法的有效性.

    仿真使用的離散傳遞函數(shù)模型如下:

    其中采樣時(shí)間為Ts=1×10-4s.仿真中采用Kp=1.05,Kd=25的PD型反饋控制器,可以寫作

    在系統(tǒng)中加入擾動(dòng)μ,它的模型為μ=H?,其中:

    ?是標(biāo)準(zhǔn)差為1×10-7的零均值白噪聲.

    仿真使用的傳遞函數(shù)的頻率響應(yīng)如圖4(a)中實(shí)線所示.根據(jù)圖中的頻率響應(yīng),可以看出系統(tǒng)主要?jiǎng)討B(tài)特性集中在0.1~100 Hz的頻段,而Kff和Ky的引入是為了補(bǔ)償系統(tǒng)的主要?jiǎng)討B(tài)特性.系統(tǒng)使用5階點(diǎn)到點(diǎn)的期望軌跡r.

    根據(jù)對(duì)仿真模型的分析和第2.2 節(jié)中定義1,建立Kff(θf(wàn)f)和Ky(θy)如下:

    對(duì)于參數(shù)向量θ,設(shè)置初始值為

    設(shè)置學(xué)習(xí)增益ρ為0.9,使用以上的Kff和Ky,運(yùn)行程序1.對(duì)運(yùn)行結(jié)果作圖:圖4(a)為對(duì)第2次迭代的Kff和Ky分別作出頻率響應(yīng),與傳遞函數(shù)P進(jìn)行對(duì)比的圖形;圖4(b)為θk隨迭代次數(shù)的變化;圖4(c)為第2次迭代的與P的對(duì)比圖.分析運(yùn)行結(jié)果,可得出以下結(jié)論:

    1) 通過(guò)圖4(a)看出,算法收斂時(shí),單獨(dú)的Kff和Ky都不足以描述系統(tǒng)的動(dòng)態(tài)特性;

    2) 由圖4(b)可以看出,程序在第2次運(yùn)行后θk變化很小,說(shuō)明調(diào)參算法收斂速度快,且可以將第3次迭代的結(jié)果看作算法收斂時(shí)得到的最優(yōu)值.此外,還可以通過(guò)調(diào)整學(xué)習(xí)增益ρ來(lái)調(diào)整算法的收斂速度;

    圖4(a) 仿真中算法收斂時(shí)Ky,Kff與P的伯德圖Fig.4(a) Simulation bode diagram of converged Ky and Kffcompared with P

    圖4(b) θ隨迭代次數(shù)變化的仿真結(jié)果Fig.4(b) Simulation result of θ during iteration

    圖4(c) 仿真中算法收斂時(shí)Ky(Kff)-1與P的伯德圖Fig.4(c) Simulation bode diagram of converged Ky(Kff)-1 compared with P

    5 實(shí)驗(yàn)

    在本章中,筆者將提出的控制方法在實(shí)驗(yàn)平臺(tái)上進(jìn)行測(cè)試,驗(yàn)證算法的有效性.

    5.1 實(shí)驗(yàn)平臺(tái)

    本文使用如圖5(a)所示的直線伺服控制平臺(tái).實(shí)驗(yàn)伺服系統(tǒng)由上位機(jī)PC機(jī),POWER PMAC控制卡,AMAC驅(qū)動(dòng)器和NEWPORT IDL225直線電機(jī)運(yùn)動(dòng)平臺(tái)組成.直線電機(jī)的位移由Heidenhain增量式編碼器測(cè)量,分辨率為20 μm,然后通過(guò)插值得到10 nm的分辨率.電機(jī)由數(shù)字伺服驅(qū)動(dòng)器驅(qū)動(dòng).算法使用C語(yǔ)言在PowerPmac控制器中進(jìn)行實(shí)現(xiàn).系統(tǒng)的采樣周期為Ts=0.0001 s.

    使用掃頻法獲得實(shí)驗(yàn)平臺(tái)的頻率響應(yīng)如圖5(b)所示.為了驗(yàn)證算法,規(guī)劃了如圖5(c)的兩種5階點(diǎn)到點(diǎn)期望軌跡r1和r2,每個(gè)期望軌跡由N=2000個(gè)點(diǎn)組成.實(shí)驗(yàn)中,使用線性PD型反饋控制器,表示為Cpd=其中:Kp=1.05,Kd=25.

    圖5(a) 直線伺服控制實(shí)驗(yàn)平臺(tái)Fig.5(a) Photograph of experimental platform

    圖5(b) 實(shí)驗(yàn)平臺(tái)的頻率響應(yīng)圖Fig.5(b) Frequency response of experimental platform

    圖5(c) 實(shí)驗(yàn)使用的兩種期望軌跡和速度、加速度示意圖Fig.5(c) Two different reference trajectorys with their velocity and acceleration

    5.2 輸入整形迭代前饋調(diào)參

    如前文所述,最優(yōu)化的Kff(θ*)和Ky(θ*)需要估計(jì)系統(tǒng)主要頻段的動(dòng)態(tài)特性.根據(jù)圖5(b)中的頻率響應(yīng)圖,系統(tǒng)主要頻段為0.1~100 Hz.基函數(shù)選擇為

    使用圖5(c)中的r1作為期望軌跡,在實(shí)驗(yàn)平臺(tái)上運(yùn)行程序1,測(cè)量每次迭代的θ和誤差并記錄.為了說(shuō)明收斂情況,對(duì)建立Kff的4個(gè)參數(shù)作圖如圖6(a).可以看出,θ在第2次迭代后波動(dòng)很小,可以認(rèn)為參數(shù)向量在第2次迭代時(shí)收斂,算法收斂速度快,因此對(duì)程序2的運(yùn)行時(shí)間影響不大.

    為了驗(yàn)證IFT的控制精度,記錄每次迭代的Vk,并作圖如圖6(b).可以發(fā)現(xiàn),第2次迭代之后,V 的波動(dòng)很小,這也與仿真相符.通過(guò)IFT,軌跡跟蹤的誤差平方和V 減小了兩個(gè)數(shù)量級(jí),說(shuō)明IFT可以對(duì)期望軌跡引入的擾動(dòng)進(jìn)行一定的補(bǔ)償,提高控制精度.算法收斂時(shí)得到的平穩(wěn)段誤差信號(hào)如圖6(c)所示,其中1ct=10 nm.

    圖6(a) 參數(shù)向量隨迭代次數(shù)的變化情況Fig.6(a) Result of θ during iteration

    圖6(b) 迭代前饋調(diào)參得到的Vk隨迭代次數(shù)k的變化情況Fig.6(b) Vk achieved from IFT

    圖6(c) 迭代前饋調(diào)參得到的平穩(wěn)段誤差信號(hào)Fig.6(c) Error signal during dwell period achieved from IFT

    5.3 本文提出的IFT-ILC算法

    為了得到最優(yōu)前饋控制器KIFT,需要計(jì)算KIFT=并希望它描述系統(tǒng)主要頻段的動(dòng)態(tài)特性.由第5.2節(jié)的實(shí)驗(yàn)結(jié)果,使用第2次迭代得到的參數(shù)向量作為最優(yōu)參數(shù)向量θ*,建立KIFT=Kff(θ*)·(Ky(θ*))-1.將系統(tǒng)頻率響應(yīng)與KIFT的幅頻特性作比較,并同時(shí)作出Kff和Ky的幅頻特性,如圖7(a).圖中:灰實(shí)線為系統(tǒng)幅頻特性,黑實(shí)線為KIFT的幅頻特性,虛線為Ky,點(diǎn)線為Kff.從圖中可以看出,KIFT可以在主要頻段描述系統(tǒng)動(dòng)態(tài)特性,而單獨(dú)的Kff(θ*)和Ky(θ*)都不足以描述系統(tǒng).因此選用KIFT作為系統(tǒng)前饋控制器是合適的.

    圖7(b)顯示了僅加入最優(yōu)前饋控制器的算法與本文提出算法的V 隨迭代次數(shù)變化情況的對(duì)比,期望軌跡為r1.由圖中結(jié)果可以看出,僅使用KIFT作為前饋控制器的算法在第1 次迭代后,Vk由不加前饋時(shí)的2.5×10-8m2降低到了7.5×10-10m2,之后的10次迭代中Vk不再有顯著的下降;而使用本文提出的IFT與ILC相結(jié)合的算法,在第1次迭代后效果與僅使用KIFT時(shí)相同,但之后的10次迭代中Vk仍隨著迭代逐漸下降,在第10次迭代后Vk降低到3.4×10-11m2.可見期望軌跡不變時(shí),本文提出的算法精度遠(yuǎn)高于無(wú)學(xué)習(xí)的精度.

    為了驗(yàn)證算法在軌跡變化時(shí)的效果,在第5次迭代后將期望軌跡由r1變?yōu)閞2,對(duì)ILC和本文提出的IFTILC算法進(jìn)行實(shí)驗(yàn),記錄Vk的變化情況如圖7(c).期望軌跡變化時(shí),PD型迭代學(xué)習(xí)控制的Vk由1.8×10-10m2變?yōu)?.2×10-9m2,發(fā)生了嚴(yán)重的效果惡化,而使用本文提出的算法,在期望軌跡變化時(shí),Vk為1.2×10-9m2,效果惡化得到減輕.可見,在期望軌跡發(fā)生變化時(shí),僅使用PD型學(xué)習(xí)律的ILC的控制效果惡化十分嚴(yán)重,而本文提出的算法能減輕效果惡化,提高控制的魯棒性.

    圖7(a) 實(shí)驗(yàn)平臺(tái)頻率響應(yīng)與,和KIFT的伯德圖Fig.7(a) Frequency response of experimental platform,,and KIFT

    圖7(b) 迭代前饋調(diào)參與本文算法的Vk的對(duì)比圖Fig.7(b) Comparison of Vk of IFT with proposed algorithm

    圖7(c) 期望軌跡發(fā)生變化時(shí)迭代學(xué)習(xí)控制與本文算法的Vk的對(duì)比圖Fig.7(c) Comparison of Vk of ILC with proposed algorithm when reference signal changes

    6 結(jié)論

    為了改善變軌跡控制任務(wù)中ILC的效果惡化問(wèn)題,同時(shí)保證高控制精度,本文提出了一種結(jié)合輸入整形濾波器,迭代前饋調(diào)參與迭代學(xué)習(xí)控制的數(shù)據(jù)驅(qū)動(dòng)控制算法.算法是數(shù)據(jù)驅(qū)動(dòng)的,不需要已知模型信息,而是使用輸入輸出數(shù)據(jù),通過(guò)迭代前饋調(diào)參在算法運(yùn)行過(guò)程中對(duì)模型進(jìn)行估計(jì).在迭代前饋調(diào)參過(guò)程中,使用基函數(shù)對(duì)控制器進(jìn)行參數(shù)化,并使用梯度下降法進(jìn)行參數(shù)優(yōu)化,從而獲得最優(yōu)前饋控制器;輸入整形濾波器的引入在保證模型估計(jì)精度的同時(shí)簡(jiǎn)化了計(jì)算.最優(yōu)前饋控制器對(duì)期望軌跡引入誤差進(jìn)行補(bǔ)償,減輕了軌跡變化帶來(lái)的效果惡化,而迭代學(xué)習(xí)控制的引入保證了控制精度.為了驗(yàn)證算法的有效性,本文將提出的算法應(yīng)用到直線電機(jī)系統(tǒng)中,并與現(xiàn)有算法進(jìn)行了對(duì)比,實(shí)驗(yàn)結(jié)果表明本算法能夠更好的處理變軌跡的軌跡跟蹤問(wèn)題.

    猜你喜歡
    控制精度濾波器軌跡
    基于無(wú)擾濾波器和AED-ADT的無(wú)擾切換控制
    基于多源異構(gòu)信息融合的采摘機(jī)械臂驅(qū)動(dòng)控制研究
    軌跡
    軌跡
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    開關(guān)電源EMI濾波器的應(yīng)用方法探討
    電子制作(2018年16期)2018-09-26 03:26:50
    MW級(jí)太空發(fā)電站微波能量波束指向控制精度分析
    軌跡
    基于安卓的智能車轉(zhuǎn)速系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)
    進(jìn)化的軌跡(一)——進(jìn)化,無(wú)盡的適應(yīng)
    成人影院久久| 亚洲精品成人av观看孕妇| 老熟妇仑乱视频hdxx| 蜜桃在线观看..| 男女之事视频高清在线观看| 少妇人妻久久综合中文| 啦啦啦中文免费视频观看日本| 水蜜桃什么品种好| av欧美777| 成在线人永久免费视频| 美女主播在线视频| 亚洲av电影在线观看一区二区三区| 热99国产精品久久久久久7| 99久久人妻综合| 91精品三级在线观看| 国产亚洲精品久久久久5区| 人人妻人人澡人人看| 91老司机精品| 日韩中文字幕欧美一区二区| 久久av网站| 欧美激情高清一区二区三区| 精品久久蜜臀av无| www日本在线高清视频| 精品人妻1区二区| 国产在线免费精品| 在线观看免费午夜福利视频| 免费一级毛片在线播放高清视频 | 男女午夜视频在线观看| 在线观看一区二区三区激情| 在线永久观看黄色视频| 国产精品久久久久久人妻精品电影 | 十八禁人妻一区二区| 国产激情久久老熟女| 不卡av一区二区三区| av在线app专区| 亚洲精品国产精品久久久不卡| 久久性视频一级片| 精品国产乱子伦一区二区三区 | 在线观看免费午夜福利视频| 国产成人精品久久二区二区免费| 女人高潮潮喷娇喘18禁视频| 精品国产乱码久久久久久男人| 久久这里只有精品19| 久久青草综合色| 亚洲成人免费电影在线观看| 十八禁网站网址无遮挡| 波多野结衣一区麻豆| 久久国产精品人妻蜜桃| 每晚都被弄得嗷嗷叫到高潮| 美女高潮到喷水免费观看| 成年人黄色毛片网站| 青春草亚洲视频在线观看| 欧美精品亚洲一区二区| 久久亚洲国产成人精品v| 最近最新免费中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 国产老妇伦熟女老妇高清| 亚洲欧美色中文字幕在线| 叶爱在线成人免费视频播放| 可以免费在线观看a视频的电影网站| 咕卡用的链子| 91成人精品电影| 亚洲人成77777在线视频| 操出白浆在线播放| 王馨瑶露胸无遮挡在线观看| 免费少妇av软件| 少妇裸体淫交视频免费看高清 | 啦啦啦啦在线视频资源| 精品人妻熟女毛片av久久网站| 亚洲美女黄色视频免费看| 亚洲av电影在线观看一区二区三区| 一边摸一边做爽爽视频免费| 亚洲成人国产一区在线观看| 一级片免费观看大全| 亚洲精品国产一区二区精华液| 老司机影院成人| 欧美97在线视频| 久久99一区二区三区| 精品欧美一区二区三区在线| 欧美日韩精品网址| 欧美老熟妇乱子伦牲交| 成年人黄色毛片网站| xxxhd国产人妻xxx| 久久久久久久国产电影| 国产精品自产拍在线观看55亚洲 | 99re6热这里在线精品视频| 久久精品亚洲av国产电影网| 久久久久久久精品精品| 久久久久久人人人人人| 免费观看人在逋| 日韩中文字幕视频在线看片| 欧美精品亚洲一区二区| 久久久水蜜桃国产精品网| tocl精华| 精品国产国语对白av| 极品少妇高潮喷水抽搐| 国产福利在线免费观看视频| 男女国产视频网站| 黑人巨大精品欧美一区二区mp4| 久久久欧美国产精品| 欧美日韩成人在线一区二区| 一级毛片精品| 一个人免费看片子| 久久精品成人免费网站| 日本av手机在线免费观看| 飞空精品影院首页| 老司机影院毛片| av片东京热男人的天堂| 一级a爱视频在线免费观看| 国产免费av片在线观看野外av| 亚洲人成电影观看| 免费久久久久久久精品成人欧美视频| 亚洲成人手机| 97在线人人人人妻| 亚洲一码二码三码区别大吗| 久久影院123| 啦啦啦啦在线视频资源| 久久中文字幕一级| 一级毛片电影观看| 国产精品一二三区在线看| 日本黄色日本黄色录像| 亚洲成国产人片在线观看| 在线十欧美十亚洲十日本专区| 人妻 亚洲 视频| 亚洲精品日韩在线中文字幕| 免费一级毛片在线播放高清视频 | 成年美女黄网站色视频大全免费| 久久久久精品国产欧美久久久 | 亚洲国产欧美一区二区综合| 精品国产乱子伦一区二区三区 | 日本猛色少妇xxxxx猛交久久| 亚洲五月色婷婷综合| 久久av网站| 国产国语露脸激情在线看| 黑人操中国人逼视频| 亚洲国产毛片av蜜桃av| 亚洲人成电影观看| 精品久久久精品久久久| 午夜福利一区二区在线看| 青草久久国产| 天天躁日日躁夜夜躁夜夜| 国产精品久久久av美女十八| 中文欧美无线码| 国产av国产精品国产| 涩涩av久久男人的天堂| 电影成人av| 高潮久久久久久久久久久不卡| 妹子高潮喷水视频| 欧美激情高清一区二区三区| 人人澡人人妻人| 亚洲av成人不卡在线观看播放网 | 亚洲欧美日韩另类电影网站| 日韩视频一区二区在线观看| 国产片内射在线| 久久久久久人人人人人| cao死你这个sao货| 欧美另类亚洲清纯唯美| 国产精品久久久久成人av| 丝袜喷水一区| 亚洲精品久久久久久婷婷小说| 高潮久久久久久久久久久不卡| 97人妻天天添夜夜摸| 在线观看一区二区三区激情| 欧美精品一区二区免费开放| 自拍欧美九色日韩亚洲蝌蚪91| 久久影院123| 侵犯人妻中文字幕一二三四区| 青春草视频在线免费观看| 欧美亚洲 丝袜 人妻 在线| 国产99久久九九免费精品| 99久久人妻综合| 人成视频在线观看免费观看| 多毛熟女@视频| 国产亚洲一区二区精品| 免费高清在线观看日韩| 亚洲欧洲精品一区二区精品久久久| 亚洲精品国产区一区二| 日韩视频一区二区在线观看| 两个人看的免费小视频| 午夜免费成人在线视频| 青草久久国产| 日韩免费高清中文字幕av| av欧美777| 午夜精品久久久久久毛片777| 丰满少妇做爰视频| 久久亚洲国产成人精品v| 波多野结衣一区麻豆| 夜夜骑夜夜射夜夜干| 考比视频在线观看| av片东京热男人的天堂| 亚洲精品中文字幕在线视频| 大型av网站在线播放| 久久影院123| 国产在线一区二区三区精| 高清黄色对白视频在线免费看| 久久午夜综合久久蜜桃| 天天影视国产精品| 一进一出抽搐动态| 亚洲欧美激情在线| 人人妻人人澡人人爽人人夜夜| 熟女少妇亚洲综合色aaa.| 亚洲精品国产区一区二| 亚洲国产欧美一区二区综合| 搡老熟女国产l中国老女人| av在线老鸭窝| 欧美日韩成人在线一区二区| 久久久久网色| 午夜福利在线免费观看网站| 美女扒开内裤让男人捅视频| av天堂在线播放| a级毛片在线看网站| 国产1区2区3区精品| 欧美av亚洲av综合av国产av| 欧美人与性动交α欧美精品济南到| 一本大道久久a久久精品| 久久av网站| 一本色道久久久久久精品综合| 侵犯人妻中文字幕一二三四区| 久久久国产精品麻豆| 欧美日韩亚洲国产一区二区在线观看 | 国产精品偷伦视频观看了| 99久久综合免费| 国产淫语在线视频| 99久久国产精品久久久| 亚洲精品在线美女| cao死你这个sao货| 午夜福利影视在线免费观看| 欧美激情极品国产一区二区三区| 欧美精品av麻豆av| 国产在线观看jvid| 久久久精品国产亚洲av高清涩受| 建设人人有责人人尽责人人享有的| 最近中文字幕2019免费版| 国产高清国产精品国产三级| 亚洲成人免费av在线播放| 国产色视频综合| 人妻人人澡人人爽人人| 国产91精品成人一区二区三区 | 欧美乱码精品一区二区三区| 1024香蕉在线观看| 色综合欧美亚洲国产小说| 亚洲精品国产区一区二| 视频在线观看一区二区三区| 我要看黄色一级片免费的| 一二三四社区在线视频社区8| 成年女人毛片免费观看观看9 | 精品第一国产精品| 成年动漫av网址| 黄色片一级片一级黄色片| 手机成人av网站| 90打野战视频偷拍视频| 亚洲 欧美一区二区三区| 精品视频人人做人人爽| 十分钟在线观看高清视频www| 性色av一级| 肉色欧美久久久久久久蜜桃| 亚洲成人手机| 亚洲精品粉嫩美女一区| 天堂8中文在线网| 久久狼人影院| 午夜精品久久久久久毛片777| 黄色视频在线播放观看不卡| 日韩欧美免费精品| 啦啦啦在线免费观看视频4| 男女无遮挡免费网站观看| 女警被强在线播放| 欧美亚洲日本最大视频资源| 日韩大码丰满熟妇| 欧美日韩国产mv在线观看视频| 久久天堂一区二区三区四区| 久久久久视频综合| 搡老熟女国产l中国老女人| 97精品久久久久久久久久精品| 涩涩av久久男人的天堂| 男女下面插进去视频免费观看| 高清av免费在线| 一边摸一边抽搐一进一出视频| 男人添女人高潮全过程视频| 母亲3免费完整高清在线观看| 熟女少妇亚洲综合色aaa.| 一级毛片女人18水好多| 99国产精品一区二区三区| 18禁观看日本| 国产亚洲精品久久久久5区| 中文精品一卡2卡3卡4更新| 国产精品熟女久久久久浪| 午夜福利乱码中文字幕| 欧美亚洲 丝袜 人妻 在线| 看免费av毛片| 日本91视频免费播放| 国产免费福利视频在线观看| 久久久精品94久久精品| 波多野结衣av一区二区av| 日韩有码中文字幕| 久久久久久人人人人人| 波多野结衣一区麻豆| 亚洲精品成人av观看孕妇| 国产成人精品无人区| 亚洲av片天天在线观看| 午夜福利免费观看在线| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 搡老乐熟女国产| 日韩视频在线欧美| 狂野欧美激情性xxxx| 亚洲av成人不卡在线观看播放网 | 免费黄频网站在线观看国产| 一进一出抽搐动态| 欧美+亚洲+日韩+国产| 国产一区二区三区av在线| 一级黄色大片毛片| 久久免费观看电影| 日韩 亚洲 欧美在线| 最近最新中文字幕大全免费视频| 免费日韩欧美在线观看| 91av网站免费观看| 免费在线观看视频国产中文字幕亚洲 | www.熟女人妻精品国产| 午夜福利一区二区在线看| 老司机福利观看| 人人妻,人人澡人人爽秒播| 99香蕉大伊视频| 蜜桃国产av成人99| 波多野结衣一区麻豆| 韩国高清视频一区二区三区| 两性夫妻黄色片| 麻豆乱淫一区二区| 最新的欧美精品一区二区| 欧美精品亚洲一区二区| 99国产精品一区二区蜜桃av | 亚洲 欧美一区二区三区| 亚洲男人天堂网一区| 中文字幕最新亚洲高清| 精品乱码久久久久久99久播| 91老司机精品| 亚洲精品美女久久av网站| 人人澡人人妻人| 国产不卡av网站在线观看| 国产成人啪精品午夜网站| 天天影视国产精品| 2018国产大陆天天弄谢| 精品人妻在线不人妻| 熟女少妇亚洲综合色aaa.| 免费少妇av软件| 可以免费在线观看a视频的电影网站| 极品少妇高潮喷水抽搐| 99久久精品国产亚洲精品| 日韩欧美国产一区二区入口| 精品国产一区二区久久| 91字幕亚洲| 国产极品粉嫩免费观看在线| 欧美+亚洲+日韩+国产| 免费在线观看黄色视频的| 亚洲精品中文字幕在线视频| 99精品久久久久人妻精品| 悠悠久久av| 69精品国产乱码久久久| 欧美 日韩 精品 国产| 亚洲天堂av无毛| 国产精品久久久久久精品古装| 嫁个100分男人电影在线观看| 秋霞在线观看毛片| 亚洲美女黄色视频免费看| 性色av一级| 久久精品亚洲av国产电影网| 免费在线观看影片大全网站| 黄色片一级片一级黄色片| 亚洲 国产 在线| 亚洲一区二区三区欧美精品| 水蜜桃什么品种好| 免费高清在线观看视频在线观看| 在线观看www视频免费| 久久久久国内视频| 欧美日韩黄片免| 波多野结衣一区麻豆| 老司机靠b影院| 日本av免费视频播放| 不卡一级毛片| 18禁黄网站禁片午夜丰满| 久久久久久人人人人人| 精品久久久久久电影网| 亚洲成人免费av在线播放| 亚洲成国产人片在线观看| 国产精品麻豆人妻色哟哟久久| 国产97色在线日韩免费| 久久狼人影院| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 久久青草综合色| 国产在线免费精品| 天天躁日日躁夜夜躁夜夜| 91av网站免费观看| 9热在线视频观看99| 欧美在线一区亚洲| 欧美成人午夜精品| 69精品国产乱码久久久| av天堂在线播放| 成人国产av品久久久| 18禁国产床啪视频网站| 啦啦啦中文免费视频观看日本| 国产一区二区激情短视频 | 国产福利在线免费观看视频| 久久久久国产精品人妻一区二区| 十八禁人妻一区二区| 日韩中文字幕欧美一区二区| 狠狠狠狠99中文字幕| 国产免费视频播放在线视频| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产色婷婷电影| a级片在线免费高清观看视频| 91大片在线观看| 国产伦理片在线播放av一区| 久久国产亚洲av麻豆专区| 免费高清在线观看日韩| 老司机亚洲免费影院| 永久免费av网站大全| 欧美日韩av久久| 色精品久久人妻99蜜桃| 妹子高潮喷水视频| 久久这里只有精品19| 国产野战对白在线观看| 91麻豆av在线| 1024香蕉在线观看| 啦啦啦视频在线资源免费观看| 性高湖久久久久久久久免费观看| 我要看黄色一级片免费的| 纵有疾风起免费观看全集完整版| 亚洲专区国产一区二区| 欧美人与性动交α欧美精品济南到| 老熟妇仑乱视频hdxx| 老司机午夜十八禁免费视频| 美女中出高潮动态图| 国产高清videossex| 亚洲男人天堂网一区| 一级毛片精品| 免费观看av网站的网址| 下体分泌物呈黄色| 国产精品偷伦视频观看了| 纯流量卡能插随身wifi吗| 老司机午夜十八禁免费视频| 搡老妇女老女人老熟妇| 日本免费a在线| 操出白浆在线播放| 欧美国产日韩亚洲一区| 日韩 欧美 亚洲 中文字幕| 性欧美人与动物交配| 精品人妻1区二区| 色综合亚洲欧美另类图片| 欧美日韩一级在线毛片| 国产精品久久久久久久电影 | 亚洲狠狠婷婷综合久久图片| 观看免费一级毛片| 亚洲国产欧洲综合997久久,| 国产1区2区3区精品| 国产精品九九99| 精品高清国产在线一区| 久久久国产精品麻豆| 亚洲午夜理论影院| 国产成人精品久久二区二区91| 国产成+人综合+亚洲专区| 最近视频中文字幕2019在线8| 韩国av一区二区三区四区| 日日干狠狠操夜夜爽| av福利片在线观看| 最近在线观看免费完整版| 男人舔奶头视频| www国产在线视频色| 国产不卡一卡二| 日本三级黄在线观看| 午夜免费观看网址| 国产1区2区3区精品| 亚洲九九香蕉| 国产精品98久久久久久宅男小说| 一本一本综合久久| av中文乱码字幕在线| 精品日产1卡2卡| 久久天堂一区二区三区四区| 日日爽夜夜爽网站| 黄频高清免费视频| 欧美成人午夜精品| 99久久国产精品久久久| 国产精品香港三级国产av潘金莲| 欧美大码av| 国产午夜精品久久久久久| 国产熟女xx| 国产精品久久久久久人妻精品电影| 亚洲熟女毛片儿| 国产成人aa在线观看| 亚洲欧美日韩无卡精品| 正在播放国产对白刺激| 久久久久久国产a免费观看| 少妇熟女aⅴ在线视频| 日韩高清综合在线| 久久精品人妻少妇| 欧美成人免费av一区二区三区| 1024香蕉在线观看| 欧美日本视频| 久久天躁狠狠躁夜夜2o2o| 麻豆久久精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 国产成人系列免费观看| 中出人妻视频一区二区| 精品乱码久久久久久99久播| 日韩欧美在线乱码| 夜夜爽天天搞| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人精品二区| 丰满的人妻完整版| 真人做人爱边吃奶动态| 最近最新中文字幕大全免费视频| 欧美三级亚洲精品| 亚洲成a人片在线一区二区| 国产高清激情床上av| 啦啦啦观看免费观看视频高清| 俄罗斯特黄特色一大片| 一级黄色大片毛片| 狂野欧美白嫩少妇大欣赏| 中文字幕人成人乱码亚洲影| 国产成人av激情在线播放| 成在线人永久免费视频| 三级男女做爰猛烈吃奶摸视频| 亚洲人成电影免费在线| www日本黄色视频网| 国产在线精品亚洲第一网站| 极品教师在线免费播放| 小说图片视频综合网站| 美女 人体艺术 gogo| 日韩精品免费视频一区二区三区| 中文字幕久久专区| 51午夜福利影视在线观看| 成人特级黄色片久久久久久久| 99热6这里只有精品| 国产主播在线观看一区二区| 成人欧美大片| 久久久久久久精品吃奶| 免费在线观看影片大全网站| 午夜福利成人在线免费观看| 亚洲男人天堂网一区| 中文资源天堂在线| 国模一区二区三区四区视频 | 国产精品九九99| 亚洲午夜精品一区,二区,三区| 精品欧美国产一区二区三| 亚洲自拍偷在线| 91大片在线观看| 香蕉丝袜av| 亚洲人成网站在线播放欧美日韩| 久久这里只有精品19| 久久这里只有精品中国| 久久 成人 亚洲| cao死你这个sao货| 成人国产一区最新在线观看| 欧洲精品卡2卡3卡4卡5卡区| av有码第一页| 国语自产精品视频在线第100页| 中国美女看黄片| 午夜福利高清视频| 精品久久久久久成人av| 亚洲国产欧洲综合997久久,| 久久久久国产精品人妻aⅴ院| 国产av又大| 午夜福利欧美成人| 夜夜躁狠狠躁天天躁| 性色av乱码一区二区三区2| 亚洲精品在线美女| e午夜精品久久久久久久| 国产精品一及| 香蕉久久夜色| 亚洲五月婷婷丁香| 精华霜和精华液先用哪个| 中文字幕久久专区| 久久亚洲真实| www.自偷自拍.com| 伦理电影免费视频| 亚洲性夜色夜夜综合| 久久香蕉精品热| 国产区一区二久久| 国产成年人精品一区二区| 亚洲自偷自拍图片 自拍| 最近在线观看免费完整版| 亚洲精品久久成人aⅴ小说| 午夜日韩欧美国产| 这个男人来自地球电影免费观看| 无遮挡黄片免费观看| 91老司机精品| а√天堂www在线а√下载| 黄色 视频免费看| 欧美成人免费av一区二区三区| 亚洲中文av在线| 变态另类成人亚洲欧美熟女| 久久久久性生活片| 美女黄网站色视频| 亚洲国产精品合色在线| 可以在线观看的亚洲视频| 香蕉av资源在线| 欧美 亚洲 国产 日韩一| 好男人在线观看高清免费视频| www日本黄色视频网| 看黄色毛片网站| 日韩大码丰满熟妇| 久久久国产欧美日韩av| 高潮久久久久久久久久久不卡| 99久久无色码亚洲精品果冻| 成人国产一区最新在线观看| xxx96com| 成人午夜高清在线视频| 18禁观看日本| 久久久久精品国产欧美久久久| 听说在线观看完整版免费高清| 日韩中文字幕欧美一区二区| 老司机午夜福利在线观看视频| 日韩三级视频一区二区三区| 欧美精品亚洲一区二区| 美女 人体艺术 gogo| 五月玫瑰六月丁香|