控制系統(tǒng)性能分析畢業(yè)論文

1、<p><b>　　摘 要</b></p><p>　　隨著計(jì)算機(jī)技術(shù)的發(fā)展，以控制系統(tǒng)為對(duì)象的計(jì)算方法和計(jì)算輔助設(shè)計(jì)問(wèn)題受到普遍的重視，同時(shí)這些理論和方法也在工程實(shí)際中得到了廣泛的應(yīng)用，控制系統(tǒng)也成為自動(dòng)控制技術(shù)和應(yīng)用中不可缺少的重要組成部分。在此當(dāng)中計(jì)算機(jī)軟件起著至關(guān)重要的作用，通過(guò)數(shù)學(xué)軟件可以解決很多有關(guān)問(wèn)題。</p><p>　　本文介紹了控制系統(tǒng)

2、的幾種性能以及對(duì)這幾種性能的分析，并應(yīng)用數(shù)學(xué)軟件Matlab、BASIC程序、以及C語(yǔ)言對(duì)線性系統(tǒng)進(jìn)行程序設(shè)計(jì)，并對(duì)系統(tǒng)性能進(jìn)行詳細(xì)分析，利用MATLAB/SIMULINK對(duì)控制系統(tǒng)進(jìn)行仿真來(lái)對(duì)非線性系統(tǒng)進(jìn)行分析，進(jìn)而全面的對(duì)系統(tǒng)性能進(jìn)行計(jì)算機(jī)的分析。</p><p>　　關(guān)鍵字：控制系統(tǒng)；性能分析；matlab；BASIC；C語(yǔ)言</p><p><b>　　Abstract&

3、lt;/b></p><p>　　along with the computer technology development, take the control system receives the universal value as the object computational method and the computation assistance design question, simu

4、ltaneously these theories and the method also actual obtained the widespread application in the project, the control system have also become in the automatic control technology and the application the essential important

5、 constituent.The computer software is playing the very important role in the middle of this,</p><p>　　This article introduced control system several performance as well as to these performance analysis, and

6、should use mathematics software Matlab, the BASIC procedure, as well as the C language carry on the programming to the linear system, and carries on the multianalysis to the system performance, carries on the simulation

7、using MATLAB/SIMULINK to the control system to come to carry on the analysis to the nonlinear system, then comprehensive carries on the computer to the system performance the </p><p>　　Key words: Control sys

8、tem; Performance analysis; matlab; BASIC; C language</p><p><b>　　目 錄</b></p><p><b>　　1 緒論1</b></p><p>　　1.1 課題選題背景及意義1</p><p>　　1.2 研究現(xiàn)狀及現(xiàn)有成果

9、1</p><p>　　2 課題來(lái)源及需求分析3</p><p>　　2.1 任務(wù)概述及可行性分析3</p><p>　　2.2 實(shí)際要求及需求分析4</p><p>　　3 設(shè)計(jì)方案論證5</p><p>　　3.1 總體方案論證5</p><p>　　3.2 應(yīng)用軟件簡(jiǎn)介5<

10、;/p><p>　　3.2.1 MATLAB語(yǔ)言5</p><p>　　3.2.2 BASIC語(yǔ)言6</p><p>　　3.2.3 C語(yǔ)言7</p><p>　　4 用數(shù)學(xué)軟件對(duì)反饋系統(tǒng)的穩(wěn)定性進(jìn)行程序設(shè)計(jì)8</p><p>　　4.1 用BASIC程序求增益裕量Kg和相角裕量γ的值8</p>

11、<p>　　4.1.1 程序流程圖8</p><p>　　4.1.2 BASIC程序設(shè)計(jì)9</p><p>　　4.1.3 計(jì)算機(jī)運(yùn)行結(jié)果10</p><p>　　4.1.4 結(jié)論10</p><p>　　4.2 用C語(yǔ)言求增益裕量Kg和相角裕量γ的值10</p><p>　　4.2.1 C

12、語(yǔ)言程序設(shè)計(jì)10</p><p>　　4.2.2 計(jì)算機(jī)運(yùn)行結(jié)果及分析12</p><p>　　4.3 具體設(shè)計(jì)分析12</p><p>　　4.3.1 相角裕量和增益裕量的意義12</p><p>　　4.3.2 增益裕量和相角裕量的計(jì)算13</p><p>　　4.3.3 波德圖上求Kg和γ的值14&l

13、t;/p><p>　　4.4 實(shí)際應(yīng)用16</p><p>　　5 用matlab對(duì)線性系統(tǒng)進(jìn)行分析17</p><p>　　5.1 線性系統(tǒng)性質(zhì)17</p><p>　　5.1.1 線性系統(tǒng)及其主要性質(zhì)17</p><p>　　5.1.2 測(cè)試系統(tǒng)的靜態(tài)特性19</p><p>　　5.

14、1.3 測(cè)試系統(tǒng)的動(dòng)態(tài)特性19</p><p>　　5.2 線性系統(tǒng)的時(shí)域分析21</p><p>　　5.2.1 程序設(shè)計(jì)21</p><p>　　5.2.2 一階系統(tǒng)的時(shí)域分析22</p><p>　　5.2.3 二階系統(tǒng)的時(shí)域分析24</p><p>　　5.3 線性系統(tǒng)頻域分析26</p>

15、;<p>　　5.3.1 頻率特性27</p><p>　　5.3.2 頻域響應(yīng)與Nyquist圖程序設(shè)計(jì)27</p><p>　　5.3.3 設(shè)計(jì)分析28</p><p>　　5.4 線性系統(tǒng)的穩(wěn)定性分析29</p><p>　　5.4.1 代數(shù)穩(wěn)定判據(jù)程序設(shè)計(jì)29</p><p>　　5.4

16、.2 Bode圖法判斷系統(tǒng)穩(wěn)定性程序設(shè)計(jì)30</p><p>　　5.4.3 系統(tǒng)穩(wěn)定性概念及定義31</p><p>　　5.4.4 穩(wěn)定條件32</p><p>　　6 MATLAB/SIMULINK仿真在非線性系統(tǒng)中的應(yīng)用33</p><p>　　6.1 程序設(shè)計(jì)33</p><p>　　6.2 非線

17、性系統(tǒng)穩(wěn)定性分析35</p><p>　　6.2.1 典型非線性特性的種類35</p><p>　　6.2.2 非線性系統(tǒng)的若干特性36</p><p>　　6.2.3 非線性系統(tǒng)的穩(wěn)定性分析37</p><p>　　6.3 非線性系統(tǒng)的應(yīng)用及發(fā)展38</p><p><b>　　結(jié) 論40<

18、;/b></p><p><b>　　致 謝41</b></p><p><b>　　參考文獻(xiàn)42</b></p><p><b>　　1 緒論</b></p><p>　　1.1 課題選題背景及意義</p><p>　　在過(guò)去的幾十年中，隨著

19、社會(huì)的不斷發(fā)展，人們的思維方式也在不斷進(jìn)步著。在實(shí)際的工程應(yīng)用中，某些系統(tǒng)數(shù)據(jù)需要進(jìn)行一些大量的計(jì)算、其要求的精度也非常高，可是，單純的筆算和腦算已經(jīng)遠(yuǎn)遠(yuǎn)達(dá)不到人們的要求，這種舊式的方法不但需要耗費(fèi)的大量時(shí)間和勞力，而且計(jì)算出的結(jié)果精度也不高，這是非常落后的，因此逐漸被人們所淘汰。在此基礎(chǔ)上，隨之而來(lái)的便是計(jì)算機(jī)時(shí)代的到來(lái)，在70年代隨著電子技術(shù)的飛速發(fā)展，大量的數(shù)學(xué)軟件及應(yīng)用程序的出現(xiàn)，使得一些復(fù)雜的問(wèn)題變的簡(jiǎn)單易操作，計(jì)算機(jī)程序的開(kāi)

20、發(fā)也給系統(tǒng)工程的計(jì)算和實(shí)施帶來(lái)了極大的方便，利用計(jì)算機(jī)程序?qū)ο到y(tǒng)進(jìn)行程序設(shè)計(jì)及仿真，則大大提高了計(jì)算的精確度并使操作簡(jiǎn)單化，結(jié)果精準(zhǔn)化，這逐漸被科學(xué)家們所利用并應(yīng)用到各個(gè)領(lǐng)域，取得了顯著的效果。同時(shí)，計(jì)算機(jī)程序也給現(xiàn)代社會(huì)的控制系統(tǒng)的設(shè)計(jì)與研究帶來(lái)了很大的方便，使社會(huì)向新時(shí)代邁進(jìn)了一大步。</p><p>　　本次畢業(yè)設(shè)計(jì)著重討論BASIC、C語(yǔ)言程序及MATLAB軟件對(duì)于線性控制系統(tǒng)和非線性控制系統(tǒng)的時(shí)域分析或

