2025/4/191控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院第3章

線性控制系統的數學模型薛定宇著《控制系統計算機輔助設計---MATLAB語言與應用》第二版，清華大學出版社2006CAI課件開發：張望舒哈爾濱工程大學薛定宇東北大學2025/4/192控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院系統的數學模型系統數學模型的重要性系統仿真分析必須已知數學模型系統設計必須已知數學模型本課程數學模型是基礎系統數學模型的獲取建模方法：從已知的物理規律出發，用數學推導的方式建立起系統的數學模型辨識方法：由實驗數據擬合系統的數學模型2025/4/193控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院系統數學模型的分類系統模型非線性線性連續離散混合單變量多變量定常時變2025/4/194控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院主要內容線性連續系統的數學模型與MATLAB表示線性離散時間系統的數學模型方框圖描述系統的化簡系統模型的相互轉換線性系統的模型降階（自學為主）線性系統的模型辨識（自學為主）飛機運動的動力學建模（補充內容）本章要點簡介2025/4/195控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.1連續線性系統的數學

模型與MATLAB表示3.1.1線性系統的狀態方程模型3.1.2線性系統的傳遞函數模型3.1.3線性系統的零極點模型3.1.4多變量系統的傳遞函數矩陣模型2025/4/196控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.1.1線性連續系統數學模型及MATLAB表示線性系統的傳遞函數模型為階次，為常數，物理可實現2025/4/197控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院傳遞函數的引入Pierre-SimonLaplace(1749--1827)，法國數學家

Laplace變換Laplace變換的一條重要性質：若則2025/4/198控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院傳遞函數表示數學方式MATLAB輸入語句2025/4/199控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院傳遞函數輸入舉例例3-1輸入傳遞函數模型MATLAB輸入語句在MATLAB環境中建立一個變量G2025/4/1910控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院另外一種傳遞函數輸入方法例3-2如何處理如下的傳遞函數？定義算子，再輸入傳遞函數2025/4/1911控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院應該根據給出傳遞函數形式選擇輸入方法例3-3輸入混合運算的傳遞函數模型顯然用第一種方法麻煩，所以2025/4/1912控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院MATLAB的傳遞函數對象2025/4/1913控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院傳遞函數屬性修改例3-4延遲傳遞函數，即若假設復域變量為，則2025/4/1914控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院傳遞函數參數提取由于使用單元數組，直接用不行有兩種方法可以提取參數這樣定義的優點：可以直接描述多變量系統第i輸出對第j輸入的傳遞函數2025/4/1915控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.1.2線性系統的狀態方程模型狀態方程模型狀態變量，階次n

，輸入和輸出非線性函數：一般非線性系統的狀態方程描述2025/4/1916控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院線性狀態方程時變模型線性時不變模型(lineartimeinvariant,LTI)2025/4/1917控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院線性時不變模型的MATLAB描述MATLAB輸入方法矩陣是方陣，為矩陣為矩陣，為矩陣可以直接處理多變量模型給出矩陣即可注意維數的兼容性2025/4/1918控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-52025/4/1919控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院帶時間延遲的狀態方程數學模型MATLAB輸入語句其他延遲屬性：ioDelay2025/4/1920控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.1.3線性系統的零極點模型零極點模型是因式型傳遞函數模型零點、極點和增益零極點模型的

MATLAB表示2025/4/1921控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-5零極點模型MATLAB輸入方法另一種輸入方法2025/4/1922控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.1.4多變量系統傳遞函數矩陣模型傳遞函數矩陣為第i輸出對第j輸入的傳遞函數可以先定義子傳遞函數，再由矩陣定義2025/4/1923控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-7多變量模型2025/4/1924控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.2線性離散時間系統的數學模型單變量系統：差分方程取代微分方程主要內容離散傳遞函數離散狀態方程2025/4/1925控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.2.1離散傳遞函數模型數學表示（Z變換代替Laplace變換）MATLAB表示（采樣周期）算子輸入方法：2025/4/1926控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-8離散傳遞函數，采樣周期MATLAB輸入方法另一種輸入方法2025/4/1927控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院離散延遲系統與輸入數學模型延遲為采樣周期的整數倍MATLAB輸入方法2025/4/1928控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院濾波器型描述方法濾波器型離散模型分子、分母除以記，則2025/4/1929控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院MATLAB表示方法例3-92025/4/1930控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.2.2離散狀態方程模型數學形式注意兼容性MATLAB表示方法2025/4/1931控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院離散延遲系統的狀態方程數學模型MATLAB表示方法2025/4/1932控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.3方框圖描述系統的化簡單環節模型前面已經介紹了實際系統為多個環節互連如何解決互連問題，獲得等效模型？主要內容控制系統的典型連接結構節點移動時的等效變換復雜系統模型的簡化2025/4/1933控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.3.1控制系統的典型連接結構系統串、并聯串聯傳遞函數并聯傳遞函數2025/4/1934控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院串、并聯狀態方程模型串聯系統的狀態方程并聯系統的狀態方程2025/4/1935控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院串、并聯系統的MATLAB求解若一個模型為傳遞函數、另一個為狀態方程，如何處理？將二者變換成同樣結構再計算基于MATLAB的計算方法串聯注意次序：多變量系統并聯優點，無需實現轉換2025/4/1936控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院系統的反饋連接反饋連接正反饋負反饋2025/4/1937控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院狀態方程的反饋等效方法其中若2025/4/1938控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院反饋連接的MATLAB求解LTI模型符號運算(置于@sym目錄)2025/4/1939控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-102025/4/1940例3-11控制器為對角矩陣2025/4/1941控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.3.2節點移動時的等效變換考慮模型難點：A點在回路間，移至輸出端2025/4/1942控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院節點移動2025/4/1943控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.3.3復雜系統模型的簡化例3-12原系統可以移動新支路模型2025/4/1944控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院得出2025/4/1945控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-13電機拖動模型

