Sobre projeto de observadores desacoplados de perturbação para sistemas descritores

IDENTIFICAÇÃO DE SISTEMAS

Sobre projeto de observadores desacoplados de perturbação para sistemas descritores

Vilemar G. da Silva^I; Eugênio B. Castelan^II; Sophie Tarbouriech^III; Germain Garcia^III

^IDE.EE-CCET-UFMA Campus universitário do Bacanga-65085-580, São Luís, MA, Brasil [e-mail]

^IIDAS-CTC-UFSC Departamento de Automação e Sistemas - 88040-900, Florianópolis, SC, Brasil [e-mail]

^IIILAAS-CNRS 7, Av. du Colonel Roche -31077, Toulouse, France [e-mail], [e-mail]

RESUMO

Condições necessárias e suficientes para a existência de observadores de ordem reduzida desacoplados de perturbação para sistemas descritores lineares contínuos são apresentadas. Essas condições são baseadas na solução de uma equação de Sylvester generalizada, sob uma restrição que descreve o desacoplamento de perturbação do estado estimado e numa condição de posto que garante a existência de matrizes para a reconstrução de todas as variáveis de estado, também com desacoplamento de perturbação. Mostra-se que as condições para a existência do observador desejado estão associadas, fundamentalmente, a certas propriedades estruturais do sistema descritor. A partir destes resultados, propõe-se um algoritmo de projeto baseado em técnica de posicionamento de autoestrutura. Os resultados apresentados destinam-se, fundamentalmente, à aplicação em síntese de leis de controle do tipo realimentação de estados estimados. Um exemplo numérico ilustrativo do procedimento de projeto é apresentado.

Palavras-chave: Sistemas descritores, Observadores, Ordem mínima, Equação de Sylvester, Posicionamento de Autoestrutura.

ABSTRACT

Necessary and sufficient conditions for the existence and design of disturbance decoupled reduced-order observers for continuous linear descriptor systems are presented. The observer existence conditions are based on the solution of a generalized Sylvester's equation, under a constraint that describes the disturbance decoupling of the estimated state, and on a rank condition that guarantees the existence of matrices for the reconstruction of all state variables, also with disturbance decoupling. Basically, the design conditions of the desired observer are associated to some structural properties of the descriptor system. From these results, we propose an observer design algorithm that is based on eigenstructure assignment technique. The presented results are useful when an estimated state feedback control law is considered. An illustrative numerical example is reported.

Keywords: Descriptor Systems, Observers, Minimal order, Sylvester's equation, Eigenstructure assignment.

1 INTRODUÇÃO

A teoria de estimação de estados exerce um papel importante em diferentes campos técnicos e científicos tais como em controle de sistemas dinâmicos (Chen, 1999) (Dai, 1989) (Liu and Patton, 1998), em detecção/diagnóstico de falhas (Liu and Patton, 1998) (Bolivar, 2001) e em filtragem (Nagpal and Khargonekar, 1991; Yaesh and Shaked, 1992). Nas aplicações em controle de sistemas dinâmicos, os observadores (ou estimadores) de estado viabilizam a implementação de leis de controle de realimentação de estados, parcial ou completamente estimados, preservando, em particular, a estabilidade do sistema em malha fechada para ambos os casos. Nesse campo de aplica cão, um observador de estados constitui-se em um subsistema dinâmico, excitado pelas entradas e saídas do sistema a controlar. Os estados do observador devem permitir recuperar, assintoticamente, as trajetórias dos estados do sistema ou de uma combinação deles.

Considerando que o sistema a controlar está sujeito a perturbações externas, é desejável que o comportamento dinâmico do observador esteja, parcialmente ou completamente, desacoplado das perturbações, ou satisfa ca algum critério de desempenho face à influência dessas perturbações no erro de observação do sistema. Portanto, um observador desacoplado de perturbações externas consiste de um sistema observador que rejeita completamente alguma perturbação considerada do erro de observação.

Embora o desacoplamento (ou rejeição completa) de perturbação possa parecer um requisito muito forte, observadores desacoplados de perturbação podem ser projetados se o sistema considerado verifica algumas condições estruturais, que estão geralmente relacionadas a alguma propriedade de detectabilidade e ao número de saídas medidas comparado ao número de entradas de perturbação. Além disso, algumas propriedades inerentes a um projeto por realimentação de estados, relativas à relação perturbação-saída, continuam válidas em malha fechada quando um observador desacoplado de perturbação é utilizado para implementar a lei de controle.

