Open-access Função energia generalizada de controle para estabilização de sistemas não lineares

Resumos

Neste trabalho, o conceito de função de Lyapunov de controle é estendido em duas direções principais. Na primeira direção, o conceito de função energia de controle (FEC) é proposto com o objetivo de não somente projetar leis de controle estabilizantes para sistemas não lineares mas também com o objetivo de fornecer uma estimativa da região de estabilidade do sistema em malha fechada. A FLC garante estabilidade local de um certo ponto de equilíbrio do sistema em malha fechada, mas usualmente não fornece estimativasótimas da região de estabilidade. O conceito de FEC é mais exigente que o conceito de FLC. Além de garantir estabilidade local do equilíbrio, a FEC fornece informações globais a respeito dos conjuntos limites e estimativas ótimas das regiões de estabilidade. A segunda direção está relacionada com o fato de que, em geral, é difícil encontrar FLC ou FEC para muitos sistemas não lineares. De forma a minimizar este problema e aplicar õ conceito de CEF em uma classe maior de sistemas, propõe-se neste artigo o conceito de função energia generalizada de controle (FEGC). A FEGC permite o projeto de controladores mesmo quando a derivada da função energia é positiva em algumas regiões limitadas do espaço de estados.

Função de Lyapunov de controle; função energia; estabilização; sistemas não lineares


extended in two main directions. In the first direction, the concept of control energy function (CEF) is proposed not only to design stabilizing feedback control laws but also to provide an estimate of the stability region of the closed-loopsystem. The CLF guarantees local stability of the closedloop- system but usually cannot provide optimal estimates of its stability region. The CEF concept is more restrict than the CLF concept. Beyond the guarantee of local stability, the CEF provides global information about limit sets and optimal estimates of stability region. The second direction is related to the fact that, in general, it is difficult to find a CLF or a CEF for many nonlinear systems. In order to overcome this problem and apply the concept of CEF for a larger class of systems, the concept of generalized control energy function (GCEF) is proposed in this paper. The GCEF allows the design of nonlinear controllers even when the derivative of the energy function is positive in some bounded regions of the state space.

control Lyapunov function; energy function; stzcabilization; nonlinear systems


CONTROLE NÃO-LINEAR

Função energia generalizada de controle para estabilização de sistemas não lineares

Flávio H.J.R. SilvaI; Renato B.L. GuedesII; Luis F.C. AlbertoIII; Newton G. BretasIII

IUniversidade de Brasília (UnB) - Faculdade do Gama (FGA), Gama, DF, Brasil, fhsilva@unb.br

IIANEEL - Agência Nacional de Energia Elétrica, Brasília, DF, Brasil, rbguedes@aneel.gov.br

IIIDepartamento de Engenharia Elétrica Escola de Engenharia de São Carlos - USP São Carlos, SP, Brasil, luis@sel.eesc.usp.br, ngbretas@sel.eesc.usp.br

RESUMO

Neste trabalho, o conceito de função de Lyapunov de controle é estendido em duas direções principais. Na primeira direção, o conceito de função energia de controle (FEC) é proposto com o objetivo de não somente projetar leis de controle estabilizantes para sistemas não lineares mas também com o objetivo de fornecer uma estimativa da região de estabilidade do sistema em malha fechada. A FLC garante estabilidade local de um certo ponto de equilíbrio do sistema em malha fechada, mas usualmente não fornece estimativasótimas da região de estabilidade. O conceito de FEC é mais exigente que o conceito de FLC. Além de garantir estabilidade local do equilíbrio, a FEC fornece informações globais a respeito dos conjuntos limites e estimativas ótimas das regiões de estabilidade. A segunda direção está relacionada com o fato de que, em geral, é difícil encontrar FLC ou FEC para muitos sistemas não lineares. De forma a minimizar este problema e aplicar õ conceito de CEF em uma classe maior de sistemas, propõe-se neste artigo o conceito de função energia generalizada de controle (FEGC). A FEGC permite o projeto de controladores mesmo quando a derivada da função energia é positiva em algumas regiões limitadas do espaço de estados.

Palavras-chave: Função de Lyapunov de controle, função energia,estabilização, sistemas não lineares.

ABSTRACT

extended in two main directions. In the first direction, the concept of control energy function (CEF) is proposed not only to design stabilizing feedback control laws but also to provide an estimate of the stability region of the closed-loopsystem. The CLF guarantees local stability of the closedloop- system but usually cannot provide optimal estimates of its stability region. The CEF concept is more restrict than the CLF concept. Beyond the guarantee of local stability, the CEF provides global information about limit sets and optimal estimates of stability region. The second direction is related to the fact that, in general, it is difficult to find a CLF or a CEF for many nonlinear systems. In order to overcome this problem and apply the concept of CEF for a larger class of systems, the concept of generalized control energy function (GCEF) is proposed in this paper. The GCEF allows the design of nonlinear controllers even when the derivative of the energy function is positive in some bounded regions of the state space.

Keywords: control Lyapunov function, energy function, stzcabilization, nonlinear systems.

1 INTRODUÇÃO

Há anos, a teoria de Lyapunov vem sendo utilizada para análise de estabilidade de sistemas não lineares. Com o passar dos anos, alguns pesquisadores perceberam que as funções de Lyapunov poderiam ser úteis para o projeto de leis de controle estabilizantes em sistemas não lineares. Como resultado desta observação, no início da década de 80, o conceito de função de Lyapunov de controle (FLC) foi introduzido (Jurdjevic and Quinn, 1978; Artstein, 1983; Sontag, 1983; Sontag, 1989; Sontag, 1990; Clarke et al., 2000). A FLC permite, para uma classe de sistemas não lineares, obter leis de realimentação de controle que garantam a estabilidade assintótica do sistema em malha fechada. Entretanto, a FLC usualmente não permite obter estimativas ótimas da região de estabilidade do sistema em malha fechada. O conhecimento da região de estabilidade é importante, já que a maioria dos sistemas físicos não lineares não apresenta um ponto de equilíbrio assintoticamente globalmente estável.

Com o objetivo de não somente encontrar uma lei de realimentação de controle, mas também de estimar a região de estabilidade dos conjuntos assintoticamente estáveis, propõese a função energia de controle (FEC). Funções energia trazem informações globais a respeito dos conjuntos limites e podem ser utilizadas para se obter estimativas ótimas da região de estabilidade. Para isto, explora-se o conceito de função energia e a caracterização da região de estabilidade propostos em (Chiang et al., 1988).

O maior obstáculo na utilização de FLC ou de FEC para o projeto de controladores em sistemas não lineares é encontrar uma função de Lyapunov ou uma função Energia para o sistema em questão. Com o objetivo de minimizar este obstáculo, propõe-se neste artigo o conceito de função energia generalizada de controle (FEGC). Este conceito é inspirado na extensão do princípio de invariância de LaSalle apresentada em (Rodrigues et al., 2000; Rodrigues et al., 2001). A generalização consiste em permitir que a derivada da FEGC possua regiões limitadas, no espaço de estados, onde a derivada possa ser positiva. Assim como a FEC, a FEGC permite obter leis de realimentação para estabilização de sistemas não lineares e também estimar a região de estabilidade do atrator. A FEGC possui a FEC como caso particular e portanto uma classe maior de sistemas não lineares pode ser tratada com esta teoria. Em algumas situações particulares, explorandose as propriedades da FEGC e particularidades do campo vetorial, demonstra-se estabilidade assintótica de equilíbrios do sistema em malha fechada.

