lag1.m

《MATLAB在控制系统中的应用》实例程序

function [ngc,dgc,k]=lag1(ng0,dg0,KK,s1,a)
ngv=polyval(ng0,s1); dgv=polyval(dg0,s1);
g=dgv/ngv; k=abs(g); beta=k/KK;
[kosi1,wn1]=s2kw(s1);
zc=-wn1*sin(a*pi/180)/sin(pi-atan(sqrt(1-kosi1^2)/kosi1)-(a*pi/180))
pc=beta*zc;
ngc=beta*[1,-zc];dgc=[1,-pc];