O desacoplamento de perturbações em projeto de observadores de estado foi tratado na literatura científica, principalmente para o caso de sistemas lineares normais (descritos por equa cões diferenciais ordinárias de primeira ordem), como por exemplo, em (Johnson, 1975) (Bhattacharyya, 1978) (Hou and Muller, 1992a) (Syrmos, 1993) (Liu and Patton, 1998). O desacoplamento de perturbações em observadores para o caso de sistemas descritores (descritos por equações diferenciais de primeira ordem e equações algébricas) foi tratado, por exemplo, em (Yang and Tan, 1989) (Paraskevopoulos et al., 1992) (Hou and Muller, 1992b) (Darouach et al., 1996) (Chu and Mehrmann, 1999a).

O objetivo deste trabalho é abordar o problema de projeto de observadores desacoplados de perturbações externas para sistemas descritores. Relativamente a trabalhos similares, a originalidade da técnica de projeto proposta reside: i) na resolução explícita de uma equação de Sylvester generalizada sujeita a duas condições complementares, e ii) na utilização de sinais auxiliares, definidos em função das entradas de controle e das saídas medidas, respectivamente. Na definição desses sinais auxiliares, são empregadas matrizes que permitem anular a ação da perturbação sobre as variáveis algébricas do sistema descritor e sobre as variáveis de saída. Já a possibilidade de solução da equação de Sylvester conjuntamente com as condições complementares, é interpretada em termos de propriedades estruturais que devem ser verificadas pelo sistema descritor para que se possa projetar o observador desacoplado de perturbação (ODP).

Este trabalho está organizado da maneira seguinte. Na se cão é feita uma apresenta cão do problema, onde são mostradas as equações do sistema descritor, do observador e dos sinais auxiliares a serem utilizados. Na seção são apresentadas condições necessárias e suficientes para a existência do observador desejado, e uma breve análise considerando a inclusão de uma dada lei de controle de realimentação de estados estimados. Na seção : i) são apresentadas e interpretadas as condições estruturais do sistema descritor que permitem o projeto do observador, ii) relacionam-se os quatro casos possíveis relativos à ação das perturba cões sobre as variáveis do sistema, indicando-se a estrutura/ordem mínima do observador mais adequada para cada caso. Na seção propõe-se um algoritmo de cálculo com base nos resultados anteriores e mostram-se os resultados numéricos e de simulação, para um exemplo de sistema descritor. Em seguida apresentam-se as conclusões.

2 APRESENTAÇÃO DO PROBLEMA

Seja um sistema descritor linear contínuo sujeito a perturbações externas descrito por:

onde x Î Âⁿ, y Î Â^p, u Î Â^m e w Î Â^r são, respectivamente, os vetores de estados, de saídas medidas, de entradas de controle e de entradas de perturbação. A matriz E Î Â^n×n é tal que posto(E) = q < n, que caracteriza um modelo de sistema descritor linear, A Î Â^n×n, B Î Â^n×m, com posto(B) = m, C Î Â^p×n, com posto(C) = p, B_w Î Â^n×r, com posto(B_w) = r e D_w Î Â^p×r, com posto(D_w) = < p¹.

Seja L Î Â^{(n - q)×n} uma matriz de posto completo que verifica

Defina B_d Î Â^{(n - q)×r} e o escalar por

Para efeito de implementação de leis de controle do tipo realimentação de estados estimados u(t) = F(t) + v(t), onde v(t) representa uma nova entrada de controle ou de referência, considera-se um observador de estados representado sob a forma:

onde z Î Â^k e Î Âⁿ são os vetores de estados estimados e de estados observados, respectivamente.

Os vetores Î e Î representam sinais auxiliares gerados a partir das entradas de controle e das saídas medidas, como segue:

com Î tal que B_d = 0, posto() = n - q - , e

com Î tal que D_w = 0 e posto() = p-.