Este artigo está organizado da seguinte forma: a seção 2 contém definições e resultados preliminares da teoria de sistemas dinâmicos não lineares. Ainda na seção 2, o conceito de prolongamento de um conjunto, inicialmente proposto por LaSalle, assim como a relação deste conceito com estabilidade e invariância são introduzidos. Na seção 3, o conceito de função de Lyapunov de controle é revisto e um novo resultado (teorema 3.2) que explora a FLC para obter estimativas da região de estabilidade do sistema em malha fechada é demonstrado. Propõe-se, na seção 4, o conceito de função energia de controle, e novos resultados (teoremas 4.2 e 4.3) que exploram a existência da função energia para obter informações a respeito dos conjuntos limites e estimativas da região de estabilidade são demonstrados. A seção 5 contém a principal contribuição deste trabalho. Nesta seção os conceitos de função energia generalizada e função energia generalizada de controle são propostos. Estudam-se as implicações da existência de funções energia generalizadas no comportamento e localização dos conjuntos limites assim como a obtenção de estimativas da região de estabilidade de conjuntos atrativos. Na seção 6, dois exemplos ilustram os conceitos e os resultados propostos e as conclusões encerram o artigo na seção 7.

2 SISTEMAS DINÂMICOS NÃO LINEARES

Nesta seção, apresenta-se uma breve revisão de algumas definições e resultados relacionados a sistemas autônomos não lineares que serão úteis para o desenvolvimento deste trabalho. Considere o sistema autônomo

onde x(t) ∈ n é o vetor de estados e f:nn é uma função de classe C1, condição esta suficiente para garantir existência e unicidade das soluções de (1). A solução de (1) iniciando em xo no instante de tempo t = 0 é denotada por φ(·, xo):n. O ponto x* é um ponto de equilíbrio de (1) se f(x*) = 0. Neste trabalho, considere a seguinte suposição sobre os pontos de equilíbrio.

(S1) Os pontos de equilíbrio de (1) são isolados.

O conceito de invariância tem um papel fundamental neste trabalho. Um conjunto n é positivamente invariante (negativamente invariante) com relação a (1) se, para todo xo, φ(t,xo) ∈ para todo t+ (t-). Um conjunto é invariante com relação a (1) se é positivamente e negativamente invariante.

Um ponto pn é um ponto limite de φ(t, xo), ou simplesmente de xo, se existir uma sequência {tn} → ∞ quando n → ∞ tal que {φ(tn, xo)} → p quando n → ∞. O conjunto de todos os pontos limites de xo é denotado por ω(xo). Se a solução φ(t, xo) for limitada para t > 0, então, além de fechado e invariante, o conjunto ω(xo) é não vazio, compacto, conexo e é o menor conjunto fechado do qual a solução φ(t, xo) se aproxima quando t → ∞, {ver pág. 29 de (LaSalle, 1976) para maiores detalhes}.

Um conjunto fechado e invariante H é atrativo se existir vizinhança U de H tal que xoU implica em d(φ(t,xo),H) → 0 quando t → ∞.

Conjuntos atrativos de sistemas não lineares não são em geral globalmente atrativos. A região de estabilidade A(H) de um certo conjunto atrativo H é o conjunto de condições iniciais cujas trajetórias tendem para o conjunto atrativo quando t → ∞, isto é, A(H): = {xn : φ(t, x) → H quando t → ∞}. A região de estabilidade de um conjunto atrativo conexo é um conjunto aberto, conexo e invariante. A fronteira topológica de A(H) é um conjunto fechado e invariante e é denotada por ∂A(H).

Um conjunto compacto H é estável se, dada uma vizinhança U de H, existe vizinhança W de H tal que xoW implica φ(t, xo) ∈ U para todo t > 0. Conjuntos atrativos e estáveis são ditos assintoticamente estáveis.

Em algumas situações, encontramos conjuntos invariantes e até mesmo atrativos que não são necessariamente estáveis. O conceito de prolongamento positivo de um conjunto H, introduzido por LaSalle (LaSalle, 1976), recupera, num certo sentido, o conjunto de pontos que são "visitados" por órbitas que se iniciam suficientemente próximas de H.

Definição 2.1 (Prolongamento)(LaSalle 1976) Um ponto zn pertence ao prolongamento positivo do conjunto H se existir sequência xn com xnyH e seqüência tn ∈ [0, τmax(xn)), onde τmax(xn) é o máximo tempo no qual a solução φ(t, xn) está definida, tal que φ(tn, xn) → z quando n → ∞.

É óbvio que H. De fato, se zH, podemos escolher uma sequência qualquer xny = z e tn = 0 para todo n. Portanto, pela continuidade das soluções com relação às condições iniciais, φ(0, xn) → z quando n → ∞ e portanto z. Os próximos resultados relacionam o conceito de prolongamento com estabilidade e invariância. Eles serão utilizados nas próximas seções para mostrar a existência de conjuntos assintoticamemte estáveis. O próximo teorema oferece uma maneira de verificar se um conjunto invariante é estável via análise de seu prolongamento.

Teorema 2.1 (Prolongamento e Estabilidade)(LaSalle 1976) Se H é um conjunto compacto positivamente invariante, então H é estável se e só se = H.

Embora o teorema 2.1 tenha sido proposto por LaSalle, a demonstração do resultado anterior não foi incluída em (LaSalle, 1976). A seguir apresenta-se uma demonstração para este teorema.

Demonstração: Suponha que z e zH. Como H é um conjunto compacto e portanto fechado, existe número real γ tal que d(z, H) > γ > 0. Portanto, podemos escolher vizinhança U de H como sendo U := {xn : d(x, H) < γ/2}. Por outro lado, como z, existem seqüências xnyH e tn ∈ [0, τmax(xn)) tal que φ(tn, xn) → z quando n → ∞. Em outras palavras, dado ε > 0, existe N > 0 tal que d(xn, y) < ε e d(φ(tn, xn), z) < ε para todo n > N. Como ε pode ser escolhido arbitrariamente pequeno, para qualquer vizinhança W de H, podemos encontar um ponto xnW tal que φ(tn, xn) ∉ U. Logo H não é estável. Provamos portanto que se H é estável, então para todo z, zH. Em outras palavras, se H é estável, então H. Como a inclusão Hé sempre verdadeira, tem-se que H estável implica = H.

Suponha que H não seja estável. Então, para alguma vizinhança limitada U de H não existe vizinhança W de H tal que xW implica φ(t, x) ∈ U para todo t > 0. Isto significa que arbitrariamente próximo de H existe xi tal que φ(t, xi) ∉ U para algum t > 0. Em particular, podemos escolher uma sequência xiyH e ti como sendo inft> 0{t:φ(t,xi) ∉ U}. Então, a seqüência φ(ti,xi) é, por construção, limitada e portanto, possui subsequência convergente, isto é, φ z. Por definição, z e zH. Como consequência, se H, então H é estável. Como a inclusão H é sempre verdadeira, tem-se que = H implica H estável.

O próximo teorema oferece condições para garantir que o prolongamento de um conjunto invariante também é invariante.

Teorema 2.2 (Prolongamento e Invariância) (LaSalle 1976) Seja H um conjunto fechado e invariante contido em um conjunto aberto G limitado e positivamente invariante, então é invariante.

