第8讲 系统分析32.pptVIP

第八讲 系统的分析 ;;num=[3 16 41 28];den=[1 14 110 528 1494 2117 112]; %系统描述 [z,p,k]=tf2zp(num,den) %求系统的零极点 ii=find(real(z)0) n1=length(ii); %检验零点的实部；求取零点实部大于零的个数 jj=find(real(p)0) n2=length(jj); %检验极点的实部；求取极点实部大于零的个数 if(n20) disp(the system is unstable) disp(the unstable pole are:) disp(p(jj)) else disp(the system is stable) end %判断系统是否稳定 if(n10) disp(the system is a nonminimal phase one) else disp(the syetem is a minimal phase one) end %判断系统是否为最小相位系统 pzmap(p,z) %绘制零极点图 p z;运行结果为： the system is stable the syetem is a minimal phase one p = -1.9474 + 5.0282i -1.9474 - 5.0282i -4.2998 -2.8752 + 2.8324i -2.8752 - 2.8324i -0.0550 z = -2.1667 + 2.1538i -2.1667 - 2.1538i -1.0000 ;8.2 系统的时域分析;num=[20];den=[1 8 36 40 0]; %开环传递函数描述 [numc,denc]=cloop(num,den); %求闭环传递函数 t=0:0.1:10; y=step(numc,denc,t); [y1,x,t1]=step(numc,denc); %绘制闭环系统的阶跃响应曲线 %对于传递函数调用，状态变量x返回为空矩阵 plot(t,y,r:,t1,y1) title(‘the step responce’)；xlabel(time-sec) disp(系统稳态值dc为：) dc=dcgain(numc,denc) %求稳态值;运行结果为： 系统稳态值dc为： dc = 1 ;2、impulse()函数的用法;wn=5;alfh=0.6;num=wn^2; den=[1 2*alfh*wn wn^2]; %系统传递函数描述 t=0:0.02:5; y=step(num,den,t); plot(t,y) %绘制闭环系统的阶跃响应曲线 title(two orders linear system step responce) xlabel(time-sec);ylabel(y(t));grid on;num=[10 25];den=[0.16 1.96 10 25]; %系统传递函数描述 t=0:0.02:5;y=step(num,den,t); plot(t,y) %绘制闭环系统的阶跃响应曲线 xlabel(time-sec); ylabel(y(t)) grid;二、常用时域分析函数;%输入典型二阶系统参数，确定系统传递函数模型 alph=input(please input alph=);wn=input(please input wn=); num=wn^2;den=[1 2*alph*wn wn^2]; [z,p,k]=tf2zp(num,den); ii=find(real(z)0);n1=length(ii); jj=find(real(p)0);n2=length(jj); if(n20) disp(the system is unstable) disp(it is no use for getting 动态参数) %判断系统是否稳定 else [y,x,t]=step(num,den); plot(t,y) [pos,tp,tr,ts]=stepchar(y,t) %调用求取二阶系统阶跃响应动态参数的函数文件 end;运行结果： please input alph=0.3 please input wn=5 pos = 37.5896 tp = 0.2577 tr = 2.2454 ts = 0.6626 ;function [pos,tr,ts,tp]=st