Exemplos
H=syslin('c',352*poly(-5,'s')/poly([0,0,2000,200,25,1],'s','c')); clf();evans(H,100);sgrid() [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) |  |  |
lugar geométrico das raízes Evans
evans(H [,kmax])
lista (sistema linear syslin)
real (ganho máximo desejado para o esboç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.
H=syslin('c',352*poly(-5,'s')/poly([0,0,2000,200,25,1],'s','c')); clf();evans(H,100);sgrid() [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) | | |