21、穩(wěn)定性分析，來(lái)說(shuō)明計(jì)算機(jī)的解算對(duì)自動(dòng)控制理論的發(fā)展起到了不可或缺的作用。數(shù)學(xué)軟件與自動(dòng)控制理論的發(fā)展緊密相關(guān)，是人們進(jìn)行科學(xué)研究及計(jì)算的有力工具，也為自動(dòng)控制理論發(fā)展提供了強(qiáng)大動(dòng)力。</p><p>　　1.2 研究現(xiàn)狀及現(xiàn)有成果</p><p>　　自動(dòng)控制系統(tǒng)是隨著現(xiàn)代大型工業(yè)生產(chǎn)自動(dòng)化的不斷興起而隨之產(chǎn)生的綜合控制系統(tǒng)，在計(jì)算機(jī)控制領(lǐng)域有著相當(dāng)重要的作用。早期的自動(dòng)控制技術(shù)僅僅局限于

22、簡(jiǎn)單的數(shù)字計(jì)算，應(yīng)用的規(guī)模很小，還沒(méi)有應(yīng)用到工程實(shí)踐當(dāng)中，隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展，開(kāi)發(fā)程序的人們已經(jīng)根據(jù)工程的需要，設(shè)計(jì)開(kāi)發(fā)出各種軟件并用相關(guān)設(shè)計(jì)程序?qū)ο到y(tǒng)進(jìn)行計(jì)算研究。對(duì)工業(yè)系統(tǒng)控制起著重要作用。</p><p>　　目前，自動(dòng)控制系統(tǒng)已被廣泛應(yīng)用于人類社會(huì)的各個(gè)領(lǐng)域。在工業(yè)方面，對(duì)于冶金、化工、機(jī)械制造等生產(chǎn)過(guò)程中遇到的各種物理量，包括溫度、流量、壓力、厚度、張力、速度、位置、頻率、相位等，都有相應(yīng)的控制系

23、統(tǒng)。在此基礎(chǔ)上通過(guò)采用數(shù)字計(jì)算機(jī)還建立起了控制性能更好和自動(dòng)化程度更高的數(shù)字控制系統(tǒng)，以及具有控制與管理雙重功能的過(guò)程控制系統(tǒng)。在農(nóng)業(yè)方面的應(yīng)用包括水位自動(dòng)控制系統(tǒng)、農(nóng)業(yè)機(jī)械的自動(dòng)操作系統(tǒng)等。在軍事技術(shù)方面，自動(dòng)控制的應(yīng)用實(shí)例有各種類型的伺服系統(tǒng)、火力控制系統(tǒng)、制導(dǎo)與控制系統(tǒng)等。在航天、航空和航海方面，除了各種形式的控制系統(tǒng)外，應(yīng)用的領(lǐng)域還包括導(dǎo)航系統(tǒng)、遙控系統(tǒng)和各種仿真器。此外，在辦公室自動(dòng)化、圖書管理 、交通 管 理乃至日常家務(wù)方面

24、，自動(dòng)控制技術(shù)也都有著實(shí)際的應(yīng)用。隨著控制理論和控制技術(shù)的發(fā)展，自動(dòng)控制系統(tǒng)的應(yīng)用領(lǐng)域還在不斷擴(kuò)大，幾乎涉及生物、醫(yī)學(xué)、生態(tài)、經(jīng)濟(jì)、社會(huì)等所有領(lǐng)域。隨著人們對(duì)控制系統(tǒng)的重視及自動(dòng)控制系統(tǒng)在現(xiàn)代社會(huì)的發(fā)展，可以預(yù)見(jiàn)，自動(dòng)控制系統(tǒng)系統(tǒng)技術(shù)會(huì)有廣闊的應(yīng)用前景。</p><p>　　2 課題來(lái)源及需求分析</p><p>　　2.1 任務(wù)概述及可行性分析</p><p>　

25、　在當(dāng)今發(fā)展社會(huì)，自動(dòng)控制系統(tǒng)已經(jīng)逐漸在人們的生活中占據(jù)越來(lái)越重要的地位，隨之而來(lái)的便是控制程序的編寫和開(kāi)發(fā)，因此，很多編程軟件也隨之被開(kāi)發(fā)出來(lái)。對(duì)于一個(gè)控制系統(tǒng)，要想實(shí)現(xiàn)其性能最優(yōu)，并在實(shí)際生產(chǎn)中發(fā)揮其最大的作用，就要選擇一個(gè)適合的阮件來(lái)對(duì)系統(tǒng)的性能進(jìn)行徹底全面的分析，通過(guò)分析結(jié)果來(lái)應(yīng)用到實(shí)際系統(tǒng)控制中，以此提高系統(tǒng)的使用價(jià)值。</p><p>　　在過(guò)去的幾十年里， BASIC等數(shù)學(xué)軟件一直應(yīng)用在對(duì)系統(tǒng)進(jìn)行分

26、析工程上，在當(dāng)時(shí)一直是數(shù)值計(jì)算領(lǐng)域所使用的主要語(yǔ)言。BASIC流行了十幾年，幾乎所有的數(shù)值計(jì)算領(lǐng)域里的許多應(yīng)用程序都是用BASIC語(yǔ)言編寫的。隨著現(xiàn)代技術(shù)的發(fā)展，數(shù)學(xué)軟件也逐步在更新?lián)Q代，雖然傳統(tǒng)的數(shù)學(xué)軟件也能進(jìn)行系統(tǒng)相應(yīng)的性能分析，但是，其方法有些落后，已經(jīng)不能滿足自控系統(tǒng)在現(xiàn)代社會(huì)的飛速發(fā)展和應(yīng)用，造成不必要的勞動(dòng)力和時(shí)間上的浪費(fèi)，也容易出現(xiàn)計(jì)算上的錯(cuò)誤。因此，目前應(yīng)用在控制系統(tǒng)領(lǐng)域出現(xiàn)的新類型的軟件，來(lái)解決這類問(wèn)題。matlab語(yǔ)

27、言延續(xù)了過(guò)去軟件的優(yōu)點(diǎn)，并更新加入了現(xiàn)有的新技術(shù)，這為人們進(jìn)行科學(xué)計(jì)算提供了極大的方便，成為自動(dòng)控制系統(tǒng)在計(jì)算機(jī)研究上的不可缺少的一部分。 MATLAB 語(yǔ)言是當(dāng)今自動(dòng)控制領(lǐng)域上最具影響力，也是最有活力的軟件，如今已經(jīng)發(fā)展成一種與控制系統(tǒng)緊密聯(lián)系的計(jì)算機(jī)語(yǔ)言。它具有運(yùn)算功能強(qiáng)、準(zhǔn)確率高、設(shè)計(jì)流程清晰、圖形質(zhì)量高，仿真方便快捷等功能，使之在處理信號(hào)、圖象、控制系統(tǒng)等問(wèn)題時(shí)顯得非常的方便快捷。因此，matlab應(yīng)用的范圍也越來(lái)越廣泛

28、，成為控制領(lǐng)域中不可不用的工具。</p><p>　　控制系統(tǒng)的程序分析只有借助于計(jì)算機(jī)才能完成 ,而MATLAB 正是控制領(lǐng)域進(jìn)行計(jì)算機(jī)輔助設(shè)計(jì)的一種非常好的工具語(yǔ)言。SIMULINK的功能也特別多，操作起來(lái)特別方便 ,所以在近幾年迅速成為國(guó)內(nèi)外系統(tǒng)工程上應(yīng)用最廣泛的軟件之一。</p><p>　　C語(yǔ)言是目前世界上流行、使用最廣泛的計(jì)算機(jī)程序設(shè)計(jì)語(yǔ)言。它既有高級(jí)語(yǔ)言的特點(diǎn)，又具有匯編語(yǔ)

