evans
lugar geométrico das raízes Evans
Seqüência de Chamamento
evans(H) evans(H, kmax)
Parâmetros
- H
lista (sistema linear
syslin
)- kmax
real (ganho máximo desejado para o esboço)
Descrição
Fornece o lugar geométrico das raízes Evans para um sistema linear
em estado-espaço ou forma de transferência H(s) (lista
syslin
). Este é o lugar geométrico das raízes de
1+k*H(s)=1+k*N(s)/D(s)
, no plano dos complexos. Para
uma amostra seleta de ganhos k <= kmax
, a parte
imaginária das raízes de D(s)+k*N(s)
é esboçada versus
a parte real.
Para obter o ganho num dado ponto do lugar geométrico, você pode
simplesmente executar a seguinte instrução:
k=-1/real(horner(h,[1,%i]*locate(1)))
e clicar no ponto
desejado sobre lugar geométrico das raízes. Se as coordenadas dos pontos
selecionados estão no vetor de reais 2 x 1 P=locate(1)
este k
resolve a equação k*N(w) + D(w)
=0
com w=P(1)+%i*P(2)=[1,%i]*P
.
Exemplos
H = syslin('c',352*poly(-5,'s')/poly([0,0,2000,200,25,1],'s','c')); clf(); evans(H,100); sgrid() // Post-tuning graphical elements ch = gca().children; curves = ch(2).children; curves.thickness = 2; asymptotes = ch(ch.type=="Segs"); asymptotes.segs_color = color("grey70"); // Completing Evans plot: [Ki,s] = kpure(H) // Gains that give pure imaginary closed loop poles plot([real(s) real(s)],[imag(s) -imag(s)],'*r') [Kr,s] = krac2(H)
Report an issue | ||
<< dtsi | Stability | g_margin >> |