資源描述:
《MATLAB經(jīng)典實(shí)例》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、MATLAB一、用MATLAB求解線性微分方程和化簡(jiǎn)系統(tǒng)方框圖1、MATLAB中數(shù)學(xué)模型的表示傳遞函數(shù):P=[1024];>>A=[1,3];>>B=[10,20,3];>>C=conv(A,B)/*C=(A*B)*/C=1050639W(s)=4(s+2)(s+3)(s+4)>>W=4*conv([1,2],conv([1,3],[1,4]))W=43610496Sys=tf(num,den)>>num=[1,2,3,4];>>den=[1,0,1];>>W=tf(num,den)Transferfunction:s^3+2s^2+3s+4---------------------s
2、^2+1傳遞函數(shù)的特征根:roots(P)>>p=[1304];>>r=roots(p)>>r=-3.35530.1777+1.0773i0.1777-1.0773i>>poly(r)>>ans=1.00003.00000.00004.0000polyval(p,a)>>n=conv([321],[14]);>>value=polyval(n,-5)>>value=-66線性系統(tǒng)的零極點(diǎn)模型:>>z=[-1-2-3];>>p=[-2-5];>>W=zpk(z,p,k)Zero/pole/gain:k(s+1)(s+2)(s+3)-------------------(s+2)(s+5)
3、2、用MATLAB求解線性微分方程[r,p,k]=residue(num,den)>>num=[2161024];>>den=[110355024];>>[r,p,k]=residue(num,den)r=-88.000091.5000-26.00003.5000p=-4.0000-3.0000-2.0000-1.0000k=2>>num=1;>>den=[conv(conv([1-1],[14]),[2-31])];>>[r,p,k]=residue(num,[den0])(R(s)=1/s)r=0.0011-0.64000.20000.8889-0.2500p=-4.00001.0
4、0001.00000.50000k=[]>>num=1;>>den=[conv(conv([1-1],[14]),[2-31])];>>[r,p,k]=residue(num,[den00])(R(s)=1/s^2)r=-0.0003-0.84000.20001.7778-0.9375-0.2500p=-4.00001.00001.00000.500000k=[]3、MATLAB在系統(tǒng)方框圖化簡(jiǎn)中的應(yīng)用串聯(lián)連接結(jié)構(gòu):[num,den]=series(num1,den1,num2,den2)[并聯(lián)連接結(jié)構(gòu):[num,den]=parallel(num1,den1,num2,den2)反饋
5、連接結(jié)構(gòu):[num,den]=feedback(num1,den1,num2,den2,sign)num1=[11];>>den1=[12];>>num2=[1];>>den2=[10];>>[num,den]=feedback(num1,den1,num2,den2,-1)num=110den=131>>printsys(num,den)num/den=s^2+s-------------s^2+3s+1二、用MATLAB進(jìn)行系統(tǒng)時(shí)域分析1、典型輸入信號(hào)的MATLAB實(shí)現(xiàn)單位脈沖響應(yīng):[y,x,t]=impulse(num,den,t)或impulse(num,den)>>t=[0:
6、0.1:10];>>num=[1];>>den=[111];>>[y,x,t]=impulse(num,den,t);>>plot(t,y);>>grid;(%繪制網(wǎng)格)>>xlabel('t');>>ylabel('y');>>title('單位脈沖響應(yīng)')單位階躍響應(yīng):[y,x,t]=step(num,den,t)或step(num,den)>>t=[0:0.1:10];>>num=[1];>>den=[10.51];>>[y,x,t]=step(num,den,t);>>plot(t,y);>>grid;(%繪制網(wǎng)格)>>xlabel('t');>>ylabel('y');>>t
7、itle('單位階躍響應(yīng)')斜坡響應(yīng):t=[0:0.1:10];num=[1];den=[10.310];c=step(num,den,t);plot(t,c);grid;xlabel('t');ylabel('y');title('單位斜坡響應(yīng)')任意函數(shù)作用下系統(tǒng)的響應(yīng):[y,x]=lsim(num,den,u,t)num=[10,20];den=[102020];%產(chǎn)生三角波函數(shù)v1=[0:0.1:2];v2=[1.9:-0.1:-2];