交通運輸自動控制原理課程設(shè)計報告

1、<p>　　《交通運輸自動控制原理》</p><p><b>　　課程設(shè)計報告</b></p><p>　　二〇一四年五月二十二日</p><p>　本組組員</p><p>　指導(dǎo)教師</p><p><b>　　目 錄 </b></p><

2、;p>　　1 自動控制系統(tǒng)建模設(shè)計1</p><p>　　1.1系統(tǒng)建模題目11</p><p>　　1.1.1基本構(gòu)成及機理分析2</p><p>　　1.1.2反饋控制結(jié)構(gòu)模型及工作原理2</p><p>　　2 自動控制系統(tǒng)仿真實驗分析7</p><p>　　2.1系統(tǒng)仿真實驗題目17<

3、;/p><p>　　2.1.1實驗步驟7</p><p>　　2.1.2實驗分析計算依據(jù)7</p><p>　　2.1.2結(jié)論圖表數(shù)據(jù)分析7</p><p>　　2.2系統(tǒng)仿真實驗題目28</p><p>　　2.3系統(tǒng)仿真實驗題目38</p><p>　　附錄A 程序源代碼40<

4、;/p><p><b>　　附錄B 42</b></p><p><b>　　參考文獻44</b></p><p>　　自動控制系統(tǒng)建模設(shè)計</p><p><b>　　系統(tǒng)建模題目1</b></p><p>　　查閱控制系統(tǒng)的相關(guān)專業(yè)資料，運用自動控

5、制系統(tǒng)的基本構(gòu)成及反饋控制原理對通控制系統(tǒng)進行機理分析，建立其反饋控制結(jié)構(gòu)模型，闡明系統(tǒng)的工作原理。</p><p>　　我們選擇ATO系統(tǒng)進行研究。ATO系統(tǒng)即列車自動駕駛子系統(tǒng)(ATO-Automatic TrainOperation)。是列車控制系統(tǒng)(automatic train control，ATC)重要的子系統(tǒng)。ATO系統(tǒng)能代替司機駕駛列車，使列車平穩(wěn)地加速到運行速度，并根據(jù)接收到的ATP限制速度命

6、令和ATO到站停車速度曲線自動調(diào)整列車速度，實現(xiàn)列車牽引加速、勻速惰行、減速制動控制和精確停車等基本駕駛功能。ATO系統(tǒng)的功能分為基本控制功能和服務(wù)功能?；究刂乒δ馨ǎ鹤詣玉{駛、自動折返、自動控制車門、屏蔽門開閉。這三個基本控制功能相互之間獨立地運行。服務(wù)及其他功能包括：故障自診斷、運行信息記錄、列車位置、允許速度、巡航/惰行功能、PTI(列車識別系統(tǒng))支持功能。 在這里我們主要分析ATO中自動駕駛（即自動調(diào)整列車速度）的功

7、能，研究其自動控制的結(jié)構(gòu)模型。 </p><p>　　1.1.1 基本構(gòu)成及機理分析</p><p>　　ATO系統(tǒng)的速度調(diào)整系統(tǒng)構(gòu)成：由運行等級模式曲線計算模塊、PID速度控制器和受控對象（列車運行模型）共同組成。運行等級模式曲線的計算模塊根據(jù)當(dāng)前區(qū)間情況離線計算運行等級模式曲線，PID速度控制器根據(jù)選定的運行等級，以該等級的運行等級模式的目標(biāo)速度一距離曲線為目標(biāo)曲線，輸出速度信號施加

8、給受控對象（列車運行模型），控制列車快速平穩(wěn)的跟隨給定的目標(biāo)速度曲線。</p><p>　　其速度控制系統(tǒng)的框圖如圖所示。</p><p>　　1.1.2 反饋控制結(jié)構(gòu)模型及工作原理</p><p>　　上圖為ATO控制系統(tǒng)的反饋結(jié)構(gòu)模型。PID速度控制器即PID（比例（proportion）、積分（integration）、微分（differentiation）

9、）控制器，由比例單元（P）、積分單元（I）和微分單元（D）組成，傳遞函數(shù)為：</p><p>　　Kp 、Ki、Kd三個參數(shù)常量通過系統(tǒng)整定得到。</p><p>　　受控對象即為列車模型，這里我們?nèi)×熊嚹Ｐ偷膫鬟f函數(shù)為：</p><p>　　所以整個閉環(huán)系統(tǒng)的傳遞函數(shù)為：</p><p>　　通過相關(guān)參數(shù)整定方法，我們可以得知：</p

