finddataerr=0;
dktnum=get(editdknum,'string');
dktden=get(editdkden,'string');
dktanum=eval(dktnum);dktaden=eval(dktden);
dhfT=get(edittime,'string');dhT=eval(dhfT);
[err1 err2]=size(dktanum);
[err3 err4]=size(dktaden);
if err2>err4,
  finddataerr=1;
else,
  if dfadeback==0,                                                       %負反饋
    dhfnum=get(editdhknum,'string');
    dhfden=get(editdhkden,'string');
    dhffnum=eval(dhfnum);dhffden=eval(dhfden);
    [dfbnum,dfbden]=feedback(dktanum,dktaden,dhffnum,dhffden,-1);
    dfknum=conv(dktanum,dhffnum),dfkden=conv(dktaden,dhffden)
  end
  if dfadeback==1,                                                       %正反饋
    dhfnum=get(editdhknum,'string');
    dhfden=get(editdhkden,'string');
    dhffnum=eval(dhfnum);dhffden=eval(dhfden);
    [dfbnum,dfbden]=feedback(dktanum,dktaden,dhffnum,dhffden,1);
    dfknum=conv(dktanum,dhffnum);dfkden=conv(dktaden,dhffden);
  end
  if dfadeback==2                                                        %無反饋
    dhfnum=get(editdhknum,'string');
    dhfden=get(editdhkden,'string');
    dhffnum=eval(dhfnum);dhffden=eval(dhfden);
    [dfbnum,dfbden]=feedback(dktanum,dktaden,dhffnum,dhffden,-1);
    dfknum=conv(dktanum,dhffnum),dfkden=conv(dktaden,dhffden);
  end
end   
[rr1 rr2]=size(dfbnum);
for i=1:1:rr2
   if dfbnum(1,i)==0,
      if dfbden(1,i)==0,
      else,
         finddataerr=0;
         break;
      end
   else,
      if dfbden(1,i)==0,
        finddataerr=1;
        break;
      end
   end
end
if finddataerr==0,
   
   [dak,dbk,dck,ddk]=tf2ss(dfknum,dfkden);
  dtak=dak;dtbk=dbk;dtck=dck;dtdk=ddk;
  [dak1,dbk1,dck1,ddk1,dtk1]=canon(dak,dbk,dck,ddk,'model');
  [dak2,dbk2,dck2,dtk2,dkk2]=ctrbf(dak,dbk,dck);
  [dak3,dbk3,dck3,dtk3,dkk3]=obsvf(dak,dbk,dck);
  [dzk,dpk,dkk]=ss2zp(dak,dbk,dck,ddk);
  [dminak,dminbk,dminck,ddmindk]=minreal(dak,dbk,dck,ddk);
  [dminnumk,dmindenk]=minreal(dfknum,dfkden);

  [da,db,dc,dd]=tf2ss(dfbnum,dfbden);
  dta=da;dtb=db;dtc=dc;dtd=dd;
  [da1,db1,dc1,dd1,dt1]=canon(da,db,dc,dd,'model');
  [da2,db2,dc2,dt2,dk2]=ctrbf(da,db,dc);
  [da3,db3,dc3,dt3,dk3]=obsvf(da,db,dc);
  [dz,dp,dk]=ss2zp(da,db,dc,dd);
  [dmina,dminb,dminc,ddmind]=minreal(da,db,dc,dd);
  [dminnum,dminden]=minreal(dfbnum,dfbden); 
   
   
   
   [da,db,dc,dd]=tf2ss(dfbnum,dfbden);
  dta=da;dtb=db;dtc=dc;dtd=dd;
  [da1,db1,dc1,dd1,dt1]=canon(da,db,dc,dd,'model');
  [d2,db2,dc2,dt2,dk2]=ctrbf(da,db,dc);
  [da3,db3,dc3,dt3,dk3]=obsvf(da,db,dc);
  [dz,dp,dk]=ss2zp(da,db,dc,dd);
  [dmina,dminb,dminc,ddmind]=minreal(da,db,dc,dd);
  [dminnum,dminden]=minreal(dfbnum,dfbden);
  test=0;
  close(changemodel3);
  set(mmodel,'enable','on');
  set(manalysis,'enable','on');
  set(msave,'enable','on');
  set(mdesign,'enable','on');
  
  
  
  
else,
  msgbox('所構造的傳遞函數分子維數大于分母維數','數據錯誤警告','help');
end
set(idwin,'enable','off');
set(idwins1,'enable','off');
set(dsmys,'enable','off');
set(msimu,'enable','off');
set(dmsimu,'enable','on');
set(dmchange,'enable','on');
set(dminfy,'enable','on');
set(mdtime,'enable','on');
set(mdcplx,'enable','on');
set(idwin,'enable','off');
set(idwins1,'enable','off');
 set(mdesign1,'enable','off');
 set(mdesign2,'enable','off');
 set(mdesign,'enable','off');