【正文】 曲線如下：　　　　與初始值 k＝1時的階躍響應比較（下圖），顯然優(yōu)化設計后，動態(tài)性能得到了很好的改善?！　　　栴}46：試分析圖573所示系統(tǒng)中死區(qū)非線性對系統(tǒng)動態(tài)性能的影響?！　　　〗獯穑阂驗殡A躍給定是10，如果c大于10，系統(tǒng)輸出響應就為零，研究它已經沒有意義，因此死區(qū)特性參數(shù)c應滿足。在Simulink環(huán)境下搭建系統(tǒng)模型如下：　　　　分別在c=0(此時為倍數(shù)為一的比例增益)，2，4，6，8的情況下，仿真得到系統(tǒng)的輸出曲線如下：　　　　由此可以看出，引入非線性環(huán)節(jié)后隨著c值的增加，系統(tǒng)動態(tài)響應速度變慢，系統(tǒng)超調減少，阻尼系數(shù)變大，當c=8時系統(tǒng)的輸出響應已經表現(xiàn)出過阻尼特性。　　問題47：如圖537所示計算機控制系統(tǒng)，試設計一最小拍控制器,并用仿真的方法分析最小拍控制器對系統(tǒng)輸入信號和對象參數(shù)變化的適應性?！　　　〗獯穑涸诖祟}中，為表述方便，設表示系統(tǒng)閉環(huán)脈沖傳遞函數(shù)：　　表示系統(tǒng)誤差脈沖傳遞函數(shù)：　　表示數(shù)字控制器脈沖傳遞函數(shù)；　　表示控制對象的z變換；表示輸入信號的z變換；表示數(shù)字控制器輸入信號的z變換。　　將被控對象離散化后可以得到　　　　本題中，T=，那么 　　也可用下列MATLAB命令語句求出：　　sysc=zpk([],[0 1],[10]):　　c2d(sysc,)　　得到結果為，可以看出，與理論計算結果相同，但顯然用MATLAB命令語句求解更加簡單方便?！　。?）首先按照快速有紋波最小拍控制系統(tǒng)設計?！　】刂茖ο蟮膫鬟f函數(shù)為，由此可見，控制對象中含有單位圓上的極點p=1,為滿足穩(wěn)定性要求，必須包含z=1的零點；因滯后一拍，閉環(huán)脈沖傳遞函數(shù)至少也要滯后一拍。系統(tǒng)針對等速輸入設計，綜上所述，設=?！　闈M足穩(wěn)態(tài)精度，列寫下列關系：　　　　解得，,即 　　由此可得，=　　建立SIMULINK仿真模型如下：　　　　用SIMULINK仿真觀察設計結果，得到系統(tǒng)的輸出響應曲線如下圖所示，　　　　數(shù)字控制器的輸出曲線如下圖所示，　　　　由上面兩個圖可見，快速有紋波系統(tǒng)其輸出值跟隨輸入值后，在非采樣時刻有明顯的紋波，而數(shù)字控制器振蕩收斂，且振蕩持續(xù)時間很長。這是由于數(shù)字控制器　　=有一個p=，其值與p=1接近，由此造成數(shù)字控制器輸出振蕩收斂。這種情況下，不僅造成非采樣時刻有誤差，而且浪費執(zhí)行機構的功率，增加機械磨損。　　因此，我們可以考慮按照無紋波最小拍控制系統(tǒng)設計?！　。?）按照無紋波最小拍控制系統(tǒng)設計?！　】刂茖ο螅纱丝梢?，控制對象中含有單位圓上的極點p=1,以及零點z=。若想令輸出無紋波，則必須包含的全部零點（不管其在z平面的位置）。同時，為滿足穩(wěn)定性要求，必須包含z=1的零點；因滯后一拍，閉環(huán)脈沖傳遞函數(shù)至少也要滯后一拍。系統(tǒng)針對等速輸入設計，綜上所述，設　　=　　為滿足穩(wěn)態(tài)精度，列寫下列關系：　　　　解得，, 　　即 　　由此可得，=　　建立SIMULINK仿真模型如下：　　　　用SIMULINK仿真觀察設計結果，得到系統(tǒng)的輸出響應曲線如下圖所示，　　　　數(shù)字控制器的輸出曲線如下圖所示，　　　　從上兩圖可見，和有紋波系統(tǒng)相比，雖然無紋波系統(tǒng)的輸出跟隨輸入函數(shù)所需時間要多一個采樣周期，但采樣點之間不存在紋波，且數(shù)字控制器輸出波形較好，控制質量較高。　　(3)按照無紋波最小拍系統(tǒng)設計方法，分析最小拍控制器對系統(tǒng)輸入信號和對象參數(shù)變化的適應性。　　?！　‘斴斎霝殡A躍信號時，在Matlab/Simulink環(huán)境下建立下圖所示的控制系統(tǒng)結構圖?！　　　》抡娴玫较到y(tǒng)輸出曲線如下：　　　　仿真得到數(shù)字控制器輸出曲線如下：　　　　由此可見，當輸入為階躍輸入時，控制器輸出正負跳動很大，而且系統(tǒng)輸出超調嚴重，達到了100%?！　‘斴斎霝榧铀俣刃盘枙r，在Matlab/Simulink環(huán)境下建立下圖所示的控制系統(tǒng)結構圖，　　　　仿真得到系統(tǒng)輸出曲線如下圖所示,為檢驗控制系統(tǒng)的性能，將輸入輸出曲線繪制在一張圖上。　　　　可以看出，系統(tǒng)對于加速度輸入時，由明顯的穩(wěn)態(tài)誤差?！　【C上所述，針對一種典型的輸入函數(shù)進行設計，得到的最小拍系統(tǒng)閉環(huán)脈沖傳遞函數(shù)，用于次數(shù)較低的輸入函數(shù)時，系統(tǒng)將出現(xiàn)超調，響應時間也將增加；而用于次數(shù)較高的輸入函數(shù)時，輸出將產生穩(wěn)態(tài)誤差。最小拍控制器對系統(tǒng)輸入信號適應性不好?！　?。　　設控制對象的傳遞函數(shù)為，本題中。為檢驗最小拍控制器對于對象參數(shù)變化的適應性，令，得到此時的控制系統(tǒng)輸出曲線如下圖所示，　　　　可以看出，所設計的最小拍控制器對于對象參數(shù)的變化，表現(xiàn)出了一定的適應性。將參數(shù)b做較大范圍的修改后，輸出仍然能夠跟蹤給定，但是恢復時間略有增加。　　問題48：為使圖575所示系統(tǒng)不產生自激振蕩，試分析a、b取值?！　　　〗獯穑壕€性部分傳遞函數(shù)為　　因此　　令線性部分頻率響應的虛部為零，即　　　　求得曲線與負實軸交點處的頻率為rad/s，將其代入到中，得到 　　因此曲線與負實軸的交點坐標為（，j0）。 　　這可以通過Matlab相關命令加以驗證，在Matlab窗口中輸入下列命令：　　num=3?！　en=conv([1,0],conv([,1],[1,1]))?！　yquist(num,den)　　得到線性部分的奈奎斯特圖如下圖所示，　　　　將圖形局部放大后，同樣可得出其與負實軸交點坐標為（，j0），與理論分析結果一致?！　D中所示系統(tǒng)為具有死區(qū)的繼電器非線性系統(tǒng)，這一非線性特性的描述函數(shù)為　　　　據(jù)此可畫出軌跡，圖中的拐點為。根據(jù)非線性控制理論，若使系統(tǒng)不產生自激振蕩，必須保證　　軌跡與軌跡不相交，因此需令　　　　可求得，可取，此時滿足不產生自激振蕩條件?！　　　栴}49：已知某一地區(qū)在有病菌傳染下的描述三種類型人數(shù)變化的動態(tài)模型為　　　　式中，X1表示可能傳染的人數(shù)；X2表示已經得病的人數(shù)；X3表示已經治愈的人數(shù)；。試用仿真方法求未來20年內三種人人數(shù)的動態(tài)變化情況?！　〗獯穑壕帉慚函數(shù)，：　　function dx=rs(t,x) 　　dx=[*x(1)*x(2)。*x(1)*x(2)*x(2)。*x(2)]。 　　在MATLAB命令窗口中輸入下列指令：　　t_end=20?！　?設定仿真終止時間為20年　　x0=[620。10。70]。　　%三種類型人數(shù)初值　　[t,x]=ode45(39。rs39。,[0,t_end],x0)?！　lot(t,x(:,1),39。b39。,t,x(:,2),39。r39。,t,x(:,3),39。k.39。)　　xlabel(39。時間（年）39。)。　　ylabel(39。三種類型人數(shù)（人）39。)?！　egend(39。X1可能傳染的人數(shù)39。,39。X2已經得病的人數(shù)39。,39。X3已經治愈的人數(shù)39。)?！　itle(39。20年內三種人人數(shù)動態(tài)變化曲線圖39。)　　x。　　%查看三種類型人數(shù)　　floor(x)　　%人數(shù)應取整數(shù)　　運行程序后可得到20年內三種人人數(shù)動態(tài)變化曲線，表格中列出了每年的三種類型人數(shù)。年數(shù) 012345678910X1620611598575539487418337254182124X21017294980124182248311359389X370707275798798113133158185　　年數(shù) 11121314151617181920X18356372517129653X2402400390375356336316297278260X3213242271299325350373396416436　　　　問題50：對于高階系統(tǒng)的設計問題，往往要進行降階近似處理，并要驗證近似效果。已知某高階系統(tǒng)模型為　　　　經簡化處理后，模型等效為 　　試比較兩個模型在單位階躍信號作用下的響應情況，并分析近似效果。 　　解答：在Matlab命令窗口中輸入下列命令：　　num1=[35,10861,13285,82402,278376,511812,422964,194480]。 　　den1=[1,33,437,3017,11870,27470,37492,28880,9600]?！　ys1=tf(num1,den1)。　　num2=[35,]?！　en2=[1,]?！　ys2=tf(num2,den2)?！　tep(sys1,39。39。,sys2,39。:r39。,15)　　legend(39。原模型39。, 39。簡化后模型39。) ?！　rid on　　程序運行結果如下：　　　　從曲線中可以看出，降階后系統(tǒng)響應無超調，調整時間縮短，系統(tǒng)響應變快，但是簡化前后終值有差異。第58頁（共5頁）