10、><p>　　工作原理：在求和單元，來自運行等級模式曲線計算模塊的目標(biāo)速度信號和實際速度信號相加，共同輸入PID速度控制器，這個控制器把收集到的實際速度信號和一個目標(biāo)速度信號進行比較，然后把這個差別用于計算新的輸入值，這個新的輸入值的目的是可以讓系統(tǒng)的信號達(dá)到或者保持在目標(biāo)速度。接著新的輸入信號傳遞到列車，通過相關(guān)設(shè)備調(diào)整列車實際運行速度，最后設(shè)備觀測的實際速度信號又反饋回求和單元。</p><p

11、>　　整個控制過程即反饋控制原理：由相關(guān)設(shè)備獲得被控對象（列車）的實際速度信號，通過PID控制器比較計算，輸出調(diào)整信號又作用于被控對象調(diào)整實際運行速度。</p><p>　　自動控制系統(tǒng)仿真實驗分析</p><p>　　2．1 系統(tǒng)仿真實驗題目1</p><p>　　使用Matlab軟件，對下圖所示交通控制系統(tǒng)的特性進行仿真實驗分析。</p>

12、<p>　?。?）繪制K1=200時的開環(huán)對數(shù)頻率特性圖(Bode圖)，仿真分析閉環(huán)系統(tǒng)穩(wěn)定裕度。</p><p>　　（2）繪制關(guān)于參數(shù)K1的閉環(huán)系統(tǒng)的根軌跡圖，仿真分析閉環(huán)系統(tǒng)的穩(wěn)定性。</p><p>　　2.1.1 實驗步驟</p><p>　　(1) 第一步：根據(jù)題目中的框圖圖求出其開環(huán)傳遞函數(shù)。 </p><p>

13、;　　第二步：由求出來的開環(huán)傳遞函數(shù)運用Matlab軟件畫出Bode圖。第三步：然后用Matlab對bode圖求解穩(wěn)定裕度。</p><p>　?。?） 第一步：由題求出其閉環(huán)系統(tǒng)的開環(huán)傳遞函數(shù)第二步：運用函數(shù)語句rlocus(sys)由開環(huán)傳遞函數(shù)畫出其閉環(huán)系統(tǒng)的根軌跡圖。第三步：分析根軌跡圖，討論系統(tǒng)穩(wěn)定性。</p><p><b>　　實驗分析計算依據(jù)</b

14、></p><p>　　由框圖計算系統(tǒng)開環(huán)傳遞函數(shù)：</p><p>　　被控對象傳遞函數(shù)為：F(s)=1/(s^2+10*s+21)</p><p>　　k1=200時系統(tǒng)開環(huán)傳遞函數(shù)為：G(s)=k1* F(s)*1/s=200/(s*(s^2+10*s+21))</p><p>　　系統(tǒng)開環(huán)傳遞函數(shù)為：G= k1*1/(s*(s^

15、2+10*s+21))</p><p><b>　　結(jié)論圖表數(shù)據(jù)分析</b></p><p>　　(1)利用matlab繪制開環(huán)對數(shù)頻率特性圖：</p><p>　　上圖為幅頻特性曲線圖，下圖為相頻特性曲線圖</p><p>　　同時利用matlab仿真分析閉環(huán)系統(tǒng)的穩(wěn)定裕度：</p><p>　

16、　由圖形可知, 幅值裕度Gm為0.424dB，相角裕度Pm為1.29deg</p><p>　　工程上推薦幅值裕度大于6dB，相角裕度30-60deg.所以該系統(tǒng)穩(wěn)定裕度較小，系統(tǒng)已經(jīng)瀕臨不穩(wěn)定，系統(tǒng)的靈敏性和快速性較強，同時系統(tǒng)的振蕩也較激烈。</p><p>　　(2) 利用開環(huán)傳遞函數(shù)1/(s*(s^2+10*s+21))，用rlocus(sys)語句繪制閉環(huán)系統(tǒng)根軌跡圖，在根軌跡圖

