《MATLAB與控制系統(tǒng)仿真》09自動(dòng)化(共43頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上MATLAB與控制系統(tǒng)仿真實(shí)驗(yàn)報(bào)告班級(jí)： 自動(dòng)化092班 學(xué)號(hào)： 姓名： 張曉梅 時(shí)間：2012 年 6 月目 錄實(shí)驗(yàn)一 MATLAB環(huán)境的熟悉與基本運(yùn)算（一）實(shí)驗(yàn)二 MATLAB環(huán)境的熟悉與基本運(yùn)算（二）實(shí)驗(yàn)三 MATLAB語(yǔ)言的程序設(shè)計(jì)實(shí)驗(yàn)四 MATLAB的圖形繪制實(shí)驗(yàn)五 基于SIMULINK的系統(tǒng)仿真實(shí)驗(yàn)六 控制系統(tǒng)的頻域與時(shí)域分析實(shí)驗(yàn)七 控制系統(tǒng)PID校正器設(shè)計(jì)法實(shí)驗(yàn)八 線性方程組求解及函數(shù)求極值南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)： 一、實(shí)驗(yàn)項(xiàng)目名稱實(shí)驗(yàn)一 MATL

2、AB環(huán)境的熟悉與基本運(yùn)算（一）二、實(shí)驗(yàn)?zāi)康模?1 熟悉MATLAB開(kāi)發(fā)環(huán)境2掌握矩陣、變量、表達(dá)式的各種基本運(yùn)算三、主要儀器設(shè)備及耗材計(jì)算機(jī)四、實(shí)驗(yàn)基本知識(shí)： 1.熟悉MATLAB環(huán)境: MATLAB桌面和命令窗口、命令歷史窗口、幫助信息瀏覽器、工作空間瀏覽器、文件和搜索路徑瀏覽器。2.掌握MATLAB常用命令表1 MATLAB常用命令clc 清除命令窗口中內(nèi)容clear 清除工作空間中變量help 對(duì)所選函數(shù)的功能、調(diào)用格式及相關(guān)函數(shù)給出說(shuō)明lookfor 查找具有某種功能的函數(shù)但卻不知道該函數(shù)的準(zhǔn)確名稱3.MATLAB變量與運(yùn)算符31變量命名規(guī)則32 MATLAB的各種常用運(yùn)算符表2 MA

3、TLAB算術(shù)運(yùn)算符操作符 功能說(shuō)明 操作符 功能說(shuō)明 + 加 矩陣左除- 減. 數(shù)組左除* 矩陣乘/ 矩陣右除.* 數(shù)組乘./ 數(shù)組右除 矩陣乘方' 矩陣轉(zhuǎn)置. 數(shù)組乘方.' 數(shù)組轉(zhuǎn)置表3 MATLAB關(guān)系運(yùn)算符操作符 功能說(shuō)明 = 等于= 不等于> 大于< 小于>= 大于等于<= 小于等于表4 MATLAB邏輯運(yùn)算符邏輯運(yùn)算符 邏輯運(yùn)算 說(shuō)明 & And 邏輯與| Or 邏輯或 Not 邏輯非Xor邏輯異或表5 MATLAB特殊運(yùn)算符號(hào) 功能說(shuō)明示例 符號(hào) 功能說(shuō)明示例 ： 1:1:4;1:2:11 . ； 分隔行 . ， 分隔列 （） % 注

4、釋 構(gòu)成向量、矩陣 ！ 調(diào)用操作系統(tǒng)命令 構(gòu)成單元數(shù)組= 用于賦值 4.MATLAB的一維、二維數(shù)組的尋訪表6 子數(shù)組訪問(wèn)與賦值常用的相關(guān)指令格式四、實(shí)驗(yàn)內(nèi)容及實(shí)驗(yàn)結(jié)果 1、新建一個(gè)文件夾（自己的名字命名，在機(jī)器的最后一個(gè)盤符） 2、啟動(dòng)MATLAB，將該文件夾添加到MATLAB路徑管理器中。3、學(xué)習(xí)使用help命令。 Help 加上所要想查詢的函數(shù)，MATLAB會(huì)自動(dòng)生成該函數(shù)功能和調(diào)用的說(shuō)明，有助于了解沒(méi)有接觸的函數(shù)。4、窗口命令l closel close alll clcl hold onl hold off了解其功能和作用，觀察command window、command histo