A figura 1 ilustra o sistema em malha fechada, contendo o sistema original, o observador e um controlador por realimentação de estados estimados. Observe que a parte dinâmica do observador, descrita por sua equação de estados (5), é similar à dos observadores de ordem reduzida para sistemas normais (ver por exemplo (Chen, 1999)). Como pode ser visto nas duas seções seguintes, a definição do sinal permite levar em conta a parte algébrica das variáveis de estado do sistema descritor, no projeto do observador.

O problema tratado na seqüência destina-se ao projeto de observadores (5)-(6) para sistemas descritores (1)-(2), sem influência das perturbações w(t). Deseja-se verificar sob quais condições é possível determinar as matrizes H, T, Z, S, e tais que o erro entre x(t) e (t) é nulo em regime estacionário, para qualquer w(t) Î Â^r. Um sistema observador de estados sob essas condições será dito um Observador Desacoplado de Perturbações (ODP).

A ordem do observador (5)-(6) é determinada não somente em função da ordem dinâmica do sistema descritor, q, e do número de saídas, p, mas também de como a perturbação atua sobre a parte algébrica das variáveis de estado (Dai, 1989) e sobre as variáveis de saída, através das matrizes e , respectivamente. Em particular, se LB_w = B_d = 0 e D_w = 0, tem-se = I_n-q e = I_p, e a ordem mínima é k = q - p (Castelan et al., 2004). Nas outras situações, a estrutura proposta também permite atingir a menor ordem possível.

3 OBSERVADORES DESACOPLADOS DE PERTURBAÇÕES

Nesta seção são apresentadas as condições para a existência de um ODP de ordem reduzida k > q - p, essencialmente descritas pela possibilidade de solução de uma equação de Sylvester conjuntamente com uma condição para o desacoplamento das perturbações. Ao final, considera-se a uma lei de controle de realimentação de estados estimados, e faz-se uma breve análise para o sistema de malha fechada correspondente.

Proposição 3.1 Dado um sistema descritor (1)-(2), existe um observador de ordem reduzida k na forma (5)-(6) tal que, para qualquer perturbação w Î Â^r:

i) o erro de estimação e(t) = z(t) - TEx(t) verifica: e(t) = e^Hte(0),"t > 0ee(t) = 0, "z(0), Ex(0); e

ii) para matrizes S, e , adequadamente calculadas, o estado estimado Î Âⁿ verifica [x(t) - (t)] = 0;

se e somente se é possível resolver a equação de Sylvester generalizada

sob as condições adicionais:

onde: H Î Â^k×k, T Î Â^k×n e Z Î Â^k×p; e s(H) denota o espectro de H.

Demonstração:

Considere o erro de estimação

Sua derivada, em relação ao tempo, é dada por:

Fazendo substituições entre as equações (1), (2), (5) e (13) obtém-se, após breve simplificação:

Observa-se que a parte (i) da proposição é satisfeita se e somente se a dinâmica do erro de estima cão for descrita pelo sistema dinâmico assintoticamente estável:

Portanto, comparando-se (14) com (15), conclui-se que as condições (9) e (10) são necessárias e suficientes para o desacoplamento de perturbação do erro de estima cão relativo à parte (i) da proposição.

Para a demonstração da parte (ii), supor que (11) também é verificada e considerar (6) escrita na forma

Multiplicando, à esquerda, a equação (1) por

L e lembrando que, por definição, LE = 0, obtém-se

Também por definição,

LB_ww(t) = 0. Então, combinando (7) com (17) resulta em

De (2) e (8), lembrando que

D_w = 0, obtém-se

Conforme a demonstração da parte (i), e(t) = 0. Então, no limite quando t ® ¥, as equações (12), (18) e (19) podem ser escritas conjuntamente como

Tomando o limite de (16), com t ® ¥, e combinando o resultado com (20), obtém-se:

em que M^† = [ S ] Î é, em geral, uma pseudo-inversa à esquerda de M = [(TE)' (LA)' (C)']'.

Logo, independentemente do valor de

x(t), o erro de observação x(t) - (t) é nulo em regime permanente, se e somente se for possível anular o termo entre parênteses na relação anterior ou, equivalentemente, se for possível calcular as matrizes S, e de modo que

Isto pode ser realizado, se e somente se a condição (11) é verificada.

A resolução de (9)-(10) está associada, fundamentalmente, ao número de entradas de perturbação e a propriedades estruturais ligadas ao sistema (E,A,C,B_w,D_w), a serem discutidas posteriormente.

Da demonstração acima, é importante salientar o papel dos sinais auxiliares e , definidos a partir das matrizes B_d = LB_w e D_w, fundamentais para a obtenção do erro de estados estimados desacoplado da perturbação. De fato, se tivéssemos B_d ¹ 0 e D_w ¹ 0, então, ao invés de (20), teríamos

o que implicaria, em particular, [x(t) - (t)] ¹ 0 se w(t) ¹ 0.

Em função da definição dos sinais auxiliares e , descrevem-se na seção seguinte diferentes casos para os quais pode-se considerar o projeto de observadores desacoplados de perturbação sob a estrutura proposta, associando-se, em particular a definição da ordem do observador à condição (11). Note que, em função de (11), a ordem do observador deve verificar k > q - p + + (veja também (21) e (22)). Então, a ordem mínima possível para o observador ocorre no caso em que a matriz M = [(TE)' (LA)' (C)']', de posto n, é quadrada e portanto inversível, de modo que M^† = M^-1.

Para finalizar esta seção, considera-se a utiliza cão do observador desacoplado de perturbação para a implementação de uma lei de controle do tipo realimenta cão de estados estimados, dada por

em que F Î Â^m×n é uma matriz de ganho de realimentação e v(t) é uma nova entrada de controle ou de referência. Supõe-se que a matriz F é tal que o par (E, A + BF) é fortemente estável, ou seja: assintoticamente estável, regular e livre de modos impulsivos (Varga, 1995).

O sistema composto em malha-fechada pode ser representado (ver Apêndice) por:

O sistema (25)-(26) é também um sistema descritor, no qual, em particular, posto = q + k. A partir das estruturas das matrizes em (3)-(23), verifica-se que:

1) pela propriedade da separação, o sistema em malha fechada possui q + k pólos finitos assintoticamente estáveis, dados por s(E, A + BF) s(H) Ì ^-, e portanto é fortemente estável;

2) devido às características do observador desacoplado de perturbação, a ação da perturbação sobre o sistema em malha fechada ocorre diretamente através de B_w e D_w e, como no caso de uma realimentação de estados medidos, u(t) = Fx(t) + v(t), o comportamento perturba cão-saída, considerando-se x(0) tal que Ex(0) = 0, pode ser representado pela matriz de transferência:

3) em regime transitório, o comportamento dinâmico do sistema em malha fechada é influenciado pela dinâmica atribuída ao observador, através das matrizes H e S obtidas pelo projeto do observador, e pelas entradas w(t) e v(t).

4 RESULTADOS COMPLEMENTARES

A seguir apresentam-se as hipóteses que permitem o cálculo de um ODP. Justificativas para a utilização dessas hipóteses são apresentadas nas subseções seguintes.

Hipótese 4.1 O sistema (1)-(2):

i) possui menos entradas de perturbação que saídas medidas, ou seja:

ii) é fortemente zero-detectável, o que requer as duas condições seguintes:

iii) e, adicionalmente:

Em geral, para a construção de observadores para sistemas descritores, é necessário que o par (E,A) seja fortemente detectável, o que requer a satisfação da condição posto = n, " l Ï , l finito, em conjunto com (30), (veja, por exemplo, (Castelan and Silva, 2005)). Então, a condição (29), sob a condição estrutural (28), constitui-se numa adapta cão que permite resolver a equação de Sylvester (9) sob a condição de desacoplamento (10). Note que a condição (30) também é necessária para que (11) possa ser verificada via a síntese da matriz T do observador. Por último, a condição adicional (31) está associada, fundamentalmente, à determinação da ordem do observador.

4.1 Condições estruturais para solução

Técnica padrão de posicionamento de autoestrutura pode ser usada para resolver a equação de Sylvester (9) sob a condição de desacoplamento (10) (Liu and Patton, 1998). Por exemplo, se considerarmos s(H) = {p₁, p₂, ..., p_k} Ì ^-, com p_i+1 = se Imag(p_i) ¹ 0², as matrizes T e Z podem ser formadas a partir dos vetores t_i Î ⁿ e z_i Î ^p que verificam:

na qual t_i+1 = e z_i+1 = se mag(p_i) ¹ 0.

