【正文】 , 39?！　en2=[1,]。)　　x。X1可能傳染的人數(shù)39。k.39?！　?三種類型人數(shù)初值　　[t,x]=ode45(39。試用仿真方法求未來20年內(nèi)三種人人數(shù)的動(dòng)態(tài)變化情況?！　　　〗獯穑壕€性部分傳遞函數(shù)為　　因此　　令線性部分頻率響應(yīng)的虛部為零，即　　　　求得曲線與負(fù)實(shí)軸交點(diǎn)處的頻率為　　　　可以看出，系統(tǒng)對于加速度輸入時(shí)，由明顯的穩(wěn)態(tài)誤差。若想令輸出無紋波，則必須包含的全部零點(diǎn)（不管其在z平面的位置）?！　】刂茖ο蟮膫鬟f函數(shù)為，由此可見，控制對象中含有單位圓上的極點(diǎn)p=1,為滿足穩(wěn)定性要求，必須包含z=1的零點(diǎn)；因滯后一拍，閉環(huán)脈沖傳遞函數(shù)至少也要滯后一拍。,[1])　　運(yùn)行結(jié)果為：　　result =　　用SIMULINK仿真驗(yàn)證所的結(jié)果，階躍響應(yīng)曲線如下：　　　　與初始值 k＝1時(shí)的階躍響應(yīng)比較（下圖），顯然優(yōu)化設(shè)計(jì)后，動(dòng)態(tài)性能得到了很好的改善?！　tt,xx,yy]=sim(39?！　　　≡诶碚撋?，Smith預(yù)估補(bǔ)償可以克服大滯后的不利影響，但是Smith預(yù)估器需要知道被控對象精確的數(shù)學(xué)模型，并且對模型的誤差非常敏感，因而在實(shí)際應(yīng)用中仍面臨很多問題?！　栴}44：試采用smith預(yù)估控制方法對題56所述系統(tǒng)進(jìn)行重新設(shè)計(jì)，并用仿真的方法分析滯后參數(shù)變化對系統(tǒng)動(dòng)態(tài)性能的影響。 　　在MATLAB命令窗口中輸入下列指令:　　global kp?！　p=x(1)。k=200*(****)sys=zpk(z,p,k)　　　　　　　　　　%求出系統(tǒng)的開環(huán)傳遞函數(shù)[gm,pm,wcg,wcp]=margin(sys)　　　%求出系統(tǒng)的幅值裕度、相位裕度及其所對應(yīng)的頻率sys_closeloop=feedback(sys,1)　　%求出系統(tǒng)的閉環(huán)傳遞函數(shù) bode(sys)　　　　　　　　　　　　 %繪制系統(tǒng)的伯德圖 grid onfigure(2)step(sys_closeloop,5)　　　　　　%繪制閉環(huán)系統(tǒng)的階躍響應(yīng) grid on運(yùn)行程序后得到校正后系統(tǒng)的剪切頻率，這些都滿足期望的性能指標(biāo)，同時(shí)系統(tǒng)伯德圖如下所示。中頻段與低頻段之間用斜率為40dB/dec的直線連接，連接線與中頻段相交的交接頻率與不宜離的太近，??；同時(shí)取，則,?！　label(39。Gain(dB)39?！　um=1。Phase(deg)39。)?！　um=1。Frequency(rad/s)39。)?！　en=conv(a,b)。　　rlocus(num,den)　　grid on　　可得到系統(tǒng)的根軌跡如下圖所示：　　　　問題39：某反饋控制系統(tǒng)的開環(huán)傳遞函數(shù)為試?yán)L制其根軌跡?！　　　=[t,K1*T]。t=0。H=KP*(1+2*TD/T)。h=?！　?　　對上式進(jìn)行Z逆變換，得到　　　　由此可編寫仿真程序。 　　　　　　　　end　　　　　　　　if(Z(i)==6)　　　　　　　　[Y(i),Ubb(i)]=backlash(Ubb(i),Y(i),Yb(i),S(i))。.*Uk+FIJ39?！　　　　　k=W*Y+W0*Y0?！　　　　　lse　　　　　　end　　　　end　　end　　Y=zeros(n,1)?！　　　　　lse　　　　　　end　　　　else　　　　　　FI(i)=exp(h*A(i)/B(i))?！　or i=1:n　　　　if(A(i)==0)。　　m=length(WIJ(:,1))?！　lot(t,y,39?！　lot(t,y,39。故取=25。10 1 1 0]。　　　　y=[y,C*X]。　　y=0。A=a(2:n+1)?！　　！　〗獯穑焊鶕?jù)等比數(shù)列求和方法，在利用matlab中的m文件，編寫程序求解。 　　解答：已知開環(huán)傳遞函數(shù)，求得閉環(huán)傳遞函數(shù)為　　　　　　綜上：當(dāng)閉環(huán)傳函形如時(shí)，可控標(biāo)準(zhǔn)型為：　　?！　　　　　= 1　　　　　　　　　顯然誤差與h2為同階無窮小，歐拉法具有一階計(jì)算精度，精度較低，但算法簡單。該方法建立的數(shù)學(xué)模型受數(shù)據(jù)量不充分，數(shù)據(jù)精度不一致，數(shù)據(jù)處理方法的不完善，很難在精度上達(dá)到更高的要求。傳遞函數(shù)是零極點(diǎn)形式和部分分式形式的基礎(chǔ)。 　　離散事件系統(tǒng)是系統(tǒng)中狀態(tài)變量的改變是由離散時(shí)刻上所發(fā)生的事件所驅(qū)動(dòng)的系統(tǒng)?！　。?）可以靈活設(shè)置仿真試驗(yàn)的時(shí)間標(biāo)尺，既可以進(jìn)行實(shí)時(shí)仿真，也可以進(jìn)行非實(shí)時(shí)仿真。由于受到理論的不完善性以及對事物認(rèn)識的不全面性等因素的影響，其應(yīng)用往往有很大局限性。　　47. 全狀態(tài)反饋控制：是以系統(tǒng)的所有狀態(tài)變量為反饋?zhàn)兞康目刂菩问??！　?1. 超前/滯后校正：如果校正部分的傳遞函數(shù)的極點(diǎn)小于零點(diǎn)，則該校正部分在高于某一頻率范圍時(shí)會給系統(tǒng)一個(gè)正相移，稱作相位超前校正；如果校正部分的傳遞函數(shù)的極點(diǎn)大于零點(diǎn)，則該校正部分會給系統(tǒng)一個(gè)負(fù)相移，稱作相位滯后校正?！　?4. FFT：即快速傅氏變換，是離散傅氏變換的快速算法，它是根據(jù)離散傅氏變換的奇、偶、虛、實(shí)等特性，對離散傅立葉變換的算法進(jìn)行改進(jìn)獲得的。　　28. 連續(xù)系統(tǒng)：系統(tǒng)的各部分的輸入和輸出信號都是連續(xù)函數(shù)的模擬量。　　20. 微分方程：表示未知函數(shù)的導(dǎo)數(shù)以及自變量之間關(guān)系的方程，就叫做微分方程?！　?2. 實(shí)體：就是存在于系統(tǒng)中的具有實(shí)際意義的物體。　　7. 數(shù)值計(jì)算：有效使用數(shù)字計(jì)算機(jī)求數(shù)學(xué)問題近似解的方法與過程。中國地質(zhì)大學(xué)（北京）繼續(xù)教育學(xué)院 2017年03課程考試詞 匯 表　　1. 解析法：就是運(yùn)用已經(jīng)掌握的理論知識對控制系統(tǒng)進(jìn)行理論上的分析、計(jì)算。數(shù)值計(jì)算主要研究如何利用計(jì)算機(jī)更好的解決各種數(shù)學(xué)問題，包括連續(xù)系統(tǒng)離散化和離散形方程的求解，并考慮誤差、收斂性和穩(wěn)定性等問題。　　13. 屬性：就是實(shí)體所具有的任何有效特征?！　?1. 紊流：當(dāng)流體的雷諾系數(shù)大于2000時(shí)，流體的流態(tài)叫紊流?！　?9. 離散系統(tǒng)：系統(tǒng)中一處或多處的信號以脈沖列或者數(shù)碼的形式傳遞的系統(tǒng)?！　?5. Z變換：離散控制系統(tǒng)的拉普拉斯變換?！　?2. 目標(biāo)函數(shù)：目標(biāo)函數(shù)是衡量控制系統(tǒng)控制性能的評判工具，它的大小就說明了控制性能的優(yōu)劣?！　?8. 輸出反饋控制：是以系統(tǒng)的輸出作為反饋?zhàn)兞康目刂菩问健?　　仿真法基于相似原理，是在模型上所進(jìn)行的系統(tǒng)性能分析與研究的實(shí)驗(yàn)方法。　　易于和實(shí)物相連。這種系統(tǒng)的輸入輸出是隨機(jī)發(fā)生的，一般采用概率模型來描述。零極點(diǎn)增益形式可用于分析系統(tǒng)的穩(wěn)定性和快速性。混合法是上述兩種方法的結(jié)合。　　問題17：用二階龍格庫塔法求解題16的數(shù)值解，并于歐拉法求得的結(jié)果比較。　　所以可控標(biāo)準(zhǔn)型是　　問題21：用matlab語言編制單變量系統(tǒng)三階龍格庫塔法求解程序，程序入口要求能接收狀態(tài)方程各系數(shù)陣（A,B,C,D）,和輸入階躍函數(shù)r(t)=R*1(t)?！　文件為　n=64。在matlab命令框中輸入　　+019　　問題26：求出y=x*sin(x) 在條件下的每個(gè)峰值 　　解答：　　問題27：繪制下面的圖形：　　解答：　　　　問題28： 　　解答：　　問題29：已知元件的實(shí)驗(yàn)數(shù)據(jù)如下，擬合這一數(shù)據(jù)，并嘗試給出其特性方程。　　A=[rot90(rot90(eye(n1,n)))。t=T0?！　　　=[t,t(i)+h]?！　　　　　　　?) 按各環(huán)節(jié)相對位置和聯(lián)接關(guān)系，有聯(lián)接矩陣如下：　　　　，所以非零元素矩陣 　　 WIJ=[1 0 1?！　=。r39。39?！　0=zeros(n,1)?！　　　　　I(i)=1?！　　　　　IM(i)=(1FI(i))*C(i)/A(i)。X=Y?！　　　　　or i=1:n　　　　　　　　if(Z(i)~=0)　　　　　　　　　　if(Z(i)==1)　　　　　　　　　　Uk(i)=satu(Uk(i),S(i))。.*Udot。 　　　　　　　　end　　　　　　end　　　　end　　end　　y=[y,Y(nout)]?！　≡贛ATLAB命令窗口中輸入下列語句：　　 KP=。Tf=10。M=KP*TD/T?！　?仿真迭代運(yùn)算：　　for K1=1:Tf/T　　　　ek=RY(1)?！　nd　　%輸出波形： 　　plot(t,yt)　　運(yùn)行結(jié)果為:　　　　此題可以用SIMULINK仿真進(jìn)行驗(yàn)證　　建立SIMULINK仿真模型:　　　　運(yùn)行結(jié)果為:　　　　問題38：設(shè)控制系統(tǒng)的開環(huán)傳遞函數(shù)為試畫出該系統(tǒng)的根軌跡?！　〗獯穑涸贛ATLAB命令窗口中輸入下列命令：　　num=1?！　=logspace(1,1,100)。　　ylabel(39。)?！　en=poly([0,1,5])。　　subplot(2,1,2)。)?！　en=poly([0 10 50 100 200])。)。Phase(deg)39?！　∽詈髞泶_定的取值。校正后閉環(huán)系統(tǒng)的階躍響應(yīng)如下圖所示，可