5、ry和workspace等窗口的變化結(jié)果。、5、工作空間管理命令l whol whosl clear6、隨機(jī)生成一個(gè)2×6的矩陣，觀察command window、command history和workspace等窗口的變化結(jié)果，實(shí)現(xiàn)矩陣左旋90°或右旋90°的功能。>> a=rand(2,6)a = 0.8462 0.2026 0.8381 0.6813 0.8318 0.7095 0.5252 0.6721 0.0196 0.3795 0.5028 0.4289>> b=rot90(a)b = 0.7095 0.4289 0.8318

6、 0.5028 0.6813 0.3795 0.8381 0.0196 0.2026 0.6721 0.8462 0.52527、求高階方程的的根，求高階多項(xiàng)式的值。 >> a=1 2 3;4 5 6;7 8 9;>> p=poly(a)p = 1.0000 -15.0000 -18.0000 -0.0000>> v=roots(p)v = 16.1168 -1.1168 -0.00008、創(chuàng)建一個(gè)二維數(shù)組A（4×8）。查詢數(shù)組A第2行、第3列的元素；查詢數(shù)組A第2行的所有元素；查詢數(shù)組A第6列的所有元素；查詢數(shù)組A按列拉長(zhǎng)形成新的數(shù)組B；查詢數(shù)組

7、A按行拉長(zhǎng)形成新的數(shù)組C。專心-專注-專業(yè)>> a=rand(4,8)a = 0.9501 0.8913 0.8214 0.9218 0.9355 0.0579 0.1389 0.2722 0.2311 0.7621 0.4447 0.7382 0.9169 0.3529 0.2028 0.1988 0.6068 0.4565 0.6154 0.1763 0.4103 0.8132 0.1987 0.0153 0.4860 0.0185 0.7919 0.4057 0.8936 0.0099 0.6038 0.7468>> a(2,3)ans = 0.4447>&

8、gt; a(2,:)ans = 0.2311 0.7621 0.4447 0.7382 0.9169 0.3529 0.2028 0.1988>> a(:,6)ans = 0.0579 0.3529 0.8132 0.0099>> b=a(:)b = 0.9501 0.2311 0.6068 0.4860 0.8913 0.7621 0.4565 0.0185 0.8214 0.4447 0.6154 0.7919 0.9218 0.7382 0.1763 0.4057 0.9355 0.9169 0.4103 0.8936 0.0579 0.3529 0.8132 0

9、.0099 0.1389 0.2028 0.1987 0.6038 0.2722 0.1988 0.0153 0.7468五、實(shí)驗(yàn)心得這是MATLAB的第一個(gè)實(shí)驗(yàn)，也是我第一次接觸這個(gè)軟件，在初步了解中可以感覺(jué)到它的功能很多而且可以解決很多工程問(wèn)題等，是我們專業(yè)很重要的應(yīng)用軟件，它的指令簡(jiǎn)單易懂但是編寫語(yǔ)句時(shí)不能出錯(cuò)，只要是一點(diǎn)不符合軟件編程要求都不能得到想要的結(jié)果，細(xì)心耐心才能做出結(jié)果，使我匪淺！南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)： 一、實(shí)驗(yàn)項(xiàng)目名稱MATLAB環(huán)境的熟悉與基本運(yùn)算（二）二、實(shí)驗(yàn)?zāi)?/p>

10、的： 1 熟悉MATLAB開(kāi)發(fā)環(huán)境2掌握矩陣、變量、表達(dá)式的各種基本運(yùn)算三、主要儀器設(shè)備及耗材計(jì)算機(jī)四、實(shí)驗(yàn)基本知識(shí)： 1.MATLAB的基本運(yùn)算表1 兩種運(yùn)算指令形式和實(shí)質(zhì)內(nèi)涵的異同表2MATLAB的常用函數(shù)五、實(shí)驗(yàn)內(nèi)容1、設(shè)有3個(gè)二維數(shù)組A（2×4），B（2×4），C（2×2），參照表1寫出所有由2個(gè)數(shù)組參與的合法的數(shù)組運(yùn)算和矩陣指令及結(jié)果。2、學(xué)習(xí)使用表2列出的常用函數(shù)。3、學(xué)習(xí)使用表3列出的數(shù)組操作函數(shù)。4、生成下列矩陣，并取出方框內(nèi)的數(shù)組元素六、實(shí)驗(yàn)結(jié)果答：1.A:數(shù)組的運(yùn)算 >> a=zeros(2,4); >> b=ones(

11、2,4);>> c=eye(2,2)>> aa = 0 0 0 0 0 0 0 0>> bb = 1 1 1 1 1 1 1 1>> cc = 1 0 0 1>> a+bans = 1 1 1 1 1 1 1 1>> a-bans = -1 -1 -1 -1 -1 -1 -1 -1>> a.*bans = 0 0 0 0 0 0 0 0>> a./bans = 0 0 0 0 0 0 0 0>> b.aans = 0 0 0 0 0 0 0 0>> exp(a)ans = 1

