evans

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'));
evans(H,100)
P=3.0548543 - 8.8491842*%i;    //P=ponto selecionado
k=-1/real(horner(H,P));
Ns=H('num');Ds=H('den');
roots(Ds+k*Ns)     //contém P como raiz particular
// outro
clf();s=poly(0,'s');n=1+s;
d=real(poly([-1 -2 -%i %i],'s'));
evans(n,d,100);
//
clf();n=real(poly([0.1-%i 0.1+%i,-10],'s'));
evans(n,d,80);

Ver Também

• kpure
• krac2
• locate