O teorema 2.2 também foi proposto por LaSalle, mas a demonstração não foi incluída em (LaSalle, 1976). A seguir apresenta-se uma demonstração para este teorema.

Demonstração: Seja U uma vizinhança de H totalmente contida em G, ou seja, HUG. Como é positivamente invariante, então . Seja z. Vamos mostrar que φ(t,z) ∈ para todo t. Como z, existem seqüências xnyH e tn ∈ [0,τmax(xn)) tal que φ(tn,xn)→ z quando t→ ∞. Para n suficientemente grande, xnG. Como G é positivamente invariante, o intervalo maximal de definição da solução φ(t,xn) contém o intervalo (0,+∞). Em particular τmax(xn) = + ∞. Da continuidade das soluções com relação as condições iniciais, tem-se que φ(t, φ(tn,xn)) → φ(t, z) quando n→ ∞. Da propriedade de fluxos, φ(t, φ(tn,xn)) = φ(t + tn, xn) = φ(tn, φ(t, xn)). Como H é um conjunto invariante, φ(t, xn) → ∞ (t,y) ∈ H quando n→ ∞ para todo t. Portanto, φ(t,z) ∈ para todo t.

3 FUNÇÃO DE LYAPUNOV DE CONTROLE

Considere o sistema (1) e suponha que a origem seja um ponto de equilíbrio deste sistema, isto é, f(0) = 0.

Definição 3.1 (Função de Lyapunov) A função :n, de classe C1, é uma função de Lyapunov de (1) se as seguintes condições forem satisfeitas: (i) (x) é localmente definida positiva; e (ii) (x) é localmente definida negativa.

A existência de uma função de Lyapunov é uma condição suficiente para a estabilidade assintótica da origem de (1) (Khalil, 1996).

Considere agora o sistema

onde xn é o vetor de estados, um é a entrada de controle e F : n × mn é uma função de classe C1. Suponha que a origem seja um ponto de equilíbrio do sistema em malha aberta (u = 0) isto é, F(0,0) = 0. O objetivo é obter uma lei de realimentação de controle u = h(x), definida na vizinhança da origem, tal que a origem do sistema em malha fechada

seja assintoticamente estável. Uma condição suficiente para garantir a estabilidade da origem do sistema em malha fechada é a existência de uma função de Lyapunov de controle (FLC) (Jurdjevic and Quinn, 1978; Artstein, 1983; Sontag, 1990).

Definição 3.2 (Função de Lyapunov de Controle) Uma função : n × m, de classe C1, é uma FLC de (2) se existir uma lei de realimentação de controle u = h(x), de classe C1, com h(0) = 0, tal que a função (x):n, definida por (x) = (x,h(x)), seja uma função de Lyapunov de (3).

O próximo teorema garante que a existência de uma FLC implica na estabilidade da origem do sistema (3) (Artstein, 1983; Sontag, 1983; Clarke et al., 2000).

Teorema 3.1 (Estabilidade do Sistema em Malha Fechada) Seja uma FLC de (2), então existe uma lei de realimentação de controle, tal que a origem do sistema em malha fechada (3) é um ponto de equilíbrio assintoticamente estável.

Demonstração: A demonstração é trivial. Pode ser obtida a partir da definição 3 e dos resultados de Lyapunov sobre estabilidade (Khalil, 1996).

O teorema 3.1 estabelece que a existência de uma FLC é uma condição suficiente para assegurar a estabilização da origem. Sontag (1983) mostra que a existência de uma FLC também é uma condição necessária, mas para isto, ele relaxa as condições sobre o funcional de Lyapunov, isto é, ele mostra que a origem do sistema é assintoticamente controlável se e somente se existe um funcional contínuo (não necessariamente de classe

1) definido positivo cuja derivada pode ser feita negativa por uma escolha apropriada de lei de realimentação. Entretanto, este resultado é existencial e, em geral, é muito difícil explicitar a FLC em termos de funções elementares conhecidas. Neste artigo, apenas a suficiência da existência dos funcionais é explorada. Particularmente, partindo-se do princípio da existência dos funcionais, estuda-se a relação entre os funcionais e a região de estabilidade.

A FLC pode fornecer alguma informação a respeito da região de estabilidade do sistema em malha fechada. Mais precisamente, prova-se neste artigo o seguinte resultado:

Teorema 3.2 (Região de Estabilidade do Sistema em Malha Fechada) Seja V uma FLC do sistema (3) e Ω um conjunto aberto contendo a origem. Se a função W(x) = V(x,h(x)) satisfaz as seguintes condições: (a) W(0) = 0 e W(x) > 0 para todo x ∈ Ω - {0}; (b) (0) = 0 e < 0 para todo x ∈ Ω - {0} então: (i) a origem é um ponto de equilíbrio assintoticamente estável do sistema em malha fechada, (ii) existe constante L tal que a componente conexa Sc(L) do conjunto de nível {x ∈ n:W(x) < L} que contém a origem está contida em Ω e (iii) Sc(L) está contida na região de estabilidade A(0).

Demonstração: A estabilidade assintótica da origem é uma conseqüência direta do teorema 3.1. Utilizando a condição (a) e a continuidade da função W, é fácil verificar que existe constante L tal que a componente conexa Sc(L) que contém a origem é limitada e está totalmente contida em Ω. Do princípio de invariância de LaSalle (Alberto, 2000), demonstra-se que Sc(L) é positivamente invariante e que toda órbita iniciando em Sc(L) está definida para todo tempo t > 0 e tende para o maior conjunto invariante contido no conjunto C = {xc(L):(x) = 0}. Utilizando a condição (b), prova-se que C = {0}. Portanto Sc(L) ⊂ A(0).

O teorema anterior, embora forneça uma estimativa da região de estabilidade da origem do sistema em malha fechada, ele não oferece indicativos de como encontrar estimativas ótimas da região de estabilidade. Em outras palavras, é difícil encontrar uma FLC que satisfaça todas as condições do teorema 3.2 em um aberto Ω grande o suficiente para que possamos obter estimativas relevantes da região de estabilidade. Além disto, o teorema 3 não oferece procedimento sistemático para uma escolha ótima do nível L. Na próxima seção, o conceito de função energia de controle será definido e algumas condições do teorema anterior serão relaxadas de forma a obter-se estimativas ótimas da região de estabilidade do sistema em malha fechada.

4 FUNÇÃO ENERGIA DE CONTROLE

Considere o sistema (1). Seja E o conjunto de todos os pontos de equilíbrio de (1), isto é, E: = {xn:f(x) = 0}.

Definição 4.1 (Função Energia)(Chiang et al., 1988) A função : n , de classe C1, é uma função energia de (1) se as seguintes condições forem satisfeitas: (i) (x) < 0 para qualquer xn; (ii) se xoE, então o conjunto {t+ : (φ(t,xo)) < 0} não possui medida finita em ; (iii) se |(φ(t,xo))| é limitada para t+, então a trajetória φ(t,xo) é também limitada para t+.

As condições (i) e (ii) implicam que a energia do sistema é estritamente decrescente ao longo de qualquer trajetória não trivial. A definição de função energia dada em (Chiang et al., 1988) exige que o conjunto {t: (φ(t, xo)) = 0} possua medida nula em se xoE. A condição (ii) da definição anterior é mais geral e contém a definição dada em (Chiang et al., 1988) como um caso particular, ou seja, a exigência de que o conjunto {t: (φ(t,xo)) = 0} possua medida nula em , se xoE, dada em (Chiang et al., 1988) é uma condição suficiente para a satisfação da condição (ii) da definição 4.1. A função não precisa ser própria1 mas a condição (iii) garante que sendo limitada ao longo da trajetória, então a trajetória φ(t) é limitada para t > 0. Maiores esclarecimentos podem ser obtidos em (Chiang et al., 1988).