29、言的特點(diǎn)。它可以作為系統(tǒng)設(shè)計(jì)語(yǔ)言，編寫工作系統(tǒng)應(yīng)用程序，也可以作為應(yīng)用程序設(shè)計(jì)語(yǔ)言，編寫不依賴計(jì)算機(jī)硬件的應(yīng)用程序。因此，它的應(yīng)用范圍非常廣泛，利用TC軟件進(jìn)行的C程序設(shè)計(jì)應(yīng)用在實(shí)際的控制系統(tǒng)工程上，在當(dāng)今的社會(huì)生產(chǎn)上有著重要地位。C語(yǔ)言對(duì)操作系統(tǒng)和系統(tǒng)使用程序以及需要對(duì)硬件進(jìn)行操作的場(chǎng)合，用C語(yǔ)言要比其它的計(jì)算型語(yǔ)言要好，很多大型應(yīng)用軟件都是用C語(yǔ)言編寫的。C語(yǔ)言具有繪圖能力強(qiáng)，可移植性，并具備很強(qiáng)的數(shù)據(jù)處理能力，因此適于編寫系統(tǒng)軟件

30、，是數(shù)值計(jì)算的高級(jí)語(yǔ)言。這為人們進(jìn)行科學(xué)計(jì)算提供了極大的方便，成為自動(dòng)控制系統(tǒng)在計(jì)算機(jī)研究上的不可缺少的一部分。</p><p>　　本課題設(shè)計(jì)就是利用這類數(shù)學(xué)軟件來(lái)分別分析系統(tǒng)的性能，為今后的系統(tǒng)在實(shí)際應(yīng)用上的研究奠定基礎(chǔ)。因此，我選擇了BASIC語(yǔ)言、c語(yǔ)言和matlab語(yǔ)言來(lái)對(duì)系統(tǒng)性能進(jìn)行分析。這三種語(yǔ)言對(duì)系統(tǒng)性能的分析在實(shí)際應(yīng)用中是非常必要的。</p><p>　　2.2 實(shí)際要求

31、及需求分析</p><p>　　在控制系統(tǒng)中將系統(tǒng)分為非線性系統(tǒng)和線性系統(tǒng)，系統(tǒng)的分析方法是隨著不同類型的的數(shù)學(xué)模型，在描述其模型的基礎(chǔ)上，用數(shù)學(xué)的方法來(lái)進(jìn)行研究討論的。因此，必須在規(guī)定的工作條件下，利用相應(yīng)的數(shù)學(xué)軟件對(duì)已知系統(tǒng)進(jìn)行系統(tǒng)性能的分析。為了使系統(tǒng)實(shí)現(xiàn)其控制目的，必須具備兩方面性能：</p><p>　?。?）使系統(tǒng)的輸出按輸入信號(hào)要求的輸出值變化；</p><

32、;p>　　（2）是系統(tǒng)的輸出盡量不受外界擾動(dòng)的影響。</p><p>　　要想達(dá)到系統(tǒng)性能的最佳狀態(tài)，則需要用相關(guān)的數(shù)學(xué)軟件對(duì)其進(jìn)行編程仿真，根據(jù)結(jié)果對(duì)系統(tǒng)進(jìn)行各種性能的分析，來(lái)優(yōu)化系統(tǒng)性能，達(dá)到本次課題的目的。同時(shí)可以將理論和實(shí)際聯(lián)系到一起，更好的理解并實(shí)踐應(yīng)用程序，使理論趨于完善，并應(yīng)用到實(shí)際生產(chǎn)生活中，在控制系統(tǒng)的發(fā)展中起到重要作用。</p><p><b>　　3

33、設(shè)計(jì)方案論證</b></p><p>　　3.1 總體方案論證</p><p>　　利用控制系統(tǒng)軟件對(duì)系統(tǒng)進(jìn)行性能分析是自控系統(tǒng)應(yīng)用和實(shí)踐的有效手段和方法。本設(shè)計(jì)對(duì)線性系統(tǒng)和非線性系統(tǒng)可以采用BASIC語(yǔ)言、matlab語(yǔ)言和C語(yǔ)言來(lái)對(duì)其進(jìn)行性能分析的方案；對(duì)線性系統(tǒng)可以采用時(shí)域分析、頻域分析和穩(wěn)定性分析的方案；對(duì)于相角裕量和增益裕量的值可以采用BASIC語(yǔ)言以及C語(yǔ)言對(duì)其進(jìn)行

34、求得，并利用運(yùn)行結(jié)果對(duì)想進(jìn)行性能分析的方案；最后通過(guò)MATLAB/SIMULINK仿真來(lái)對(duì)非線性系統(tǒng)在控制系統(tǒng)中的應(yīng)用進(jìn)行分析與計(jì)算。</p><p>　　本設(shè)計(jì)方案通過(guò)實(shí)際實(shí)驗(yàn)仿真可以證明，利用系統(tǒng)軟件語(yǔ)言可以很方便快捷并且準(zhǔn)確的對(duì)控制系統(tǒng)進(jìn)行性能分析，無(wú)論是系統(tǒng)數(shù)值計(jì)算上還是工程實(shí)際應(yīng)用上，這類數(shù)學(xué)軟件都起著相當(dāng)重要的作用，將系統(tǒng)的分析結(jié)果應(yīng)用在實(shí)際控制系統(tǒng)工程中，對(duì)日后的控制系統(tǒng)的發(fā)展有很大的幫助。因此，我

35、所采用的課題方案可以達(dá)到任務(wù)的要求。但在實(shí)際中，由于BASIC語(yǔ)言出現(xiàn)時(shí)間較早，應(yīng)用起來(lái)不太方便，因此需要將需要新的軟件將其代替。</p><p>　　3.2 應(yīng)用軟件簡(jiǎn)介</p><p>　　3.2.1 MATLAB語(yǔ)言</p><p>　　MATLAB 語(yǔ)言是當(dāng)今自動(dòng)控制領(lǐng)域上最具影響力，也是最有活力的軟件，如今已經(jīng)發(fā)展成一種與控制系統(tǒng)緊密聯(lián)系的計(jì)算機(jī)語(yǔ)言。它

36、具有運(yùn)算功能強(qiáng)、準(zhǔn)確率高、設(shè)計(jì)流程清晰、圖形質(zhì)量高，仿真方便快捷等功能，使之在處理信號(hào)、圖象、控制系統(tǒng)等問(wèn)題時(shí)顯得非常的方便快捷。因此，matlab應(yīng)用的范圍也越來(lái)越廣泛，成為控制領(lǐng)域中不可不用的工具。</p><p>　　Matlab是一種解釋性執(zhí)行語(yǔ)言，具有強(qiáng)大的計(jì)算、仿真、繪圖等功能。由于它使用簡(jiǎn)單，擴(kuò)充方便，尤其是世界上有成千上萬(wàn)的不同領(lǐng)域的科研工作者不停的在自己的科研過(guò)程中擴(kuò)充matlab的功能，使其成

37、為了巨大的知識(shí)寶庫(kù)。可以毫不夸張的說(shuō)，哪怕是你只真正理解了一個(gè)工具箱，那么就是理解了一門非常重要的科學(xué)知識(shí)?？蒲泄ぷ髡咄ǔ？梢酝ㄟ^(guò)matlab來(lái)學(xué)習(xí)某個(gè)領(lǐng)域的科學(xué)知識(shí)，這就是matlab真正在全世界推廣開(kāi)來(lái)的原因。目前的matlab版本已經(jīng)可以方便的設(shè)計(jì)漂亮的界面，它可以象vb等語(yǔ)言一樣設(shè)計(jì)漂亮的用戶接口，同時(shí)因?yàn)橛凶钬S富的函數(shù)庫(kù),所以計(jì)算的功能實(shí)現(xiàn)也很簡(jiǎn)單，進(jìn)一步受到了科研工作者的歡迎。另外，,matlab和其他高級(jí)語(yǔ)言也具有良好的接

38、口，可以方便的實(shí)現(xiàn)與其他語(yǔ)言的混合編程，進(jìn)一步拓寬了matlab的應(yīng)用潛力?？梢哉f(shuō)，matlab已經(jīng)稱為實(shí)際系統(tǒng)工程上一個(gè)不可缺少的計(jì)算機(jī)語(yǔ)言，也使大量的控制系統(tǒng)在實(shí)際生產(chǎn)上得以最大的發(fā)揮，對(duì)以后的研究具有很大的推動(dòng)作用。</p><p>　　3.2.2 BASIC語(yǔ)言</p><p>　　BASIC語(yǔ)言是世界上最早出現(xiàn)并第一個(gè)被正式推廣使用的程序設(shè)計(jì)語(yǔ)言，至今已有三十多年的歷史，始終是

