test1.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 9 行

G=tf([1],[10 1]);set(G,'inputdelay',2);
syms s nums dens denm numm den num w;m=0.221;
G=tf(G);num=poly2sym(G.num{1},s);den=poly2sym(G.den{1},s);
[MAG,PHASE,ww] = bode(G);DT=get(G,'inputdelay');
wM=max(ww);wm=min(ww);w=logspace(log10(wm),log10(wM),100);
numm=subs(num,s,(j-m)*w);denm=subs(den,s,(j-m)*w);
Gg=numm./denm.*exp(m*w*DT);
Gm1=Gg./(1+Gg);sysm=frd(Gm1,w);set(sysm,'inputdelay',DT);
[Gmm1,Pmm,Wcgm,Wcpm]=margin(sysm);k=Gmm1