A existência de vetores não-nulos que permitem a escolha arbitrária dos autovalores p_i está relacionada às propriedades estruturais do sistema representado por (E,A,C,B_w,D_w) e que podem ser descritas a partir da matriz do sistema P(l), de dimensão (n + p) × (n + r), dada por:

Para todo valor arbitrário de p_iÎ , a existência de um vetor não-nulo que satisfaz (32) é garantida se o número de entradas de perturbações é menor que o número de saídas medidas, ou seja, se a condição (i) da Hipótese 4 é satisfeita. Observe, entretanto, que se para algum valor particular de l tivermos posto(P(l)) < n + r, este valor deverá fazer parte obrigatoriamente do espectro da matriz H. Quando r > p, a matriz P(l) perde posto somente para valores particulares de l e a resolução de (9)-(10) não é geralmente possível.

A condição (30) corresponde à detectabilidade dos pólos infinitos do sistema representado por (E,A,C) (Castelan and Silva, 2005). Além disso, utilizando-se a desigualdade de Sylvester (Chen (1999), pag.207) sobre o posto de produto de matrizes, pode-se verificar que a condição (30) é necessária para que a condição (11) possa ser satisfeita.

4.2 Ordem do Observador

Na classificação apresentada a seguir, consideram-se certas propriedades das matrizes B_d = LB_w e D_w e, a partir delas, identificam-se a ordem mínima do ODP em quatro casos possíveis.

Observe, inicialmente, que a ordem mínima do ODP é determinada pelo número mínimo de linhas linearmente independentes que deverá ter o produto TE, para que a condição de posto (11) seja satisfeita. Então, ela é determinada pela relação seguinte:

Por hipótese, a condição (31) é verificada. Então, utilizando-se a desigualdade de Sylvester (Chen (1999), pag.207) sobre o posto de produto de matrizes, tem-se que

Assim, de (34) e (35), resulta:

Levando em conta a estrutura das matrizes B_d = LB_w e D_w, quatro casos distintos podem ser identificados. Então, utilizando (36), a ordem mínima do ODP associado à cada caso é determinada como segue:

Caso 1 Se B_d ¹ 0 e D_w ¹ 0 utiliza-se um observador de ordem mínima k = q - p + + .

Caso 2 Se B_d = 0 e D_w = 0 utiliza-se um observador de ordem mínima k = q - p, no qual Î Â^{(n-q) ×(n-q)} e Î Â^p×p são matrizes de posto completo. Por simplicidade, pode-se escolher = I_n-q e = I_p.

Caso 3 Se B_d = 0 e D_w ¹ 0 utiliza-se um observador de ordem k = q - p + , no qual = I_n-q.

Caso 4 Se B_d ¹ 0 e D_w = 0 utiliza-se um observador de ordem reduzida k = q - p + , para o qual = I_p.

4.3 Discussão complementar

Como comentado anteriormente, na Hipótese 4.1, a condição (29) se reduz à condição de detectabilidade dos pólos finitos do sistema descritor se o desacoplamento de perturbações não é considerado (Castelan and Silva, 2005). Além disso, a verificação da condição (30), que corresponde à detectabilidade dos pólos infinitos do sistema representado por (E, A, C), é fundamental para que os valores calculados para as matrizes T, e satisfa cam a condição (11).

Em relação à condição adicional (31), a partir da qual define-se a ordem do observador, cabe ressaltar que ela pode ser relaxada para posto = n - q + p - d. Neste caso, a deficiência de posto representada por d deve ser compensada pelo acréscimo correspondente na ordem mínima do observador, em cada um dos casos descritos anteriormente.

Quanto à estrutura utilizada na definição do observador, destaca-se como original a utilização dos sinais auxiliares definidos em função das entradas de controle e das saídas medidas, conforme equações (7) e (8) respectivamente, o que se constitui na idéia chave para anular os termos decorrentes das perturba cões no erro de observação. Com isto, a ordem do observador é determinada, integradamente entre os casos 1 a 4 da Seção 4.2, em função das matrizes que descrevem a dinâmica das perturbações, o que é, de fato, um resultado singular.