12、 1 1 1 1 1 1 1>> log(b)ans = 0 0 0 0 0 0 0 0>> sqrt(c)ans = 1 0 0 1>> c'ans = 1 0 0 1B：矩陣的運(yùn)算>> a=1 1;2 2a = 1 1 2 2>> a'ans = 1 2 1 2>> 5*aans = 5 510 10>> a2ans = 3 3 6 6>> a1.2ans = 1.2457 1.2457 2.4915 2.4915>> b=1 2;3 4b = 1 2 3 4>

13、> 2*inv(b)ans = -4.0000 2.0000 3.0000 -1.0000>> 4*inv(b)ans = -8.0000 4.0000 6.0000 -2.0000>> expm(b)ans = 51.9690 74.7366 112.1048 164.0738> b=1 2 ;3 4;>> logm(b)Warning: Principal matrix logarithm is not defined for A with nonpositive real eigenvalues. A non-principal matri

14、x logarithm is returned.> In funm at 153 In logm at 27ans = -0.3504 + 2.3911i 0.9294 - 1.0938i 1.3940 - 1.6406i 1.0436 + 0.7505i>> sqrtm(b)ans = 0.5537 + 0.4644i 0.8070 - 0.2124i 1.2104 - 0.3186i 1.7641 + 0.1458i2、表2常用函數(shù)的使用>> V=1 2 3 4;>> A=diag(V)A = 1 0 0 0 0 2 0 0 0 0 3 0 0 0

15、 0 4>> a=eye(3)a = 1 0 0 0 1 0 0 0 1>> b=magic(4)b = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1>> c=ones(3)c = 1 1 1 1 1 1 1 1 1>> a=rand(3)a = 0.9501 0.4860 0.4565 0.2311 0.8913 0.0185 0.6068 0.7621 0.8214>> a=randn(3)a = -0.4326 0.2877 1.1892 -1.6656 -1.1465 -0.0376 0.125

16、3 1.1909 0.3273>> a=zeros(3)a = 0 0 0 0 0 0 0 0 0>> a=1 2 3,4 5 6;>> b=size(a)b = 1 63、表3函數(shù)的使用>> a=1 2 3 4;>> b=diag(a)b = 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4>> c=flipud(a)c = 1 2 3 4>> d=fliplr(a)d = 4 3 2 1>> a=1 2 3;>> rot90(a)ans = 3 2 1>>

17、q=1 2 3;1 2 3;1 2 3;>> det(q)ans = 0>> a=1 2;2 1;>> rank(a)ans = 2>> a=1 1;2 2;>> rank(a)ans = 14、取出方塊中的元素>> A=1:25;A=reshape(A,5,5)'>> A(2,:)=fliplr(A(2,:);>> AA = 1 2 3 4 5 10 9 8 7 6 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25>> A(2,2:3)a

18、ns = 9 8>> A(2:4,4)ans = 7 14 19>> A(4:5,1:3)ans = 16 17 18 21 22 23七、實(shí)驗(yàn)心得：這個(gè)實(shí)驗(yàn)室考察矩陣和向量的運(yùn)算，在其中有很多細(xì)節(jié)要注意，例如，在矩陣的運(yùn)算中，最容易出現(xiàn)的問(wèn)題是求逆矩陣的運(yùn)算，所要求的矩陣必須是非奇異的矩陣，否則會(huì)產(chǎn)生錯(cuò)誤的結(jié)果。南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)： 一、實(shí)驗(yàn)項(xiàng)目名稱實(shí)驗(yàn)三 MATLAB語(yǔ)言的程序設(shè)計(jì)二、實(shí)驗(yàn)?zāi)康模?1、熟悉MATLAB程序編輯與設(shè)計(jì)環(huán)境2、掌握各種編程語(yǔ)句語(yǔ)

