1、 廣西大學實驗報告紙 姓名： 成績：學院：電氣工程學院專業(yè)：班級：實驗內容：Lyapunov方程求解2015年 月 日【實驗時間】2015年6月22日【實驗目的】掌握求解Lyapunov方程的一種方法，了解并使用MATLAB中相應函數(shù)?！緦嶒炘O備與軟件】硬件：PC機一臺；軟件：MATLAB/Simulink。【實驗原理】線性定常系統(tǒng)漸進穩(wěn)定的Lyapunov方程判據(jù)線性定常連續(xù)系統(tǒng)為漸進穩(wěn)定的充要條件是：對給定的任一個正定對稱陣Q，都存在唯一的對稱正定陣P，滿足如下矩陣Lyapunov方程：該條件在傳遞函數(shù)最小實現(xiàn)下等價于：全部特征根都是負實數(shù)或實部為負的復數(shù)，亦即全部根都位于左半復平面。線性

2、定常離散系統(tǒng)為漸進穩(wěn)定的充要條件：對給定的任一個正定對稱陣Q，都存在唯一的對稱正定陣P，滿足如下矩陣Lyapunov方程：該條件在傳遞函數(shù)最小實現(xiàn)下等價于：全部特征根的摸均小于1，即都在單位圓內。在MATLAB控制工具箱中，函數(shù)lyap和dlyap用來求解lyapunov方程。P=lyap（,Q)可解連續(xù)時間系統(tǒng)的lyapunov方程，其中，Q和A為具有相同維數(shù)的方陣（A是系統(tǒng)矩陣）。如果Q是對稱的，則解P也是對稱的。P=dlyap（，Q）可解離散時間系統(tǒng)的lyapunov方程，其中，Q和G為具有相同維數(shù)的方陣(G是系統(tǒng)矩陣)。如果Q是對稱的，則解P也是對稱的。3、連續(xù)情況下的最小相位系統(tǒng)：系

3、統(tǒng)的零點均在左半復平面，但系統(tǒng)首先是穩(wěn)定的，其他情況為非最小相位系統(tǒng)?！緦嶒瀮热荨⒎椒?、過程與分析】題目1實驗內容：輸入連續(xù)狀態(tài)空間模型:選取正定矩陣,求穩(wěn)定性判別矩陣P，判定系統(tǒng)是否穩(wěn)定。求線性系統(tǒng)階躍響應曲線，并判定是否為最小相位系統(tǒng)，求系統(tǒng)的實現(xiàn)，判定是否是最小實現(xiàn)并比較。題目1實驗過程及結果分析：根據(jù)題意，在實驗中，先通過運算可以得出結果，根據(jù)結果做出如下的c文件程序： = 1 * GB3 * MERGEFORMAT 、由實驗c文件程序運行后結果：得到正定矩陣P： = 2 * GB3 * MERGEFORMAT 、由題意得出系統(tǒng)的響應曲線：由圖可知：該系統(tǒng)是漸進穩(wěn)定的。求特征根由結果可

4、以得出，此系統(tǒng)特征根的實部全部都為負數(shù)，亦全部的根都在左邊平面。所以該系統(tǒng)為最小相位系統(tǒng)。所以，根據(jù)題意，更改A矩陣，求其階躍響應曲線，并進行比較得：之前的A矩陣：更改之前的特征值：更改前的階躍響應：更改之后的A矩陣：更改之后的特征值：更改后的階躍響應：對比特征值可知，更改矩陣A后特征根有一個為正數(shù)，即在右半平面；對比階躍響應圖可知，更改矩陣A后，其階躍響應為發(fā)散的。題目2實驗內容：輸入離散狀態(tài)空間模型選定正定矩陣,求穩(wěn)定性判別矩陣P。請定義離散情況下的最小相位系統(tǒng)。求線性系統(tǒng)階躍響應曲線，并按你所定義的判別矩陣是否為最小相位系統(tǒng)。題目2實驗過程及結果分析：根據(jù)題意，在實驗中，先通過運算可以得

5、出結果，根據(jù)結果做出如下的c文件程序： = 1 * GB3 * MERGEFORMAT 、由實驗c文件程序運行后結果：得到矩陣P：由圖可知：得到的矩陣P是非正定的。 = 2 * GB3 * MERGEFORMAT 、由題意得出系統(tǒng)的響應曲線： 由圖可知：階躍響應發(fā)散，該系統(tǒng)是不穩(wěn)定的。系統(tǒng)不穩(wěn)定，所以不是最小相位系統(tǒng)。之前的G矩陣：更改之前的特征值：更改前的階躍響應：更改之后的G矩陣：更改之后的特征值：更改后的階躍響應：對比更改前后特征值可知，更改矩陣G后特征根全部為負數(shù)，即在右半平面。對比階躍響應圖可知，更改矩陣G后，其階躍響應為收斂的?！緦嶒灴偨Y】通過本次實驗了解并掌握了Lyapunov方

