Exercícios Resolvidos de Transformada de Laplace

Ver Teoria

Enunciado

O circuito RLC abaixo pode ser usado como um filtro de banda estreita. Se a tensão de entrada v i ( t ) consiste na soma de tensões de variação sinusoidal com diferentes frequências, o filtro de banda estreita permitirá passar somente tensões cujas frequências estão dentro de uma faixa estreita. Esses filtros são usados em circuitos de sintonia, como os usados em rádios AM, para permitir a recepção apenas do sinal portador da estação de rádio desejada.

O circuito RLC tem a seguinte equação diferencial:

L C v ¨ 0 + R C v ˙ 0 + v 0 = R C v ˙ i ( t )

Onde v 0 ( t ) é a tensão de saída.

Use o método da transformada de Laplace para encontrar a resposta ao degrau unitário de v 0 ( t ) para condições iniciais nulas, onde C = 10 ∙ 10 - 5   F e L = 5 ∙ 10 - 3   H . Para modelar um filtro de banda larga, use R = 1000   Ω . Já para modelar um filtro de banda estreita, use R = 10   Ω . Compare a resposta ao degrau dos dois casos.

Passo 1

O exercício pede para encontrarmos a resposta de v 0 ( t ) a v i ( t ) igual a um degrau unitário, usando a transformada de Laplace. Em outras palavras, vamos encontrar a solução v 0 da equação diferencial de segunda ordem dada pelo enunciado. Primeiramente, aplicaremos a função l a p l a c e do MATLAB na equação diferencial para encontrarmos a transformada de v 0 . A transformada de Laplace de v 0 ( t ) será indicada por V 0 ( s ) .

Após encontrarmos V 0 ( s ) , vamos aplicar i l a p l a c e para encontrar a inversa da transformada, isto é, v 0 ( t ) , que é a solução da equação diferencial.

Vamos fazer esse procedimento considerando C , L e R simbólicos, substituindo-os pelos valores dados pelo enunciado somente no final. Para compararmos a resposta ao degrau, substituiremos em v 0 ( t ) os dois valores de R e veremos como ficará a expressão simbólica para ambos os casos.

Passo 2

Comece entrando com a equação diferencial no MATLAB, termo a termo. Em um script novo, digite as linhas de código:

Na primeira linha, criamos as variáveis simbólicas L , C , R e v _ 0 ( t ) , indicando que v _ 0 é em função de t .

Na segunda linha, calculamos a derivada v ˙ 0 usando a função d i f f ( v _ 0 , t ) do MATLAB. Como estamos trabalhando com variáveis simbólicas, d i f f ( v _ 0 , t ) apenas representa a derivada de v _ 0 em relação a t e é armazenado em D v _ 0 .

Na terceira linha, calculamos a derivada v ¨ 0 , pedindo ao MATLAB para calcular a derivada de D v _ 0 em relação a t ( d i f f ( D v _ 0 , t ) ).

Na quarta linha, criamos um degrau unitário v _ i em função de t usando a função h e a v i s i d e ( t ) , representando a tensão de entrada do circuito.

Na quinta linha, definimos a primeira derivada do degrau unitário v ˙ i usando d i f f ( v _ i , t ) . A derivada é armazenada em D v _ i .

Tendo definido os termos da equação, vamos criá-la e armazená-la na variável e q :

Repare que joguei todos os termos da equação da para a esquerda, para fazer com que ela seja igual a 0 .

Passo 3

Com a equação diferencial no MATLAB, podemos usar a função l a p l a c e para calcular a transformada de Laplace dela. Lembrando que l a p l a c e tem como estrutura:

F = l a p l a c e ( f )

Ela retorna a transformada de Laplace da função f , armazenando-a na variável F .

Acrescente em seu programa o comando abaixo para calcular a transformada de Laplace de e q :

O resultado é armazenado em E . Vamos ver como ele ficou? Rode o programa que fizemos até agora e imprima E no Command Window:

Ficou uma expressão bem comprida, tive que deixar a letra minúscula para caber aqui. Mas podemos dar uma enxugada nessa expressão: note que aparecem uns termos l a p l a c e ( v _ 0 ( t ) , t , s ) em E , e eles têm exatamente a mesma estrutura da função l a p l a c e . Essa expressão representa a transformada de Laplace de v 0 ( t ) e pode ser substituída por V 0 ( s ) . No caso, chamaremos V 0 ( s ) só de V 0 .

Digite os comandos abaixo para fazer a substituição:

Primeiramente, devemos criar mais duas variáveis simbólicas: s e V 0 . Em seguida, usamos a função s u b s para pedir ao MATLAB que substitua em E a expressão l a p l a c e ( v _ 0 ( t ) , t , s ) por V 0 . O resultado da substituição é armazenado em E , atualizando a expressão anterior. Após a susbtituição, E fica assim:

Já deu uma boa enxuagada. Agora vamos substituir as condições iniciais em E . O enunciado disse que elas são nulas, logo, v 0 0 = v ˙ 0 0 = 0 . Na expressão acima, é fácil de identificar quem é v 0 ( 0 ) . Já v ˙ 0 0 é representado por D ( v _ 0 ) ( 0 ) , um pouco mais sutil.

Para fazer as substituições, adicione ao seu programa a linha de código:

Com este comando, substituímos em E v _ 0 ( 0 ) e D ( v _ 0 ) ( 0 ) por 0 e armazenamos a nova expressão em E , no lugar da antiga. Note que usei aspas simples em ' D ( v _ 0 ) ( 0 ) ' , para indicar que os parênteses pertencem a uma expressão só. Se você tentar rodar o programa sem elas, o MATLAB acusará erro de sintaxe.

Agora E está assim:

Passo 4

Agora podemos encontrar v 0 ( t ) . Para isto, precisamos encontrar V 0 ( s ) para daí aplicar nele i l a p l a c e , para termos a inversa da tranformada de Laplace. Em E , a variável V 0 representa V 0 ( s ) e podemos usar a função s o l v e do MATLAB para encontrá-la. Adicione em seu código:

Com este comando, a função s o l v e resolve a equação E para a variável V 0 e armazena a solução em V 0 , atualizando seu valor. Na prática, isolamos V 0 na equação e V 0 passa a representar a expressão algébrica de V 0 ( s ) . Rodando o programa, temos que V 0 é igual a:

Para encontrarmos v 0 ( t ) , devemos aplicar a inversa da transformada de Laplace em V 0 ( s ) , isto é, usar a função i l a p l a c e do MATLAB em V 0 . A função i l a p l a c e tem a estrutura:

f = i l a p l a c e ( F )

Ela retorna a inversa da transformada de Laplace de F , armazenando-a na variável f .

Logo, em seu código, entre o comando:

Com este comando, o MATLAB calcula a inversa da transformada de Laplace de V 0 e armazena o resultado em v _ 0 , atualizando seu valor. v _ 0 é a solução da equação diferencial, representadno a resposta do sistema ao degrau unitário. Rodando o programa, temos que v _ 0 é igual a:

A expressão acima ficou grandinha, mas ela vai encolher quando substituirmos R , L e C pelos valores fornecidos no enunciado.

Passo 5

Vamos comparar as respostas ao degrau de um filtro de banda larga, com R = 1000   Ω , e de um filtro de banda estreita, com R = 10   Ω . O primeiro passo é substituir em v _ 0 C e L por C = 10 ∙ 10 - 5   F e L = 5 ∙ 10 - 3   H , usando o comando s u b s :

Com este comando, v _ 0 foi atualizado com valores numéricos de C e L .

Em seguida, vamos substituir R por R = 1000   Ω , para ver a resposta ao degrau do filtro de banda larga:

Armazenamos em v _ 0 _ b l a expressão de v _ 0 com R = 1000 . Vamos ver como ela ficou?

Podemos transformar as frações e as exponenciações em números decimais, usando o comando v p a do MATLAB. Este comando avalia cada elemento de uma expressão simbólica e identifica as frações ou as exponenciações com números reais, e as transforma em números “quebrados”. Por isso, acrescente no seu código o comando:

Este comando aplica a função v p a em v _ 0 _ b l e armazena o resultado em v _ 0 _ b l , substituindo a expressão anterior. Ele vai nos dar os coeficientes de v _ 0 _ b l com 7 algarismos significativos. Ao rodar programa, podemos ver v _ 0 _ b l com seus coeficientes transformados em números decimais:

Repare que v _ 0 _ b l tem números complexos, identificados por multiplicarem i .

Vamos agora ver como é v 0 ( t ) em um filtro de banda estreita, isto é, vamos substituir R por R = 10   Ω :

Armazenamos em v _ 0 _ b e a expressão de v _ 0 com R = 10 . Rode o programa e imprima-a no Command Window para ver v _ 0 _ b e :

Bem mais simpática que v _ 0 _ b l , não?

Resposta

A resposta ao degrau do filtro de banda larga é:

v 0 B L ( t ) = - 2 i ∙ e - 10 5 t sin ⁡ ( 99990 i ∙ t )

A resposta ao degrau do filtro de banda estreita é:

v 0 B E ( t ) = 2 e - 10 3 t sin ⁡ ( 1000 t )