dmystep.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 106 行

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(dfbnum);
drr2=size(dfbden);
dfinderr=0;
%判断是否为比例环节
if ((drr1(1,2)==1)&(drr2(1,2)==1)),
   dfinderr=1;
   dsysb=tf(dfbnum,dfbden,dhT);
   ddcg=dcgain(dsysb);
   subplot(1,1,1);
   step(dsysb);
   grid on;hold on
   strtext=['这是比例环节，系统指标: 稳态值=' num2str(ddcg) ' ; tr,tp,Mp,ts 无意义 '];
   title(strtext);
else,
n=length(da);
for i=1:n
   de(i,i)=1;
end  
dsysb=dss(da,db,dc,dd,de,dhT);       
ddpb=pole(dsysb);
dii=find(abs(ddpb)>=1);dn1=length(dii);
if (dn1>0)
ddcg=dcgain(dsysb);
[dyc 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');
%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