2025/4/1946控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院

信號單獨輸入得出另一個傳遞函數2025/4/1947控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院最終得出傳遞函數矩陣2025/4/1948控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.4系統模型的相互轉換前面介紹的各種模型之間的相互等效變換主要內容連續模型和離散模型的相互轉換系統傳遞函數的獲取控制系統的狀態方程實現狀態方程的最小實現傳遞函數與符號表達式的相互轉換2025/4/1949控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.4.1連續模型和離散模型的相互轉換連續狀態方程的解析階采樣周期選擇2025/4/1950控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院這樣可以得出離散模型記則可以得出離散狀態方程模型MATLAB函數直接求解2025/4/1951控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院還可以采用Tustin變換（雙線性變換）例3-14雙輸入模型，2025/4/1952控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院輸入模型、變換模型2025/4/1953控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-15時間延遲系統的離散化MATLAB求解零階保持器變換變換結果2025/4/1954控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院Tustin變換數學表示其他轉換方法FOH一階保持器matched單變量系統零極點不變imp脈沖響應不變準則2025/4/1955控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院離散模型連續化對前面的變換求逆Tustin反變換MATLAB求解(無需)2025/4/1956控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-16對前面的連續狀態方程模型離散化，對結果再連續化，則

可以基本上還原連續模型2025/4/1957控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.4.2系統傳遞函數的獲取已知狀態方程兩端Laplace變換則2025/4/1958控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院因此可以得出傳遞函數難點基于Fadeev-Fadeeva算法能得出更好結果由零極點模型，直接展開分子分母用MATLAB統一求解2025/4/1959控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-17多變量模型，求傳遞函數矩陣2025/4/1960控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.4.3控制系統的狀態方程實現由傳遞函數到狀態方程的轉換不同狀態變量選擇，結果不唯一默認變換方式，采用MATLAB函數G可以是傳遞函數、狀態方程和零極點模型適用于有延遲的、離散的或多變量模型2025/4/1961控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-18連續多變量模型狀態方程獲取2025/4/1962控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院得出的狀態方程模型ioDelay矩陣2025/4/1963控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院該模型可以轉換回傳遞函數矩陣得出的轉換結果2025/4/1964控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院均衡實現（banlancedrealization）由一般狀態方程輸入輸出關系顯著程度不明顯，需要進一步變換均衡實現是一種很有用的方式用MATLAB直接求解得出均衡實現的模型得出排序的Gram矩陣2025/4/1965控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例原系統模型引入(內部坐標變換)2025/4/1966控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.4.4狀態方程的最小實現例3-19觀察傳遞函數模型未見有何特殊求取零極點模型2025/4/1967控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院得出結果相同位置的零極點，可以對消問題：狀態方程如何處理？MATLAB解決方法2025/4/1968控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-20多變量模型不能直接看出是否最小實現2025/4/1969控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院MATLAB求解2025/4/1970控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.4.5傳遞函數與符號表達式

的相互轉換傳遞函數到符號表達式表達式到傳遞函數置于@sym目錄下2025/4/1971控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.5線性系統模型降階用低階模型近似高階模型和最小實現不同最早由EdwardJ.Davison提出(1966)主要內容與Routh算法時間延遲模型的近似帶有延遲的最優降階算法狀態空間的降階算法2025/4/1972控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.5.1降階算法

與Routh降階算法原始模型尋求降階模型假設2025/4/1973控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院展開原模型其中時間矩量可以遞推求出若已知狀態方程模型2025/4/1974控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院時間矩量的MATLAB求解降階思想：保留前時間矩量2025/4/1975控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院對比系數，則2025/4/1976控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院這樣可以得出2025/4/1977控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院

降階求解函數2025/4/1978控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-21原始模型Padé

近似結果2025/4/1979控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-22反例零極點模型求取穩定模型2025/4/1980控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院Padé

近似不穩定降階模型Padé

不能保證降階模型的穩定性不穩定降階模型可能得出穩定降階模型2025/4/1981控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院Routh降階方法與實例Routh算法（較煩瑣，從略）2025/4/1982控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院Routh算法的最大特色：穩定系統降階后能保證降階模型穩定性例3-23仍考慮穩定模型2025/4/1983控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.5.3時間延遲模型的Padé

