(哈工大)系統(tǒng)辨識(shí)與自適應(yīng)控制——第一講

1、Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第1頁系統(tǒng)辨識(shí)與自適應(yīng)控制系統(tǒng)辨識(shí)與自適應(yīng)控制黃顯林、班曉軍黃顯林、班曉軍控制理論與制導(dǎo)技術(shù)研究中心控制理論與制導(dǎo)技術(shù)研究中心哈爾濱工業(yè)大學(xué)哈爾濱工業(yè)大學(xué)Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第2頁第一講第一講系統(tǒng)辨識(shí)的基本概念系統(tǒng)辨識(shí)的基本概念一、什么是系統(tǒng)辨識(shí)？一、什么是系統(tǒng)辨識(shí)？1. 機(jī)理分析建模方法 (白箱法)圖1 單級(jí)倒立擺實(shí)驗(yàn)裝置Harbin Institute of Technology

2、HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第3頁圖2 單級(jí)倒立擺示意圖 Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第4頁圖中所示變量名的物理含義如表圖中所示變量名的物理含義如表1所示。所示。 Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第5頁MruFPN步驟一：對(duì)小車進(jìn)行受力分析，小車的受力分析如圖步驟一：對(duì)小車進(jìn)行受力分析，小車的受力分析如圖3所所示。示。 圖3 小車受力分析圖圖中，圖中，P表示擺桿對(duì)小車水平方向上的作用力，單位表示擺桿

3、對(duì)小車水平方向上的作用力，單位N；N 表示擺桿對(duì)小車垂直方向上的作用力，單位（表示擺桿對(duì)小車垂直方向上的作用力，單位（N）。）。根據(jù)牛頓定律，小車水平方向上的力平衡方程為：根據(jù)牛頓定律，小車水平方向上的力平衡方程為：Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第6頁步驟二：對(duì)擺桿進(jìn)行受力分析，擺桿的受力如圖步驟二：對(duì)擺桿進(jìn)行受力分析，擺桿的受力如圖4所示。所示。圖圖4 擺桿受力分析圖擺桿受力分析圖擺桿水平方向上的力平衡方程如下，擺桿水平方向上的力平衡方程如下，22tdrdMPFutdrdFmgNPHarbin Instit

4、ute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第7頁sincos)sincos()cos()sin(2222 mlmlrmllrmlrtddmlrtddmP將式(1-3)合并可得下式,sincos)(sincossincos222 mlmlrrmMurmlmlrmrMumlmlrmrMFurmlumlrmM sincos)(2Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第8頁擺桿垂直方向上的力平衡方程式如下,)sincos()cos(222 mlltddmmgN)sincos(2

5、 mlmgN JwfPlNlcossinwfJmlrmlmgl 2cossinwfmglJmlrml sin)(cos2擺桿的轉(zhuǎn)矩平衡方程式如下,將3、7式代入8式并化簡(jiǎn)得化簡(jiǎn)得Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第9頁步驟三：由5式與10式連列即得到單級(jí)倒立擺動(dòng)力學(xué)非線性方程組。rmlumlrmMwfmglJmlrml sincos)(sin)(cos22步驟四：化成狀態(tài)空間描述。1222221222421122124431222221141112222221cos)()(sin)()(cossincoscos)

6、()(sin)(coscossincosxlmmlJmMumlJxmlxmlJxmlJxxglmxmlfxxxxxlmmlJmMfxmMxmglmMxmlxxmluxxxlmxxxHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第10頁2222422122443222421221)()()()()()(lmmlJmMumlJxmlJmlfxgxlmxxxlmmlJmMmlumlxfxmMmglxmMxxx (12)Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心

7、 第11頁問題：?jiǎn)栴}：(1). 效率低：隨著系統(tǒng)復(fù)雜程度的增加，建模過程愈加復(fù)效率低：隨著系統(tǒng)復(fù)雜程度的增加，建模過程愈加復(fù)雜；雜；(2). 不方便不方便“計(jì)算機(jī)計(jì)算機(jī)”在線決策。在線決策。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第12頁2. 系統(tǒng)辨識(shí)法系統(tǒng)辨識(shí)法(黑箱法黑箱法)能否根據(jù)能否根據(jù)“輸入、輸出數(shù)據(jù)輸入、輸出數(shù)據(jù)”獲取獲取“對(duì)象對(duì)象”的數(shù)學(xué)模型的數(shù)學(xué)模型呢？呢？例：原被控對(duì)象的差分形式為：例：原被控對(duì)象的差分形式為：Y(t) - 2.85y(t-1) + 2.717y(t-2) - 0.865y(t-3)

