bodelead.m

xnsy算法的matlab实现

function [nc,dc]=bodelead(ng,dg,kc,xpm)
[mag,phase,w]=bode(kc*ng,dg);
[gm,pm,wcg,wcp]=margin(mag,phase,w);
dpm=xpm-pm+6;
smo=length(mag);
phi=dpm*pi/180;
a=(1+sin(phi))/(1-sin(phi));
mag_db=20*log10(mag);mm=-10*log10(a);
wgc=spline(mag_db,w,mm);
T=1/(wgc*sqrt(a));
z=a*T;
p=T;
nc=[z,1];
dc=[p,1];