39、數(shù)值計(jì)算領(lǐng)域所使用的主要語(yǔ)言。BASIC是一種高級(jí)語(yǔ)言，它的英文含義是“初學(xué)者通用符號(hào)指令代碼”，是在1965年5月，由美國(guó)科學(xué)家托馬斯·庫(kù)爾茲研制出來(lái)的。10多年后，比爾·蓋茨把它移植到PC上。三十多年來(lái)，BASIC語(yǔ)言一直是初學(xué)計(jì)算機(jī)語(yǔ)言者使用最廣泛的一種高級(jí)語(yǔ)言。它能進(jìn)行數(shù)值計(jì)算、畫圖、演奏音樂(lè)，功能十分強(qiáng)大，而學(xué)起來(lái)又是非常容易。 </p><p>　　BASIC語(yǔ)言的主要特點(diǎn)是：&l

40、t;/p><p>　?。?）構(gòu)成簡(jiǎn)單。 BASIC語(yǔ)言的最基本語(yǔ)句只有17種，而且它們都是常見(jiàn)的英文單詞或其變形，很容易學(xué)習(xí)和掌握。 </p><p>　　（2）是一種“人機(jī)會(huì)話”式的語(yǔ)言。通過(guò)鍵盤操作，用BASIC語(yǔ)言編寫完的程序，可以在計(jì)算機(jī)上邊編寫、邊修改、邊運(yùn)行。而且還可以在運(yùn)行中向人們提示信息的指出錯(cuò)誤，要求人去改正，即實(shí)現(xiàn)了人和機(jī)器的對(duì)話。 </p><p>

41、;　?。?）BASIC語(yǔ)言應(yīng)用廣泛。在許多的控制系統(tǒng)的程序編譯仿真上都起著很重要的作用。</p><p>　　BASIC是為科學(xué)、工程問(wèn)題中的那些能用數(shù)學(xué)公式表達(dá)的問(wèn)題而設(shè)計(jì)的，其數(shù)值計(jì)算的功能較強(qiáng)。在對(duì)矩陣和復(fù)數(shù)進(jìn)行運(yùn)算上也跟matlab很像。BASIC語(yǔ)言自誕生以來(lái)廣泛地應(yīng)用于數(shù)值計(jì)算領(lǐng)域，并以其特有的功能在數(shù)值、科學(xué)和工程計(jì)算領(lǐng)域發(fā)揮著重要作用。</p><p>　　3.2.3 C

42、語(yǔ)言</p><p>　　C語(yǔ)言是國(guó)際上廣泛流行的計(jì)算機(jī)高級(jí)語(yǔ)言。早期的C語(yǔ)言主要是用于UNIX系統(tǒng)。由于C語(yǔ)言的強(qiáng)大功能和各方面的優(yōu)點(diǎn)逐漸為人們認(rèn)識(shí)，到了八十年代，C開(kāi)始進(jìn)入其它操作系統(tǒng)，并很快在各類大、中、小和微型計(jì)算機(jī)上得到了廣泛的使用。成為當(dāng)代最優(yōu)秀的程序設(shè)計(jì)語(yǔ)言之一。 </p><p>　　C語(yǔ)言是一種結(jié)構(gòu)化語(yǔ)言。它層次清晰，便于按模塊化方式組織程序，易于調(diào)試和維護(hù)。C語(yǔ)言的表現(xiàn)

43、能力和處理能力極強(qiáng)。它不僅具有豐富的運(yùn)算符和數(shù)據(jù)類型，便于實(shí)現(xiàn)各類復(fù)雜的數(shù)據(jù)結(jié)構(gòu)。它還可以直接訪問(wèn)內(nèi)存的物理地址，進(jìn)行位(bit)一級(jí)的操作。由于C語(yǔ)言實(shí)現(xiàn)了對(duì)硬件的編程操作，因此C語(yǔ)言集高級(jí)語(yǔ)言和低級(jí)語(yǔ)言的功能于一體。既可用于系統(tǒng)軟件的開(kāi)發(fā)，也適合于應(yīng)用軟件的開(kāi)發(fā)。此外，C語(yǔ)言還具有效率高，可移植性強(qiáng)等特點(diǎn)。因此廣泛地移植到了各類各型計(jì)算機(jī)上，從而形成了多種版本的C語(yǔ)言。C語(yǔ)言發(fā)展迅速, 而且成為最受歡迎的語(yǔ)言之一。用C語(yǔ)言加上一些匯

44、編語(yǔ)言子程序, 就更能顯示C語(yǔ)言的優(yōu)勢(shì)了,象PC- DOS 、WORDSTAR等就是用這種方法編寫的。目前，用C語(yǔ)言來(lái)編寫控制系統(tǒng)的應(yīng)用程序是現(xiàn)代化控制系統(tǒng)的軟件設(shè)計(jì)實(shí)施中也處于重要地位，使我們的生活向現(xiàn)代化更進(jìn)了一步。</p><p>　　4 用數(shù)學(xué)軟件對(duì)反饋系統(tǒng)的穩(wěn)定性進(jìn)行程序設(shè)計(jì)</p><p>　　4.1 用BASIC程序求增益裕量Kg和相角裕量γ的值</p><

45、;p>　　單位反饋系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為G（S）=</p><p>　　4.1.1 程序流程圖</p><p>　　圖4-1 程序求Kg和γ值的流程圖</p><p>　　流程圖圖4-1的編碼符號(hào)為</p><p>　　W=ω，G2=,P=G(S)的相角，PM=γ相角裕量,GOSUB計(jì)算G2和P。</p><p>　

46、　4.1.2 BASIC程序設(shè)計(jì)</p><p><b>　　1.REM</b></p><p>　　10.LET W=0.01</p><p>　　20.GOSUB 170</p><p>　　30.IF G2<1 THEN 90</p><p>　　40.IF G2<100 TH

47、EN 70</p><p>　　50.LET W=2*W</p><p>　　60.GOTO 20</p><p>　　70.LET W=1.01*W</p><p>　　80 GOTO 20</p><p>　　90 IF P≧180° THEN 210 </p><p>　　100

48、 PRINT “UNITY GAIN(峰位增益)”“W=”W,“P=”P</p><p>　　110 LET W=1.01*W</p><p>　　120 GOSUB 170</p><p>　　130 IF P≧180° THEN 150</p><p>　　140 GOTO 110</p><p>　　1

49、50 PRINT “W=”W,“GAIN MARGIN(增益裕量)=”4.3429448*LOG(1/G2)</p><p>　　160 GOTO 990</p><p>　　170 LET P=57.29578*（2*ATN（0.01*W）+ATN（0.2*W）+ATN（1.5*W）-ATN（0.1*W））</p><p>　　18O LET X=W*W</

50、p><p>　　190 LET G2=99*99*(1+0.01*X)/((1+0.0001*X)^ 2*(1+0.04*X)*(1+2.25*X))</p><p>　　200 RETURN</p><p>　　210 PRINT “W=”W,“SYSTEM UNSTABLE(系統(tǒng)不穩(wěn)定)”</p><p><b>　　990 END

51、</b></p><p>　　4.1.3 計(jì)算機(jī)運(yùn)行結(jié)果</p><p><b>　　RUN</b></p><p>　　UNITY GAIN(單位增益)W=31.21</p><p>　　P=-132.107</p><p>　　W=96.0747 GAIN MARGIN(增益裕量)

52、=14.9268</p><p><b>　　4.1.4 結(jié)論</b></p><p>　　根據(jù)basic程序數(shù)據(jù)，求得，</p><p>　　γ=180+G(j)=180°-132.107°</p><p><b>　　=47.893°</b></p>

53、<p>　　Kg=14.9268（db）</p><p>　　因此，系統(tǒng)是穩(wěn)定的，而且儲(chǔ)備了適當(dāng)?shù)姆€(wěn)定裕度。</p><p>　　4.2 用C語(yǔ)言求增益裕量Kg和相角裕量γ的值</p><p>　　單位反饋系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為G（S）=</p><p>　　4.2.1 C語(yǔ)言程序設(shè)計(jì) </p><p>　

54、　#include <stdio.h></p><p>　　#include <math.h></p><p>　　int main()</p><p><b>　　{</b></p><p>　　double w = 0.01;</p><p>　　A:double p

