【文章內(nèi)容簡介】 4,8,4。3,34,3。3,6,8]。b=[4。6。2]。X=linsolve(A,b)% 調(diào)用 linsolve 函數(shù)求解 Ab % 用另一種方法求解 diff(f)%未指定求導(dǎo)變量和階數(shù)，按缺省規(guī)則處理 ans = 6 ． 對方程組求解 a11*x1+a12*x2+a13*x3=b1 a21*x1+a22*x2+a23*x3=b2 a31*x1+a32*x2+a33*x3=b3 syms a11 a12 a13 a21 a22 a23 a31 a32 a33 b1 b2 b3。A=[a11,a12,a13。a21,a22,a23。a31,a32,a33]。b=[b1。b2。b3]。XX=Ab %用左除運算求解(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。f=sqrt(1+exp(x))。%求 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)。diff(f,x,2)diff(f,x,3)f1=a*cos(t)。f2=b*si n(t)。diff(f2)/diff(f1)diff(f)%未指定求導(dǎo)變量和階數(shù)，按缺省規(guī)則處理仿真圖: 波形圖: 其中: G i(R實驗五 MATLAB 在控制系統(tǒng)分析中的應(yīng)用 實驗任務(wù) MATLAB 在控制系統(tǒng)時間響應(yīng)分析中的應(yīng)用。 MATLAB 在系統(tǒng)根軌跡分析中的應(yīng)用； MATLAB 控制系統(tǒng)頻率分析中的應(yīng)用； MATLAB 在控制系統(tǒng)穩(wěn)定性分析中的應(yīng)用 基本命令 itial nd 程序舉例 G(s)t num=[4]。den=[1 , 1 , 4]。step(num , den)[y , x , t]=step(num , den)。tp=spli ne(y , t , max(y))% 計算峰值時間 max(y)% 計算峰值 tp = Step Response Time(sec)p m ans = X i 0 1 X i 0 u X 2 6 5 X 2 1 a=[0,1。6,5]。b=[0。1]。c=[1,0]。d=0。[y,x]=step(a,b,c,d)。Plot(y)1,0 X i X 2 : G(s)兒 num=[4]。den=[1 , 1 ,4]。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]。xO=[1 ,0] ； subplot(1 , 2,1)。in itial(a , b , c ,d,x0)subplot(1 , 2,2)。impulse(a , b , c , d):系統(tǒng)傳遞函數(shù)為: G(s)占 輸入正弦信號時，觀察輸出 信號的相位差。 num=[1]。den=[1 ,1]。t=0 : : 10。u=s in(2*t)。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]。den=[1 2 3]。t=(0:.1:10)。period=4。u=(rem(t,period)=period./2)。% 看 rem 函數(shù)功能 lsim(nu m,de n,u,t)。，繪制系統(tǒng)的根軌跡，并分析其穩(wěn)定性G(s)k(s 2)(s 24s 3)2num=[1 2]。den 仁 [1 4 3]。den=conv(de n1,de n1)。figure(1)6 rlocus(num,den)[k,p]= rlocfi nd(nu m,de n)l L l ■ 匚一 L — Cr t r r108642 4 6 Root Locus 0 2 2 20 0 2 2CPXA vyanma n8 G(s)Lin ear Simulati on e Time(sec)p mimpulse response(k=55)5figure(2)k=55。num 仁 k*[1 2]。den=[1 4 3]。den 1=c onv(de n,den)。[nu m,de n]=cloop(nu m1,de n1,1)。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。num 仁 k*[1 2]。den=[1 4 3]。den 1=c onv(de n,den)。[nu m,de n]=cloop(nu m1,de n1,1)。impulse(nu m,de n)2468 L0 500 1000 Time(sec)title(”impulse resp on se(k=56)“)1500 2000 2500Select a point in the graphics win dow selected_po int = + Bode Diagram bode 圖 G(s) n=[1 , 1]。d=[1 , 4,11 , 7]。bode(n , d),grid on Frequency(rad/sec)s 1 s 34s 211s 7 G(s)s 1 e(s 2)3 求有理傳函的頻率響應(yīng)，然后在同一張圖上繪出以四階伯德近似表示的系統(tǒng)頻 率響應(yīng) num=[1]。de n=co nv([1 2],co nv([1 2],[1 2]))。w=logspace(1,2)。t=。ylabel(”gai n“)。subplot(2,1,2)。semilogx(w,p1,w,p2,”g“)。grid on。xlabel(”freque ncy“)。ylabel(”phase“)。求它的幅值裕度和相角裕度 n=[]。d=[1 2 3 2]。[Gm,Pm,Wcg,Wcp]=margi n(n,d)G(s) s 32s 23s 2 bode plot [m1,p1]=bode(nu m,de n, 2)。p1=p1t*w”*180/pi。[n 2,d2]=pade(t,4)。nu mt=c onv(n2,nu m)。den t=(c onv(de n, d2))。2 freque