dnum=[3.4 -1.5];dden=[1 -1.6 0.8]; Ts=0.05;sysd=tf(dnum,dden,Ts);
num1=0.25;den1=[1 1];num2=1;den2=[1 2];
sys1=tf(num1,den1);sys2=tf(num2,den2);sysc2=sys1*sys2;
Ts=0.05:0.01:100;m=length(Ts);
for i=1:m
    sysd=tf(dnum,dden,Ts(i));syscd=c2d(sysc2,Ts(i),'zoh');     %修改 
    dsys=sysd*syscd;dsysb=feedback(dsys,1);            %組成新的閉環脈沖傳遞函數
    [p,z]=pzmap(dsysb);pp=abs(p);                     %求出系統極點向量的模向量
    n=find(pp>=1);
if length(n)>=1
    disp('系統臨界穩定的采樣周期Ts='),Ts=Ts(i),
    disp('系統臨界穩定的零極點='),p,z,break
end;end