19、法規(guī)則及程序設(shè)計(jì)方法 3、函數(shù)文件的編寫和設(shè)計(jì) 4、了解和熟悉跨空間變量傳遞和賦值 三、主要儀器設(shè)備及耗材計(jì)算機(jī)四、實(shí)驗(yàn)基本知識(shí)：1、程序流程控制語(yǔ)句 for循環(huán)結(jié)構(gòu)語(yǔ)法：for i初值：增量：終值 語(yǔ)句1 語(yǔ)句n end 說(shuō)明：1i初值：終值，則增量為1。 2初值、增量、終值可正可負(fù)，可以是整數(shù)，也可以是小數(shù)，只須符合數(shù)學(xué)邏輯。while 循環(huán)結(jié)構(gòu) 語(yǔ)法：while 邏輯表達(dá)式 循環(huán)體語(yǔ)句 end 說(shuō)明：1、whiIe結(jié)構(gòu)依據(jù)邏輯表達(dá)式的值判斷是否執(zhí)行循環(huán)體語(yǔ)勾。若表達(dá)式的值為真，執(zhí)行循環(huán)體語(yǔ)句一次、在反復(fù)執(zhí)行時(shí)，每次都要進(jìn)行判斷。若表達(dá)式的值為假，則程序執(zhí)行end之后的語(yǔ)句。2、為了避免因

20、邏輯上的失誤，而陷入死循環(huán)，建議在循環(huán)體語(yǔ)句的適當(dāng)位置加break語(yǔ)句、以便程序能正常執(zhí)行。（執(zhí)行循環(huán)體的次數(shù)不確定；每一次執(zhí)行循環(huán)體后，一定會(huì)改變while后面所跟關(guān)系式的值。） 3、while循環(huán)也可以嵌套、其結(jié)構(gòu)如下： while邏輯表達(dá)式1 循環(huán)體語(yǔ)句1 while邏輯表達(dá)式2 循環(huán)體語(yǔ)句2 end 循環(huán)體語(yǔ)句3 end if-else-end分支結(jié)構(gòu)if 表達(dá)式 1 語(yǔ)句1 else if 表達(dá)式 2(可選) 語(yǔ)句2 else(可選) 語(yǔ)句3 end end 說(shuō)明：1.if結(jié)構(gòu)是一個(gè)條件分支語(yǔ)句，若滿足表達(dá)式的條件，則往下執(zhí)行；若不滿足，則跳出if結(jié)構(gòu)。2else if表達(dá)式2與el

21、se為可選項(xiàng)，這兩條語(yǔ)句可依據(jù)具體情況取舍。3注意：每一個(gè)if都對(duì)應(yīng)一個(gè)end，即有幾個(gè)if，記就應(yīng)有幾個(gè)end。switch-case結(jié)構(gòu)語(yǔ)法：switch表達(dá)式case常量表達(dá)式1 語(yǔ)句組1 case常量表達(dá)式2 語(yǔ)句組2 otherwise 語(yǔ)句組n end 說(shuō)明：1switch后面的表達(dá)式可以是任何類型，如數(shù)字、字符串等。2當(dāng)表達(dá)式的值與case后面常量表達(dá)式的值相等時(shí)，就執(zhí)行這個(gè)case后面的語(yǔ)句組如果所有的常量表達(dá)式的值都與這個(gè)表達(dá)式的值不相等時(shí)，則執(zhí)行otherwise后的執(zhí)行語(yǔ)句。程序流程控制指令break、return、pause 2、函數(shù)文件的結(jié)構(gòu)和編寫方法 函數(shù)只能在函數(shù)

22、體內(nèi)對(duì)變量進(jìn)行操作，也就是只能訪問(wèn)函數(shù)本身工作空間中的變量。M函數(shù)文件的結(jié)構(gòu)（1） 函數(shù)定義行（function） （2） H1行（函數(shù)幫助文本的第一行） （3） 函數(shù)幫助文本（4） 函數(shù)體（5） 注釋函數(shù)文件編寫后，保存時(shí)，其文件名必須與函數(shù)名相同。注意：函數(shù)名不要與MATLAB自身的函數(shù)命令相同。五、實(shí)驗(yàn)內(nèi)容： 練習(xí)A 1、熟悉MATLAB程序編輯與設(shè)計(jì)環(huán)境2、用for循環(huán)語(yǔ)句實(shí)現(xiàn)求1100的和3、用for循環(huán)語(yǔ)句實(shí)現(xiàn)編寫一個(gè)求階乘n!的函數(shù)文件4、判斷y向量的元素屬于3n 3n+1或3n+2，設(shè)y=3,5,7,9,11,16。(使用函數(shù)rem求余數(shù))7、已知數(shù)組A=2,4,6,8,10，