A aplicação da equação de Sylvester generalizada (9), do observador com dinâmica descrita sob a forma (5) e a ordem reduzida do observador são características centrais deste trabalho encontradas também em trabalhos similares: i) não considerando o desacoplamento de perturbações (Verhaegen and dooren, 1986; Muller and Hou, 1993; Darouach and Boutayeb, 1995; Castelan and Silva, 2005), e ii) considerando o desacoplamento de perturbações (Yang and Tan, 1989; Darouach et al., 1996; Chu and Mehrmann, 1999b; Koenig and Mammar, 2002; Castelan et al., 2004).

É importante salientar que, sob o ponto de vista das hipóteses que indicam como as perturbações externas podem agir sobre os atuadores e sensores de um sistema descritor, os resultados propostos no presente trabalho são mais abrangentes que seus similares citados no item (ii) anterior. Esta abrangência refere-se ao fato de que no presente trabalho são previstos quatro casos para o par de matrizes de perturba cão B_w, D_w, relacionados ao final da seção anterior.

5 ALGORITMO E EXEMPLO

Com base nos resultados anteriores é proposto na seqüência um algoritmo para projeto de ODP de ordem mínima para sistemas descritores que satisfazem a Hipótese 4.1. Um exemplo numérico é apresentado na Subseção 5.2.

5.1 Algoritmo

Passo 1 Determinar a ordem mínima do observador, k_min, utilizando a equação (36) dentre os casos possíveis 1 a 4. Definir o espectro de autovalores desejados para o observador: s(H) = P = {p₁,p₂,..., }.

Passo 2 Para cada p_i, encontrar t_i Î ⁿ e z_iÎ ^p que verificam (32). As linhas da matriz T Î Â^k×n, denotadas por T_i, são formadas a partir dos vetores t_i como segue: se p_iÎ Â, então T_i = ; se p_i Î , considera-se p_i+1 = e T_i = e(), T_i+1 = mag(). De maneira similar, encontra-se a matriz Z. A matriz H correspondente, na forma de Jordan, é H = diag{D_i} onde D_i = p_i se p_i Î Â, e D_i = se p_i Î .

Passo 3 Encontrar as matrizes e que verificam B_d = 0 e D_w = 0, respectivamente.

Passo 4 Calcular as matrizes S, e utilizando a equação (22), com M^† = M^-1.

É importante lembrar que os vetores t_i, determinados no Passo 2, devem formar um conjunto de vetores linearmente independentes, além de garantir rank(TE) = k_min. Este tipo de requisito, comumente encontrado em técnicas de projeto via posicionamento de autoestrutura, pode ser em geral atendido escolhendo-se valores distintos para os elementos do conjunto P (veja, por exemplo, (Castelan, 2005)). Além disso, pode-se fazer uma verificação numérica do atendimento deste requisito ao longo dos cálculos realizados no Passo 2. Por construção, se posto(TE) = k_min, então a condição (11) é verificada e M, em (22), é uma matriz quadrada e inversível.

Deve-se destacar o fato de que a partir da definição de um ODP, é possível determinar novos parâmetros para o observador com o objetivo, por exemplo, de melhorar o condicionamento da matriz não-singular envolvida na condição (11), como proposto em (Silva, 2005), a partir dos graus de liberdade existentes na escolha das matrizes e . Além disso, no Caso (ii), em que B_d = 0 e D_w = 0, a síntese do ODP pode ser realizada via uma técnica de posicionamento regional de pólos em regiões do tipo LMI (Castelan and Silva, 2005; Silva, 2005).

5.2 Exemplo

Considerar o sistema descritor sob perturbações, que verifica a Hipótese 4.1, representado pelas matrizes:

Para ilustrar a utilização do algoritmo para os quatro casos enunciados na seção 4.2, supõe-se que B_d = LB_w pode assumir os valores = [1 1] ou = [0 0]. Da mesma forma, = [0.7382 0.1763] ou = [0 0]. Os autovalores utilizados foram: caso 1) {-7 ± 2i, -8 ± 3i}, caso 2) {-7 ± 2i}, casos 3 e 4) {-7, -8 ± 2i}. Para os dados acima, foram obtidas as matrizes do ODP, mostradas nas tabelas 1 e 2. As matrizes S, e podem ser calculadas diretamente da equação (22).