A existência de uma função energia possui implicações importantes a respeito dos conjuntos limites. O resultado a seguir é uma adaptação dos resultados obtidos em (Chiang et al., 1988).

Teorema 4.1 (Função Energia e Conjuntos Limites) Seja φ(t, xo) uma solução de (1) e uma função energia de (1) tal que |( φ(t, xo))| é limitada para t > 0. Então, o conjunto limite ω(xo) é composto por pontos de equilíbrio.

Demonstração: Sendo |(φ(t,xo))| limitada para t> 0, a condição (iii) da definição 4.1 implica que φ(t,xo) é uma solução limitada de (1), para t > 0. Portanto o conjunto ω-limite de φ(t,xo) é não vazio. A condição (i) da definição 4.1 implica que (φ(t,xo)) é uma função decrescente em t. Portanto, existe α tal que (φ(t,xo)) → α quando t → ∞. Seja pω (xo). Então, pela continuidade de , (p) = αpω(xo). A invariância de ω(xo) garante que (p) = 0 para todo p ∈ Ω. A condição (ii) da definição 4.1 implica que p é um ponto de equilíbrio.

Corolário 4.2 Se a suposição (S1) for satisfeita, for uma função energia de (1) e φ(t, xo) uma solução de (1) tal que |(φ(t, xo))| seja limitada para t > 0, então o conjunto ω(xo) é composto por um único ponto de equilíbrio.

Demonstração: O teorema 4 garante que todo ponto p∈ ω(xo) é um ponto de equilíbrio. Como os pontos de equilíbrio são isolados por hipótese e o conjunto ω(xo) é um conjunto conexo, então ω(xo) é constituído por um único ponto de equilíbrio.

O conceito de função energia, apresentado anteriormente, pode ser utilizado com a finalidade de projetar leis de realimentação. Para isto, define-se o conceito de função energia de controle (FEC).

Definição 4.3 (Função Energia de Controle) Uma função : n × m, de classe C1, é uma FEC de (2) se existir uma lei de realimentação de controle u = h(x), de classe C1, tal que = (x, h(x)) seja uma Função Energia de (3).

O próximo teorema garante que a existência de uma FEC em conjunto com o conceito de prolongamento de conjuntos asseguram a estabilidade do sistema (3) e fornece também uma estimativa de sua região de estabilidade.

Teorema 4.2 (Estabilidade e Região de Estabilidade do Sistema em Malha Fechada) Seja : n × m uma FEC de (2). Se existir uma componente conexa limitada Sc(L) do conjunto {xn : (x, h(x)) < L}, tal que xs seja o único ponto de equilíbrio contido no interior de Sc(L), então xs é um ponto de equilíbrio assintoticamente estável do sistema em malha fechada (3). Além disso, Sc(L) é uma estimativa da região de estabilidade.

Demonstração: Vamos demonstrar, em primeiro lugar, que o equilíbrio é atrativo. Seja (x) uma FEC de (2). Então, de acordo com a definição 4.3, existe uma lei de realimentação de controle u = h(x), de classe C1, tal que (x) = (x,h(x)) é uma função energia de (3). Por hipótese, Sc(L) é uma componente conexa limitada do conjunto {xn :(x,h(x)) < L} e xs é o único ponto de equilíbrio contido em Sc(L). O princípio de invariância de LaSalle (Alberto, 2000) garante que Sc(L) é um conjunto positivamente invariante, isto é, todas as trajetórias de (3) partindo do interior de Sc(L) são limitadas e convergem para o maior conjunto invariante contido em Sc(L). O Teorema 4.1 garante que os conjuntos limites são compostos somente por pontos de equilíbrio. Sendo xs o único ponto de equilíbrio contido em Sc(L), então toda trajetória partindo em Sc(L) converge para xs quando t → ∞. Como consequência, este ponto de equilíbrio é atrativo. Vamos mostrar que o equilíbrio xs é estável. O ponto de equilíbrio xs é um conjunto compacto e invariante, logo, em acordo com o teorema 2.1, H = {xs} é estável se e só se = H. O conjunto H = {xs} é um conjunto fechado e invariante contido em Sc(L) que, por sua vez, é um conjunto aberto, limitado e positivamente invariante, portanto, segue do teorema 2.2 que é um conjunto invariante. Suponha existir z tal que zxs. Então, existem seqüências xn xs e tn ∈ [0,τmax(xn)) tal que φ(tn,xn)→ z quando n→ ∞. Portanto, segue da suposição (i) da definição 4.1 que V(xn) > V(φ(tn,xn)). Tomando o limite quando n→ ∞, conlcui-se que V(xs) > V(z). Por outro lado, xs é um ponto de equilíbrio isolado que pertence ao conjunto aberto Sc(L). Logo, existe ε > 0 suficientemente pequeno tal que xsSc(L-ε): = {xSc(L):V(x) < L-ε}. Por construção, Sc(L-ε) é um conjunto positivamente invariante e portanto . Conclui-se portanto que se z, então zSc(L) e φ(t,z)→ xs quando t→ ∞. Novamente, utilizando a suposição (i) da definição 4.1, conclui-se que V(z) > V(xs). Isto significa que V é constante e igual a V(xs) em . Como é invariante e z não é ponto de equilíbrio, tem-se que V(φ(t,z)) = V(xs) para todo t o que implica que (φ(t,z)) = 0 para todo t violando a propriedade (ii) da definição 4.1. Portanto não existe z tal que zxs, o que implica em = H.

A vantagem em encontrar uma FEC quando comparada a FLC é que podemos utilizar a teoria de caracterização da fronteira da região de estabilidade para obter estimativas ótimas da região de estabilidade. Mais precisamente, tem-se o seguinte resultado:

Teorema 4.3 (Estimativa Ótima da Região de Estabilidade do Sistema em Malha Fechada) Seja :n× m uma FEC e seja xs um ponto de equilíbrio assintoticamente estável do sistema em malha fechada. Defina o nível de energia

Então:

(i) a componente conexa Sc(L) do conjunto de nível {xn:W(x) = V(x,h(x)) < L} contendo xs está contida na região de estabilidade A(xs) do sistema em malha fechada,

(ii) a componente conexa Sc(M) do conjunto de nível {xn:W(x)) < M} contendo xs tem interseção não vazia com o complemento da região de estabilidade Ac(xs) do sistema em malha fechada para qualquer M > L.