17、上任意選擇軌跡上某點就可以得到對應(yīng)的參數(shù)K1=Gain.</p><p>　　由根軌跡圖我們可以發(fā)現(xiàn)，在虛軸上有3個極點，粗略估計分別當(dāng)K1=209,0,209時取得。說明系統(tǒng)穩(wěn)定的范圍在K1取0—209左右 ，當(dāng)K1大于209后系統(tǒng)不穩(wěn)定。</p><p>　　2．2 系統(tǒng)仿真實驗題目2</p><p>　　司機駕駛汽車的模型如下圖所示，其中，K=5.3。&

18、lt;/p><p>　?。╝）若反應(yīng)時間T=0,求出系統(tǒng)的閉環(huán)頻率響應(yīng)、增益裕度和相角裕度；</p><p>　?。╞）當(dāng)T=0.1s時，估算系統(tǒng)的相角裕度；</p><p>　?。╟）為使系統(tǒng)臨界穩(wěn)定（即相角裕度為00），計算所需的反應(yīng)時間。</p><p>　　2.2.1 實驗步驟</p><p>　?。╝）根據(jù)題

19、目框圖求出系統(tǒng)開環(huán)傳遞函數(shù)；用matlab繪制開環(huán)bode圖求閉環(huán)頻率特性、增益裕度、相角裕度。</p><p>　?。╞）根據(jù)題目框圖求出系統(tǒng)開環(huán)傳遞函數(shù)；并且用matlab繪制出相應(yīng)的Bode圖求出相角裕度。</p><p>　?。╟）利用matlab求出系統(tǒng)相角裕度的具體數(shù)值，令反應(yīng)時間為一系列連續(xù)取值，編寫循環(huán)程序，找到相角裕度為零時的反應(yīng)時間。</p><p&

20、gt;　　2.2.2 實驗分析計算依據(jù)</p><p>　　T=0時系統(tǒng)開環(huán)傳遞函數(shù)：G=(5.3*(s^2+0.8*s+0.32))/(s^3)</p><p>　　T=0時系統(tǒng)閉環(huán)傳遞函數(shù)：G=5.3/(s^3+5.3*s^2+5.3*0.8*s+5.3*0.32)</p><p>　　T=0.1時系統(tǒng)開環(huán)傳遞函數(shù)：G=exp(-0.1*s)*5.3*(s^

21、2+0.8*s+0.32)/s^3</p><p>　　開環(huán)G=exp(-t*s)*5.3*(s^2+0.8*s+0.32)/s^3</p><p>　　2.2.3 結(jié)論圖表數(shù)據(jù)分析</p><p>　　（a）系統(tǒng)的閉環(huán)頻率響應(yīng)、增益裕度和相角裕度</p><p>　　由圖可知幅值裕度Gm為-22.4dB，相角裕度Pm為81.3deg&l

22、t;/p><p><b>　　系統(tǒng)不穩(wěn)定。</b></p><p>　?。╞）T=0.1時估算系統(tǒng)相角裕度</p><p>　　由上圖可知，當(dāng)T=0.1時估算系統(tǒng)相角裕度Pm=51deg</p><p>　?。╟）為使系統(tǒng)臨界穩(wěn)定（即相角裕度為00），計算所需的反應(yīng)時間。</p><p>　　通過前期

23、工作循環(huán)迭代工作，大概確定滿足要求的最優(yōu)反應(yīng)時間在0.1到0.5之間，將T值范圍確定在0.1-0.5步長取0.0001，帶入程序求最優(yōu)值。</p><p>　　由循環(huán)程序求出臨界反應(yīng)時間為T=0.2677</p><p>　　此時系統(tǒng)Pm約等于0，系統(tǒng)處于臨界穩(wěn)定。</p><p>　　2．3 系統(tǒng)仿真實驗題目3</p><p>　　人

