bodelag.m

xnsy算法的matlab实现

function [nc,dc]=bodelag(ng,dg,kc,xpm)
[mag,phase,w]=bode(kc*ng,dg);
[gm,pm,wcg,wcp]=margin(mag,phase,w);
wgc=spline(phase,w,xpm+5-180);
mag_db=20*log10(mag);
gr=-spline(w,mag_db,wgc); 
alpha=10^(gr/20);
T=10/(alpha*wgc);
nc=[alpha*T,1];
dc=[T,1];