《基于MATLAB的连续时间系统分析_马俊兴》由会员分享，可在线阅读，更多相关《基于MATLAB的连续时间系统分析_马俊兴（3页珍藏版）》请在金锄头文库上搜索。
1、143河南教育学院学报(自然科学版) Vol.14No.320 05M9Journal of Henan Institute of Education(NaturalScience) Sep.2005l:2005-07-06Te:.(1965),3,2,2q.MATLABHWds马俊兴,刘晓艳(河南教育学院物理系,河南郑州450014)K1:由人工完成信号与系统的高阶系统分析几乎是不可能实现的,本文结合实际的高阶微分方程给出了如何应用MATLAB工程计算语言简便解决这些问题的方法.1oM:信号与系统;MATLAB;系统分析ms|:O411DSM:AcI|:1007-0834(2005)03-0
2、060-030|dsY、1e、_、Se51T.|dsB15,)d&S#|Y.dBVsZsZ,HWdsZ,HWdsZ,sZsZ9Tv,yN,ds+V?LC.kmBASIC、FORTRANCIM19,E1I/F,yNK.SMATHWORKS1967Mw“Matrix laboratory”(MATLAB),rs、!|m)、?QB8,BZLOz,Baq9.1+YMATLAB|)q,VZL|ds9、Vjy#d!9、_k.MATLABSv*X1|d.MATLABC|ds4vZL.L=sZMATLAB9eL%t5ZE/F.1dsBHWds/+:pd,Y;pd,Y;d;d,sm./sZMATLABqsd8Z
3、E.2MATLABdssZ:4 d4ydt4 +3d3ydt3 +2d2ydt2+1=U(t),kpd,Y,Y,isd,d,sm.!S.y (0 )=0, y(0 )=1, y(0 )=0, y(0 )=0.pV/.2.1p,Y(dS)S, y (0 )=0, y(0 )=1, y(0 )=0,y(0 )=0MATLAB:TY0=(0,1, 0,0).5dsZ7Z,1XMAT-LAB,pfinitial(A,B,C,D, x0,t),A5sZpdA,B,C, D.f A,B,C,D=tf2ss(num,den)Vp,num,densYdsZfYf.1p/:B:u7MATLABq,y,/:cle
4、ar, close all%use num anddenfindthe initialmap,firstusenum and denmake A, B, C, D,thencanuse the“initial”function.num=1 ;den=43201 ; A,B, C,D =tf2ss(num,den);x0= 0100 ;t=0:0.1:10;initial(A,B,C, D, x0,t);=:ki,Vd,YH_mm1U.60m1HW30s,Ym:9,AU,YmV,dd,Ff/,?vt,79v.2.2s,Y2.2.1pdYfimpluse(num, den,t)fimpluse(A
5、, B, C, D)2Vp.B:u7MATLABq,y,/:clear all%use impluse(num,den)or(A,B, C,D)num=input(|_num= );den=input(|-_den=);t=input(_HWt= ); y,x =impulse(num,den,t);impulse(num,den, t),grid;=:ki,v74U1pnumden_/:|_num= 10|-_den=432 01_HWt=30:dYH_mm2U.2.2.2pYfstep(num, den, t)step(A, B, C, D),fpBVY_m.3B:u7MATLABq,y,
6、/:clear all%use step(num,den)or(A, B,C,D)num=input(|_num= );den=input(|-_den=);t=input(_HWt= ); y,x =step(num,den, t);step(num,den, t);=:ki,v74U1pnumden_/:|-_den=432 01_HWt=30|_num= 10:dYH_mm3U.m2Y_mm3dY_m2.3dd1dxP5d.pd,VMATLABf z, p,k =tf2zp(num, den)1p,realfdl,.B:u7MATLABq,y,/:clear all;%decide wh
7、ether the systerm issteadynum=input(|_num= );den=input(|-_den=); z, p,k =tf2zp(num, den);display(p);n=find(real(p)0;if(length(n)0)display( The systemisnotsteady.),else display( the system is steady. )end;=:ki,v74Unumden_/:num=161den= 43201:7VAT/.p = -0.6205 +0.6402i-0.6205 -0.6402i0.2455 +0.5042i0.2
8、455 -0.5042iThe system isnot steady.2.4,sm4,smZE,Vp,mf,9VMATLABmap(p, z),BZT.B:u7MATLABq,y,/:clear,close all;%draw the p,zmap.num=input(|_num= );den=input(|-_den=); z, p,k =tf2zp(num, den);display(p);display(z);%draw a ciclezplane(z,p);x=-1:0.001:1;y1=(1-x. 2). 0.5;y2=-(1-x. 2). 0.5;plot(x,y1,x,y2);hold;zplane(z, p);=:7VAT/,iVm