Visando ilustrar a validade do método proposto, considere o Caso 1, em que B_d ¹ 0, D_w ¹ 0. Além do ODP correspondente a este caso, relacionado nas Tabelas 1 e 2, é utilizado, para comparação, um outro observador com estrutura similar ao ODP, mas que atende apenas ao requisito de desacoplamento parcial da perturbação relativamente ao erro de estimação z(t) - TEx(t). Assim, este observador compartilha as mesmas matrizes H, T e Z do ODP, o que implica na satisfação da condição de desacoplamento TB_w + ZD_w = 0, mas as matrizes e são tais que B_d ¹ 0, D_w ¹ 0, ou seja, elas não satisfazem as condições para a obtenção do desacoplamento da perturbação relativamente ao erro de observação, x(t) - (t), no caso em consideração. Na obtenção das curvas de simulação temporal, as condições iniciais utilizadas para o sistema descritor e o observador foram x(0) = [0 1 0 0 0 0]' e z(0) = 0, respectivamente. Na figura 2 mostram-se as curvas dos erros de estimação e de observação, para o ODP e para o observador com desacoplamento parcial. Observa-se que os erros de estima cão e observação de um ODP convergem aos correspondentes estados nulos, conforme Figura 2 (a), enquanto apenas o erro de estimação converge ao estado nulo (o mesmo não ocorrendo com o erro de observação) para o observador com desacoplamento parcial, conforme Figura 2 (b).

6 CONCLUSÕES

Foram apresentadas condições para a existência de observadores de ordem reduzida para sistemas descritores, desacoplados de pertubações externas. Tais condi cões se basearam na solução de uma equação de Sylvester generalizada sob duas restrições adicionais. A primeira é uma restrição matricial que descreve o desacoplamento de perturbações do erro de estima cão, podendo ser resolvida conjuntamente com a equação de Sylvester. A segunda é uma restrição de posto necessária para a reconstrução das variáveis de estado, com desacoplamento de perturbações do erro de observação. Foram definidos sinais auxiliares em fun cão das entradas de controle e das saídas medidas, respectivamente, sendo utilizadas matrizes ortogonais às matrizes que descrevem a dinâmica das perturbações. Esta idéia foi o fundamento para anular termos decorrentes das perturbações, que estariam presentes no erro de observa cão. Uma realimentação de estados estimados foi considerada como ilustração de uma aplicação para o ODP proposto. Resultados numéricos foram apresentados para um exemplo de sistema descritor, onde foi utilizado o algoritmo de cálculo do ODP, elaborado com base nos resultados teóricos anteriores.

AGRADECIMENTOS:

Parte deste trabalho foi desenvolvido no LAAS-CNRS, Toulouse-França, durante doutorado sanduíche e pós-doutorado dos dois primeiros autores, respectivamente, com suporte financeiro da CAPES, Brasil. Este trabalho também contou com o suporte financeiro do CNPq, Brasil. Uma versão preliminar do presente artigo foi apresentada no CBA 2004, Gramado-R.S. Os autores agradecem ao editor associado e aos revisores por seus comentários, decorrentes do processo de revisão do artigo.

Artigo submetido em 03/08/2004

1a. Revisão em 03/08/2005

2a. Revisão em 21/09/2006

3a. Revisão em 21/05/2007

Aceito sob recomendação do Editor Associado

Prof. Liu Hsu

APÊNDICE: EXPRESSÃO PARA O SISTEMA EM MALHA FECHADA SOB UMA REALIMENTAÇÃO DE ESTADOS ESTIMADOS

Combinando (6) com (24), obtém-se

Substituindo (18) e (19) em (37), obtém-se

A equação (22) pode ser reescrita como

Substituindo (39) em (38), após breve simplificação, chega-se a

Substituindo (40) em (1) e realizando alguns passos algébricos, encontra-se:

A partir de (2), (15) e (41) é obtido o modelo do sistema em malha fechada dado por (25)-(26).