8、= u(t-1) + u(t-2) + u(t-3); 傳遞函數(shù)形式：傳遞函數(shù)形式：321321865. 0717. 285. 21)(zzzzzzzGHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第13頁給定輸入信號(hào)：給定輸入信號(hào)：10階階M序列。序列。輸入為為輸入為為0.002秒的秒的10階階M序列序列(周期周期 s )：10(21) 0.0022.04600.511.52-1-0.500.51 Time in Seconds Output圖5. 10階M序列Harbin Institute of Technology

9、HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第14頁00.050.10.150.2-1-0.500.51 Time in Seconds Output圖6. 10階M序列局部放大圖Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第15頁系統(tǒng)的響應(yīng)：00.511.52-500005000 Time in Seconds Output圖7. 系統(tǒng)對(duì)10階M序列的響應(yīng)曲線Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第16頁辨識(shí)過程：辨識(shí)過程：分別對(duì)應(yīng)的

10、辨識(shí)結(jié)果：分別對(duì)應(yīng)的辨識(shí)結(jié)果：給定階數(shù)給定階數(shù) 3，3，1根據(jù)不受噪聲干擾時(shí)的數(shù)據(jù)辨識(shí)出來的結(jié)果：根據(jù)不受噪聲干擾時(shí)的數(shù)據(jù)辨識(shí)出來的結(jié)果：Discrete-time IDPOLY model: A(q)y(t) = B(q)u(t) + e(t)A(q) = 1 - 2.85 q-1 + 2.717 q-2 - 0.865 q-3 B(q) = q-1 + q-2 + q-3 Estimated using ARX from data set mydata Loss function 6.25668e-024 and FPE 6.33214e-024 Sampling interval: 0.

11、002 321321865. 0717. 285. 21)(zzzzzzzGHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第17頁問題：?jiǎn)栴}：(1). 輸入信號(hào)為什么要選輸入信號(hào)為什么要選M序列，正余弦函數(shù)行不行，階序列，正余弦函數(shù)行不行，階躍信號(hào)行不行？躍信號(hào)行不行？(2). 預(yù)定的模型階次怎么確定？預(yù)定的模型階次怎么確定？(3). 具體的參數(shù)怎么確定？具體的參數(shù)怎么確定？Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第18頁3. 機(jī)理分析法機(jī)理分析法+

12、系統(tǒng)辨識(shí)法系統(tǒng)辨識(shí)法 （工程常用，灰箱法）（工程常用，灰箱法）電機(jī)系統(tǒng)：電機(jī)系統(tǒng)：Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第19頁圖8. 電機(jī)系統(tǒng)示意圖Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第20頁圖9. 電機(jī)系統(tǒng)的傳遞函數(shù) Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第21頁二、系統(tǒng)辨識(shí)方法的基本分類二、系統(tǒng)辨識(shí)方法的基本分類1. 參數(shù)辨識(shí)方法參數(shù)辨識(shí)方法a. 經(jīng)典辨識(shí)

13、方法經(jīng)典辨識(shí)方法 階躍響應(yīng)法；脈沖響應(yīng)法；頻域響應(yīng)法；相關(guān)分析法；譜分析法。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第22頁b. 最小二乘類參數(shù)辨識(shí)方法最小二乘類參數(shù)辨識(shí)方法 最小二乘一次性算法；最小二乘遞推算法增廣最小二乘算法；廣義最小二乘算法。c. 極大似然法和預(yù)報(bào)誤差方法極大似然法和預(yù)報(bào)誤差方法d. Bayes方法方法e. 模型參考自適應(yīng)方法模型參考自適應(yīng)方法Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第23頁2. 結(jié)構(gòu)辨識(shí)方法結(jié)構(gòu)辨識(shí)方法a

14、. 根據(jù)Hankel 矩陣的秩估計(jì)模型的階次； b. 利用行列式比估計(jì)模型的階次；c. 利用殘差的方差估計(jì)模型的階次；d. 利用Akaike準(zhǔn)則估計(jì)模型的階次；e. 利用最終預(yù)報(bào)誤差準(zhǔn)則估計(jì)模型的階次。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第24頁三、辨識(shí)的基本要素三、辨識(shí)的基本要素1. 輸入輸出數(shù)據(jù)（辨識(shí)的基礎(chǔ)）輸入輸出數(shù)據(jù)（辨識(shí)的基礎(chǔ)） 必須包含有關(guān)系統(tǒng)特性的足夠信息必須包含有關(guān)系統(tǒng)特性的足夠信息 時(shí)域的角度：信號(hào)變化劇烈，且呈現(xiàn)非周期性；時(shí)域的角度：信號(hào)變化劇烈，且呈現(xiàn)非周期性； 頻域的角度：頻譜寬。頻域的角度