近似純延遲的Padé近似方法近似函數純滯后逼近2025/4/1984控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院編寫MATLAB函數其中r/m任意選擇可以選擇0/m，以避免非最小相位模型2025/4/1985控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-24純延遲模型MATLAB求解擬合結果2025/4/1986控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-25已知帶有延遲的線性模型可以得出近似模型2025/4/1987控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.5.4帶有時間延遲系統的

次最優降階算法降階模型的降階效果誤差定義ISE準則2025/4/1988控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院原模型降階模型降階誤差定義2025/4/1989控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院參數向量誤差MATLAB實現（從略）調用格式2025/4/1990控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-26對給出的傳遞函數進行降階研究可以給出下面的語句得出的降階模型為2025/4/1991控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-27已知高階模型可以給出如下命令得出的降階模型2025/4/1992控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.5.4狀態方程模型的降階算法均衡實現模型的降階算法2025/4/1993控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院MATLAB求解函數例3-282025/4/1994控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院基于Schur均衡實現模型的降階算法MATLAB求解函數例3-29高階傳遞函數思路：先轉換成狀態方程，再降階2025/4/1995控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院模型輸入與降階Schur降階模型2025/4/1996控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院最優Hankel范數的降階模型近似MATLAB求解函數例3-30仍采用前面模型2025/4/1997控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院降階算法綜述狀態方程方法不能任意選擇分母分子階次，而很多傳遞函數方法可以降階效果比較，下章給出時域響應比較頻域響應比較降階模型的應用仿真應用（用途越來越小）控制器設計應用2025/4/1998控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.6線性系統的模型辨識模型辨識由已知實測數據獲得系統模型的方法實測數據時域響應數據、頻率響應數據主要內容離散系統辨識方法辨識信號生成多變量系統辨識離散系統在線辨識2025/4/1999控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.6.1離散系統的模型辨識離散傳遞函數模型對應的差分方程模型2025/4/19100控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院已知實測信號輸入輸出由數據可以得出2025/4/19101控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院矩陣形式定義殘差最小指標最小二乘解2025/4/19102控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院系統辨識工具箱求解T為結構體變量，T.a,T.b,tf(T)當然由前面的公式也能直接求解2025/4/19103控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-31實測數據2025/4/19104控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院基于MATLAB的求解2025/4/19105控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院數學形式辨識模型的提取還可以寫成2025/4/19106控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院還可以由下面語句求解辨識結果2025/4/19107控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院直接辨識方法辨識結果辨識界面：ident2025/4/19108控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.6.2離散系統辨識信號的生成問題：什么樣信號激勵系統，辨識效果最好？有豐富頻率信息的信號最好，如PRBS偽隨機二進制序列pseudo-randombinarysequence頻率豐富值為可重復構建MATLAB直接生成2025/4/19109控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院例3-32生成63個點的PRBS信號辨識效果殘差明顯減小2025/4/19110控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院連續系統的辨識可以考慮的方法連續頻率擬合方法，不惟一離散方法，再轉換成連續模型例3-332025/4/19111控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院正弦信號激勵辨識結果問題原因：輸入頻率單一2025/4/19112控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.6.3多變量離散系統的辨識離散傳遞函數矩陣模型其中例3-342025/4/19113控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院MATLAB求解2025/4/19114控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院得出的高階模型應該最小實現辨識結果2025/4/19115控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院3.6.4離散系統的遞推最小二乘辨識在系統運行過程中實時獲取系統參數，而不是象前面介紹的方法那樣一次性獲得模型，適合于變參數模型的實時控制廣泛應用于自適應控制漸近地逼近參數真值這里介紹算法，仿真研究將在后面介紹2025/4/19116控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院遞推最小二乘辨識傳遞函數模型差分方程模型待辨識參數2025/4/19117控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院遞推初值和加權矩陣輸入輸出數據向量遞推公式2025/4/19118控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院補充內容飛機的動力學模型坐標系地心慣性坐標系地心固定坐標系航跡坐標系機體坐標系2025/4/19119控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院坐標系之間的轉換地心慣性坐標系到地心固定坐標系以上公式使用mathematics軟件推導，matlab也具備這個功能2025/4/19120控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院坐標系之間的轉換地心固定坐標系到航跡坐標系（L維度，λ經度）2025/4/19121控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院坐標系之間的轉換航跡坐標系到機體坐標系（三個歐拉角）2025/4/19122控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程基本理論牛頓第二定理

動軸定理2025/4/19123控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程xyzP飛機的質心在地心慣性系中2025/4/19124控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程相對航跡坐標系的加速度隨地球自轉產生的向心加速度科氏加速度2025/4/19125控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程2025/4/19126控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程對于飛機2025/4/19127控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程忽略了地球自轉，球形地球2025/4/19128控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機的質心動力學方程質心動力學方程與導航方程是不耦合的質心動力學方程質心運動方程（導航方程）2025/4/19129控制系統計算機輔助設計---MATLAB語言與應用東北大學信息學院飛機繞質心動力學方程基本理論動量矩定理慣性系中很不方便體坐標系中注意：忽略了地球自轉202