Demonstração: Se é uma FEC, então existe uma lei de realimentação u = h(x) tal que é uma função energia do sistema em malha fechada. A região de estabilidade A(xs) do ponto de equilíbrio assintoticamente estável xs do sistema em malha fechada é um conjunto aberto, conexo e invariante. A fronteira de todo conjunto invariante também é invariante. Portanto ∂A(xs) é um conjunto invariante. Da continuidade da função e da condição (i) da definição 4.1, prova-se que W(xs) é um limitante inferior da função na fronteira ∂A(xs) da região de estabilidade, ou seja, W(x) > W(xs) para qualquer x ∈ ∂A(xs). Em seguida, vamos provar que toda solução na fronteira tende para um ponto de equilíbrio na fronteira. Seja xo ∈ ∂A(xs) um ponto genérico na fronteira da região de estabilidade. Como ∂A(xs) é invariante, φ(t,xo) ∈ ∂A(xs) para todo t. Além disto, da condição (i) da definição 4.1, tem-se que W(xs) < W(φ(t,xo)) < W(xo) para todo t > 0. Portanto, como conseqüência do corolário 4.2, conclui-se que todas trajetórias na fronteira da região de estabilidade tendem para um equilíbrio na fronteira ∂A(xs). Para completar a prova, considere xj ∈ ∂A(xs) como sendo o equilíbrio com menor energia na fronteira da região de estabilidade, ou seja, L = W(xj). Como conseqüência das conclusões anteriores, W(x) > W(xj) para qualquer x ∈ ∂A(xs). Isto significa que qualquer caminho conexo saindo de xs e cruzando a fronteira da região de estabilidade tem que antes cruzar a fronteira do conjunto Sc(L) antes de sair da região de estabilidade. Portanto Sc(L) ⊂ A(xs). Isto prova (i). Para provar (ii), seja M > L. Então existem pontos na fronteira que pertencem a Sc(M) e portanto Sc(M) ∩ Ac(xs) é não vazio.

O teorema 4.3 utiliza uma curva de nível da função energia para obter uma estimativa ótima da região de estabilidade. Deste teorema conclui-se que o valor da energia do ponto de equilíbrio instável de menor energia na fronteira da região de estabilidade oferece a melhor estimativa da região de estabilidade na forma de um conjunto de nível da função . Em outras palavras Sc(L) é o maior conjunto de nível da função que está contido na região de estabilidade do sistema em malha fechada. Além disto, o teorema 4.3 oferece um procedimento sistemático para o cálculo do nível L.

Embora a FEC forneça estimativas ótimas da região de estabilidade, a FEC, assim como a FLC, exige que se tenha ( φ(t)) < 0. Em vários sistemas dinâmicos, é difícil encontrar função V satisfazendo esta condição. Com o próposito de relaxar tal condição, na próxima seção será apresentado o conceito de função energia generalizada de controle.

5 FUNÇÃO ENERGIA GENERALIZADA DE CONTROLE

Nesta seção, generaliza-se a teoria da seção anterior permitindo que a função possua, no espaço de estados, regiões limitadas de derivada positiva. Para isto define-se o conceito de função energia generalizada de controle (FEGC).

Definição 5.1 (Função Energia Generalizada) Uma função : n, de classe C1, é uma função energia generalizada (FEG) de (1) se: (i) todo conjunto compacto intercepta um número finito de componentes conexas limitadas e isoladas i's do conjunto : = {xn:(x) > 0} e (ii) supt> 0|(φ(t,xo))| < ∞ implica φ(t,xo) é limitado para t > 0.

As condições exigidas para que uma função seja uma FEG são menos restritivas do que as condições da função energia, ainda assim, não é trivial encontrar uma função energia generalizada para um sistema não-linear. Condições que garantem a existência de uma função energia generalizada para uma classe de sistemas não-lineares com aplicações em sistemas elétricos de potência foram estudadas em (Alberto, 2006). O artigo (Martins et al., 2006) estuda condições semelhantes para uma classe de sistemas não lineares que é uma generalização dos sistemas na forma de Lur'e. Neste artigo, na seção de exemplos, propriedades particulares dos sistemas não-lineares em estudo serão utilizadas para provar a satisfação de tais condições.

É interessante ressaltar que a condição (i) da definição 5 poderia ser substituída pela condição de que o conjunto C é composto por um número finito de componentes conexas compactas e isoladas; entretanto, esta alternativa é menos geral e restringe a aplicação dos resultados para uma classe de problemas que satisfazem a condição (i) da definição 5.1, mas contêm um número infinito de componentes compactas e isoladas. Este é o caso do exemplo de sistemas elétricos de potência discutido em (Alberto, 2006).

A existência de uma função energia, garante, conforme demonstrado no teorema 4.1, que os conjuntos limites são constituídos exclusivamente por pontos de equilíbrio. Embora a função energia generalizada não possua a derivada ao longo das trajetórias semi-definida negativa, ainda assim, informações importantes a respeito dos conjuntos limites podem ser obtidas conforme mostra o teorema a seguir. Estas informações, em conjunto com o conceito de prolongamento de conjuntos, serão exploradas no teorema 5.2 para provar a existência e estabilidade de um conjunto atrativo.

Teorema 5.1 (Função Energia Generalizada e Conjuntos Limites) Considere que o sistema (1) possua uma FEG e que φ(t,xo) seja uma solução do sistema (1) tal que |(φ(t,xo))| é limitado para t > 0. Então φ(t,xo) é uma solução limitada, seu conjunto ω-limite é não vazio e intercepta pelo menos uma componente conexa i do conjunto .

Demonstração: Seja (x) uma FEG de (1). De acordo com a condição (ii) da definição 5.1, sendo |(φ(t,xo))| limitada, para t > 0, tem-se que φ(t,xo) é uma solução limitada de (1). Portanto, o conjunto ω-limite é não vazio. Suponha que φ(t,xo) não intercepte o conjunto para todo t > 0. Então usando argumentos similares aqueles utilizados na demonstração do princípio de invariância de LaSalle (Khalil, 1996), demonstra-se que φ(t,xo) tende para o maior conjunto invariante contido em {xn : (x) = 0} quando t→ ∞. Neste caso, a conexidade do conjunto ω-limite e o fato das componentes i's serem isoladas garante que existe inteiro j tal que φ(t,xo) tende para o maior conjunto invariante contido em {xn : (x) = 0}∩j. Suponha que φ(t,xo) intercepte . Então existe um tempo t1 e uma componente i1 tal que φ(t1,xo) ∈ i1. Se a solução permanecer em i1 para todo t > t1, então demonstra-se que φ(t,xo) tende para o maior conjunto invariante contido em {xn : (x) = 0}∩i1 quando t→ ∞. Se por outro lado, a solução abandonar Ci1, ou seja, existir tal que φ(t,xo) ∈i1 para t1 < t < e φ(t,xo) ∉i1 para t ligeiramente maior que , então duas possibilidades podem ocorrer: (i) a solução não intercepta para todo t > ou (ii) ∃ t2 e componente i2 tal que φ(t2,xo) ∈i2. No caso (i), a demonstração segue da primeira parte. No caso (ii) novamente temos duas possibilidades, ou seja, ou φ(t,xo) ∈ 2 para todo t > t2 ou existe t3 > tal que φ(t3,xo) ∈ i3. Se houver uma sequência infinita de tempos {tr} com tr → + ∞ tal que φ(tr,xo) ∈ ir, então, como o número de componentes ir com interseção não vazia com conjuntos compactos é finito, existe pelo menos uma componente j que é visitada infinitas vezes. Como toda sequência infinita num conjunto compacto possui subsequência convergente, então j possui pelo menos um ponto pertencente ao conjunto ω-limite de φ(t,xo). Logo a intersecção do conjunto ω-limite com a componente j é não vazia.

Como se exige menos da FEG se comparada a função energia, informações mais fracas a respeito dos conjuntos limites são obtidas no teorema anterior. Em compensação, o teorema anterior permite que os conjuntos limites sejam mais complexos incluindo ciclos limites e comportamentos caóticos.