6、程的一種用MATLAB求解的方法，并熟悉了線性定常系統(tǒng)漸進穩(wěn)定的Lyapunov方程判據(jù)和求解lyapunov方程的一些函數(shù)。附錄資料：MATLAB的30個方法1 內部常數(shù)pi 圓周率 exp(1)自然對數(shù)的底數(shù)ei 或j 虛數(shù)單位Inf或 inf 無窮大 2 數(shù)學運算符a+b 加法a-b減法a*b矩陣乘法a.*b數(shù)組乘法a/b矩陣右除ab矩陣左除a./b數(shù)組右除a.b數(shù)組左除ab 矩陣乘方a.b數(shù)組乘方-a負號 共軛轉置.一般轉置3 關系運算符=等于大于=大于或等于=不等于4 常用內部數(shù)學函數(shù) 指數(shù)函數(shù)exp(x)以e為底數(shù)對數(shù)函數(shù)log(x)自然對數(shù)，即以e為底數(shù)的對數(shù)log10(x)常用

7、對數(shù)，即以10為底數(shù)的對數(shù)log2(x)以2為底數(shù)的x的對數(shù)開方函數(shù)sqrt(x)表示x的算術平方根絕對值函數(shù)abs(x)表示實數(shù)的絕對值以及復數(shù)的模三角函數(shù)（自變量的單位為弧度）sin(x)正弦函數(shù)cos(x)余弦函數(shù)tan(x)正切函數(shù)cot(x)余切函數(shù)sec(x)正割函數(shù)csc(x)余割函數(shù)反三角函數(shù) asin(x)反正弦函數(shù)acos(x)反余弦函數(shù)atan(x)反正切函數(shù)acot(x)反余切函數(shù)asec(x)反正割函數(shù)acsc(x)反余割函數(shù)雙曲函數(shù) sinh(x)雙曲正弦函數(shù)cosh(x)雙曲余弦函數(shù)tanh(x)雙曲正切函數(shù)coth(x)雙曲余切函數(shù)sech(x)雙曲正割函數(shù)cs

8、ch(x)雙曲余割函數(shù)反雙曲函數(shù) asinh(x)反雙曲正弦函數(shù)acosh(x)反雙曲余弦函數(shù)atanh(x)反雙曲正切函數(shù)acoth(x)反雙曲余切函數(shù)asech(x)反雙曲正割函數(shù)acsch(x)反雙曲余割函數(shù)求角度函數(shù)atan2(y,x)以坐標原點為頂點，x軸正半軸為始邊，從原點到點（x，y）的射線為終邊的角，其單位為弧度，范圍為（ ， 數(shù)論函數(shù)gcd(a,b)兩個整數(shù)的最大公約數(shù)lcm(a,b)兩個整數(shù)的最小公倍數(shù)排列組合函數(shù)factorial(n)階乘函數(shù)，表示n的階乘 復數(shù)函數(shù) real(z)實部函數(shù)imag(z)虛部函數(shù)abs(z)求復數(shù)z的模angle(z)求復數(shù)z的輻角，其范

