【正文】 定性. G=tf(1000,co nv([1,3,2],[1,5]))。d4=[1,1,1]。d 仁 [1 , 4,1]。Nyquist Diagram 80 |[[: G(s)n s 2n S S 令 wn=1, 分別作出 E =2,1 , , C A 5 ya 時的 nyquist 曲線。nyq uist(n, d3)。2 freque xlabel(“frequeGm = Pm = Wcg = Wcp =nyq uist(n, d1)。nu mt=c onv(n2,nu m)。p1=p1t*w”*180/pi。d=[1 2 3 2]。ylabel(”phase“)。grid on。subplot(2,1,2)。t=。de n=co nv([1 2],co nv([1 2],[1 2]))。d=[1 , 4,11 , 7]。[nu m,de n]=cloop(nu m1,de n1,1)。den=[1 4 3]。impulse(nu m,de n)title(”impulse resp onse(k=55)“) 0 200 400 600 Time(sec)800 1000 1200 5 50 050 0ox 10 8!impulse resp on se(k=56)figure(3)k=56。den 1=c onv(de n,den)。num 仁 k*[1 2]。den=conv(de n1,de n1)。繪制系統(tǒng)的根軌跡，并分析其穩(wěn)定性G(s)k(s 2)(s 24s 3)2num=[1 2]。u=(rem(t,period)=period./2)。t=(0:.1:10)。hold on plot(t,u, ”r“)lsim(nu m,de n,u,t)Real Axis，求出周期為 4 秒的方波的輸出響應(yīng) 2s 25s 1 ~2 s 2s 3 num=[2 5 1]。t=0 : : 10。 num=[1]。in itial(a , b , c ,d,x0)subplot(1 , 2,2)。impulse(nu m,de n)Time(sec)Response to Initial Conditions : e u p m A L I 1r廠、/ Impulse Responsec=[1 , 0] ；d=[0]。Plot(y)1,0 X i X 2 : G(s)兒 num=[4]。d=0。1]。6,5]。step(num , den)[y , x , t]=step(num , den)。 MATLAB 在系統(tǒng)根軌跡分析中的應(yīng)用； MATLAB 控制系統(tǒng)頻率分析中的應(yīng)用； MATLAB 在控制系統(tǒng)穩(wěn)定性分析中的應(yīng)用 基本命令 itial nd 程序舉例 G(s)t num=[4]。f2=b*si n(t)。%求 f 對 x 的二階導(dǎo)數(shù) %求 f 對 x 的三階導(dǎo)數(shù) %按參數(shù)方程求導(dǎo)公式求 y 對 x 的導(dǎo)數(shù) ans = 1/2/(1+exp(x))A(1/2)*exp(x)ans =2*si n(x)x*cos(x)ans =3*cos(x)+x*si n(x)ans =b*cos(t)/a/si n(t)三、SIMULINK 的使用f=x*cos(x)。XX=Ab %用左除運(yùn)算求解(X=linsolve(A,b)%調(diào)用 linsolve 函數(shù)求的解)XX =(a12*a23*b3a12*b2*a33+a13*a32*b2a13*a22*b3+b1*a22*a33b1*a32*a23)/(a11*a22*a33a11*a32*a23a12*a21*a33+a32*a21*a13a22*a31*a13+a31*a12*a 23)(a11*a23*b3a11*b2*a33a21*a13*b3a23*a31*b1+b2*a31*a13+a21*b1*a33)/(a11*a22*a33a11*a32*a23a12*a21*a33+a32*a21*a13a22*a31*a13+a31*a12* a23)(a32*a21*b1a11*a32*b2+a11*a22*b3a22*a31*b1a12*a21*b3+a31*a12*b2)/(a11*a22*a33a11*a32*a23a12*a21*a33+a32*a21*a13a22*a31*a13+a31*a12*a 23)7 ． syms a b t x y z。b2。a31,a32,a33]。A=[a11,a12,a13。2]。b=[4。3,34,3。s=(7*x A 28*y A 2)*(x A 2+3*y A 2)。%定義 4 個符號變量 %定義 4 個數(shù)值變量 %建立符號矩陣 A %建立數(shù)值矩陣 B %計算符號矩陣 A 的行列式 %計算數(shù)值矩陣 B 的行列式[ c, d]X108 11 ans = a*db*c ans = 150 4.187。y,z] det(A)det(B)A = [ a, b] ”c“)。A=[a,b。y=8。c=sym(w=10。da=det(a)ea=eig(a)da = 7 ea = 7/2+1/2*21A(1/2)7/21/2*21A(1/2)(包括精確解和一定精度的解) r1=solve(”xA2+x1“)rv=vpa(r1)rv4=vpa(r1,4)rv20=vpa(r1,20)r1 = 1/2*5 八(1/2)1/21/2*5 八(1/2)1/2 rv =. rv4 =. rv20 =. 3 ． a=sym(”a“)。da=det(a)ea=eig(a)da = a11*a22a12*a21 ea = 1/2*a11+1/2*a22+1/2*(a11 A 22*a11*a22+a22 A 2+4*a12*a21)A(1 /2)1/2*a11+1/2*a221/2*(a11 A 22*a11*a22+a22 A 2+4*a12*a21)A(1 /2)a=sym(”[2 3。程序舉例 a=sym(”[a11 a12。disp([num2str(num3),”=“,num2str(num1),”x“,num2str(num2)])End End Endend y=input(” 是否繼續(xù)判斷？繼續(xù)請按 1 ,按任意鍵退出 ：“)； if y~=1 break end end 實驗四 MATLAB 的符號計算與 SIMULINK 的使用 實驗任務(wù) MATLAB 符號計算的特點和常用基本命令。num1=n。function sushu while 1 x=input（” 請輸入一個自然數(shù) “）。例：輸入 6，因為 6 不是素數(shù)。如果操作者希望停止這種判斷，則可以退出程序。最好能具有如下功能：①設(shè)計較好的人機(jī)對話界面，程序中 含有提示性的輸入輸出語句。pp =3 e2(500)ans =1 2 3 5 8 13 21 34 55 89 144 233377 lookfor ffibno e2ffibno 計算斐波那契亞數(shù)列的函數(shù)文件 help e2 ffibno 計算斐波那契亞數(shù)列的函數(shù)文件 n 可取任意自然數(shù) 程序如下（用法：lookfor 關(guān)鍵詞 在所有 M 文件中找 “關(guān)鍵詞 ”，比如：lookfor max（即尋找關(guān)鍵詞 “ max”）其實就和我們平時用 CTRL+F 來查找關(guān)鍵詞”是一樣的 而 help 是顯示 matlab 內(nèi)置的幫助信息 用法：help 命令，比如 help inv，作用就是調(diào)用 inv 這個命令的幫助）程序設(shè)計題 用一個 MATLAB 語言編寫一個程序：輸入一個自然數(shù)，判斷它是否是素數(shù)，如果是，輸出“ It is one prime ”，如果不是，輸出“ It is n ot o ne pr