55、 = </p><p>　　57.29578*(2*atan(0.01*w)+atan(0.2*w)+atan(1.5*w)-atan(0.1*w));</p><p>　　double x = w*w;</p><p>　　double g2 = </p><p>　　99*99*(1+0.01*x)/(pow(1+0.0001*x,2)

56、*(1+0.04*x)*(1+2.25*x));</p><p>　　if (g2 < 1)</p><p><b>　　{</b></p><p>　　if (p >= 180)</p><p><b>　　{</b></p><p>　　printf(&quo

57、t;w= %lf\n",w,"system unstable(系統(tǒng)不穩(wěn)定)");</p><p>　　return -1;</p><p><b>　　}</b></p><p>　　printf("unity gain (峰位增益)""w= %lf, p= %lf\n",w

58、,p);</p><p><b>　　while (1)</b></p><p><b>　　{</b></p><p>　　w = 1.01*w;</p><p><b>　　p = </b></p><p>　　57.29578*(2*atan

59、(0.01*w)+atan(0.2*w)+atan(1.5*w)-atan(0.1*w));</p><p><b>　　x = w*w;</b></p><p><b>　　g2 = </b></p><p>　　99*99*(1+0.01*x)/(pow(1+0.0001*x,2)*(1+0.04*x)*(1+2.25

60、*x));</p><p>　　if (p >= 180) </p><p><b>　　{</b></p><p>　　printf("w= %lf GAIN MARGIN(增益裕量)= </p><p>　　%lf\n",w,4.3429448*log(1/g2));</p>

61、<p><b>　　return 0;</b></p><p><b>　　}</b></p><p><b>　　}</b></p><p><b>　　}</b></p><p>　　else if ((g2 < 100) &&

62、amp; (g2 >= 1))</p><p><b>　　{</b></p><p>　　w = 1.01*w;</p><p><b>　　goto A;</b></p><p><b>　　}</b></p><p><b>　　e

63、lse</b></p><p><b>　　{</b></p><p><b>　　w = 2*w;</b></p><p><b>　　goto A;</b></p><p><b>　　}</b></p><p>&

64、lt;b>　　return 0;</b></p><p><b>　　}</b></p><p>　　4.2.2 計(jì)算機(jī)運(yùn)行結(jié)果及分析</p><p><b>　　程序運(yùn)行結(jié)果為：</b></p><p>　　Unity gain(峰位增益)W=31.210004,P=132.10

65、6852</p><p>　　W=96.074705 ,gain margin（增益裕量）=14.926804</p><p>　　結(jié)果表明γ=47.893°，Kg=14.9268（db）。所以，系統(tǒng)是穩(wěn)定的，儲(chǔ)備裕量充分。</p><p>　　4.3 具體設(shè)計(jì)分析</p><p>　　4.3.1 相角裕量和增益裕量的意義</

66、p><p>　　對(duì)于控制系統(tǒng)，可以根據(jù)系統(tǒng)開(kāi)環(huán)幅相特性曲線與（—1，j0）點(diǎn)的關(guān)系，來(lái)判定閉環(huán)系統(tǒng)是否穩(wěn)定，即奈奎斯特判據(jù)，其判別式Z=P-2N，只有閉環(huán)系統(tǒng)在右半S平面上的極點(diǎn)數(shù)Z=0，閉環(huán)系統(tǒng)才是穩(wěn)定的。</p><p>　　在實(shí)際的控制系統(tǒng)工程中，穩(wěn)定性是非常重要的。圖4-2中三條幅相特性曲線，都不包圍（—1，j0）點(diǎn)，所以閉環(huán)系統(tǒng)都是穩(wěn)定的。曲線3離（—1，j0）點(diǎn)很近，系統(tǒng)受到外界擾

67、動(dòng)后，其穩(wěn)定性容易遭到破壞。曲線1的穩(wěn)定性很好，但快速性很差，所以，系統(tǒng)過(guò)于穩(wěn)定也是不好的。為了使控制系統(tǒng)的穩(wěn)定性和快速性都保持良好的狀態(tài)，因此增益裕量和相角裕量的研究十分必要。</p><p>　　圖4-2 三條幅相特性曲線</p><p>　　4.3.2 增益裕量和相角裕量的計(jì)算</p><p>　　圖4-3相位裕量和增益裕量</p><p&

68、gt;　　增益裕量如圖4-3所示，幅相特性曲線與負(fù)實(shí)軸的交點(diǎn)頻率Wg為相位穿越頻率， 原點(diǎn)與wg之間的距離為G（jwg）H(jwg)，則距離的倒數(shù)Kg=為系統(tǒng)的增益裕量。從圖中可以看出， Kg>1時(shí)系統(tǒng)穩(wěn)定，且Kg越大系統(tǒng)越穩(wěn)定， Kg<1時(shí)，系統(tǒng)將不穩(wěn)定，如果系統(tǒng)的開(kāi)環(huán)傳遞系數(shù)增大Kg倍時(shí)，系統(tǒng)為臨界穩(wěn)定。</p><p>　　相角裕量如圖4-3所示，幅相特性曲線與單位圓的交點(diǎn)Wc為稱幅值穿越頻率，

69、則直線OWc與負(fù)實(shí)軸的夾角γ=+| G(jwc)H(jwc)|為相角裕量。從圖中可以看出，γ越大，系統(tǒng)越穩(wěn)定。如果系統(tǒng)推遲γ角，則系統(tǒng)為臨界穩(wěn)定。</p><p><b>　　其中：</b></p><p>　　（1）增益裕量和相角裕量，通常只用來(lái)分析最小相位系統(tǒng)；</p><p>　?。?）增益裕量和相角裕量，同時(shí)給出才能說(shuō)明系統(tǒng)的穩(wěn)定程度；

70、</p><p>　?。?）選擇多大的增益裕量和相角裕量，目前尚無(wú)統(tǒng)一標(biāo)準(zhǔn)，一般而言Kg>6(db),選在Kg=10~20(db)；γ>,而選在γ=，是適宜的；</p><p>　　（4）在實(shí)際工程中，應(yīng)避免使用條件穩(wěn)定系統(tǒng)。</p><p>　　4.3.3 波德圖上求Kg和γ的值</p><p>　　通過(guò)幅相特性曲線可以求出增益

71、裕量和相角裕量，但是在實(shí)際上卻需要計(jì)算幅相曲線許多點(diǎn)的值，才能作出幅相特性曲線，所以這種方法在實(shí)際應(yīng)用上是有困難的，而工程上經(jīng)常是利用對(duì)數(shù)頻率特性來(lái)求Kg和γ的值。</p><p>　　圖4-4伯德圖上Kg和γ</p><p>　　由圖可知，極坐標(biāo)圖和伯德圖間存在如下對(duì)應(yīng)關(guān)系：</p><p>　　第一、伯德圖上的O（db）線對(duì)應(yīng)極坐標(biāo)圖上的單位圓；</p&g

72、t;<p>　　第二、伯德圖上的-線對(duì)應(yīng)極坐標(biāo)圖上的負(fù)實(shí)軸；</p><p>　　第三、伯德圖O（db）以上區(qū)域?qū)?yīng)極坐標(biāo)圖上的單位圓外部分，伯德圖O（db）以下區(qū)域?qū)?yīng)極坐標(biāo)圖上的單位圓內(nèi)部分；</p><p>　　第四、伯德圖上隨ω增加，φ（ω）從上向下穿越-180線，對(duì)應(yīng)幅相特性曲線逆時(shí)針包圍（-1，j0）點(diǎn)；</p><p>　　第五、伯德圖與

73、幅相特性曲線的穩(wěn)定判據(jù)，都是</p><p><b>　　Z=P-2N</b></p><p><b>　　N=—</b></p><p>　　顯然，Kg和γ值由圖4-4可得。</p><p>　　綜上所述，奈奎斯特穩(wěn)定判據(jù)表現(xiàn)為兩種形式，一種是幅相特性曲線形式，另一種是伯德圖形式，它們即能判定系統(tǒng)

74、的穩(wěn)定性，又能確定系統(tǒng)的相角裕量和增益裕量。工程上通常采用后一種形式確定Kg和γ值。</p><p><b>　　4.4 實(shí)際應(yīng)用</b></p><p>　　在控制系統(tǒng)的實(shí)際應(yīng)用中，增益裕量和相角裕量的研究是非常有必要的，可以應(yīng)用在判斷系統(tǒng)的穩(wěn)定性，也能更好更有效的提高系統(tǒng)的性能，在實(shí)際的工程操作中處于非常重要的地位，并在當(dāng)今人們的生產(chǎn)生活中起著非常重要的角色作用。

