dmystepbb.m

详细介绍了Mathlab的使用方法与技巧

%if dfadeback==2
%      dmystepbb
%   else
%for i=1:n
%  if db(i)==0
%      db(i);
%   else
%      figure(i+1)
%      dstep(da,db,dc,dd,i);
%   end
%end

 dwindstep=figure('color',[0.8,0.7,0.8],'position',[0.1*winw,0.1*winh,0.8*winw,0.8*winh],...
       'name','开环离散系统输出Y(z)和阶跃响应','numbertitle','off','menubar','figure',...
       'resize','off');
drr1=size(dfknum);
drr2=size(dfkden);
dfinderr=0;
%判断是否为比例环节
if ((drr1(1,2)==1)&(drr2(1,2)==1)),
   dfinderr=1;
   dsysk=tf(dfknum,dfkden,dhT);
   ddcg=dcgain(dsysk);
   subplot(1,1,1);
   step(dsysk);
   grid on;hold on
   strtext=['这是比例环节，系统指标: 稳态值=' num2str(ddcg) ' ; tr,tp,Mp,ts 无意义 '];
   title(strtext);
else,
n=length(dak);
for i=1:n
   dek(i,i)=1;
end  
dsysk=dss(dak,dbk,dck,ddk,dek,dhT);       
ddpk=pole(dsysk);
dii=find(abs(ddpk)>=1);dn1=length(dii);
if (dn1>0)
ddcg=dcgain(dsysk);
[dyc dx]=step(dsysk);
dnsn=length(dyc);dt=1:dnsn;
dyc1=dyc/ddcg;
subplot(2,1,1);
plot(dt,dx,'r');
grid on;hold on
xlabel('s/t');
ylabel('X(z)');
legend('红线:X(z)状态曲线',0);
title('X(z)状态曲线');
subplot(2,1,2);
grid on;hold on
step(dsysb);
 %  plot(t,yc1,'r',t,yc,'b');
grid on;hold on
xlabel('k*T');
ylabel('Y(z)和Y(z)/dcg');
%legend('红线:Y(z)/dcg阶跃响应曲线','蓝线:Y(z)阶跃响应曲线',0);
strtext=[' 不稳定系统，无稳态值且tr,tp,Mp,ts 无意义 '];
title(strtext); 

else
   
ddcg=dcgain(dsysb);
[dyc dty dx]=step(dsysb);
dnsn=length(dyc);dt=1:dnsn;
dyc1=dyc/ddcg;
subplot(2,1,1);
plot(dt,dx,'r');
grid on;hold on
xlabel('s/t');
ylabel('X(z)');
legend('红线:X(z)状态曲线',0);
title('X(z)状态曲线');
subplot(2,1,2);
grid on;hold on
step(dsysb);
 %  plot(t,yc1,'r',t,yc,'b');
grid on;hold on
xlabel('k*T');
ylabel('Y(z)和Y(z)/dcg');
 
   
  dycn=step(dsysb);
  dtn=length(dycn);
  
  [dycnn dty]=step(dsysb);dcin=length(dycnn);
  %[dz,dp]=ss2zp(dsysb);
  dp=pole(dsysb);  
  uicontrol(dwindstep,'style','push','string','稳态值',...
       'position',[0.06,0.01,0.11,0.05],'units','normalized',...
      'callback','dmysteppre_dcg');
uicontrol(dwindstep,'style','push','string','上升时间tr',...
       'position',[0.22,0.01,0.11,0.05],'units','normalized',...
      'callback','dmysteppre_tr');
uicontrol(dwindstep,'style','push','string','峰值时间tp和最大超调量Mp',...
       'position',[0.37,0.01,0.25,0.05],'units','normalized',...
      'callback','dmysteppre_tpmp');
uicontrol(dwindstep,'style','push','string','调整时间ts',...
       'position',[0.66,0.01,0.11,0.05],'units','normalized',...
      'callback','dmysteppre_ts');
uicontrol(dwindstep,'style','push','string','结束',...
       'position',[0.83,0.01,0.11,0.05],'units','normalized',...
       'callback','dclosestep');
   
end
end
%end