Use o critério de Nyquist para determinar a estabilidade de um sistema com realimentação unitária e função de transferência mostrada abaixo:

A função serve para criar funções de transferência e tem como sintaxe:

Ela tem como entrada dois vetores, o primeiro contendo os coeficientes do polinômio do numerador da função de transferência (representado por ) e outro contendo os coeficientes do polinômio do denominador (representado por ). retorna um objeto representando a função de transferência, que será armazenada em .

Para , o polinômio do numerador é . Já o polinômio do denominador é .

Em um novo script, vamos as entradas de . Como o polinômio do numerador tem somente um elemento correspondendo a , a variável contendo seus coeficientes pode ser um número real. Chamaremos essa variável de e ela é igual a .

A outra variável é um vetor chamado , contendo os coeficientes de , indo do elemento com maior expoente ao de menor expoente:

Logo abaixo das linhas onde e foram criados, já digitamos o comando para criar , que será armazenado em .

Antes de continuar, verifique pra ver se não digitou nada de errado. Para isto, rode o programa que você tem até agora e imprima no Command Window:

Estando tudo OK, podemos plotar o diagrama de Nyquist.

Use a função para gerar o diagrama de Nyquist, que tem como sintaxe:

Ela tem como entrada uma função de transferência, criada usando o comando do MATLAB (representada por acima), e ela gera um gráfico com o diagrama de Nyquist para .

Adicione em seu código o comando:

Execute o programa para ver o diagrama de Nyquist de :

Meio esquisito, né? Talvez você se assuste com esse formato diferente do que vimos nos exercícios anteriores, mas você vai ver no próximo passo que a análise do diagrama pelo critério de Nyquist não tem nada de diferente do que vimos até aqui.

Vamos utilizar o diagrama de Nyquist plotado para avaliar a estabilidade do sistema por meio do critério de Nyquist, que é dado pela equação:

Onde:

• é o número de polos positivos da função de transferência em malha fechada;
• é o número de voltas que o diagrama de Nyquist dá em torno de . é positivo se a volta for no sentido horário. Se for no sentido anti-horário, é negativo.
• é o número de polos positivos da função de transferência em malha aberta.

Pelo critério de Nyquist, para que um sistema seja estável, deve ser igual a .

O é obtido calculando as raízes do denominador da função de transferência em malha aberta, no caso, :

O denominador de é , que tem como raízes e . Logo, não tem polos positivos e .

é obtido pelo diagrama de Nyquist. Vamos rever o diagrama de Nyquist gerado no passo anterior:

Note a cruzinha vermelha na lateral esquerda do gráfico. Essa cruzinha marca o ponto de coordenadadas . Como a cruzinha vermelha está fora do diagrama de Nyquist (representado pela curva azul), o diagrama não dá volta alguma em torno dela, por isso, .

Tendo e , calculamos . Temos que . Como é igual a , o sistema é estável.