É interessante observar que a interseção do conjunto ω-limite com o conjunto não garante que o conjunto limite esteja contido em . Em verdade, as soluções podem entrar e sair da região onde a derivada é positiva. Apesar disto, o teorema 5.2, que será apresentado a seguir, oferece condições sobre a FEG e explora o conceito de prolongamento de um conjunto para provar a existência de um conjunto assintoticamente estável e obter estimativas da região de estabilidade deste conjunto. Mostra-se que este conjunto assintoticamente estável intercepta o conjunto C embora não esteja necessariamente contido em C, ou seja, a solução de "regime permanente" pode ficar entrando e saindo do conjunto C indefinidamente. Para maiores eslcarecimentos a respeito destas idéias vide (Rodrigues et al., 2000) e (Rodrigues et al., 2001).

Embora o teorema 5.1 não possa ser utilizado para provar estabilidade do conjunto limite, ele em conjunto com o conceito de prolongamento são explorados para demonstrar a existência de um conjunto assintoticamente estável. Com esta finalidade, propõe-se a função energia generalizada de controle (FEGC).

Definição 5.2 (Função Energia Generalizada de Controle) Uma função : n ×mde classe C1, é uma FEGC de (2), se existir uma lei de realimentação de controle u = h(x), de classe C1, tal que (x) = (x,h(x)) seja uma FEG do sistema em malha fechada (3).

Teorema 9 (Estabilidade e Região de Estabilidade do Sistema em Malha Fechada) Considere que o sistema (2) possua uma FEGC . Seja um número real tal que a componente conexa Sc(L) de {xn :(x) < L} seja limitada. Suponha que . Então, Sc(l): = {xSc(L):(x) < l} contém um conjunto H invariante e assintoticamente estável e Sc(L) é uma estimativa da região de estabilidade de H. Além disto, H tem interseção não vazia com o conjunto CSc(l).

Demonstração: Seja (x) uma FEGC de (2). Então, existe uma lei de realimentação de estados u = h(x), de classe C1, tal que (x) = (x,h(x)) é uma FEG de (3). Por hipótese Sc(L) é uma componente conexa e limitada de {xn :(x)) < L}. Se então CSc(l) ⊂ Sc(L). Pela extensão do princípio de invariância (Rodrigues et al., 2000; Rodrigues et al., 2001), Sc(L) e Sc(l) são conjuntos positivamente invariantes e todas as trajetórias de (3) partindo do interior de Sc(L) convergem para o maior conjunto fechado e invariante H contido em Sc(l). O teorema 5.1 garante que H tem interseção não vazia com o conjunto CSc(l). Desde que d(φ(t,xo),H) → 0 quando t → + ∞ para todo xoSc(L), então H é um conjunto atrativo e Sc(L) é uma estimativa da região de estabilidade de H. Falta demonstrar que o conjunto H é estável. Como H é um conjunto compacto e invariante, então, pelo teorema 2.1, H é estável se e só se = H. Provaremos a seguir que Sc(l). Como H é um conjunto fechado e invariante contido em Sc(L) que é um conjunto aberto, limitado e positivamente invariante, então, pelo teorema 2.2, é um conjunto invariante. Por outro lado, Sc(l) é um conjunto estável, compacto e positivamente invariante. Então, novamente pelo teorema 2.1, c(l) = Sc(l). Como HSc(l), então c(l) = Sc(l). Mas é um conjunto invariante contido em Sc(l). Portanto H o que implica = H.

O teorema proposto anteriormente permite que se tenha, em algumas regiões limitadas do espaço de estados, (φ(t)) > 0. Esta flexibilidade no conceito de FEGC permite o tratamento de problemas que não apresentam função energia ou em que a função energia seja de difícil obtenção. Mesmo sob esta generalização, demonstra-se a existência de um conjunto assintoticamente estável e obtêm-se estimativas de sua localização e de sua região de estabilidade. Estimativas ótimas da região de estabilidade na forma de curvas de nível da FEGC também podem ser encontradas. Maiores detalhes a este respeito podem ser encontrados em (Alberto and Chiang, 2008).

6 EXEMPLOS

6.1 Exemplo 1

Considere o sistema não-linear:

onde kd = 1, g = 6 e u é a entrada de controle. Este sistema tem origem no problema de controle de velocidade de um motor de corrente contínua via tensão de campo quando não linearidades referentes ao efeito de reação de armadura são levados em consideração (Fallside and Patel, 1965).

Linearizando este sistema, verifica-se que a origem é um ponto de equilíbrio assintoticamente estável do sistema em malha aberta para todo kd > 0. Entretanto, para g > 4, o sistema apresenta 3 pontos de equilíbrio indicando que a origem não é um ponto de equilíbrio globalmente assintoticamente estável do sistema em malha aberta.

Neste exemplo, utilizaremos o conceito de função energia generalizada e exploraremos as peculiaridades do sistema (4) para projetar uma lei de realimentação de estados u = h(x1, x2) que assegure estabilidade global para a origem. Com este propósito, considere a seguinte candidata a função energia generalizada de controle:

Escolhendo e k5 = 2k2 e calculando a derivada da função V ao longo das trajetórias de (4) obtém-se:

Podemos escolher as constantes da função energia de tal forma que 2k2kd > k3. Neste exemplo consideraremos k3 = Kd e k2 = 1. Com esta escolha o termo quadrático em x2 da derivada fica definido negativo e não é difícil verificar que V é uma função energia para o sistema (4) em malha aberta se g < 4. Entretanto, para g > 4 não conseguimos provar a condição (i) da definição 4.1 que exige < 0.

A seguir, iremos mostrar que V é uma função energia generalizada de controle para o sistema (4). Para isto, escolha u = h(x1,x2) = . Observe que h(0,0) = 0 e portanto a origem é um ponto de equilíbrio do sistema em malha fechada. Em malha fechada é fácil verificar que a origem é o único ponto de equilíbrio deste sistema. Para mostrar que V é uma função energia generalizada de controle, precisamos mostrar que W(x) = V(x,h(x)) é uma função energia generalizada do sistema em malha fechada. Substituindo a lei de realimentação u na expressão de obtém-se:

Escolhendo as constantes de tal forma que 2k2kd-k3 > 0, todos os termos da expressão acima são não positivos exceto o termo . As regiões de derivada positiva da função W são apresentadas na figura 1. O conjunto C neste caso é constituído por três componentes conexas limitadas, são elas os conjuntos C1 e C2 e a origem, logo a condição (i) da definição 5.1 está satisfeita. Escolhendo os parâmetros de tal forma que 2k2kd > k3, a função V fica radialmente ilimitada e consequentemente a condição (ii) da definição 5.1 também é satisfeita. Portanto, V é uma função energia generalizada de controle associada ao sistema (4).


Escolhendo l = maxx C V = 0,9406 e sabendo que W é uma função radialmente ilimitada, prova-se, com o auxílio do teorema 5.2, que as soluções de (4) em malha fechada entram no conjunto limitado Sc(0,9406) para tempos positivos. O conjunto Sc(0,9406), em acordo com o teorema 5.2, é positivamente invariante e possui um conjunto H invariante não vazio e assintoticamente estável. Este conjunto H tem interseção não vazia com o conjunto C em acordo com o teorema 5.1