23、B=4,9,6,7,4,求,其中n=5。六、實(shí)驗(yàn)結(jié)果：1、列出實(shí)驗(yàn)程序代碼及實(shí)驗(yàn)結(jié)果。l 用for循環(huán)語(yǔ)句實(shí)現(xiàn)求1100的和>> %ex1.m>> mysum=0;i=1;>> for i=1:100 mysum=mysum+i;end>> mysummysum = 5050l 用for循環(huán)寫n！函數(shù)函數(shù)是：%at=1;for i=1:n; t=i*t;endt運(yùn)行結(jié)果為：>> n=3;at = 6>> n=4;at = 24l 判斷y向量的元素屬于3n 3n+1或3n+2，設(shè)y=3,5,7,9,11,16。(使用函數(shù)re

24、m求余數(shù))>> y=3,5,7,9,11,16y = 3 5 7 9 11 16>> a=rem(y,3)a = 0 2 1 0 2 1l 已知數(shù)組A=2,4,6,8,10，B=4,9,6,7,4,求,其中n=5程序文件如下：%wA=2,4,6,8,10;B=4,9,6,7,4;m=0;for i=1:5 m=m+A(i)*B(6-i)end程序調(diào)試結(jié)果：>> wm = 8m = 36m = 72m = 144m = 184七寫出實(shí)驗(yàn)體會(huì)。這個(gè)實(shí)驗(yàn)是關(guān)于小程序的編寫，主要考察語(yǔ)句的使用和整體程序的理解，需要有語(yǔ)言編寫的能力，而且需要耐心調(diào)試修改。 南昌大學(xué)實(shí)

25、驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)： 一、實(shí)驗(yàn)項(xiàng)目名稱實(shí)驗(yàn)四 MATLAB的圖形繪制二、實(shí)驗(yàn)?zāi)康模?1、學(xué)習(xí)MATLAB圖形繪制的基本方法； 2、熟悉和了解MATLAB圖形繪制程序編輯的基本指令； 3、熟悉掌握利用MATLAB圖形編輯窗口編輯和修改圖形界面，并添加圖形的各種標(biāo)注； 4、掌握plot、subplot的指令格式和語(yǔ)法。三、主要儀器設(shè)備及耗材計(jì)算機(jī)四、實(shí)驗(yàn)基本知識(shí)：1、基本的繪圖命令plot（x,y,）2、建立圖形窗口命令figure（1）；figure（2）；figure(n)打開(kāi)不同的圖形窗口，

26、以便繪制不同的圖形。3、grid on：在所畫出的圖形坐標(biāo)中加入柵格grid off：除去圖形坐標(biāo)中的柵格。4、hold on：把當(dāng)前圖形保持在屏幕上不變，同時(shí)允許在這個(gè)坐標(biāo)內(nèi)繪制另外一個(gè)圖形。hold off：使新圖覆蓋舊的圖形。5、axis設(shè)定軸的范圍axis（xmin xmax ymin ymax）設(shè)定x 軸與y軸的最大、最小坐標(biāo)。axis(equal)：將x坐標(biāo)軸和y坐標(biāo)軸的單位刻度大小調(diào)整為一樣。6、文字標(biāo)示l text(x,y,字符串)在圖形的指定坐標(biāo)位置(x,y)處，標(biāo)示單引號(hào)括起來(lái)的字符串。l gtext(字符串)利用鼠標(biāo)在圖形的任意位置標(biāo)示字符串。l title(字符串)在所

27、畫圖形的最上端顯示說(shuō)明該圖形標(biāo)題的字符串。l xlabel(字符串)，ylabel(字符串)設(shè)置x，y坐標(biāo)軸的名稱。l 輸入特殊的文字需要用反斜杠（）開(kāi)頭7、subplot（m，n，k）：分割圖形顯示窗口，在同一個(gè)窗口中顯示多個(gè)圖形。l m：上下分割個(gè)數(shù)l n：左右分割個(gè)數(shù)l k:：子圖編號(hào)8、semilogx：繪制以x軸為對(duì)數(shù)坐標(biāo)（以10為底），y軸為線性坐標(biāo)的半對(duì)數(shù)坐標(biāo)圖形。semilogy：繪制以y 軸為對(duì)數(shù)坐標(biāo)（以10 為底），x 軸為線性坐標(biāo)的半對(duì)數(shù)坐標(biāo)圖形。五、實(shí)驗(yàn)內(nèi)容1、繪制下列曲線，要求盡可能多地使用二中所列命令。 y=sin(t)sin(9t)2、繪制二階系統(tǒng)階躍響應(yīng)，綜合演