• Bhattacharyya, S. O. (1978). Observer design for linear systems with unknown inputs, IEEE Trans. on Automatic Control 23: 483484.
• Bolivar, A. R. (2001). Sur la synthèse de filtres de detection de defaillances, These doctorat de lUniversite Toulouse III, Spécialité automatique, LAAS du CNRS.
• Castelan, E. B. (2005). Estabilização de sistemas descritores por realimentação de saídas via subespaços invariantes (aceito), Controle e Automação - SBA (5).
• Castelan, E. B. and Silva, V. G. (2005). On the solution of a sylvester equation appearing in descriptor systems control theory, Systems & Control Letters, accepted for publication 54: 109117.
• Castelan, E. B., Silva, V. G., Tarbouriech, S. and Garcia, G. (2004). Disturbance decoupled minimal-order observers for a class of linear descriptor systems, Proc. of IFAC Symposium on Systems Structure and Control - SSSC04
• Chen, C. T. (1999). Linear System Theory and Design, Oxford University Press.
• Chu, D. and Mehrmann, V. (1999a). Disturbance decoupled observer design for descriptor systems, Systems & Control Letters 38: 3748.
• Chu, D. and Mehrmann, V. (1999b). Disturbance decoupled observer design for descriptor systems, Systems & Control Letters 38: 3748.
• Dai, L. (1989). Singular Control Systems. Lecture Notes in Control and Information Sciences, Springer-Verlag.
• Darouach,M. and Boutayeb,M. (1995). Design of observers for descriptor systems, IEEE Trans. on Automatic Control 40(7): 13231327.
• Darouach, M., Zasadzinski, M. and Hayar, M. (1996). Reduced-order observer design for descriptor systems with unknown inputs, IEEE Trans. on Automatic Control 41(7): 10681072.
• Hou, M. and Muller, P. C. (1992a). Design of observers for linear systems with unknown inputs, IEEE Trans. on Automatic Control 37: 871875.
• Hou, M. and Muller, P. C. (1992b). Observer for linear descriptor systems with unknown inputs control, Auto matisierungstechnic 40: 220227.
• Johnson, C. D. (1975). On observers for linear systems with unknown and inaccessible inputs, International Journal of Control 21: 825831.
• Koenig, D. and Mammar, S. (2002). Design of proportionalintegral observer for unknown input descriptor systems, IEEE Trans. on Automatic Control 47(12): 20572062.
• Liu, G. P. and Patton, R. J. (1998). Eigenstructure Assignment for Control System Design, JohnWiley and Sons.
• Muller, P. and Hou, M. (1993). On the observer design for descriptor systems, IEEE Trans. on Automatic Control 38(11): 16661671.
• Nagpal, K. M. and Khargonekar, P. P. (1991). Filtering and smoothing in an H_¥" setting, IEEE Trans. on Automatic Control 36: 152166.
• Paraskevopoulos, P. N., Koumboulis, F. N., Tzierakis, K. G. and Panagiotakis, G. E. (1992). Observers design for generalized state space systems with unknown inputs, Systems & Control Letters 18: 309321.
• Silva, V. G. (2005). Sobre síntese de observadores para sistemas descritores, Programa de pós graduação em engenharia elétrica, ppgeel, Universidade Federal de Santa Catarina-UFSC, Campus Universitário, Trindade, Florianópolis-SC. Tese de doutorado.
• Syrmos, V. L. (1993). Computational observer design techniques for linear systems with unknown inputs using the concept of transmission zeros, IEEE Trans. on Automatic Control 38(5): 790794.
• Varga, A. (1995). On stabilization methods of descriptor systems, Systems & Control Letters 24: 133138.
• Verhaegen, M. H. and dooren, P. V. (1986). A reduced observer for descriptor systems, Systems & Control Letters 8(5): 2937.
• Yaesh, I. and Shaked, U. (1992). Game theory approach to optimal linear state estimation and its relation to H_¥ norm estimation, IEEE Trans. on Automatic Control 37: 828831.
• Yang, C. W. and Tan, H. L. (1989). Observer design for singular systems with unknown inputs, International Journal of Control 49: 19371946.

1

O conjunto de pólos (finitos e infinitos) de um sistema descritor corresponde ao conjunto de autovalores generalizados do par (

E,A). Se o sistema é regular, o conjunto dos pólos finitos, denotado s(

E,A), é formado pelas raízes da equação característica det(l

E - A) = 0 (Castelan, 2005; Dai, 1989).

2

Imag(p

_i) denota a parte imaginária de p

_i

Datas de Publicação

• Publicação nesta coleção
  21 Jan 2008
• Data do Fascículo
  Dez 2007

Histórico

• Revisado
  03 Ago 2005
• Recebido
  03 Ago 2004
• Aceito
  21 Maio 2007