A seguir, vamos explorar as particularidades deste sistema para mostrar que H é a origem. Concluiu-se anteriormente que os conjuntos ω-limite do sistema (4) em malha fechada estão contidos em Sc(0,9406). De acordo com o teorema de Poincaré-Bendixson (Guckenheimer and Holmes, 1993), se o conjunto ω-limite não contém pontos de equilíbrio, ele é uma órbita fechada. Utilizando o critério de Bendixson (Guckenheimer and Holmes, 1993), verifica-se que < 0 para kd > 0, logo, não existem órbitas fechadas em Sc(0,9406). Os conjuntos limites portanto contêm o único equilíbrio do sistema em malha fechada que é a origem.

Escolhendo L = 0,04 e utilizando o teorema 4.2, prova-se que a origem é um ponto de equilíbrio assintoticamente estável, logo, como os conjuntos limites são conexos, temos que a origem é um conjunto limite isolado, ou seja, não existem órbitas homoclínicas conectando a origem a ela mesma. Portanto, a origem é um ponto de equilíbrio globalmente assintoticamente estável do sistema em malha fechada.

6.2 Exemplo 2

Considere o seguinte sistema não linear:

onde x = [x1; x2] é a variável de estado, u é a entrada de controle e

são funções de classe C1. Este sistema é uma versão simplificada de um sistema elétrico de potência constituído de um gerador e um barramento infinito com um dispositivo FACTS (Flexible AC Transmission System) para controle da impedância série da linha de transmissão. Os parâmetros, do sistema, são: a = 1, b = 2, c = 0,1 e d = 0,3. O ponto xs = [0,4728;0] e xu = [2,5689;0] são respectivamente um ponto de equilíbrio assintoticamente estável e um ponto de equilíbrio hiperbólico do tipo-1 do sistema em malha aberta (u = 0 em (4)). O objetivo é encontrar uma lei de realimentação u = h(x) tal que o equilíbrio xs seja assintoticamente estável com maior eficiência no amortecimento de oscilações e que a região de estabilidade seja suficientemente grande.

6.2.1 Solução utilizando uma FEC

Inicialmente, a lei de realimentação de controle u = h(x) será projetada utilizando uma FEC. Para isto, considere a função

como candidata a FEC do sistema (5). A constante α pode ser escolhida de maneira arbitrária, entretanto ela pode ser convenientemente escolhida de tal forma que o valor da função 1 calculada em xs seja igual a zero. Neste caso, temos α = -2,208. Considerando u = h(x1, x2), tem-se a seguinte candidata à função energia para o sistema em malha fechada:

Calculando a derivada de W1 ao longo das trajetórias do sistema (5) em malha fechada, obtém-se:

onde g(x) é dada em (7). Para que W1 seja uma função energia do sistema em malha fechada (5), deve-se escolher a lei de realimentação de controle h(x), tal que 1< 0. Com a escolha

tem-se:

onde k > 0 , é um parâmetro do controlador a ser escolhido. A escolha de k é feita conforme o nível de amortecimento desejado para o sistema. Neste exemplo, considere k = 0,5. Observe que h(xs) = h(xu) = 0, logo, a realimentação de estados h(x) não altera a posição dos equilíbrios xs e xu. Vamos mostrar que W1 é uma função energia do sistema em malha fechada. Para todo k > 0, mostrou-se que a função W1 satisfaz a condição (i) da definição 4.1. Suponha que a condição (ii) da definição 4.1 não seja atendida, então existe xoE e intervalo I = [t1,t2] ⊂ de medida não nula, com t2 > t1 tal que 1 (φ(t,xo)) = 0 para todo tI, isto é, [1 + kg2(x)] = 0 para todo tI. Como o termo entre colchetes na igualdade anterior é sempre positivo para k > 0, obrigatoriamente x2 = 0 para todo tI. Consequentemente, 1 = 0 para todo tI o que significa que no intervalo I a solução é um ponto de equilíbrio. Da unicidade de soluções, tem-se que xo é um ponto de equilíbrio, mas isto contraria o fato de que xoE. Portanto, a condição (ii) da definição 4.1 é atendida para todo k > 0. Para provar a condição (iii), observa-se que x2(t) é sempre limitado para t > 0. Isto é uma consequência direta da fórmula da variação das contantes aplicada a segunda equação diferencial do sistema (5). Então, se W1 é limitada para todo t > 0, então, necessariamente x1(t) é limitada para todo t > 0. Portanto a condição (iii) é atendida. Para estimar a região de estabilidade do ponto de equilíbrio xs, deve-se encontrar o maior número L tal que as condições do Teorema 4.2 sejam satisfeitas. Neste exemplo, a melhor estimativa da região de estabilidade na forma de um conjunto de nível da função 1 é obtida escolhendo-se, em acordo com o teorema 4.3, L = 1(xu) = 1,3740, ou seja, L é a energia do único ponto de equilíbrio na fronteira da região de estabilidade de xs. A figura 2, ilustra a estimativa da região de estabilidade. Observe que xs é o único ponto de equilíbrio contido em L = {x2: L < 1,3740}. Portanto, (8) é uma FEC do sistema (5). As trajetórias do sistema, sem e com a ação da lei de realimentação (11), podem ser observadas na figura 2. Note que (11) é eficaz no amortecimento de oscilações, sendo capaz de estabilizar o sistema em xs, para qualquer que seja xoL. Além disso, com a função energia proposta, pode-se mostrar que a estimativa da região de estabilidade do sistema em malha fechada é pelo menos tão grande quanto a estimativa da região do sistema em malha aberta, ou seja, consegue-se melhor desempenho em amortecimento das oscilações sem prejuízo a região de estabilidade do sistema.


6.2.2 Solução utilizando uma FEGC

Na subseção anterior, mostrou-se que (5) possui uma FEC. Apesar disto, uma FEGC será proposta com o objetivo de ilustrar o teorema 5.2. Para isto, considere a função candidata a FEGC:

onde z(x) = f(x) + dx2, β é um parâmetro a ser determinado e α é o valor da função 2 calculada em xs. Considerando u = h(x), a derivada de (8) ao longo das trajetórias de (5) é dada por:

onde

Observe que

2 é composta por um termo quadrático Q dependente do parâmetro β. Este parâmetro deve ser determinado de tal forma que a matriz Q seja definida positiva. Aplicando o critério de Sylvester pode-se facilmente demonstrar que isto é garantido se

Neste exemplo, considere β = -0,0339. O sinal do terceiro termo de (14), '(g(x)x2 + βz(x) g(x)) h(x)', é dependente da lei de realimentação de controle h(x). O sinal do segundo termo de (14) é indefinido. Com o objetivo de aumentar o número de termos definidos negativos em (14), escolhe-se

onde k > 0 é um parâmetro do controlador a ser determinado. Neste exemplo, considera-se k = 0,5. Desta forma, tem-se

