[工學]第四章 面向結構數(shù)字仿真法

1、第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法 對一個控制系統(tǒng)進行研究，其中一個很重要的問題就是考察系統(tǒng)中一些參數(shù)改變對系統(tǒng)動態(tài)性能的影響，面向微分方程的仿真方法很難得到這一點。這主要是由小回路的傳遞函數(shù)得到的全系統(tǒng)大回路的傳遞函數(shù)之間的參數(shù)對應關系將變得非常復雜。其次，將復雜系統(tǒng)中諸多小回路化簡求出總的系統(tǒng)模型也是十分麻煩的，更何況對于非線性系統(tǒng)，或難以用非數(shù)學模型描述的系統(tǒng)，則無法找到系統(tǒng)的總的閉環(huán)模型。第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法 本章介紹兩種由一些典型環(huán)節(jié)構成的復雜系統(tǒng)仿真的方法。在這類仿真程序中，先將仿真這些典型環(huán)節(jié)特性的仿真子程序編制好；在

2、仿真時，只要輸入各典型環(huán)節(jié)的參數(shù)以及環(huán)節(jié)間的連接關系的參數(shù)便可以作系統(tǒng)的仿真。這就是面向結構圖的數(shù)字仿真法，它可以解決上述困難，且具有一些優(yōu)點：第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法1.很容易改變某些參數(shù)環(huán)節(jié)，便于研究各環(huán)節(jié)參數(shù)對系統(tǒng)的影響。2.不需要計算出總的傳遞函數(shù)，并且可以直接得到各個環(huán)節(jié)的動態(tài)性能。3.系統(tǒng)中含有非線性環(huán)節(jié)時也比較容易處理。 第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法本章內容： 第一節(jié)介紹面向結構圖仿真各典型環(huán)節(jié)仿真模型的確定。 第二節(jié)介紹面向結構圖模型離散相似法仿真的方法。 第三節(jié)介紹對于含有典型非線性環(huán)節(jié)的處理方法。 第四節(jié)介紹連

3、續(xù)系統(tǒng)結構圖仿真方法、程序的編制及應用。4.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定n在第二章第四節(jié)中已經介紹了狀態(tài)方程離散化的方法，即對一個狀態(tài)方程加入虛擬的采樣器和保持器，當采樣頻率合適時則可實現(xiàn)信號重構。面向結構圖仿真方法其基本思想就是將結構圖化簡為各個典型環(huán)節(jié)組成，然后在各個典型環(huán)節(jié)前加入虛擬的采樣器和保持器使各環(huán)節(jié)獨自構成一個便于計算機仿真的差分方程。本節(jié)就是求出個典型環(huán)節(jié)對應的離散狀態(tài)方程得系數(shù)矩陣 ，即 。 1.積分環(huán)節(jié)n積分環(huán)節(jié)如圖4.1.1所示，其傳遞函數(shù)可寫為 (4.1.1) )(),(ttm)(tmsasusysg0)()()(4.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)

4、節(jié)仿真模型的確定狀態(tài)方程為 (4.1.2)根據(jù)(2.4.9)式可得 圖4.1.1 積分環(huán)節(jié)結構圖 其中 。離散狀態(tài)方程為 (4.1.3) xyuax02000)(0000)()()(1)(tabdbdettadabdetetttttamttttamat) 1() 1()0()0()()()() 1(200nxnyyxnutantuanxnx0; 0aba4.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定 2.比例積分環(huán)節(jié)n比例積分環(huán)節(jié)如圖4.1.2所示。顯見，狀態(tài)方程與積分環(huán)節(jié)一致，不同的是輸出方程、傳遞函數(shù)可寫為 (4.1.4)n其中： ， 。n根據(jù)(2.4.9)式，比例積分環(huán)節(jié)的狀態(tài)方程

5、和輸出方程可寫為 (4.1.4) 圖4.1.2 比例積分環(huán)節(jié)結構圖 100)(aasabsdscsgbca 0cda 1uaaxyuax1004.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定n顯見， ， 同積分環(huán)節(jié)一樣，僅離散狀態(tài)方程中的輸出方程與(4.1.3)是不一樣。即(4.1.6) 3.慣性環(huán)節(jié)n慣性環(huán)節(jié)的結構圖如圖4.1.3所示，其傳遞函數(shù)可寫為 (4.1.7) )(),(ttm)(tm) 1() 1() 1()0()0()(21)()() 1(10200nuaanxnyyxnutantuanxnx10)(asabsacsg4.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定其中:

6、 , 環(huán)節(jié)的狀態(tài)方程和輸出方程為 (4.1.8) 圖4.1.3 慣性環(huán)節(jié)結構圖根據(jù)(2.4.9)式，其差分方程的各項系數(shù)為bca 0baa 1xyuaxax01taet1)()1)()(11100)(0tattameaadeat4.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定離散狀態(tài)方程為 (4.1.9)4.比例慣性環(huán)節(jié)n比例慣性環(huán)節(jié)的結構圖如圖4.1.4所示。傳遞函數(shù)可寫為 (4.1.10)其中: , , 。 ) 1)()()(11210100)(0tattameaataadeat)0()0()1()1()()()()()()()1(yxnxnynutnutnxtnxmm11200)()

