ex513.m

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

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