Inicialmente, deve-se carregar os dados fornecidos pelo professor, para isso, basta
salvar o arquivo [Link] na pasta onde será implementado o projeto e construir o código
dentro da mesma pasta. Dessa forma, ao utilizar a função load() do Matlab, basta indicar o no
do arquivo sem precisar especificar o caminho.
Pode-se ver que o arquivo [Link] e o arquivo TunelVento.m desenvolvido estão na
mesma pasta chamada TunelVentoMatlab. Ao executar o código serão carregadas as variáveis
t referente ao tempo, u referente a entrada e y referente a saída.
Com os dados em memória, pode-se abrir o App – System Identification.
Ao abrir, inicialmente devemos importar os dados, em import data – time domain data.
Indicar o nome das variáveis de entrada e de saída, nome dos dados e indicar que o tempo
iniciará em 0 com passo de 0.1 e clicar em Import.
Os dados serão apresentados em Data Views. Para visualizar basta clicar em time plot.
Com os dados, agora é hora de estimar a função de transferência, para isso, clica-se
em Estimate – Transfer Functions Models
Indica-se a quantidade de polos e zeros, com a quantidade de polos maior que a de
zeros e estimate.
Para visualizar, clica-se na modelo criado e em Model Output. Serão indicados o sinal
real, o criado pelo modelo e a aproximação do modelo.
Criou-se modelos aumentando a quantidade de polos e zeros em 1 até que a
aproximação seja maior que 95%. Dessa forma, foram criados o modelo TV_FT1 com 0 zeros
e 1 polo, TV_FT2 com 1 zero e 2 polos, TV_FT3 com 2 zeros e 3 polos, TV_FT4 com 3 zeros
e 4 polos, TV_FT5 com 4 zeros e 5 polos, TV_FT6 com 5 zeros e 6 polos, TV_FT7 com 6 zeros
e 7 polos.
O modelo TV_FT7 possui o seguinte formato:
Para levar o modelo até as variáveis do MATLAB arrasta-se o modelo até “To Workspace”
Com o modelo na memória, pode-se plotar o lugar das raízes do modelo, para isso,
utiliza-se da função rlocus() do MATLAB:
Pelo lugar das raízes pode-se notar o lugar dos polos, dos zeros e um dado importante
é a passagem dos polos para o semi-plano direito indicando instabilidade. Para observar os
valores de ganho que levam a instabilidade, basta verificar a fronteira:
O sistema passa a ter polos com parte real positiva quando o ganho se aproxima de
0,995, dessa forma, esse sistema é estável para 0 < 𝐾 < 0,995
Outra informação importante é indicar quais são os polos, zeros e o ganho da função
de transferência do modelo. Para isso encontra-se o numerador, o denominador e a partir deles
encontra-se o zpk do sistema.
É possível notar que todos os polos possuem parte real negativa, indicando que o
sistema em malha aberta é estável.