75、因此，自動(dòng)控制系統(tǒng)的分析和計(jì)算，摻進(jìn)計(jì)算機(jī)的幫助，前景將是美好的。</p><p>　　5 用matlab對(duì)線性系統(tǒng)進(jìn)行分析</p><p>　　5.1 線性系統(tǒng)性質(zhì)</p><p>　　5.1.1 線性系統(tǒng)及其主要性質(zhì) </p><p>　　我們把外界對(duì)系統(tǒng)的作用稱之為系統(tǒng)的輸入或激勵(lì)，而將系統(tǒng)對(duì)輸入的反應(yīng)稱為系統(tǒng)的輸出或響應(yīng)。</p

76、><p>　　圖5-1 測(cè)試系統(tǒng)方框圖</p><p>　　當(dāng)系統(tǒng)的輸入x(t)和輸出y(t)之間的關(guān)系可用常系數(shù)線性微分方程式（5.1）來(lái)描述時(shí)，則稱該系統(tǒng)為定常線性系統(tǒng)或時(shí)不變線性系統(tǒng)。 </p><p>　　= 式（5.1）</p><p>　　式中：t為時(shí)間自變量；系數(shù),，…，和，，…, 均為不

77、隨時(shí)間而變化的常數(shù)。 </p><p>　　若以x(t)y(t)表示定常線性系統(tǒng)輸入與輸出的對(duì)應(yīng)關(guān)系，則定常線性系統(tǒng)具有以下主要性質(zhì)： </p><p><b>　　(1)疊加原理 </b></p><p>　　當(dāng)幾個(gè)輸入同時(shí)作用于線性系統(tǒng)時(shí)，則其響應(yīng)等于各個(gè)輸入單獨(dú)作用于該系統(tǒng)的響應(yīng)之和。即 </p><p><

78、b>　　若 ;</b></p><p><b>　　則 </b></p><p>　　疊加原理表明，對(duì)于線性系統(tǒng)，一個(gè)輸入的存在并不影響另一個(gè)輸入的響應(yīng)，各個(gè)輸入產(chǎn)生的響應(yīng)是互不影響的。因此，對(duì)于一個(gè)復(fù)雜的輸入，就可以將其分解成一系列簡(jiǎn)單的輸入之和，系統(tǒng)對(duì)復(fù)雜激勵(lì)的響應(yīng)便等于這些簡(jiǎn)單輸入的響應(yīng)之和。</p><p><b

79、>　　(2)比例特性</b></p><p>　　若線性系統(tǒng)的輸入擴(kuò)大 k 倍，則其響應(yīng)也將擴(kuò)大k倍，即對(duì)于任意常數(shù)k，必有kx(t)ky(t)</p><p><b>　　(3)微分特性</b></p><p>　　線性系統(tǒng)對(duì)輸入導(dǎo)數(shù)的響應(yīng)等于對(duì)該輸入響應(yīng)的導(dǎo)數(shù)，即 </p><p><b

80、>　　(4)積分特性 </b></p><p>　　若線性系統(tǒng)的初始狀態(tài)為零（即當(dāng)輸入為零時(shí)，其響應(yīng)也為零）。則對(duì)輸入積分的響應(yīng)等于對(duì)該輸入響應(yīng)的積分，即 </p><p><b>　　(5)頻率保持特性</b></p><p>　　若線性系統(tǒng)的輸入為某一頻率的簡(jiǎn)諧信號(hào)，則其穩(wěn)態(tài)響應(yīng)必是同一頻率的簡(jiǎn)諧信號(hào)。&l

81、t;/p><p>　　線性系統(tǒng)的頻率保持性，在測(cè)試工作中具有非常重要的作用。因?yàn)樵趯?shí)際測(cè)試中，測(cè)試得到的信號(hào)常常會(huì)受到其他信號(hào)或噪聲的干擾，這時(shí)依據(jù)頻率保持特性可以認(rèn)定測(cè)得信號(hào)中只有與輸入信號(hào)相同的頻率成分才是真正由輸入引起的輸出。同樣，在故障診斷中，根據(jù)測(cè)試信號(hào)的主要頻率成分，在排除干擾的基礎(chǔ)上，依據(jù)頻率保持特性推出輸入信號(hào)也應(yīng)包含該頻率成分，通過(guò)尋找產(chǎn)生該頻率成分的原因，就可以診斷出故障的原因。</p>

82、;<p>　　5.1.2 測(cè)試系統(tǒng)的靜態(tài)特性 </p><p>　　在式（5.1）描述的線性系統(tǒng)中，當(dāng)系統(tǒng)的輸入為常數(shù)，即輸入信號(hào)的幅值不隨時(shí)間變化或其隨時(shí)間變化的周期遠(yuǎn)遠(yuǎn)大于測(cè)試時(shí)間時(shí)，式（5.1）變成： </p><p><b>　　式(5.2)</b></p><p>　　也就是說(shuō)，理想線性系統(tǒng)其輸出與輸入之間是呈單調(diào)、線性

83、比例的關(guān)系，即輸入、輸出關(guān)系是一條理想的直線，斜率為常數(shù)。 </p><p>　　但是實(shí)際測(cè)試系統(tǒng)并非是理想定常線性系統(tǒng)，輸入、輸出曲線并不是理想的直線，式（5.2）實(shí)際上變成 </p><p>　　測(cè)試系統(tǒng)的靜態(tài)特性就是在靜態(tài)測(cè)量情況下描述實(shí)際測(cè)試裝置與理想定常線性系統(tǒng)的接近程度。 </p><p>　　5.1.3 測(cè)試系統(tǒng)的動(dòng)態(tài)特性 </p>&l

84、t;p>　　測(cè)試系統(tǒng)的動(dòng)態(tài)特性是指，輸入量隨時(shí)間變化時(shí)其輸出隨輸入而變化的關(guān)系。一般在所考慮的測(cè)量范圍內(nèi)，測(cè)試系統(tǒng)都可以認(rèn)為是線性系統(tǒng)，因此就可以用式（5.1）這一定常線性系統(tǒng)微分方程來(lái)描述測(cè)試系統(tǒng)以及和輸入x(t)、輸出y(t)之間的關(guān)系，通過(guò)拉普拉斯變換建立其相應(yīng)的“傳遞函數(shù)”，該傳遞函數(shù)就能描述測(cè)試裝置的固有動(dòng)態(tài)特性，通過(guò)傅里葉變換建立其相應(yīng)的“頻率響應(yīng)函數(shù)”，以此來(lái)描述測(cè)試系統(tǒng)的特性。 </p><p&

85、gt;<b>　　[1] 傳遞函數(shù) </b></p><p>　　系統(tǒng)的傳遞函數(shù)H(s)為輸出量和輸入量的拉普拉斯變換之比，即 </p><p>　　其中s是復(fù)變量，即。</p><p>　　傳遞函數(shù)包含了瞬態(tài)、穩(wěn)態(tài)時(shí)間響應(yīng)和頻率響應(yīng)的全部信息。有以下幾個(gè)特點(diǎn)： </p><p>　?。?）H(s)描述了系統(tǒng)本身的動(dòng)態(tài)特

86、性，而與輸入量x(t)及系統(tǒng)的初始狀態(tài)無(wú)關(guān)。 </p><p>　?。?）H(s)是對(duì)物理系統(tǒng)特性的一種數(shù)學(xué)描述，而與系統(tǒng)的具體物理結(jié)構(gòu)無(wú)關(guān)。H(s)是通過(guò)對(duì)實(shí)際的物理系統(tǒng)抽象成數(shù)學(xué)模型式（5.1）后，經(jīng)過(guò)拉普拉斯變換后所得出的，所以同一傳遞函數(shù)可以表征具有相同傳輸特性的不同物理系統(tǒng)。 </p><p>　　（3）H(s)中的分母取決于系統(tǒng)的結(jié)構(gòu)，而分子則表示系統(tǒng)同外界之間的聯(lián)系，如輸入點(diǎn)

87、的位置、輸入方式、被測(cè)量以及測(cè)點(diǎn)布置情況等。分母中s的冪次n代表系統(tǒng)微分方程的階數(shù)，當(dāng)n=1時(shí)，系統(tǒng)為一階系統(tǒng)；n=2 時(shí)，系統(tǒng)為二階系統(tǒng)。 </p><p>　　一般測(cè)試系統(tǒng)都是穩(wěn)定系統(tǒng)，其分母中s的冪次總是高于分子中s的冪次（n>m ）。 </p><p>　　[2]頻率響應(yīng)函數(shù) </p><p>　　傳遞函數(shù)與頻率響應(yīng)函數(shù)的關(guān)系為： </p>

