Um sistema LIT de tempo discreto é descrito por
- Obtenha a função de transferência do sistema, tendo como saída e como entrada;
- Obtenha a resposta do sistema a um degrau unitário.
Passo 1
Uma função de transferência é a razão entre a saída e a entrada de um sistema. No tempo discreto, ela é obtida aplicando a transformada Z em uma equação das diferenças, considerando as condições iniciais nulas, e dividindo a saída pela entrada. Logo, a função de transferência será:
Onde é a tranformada Z de e , de .
Para resolver a letra a o exercício, primeiramente aplicaremos a função do MATLAB na equação de diferenças para encontrarmos a transformada de e , que indicaremos por e , respectivamente. Em seguida, manipularemos a equação para obtermos , que será a função de transferência .
Para resolver a letra b do exercício, calcula-se produto entre e a transformada Z do degrau do unitário, que é criado usando a função . Em seguida, usamos a função para calcular sua transformada inversa e encontrar a resposta ao degrau.
Passo 2
Para resolver a letra a, vamos rearranjar a equação, jogando todos os termos para a esquerda:
No MATLAB, vamos aplicar a função em cada termo da equação individualmente, pois dessa forma conseguiremos fazer a manipulação algébrica que nos dará a função de transferência. Lembrando que tem a sintaxe:
Ela retorna a transformada Z da função , armazenando-a na variável .
A transformada Z dos termos da equação das diferenças será identificada da seguinte maneira:
Onde indica a transformada Z. Logo:
Considerando tudo isso, em um novo script, digite as linhas de código:
Na primeira linha, definimos as variáveis simbólicas , e (indicamos que e são em função de ).
Na segunda linha, assumimos que , para que o degrau unitário , que será usado para resolver a letra b, seja igual a .
Na terceira, quarta, quinta, sexta e sétima linha, a transformada Z de cada um dos termos da equação das diferenças rearranjada é calculada, sendo armazenada em , , , e .
Após definir os termos da equação, construa a transformada Z dela a partir de , , , e , com o comando:
Desta forma, a transformada da equação é armazenada na variável .
Vamos ver como ficou a equação após a aplicação da transformada Z, executando o programa e imprimindo no Command Window:
Passo 3
Nossa, tá enorme, nem dá pra ler. Mas podemos deixar mais enxuto: ele tem termos que têm exatamente a mesma estrutura da função e que representam transformadas que podem ser substituídas por uma única letra, indicando que é a transformada de uma função de . representa a transformada Z de , podendo ser substituído por . No caso, chamaremos só de .
Já representa a transformada Z de , podendo ser substituído por . Como fizemos para , chamaremos só de .
Adicione os comandos abaixo em seu programa para fazer a substituição:
Primeiramente, criamos mais três variáveis simbolicas: , e . Em seguida, usamos a função para pedir ao MATLAB que substitua em a expressão por . Logo abaixo, fazemos a mesma coisa para . O resultado da substituição é armazenado em , ficando no lugar da expressão antiga:
Nossa, bem melhor. E tem mais coisa que podemos substituir em . Precisamos encontrar a função de transferência do sistema, que é determinada considerando as condições iniciais nulas. Logo, . Para fazer a substiuição, adicione ao seu programa a linha de código:
Com este comando, substituímos e por em e armazenamos a nova expressão em , atualizando a equação. Executando o programa, podemos ver como está agora:
Bem mais agradável de se olhar!
Passo 4
Para calcular , precisamos dividir por . Em , a variável representa e podemos usar a função do MATLAB para isolá-la na equação. Adicione em seu código:
Com este comando, a função resolve a equação para a variável e armazena a solução em , atualizando seu valor. Na prática, isolamos na equação com este comando e passa a representar a expressão algébrica de . Rodando o programa, temos que é igual a:
Para definir , divida por , já aplicando a função para eliminar o do numerador de :
E temos que a função de transferência é igual a:
Está aí a resposta da letra a ;)
Passo 5
Com , podemos resolver a letra b do exercício e calcular a resposta ao degrau unitário do sistema, que é dada pelo produto entre e a transformada Z do degrau do unitário. Para criar o degrau unitário, usaremos a função , que tem a estrutura:
Ela cria um degrau unitário e o armazena em : quando , . Se , . Finalmente, para , . Como dissemos que , temos .
Após criar o degrau com , aplicamos nele para obter sua transformada Z. Por isso, digite em seu programa os comandos:
Na primeira linha, criamos o degrau unitário e o armazenamos em . Em seguida, calculamos sua transformada Z e a registramos em .
Para encontrar a resposta ao degrau, inclua no seu programa o comando:
é o produto entre e e representa a transformada Z da resposta ao degrau unitário:
Passo 6
Para finalizar, vamos calcular a transformada inversa de para encontrar , a resposta do sistema ao degrau unitário no domínio de tempo discreto.
Para encontrarmos , vamos aplicar a função do MATLAB em . A função tem a estrutura:
Ela retorna a inversa da transformada Z de , armazenando-a na variável .
Para encontrar , digite os comandos:
Com o primeiro comando, o MATLAB calcula a inversa da transformada Z de e armazena o resultado em . O segundo comndo simplifica , para deixá-lo mais enxuto.
Rodando o programa, é igual a:
Resposta
- A função de transferência é:
- A resposta ao degrau unitário é: