第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 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) %求稳态值 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 function [pos,tr,ts,tp]=stepchar(y,t) %find pos and tp [mp,ind]=max(y); %返回阶跃响应输出y列向量的最大值mp及对应下标值ind dimt=length(t); %求取时间向量的长度dimt yss=y(dimt); %确定最终的稳定值yss pos=100*(mp-yss)/yss;tp=t(ind); % find rise time tr i=1;j=1;k=1;q=1; while y(i)0.1 %确定输出为0.1时的时刻 i=i+1; end t1=t(i); while y(j)0.9 %确定输出为0.9时的时刻 j=j+1; end t2=t(j);tr=t2-t1; 举例： wn=6;kosi=[0.1:0.1:1.0]; %在对数空间上生成从10^(-1)到10^1共100个数据的横坐标 w=logspace(-1,1,100);num=wn^2; for kos=kosi den=[1 2*kos*wn wn^2]; [mag,pha,w]=bode(num,den,w); % 注意mag的单位不是分贝，若需要分贝表示，可以通过20*log10(mag)进行转换 subplot(211); hold on; semilogx(w,mag) % 注意在所绘制的图形窗口会发现x轴并没有取对数分度 subplot(212); hold on; semilogx(w,pha) end subplot(211);grid on title(bode plot);xlabel(frequency(w));ylabel(amplitude) text(6.2,5,kosi=0.1);text(2,0.5,kosi=1.0) subplot(212);grid on xlabel(frequency(w));ylabel(phase deg) text(5,-20,kosi=0.1);text(2,-85,kosi=1.0) num=1; den=[0.5 1 0]; rlocus(num,den); text(0.1,0.1,k=0); t