7、(asaaaabsadscsgbda 0baa 1dca 24.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定n狀態(tài)方程為 (4.1.11) 圖4.1.4 比例慣性環(huán)節(jié)結構圖n顯見狀態(tài)方程與慣性環(huán)節(jié)一樣，故 （t）， m（t）， m（t）的計算也一樣，僅輸出方程不一樣，故得離散狀態(tài)方程為 (4.1.12)uaxaayuaxax01201)() 0()() 0() 1() 1()() 1()()()()()()() 1(012ydcxnuanxaanynutnutnxtnxmm4.1 典型環(huán)節(jié)仿真模型的確定典型環(huán)節(jié)仿真模型的確定n處理上述幾種典型環(huán)節(jié)外，常用的還有二階環(huán)節(jié) ，它可由圖4.1.5

8、所示結構組成。圖4.1.5 二階環(huán)節(jié)等效結構圖n可見高階環(huán)節(jié)均可用前述幾種典型環(huán)節(jié)獲得。2120)(asasabsg4.2 結構圖離散相似法仿真結構圖離散相似法仿真n面向結構圖模型的離散相似法仿真除了需要建立典型環(huán)節(jié)的差分式外，還需要建立能描述系統(tǒng)連接方式的方程。在上一節(jié)的基礎上，本節(jié)將進一步介紹系統(tǒng)連接矩陣的建立和面向結構圖模型的離散相似法仿真方法以及計算程序的實現(xiàn)。一、連接矩陣n上一節(jié)，介紹了環(huán)節(jié)離散化方法以及所得到的差分方程模型的形式。但這僅僅表示了各個單獨環(huán)節(jié)輸入和輸出之間的關系。為了實現(xiàn)面向結構圖離散相似法仿真，還必須把這些環(huán)節(jié)按照系統(tǒng)結構圖的要求連接起來、以保證正確的計算次序。設系

9、統(tǒng)的第 個環(huán)節(jié)輸入、輸出分別用 表示， 為系統(tǒng)的外部輸入量，則 i),.,2 , 1(,niyuii0y4.2 結構圖離散相似法仿真結構圖離散相似法仿真 (4.2.1)n可把(4.2.1)式寫成 (4.2.2) 式中 是一個 維長方矩陣。這是把表示輸入信號與系統(tǒng)連接情況的 矩陣放在原連接矩陣的第一列，也就是n 001ywywuwyyywwu010nnnnnnwwwwwwwwww1022120111100w) 1( nnw4.2 結構圖離散相似法仿真結構圖離散相似法仿真 表示第 個環(huán)節(jié)輸入之間的連接方式。n 而 是一個 的列矢量， 。例如，有一系統(tǒng)如圖4.2.1所示。如果已知各環(huán)節(jié)的傳遞函數(shù)，則

10、很容易將其離散化，各個環(huán)節(jié)的輸入輸出關系為 (4.2.3) (4.2.4)ijwjy1) 1(ntnyyyyy210432104321 00100001001001001001yyyyyuuuuywu4.2 結構圖離散相似法仿真結構圖離散相似法仿真 圖4.2.1 系統(tǒng)結構圖 二、仿真程序的設計n把不同類型環(huán)節(jié)的離散系數(shù)的計算分別編成子程序。在程序中引入一個標志參數(shù) ，表示該典型環(huán)節(jié)的類型，假設一個通用程序只包括下列四種典型環(huán)節(jié)，且 與典型環(huán)節(jié)對應關系如下n當 =0時，表示第i個環(huán)節(jié)為積分環(huán)節(jié) 。n當 =1時，表示第i個環(huán)節(jié)為比例積分環(huán)節(jié) 。n當 =2時，表示第i個環(huán)節(jié)為慣性環(huán)節(jié) 。n當 =3時

11、，表示第i個環(huán)節(jié)為比例慣性環(huán)節(jié) 。)(ih)(ih)(ih)(ih)(ih)(ihscbsdsc bsacbsadsc4.2 結構圖離散相似法仿真結構圖離散相似法仿真n由前述可知，對于 =0和 =1兩種典型環(huán)節(jié)，計算狀態(tài)變量的公式相同，只是它們輸出變量計算公式不同。而同樣對于 =2和=3的典型環(huán)節(jié)，也是計算狀態(tài)變量的公式相同，僅僅是輸出方程不同。在步長取定后，典型環(huán)節(jié)的離散（t）， m（t）， m（t）， 圖4.2.2 仿真流程圖 )(ih)(ih)(ih)(ih4.2 結構圖離散相似法仿真結構圖離散相似法仿真n就僅是典型環(huán)節(jié)的參數(shù)（時間常數(shù)、放大增益）的函數(shù)，可以預先根據(jù)典型環(huán)節(jié)的類型分別編

12、成子程序，仿真時即可根據(jù) 方便地調用。n系統(tǒng)的連接情況，仍用連接矩陣w來描述。n面向系統(tǒng)結構圖離散化仿真的工作流程圖如圖4.2.2所示。n按系統(tǒng)的典型環(huán)節(jié)離散化仿真，其主要優(yōu)點是：（1）各個環(huán)節(jié)的離散狀態(tài)方程系數(shù)計算簡單，而且可以一步求出，不像龍格-庫塔法那樣，每一步都要重新計算龍格-庫塔系數(shù)，因而計算量相對來說較小。（2）由于各個環(huán)節(jié)的輸入量 ，輸出量 每一步都可求出，所以很容易推廣到含有非線性環(huán)節(jié)的系統(tǒng)仿真中去。)(ihiuiy4.2 結構圖離散相似法仿真結構圖離散相似法仿真n 該方法的主要缺點是計算精度低。因為每個環(huán)節(jié)的輸入實際上都是使用了它們的近似值（舉行近似或梯形近似），故僅有一階或

13、二階精度，這會帶來計算誤差，而且環(huán)節(jié)越多，誤差越大。這一點下面還將進一步分析。另外，需要指出的是，當輸入采用梯形近似法時，需要用到 來求取 ， ，這通常是難以辦到的。于是在仿真中有時只得采用簡單的向后差分的方法來計算 即 。由于 本來的定義是表示在 區(qū)間輸入信號的平均變化速度，所以用向后差分的方法來計算 實際上使用前一個周期 的輸入信號的平均變化速度來近似代替周期 的輸入信號變化速度，相差一個采樣周期。這顯然會使計算誤差增大。 ) 1( nu)(nu/ )() 1()(tnununu)(nu tnununu/)1()()()(nu nttn) 1( )(nu tn) 1( ntnttn) 1(

14、 4.2 結構圖離散相似法仿真結構圖離散相似法仿真三、仿真算例及分析n用該程序求某四階系統(tǒng)（結構圖見圖4.2.3）在階躍函數(shù)作用下的過渡過程。 圖4.2.3 四階系統(tǒng)結構圖n首先，確定典型環(huán)節(jié)類型和環(huán)節(jié)編號，本例從左到右順序排號 ，第一塊類型號 ，第二塊類型號 ，第三塊類型號 ，第四塊類型號 根據(jù)圖4.2.3所示可寫出連接矩陣為3) 1 (h0) 2 (h2) 3 (h2) 4 (h4.2 結構圖離散相似法仿真結構圖離散相似法仿真n根據(jù)經驗公式 ，大約可達到0.5%左右的精度， 為系統(tǒng)開環(huán)頻率特性的剪切頻率。在此例中， ， ，因此可選。輸入數(shù)據(jù)有01000001000001010001wct)