24、們常用電子電路和計算機來控制汽車。下圖給出了一個汽車駕駛控制系統(tǒng)的例子，其中，控制桿負(fù)責(zé)操縱車輪。假定司機的反應(yīng)時間T=0.2s。</p><p>　　（a）用Bode圖確定增益K的取值，使閉環(huán)系統(tǒng)的諧振峰Mp不大于2dB； </p><p>　　（b）根據(jù)Mp或相角裕度分別估計系統(tǒng)的阻尼系數(shù)，若所得結(jié)果不同，請解釋原因；</p><p>　?。╟）估計閉環(huán)系統(tǒng)的帶

25、寬。</p><p>　　2.3.1 實驗步驟</p><p>　?。╝）使用循環(huán)求解法。前期試湊出K值的大概范圍，然后在此范圍內(nèi)連續(xù)取一系列的K值，找出符合條件的最優(yōu)K值。求出K的取值范圍。 </p><p>　?。╞）取K=2，求出此時系統(tǒng)的Mp與Pm分別利用兩種公式求出系統(tǒng)的阻尼系數(shù)，比較差異，分析誤差原因。</p><p>　?。?/p>

26、c）利用系統(tǒng)閉環(huán)傳遞函數(shù)，求出系統(tǒng)帶寬?？刂葡到y(tǒng)的閉環(huán)頻率響應(yīng)包括帶寬和諧振峰值，帶寬指閉環(huán)系統(tǒng)的bode圖中幅頻特性曲線下降到－3分貝所對應(yīng)的頻率，諧振峰值指幅頻特性曲線圖中最大幅頻值。bandwidth(sys)可以直接求出帶寬</p><p>　　2.3.2 實驗步驟實驗分析計算依據(jù)</p><p>　　系統(tǒng)閉環(huán)傳遞函數(shù)：G=k*exp(-0.2*s)/(0.1*s^2+s+k*

27、exp(-0.2*s))</p><p>　　系統(tǒng)t=2時閉環(huán)傳遞函數(shù)：G=2*exp(-0.2*s)/(0.1*s^2+s+2*exp(-0.2*s))</p><p>　　系統(tǒng)t=2時開環(huán)傳遞函數(shù)：G=2*exp(-0.2*s)/ (0.1*s^2+s)</p><p>　　2.3.3 結(jié)論圖表數(shù)據(jù)分析</p><p>　?。╝）用B

28、ode圖確定增益K的取值，使閉環(huán)系統(tǒng)的諧振峰Mp不大于2dB；</p><p>　　由循環(huán)程序得最優(yōu)K值為2.4780，此時諧振峰值為2</p><p>　　所以K的取值范圍為（0，2.4780）</p><p>　　（b）根據(jù)Mp或相角裕度分別估計系統(tǒng)的阻尼系數(shù)，若所得結(jié)果不同，請解釋原因；</p><p>　　由程序求出此時系統(tǒng)Mp=0.

29、4837dB，Pm = 3.0576deg</p><p><b>　　阻尼系數(shù)</b></p><p>　　由諧振峰值求得阻尼系數(shù)為：=</p><p>　　由相角裕度求得阻尼系數(shù)為：=</p><p><b>　　誤差原因：</b></p><p>　　（c）估計閉環(huán)系統(tǒng)

30、的帶寬。設(shè)K=2。此時閉環(huán)系統(tǒng)帶寬用bandwidth(sys)語句求出。</p><p>　　同時繪制bode圖，由圖形求帶寬，驗證結(jié)果。</p><p>　　由函數(shù)語句求出系統(tǒng)帶寬為（0，4.2744）rad/s</p><p>　　由閉環(huán)系統(tǒng)bode圖可知，幅值-3dB對應(yīng)的頻率為帶寬頻率4.27rad/s</p><p>　　與函數(shù)求

31、出的帶寬值符合</p><p>　　附錄A 程序源代碼</p><p>　　系統(tǒng)仿真實驗題目1程序源代碼：</p><p>　　1.Bode圖繪制代碼:</p><p>　　>> s=tf('s')</p><p>　　>> G=200/(s*(s^2+10*s+21))<

