Conversão de componentes do controlador PID com realimentação de estado em função de transferência única e forma discreta de espaço de estado

9

Estou lutando contra esse problema há cerca de uma semana, como parte de um projeto de um ano. Estamos projetando um controlador para um reator específico com base em um modelo. Depois de analisar isso por um tempo, ainda não consigo fazer isso funcionar - então eu realmente aprecio isso se conseguir alguma ajuda.

Uma das revisões de literatura publicadas em que nos baseamos fortemente lista um controlador PID em cada componente separado em vez de uma equação combinada, como segue:

{P(n)=Kp[G(n)-tumarget]Eu(n)=Eu(n-1 1)+KpTEu[G(n)-tumarget]D(n)=KpTDdGdt(n)

Simplesmente combinando os três componentes na saída do controlador PID:

PEuD(n)=P(n)+Eu(n)+D(n)

E a partir daqui, o autor adiciona uma camada adicional de feedback de estado sobre o sinal PID para obter a saída final do controlador aplicada no sistema.

{Q(n)=K0R(n-1 1)+K1 1Q(n-1 1)-K2Q(n-2)R(n)=(1 1+γ)PEuD(n)-γQ(n-1 1)

E R é a "saída do controlador" final. Aqui, é o ganho do processo, e são os ganhos integrais e derivativos, , e são "ganhos" ajustados para o feedback de estado (imutável) e é uma constante 0,5. é o estado do sistema, é um estado estimado que afeta a dinâmica do modelo e é a saída final real que é enviada para a planta.T I T D K 0 , K 1 K 2 γ G ( n ) Q ( n ) R ( n )KpTEuTDK0 0,K1 1K2γG(n)Q(n)R(n)

Eu estava tentando primeiro converter tudo em uma única função de transferência de controlador, mas me disseram que simplesmente adicioná-los não funcionaria.

Também fui encarregado de encontrar uma representação discreta do espaço de estado desse controlador. Para isso, tentei alterar o para para me livrar desse problema.G(n)-G(n-1)dGdt(n)G(n)-G(n-1 1)

Em seguida, tentei definir uma nova variável de estado para para que e pudessem ser convertidos em primeira ordem.Q ( n - 1 ) Q ( n - 2 )Q(n)Q(n-1 1)Q(n-2)

Tentei então substituir os valores no controlador PID, para obter o como a variável de estado. Esses esforços foram todos baseados em recomendações do meu professor.G(n)

No entanto, ainda estou muito preso, pois tenho seguido cegamente as instruções dele, sem uma visão geral para trabalhar nisso. Eu pensei que seria uma questão simples de transformação de Tustin - oh, como eu estava muito errado ...

Estou bastante frustrado, pois após um esforço de uma semana, ainda estou perplexo com o que parece ser um problema simples.

Se possível, posso humildemente pedir sua ajuda nessas duas perguntas específicas?

  1. Converta esse controlador em uma única função de transferência de controlador (como normalmente visto em qualquer representação de função de transferência, ou seja, )G(s)=1 1s+1 1
  2. Converter este controlador em uma representação discreta do espaço de estado, deixando a taxa de amostragem como uma variável?
John Galt
fonte
MATLAB e Maple podem resolver esses problemas. Eu tenho os dois programas. Imprimi sua postagem e tentarei trabalhar com eles. Eu fiz um pouco disso na faculdade.
Wesley Wortman
Você pode dar o título da publicação?
Hazem em Hazem

Respostas:

1

Não é uma resposta completa, mas espero que possa ser de alguma ajuda.

Você pode reescrever o primeiro sistema como

{P(n)=KPE(n)Eu(n)=Eu(n-1 1)+KPTEuE(n)ΔtD(n)=KPTDE(n)-E(n-1 1)Δt

Onde e é o seu intervalo de amostragem. Observe que e não são definidos como ganhos. e são respectivamente o ganho integral e o ganho derivado.E(n)=G(n)-tumarget(n)ΔtTDTEuKEu=KPTEuKD=KPTEu

Agora você pode reescrever o sistema como uma única função do erro.

PEuD(n)=P(n)+Eu(n)+D(n)

Eu(n-1 1)=PEuD(n-1 1)-P(n-1 1)-D(n-1 1)=PEuD(n-1 1)-KPE(n-1 1)-KPTDE(n-1 1)-E(n-2)Δt

PEuD(n)=KPE(n)+PEuD(n-1 1)-KPE(n-1 1)-KPTDE(n-1 1)-E(n-2)Δt+KPTEuE(n)Δt+KPTDE(n)-E(n-1 1)Δt=PEuD(n-1 1)+KP((1 1+ΔtTEu+TDΔt)E(n)-(1 1+2TDΔt)E(n-1 1)+TDΔtE(n-2))

O segundo é um pouco mais complexo para reescrever como uma única equação, mas você pode fazê-lo de maneira semelhante. O resultado deve ser

R(n)=K1 1R(n-1 1)-(γK0 0+K2)R(n-2)+(1 1+γ)(PEuD(n)-K1 1PEuD(n-1 1)+K2PEuD(n-2))

Agora você só precisa substituir a equação do PID para obter a equação do regulador em função do erro.

gvgramazio
fonte