15、5030(1c1cstc02. 0033. 0)5030(14.2 結構圖離散相似法仿真結構圖離散相似法仿真n環(huán)節(jié)序號 初始值 初始值 1 1 1 0 0 2 0 1 1 0 0 0 3 1 0 0 0 4 1 0 0 0n連接矩陣abcdxy1a1b3a3a4a4a01000001000001010001w4.2 結構圖離散相似法仿真結構圖離散相似法仿真n仿真參數(shù) 采樣周期 仿真時間 打印、顯示時間 間隔輸出環(huán)節(jié)號1 0.01 10 1 1n輸入以上三組參數(shù)后，便可在計算機上仿真。四、采用補償器提高模型精度和穩(wěn)定性的方法n系統(tǒng)的離散化過程，就是在連續(xù)系統(tǒng)中加入虛擬的采樣開關和保持器。由于保持

16、器不可能完整無誤地將連續(xù)信號重構出來，因此必然會產生仿真誤差。一般來講，采樣間隔越大，仿真的誤差也就越大。為了減少誤差，很自然地就想到是否能在這個仿真器中模型中加進校正補償環(huán)節(jié)。一般所加入的補償 4.2 結構圖離散相似法仿真結構圖離散相似法仿真n器應盡可能好地抵消經過采樣-保持器所造成的失真，補償器常常采用超前的 的形式，其中 可以根據(jù)實際情況選取。整個仿真模型如圖4.2.4所示。n下面以積分環(huán)節(jié)為例來說明這種方法的基本原理。n假定 ，則按圖4.2.4所示構成的仿真模型的 為 (4.2.5) 圖4.2.4 加校正的數(shù)字仿真模型 ste,ssg1)()(zg)1(1)(2sezzzsesezzg

17、ststst4.2 結構圖離散相似法仿真結構圖離散相似法仿真n對 做一次近似，即去 (4.2.6)n則(4.2.5)是變成 (4.2.7)n寫成差分方程 (4.2.8)n選擇不同的 ，可得各種不同數(shù)值及分公式。比如stetsest 11)1 (1) 1()1(1)1()(22zztztzztzzzsstzzzzg)1 (11nnnnuutyy和4.2 結構圖離散相似法仿真結構圖離散相似法仿真n選 ， ，則有 歐拉公式n選 ， ，則有 梯形公式n選 ， ，則有 超前歐拉公式。 n在梯形公式及超前歐拉公式都有 項，一般它是未知的，在計算 時只知道 。為此，可以先對輸入信號加一拍延滯，然后再加大 ，

18、補償這種延滯所造成的誤差。如圖4.2.5所示，則有 (4.2.9) (4.2.10) 111021111nnntuyy)(211nnnnuutyynnntuyy1nyny1ny1)1 ()(1zztsg)1 (11nnnnuutyy4.2 結構圖離散相似法仿真結構圖離散相似法仿真n選 ， ，根據(jù)(4.2.10)式可得 (4.2.11)n這就是亞當斯公式。 圖4.2.5 補償延滯造成的誤差n由于 ， 可調，故將(4.2.8)及(4.2.10)式稱為可調整的數(shù)值計分共識。把這種方法用于復雜系統(tǒng)的快速仿真，就可以得出允許較大步距、又有一定精度的仿真模型。通常將這種方法稱為可調的數(shù)值積分法。132)3

19、(2211nnnnuutyy4.2 結構圖離散相似法仿真結構圖離散相似法仿真n當將這種方法用于復雜系統(tǒng)時，為獲得仿真模型，其基本步驟是：（1）在系統(tǒng)的輸入端加虛擬的采樣器及保持器，然后加 上補償環(huán)節(jié)，如圖4.2.4所示。（2）求出該圖所示的離散化系統(tǒng)的脈沖傳遞函數(shù) ，并列出它的差分方程，這就是仿真模型。（3） 用高階的龍格-庫塔法計算該系統(tǒng)的相應，將它作為一個標準解，然后給出不同的 ， ，計算仿真模型的響應，并將它與標準解進行比較，直到誤差達到最小為止。)(zg4.2 結構圖離散相似法仿真結構圖離散相似法仿真n利用上述步驟僅僅是計算出了系統(tǒng)的輸入量y。如果不僅對y感興趣，而且對于系統(tǒng)中的其他變

20、量也有興趣，那么就必須將系統(tǒng)分成幾個部分，每部分都要加虛擬的采樣器及保持器。至于校正補償環(huán)節(jié)則按一般系統(tǒng)的校正原則，可以對每一個小閉環(huán)加一個 。調整時，一般是先調外環(huán)的。調整的目標是要求所獲得的仿真模型在較大的計算步距時仍能最好地與實際模型相接近。 ， 地選取可采用第六章所介紹的參數(shù)尋優(yōu)程序來確定。ste4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真n在本章的上一節(jié)中曾提到，利用離散相似法編制的仿真程序雖然精度低，但是卻可以十分方便地推廣應用到這類非線性系統(tǒng)中去，其主要原因是在仿真計算程序中，每走一步，各個環(huán)節(jié)的輸入量及輸出量都將重新計算一次。因此非線性環(huán)節(jié)子程序很容易加入到仿真程序中去。下面

21、首先介紹典型的非線性環(huán)節(jié)的仿真。一、非線性環(huán)節(jié)仿真子程序一、非線性環(huán)節(jié)仿真子程序 1.飽和非線性n完成圖4.3.1所示飽和非線性特性輸入-輸出之間的仿真程序，可采用圖4.3.2所示的仿真流程圖，并相應地編制子程序在使用中調用。4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 圖4.3.1 飽和非線性特性 圖4.3.2 飽和非線性仿真程序流程圖4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 2.失靈區(qū)非線性n圖4.3.3所示的失靈區(qū)非線性特性輸入-輸出之間的仿真流程圖如圖4.3.4所示。 圖4.3.3 失靈區(qū)非線性特性 圖4.3.4 失靈區(qū)非線性仿真程序流程圖 4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)

