matlab控制系統(tǒng)仿真(已改無錯字)
2022-09-01 23:07:19 本頁面
【正文】 ck()函數(shù)的應(yīng)用 例 一個多環(huán)的反饋系統(tǒng)如圖 249所示 ， 給定各環(huán)節(jié)的傳 遞函數(shù)為： 11()10Gs s? ? 21()1Gs s? ?23 21()44sGsss???? 4 1() 6sGs s ?? ?1 ( ) 1Hs ? 2 ( ) 2Hs ?31()2sHss???試求閉環(huán)傳遞函數(shù) GB(s)=C(s)/R(s)。 解 求解步驟如下： 步驟 1：輸入系統(tǒng)各環(huán)節(jié)的傳遞函數(shù)； 步驟 2：將 H2的綜合點移至 G2后； 步驟 3：消去 G3， G2， H2環(huán)； 步驟 4：消去包含 H3的環(huán)； 步驟 5：消去其余的環(huán) ， 計算 GB (s)。 根據(jù)上述步驟的 MATLAB指令以及計算結(jié)果在圖 272中 。 ng1=[1]。dg1=[1 10]。 ng2=[1]。dg2=[1 1]。 ng3=[1 0 1]。dg3=[1 4 4]。 ng4=[1 1]。dg4=[1 6]。 nh1=[1]。dh1=[1]。 nh2=[2]。dh2=[1]。 nh3=[1 1]。dh3=[1 2]。 [n1,d1]=series(ng2,dg2,nh2,dh2)。 [n2,d2]=feedback(ng3,dg3,n1,d1,1)。 [n3,d3]=series(n2,d2,ng4,dg4)。 [n4,d4]=feedback(n3,d3,nh3,dh3,1)。 [n5,d5]=series(ng1,dg1,ng2,dg2)。 [n6,d6]=series(n5,d5,n4,d4)。 [n7,d7]=cloop(n6,d6,1)。 printsys(n7,d7) num/den= s^4+ 3s^3+ 3s^2+3s+2 2s^6+38s^5+261s^4+1001s^3+1730s^2+1546s+732 圖 272 多環(huán)結(jié)構(gòu)圖簡化 通過 pzmap()或 roots()函數(shù)可查看傳遞函數(shù)是否有相同的零極 點 ， 還可使用 minreal()函數(shù)除去傳遞函數(shù)共同的零極點因子 。 如圖 273所示 。 numg=[1 6 11 6]。deng=[1 7 12 11 5]。 printsys(numg,deng) numg/deng= s^3+6s^2+11s+6 s^4+7s^3+12s^2+11s+5 [num,den]=minreal(numg,deng)。 printsys(num,den) 1 polezeros cancelled num/den= s^2+4s+3 s^3+6s^2+6s+5 圖 273 minreal()函數(shù)的應(yīng)用 例 ，在給定各元件參數(shù)并忽略 La和令 ML = 0的情況下，其結(jié)構(gòu)圖如圖 274所示： 圖 274 位置隨動系統(tǒng)的結(jié)構(gòu)圖 第一步求閉環(huán)傳遞函數(shù) GB (s)=? c(s) /? r(s)，求解過程及結(jié)果 如圖 275所示。第二步利用 step()函數(shù)計算參考輸入 ? r (t)為單位 階躍信號時輸出 ? c (t)的響應(yīng)。 num1=[200]。 den1=[20]。num2=[1]。 den2=[2 0]。 num3=[ 0]。 den3=[1]。 num
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1