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

				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