88、<p>　　對(duì)于定常線性系統(tǒng)，有頻率響應(yīng)函數(shù)H（w）：</p><p><b>　　式中。</b></p><p>　　若在t=0時(shí)刻將輸入信號(hào)接入定常線性系統(tǒng)時(shí)，將代入拉普拉斯變換中，實(shí)際上是將拉普拉斯變換變成傅里葉變換。由于系統(tǒng)的初始條件為零，因此系統(tǒng)的頻率響應(yīng)函數(shù)H(w)就成為輸出y(t) 、輸入x(t)的傅里葉變換Y(ω)、X(ω)之比，即

89、 </p><p>　　因此，在測(cè)得輸出y(t)和輸入x(t)后，由其傅里葉變換Y(ω)和X(ω)可求得頻率響應(yīng)函數(shù)H(w)= Y(ω)/X(ω)。</p><p>　　頻率響應(yīng)函數(shù)是描述系統(tǒng)的簡(jiǎn)諧輸入和其穩(wěn)態(tài)輸出的關(guān)系，在測(cè)量系統(tǒng)頻率響應(yīng)函數(shù)時(shí)，必須在系統(tǒng)響應(yīng)達(dá)到穩(wěn)態(tài)階段時(shí)才測(cè)量。 </p><p>　　頻率響應(yīng)函數(shù)是

90、復(fù)數(shù)，因此可寫為</p><p>　　式中：－系統(tǒng)的幅頻特性；－系統(tǒng)的相頻特性。 </p><p>　　因此，系統(tǒng)的頻率響應(yīng)函數(shù)或其幅頻特性 、相頻特性，都是簡(jiǎn)諧輸入頻率ω的函數(shù)。 </p><p>　　[3] 脈沖響應(yīng)函數(shù) </p><p>　　若輸入為單位脈沖，即時(shí)，則x(s)=1。因此，有 H(s)=Y(s)。經(jīng)拉普拉斯反變換為 y(t

91、)=h(t)。h(t)為系統(tǒng)的脈沖響應(yīng)函數(shù)。脈沖響應(yīng)函數(shù)可作為系統(tǒng)特性的時(shí)域描述。 </p><p>　　由上可以得出，系統(tǒng)特性在時(shí)域可以用h(t)來(lái)描述，在頻域可以用來(lái)描述，在復(fù)數(shù)域可以用H(s)來(lái)描述。三者的關(guān)系也是一一對(duì)應(yīng)的。</p><p>　　5.2 線性系統(tǒng)的時(shí)域分析</p><p>　　在確定系統(tǒng)的數(shù)學(xué)模型后，便可以用幾種不同的方法去分析控制系統(tǒng)的性能

92、。在經(jīng)典控制理論中，常用時(shí)域法、根軌跡法和頻域法來(lái)分析線形系統(tǒng)的性能。時(shí)域分析法是一種直接在時(shí)間域中對(duì)系統(tǒng)進(jìn)行分析的方法，具有直觀、準(zhǔn)確的優(yōu)點(diǎn)，并且可以提供系統(tǒng)時(shí)間響應(yīng)的全部信息。</p><p>　　5.2.1 程序設(shè)計(jì)</p><p>　　典型二階系統(tǒng)，繪制=0.5，分別為2、4、6、8、10、12時(shí)的系統(tǒng)單位階躍響應(yīng)。</p><p>　　編寫M文件：%求系統(tǒng)

93、的單位階躍響應(yīng)</p><p><b>　　w=2:2:12;</b></p><p>　　kosai=0.5;</p><p><b>　　figure(1)</b></p><p><b>　　hold on</b></p><p><b>

94、;　　for Wn=w</b></p><p><b>　　num=Wn^2;</b></p><p>　　den=[1,2*kosai*Wn,Wn^2];</p><p>　　step(num,den);</p><p><b>　　end</b></p><p&g

95、t;<b>　　hold off</b></p><p><b>　　grid on;</b></p><p>　　title(‘單位階躍響應(yīng)’)</p><p>　　xlabel(‘時(shí)間’)</p><p>　　ylabel(‘振幅’)</p><p><b>　

96、　得到運(yùn)行結(jié)果：</b></p><p>　　從圖可以看出，隨著的逐漸增大，系統(tǒng)的響應(yīng)速度明顯加快。</p><p>　　5.2.2 一階系統(tǒng)的時(shí)域分析</p><p>　　一階系統(tǒng)，即可以用一階微分方程來(lái)描述的系統(tǒng)。常常被應(yīng)用在控制工程實(shí)踐中，有些高階系統(tǒng)的特性也可用一階系統(tǒng)的特性近似表征。</p><p>　　一階系統(tǒng)傳遞函數(shù)

97、為 </p><p>　　式中T為慣性環(huán)節(jié)實(shí)踐常數(shù),反應(yīng)系統(tǒng)慣性。</p><p>　　一階系統(tǒng)的方框圖如圖所示</p><p><b>　　圖5-2 一階系統(tǒng)</b></p><p>　　圖5-3 單位反饋一階系統(tǒng)</p><p>　　當(dāng)輸入信號(hào)r(t)=1(t)時(shí)，系統(tǒng)單位

98、階躍響應(yīng)c(t)=1- (t≥0)</p><p>　　圖5-4 系統(tǒng)單位階躍響應(yīng)曲線</p><p>　　時(shí)間常數(shù)T反應(yīng)系統(tǒng)的慣性，慣性越小，系統(tǒng)響應(yīng)過(guò)程越快；反之T越大，響應(yīng)越慢。</p><p>　　當(dāng)輸入信號(hào)r(t)=(t)時(shí)，系統(tǒng)單位脈沖響應(yīng)c(t)= (t≥0)</p><p>　　圖5-5 系統(tǒng)單位脈沖響應(yīng)曲線&l

99、t;/p><p>　　當(dāng)t=0時(shí)，響應(yīng)取最大值c(t)=1/T；當(dāng)t時(shí)，響應(yīng)幅值衰減為零。</p><p>　　在單位斜坡輸入作用下，輸入信號(hào)r(t)=t1(t)時(shí)，系統(tǒng)的單位斜坡響應(yīng)c(t)=t-T+T (t≥0)</p><p>　　圖5-6系統(tǒng)的單位斜坡響應(yīng)曲線</p><p>　　脈沖函數(shù)(t)和斜坡函數(shù)t1(t)分別是階躍函數(shù)1(

100、t)的對(duì)時(shí)間t的一階微分和積分，而系統(tǒng)的單位脈沖響應(yīng)和單位斜坡響應(yīng)分別是系統(tǒng)的單位階躍響應(yīng)對(duì)時(shí)間t的一階微分和積分。這一關(guān)系表明，系統(tǒng)對(duì)輸入信號(hào)導(dǎo)數(shù)的響應(yīng)等于系統(tǒng)對(duì)該輸入信號(hào)響應(yīng)的導(dǎo)數(shù)；系統(tǒng)對(duì)輸入信號(hào)積分的響應(yīng)等于系統(tǒng)對(duì)該輸入信號(hào)響應(yīng)的積分。</p><p>　　5.2.3 二階系統(tǒng)的時(shí)域分析</p><p>　　二階系統(tǒng)通常表示為圖示的單位負(fù)反饋的結(jié)構(gòu)形式。</p><

101、;p>　　圖5-7 標(biāo)準(zhǔn)化二階系統(tǒng)</p><p>　　二階系統(tǒng)閉環(huán)傳遞函數(shù)為 </p><p>　　其中為阻尼比，無(wú)阻尼振蕩頻率=1/T </p><p>　　則 </p><p>　　在單位階躍函數(shù)作用下，二階系統(tǒng)的單位階躍響應(yīng)有四種形式:</p><p>

102、;　　(1).無(wú)阻尼響應(yīng)：c(t)=1-cost </p><p>　　(2).欠阻尼響應(yīng)：c(t)= (其中系統(tǒng)阻尼振蕩頻率，為時(shí)間響應(yīng)的初始相位角)</p><p>　　(3).臨界阻尼響應(yīng)：c(t)=1-</p><p>　　(4).過(guò)阻尼響應(yīng)：c(t)=</p><p>　　二階系統(tǒng)單位脈沖響應(yīng)g(t)形式有四種：</p&

