fun1.m

文件是编辑好的m文件

function sol=fun1(G);
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);
Gm=Gg*y./(1+y*Gg);sysm=frd(Gm,w);set(sysm,'inputdelay',DT);[Gmm,Pmm,Wcgm,Wcpm]=margin(sysm);
sol=Gmm-1;