28、示圖形標(biāo)識(shí)。clf;t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t); tt=t(find(abs(y-1)>0.05);ts=max(tt); plot(t,y,'r-','LineWidth',3) axis(-inf,6*pi,0.6,inf) set(gca,'Xtick',2*pi,4*pi,6*pi,'Ytick',0.95,1,1.05,max(y) grid on title('it y = 1 - e -alphatcosomegat') text

29、(13.5,1.2,'fontsize12alpha=0.3') text(13.5,1.1,'fontsize12omega=0.7') hold on;plot(ts,0.95,'bo','MarkerSize',10);hold off cell_string1='fontsize12uparrow' cell_string2='fontsize16 fontname隸書鎮(zhèn)定時(shí)間' cell_string3='fontsize6 ' cell_string4='font

30、size14rmt_s = ' num2str(ts); text(ts,0.85,cell_string) xlabel('fontsize14 bft rightarrow')ylabel('fontsize14 bfy rightarrow')六、實(shí)驗(yàn)結(jié)果1列出實(shí)驗(yàn)程序（命令）代碼及實(shí)驗(yàn)結(jié)果。試驗(yàn)程序如下：>> x=0:0.1:2*pi;y=x+(x.3)/3;plot(x,y)試驗(yàn)程序如下：x=0:0.1:2*pi;y=(1/2.*pi).*exp(x.2/2);plot(x,y) 試驗(yàn)程序如下：t=0:0.1:2*pi;r=2.5;

31、x=r*sin(t);y=r*cos(t);plot(x,y) y=sin(t)sin(9t)試驗(yàn)程序如下：x=0:0.1:2*pi;y=sin(x).*sin(9*x);plot(x,y)實(shí)驗(yàn)圖一實(shí)驗(yàn)圖二：實(shí)驗(yàn)圖三：實(shí)驗(yàn)圖四七、心得體會(huì)： 這個(gè)實(shí)驗(yàn)是首次運(yùn)用MATLAB的繪圖功能，需要先寫出原函數(shù)，在窗口中調(diào)用，主要是考察plot函數(shù)的使用，難點(diǎn)在于源文件的而編寫，在整個(gè)程序中不能出現(xiàn)一點(diǎn)差錯(cuò)就不能得到想要的圖，考驗(yàn)編程能力和耐心細(xì)心程度。 南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)： 一、實(shí)驗(yàn)項(xiàng)目名稱實(shí)

32、驗(yàn)五 基于SIMULINK的系統(tǒng)仿真二、實(shí)驗(yàn)?zāi)康模?1、熟悉SIMULINK 工作環(huán)境及特點(diǎn)2、掌握線性系統(tǒng)仿真常用基本模塊的用法3、掌握SIMULINK 的建模與仿真方法4、子系統(tǒng)的創(chuàng)建和封裝設(shè)計(jì)三、主要儀器設(shè)備及耗材計(jì)算機(jī)四、實(shí)驗(yàn)基本知識(shí)：1、了解SIMULINK模塊庫(kù)中各子模塊基本功能 2、SIMULINK 的建模與仿真方法（1）打開(kāi)模塊庫(kù)，找出相應(yīng)的模塊。鼠標(biāo)左鍵點(diǎn)擊相應(yīng)模塊，拖拽到模型窗口中即可。（2）創(chuàng)建子系統(tǒng)：當(dāng)模型大而復(fù)雜時(shí)，可創(chuàng)建子系統(tǒng)。（3）設(shè)置仿真控制參數(shù)。五、實(shí)驗(yàn)內(nèi)容及實(shí)驗(yàn)結(jié)果 1、SIMULINK仿真實(shí)際應(yīng)用（1） 雙環(huán)調(diào)速的電流環(huán)系統(tǒng)的方框圖模型為： 圖中參數(shù)設(shè)為K