9、圍是（ ， conj(z)求復數(shù)z的共軛復數(shù)求整函數(shù)與截尾函數(shù)ceil(x)表示大于或等于實數(shù)x的最小整數(shù)floor(x)表示小于或等于實數(shù)x的最大整數(shù)round(x)最接近x的整數(shù)最大、最小函數(shù)max(a，b，c，)求最大數(shù)min(a，b，c，)求最小數(shù)符號函數(shù) sign(x)5 自定義函數(shù)-調用時：“返回值列=M文件名（參數(shù)列）”function 返回變量=函數(shù)名（輸入變量） 注釋說明語句段（此部分可有可無）函數(shù)體語句 6進行函數(shù)的復合運算compose(f,g) 返回值為f(g(y)compose(f,g,z) 返回值為f(g(z)compose(f,g,x,.z) 返回值為f(g(z)

10、compose(f,g,x,y,z) 返回值為f(g(z)7 因式分解syms 表達式中包含的變量 factor(表達式) 8 代數(shù)式展開syms 表達式中包含的變量 expand(表達式)9 合并同類項syms 表達式中包含的變量 collect(表達式，指定的變量)10 進行數(shù)學式化簡syms 表達式中包含的變量 simplify(表達式)11 進行變量替換syms 表達式和代換式中包含的所有變量 subs(表達式，要替換的變量或式子，代換式)12 進行數(shù)學式的轉換調用Maple中數(shù)學式的轉換命令，調用格式如下：maple(Maple的數(shù)學式轉換命令) 即：maple(convert(表達

11、式，form)將表達式轉換成form的表示方式 maple(convert(表達式，form, x) 指定變量為x，將依賴于變量x的函數(shù)轉換成form的表示方式（此指令僅對form為exp與sincos的轉換式有用） 13 解方程solve(方程，變元) 注：方程的等號用普通的等號： = 14 解不等式調用maple中解不等式的命令即可，調用形式如下： maple(maple中解不等式的命令)具體說，包括以下五種：maple( solve（不等式）) maple( solve（不等式，變元） ) maple( solve（不等式，變元） ) maple( solve（不等式，變元） ) map

12、le( solve（不等式，變元） )15 解不等式組調用maple中解不等式組的命令即可，調用形式如下： maple(maple中解不等式組的命令) 即：maple( solve（不等式組，變元組） )16 畫圖方法：先產生橫坐標的取值和相應的縱坐標的取值，然后執(zhí)行命令： plot(x,y) 方法2：fplot(f(x)，xmin,xmax) fplot(f(x)，xmin,xmax,ymin,ymax) 方法3：ezplot(f(x) ezplot(f(x) ,xmin,xmax) ezplot(f(x) ,xmin,xmax,ymin,ymax) 17 求極限（1）極限：syms x l

13、imit(f(x), x, a) （2）單側極限：左極限：syms x limit(f(x), x, a,left)右極限：syms x limit(f(x), x, a,right) 18 求導數(shù)diff(f(x) diff(f(x),x) 或者：syms x diff(f(x) syms x diff(f(x), x) 19 求高階導數(shù) diff(f(x)，n) diff(f(x),x，n)或者：syms x diff(f(x)，n)syms x diff(f(x), x，n) 20 在MATLAB中沒有直接求隱函數(shù)導數(shù)的命令，但是我們可以根據(jù)數(shù)學中求隱函數(shù)導數(shù)的方法，在中一步一步地進行推

14、導；也可以自己編一個求隱函數(shù)導數(shù)的小程序；不過，最簡便的方法是調用Maple中求隱函數(shù)導數(shù)的命令，調用格式如下： maple(implicitdiff(f(x,y)=0,y,x) 在MATLAB中，沒有直接求參數(shù)方程確定的函數(shù)的導數(shù)的命令，只能根據(jù)參數(shù)方程確定的函數(shù)的求導公式 一步一步地進行推導；或者，干脆自己編一個小程序，應用起來會更加方便。21 求不定積分 int(f(x) int (f(x),x)或者：syms x int(f(x) syms x int(f(x), x) 22 求定積分、廣義積分 int(f(x),a,b) int (f(x),x,a,b)或者：syms x int(f

15、(x),a,b) syms x int(f(x), x,a,b) 23 進行換元積分的計算自身沒有提供這一功能，但是可以調用Maple函數(shù)庫中的changevar命令，調用方法如下：maple( with(student) ) 加載student函數(shù)庫后，才能使用changevar命令maple( changevar( m(x)=p(u), Int(f(x),x) ) ) 把積分表達式中的m(x)代換成p(u)24 進行分部積分的計算自身沒有提供這一功能，但是可以調用Maple函數(shù)庫中的intparts命令，調用方法如下： maple( with(student) ) 加載student函數(shù)庫

16、后，才能使用intparts命令maple(intparts(Int(f(x),x),u) ) 指定u，用分部積分公式 進行計算25 對數(shù)列和級數(shù)進行求和 syms n symsum(f(n), n a ,b )26 進行連乘 maple(product(f(n),n=a.b)27 展開級數(shù)syms x taylor(f(x), x, n, a )28 進行積分變換syms s t laplace( f(t), t, s ) 拉普拉斯變換 ilaplace( F(s), s, t ) 拉普拉斯變換的逆變換 syms t fourier( f(t), t, ) 傅立葉變換 ifourier( F(), , t ) 傅立葉變換的逆變換 syms n z ztrans(