Exercícios Resolvidos de Funções de Transferência

Ver Teoria

Enunciado

Plotar a resposta ao impulso de

G s = 3 e - 2 s 10 s + 1

Passo 1

O exercício pede para gerarmos o gráfico da resposta ao impulso de G ( s ) . Para isto, usaremos a função t f do MATLAB para criarmos a função de transferência G ( s ) e, em seguida, usaremos a função i m p u l s e para ver sua resposta ao impulso.

Contudo, você deve ter reparado que G ( s ) tem um formato um pouco diferente do que estamos acostumados: existe um exponencial no numerador. Um exponencial no numerador de uma função de transferência representa matematicamente um atraso entre o sinal de entrada e o de saída. É uma forma de modelar matematicamente atrasos em sistemas de controle reais. Por exemplo, em um reator químico ou em aeronaves, pode haver uma demora de alguns segundos do sinal de saída, devido a atuadores ou sensores.

Normalmente, funções de transferência que modelam atrasos no sinal tem o seguinte formato:

G s = H ( s )   e - τ s

Onde H ( s ) seria uma função de transferência no formato que estamos acostumados (polinômios no numerador e no denominador) e e - τ s modela um atraso de τ segundos entre o sinal de entrada e o sinal de saída.

No exercício, temos τ = 2 segundos e

H s = 3 10 s + 1

Veremos que, para incluirmos e - τ s em uma função de transferência usando o f t do MATLAB, basta adicionarmos mais uma entrada em f t indicando o tempo de atraso τ .

Passo 2

Para entendermos bem como funciona esse conceito de tempo de atraso entre os sinais de entrada e de saída, vamos primeiro trabalhar somente com H ( s ) e ver sua resposta ao impulso.

Usaremos a função f t para criarmos H ( s ) no MATLAB, que tem a sintaxe:

s i s t e m a = t f ( n u m e r a d o r , d e n o m i n a d o r )

A função t f tem como entrada dois vetores, o primeiro contendo os coeficientes do polinômio do numerador da função de transferência (representado por n u m e r a d o r ) e outro contendo os coeficientes do polinômio do denominador (representado por d e n o m i n a d o r ). t f retorna um objeto representando a função de transferência, que será armazenada em s i s t e m a .

O polinômio do numerador de H ( s ) é 3 . Já o polinômio do denominador é 10 s + 1 .

No nosso script, vamos criar duas variáveis representando os coeficientes dos polinômios. Como o polinômio do numerador tem somente um elemento correspondendo a 3 ∙ s 0 = 3 ∙ 1 = 3 , a variável representando seus coeficientes pode ser um número real. Chamaremos essa variável de n u m e ela é igual a 3 .

A outra variável é um vetor chamado d e n , contendo os coeficientes de 10 s + 1 , indo do elemento com maior expoente ( s ) ao de menor expoente ( 1 , pois ele multiplica s 0 = 1 ):

Tendo criado as entradas de t f , podemos incluir o comando abaixo para criar a função de transferência H ( s ) :

A função de transferência é armazenada na variável H , isto é, H é igual a H ( s ) . Podemos conferir isso rodando o programa e pedindo ao MATLAB para imprimir X no Command Window:

Passo 3

Após criarmos H ( s ) , vamos ver sua resposta ao impulso, sem o atraso de 2 segundos, isto é, desconsiderando o componente e - 2 s . A resposta ao impulso de uma função de transferência é calculada usando o comando i m p u l s e :

i m p u l s e ( s i s t e m a )

Onde s i s t e m a é a função de transferência criada usando o comando t f . Usando a sintaxe acima, plota-se automaticamente o gráfico da resposta ao impulso de s i s t e m a .

Inclua em seu script os comandos:

Na primeira linha, dizemos que o gráfico que será gerado por i m p u l s e ( H ) será a primeira figura. Incluimos esse comando para que os gráficos sem e com atraso sejam gerados em janelas diferentes e possamos compará-los.

Na segunda linha, usamos a função i m p u l s e para vermos a resposta ao impulso de H ( s ) .

Vamos rodar o programa para vermos a resposta ao impulso sem o atraso de 2 segundos:

Perceba que a curva azul já começa no instante 0 , que é o instante quando ocorre o impulso.

Passo 4

Agora vamos incluir no sistema o atraso de 2 s entre a entrada e a saída, isto é, vamos criar a função de transferência G ( s ) . Lembrando que um atraso entre o sinal de entrada e de saida pode ser modelado assim:

G s = H ( s )   e - τ s

Para modelar um atraso de tempo em uma função de transferência, devemos incluir mais uma entrada na função t f , que terá o seguinte formato:

s i s t e m a = t f ( n u m e r a d o r , d e n o m i n a d o r ,   ' I n p u t D e l a y ' , τ )

Os vetores n u m e r a d o r e d e n o m i n a d o r contêm os coeficientes dos polinômios do numerador e do denominador de H ( s ) .

Além dos vetores com os coeficientes dos polinômios do numerador e do denominador de H ( s ) , existe uma terceira entrada: ' I n p u t D e l a y ' , τ . Essa terceira entrada indica que há um atraso de τ segundos entre o sinal de entrada e o de saída.

Para o exercício, como já criamos n u m e d e n com os coeficientes de H ( s ) no passo 2, para criarmos G ( s ) , vamos só adicionar ao nosso script o comando:

Com este comando, criamos a função de transferência G ( s ) do enunciado e a armazenamos em G . Ao incluirmos ‘ I n p u t D e l a y ’ , 2 como entrada, indicamos que τ = 2 em G ( s ) .

Se você rodar o programa que tem até agora e pedir para ver G , você verá que G é exatamente a função de transferência do enunciado:

Passo 5

Após criar G ( s ) , vamos ver sua resposta ao impulso. Para isto, adicione em seu script os comandos:

Na primeira linha, dizemos que o gráfico que será gerado por i m p u l s e ( G ) será a segunda figura. Na segunda linha, usamos a função i m p u l s e para vermos a resposta ao impulso de G ( s ) .

Ao executarmos nosso programa, teremos dois gráficos gerados. Você pode colocá-los lado a lado para compará-los:

À esquerda, temos a Figure 1, que tem a resposta ao impulso de H ( s ) , isto é, sem o atraso de 2 s.

À direita, temos a Figure 2, que tem a resposta ao impulso de G ( s ) , isto é, com o atraso de 2 s. Perceba que há uma diferença de em torno de 2 segundos entre as curvas azuis da esquerda e da direita. Isso ocorre devido ao exponencial de G ( s ) , que representa esse tempo de atraso.

Resposta