33、s=44；Ts=0.00167；Ta=0.017；R=1；Tm=0.075；Ce=0.1925；Kt=0.01178；T1=0.049；T2=0.088(1)在Simulink集成環(huán)境下建立模型，在給定信號(hào)作用點(diǎn)處輸入單位給定階躍響應(yīng)信號(hào)，0.3秒后在擾動(dòng)信號(hào)點(diǎn)輸入單位階躍響應(yīng)信號(hào)。并繪制相應(yīng)的響應(yīng)曲線Ks=44;Ts=0.00167;Ta=0.017;R=1;Tm=0.075;Ce=0.1925;Kt=0.01178;T1=0.049;T2=0.088;>> qingchundou2、用Simulink對(duì)以下系統(tǒng)進(jìn)行仿真 其中為系統(tǒng)輸入，y(t)為系統(tǒng)輸出，仿真當(dāng)輸入為正弦信號(hào)時(shí)

34、，輸出的信號(hào)的波形，仿真時(shí)間0<=0t<=100。仿真的系統(tǒng)模型：結(jié)果顯示：3、在滑艇的運(yùn)行過(guò)程中，滑艇主要受到如下作用力的控制：滑艇自身的牽引力F，滑艇受到的水的阻力f。其中水的阻力，u為滑艇的運(yùn)動(dòng)速度。由運(yùn)動(dòng)學(xué)的相關(guān)定理可知，整個(gè)滑艇系統(tǒng)的動(dòng)力學(xué)方程為：其中，m為滑艇的質(zhì)量。假設(shè)滑艇的質(zhì)量為1000kg，建立此系統(tǒng)的Simulink模型并進(jìn)行分析。解：（1）滑艇速度控制系統(tǒng)的模型建立與仿真仿真模型建立如下圖中的step模塊，用來(lái)產(chǎn)生滑艇的牽引力。將其終值設(shè)為1000，起始時(shí)間設(shè)置為0,增益gain設(shè)為1/1000，函數(shù)模塊用于取水的阻力，設(shè)置為u2-u，scope模塊用于顯示滑

35、艇的速度。根據(jù)模型仿真出來(lái)的圖形如下：由仿真曲線可知，滑艇在牽引力F的作用下，在經(jīng)過(guò)80s左右的時(shí)間后，速度由0上升并穩(wěn)定在33km/h.六、實(shí)驗(yàn)心得; 這個(gè)試驗(yàn)中我學(xué)習(xí)到了matlab內(nèi)很多的函數(shù)構(gòu)造和系統(tǒng)模型的建立，并用示波器模擬出實(shí)驗(yàn)結(jié)果的波形，使我了解到用計(jì)算機(jī)進(jìn)行模擬的方法，同時(shí)可以看出matlab這個(gè)軟件的應(yīng)用廣泛和功能的多樣強(qiáng)大南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)： 一、實(shí)驗(yàn)項(xiàng)目名稱實(shí)驗(yàn)六 控制系統(tǒng)的頻域與時(shí)域分析二、實(shí)驗(yàn)?zāi)康模?1、掌握控制系統(tǒng)數(shù)學(xué)模型的基本描述方法和相互轉(zhuǎn)化2、了解控

36、制系統(tǒng)的穩(wěn)定性分析方法3、掌握控制系統(tǒng)頻域與時(shí)域分析基本方法三、主要儀器設(shè)備及耗材計(jì)算機(jī)四、實(shí)驗(yàn)基本知識(shí)：1、系統(tǒng)數(shù)學(xué)模型的幾種表示方法 2、在MATLAB工具箱中，提供了子系統(tǒng)的連接處理函數(shù)：1） series()函數(shù)：系統(tǒng)串聯(lián)實(shí)現(xiàn)。2） parallel()函數(shù)：系統(tǒng)并聯(lián)實(shí)現(xiàn)。3） feedback()函數(shù)：系統(tǒng)反饋連接。3、控制系統(tǒng)根軌跡繪制1）rlocfind():計(jì)算給定根的根軌跡增益 2）rlocus() 函數(shù)：功能為求系統(tǒng)根軌跡。 4、控制系統(tǒng)頻域分析基本方法1）Bode圖：bode()函數(shù)2）Nyquist圖：nyquist()函數(shù)3）穩(wěn)定裕度計(jì)算：margin( )函數(shù)5、線

