%   if dbodeabcd==1
dwindbodek=figure('color',[.8,.8,.9],'position',[0.1*winw,0.1*winh,0.8*winw,0.8*winh],...
       'name','離散系統Z域開環Bode圖分析','numbertitle','off','menubar','figure',...
       'resize','off');
% end
% if bodeabcd==0;
%    figure(windbode)
% end
% bodeabcd=0;
dnk=length(dak);
for i=1:dnk
   dek(i,i)=1;
end  
dsysk=dss(dak,dbk,dck,ddk,dek,dhT); 
grid on
margin(dsysk)
uicontrol(dwindbodek,'style','push','string','結束',...
       'position',[0.83,0.01,0.11,0.05],'units','normalized',...
       'callback','dclosebodek');
%   abcdefg=0;
 %  mfksys=tf(knum,kden);
 %  ltiview('bode',mfksys)