15、：頻譜寬。2. 模型類模型類3. 等價(jià)準(zhǔn)則等價(jià)準(zhǔn)則評(píng)判評(píng)判“辨識(shí)得到的模型辨識(shí)得到的模型”是否滿足是否滿足“實(shí)際需要實(shí)際需要”的一的一個(gè)個(gè)“準(zhǔn)準(zhǔn) 則則”。 辨識(shí)就是按照一定的準(zhǔn)則從某一類模型中找出一個(gè)與辨識(shí)就是按照一定的準(zhǔn)則從某一類模型中找出一個(gè)與輸入輸出數(shù)據(jù)擬合得最好的模型。輸入輸出數(shù)據(jù)擬合得最好的模型。 Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第25頁例子：一個(gè)熱交換過程，如下圖所示。預(yù)建立例子：一個(gè)熱交換過程，如下圖所示。預(yù)建立T/Q模型。模型。經(jīng)觀測(cè)得到一組輸入輸出數(shù)據(jù)，記為經(jīng)觀測(cè)得到一組輸入輸出數(shù)據(jù)，記為Q(

16、k)，T(k)， 。 1,2,3,kLHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第26頁冷水熱水T 出口溫度 （輸出量）Q 蒸汽流量 （輸入量）圖10. 熱交換系統(tǒng)示意圖 Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第27頁選定模型類：選定模型類：選定等價(jià)準(zhǔn)則：選定等價(jià)準(zhǔn)則：最小化最小化 J，)()() 1()() 1()(11kenkQbkQbnkTakTakTnnLkkeJ12)(LknnnkQbkQbnkTakTakT1211)() 1()()

17、 1()(Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第28頁四、辨識(shí)問題的表達(dá)形式四、辨識(shí)問題的表達(dá)形式 (最小二乘格式最小二乘格式)圖11. 辨識(shí)問題的表達(dá)形式Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第29頁其中，其中，h(k): 系統(tǒng)的輸入變量；系統(tǒng)的輸入變量；z(k): 系統(tǒng)的量測(cè)輸出變量；系統(tǒng)的量測(cè)輸出變量；e(k): 系統(tǒng)的模型噪聲；系統(tǒng)的模型噪聲； ： 未知參數(shù)。未知參數(shù)。Harbin Institute of Technology

18、HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第30頁TNkhkhkhkh)(,),(),()(21TN,21)()()(kekhkzTHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第31頁例例 2. 對(duì)于給定質(zhì)量的氣體，不同的體積對(duì)于給定質(zhì)量的氣體，不同的體積V對(duì)應(yīng)著不同的對(duì)應(yīng)著不同的壓力壓力P。根據(jù)熱力學(xué)原理，壓力和體積之間存在如下的關(guān)。根據(jù)熱力學(xué)原理，壓力和體積之間存在如下的關(guān)系系PVrc其中，r 和 c 為待定常數(shù)；P和V 在各點(diǎn)上都是可觀測(cè)的。Harbin Institute of Technology H

19、IT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第32頁ckVrkPlg)(lg)(lgckVrkPlg)(lg)(lg)(lg)(kPkyTkVkh 1)(lg)(Tcrlg)()()(kekhkyT其中， Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第33頁例例2總結(jié)：最小二乘格式中的輸入、輸出量可以不是原總結(jié)：最小二乘格式中的輸入、輸出量可以不是原過程的輸入輸出量。過程的輸入輸出量。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第34頁圖12.

20、 辨識(shí)問題的表達(dá)形式Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第35頁五、遞推辨識(shí)算法的基本原理五、遞推辨識(shí)算法的基本原理 (在線辨識(shí)算法的基本原在線辨識(shí)算法的基本原理；遞推辨識(shí)算法的基本原理；類似神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過理；遞推辨識(shí)算法的基本原理；類似神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程程)Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第36頁圖13. 遞推辨識(shí)算法的基本原理Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)