22、的數(shù)字仿真3.齒輪間隙(磁滯回環(huán))非線性(圖4.3.5)n設 為上一次的輸入， 為上一次的輸出。若當 0時，且 ，則 = 即，若只滿足前一個條件，而不滿足后一個條件，則是工作在由左邊的特性過渡到右邊的特性上。 圖4.3.5 齒輪間隙非線性特性）（入mu0）（出mu0）（入mu0）（出mu0）（出mu0）（入mu0）（出mu0）（入mu01c1c4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 圖4.3.6 齒輪間隙非線性仿真程序流程圖 4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真n若當 0時，且 + ，則 = + 其他情況， = ，即輸出維持不變，正好在走間隙這一段。程序流程圖見4.3.6。二

23、、含有非線性環(huán)節(jié)的離散相似法仿真程序的計算方法二、含有非線性環(huán)節(jié)的離散相似法仿真程序的計算方法n當系統(tǒng)中有上述典型環(huán)節(jié)時，本章第二節(jié)講的離散相似法仿真程序要作如下修改： (1) 對每個環(huán)節(jié)要增設一個參數(shù) ，表示第 個環(huán)節(jié)的入口或出口有那種類型的非線性環(huán)節(jié)。 (2) 對每個環(huán)節(jié)要增設一個參數(shù) ，表示第 個環(huán)節(jié)的入口的那個非線性環(huán)節(jié)的參數(shù) ，當?shù)?個環(huán)節(jié)入口沒有非線性時， =0。）（入mu0）（出mu0）（出mu0）（入mu0）（入mu0）（出mu0）（出mu0）（出mu1c1c)(iziii)(icic)(ic4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真n因此在輸入數(shù)據(jù)時，對于每一個非線性環(huán)節(jié)

24、都要同時送 ， ， ， ， ， ， ， 8個數(shù)據(jù)。其中 ， ， ， 為線性環(huán)節(jié) 的系數(shù)， 的含義為： =0 表示該環(huán)節(jié)前、后無非線性環(huán)節(jié)。 =1 表示該環(huán)節(jié)前有飽和非線性環(huán)節(jié)。 =2 表示該環(huán)節(jié)前有失靈區(qū)非線性環(huán)節(jié)。 =3 表示該環(huán)節(jié)前有齒輪間隙非線性環(huán)節(jié)。 =4 表示該環(huán)節(jié)后有飽和非線性環(huán)節(jié)。)(ia)(ib)(ic)(id)(iy)(ix)(iz)(is)(ia)(ib)(ic)(idsbasdcxyiiiiii)(iz)(iz)(iz)(iz)(iz)(iz4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 =5 表示該環(huán)節(jié)后有失靈區(qū)非線性環(huán)節(jié)。 =6 表示該環(huán)節(jié)后有齒輪間隙非線性環(huán)節(jié)。n

25、的意義可以參見圖4.3.1，圖4.3.5所示的非線性參數(shù)。 (3) 一個完整的面向結構圖的離散相似法仿真程序框圖如圖4.3.7所示：)(iz)(iz)(is4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真三、非線性系統(tǒng)仿真舉例三、非線性系統(tǒng)仿真舉例n有一個四階非線性系統(tǒng)仿真結構圖如圖4.3.8所示。試分析：當階躍輸入 時：無非線性環(huán)節(jié)；非線性環(huán)節(jié)為飽和特性（見圖4.3.1），且 ；非線性環(huán)節(jié)為失靈區(qū)特性（見圖3.3.3）,且 ，這三種情況下的系統(tǒng)輸出響應，并分析結果。n第一步：確定系統(tǒng)各個環(huán)節(jié)號。本例除第一個環(huán)節(jié)前有非線性環(huán)節(jié)外，其余都為線性環(huán)節(jié)。n第二步：根據(jù)圖4.3.9所示寫出連接矩陣為10

26、u51c11c432104321yyy 01000001000001010001yyuuuu4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 圖4.3.7 離散相似法仿真程序框圖 4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 圖4.3.8 四階非線性系統(tǒng)結構圖 圖4.3.9 四階非線性系統(tǒng)仿真框圖4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真n第三步：運行程序根據(jù)提示輸入數(shù)據(jù)。輸入的數(shù)據(jù)有：（1）各環(huán)節(jié)參數(shù)，按第二種情況考慮，非線性環(huán)節(jié)為飽和特性，即 0.1 1 0.5 1 0 0 1 0 1 1 0 0 0 0 2 1 2 0 0 0 0 10 1 10 0 0 0 0（2）輸入連接矩陣數(shù)據(jù)

27、，即)(ia)(ib)(ic)(id)(0iu)(0iy)(iz4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真（3）輸入仿真參數(shù)。根據(jù)經驗公式 ，采樣周期按各環(huán)節(jié)最小的時間常數(shù)的選取，本例中最小時間常數(shù)為0.1，故采樣周期選為0.01。仿真時間去10s，且觀察第1號、4號環(huán)節(jié)輸出，因此仿真參數(shù)輸入時如下：01000001000001010001w1014.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真 采樣周期 仿真時間 打印、顯示時間間隔 輸出環(huán)節(jié)號 0.01 10 1 1，4n第五步：結果分析。將以上三組數(shù)據(jù)輸入到仿真程序中，運行后可得到數(shù)據(jù)結果。通過該程序的仿真結果，可以分析系統(tǒng)中典型非線性

28、環(huán)節(jié)對系統(tǒng)的影響。1.飽和非線性對系統(tǒng)過渡過程的影響n當自動控制系統(tǒng)中（非條件穩(wěn)定系統(tǒng)）存在飽和元件時，此時系統(tǒng)的穩(wěn)定性將變好，而快速性將變壞，也即超調量將減小，而過渡過程時間增加。這與自動控制原理理論分析結果相同。 4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真2.失靈區(qū)非線性對系統(tǒng)過渡過程的影響n根據(jù)調解原理分析可知：若系統(tǒng)中具有失靈區(qū)非線性環(huán)節(jié)，那么系統(tǒng)的動態(tài)品質將變壞，而對穩(wěn)定性影響不大。其原因是： (1) 由于有失靈區(qū)，在過渡過程的起始段，相當于減小了系統(tǒng)的開環(huán)增益，故過度過程變緩。 (2) 當輸入量接近穩(wěn)定時，放大器處于穩(wěn)定區(qū)，系統(tǒng)處于失控狀態(tài)，控制作用為零，故超調量將略微增大。

29、(3) 由于放大器有失靈區(qū)，故在過渡過程中有尾部，系統(tǒng)也處于失控狀態(tài)，因此將出現(xiàn)一個很長的尾巴。即從系統(tǒng)進入失靈區(qū)到輸出量進入穩(wěn)態(tài)值區(qū)（5%），輸出量變化十分緩慢。4.3非線性系統(tǒng)的數(shù)字仿真非線性系統(tǒng)的數(shù)字仿真3.齒輪間隙非線性對系統(tǒng)的影響n由于存在齒輪間隙，當系統(tǒng)的輸出值超過穩(wěn)態(tài)值時，因系統(tǒng)有反向調節(jié)的趨勢，輸出將維持不變一直要等非線性環(huán)節(jié)的輸入 走完間隙時輸出才能下降。而當輸出值反向偏離穩(wěn)態(tài)值時，系統(tǒng)郵政向調節(jié)的趨勢，輸出又將維持不變，一直要等非線性環(huán)節(jié)的輸入 走完間隙時輸出才能回升。其結果，系統(tǒng)將會在穩(wěn)態(tài)值附近以某一幅度和頻率進行震蕩，即系統(tǒng)始終在一個極限環(huán)內運動，而無法穩(wěn)定下來。uu4

30、.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序4.4.1 cssf程序包簡單介紹程序包簡單介紹n本節(jié)將介紹一個面向結構圖的數(shù)字仿真程序包css(continuons system simulation)。該程序包是1981年5月引入我國的。該程序包首先被移植在國產的djs-130機上，后因basic語言運行速度太慢，國內有關單位將該程序包翻譯成fortran語言，并對原程序作了大量的修改和補充，取名為cssf，zfx等。ncssf仿真程序包的主要特點： (1) 配備了多種積分方法。例如，定步長龍格-庫塔法、邊步長梅森法、定步長漢明法等。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖

31、仿真及程序 (2) 可以面向多種形式的數(shù)學模型。例如，狀態(tài)方程，一階及二階傳遞函數(shù)，n階傳遞函數(shù)或微分方程，對所有這些原始方程，只需直接輸入系統(tǒng)，不必進行交換。 (3) 增加了延遲、微分、一元及二元函數(shù)發(fā)生器等。 (4) 增加了采樣系統(tǒng)仿真功能，并且有參數(shù)優(yōu)化模塊。 (5) 具有非線性兩點邊值自動求解功能。ncssf程序包提供了約40種不同類型的運算塊，包括模擬機中的積分器、比例器、三角函數(shù)、對數(shù)、乘法、除法以及各種非線性函數(shù)運算模塊，故一套cssf程序包的計算功能相當于一臺的大型模擬機。除此之外cssf程序包還包括幾塊用戶自定義塊意識和用戶的特殊需要。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)