O termo - c cos(x1)x2 será responsável por gerar regiões nas quais a derivada de 2 é positiva. O conjunto , da definição 5.1, é composto pela união de componentes conexas localizadas nas proximidades dos pontos de equilíbrio. A figura 3 mostra as duas componentes que serão de nosso interesse considerando o sistema com e sem controle. A componente 2 está próxima ao ponto de equilíbrio instável xu ao passo que a componente 1 está próxima ao ponto de equilíbrio estável xs. Claramente, a condição (i) da definição 5.1 é satisfeita. Utilizando argumentos similares aqueles utilizados na seção 6.2.1, demonstra-se que a condição (ii) da definição 5.1 também é satisfeita e portanto, a função (13) é uma FEGC do sistema (5). Além disso, pode-se mostrar que (13) é uma FEG de (5) do sistema (5) em malha aberta. O máximo valor de 2 em 1 define o conjunto l = {x2: 2 < l}. Neste exemplo, obteve-se lc = 0,0605 para o sistema controlado e lsc = 0,1286 para o sistema sem controle. Para estimar a região de estabilidade, procura-se pelo maior número L tal que as condições do teorema 5.2 sejam satisfeitas. Na prática, deve-se garantir que L = {x2: 2 < L} não intercepte a componente 2, próxima ao ponto de equilíbrio instável. Neste exemplo, obteve-se respectivamente Lc = 1,3486 e Lsc = 1,3253 para o sistema com e sem controle. Observando a figura 3 verifica-se que e que .


As trajetórias do sistema, sem e com a ação da lei de realimentação (15), podem ser observadas na figura 3. A lei de realimentação (15) é eficaz no amortecimento de oscilações, sendo capaz de estabilizar o sistema no maior conjunto invariante H contido em L, que neste caso é o ponto de equilíbrio xs. Na figura 3, comparando os conjuntos e (sistema com o controle (15)), com os conjuntos e (sistema sem controle, k = 0 em (15)), pode-se observar que a lei de realimentação de controle (15) é eficiente na redução das regiões de derivada positiva da função (13). Esta redução implica em uma estimativa um pouco menos conservadora da região de estabilidade se comparada a estimativa obtida em malha aberta.

7 CONCLUSÕES

Neste trabalho, os conceitos de função energia de controle (FEC) e função energia generalizada de controle (FEGC) foram propostos. As implicações da existência de uma função energia e uma função energia generalizada em termos dos conjuntos limites foram estudadas. Mostra-se que os conjuntos ω-limites de sistemas com funções energia são pontos de equilíbrio e que os conjuntos limites de sistemas com funções energia generalizada interceptam o conjunto de derivada não negativa. A FEC foi proposta utilizando o conceito de função energia apresentado em (Chiang et al., 1988). A principal vantagem da teoria proposta, quando comparada com a tradicional FLC, é a possibilidade de obter uma lei de realimentação de controle, para a estabilização do sistema, e também estimativas ótimas da região de estabilidade do sistema em questão. A FEGC é proposta utilizando as idéias da extensão do princípio de invariância de LaSalle apresentado em (Rodrigues et al., 2000; Rodrigues et al., 2001). A função energia generalizada de controle permite que a derivada da FEGC possua regiões limitadas, no espaço de estados, de derivada positiva. Com isto, uma classe maior de sistemas não lineares pode ser tratada via esta teoria. A FEGC permite obter uma lei de realimentação de controle para estabilização de conjuntos atrativos de sistemas não lineares, e também uma estimativa de sua região de estabilidade. Em situações particulares, pode-se provar, utilizando funções energia generalizadas de controle, estabilidade assintótica de pontos de equilíbrio em malha fechada.

AGRADECIMENTOS

Agradecemos à FAPESP e ao CNPq pelo apoio financeiro.

Artigo submetido em 14/06/2007 (Id.: 00806)

Revisado em 15/01/2008, 30/07/2008, 01/11/2008, 12/02/2009

Aceito sob recomendação do Editor Associado Prof. Takashi Yoneyama

Referências bibliográficas

  • Alberto, L. F. C. (2000). O Princípio de Invariância de LaSalle Estendido Aplicado ao Estudo de Coerência de Geradores e à Análise de Estabilidade Multi-Swing, Tese de Doutorado da Escola de Engenharia de São Carlos -Universidade de São Paulo.
  • Alberto, L. F. C. (2006). Caracterização e Estimativas da Área de Atração de Sistemas Dinâmicos Não Lineares, Tese de Livre-Docência da Escola de Engenharia de São Carlos -Universidade de São Paulo.
  • Alberto, L. F. C. and Chiang, H. D. (2008). Generalized energy functions theory and stability region characterization of complex nonlinear systems, submetido para possível publicação no IEEE Trans. on Automatic Control
  • Artstein, Z. (1983). Stabilization with relaxed controls, Nonlinear Analysis, Theory, Methods and Application 7(11): 11631173.
  • Chiang, H., Hirsch, M. and Wu, F. (1988). Stability regions of nonlinear autonomous dynamical systems, IEEE Transactions on Automatic Control 33(1): 1627.
  • Clarke, F. H., Ledyaev, Y. S., Rifford, L. and Stern, R. J. (2000). Feedback stabilization and lyapunov functions, SIAM J. Control Optim. 39(1): 2548.
  • Fallside, F. and Patel, M. R. (1965). Step-response behaviour of a speed-control system with a back-e.m.f. nonlinearity, Proc. of the IEE 112(10): 19791984.
  • Guckenheimer, J. and Holmes, P. (1993). Nonlinear Oscilations, Dynamical Systems, and Bifurcations of Vector Fields, Springer-Verlag, New York.
  • Jurdjevic, V. and Quinn, J. P. (1978). Controllability and stability, Journal of Differential Equations 28: 381389.
  • Khalil, H. K. (1996). Nonlinear Systems, Prentice-Hall New Jersey.
  • LaSalle, J. P. (1976). The Stability of Dynamical Systems, Society for Industrial and Applied Mathematics.
  • Martins, A. C. P., Alberto, L. F. C. and Bretas, N. G. (2006). Uniform estimates of attractor sets of extended lure systems using lmi's, IEEE Transactions on Automatic Control 51: 16751678.
  • Rodrigues, H. M., Alberto, L. F. C. and Bretas, N. G. (2000). On the invariance principle. generalizations and applications to synchronism, IEEE Transactions on Circuits and Systems 47(5).
  • Rodrigues, H. M., Alberto, L. F. C. and Bretas, N. G. (2001). Uniform invariance principle and synchronization. robustness with respect to parameter variation, Journal of Differential Equations 169(1): 228254.
  • Sontag, E. (1983). A lyapunov-like characterization of asymptotic controllability, SIAM J. Control Optim. 21(3): 462471.
  • Sontag, E. (1989). A universal construction of artstein's theorem on nonlinear stabilization, System & Control Letters 13: 117123.
  • Sontag, E. (1990). Mathematical control theory. Deterministic finite dimensional Systems, Springer-Verlag -New York.
  • 1
    Uma aplicação
    f:
    X
    Y é própria se para cada conjunto compacto
    D
    Y, o conjunto
    F
    -1(
    D) é compacto em
    X.
  • Datas de Publicação

    • Publicação nesta coleção
      28 Maio 2009
    • Data do Fascículo
      Jun 2009

    Histórico

    • Revisado
      15 Jan 2008
    • Recebido
      14 Jun 2007
    • Aceito
      12 Fev 2009
    location_on
    Sociedade Brasileira de Automática Secretaria da SBA, FEEC - Unicamp, BLOCO B - LE51, Av. Albert Einstein, 400, Cidade Universitária Zeferino Vaz, Distrito de Barão Geraldo, 13083-852 - Campinas - SP - Brasil, Tel.: (55 19) 3521 3824, Fax: (55 19) 3521 3866 - Campinas - SP - Brazil
    E-mail: revista_sba@fee.unicamp.br
    rss_feed Stay informed of issues for this journal through your RSS reader
    Acessibilidade / Reportar erro