37、性系統(tǒng)時(shí)間響應(yīng)分析1）step( )函數(shù)-求系統(tǒng)階躍響應(yīng)2）impulse( )函數(shù)：求取系統(tǒng)的脈沖響應(yīng)3）lsim( )函數(shù)：求系統(tǒng)的任意輸入下的仿真五、實(shí)驗(yàn)內(nèi)容及實(shí)驗(yàn)結(jié)果1、表示下列傳遞函數(shù)模型,并轉(zhuǎn)化成其他的數(shù)學(xué)模型(1)>> K=4;Z=-2;-1;-5;-1;-5;P=0;-1;-1;-1;-0.0479 + 1.3112i;-0.0479 - 1.3112iP = 0 -1.0000 -1.0000 -1.0000 -0.0479 + 1.3112i -0.0479 - 1.3112i>> A,B,C,D=zp2ss(Z,P,K),num,den=zp2tf

38、(Z,P,K)A = -1.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 1.0000 -0.0958 -1.3121 0 0 0 0 1.3121 0 0 0 1.0000 1.0000 9.9042 17.7417 -2.0000 -1.0000 0 0 0 0 1.0000 0B = 1 0 0 0 0 0C = 4.0000 4.0000 39.6168 70.9669 4.0000 4.0000D = 0num = 0 4 56 280 608 580 200den = 1.0000 3.0958 5.0089 6.4520 5.2604 1.7215

39、 0(2)>> num=0 4 -2;den=1 0 2 5;printsys(num,den) num/den = 4 s - 2 - s3 + 2 s + 5>> num=0 4 -2;den=1 0 2 5;Z,P K=tf2zp(num,den)Z = 0.5000P = 0.6641 + 1.8230i 0.6641 - 1.8230i -1.3283 K = 4（3）>> k=1,p=-1/2;-2k = 1p = -0.5000 -2.0000>> K=1;P=-1/2;-2;Z=0;>> A,B,C,D=zp2ss(Z

40、,P,K),num,den=zp2tf(Z,P,K)A = -2.5000 -1.0000 1.0000 0B = 1 0C = 1 0D = 0num = 0 1 0den = 1.0000 2.5000 1.0000(4)>> A=0 1;-1 -2A = 0 1 -1 -2>> B=0;1B = 0 1>> C=0 1, D=zeros(2,2)C = 0 1D = 0 0 0 02、一個(gè)單位負(fù)反饋開(kāi)環(huán)傳遞函數(shù)為 試?yán)L出系統(tǒng)閉環(huán)的根軌跡圖；并在跟軌跡圖上任選一點(diǎn)，試計(jì)算該點(diǎn)的增益K及其所有極點(diǎn)的位置。解：試驗(yàn)程序是：%gsnum=1;den=conv(

41、1,0,conv(0.5,1,4,1);rlocus(num,den),k,poles=rlocfind(num,den)得到的根軌跡圖如下：在圖中計(jì)算出幾點(diǎn)的位置是：>> gsSelect a point in the graphics windowselected_point = 0 + 0.6832ik = 2.1011poles = -2.2365 -0.0068 + 0.6853i -0.0068 - 0.6853i在根軌跡上任意取了一點(diǎn)的值>> gsSelect a point in the graphics windowselected_point = 0

42、.6825 + 2.2236ik = 39.0822poles = -3.6105 0.6803 + 2.2247i 0.6803 - 2.2247i3、求下面系統(tǒng)在階躍信號(hào)為0.11(t)時(shí)系統(tǒng)的響應(yīng)。并求系統(tǒng)性能指標(biāo)：穩(wěn)態(tài)值、上升時(shí)間、調(diào)節(jié)時(shí)間、超調(diào)量。% gnum=20;den=1 8 36 40 20;numc,denc=cloop(num,den);t=0:0.1:10;y,x,t=step(numc,denc,t);plot(t,y)六、實(shí)驗(yàn)心得傳遞函數(shù)是現(xiàn)代控制理論中很重要的數(shù)學(xué)模型，我以前只學(xué)習(xí)到用數(shù)學(xué)和物理的方法分析它的極點(diǎn)和特征根，根據(jù)所計(jì)算出的數(shù)據(jù)描繪出根軌跡然后在進(jìn)行系統(tǒng)的分析，而學(xué)習(xí)了matlab這個(gè)軟件后可以用軟件編程的方法實(shí)現(xiàn)以前計(jì)算的內(nèi)容，簡(jiǎn)單直觀效率高，是現(xiàn)代控制中很重要的工具。 南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名： 張曉梅 學(xué) 號(hào)： 專業(yè)班級(jí)： 自動(dòng)化092班 實(shí)驗(yàn)類型： 驗(yàn)證 綜合 設(shè)計(jì) 創(chuàng)新 實(shí)驗(yàn)日