32、的結構圖仿真及程序n用cssf進行仿真的最大優(yōu)點就是使用方便。用戶只需將仿真系統(tǒng)轉化成程序包所含有的運算功能塊所組成的仿真結構圖，啟動程序包，按指定的方式輸入個功能塊的編號、類型參數(shù)、連接方式，以及積分步長、仿真時間等參數(shù)，就可以得到仿真結果。n4.4.2 micro-css仿真程序仿真程序n可以說，cssf仿真程序是一種功能全、性能高的程序包，可用于大型、復雜的連續(xù)或采樣系統(tǒng)的仿真，但該程序包的結構復雜。micro-css（mcss）是css的微型化，它是用basic語言編寫的?；镜膍css約有80條語句，而擴展mcss約有200條語句，并采用c語言編寫。它包括了css程序包的主要內容，并

33、且參數(shù)的輸入風格也與css相同，故通過mcss程序的學習，也可解剖css的結構，并掌握css的使用方法。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序nmcss仿真程序共有十七種函數(shù)功能塊，且為用戶提供了方便的擴展方式。仿真最大塊數(shù)定義為200塊，也可根據(jù)機器內存將其擴展。mcss仿真程序具有自動排序功能。其功能塊名稱見表4.4.1。每一功能塊用一種型號表示，例如，1型為常數(shù)快；2型為比例加法塊；3型為積分器；13型為采樣控制器模塊等。用戶使用時，首先要設定系統(tǒng)中每一功能塊的順序號（不是類型號），順序號的編號可由用戶自定，但最大號數(shù)不能超過200。實際運行次序有仿真程序在運行前自動