103、gt;<p>　　(1). 無(wú)阻尼響應(yīng)（=0）：g(t)=</p><p>　　(2). 欠阻尼響應(yīng)（1>>0）：g(t)=</p><p>　　(3). 臨界阻尼響應(yīng)(=1)：g(t)=</p><p>　　(4). 過(guò)阻尼響應(yīng)(>1)：g(t)=</p><p>　　圖5-8 二階系統(tǒng)的脈沖響應(yīng)曲線<

104、/p><p>　　由r(t)=t1(t)，二階系統(tǒng)的單位斜坡響應(yīng)的時(shí)間響應(yīng)為</p><p>　　其中為阻尼比，為自然頻率，為阻尼振蕩頻率。</p><p>　　圖5-9 二階系統(tǒng)的單位斜坡響應(yīng)曲線</p><p>　　5.3 線性系統(tǒng)頻域分析</p><p>　　頻域分析是應(yīng)用頻率特性研究線性系統(tǒng)的一種經(jīng)典方法。它以控制

105、系統(tǒng)的頻率特性作為數(shù)學(xué)模型，以伯德圖或其他圖表作為分析工具，來(lái)研究、分析控制系統(tǒng)的動(dòng)態(tài)性能與穩(wěn)態(tài)性能。在當(dāng)前的控制系統(tǒng)的分析應(yīng)用中起著重要作用。</p><p>　　5.3.1 頻率特性</p><p>　　穩(wěn)定的線性定常系統(tǒng)，在正弦信號(hào)作用下，系統(tǒng)輸出的穩(wěn)態(tài)分量為同頻率的正弦信號(hào)，其振幅與輸入正弦信號(hào)的振幅之比相對(duì)于正弦信號(hào)角頻率間的關(guān)系為幅頻特性；其相位與輸入正弦信號(hào)的相位之差相對(duì)于正

106、弦信號(hào)角頻率間的關(guān)系為相頻特性。系統(tǒng)頻率響應(yīng)與輸入正弦信號(hào)的復(fù)數(shù)比即為系統(tǒng)的頻率特性。即=</p><p>　　5.3.2 頻域響應(yīng)與Nyquist圖程序設(shè)計(jì)</p><p><b>　　系統(tǒng)的開(kāi)環(huán)函數(shù)為，</b></p><p>　　利用頻率響應(yīng)函數(shù)freqs()編寫程序：</p><p>　　num=4;den=[1

107、 2 4]</p><p>　　w=0:0.01:3;</p><p>　　g= freqs(num,den,w);</p><p>　　mag=abs(g)</p><p>　　plot(w,mag)</p><p>　　xlabel(‘頻率(rad./s)’);</p><p>　　ylab

108、el(‘幅值’)；</p><p><b>　　grid； </b></p><p>　　axis([0 3 0.5 1.2])</p><p>　　title(‘幅頻特性曲線’)</p><p>　　運(yùn)行后系統(tǒng)的幅頻特性曲線為：</p><p>　　Nyquist圖對(duì)頻率特性函數(shù)給出從到的一系列

109、數(shù)值，分別求出Im()和Re()，以Re()為橫坐標(biāo)、Im()為縱坐標(biāo)繪制出的極坐標(biāo)頻率特性圖。</p><p>　　5.3.3 設(shè)計(jì)分析</p><p>　　將傳遞函數(shù)按實(shí)部和虛部改寫為 </p><p>　　則和都是的實(shí)函數(shù)，曲線-和-分別為系統(tǒng)的實(shí)頻特性和虛頻特性曲線。如果將的虛部和實(shí)部分別作為縱橫坐標(biāo)，則曲線即為-奈魁斯特圖（Nyquist圖）。則

110、</p><p><b>　　,</b></p><p>　　由圖，奈奎斯特曲線逆時(shí)針包圍（-1，j0）點(diǎn)一圈，而開(kāi)環(huán)系統(tǒng)在右半平面有一個(gè)極點(diǎn)，所以系統(tǒng)穩(wěn)定。</p><p>　　在頻域分析中，此程序所應(yīng)用的頻率響應(yīng)法由于物理意義鮮明、計(jì)算量較小、直觀性較強(qiáng)等優(yōu)點(diǎn)，在實(shí)際應(yīng)用上非常廣泛。</p><p>　　5.4 線性

111、系統(tǒng)的穩(wěn)定性分析</p><p>　　穩(wěn)定是對(duì)控制系統(tǒng)提出的基本要求，也是保證系統(tǒng)正常工作的基本條件。</p><p>　　5.4.1 代數(shù)穩(wěn)定判據(jù)程序設(shè)計(jì)</p><p>　　求解控制系統(tǒng)閉環(huán)特征方程的根并判斷所有根的實(shí)部是否小于零，在matlab中可以調(diào)用roots()函數(shù)來(lái)實(shí)現(xiàn)。Root(P)函數(shù)輸入?yún)⒘縋是降冪排列多項(xiàng)式系數(shù)向量，在線性系統(tǒng)的穩(wěn)定性分析中，P

112、就是系統(tǒng)閉環(huán)特征多項(xiàng)式降冪排列的系數(shù)向量，如果能求出P，則其根也可以求出。</p><p>　　若系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為：</p><p>　　調(diào)用roots()函數(shù)編寫程序：</p><p><b>　　k=100;</b></p><p><b>　　z=-2;</b></p><

113、;p>　　p=[0,-1,-20];</p><p>　　[n1,d1]=zp2tf(z,p,k);</p><p>　　G=tf(n1,d1);</p><p><b>　　P=n1+d1</b></p><p><b>　　roots(P)</b></p><p>

114、　　運(yùn)行程序后可得多項(xiàng)式系數(shù)向量P及其根：</p><p><b>　　P =</b></p><p>　　1 21 120 200</p><p><b>　　ans =</b></p><p><b>　　-12.8990</b></p><

115、;p><b>　　-5.0000</b></p><p><b>　　-3.1010</b></p><p>　　由數(shù)據(jù)可看出所有特征根的實(shí)部均為負(fù)值，所以該閉環(huán)系統(tǒng)是穩(wěn)定的。</p><p>　　5.4.2 Bode圖法判斷系統(tǒng)穩(wěn)定性程序設(shè)計(jì)</p><p>　　立用求系統(tǒng)模值裕度和相位裕

116、度的函數(shù)margin()可以繪制處系統(tǒng)的Bode圖，同時(shí)根據(jù)Bode圖計(jì)算系統(tǒng)頻域性能指標(biāo)。</p><p>　　單位負(fù)反饋系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為：</p><p>　　調(diào)用margin（）函數(shù)編寫MATLAB程序如下：</p><p>　　num=[0 0 0 3.6];</p><p>　　den=[2 8 5 0];</p>

117、<p>　　s1=tf(num,den);</p><p>　　[Gm,Pm,Wcp,Wcg]=margin(s1)</p><p>　　margin(s1)</p><p>　　title('Bode圖')</p><p>　　ylabel('幅值');</p><p>&l

118、t;b>　　grid on</b></p><p>　　xlabel('頻率');</p><p>　　ylabel('相位');</p><p>　　執(zhí)行文件后，得到頻率性能指標(biāo)，如下：</p><p>　　模值穩(wěn)定裕度：Gm = 5.5556 即20*log10（5.5556）=14.89

119、46dB</p><p>　　相位穩(wěn)定裕度：Pm = 43.5933°</p><p>　　-π穿越頻率：Wcp = 1.5811 1/s</p><p>　　剪切頻率：Wcg = 0.5709 1/s</p><p>　　系統(tǒng)的Bode圖，如圖所示：</p><p>　　從上面的頻域性能指標(biāo)數(shù)據(jù)可以看出該系

120、統(tǒng)閉環(huán)是穩(wěn)定的，而且穩(wěn)定裕度也不錯(cuò)。</p><p>　　5.4.3 系統(tǒng)穩(wěn)定性概念及定義</p><p>　　（1）在經(jīng)典控制分析中，對(duì)線性定常系統(tǒng)穩(wěn)定性的概念是這樣定義的：若系統(tǒng)由于受到擾動(dòng)作用而偏離了原來(lái)的平衡狀態(tài)，當(dāng)擾動(dòng)去除后，如果能恢復(fù)到原來(lái)的平衡狀態(tài)，則稱該系統(tǒng)是穩(wěn)定的，否則該系統(tǒng)就是不穩(wěn)定的。</p><p>　?。?）穩(wěn)定性是系統(tǒng)自身的固有特性。穩(wěn)定