lag1.m

来自「matlab在控制系统中的应用（内附使用方法及每章内容、程序）」· M 代码 · 共 8 行

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];