32、;/p><p>　　>> bode(G)</p><p><b>　　>> grid</b></p><p>　　2.Bode圖穩(wěn)定裕度代碼:</p><p>　　>> s=tf('s')</p><p>　　>> G=200/(s*(s

33、^2+10*s+21))</p><p><b>　　margin(G)</b></p><p>　　3.根軌跡圖源代碼:</p><p>　　>> s=tf('s')</p><p>　　G=1/(s*(s^2+10*s+21))</p><p><b>　

34、　rlocus(G)</b></p><p>　　系統(tǒng)仿真實驗題目2程序源代碼：</p><p>　　T=0求閉環(huán)頻率響應(yīng)、增益裕度和相角裕度：</p><p>　　>> s=tf('s');</p><p>　　>> G=(5.3*(s^2+0.8*s+0.32))/(s^3);</

35、p><p>　　>> margin(G);</p><p>　　T=0.1求相角裕度：</p><p>　　>> s=tf('s');</p><p>　　>> G=exp(-0.1*s)*5.3*(s^2+0.8*s+0.32)/s^3;</p><p>　　>

36、> margin(G);</p><p>　　系統(tǒng)臨界穩(wěn)定，計算所需反應(yīng)時間：</p><p><b>　　clc;</b></p><p><b>　　clear;</b></p><p>　　t=0.1:0.0001:0.5; %t的取值范圍</p><p>　　s

37、=tf('s');</p><p>　　for i=1:4000</p><p>　　G=exp(-t(i)*s)*5.3*(s^2+0.8*s+0.32)/s^3 %傳遞函數(shù)</p><p>　　[Gm,Pm,Wcg,Wcp] = margin(G);</p><p>　　if 0<Pm<0.0001 %相角裕度

38、辨別</p><p><b>　　break</b></p><p><b>　　end</b></p><p>　　T=t(i); %最優(yōu)反應(yīng)時間</p><p><b>　　end</b></p><p>　　margin(G); %畫出bode圖&

39、lt;/p><p>　　系統(tǒng)仿真實驗題目3程序源代碼：</p><p>　　確定增益K的取值，使閉環(huán)系統(tǒng)的諧振峰Mp不大于2dB：</p><p><b>　　clear</b></p><p><b>　　clc</b></p><p>　　k=2:0.0001:3; %k的取

40、值范圍</p><p>　　s=tf('s');</p><p>　　for i=1:10000</p><p>　　G=k(i)*exp(-0.2*s)/(0.1*s^2+s+k(i)*exp(-0.2*s)); %系統(tǒng)閉環(huán)傳遞函數(shù)</p><p>　　[m,p,w,]=bode(G);</p><p&g

41、t;　　mr=max(m);</p><p>　　lwr=20*log10(mr); %將諧振峰值換算單位</p><p>　　if lwr>=2; %諧振峰判斷條件</p><p><b>　　break</b></p><p><b>　　end</b></p><p&

42、gt;　　K=k(i); %最優(yōu)K值</p><p><b>　　end</b></p><p>　　bode(G) %畫出bode圖</p><p>　　2. 求K=2時 閉環(huán)系統(tǒng)諧振峰值</p><p>　　>> s=tf('s');</p><p>　　>

43、> G=2*exp(-0.2*s)/(0.1*s^2+s+2*exp(-0.2*s));</p><p>　　>> [m,p,w,]=bode(G);</p><p>　　mr=max(m);</p><p>　　Mp=20*log10(mr); </p><p><b>　　>> Mp</b&g

44、t;</p><p>　　求K=2時 閉環(huán)系統(tǒng)諧振相角裕度</p><p>　　>> s=tf('s');</p><p>　　>> G=2*exp(-0.2*s)/ (0.1*s^2+s);</p><p>　　>> [Pm] = margin(G);</p><p&g

45、t;<b>　　>> Pm</b></p><p>　　3.求K=2時閉環(huán)系統(tǒng)帶寬：</p><p>　　>> s=tf('s');</p><p>　　G=2*exp(-0.2*s)/(0.1*s^2+s+2*exp(-0.2*s));</p><p>　　bandwidth(G