21、技術(shù)研究中心 第37頁k 時(shí)刻的輸出值預(yù)測(cè)：時(shí)刻的輸出值預(yù)測(cè)：k 時(shí)刻的輸出誤差，或稱為時(shí)刻的輸出誤差，或稱為“新息新息”) 1()()( kkhkzT)( )()(kzkzkzHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第38頁原理：原理：將新息將新息(Innovation)“反饋反饋”到辨識(shí)算法中去，依據(jù)該值到辨識(shí)算法中去，依據(jù)該值修修正正“下一時(shí)刻下一時(shí)刻”模型參數(shù)的估計(jì)值。此迭代過程不斷進(jìn)模型參數(shù)的估計(jì)值。此迭代過程不斷進(jìn)行行下去，直至對(duì)應(yīng)的準(zhǔn)則函數(shù)取得最小值。下去，直至對(duì)應(yīng)的準(zhǔn)則函數(shù)取得最小值。反饋的又一功能。反

22、饋的又一功能。與神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)算法（與神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)算法（Bp）算法相似。）算法相似?！氨孀R(shí)辨識(shí)”的過程就是的過程就是“學(xué)習(xí)學(xué)習(xí)”的過程。的過程。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第39頁六、辨識(shí)的精度問題六、辨識(shí)的精度問題“時(shí)域評(píng)價(jià)結(jié)果時(shí)域評(píng)價(jià)結(jié)果”與與“頻域評(píng)價(jià)結(jié)果頻域評(píng)價(jià)結(jié)果”不一致。不一致。設(shè)對(duì)象具有如下傳遞函數(shù)：設(shè)對(duì)象具有如下傳遞函數(shù)：辨識(shí)得到的模型為：辨識(shí)得到的模型為： 11332201)(23ssssG16 .124 .321)(2sssGHarbin Institute of Technology

23、 HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第40頁以階躍響應(yīng)為評(píng)價(jià)指標(biāo)：精度較高。以階躍響應(yīng)為評(píng)價(jià)指標(biāo)：精度較高。)(/)()()(hththt 2 % Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第41頁010203040506070809010000.10.20.30.40.50.60.70.80.91Step ResponseTime (sec)Amplitude圖14. 階躍響應(yīng)的對(duì)比Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第

24、42頁0510152025303540455000.010.020.030.040.050.060.07Impulse ResponseTime (sec)Amplitude圖15. 脈沖響應(yīng)的對(duì)比Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第43頁-150-100-500Magnitude (dB)10-310-210-1100101102-270-180-900Phase (deg)Bode DiagramFrequency (rad/sec)圖16. 頻率特性的對(duì)比Harbin Institute of Techno

25、logy HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第44頁結(jié)論：結(jié)論：1. 辨識(shí)得到的模型只是實(shí)際過程的近似，需要有明確的辨識(shí)得到的模型只是實(shí)際過程的近似，需要有明確的評(píng)價(jià)指標(biāo)；評(píng)價(jià)指標(biāo)；2. 不同的評(píng)價(jià)指標(biāo)會(huì)得出不同的不同的評(píng)價(jià)指標(biāo)會(huì)得出不同的“精度評(píng)價(jià)精度評(píng)價(jià)”結(jié)果。結(jié)果。提示：不必要一味追求“精確”的模型。評(píng)價(jià)標(biāo)準(zhǔn)：實(shí)際應(yīng)用的效果。評(píng)價(jià)標(biāo)準(zhǔn)：實(shí)際應(yīng)用的效果。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第45頁八、辨識(shí)的應(yīng)用八、辨識(shí)的應(yīng)用1. 用于控制系統(tǒng)的設(shè)計(jì)與分析；（線性、用于控制系統(tǒng)的設(shè)計(jì)與分析；（

26、線性、T-S模糊模型）模糊模型）2. 用于在線控制；（非線性模型：模糊語言規(guī)則模型，用于在線控制；（非線性模型：模糊語言規(guī)則模型， 神經(jīng)網(wǎng)絡(luò)模型，模糊神經(jīng)網(wǎng)絡(luò)模型，神經(jīng)網(wǎng)絡(luò)模型，模糊神經(jīng)網(wǎng)絡(luò)模型，T-S模型）模型）3. 用于天氣、水文、人口、能源、客流量甚至股票走勢(shì)用于天氣、水文、人口、能源、客流量甚至股票走勢(shì) 等問題的預(yù)報(bào)；（非線性模型）等問題的預(yù)報(bào)；（非線性模型）4. 用于監(jiān)視過程參數(shù)并實(shí)現(xiàn)故障診斷（與故障診斷相結(jié)用于監(jiān)視過程參數(shù)并實(shí)現(xiàn)故障診斷（與故障診斷相結(jié) 合）。合）。Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第

27、46頁九、模型簡(jiǎn)介九、模型簡(jiǎn)介 (書目書目3 Page. 81)Equation Error Model Structure (ARX模型結(jié)構(gòu)模型結(jié)構(gòu))Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第47頁Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第48頁圖17. The ARX model structure.Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第49頁2. ARMAX

28、Model Structure (ARMAX 模型結(jié)構(gòu)模型結(jié)構(gòu))MA: 滑動(dòng)平均項(xiàng)(the moving average part) )()(keqCHarbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第50頁圖18. The ARMAX model structure.Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第51頁3. ARARX Model Structure (ARARX 模型結(jié)構(gòu)模型結(jié)構(gòu))()(1)()()()(keqDkuqBkyqAndndqdqdqD111)(Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第52頁圖19. The ARARX model structure.Harbin Institute of Technology HIT2010-02-20控制理論與制導(dǎo)技術(shù)研究中心 第53頁4. ARARMAX Model Structure (ARARMAX 模型結(jié)構(gòu))()()()()()()(keqDqCkuqBkyqAHarbin Institute