MATLAB经典实例

《MATLAB经典实例》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、MATLAB一、用MATLAB求解线性微分方程和化简系统方框图1、MATLAB中数学模型的表示传递函数：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传递函数的特征根：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线性系统的零极点模型：>>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在系统方框图化简中的应用串联连接结构：[num,den]=series(num1,den1,num2,den2)[并联连接结构：[num,den]=parallel(num1,den1,num2,den2)反馈

5、连接结构：[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进行系统时域分析1、典型输入信号的MATLAB实现单位脉冲响应：[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;(%绘制网格)>>xlabel('t');>>ylabel('y');>>title('单位脉冲响应')单位阶跃响应：[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;(%绘制网格)>>xlabel('t');>>ylabel('y');>>t

7、itle('单位阶跃响应')斜坡响应：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,x]=lsim(num,den,u,t)num=[10,20];den=[102020];%产生三角波函数v1=[0:0.1:2];v2=[1.9:-0.1:-2];