34、排列。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 表4.4.1 mcss仿真程序模塊名稱和功能 類型 名稱 符號輸入輸出之間的函 數(shù)關系1常數(shù)塊2比例加法器3比例微分器4延遲環(huán)節(jié)塊11pc 3423121pcpcpccdtpcpccpct)(34230211)(121ptcc4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序5正弦函數(shù)塊6余弦函數(shù)塊7乘法器8反號器9反正切函數(shù)10除法器)sin(21cc )cos(21cc 321ccc 21cc21arctancc 0,/3321cccc4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序11開關模塊12時間輸入

35、13采樣控制器模塊控制器系數(shù)采樣周期81802881188110,1)()()(bbaapzbzbzazaazdzuzytc 10024231ccccc4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序14理想繼電器15失靈區(qū)特性16包和特性17間隙特性4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序4.4.3 mcss仿真程序的使用方法仿真程序的使用方法n下面結合具體例子說明構成mcss仿真圖的方法、系統(tǒng)參數(shù)、連接情況，以及仿真參數(shù)的輸入方式。n例例4.4.1 有一個二階系統(tǒng)，其微分方程如下： 初值： ， 。試畫出仿真結構圖，并寫出仿真數(shù)據(jù)。n解解 1.仿真圖構成方法 由于

36、原微分方程可改寫為 1)(1 . 0)()(tytyty 0)0(y0)0(y )(1 . 0)(1)(tytyty 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 故可由積分器、常數(shù)塊構成如圖4.4.1所示的仿真圖，并在仿真圖上方標出順序號。 2.各功能塊連接方式輸入數(shù)據(jù)n各功能塊連接方式輸入數(shù)據(jù)的格式為 data i，16,j，k,l 其中 i-功能號順序號； 16-功能號類型號； j-同第1輸入端 連接的模塊號； k-同第2輸入端 連接的模塊號； l-同第3輸入端 連接的模塊號。)(2c)(3c)(4c4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 這部分數(shù)據(jù)輸入

37、結束時，用全零行表示，即 data 0,0,0,0,0 本例的數(shù)據(jù)格式為10 data 1，1，0，0，012 data 2，3，1，3，214 data 3，3，2，0，016 data 0，0，0，0，0 3.參數(shù)部分數(shù)據(jù)輸入 n參數(shù)部分數(shù)據(jù)輸入的格式為 圖4.4.1 例4.4.1的仿真結構圖 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 data i， ， ， n其中： i功能塊順序號 ， ， 各功能塊的參數(shù) ， ， 這部分數(shù)據(jù)輸入結束時，同樣用全零行表示。n注意：積分塊（3型功能塊）的第一參數(shù)輸入為初始條件數(shù)據(jù)。當某一功能塊中 ， ， 三個參量都沒有時，可以不輸入該塊的數(shù)

38、據(jù)。例如，正弦函數(shù)塊、余弦函數(shù)塊、乘法器等。n本例的數(shù)據(jù)格式為) 1 (y) 1 (y)2(y)2(y)3(y)3(y1p2p3p1p2p3p4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 60 data 1，1， 0， 0 62 data 2，0， -0.1， -1 64 data 3，0， 0， 0 66 data 0，0， 0， 04. 仿真控制數(shù)據(jù)輸入n仿真控制數(shù)據(jù)輸入的格式如下： 80 data d，t1，t2，k1，k2，k3，k44.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n其中 d積分步長； t1仿真時間； t2打印、顯示時間間隔； k1k4輸出量所在

39、模塊。n本例中，如果設計步長為0.1s，仿真時間為10s，每秒顯示一次，輸出量為，則數(shù)據(jù)格式為90 data 0.1，1，3，2，0，0n 其輸出結果如下：4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 time output3 output2 output0 output0 0 0 0 0 00.999 996 0.365 040 3 0.721 808 0 01.999 999 1.103 401 0.811 573 8 0 02.999 998 1.932 805 0.830 502 8 0 03.999 997 2.743 704 0.785 397 5 0 04.999

40、996 3.496 954 0.719 563 2 0 05.999 995 4.181 771 0.650 230 8 0 06.999 994 4.798 479 0.583 958 2 0 07.999 993 5.351 465 0.522 966 2 0 09.000 016 5.846 339 0.467 740 8 0 04.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n例例4.4.2 設有一非線性系統(tǒng)如圖4.4.2所示，求其響應 。 圖4.4.2 例4.2.2系統(tǒng)結構圖 )(ty4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n解解 首先畫出其仿真結構圖，如

41、圖4.4.3所示。 圖4.4.3 例4.2.2系統(tǒng)仿真結構圖 1 . 04 . 011 . 05 . 0sss4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n數(shù)據(jù)文件為 10 data 1，1，0，0，0 12 data 2，2，1，8，1 14 data 3，16，2，0，0 16 data 4，3，0，3，0 18 data 5，3，4，5，0 20 data 6，3，5，6，0 22 data 7，3，0，6，7 24 data 8，2，6，7，0 26 data 0，0，0，0，0 50 data 1，1，0，04.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 5

42、2 data 2，1，-1，0 54 data 3，1，-1，0 56 data 4，0，4，0 58 data 5，0，-2，0 60 data 6，0，-10，0 62 data 7，0，0.4，-0.1 64 data 8，1，1，0 66 data 0，0，0，0，0 90 data 0.01，5，0.5，8，0，0，04.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n其輸出結果為 time output3 output2 output0 output0 0 0 0 0 00.499 999 8 2.753 746e-02 0 0 00.999 999 3 0.110 005

43、 2 0 0 01.499 999 0.226 193 0 0 01.999 998 0.366 024 2 0 0 02.499 998 0.525 232 6 0 0 02.999 998 0.701 728 2 0 0 03.499 997 1.101 851 0 0 04.500 008 1.323 749 0 0 04.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序4.4.4 mcss仿真程序分析仿真程序分析n為了方便讀者進一步該程序的使用方法和應用范圍，并熟悉變成技巧，本節(jié)對該程序的組成及功能作更深入的分析和說明。nmcss程序共有以下四部分組成： (1) 輸入數(shù)據(jù)塊； (

44、2) 自動安裝模塊順序塊； (3) 運行程序快； (4) 輸出打印塊； 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n該程序編排緊湊，有些塊又是相互交叉進行的。圖4.4.4為mcss仿真程序流程圖。n下面分別說明以上四部分程序。1.輸入數(shù)據(jù)塊n這部分主要有讀語句組成，輸入的數(shù)據(jù)有三部分，按先后順序為 (1) 連接方式數(shù)據(jù)：i，16,j，k,l (2) 典型塊參數(shù)數(shù)據(jù)：i， (3) 仿真控制數(shù)據(jù)：d，t1，t2，k1，k2，k3，k4 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n其中（1），（2）兩類輸入數(shù)據(jù)結束標志為零。數(shù)據(jù)段可放在程序的前部或后部，但三部分數(shù)據(jù)的

45、次序不能改變。n在編程時考慮到程序的緊湊，故將這部分插入在其他部分之中。在讀數(shù)據(jù)的同時為了后一部分排序的要求，計算機自動檢查仿真系統(tǒng)中積分塊、常數(shù)塊、延時塊的個數(shù)即編號。安排方式如表4.4.1所示。 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 表4.4.1 三類特殊模塊的編號類型數(shù)目編號延遲塊n4id(n4)積分塊n9v(n9)常數(shù)塊n8n(n8)4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序2.自動安排各模塊計算順序n這一段程序的主要功能是編排各模塊計算順序并賦予順序號。用戶在完成仿真結構圖后，先要給你一個計算順序，即編號。至于這一順序號是否合理，程序會自動進行判

46、斷，并作出合理安排。前排的原則是，該模塊的三個輸入信號是否已有確定值，若無，則這一模塊必須排在以后計算。但有三種類型塊例外：常數(shù)塊(1)、積分塊(3)、和延遲器(4)。它們的輸出值可由給定的參數(shù)和初始值來確定，不必依賴于輸入值。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序圖4.4.4 mcss流程圖 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n在進入這段程序以前，除了常值信號源外(因為它們不必計算，故沒有安排計算順序的必要)，其他運算塊的類型代碼(程序中用u(i，1)表示)，均為負值，以表示尚未安排好的順序。程序根據(jù)用戶編號順序從小到大檢查類型代碼u(i，1)的符

47、號，如果u(i，1)為負，就開始檢查該塊的三個輸入端信號是否都有確定值，即檢查三個輸入信號是否來自延遲、積分、常數(shù)或已安排好的功能模塊的輸出。如果三個輸入均滿足上述情況之一，即可安排該塊的計算順序并編號，且將類型代碼u(i，1)取正值。該塊的流程圖如圖4.4.5所示。其中，n7為環(huán)節(jié)總數(shù)目；n()位數(shù)組存放排好的計算順序；i，k，l1為該子程序用到的變量。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 圖4.4.5 自動排序模塊流程圖 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n下面用一個例子來說明mcss程序中各模塊計算順序自動安排的過程。n例例 4.4.3 有一

48、方程如下 設初始條件為 試分析計算機排序的方法和結果。n解：解：首先將方程改寫為 n用上式可畫出仿真結構圖。先對仿真結構圖中各模塊任意編號。如圖4.4.6 所示。0)()()1()(2tytyypty 0)0(y0.1)0(y )()()()()(2tytytpytypty 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 圖4.4.6 例4.4.3仿真結構圖 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n程序先按用戶編號1號開始檢查，因為1號比例器只有第3輸入端是確定的，其余兩個未定，所以不能列為第1號。接著檢查第2號，由于2號積分器的輸入端需要1號提供，故也不能列為

49、第1號。而3號積分器的輸入端由2號積分器的初始值便可確定，故列為第1號。這樣反復檢查，最后得到的計算順序為 345612n顯然這種計算順序并非唯一，它和用戶編號順序有關。若把原用戶編號作如下改動，即 13，24，35，41，52，664.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n則計算順序為 125634n雖然順序不同，只要確定正確無誤，那么其計算結果完全相同。3.運行程序塊n該段主要包括（1）運行準備：指針 a 龍格-庫塔法系數(shù)； b 典型塊運行控制。（2）典型塊計算，目前mcss仿真程序共有17種典型塊，且程序為用戶提供了方便的備用模塊入口地址，以適應特殊需要。另外程序為了突

50、出簡單，積分方法只選擇了四階龍格-庫塔法。 4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序4.輸出數(shù)據(jù)塊nmcss仿真程序追求的是簡單，故輸出數(shù)據(jù)的方式很簡單，只有三條語句為數(shù)據(jù)輸出語句。每次運行最多可將4個模塊的數(shù)據(jù)輸出。4.4.5 代數(shù)環(huán)問題n在連續(xù)時間系統(tǒng)面向結構圖仿真中，完全可能出現(xiàn)這種情況，即計算機無法對所有的環(huán)節(jié)進行排序。換句話說，在系統(tǒng)的仿真結構圖中，可能存在一些環(huán)節(jié)無法安排計算順序。這種情況我們稱系統(tǒng)的仿真結構圖出現(xiàn)了代數(shù)環(huán)。這是需要特別注意的問題，因為這時，系統(tǒng)將不可能進行仿真計算。所謂的代數(shù)環(huán)就是在系統(tǒng)的仿真結構圖中出現(xiàn)了純粹是有代數(shù)環(huán)節(jié)構成的閉環(huán)回路。 4.4

51、連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n例例 4.4.4 考慮下面方程所描述系統(tǒng)的仿真 n解：我們首先畫出系統(tǒng)的仿真結構圖，如圖4.4.7所示。22xyyyxyx4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序 圖4.4.7 例4.4.4的仿真結構圖4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n對圖4.4.7分析，可以發(fā)現(xiàn)由于模塊2無法確定計算順序，導致模塊3、模塊4、模塊5和模塊7也無法確定計算順序，也這是模塊7無法確定計算順序，有導致模塊2無法確定計算順序。這就是說仿真結構圖出現(xiàn)代數(shù)環(huán)，即由模塊2、模塊4和模塊7形成一個封閉環(huán)，在該封閉環(huán)內沒有動態(tài)環(huán)節(jié)

52、，而僅有一些代數(shù)環(huán)節(jié)組成?；蛘咭部梢赃@樣說，如果某環(huán)節(jié)的輸出僅僅經過代數(shù)運算就反饋到該環(huán)節(jié)的某個輸入，則就構成一個代數(shù)環(huán)。n代數(shù)環(huán)的解決方法很多，可以根據(jù)被仿真對像的具體情況來定。在simulink里面，采用了newton-raphson技術來求解代數(shù)環(huán)。由于這種方法利用的是迭代優(yōu)化，因此存在計算時間長，甚至找不到解的可能性。在此我們介紹兩種比較實用的方法。4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n 如果系統(tǒng)的仿真結構圖比較復雜，并且計算步長也較小，可以在代數(shù)環(huán)上的任意一個地方增加一個一步延時環(huán)。這種方法比較簡單，并且有效。n 如果系統(tǒng)的模型不是特別復雜，可以將模型重寫為一個

53、標準的一階微分方程組，例如例4.4.4的系統(tǒng)可以重寫為)2(12)2(1122xyyyyxyyx4.4連續(xù)系統(tǒng)的結構圖仿真及程序連續(xù)系統(tǒng)的結構圖仿真及程序n即按上式重新畫仿真結構圖，就可以避免代數(shù)環(huán)的出現(xiàn)。因此，一般按標準的一階微分方程組，即方程右邊不出現(xiàn)變量的導數(shù)，畫出的仿真結構圖不會有代數(shù)環(huán)。第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法本章小結本章小結n本章是全書的重點之一，所介紹的兩種面向結構圖仿真方法是目前科學研究和核工程實踐中常用的仿真方法。 (1) 當連續(xù)系統(tǒng)用結構圖形式給定后，離散相似法是一種較為簡單的方法，該方法的實質就是在系統(tǒng)必要環(huán)節(jié)的輸入和輸出端加入虛擬采樣器

54、和保持器，將連續(xù)系統(tǒng)離散化，然后分別計算分割開的各個環(huán)節(jié)的輸出量，并按結構圖上的關系把相應的輸入與輸出連接起來，順序求解計算。由于環(huán)節(jié)的離散化方程可離線計算，因此該方法突出的特點是運算速度快，但精度低。第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法（2）相對于離散相似法而言在本章第四節(jié)敘述的連續(xù)系統(tǒng)結構圖仿真方法應用更廣泛。它們的基本原理是一樣的，但后者對各環(huán)節(jié)的計算是在仿真運行中分別計算的，因此它可以提供多種積分解法供用戶選擇，并且有較高的精度，但帶來的問題是計算速度較慢，并且需要排序，以避免結構圖中出現(xiàn)的代數(shù)環(huán)。第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法習習 題

55、題n4-1 有一閉環(huán)系統(tǒng)如下圖所示。 （1）求出的，列出求解的差分方程。（2）求出閉環(huán)系統(tǒng)的，列出求解的差分方程。第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法 （3）在閉環(huán)入口處加虛擬采樣器及零階保持器，求出開環(huán)的脈沖傳遞函數(shù)，并列出求解的差分方程。 （4）在系統(tǒng)入口處加虛擬采樣器及零階保持器，求出開環(huán)的脈沖傳遞函數(shù)，并列出求解的差分方程。n4-2已知系統(tǒng)： （ ）的 = ， ，故差分方程為n試分析步距t應如何選擇。若選擇的過大，計算時是否會發(fā)生不穩(wěn)定？ ask0a)(tate)1 ()(atmeaktnatnatnueakyey)1 (1第四章第四章 面向結構圖的數(shù)字仿真法面向

56、結構圖的數(shù)字仿真法n4-3 試編出下圖所示的非線性環(huán)節(jié)的仿真程序。 第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n4-4 設有系統(tǒng)狀態(tài)方程如下：（1） （2） 試編制程序，求出這兩組狀態(tài)方程轉移矩陣和。已知采樣周期t=0.5。uxxxx1014. 335.16102121uxxxxxx10011005 . 00015 . 0321321第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n4-5 試用mcss仿真程序對下圖所示系統(tǒng)進行仿真。其中，輸入u為單位階躍函數(shù)；參數(shù) ， 。分析在離散化采樣周期（步長）t=0.02和t=0.05兩種情況下，仿真結果是否相同。 11c05

57、. 02c第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法4-6 設有一系統(tǒng)結構如下圖所示。當輸入為單位階躍函數(shù)時，輸出y的精確解為 n試用離散仿真程序對上述系統(tǒng)在以下幾種條件下進行仿真，并與精確解y(t)相比較。tetetetetyttttsin813cos8113sin24173cos831)(33第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法（1）在離散化采樣周期 0.02，0.005，0.001，0.002四種情況下的仿真結果與精確解的誤差。（2）在離散化采樣周期 0.005條件下，對于把環(huán)節(jié) 與 合在一起離散化以及兩個環(huán)節(jié)分開離散化 兩種情況，若將兩者所得結果同

58、精確解比較，那一種誤差更小一些？0d59. 420sss50d第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n4-7 試用mscss程序仿真一個具有延遲環(huán)節(jié)的系統(tǒng)，其結構圖如下圖所示。離散化采樣周期 0.01。 0d第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n附錄：附錄：mcss仿真程序原碼仿真程序原碼n500 data 1,1,0,0,0n510 data 2,2,1,8,0n520 data 3,16,2,0,0n530 data 4,3,0,3,0n540 data 5,3,4,5,0n550 data 6,3,5,6,0n560 data 7,3,0,6,7

59、n570 data 8,2,6,7,0第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n580 data 0,0,0,0,0n590 data 1,1,0,0n600 data 2,1,-1,0n610 data 3,1,1,0n620 data 4,0,4,0n630 data 5,0,-2,0n640 data 6,0,-10,0n650 data 7,0,.4,-.1n660 data 8,1,1,0n670 data 0,0,0,0第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n670 data 0,0,0,0n680 data .05,30,.2,8,2,6,7

60、n690 dim u(200,5),v(60),p(200,3),c(200),e(60,4),y(200),tt(200),ck4(200,5)n700 dim s(20,200),n(200),id(20),ad(9),bd(9),cy(9),cu(9)n710 n9=0:t=0:c(0)=0:jnp=1n720 print input i,i6,j,k,ln730 input j,u(j,1),u(j,2),u(j,3),u(j,4)第四章第四章 面向結構圖的數(shù)字仿真法面向結構圖的數(shù)字仿真法n740 if j=0 goto 830n750 n7=n7+1:i6=u(j,1):if i64