Apostila Prof Trofino
Apostila Prof Trofino
APOSTILA DA DISCIPLINA
Sistemas Lineares
Conte
udo
1 Introdu
c
ao Geral
1.1
Definicoes Gerais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.2
10
1.3
10
1.4
12
1.5
Classificacao de Sistemas . . . . . . . . . . . . . . . . . . . . . . . . . . .
13
1.5.1
13
1.5.2
15
1.5.3
16
1.5.4
17
1.5.5
17
1.5.6
18
1.5.7
Sistemas inversveis . . . . . . . . . . . . . . . . . . . . . . . . . .
18
1.5.8
Sistemas estaveis . . . . . . . . . . . . . . . . . . . . . . . . . . .
19
19
1.6.1
Sinal exponencial . . . . . . . . . . . . . . . . . . . . . . . . . . .
20
1.6.2
Degrau unitario . . . . . . . . . . . . . . . . . . . . . . . . . . . .
22
1.6.3
Impulso unitario . . . . . . . . . . . . . . . . . . . . . . . . . . .
23
24
1.7.1
24
1.6
1.7
Compressao no tempo: . . . . . . . . . . . . . . . . . . . . . . . .
Sistemas Lineares
1.8
1.7.2
Expansao no tempo: . . . . . . . . . . . . . . . . . . . . . . . . .
25
1.7.3
Deslocamento no tempo: . . . . . . . . . . . . . . . . . . . . . . .
25
1.7.4
Inversao no tempo: . . . . . . . . . . . . . . . . . . . . . . . . . .
26
27
1.8.1
27
1.8.2
28
1.8.3
30
1.8.4
33
1.8.5
34
2 S
erie de Fourier
2.1
37
Introducao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
37
2.1.1
Sinais periodicos . . . . . . . . . . . . . . . . . . . . . . . . . . .
38
2.1.2
38
39
2.2.1
41
2.2.2
Observacoes Adicionais . . . . . . . . . . . . . . . . . . . . . . . .
47
2.3
48
2.4
49
2.4.1
52
O Espectro de Freq
uencias de Sinais Periodicos . . . . . . . . . . . . . .
53
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
59
2.2
2.5
3 Transformada de Fourier
61
3.1
Energia de sinais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
62
3.2
63
3.2.1
63
3.2.2
Sinal Porta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
64
Sistemas Lineares
3.2.3
Sinal Impulso . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
64
3.2.4
65
3.2.5
Sinais Senoidais . . . . . . . . . . . . . . . . . . . . . . . . . . . .
67
3.2.6
Exponencial Eterna . . . . . . . . . . . . . . . . . . . . . . . . . .
67
3.2.7
Funcoes Periodicas . . . . . . . . . . . . . . . . . . . . . . . . . .
67
Propriedades da transformada . . . . . . . . . . . . . . . . . . . . . . . .
70
3.3.1
Linearidade . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
71
3.3.2
Simetria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
71
3.3.3
Escalonamento . . . . . . . . . . . . . . . . . . . . . . . . . . . .
71
3.3.4
Deslocamento em freq
uencia e Modulacao . . . . . . . . . . . . .
72
3.3.5
Deslocamento no Tempo . . . . . . . . . . . . . . . . . . . . . . .
73
3.3.6
74
3.3.7
Diferenciacao em freq
uencia . . . . . . . . . . . . . . . . . . . . .
75
3.3.8
Convolucao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
75
3.4
Amostragem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
76
3.5
82
3.5.1
Resposta em freq
uencia de sistemas lineares . . . . . . . . . . . .
82
3.5.2
83
3.5.3
85
3.5.4
86
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
89
3.3
4 Diagramas de Bode
91
4.1
Introducao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
91
4.2
. . . . . . . .
93
4.3
Diagramas basicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
94
4.3.1
94
Constante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Sistemas Lineares
4.3.2
Polo na origem . . . . . . . . . . . . . . . . . . . . . . . . . . . .
95
4.3.3
Zero na origem . . . . . . . . . . . . . . . . . . . . . . . . . . . .
96
4.3.4
97
4.3.5
99
106
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
109
4.4
5 Transformada de Laplace
111
5.1
111
5.2
114
5.3
Propriedades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
116
5.3.1
Operacao Linear . . . . . . . . . . . . . . . . . . . . . . . . . . .
116
5.3.2
117
5.3.3
117
5.3.4
119
5.3.5
120
5.3.6
120
5.3.7
122
5.3.8
122
5.3.9
123
123
124
Transformada Inversa . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
127
5.4.1
127
5.4.2
129
5.4.3
131
132
5.4
5.5
Sistemas Lineares
5.6
134
5.7
137
5.8
138
5.9
139
143
145
147
149
150
151
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
153
6 Resposta ao Degrau
155
6.1
Introducao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
155
6.2
157
6.3
159
6.3.1
160
6.3.2
160
6.3.3
161
6.3.4
161
6.3.5
162
162
6.5
165
6.6
Problemas complementares . . . . . . . . . . . . . . . . . . . . . . . . . .
174
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
177
Introducao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
179
179
Sistemas Lineares
7.1.1
Conversao A/D . . . . . . . . . . . . . . . . . . . . . . . . . . . .
179
7.1.2
180
7.2
185
7.3
Transformada Z . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
188
7.3.1
Definicao e exemplos . . . . . . . . . . . . . . . . . . . . . . . . .
189
7.3.2
191
Propriedades da Transformada Z . . . . . . . . . . . . . . . . . . . . . .
192
7.4.1
Linearidade . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
192
7.4.2
196
7.4.3
196
7.4.4
197
7.4.5
Convolucao Discreta . . . . . . . . . . . . . . . . . . . . . . . . .
198
Transformada Z Inversa . . . . . . . . . . . . . . . . . . . . . . . . . . .
199
7.5.1
200
7.5.2
201
7.6
201
7.7
205
7.7.1
. . . . . . . . . . . . .
205
7.7.2
207
7.8
Sistemas Amostrados . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
208
7.9
Sistemas Realimentados . . . . . . . . . . . . . . . . . . . . . . . . . . .
216
220
220
222
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
225
7.4
7.5
Captulo 1
Introduc
ao Geral
1.1
Definic
oes Gerais
Para se iniciar o estudo de sistemas lineares, primeiramente deve-se definir o que sao
sinais e o que sao sistemas.
uma funcao ou seq
Sinal: E
uencia de valores que representa informacao. Eles descrevem grandezas cujos valores se alteram em funcao de uma variavel independente,
geralmente o tempo.
um conjunto de dispositivos interconectados que realizam uma tarefa esSistema: E
pecfica, estabelecendo uma relacao entre um certo n
umero de sinais. Um sistema
pode representar um equipamento ou fenomeno fsico.
Normalmente, os sinais relacionados ao sistema podem ser classificados como sinais de
entrada e sinais de sada. Os sinais de entrada existem independentemente do sistema
e nao sao afetados por ele. Ja os sinais de sada representam informacoes geradas pelo
sistema, geralmente em resposta aos sinais de entrada. Os sistemas sao geralmente representados atraves de diagramas de blocos. A forma mais basica de representacao de
um sistema e apresentada na Figura 1.1
Sistemas Lineares
10
Processo: Fenomenos (naturais ou criados artificialmente) que evoluem progressivamente segundo dinamicas que lhe sao proprias. (Fenomeno a ser controlado: processos qumicos, economicos, biologicos,...).
Perturbac
ao ou Dist
urbio: Sinal indesejado (interno ou externo).
Controle Realimentado: Operacao que visa corrigir (automaticamente ou manualmente) certas variaveis (grandezas fsicas) de um sistema. Diminui o efeito de
fenomenos indesejaveis.
um sistema realimentado para controle automatico de posicao,
Servomecanismo: E
velocidade ou aceleracao. Muito frequente na ind
ustria.
Sistemas Reguladores Autom
aticos: Sistema de controle cujo principal objetivo e
manter constante algumas variaveis do mesmo. (Controle de nvel constante, posicao constante, velocidade, aceleracao, ...). Exemplos: robos, elevadores, estufas,...
A teoria de sistemas e uma ferramenta matematica para seu estudo e projeto. Inicialmente, a teoria e apresentada de uma forma generica, abstrata, sem vnculo com alguma
aplicacao especfica. Desta forma ela pode ser utilizada em diferentes areas da engenharia
eletrica (bem como de outras engenharias, fsica, economia, etc.).
1.2
Sistemas onde a variavel a ser controlada (sada) nao interfere na acao de controle
(variavel de entrada) sao conhecidos como Sistemas de malha aberta.
Perturbac
oes
Entrada
SISTEMA
Sada
1.3
Sistemas Lineares
11
perturbac
ao
Vari
avel
Observada
Ref.
Comparador
Controlador
Atuador
SISTEMA
Medidor
sinal de medic
ao
rudo de medic
ao
Controlador
Ref.
Comparador
A/D
D/A
Computador
Atuador
Sada
SISTEMA
Medidor
posica
o
da antena
c(t)
potenci
ometro
refer
encia
r(t)
comparador
Vr (t)
Vc (t)
potenci
ometro
erro e(t)
engrenagem
amplificador
de pot
encia
Ea (t)
motor DC
Sistemas Lineares
12
1.4
TEMPO DISCRETO: t e uma variavel discreta que assume valores apenas em instantes
discretos do tempo. Por exemplo, t = kT onde k e uma variavel k = 0, 1, 2, . . . e T e uma
constante. Nesse caso um sinal f (kT ) sera uma seq
uencia de pulsos, isto e, um sinal de
tempo discreto, Figura 1.7.
Sistemas Lineares
13
f(kT)
Ref.
t = kT
1.5
Classificac
ao de Sistemas
Os sistemas podem ser classificados de acordo com diversas categorias. Esta classificacao e importante na escolha das tecnicas de analise e projeto de sistemas.
1.5.1
Sistemas lineares e n
ao-lineares
A propriedade da linearidade e uma das mais importantes que serao consideradas neste
curso. As tecnicas de solucao de sistemas lineares estao bem definidas e consolidadas.
Estas tecnicas tambem sao geralmente utilizadas em sistemas nao-lineares fazendo-se
aproximacoes em torno de um determinado ponto de operacao.
Um sistema e dito linear se obedecer as regras abaixo, considerando x(t) como um sinal
de entrada, y(t) o sinal de sada e representando a transformacao que o sistema aplica
ao sinal de entrada.
1. Princpio da homegeneidade: Se uma entrada e multiplicada por um fator k, a sada
tambem e multiplicada pelo mesmo fator k:
k.x(t) k.y(t)
2. Princpio da aditividade: Se duas entradas sao aplicadas simultaneamente, a resposta total sera a soma das respostas individuais a cada uma das entradas separadamente. Assim, se x1 (t) y1 (t) e x2 (t) y2 (t), entao:
x1 (t) + x2 (t) y1 (t) + y2 (t)
Estes dois princpios, se combinados, resultam no princpio da superposicao. Um sistema e classificado como linear se e somente se atender a este princpio, ou seja, se a
Sistemas Lineares
14
Nota-se que a linearidade nao necessariamente implica que a funcao de tempo da sada
tem a mesma forma da entrada, por exemplo:
Quando se trata de sistemas reais, muitas vezes a linearidade se aplica somente para
uma faixa restrita de valores das constantes k1 , k2 e das entradas. Por exemplo, a curva
de deformacao de um certo material em funcao da tensao mecanica aplicada sobre o
corpo de prova. Neste caso, a medida em que a tensao aplicada aumenta, a deformacao
aumenta proporcionalmente em uma certa faixa de valores de entrada. Para esta faixa de
valores, o sistema se comporta como linear. Se a tensao continuar aumentando o sistema
entrara em um regiao de nao-linearidade, podendo ate haver o rompimento do material.
Um grafico exemplificando este comportamento pode ser visto na Figura 1.10.
Sistemas Lineares
15
1.5.2
A variancia/invariancia no tempo e a segunda caracterstica mais importante dos sistemas a serem estudados nesta disciplina.
Os sistemas invariantes no tempo podem ser definidos da seguinte forma:
Sistema Invariante no Tempo: Um sistema que responde a um sinal de entrada deslocado no tempo com um sinal de sada deslocado no tempo correspondente e chamado
de sistema invariante no tempo.
Descrevendo de outra forma, um sistema variante e aquele onde seus parametros variam com o tempo, como por exemplo, um circuito RLC onde se considera o efeito da
temperatura. Ja o sistema invariante pode ser tambem descrito como aquele em que os
parametros nao variam com o tempo. Exemplo, um circuito RLC onde nao se considera
o efeito da temperatura.
A Figura 1.11 apresenta dois sistemas S1 e S2. Nestes dois sistemas e aplicada uma
entrada x(t) e a mesma entrada deslocada de t1 unidades de tempo. Pela figura pode-se
perceber que o sistema S1 e invariante no tempo, pois a sada deslocada s1 (t) tem a
mesma forma da sada nao deslocada s1(t). Ja para o sistema S2, a mesma afirmativa
nao e verdadeira, pois a sada deslocada tem uma forma diferente da original.
Os sistemas invariantes no tempo podem ser modelados por equacoes diferenciais de
coeficientes constantes. Por exemplo, um circuito RLC sem considerar as variacoes dos
parametros com o tempo:
V (t) = R.i(t) + L.
di(t)
1
+
dt
C
i(t).dt
V (t) = R(t).i(t) + L.
1
di(t)
+
dt
C
i(t).dt
Sistemas Lineares
16
1.5.3
Sistemas est
aticos e din
amicos (sem mem
oria e com mem
oria)
x( )d
y(t) = K
onde a tensao de sada y(t) depende de todos os valores anteriores da entrada de tensao
x(t).
Sistemas Lineares
17
1.5.4
Sistemas causais e n
ao causais
1.5.5
Sistemas Determinsticos e N
ao Determinsticos
Sistemas Lineares
1.5.6
18
Sistemas a par
ametros concentrados ou distribudos
Um sistema e a par
ametros concentrados quando seus elementos sao discretos, nao
se considerando o seu tamanho fsico. Por exemplo, um resistor: R = 2k. Estes sistemas
sao modelados por equacoes diferenciais ordinarias.
Quando um sistema e a par
ametros distribudos e considerada a distribuicao dos
seus parametros do espaco. Sua modelagem e dada por equacoes diferenciais parciais
onde esta presente a variavel espaco.
Um exemplo de sistema a parametros distribudos sao as linhas de transmissao, Figura
1.13, onde a linha e separada em diversos segmentos, cada qual com seu modelo equivalente. Faz se entao o tamanho de cada segmento tender a zero e com isso chega-se a uma
equacao diferencial parcial que descreve a linha em funcao da distancia.
OBS.: Nos sistemas a parametros concentrados, a dimensao do elemento e muito pequena comparada com o comprimento de onda dos sinais a serem transmitidos. A modelagem de um sistema por parametros concentrados ou distribudos deve ser escolhida de
acordo com a aplicacao. Muitas vezes, uma modelagem detalhada nao se faz necessario
para o tipo de estudo que se deseja realizar.
1.5.7
Sistemas inversveis
1
w(t) = .y(t)
3
Sistemas Lineares
19
Seja o sistema y(t) = x2 (t). Neste caso, diferentes entradas resultam na mesma sada.
Por exemplo, y(t) = 4 e obtido se x(t) = 2 ou x(t) = 2. Portanto, este sistema nao e
inversvel.
Para um sistema ser inversvel, diferentes entradas devem levar a diferentes sadas.
1.5.8
Sistemas est
aveis
Os sistemas tambem podem ser estaveis ou nao. Existem diversas maneiras de se definir
a estabilidade de um sistema. Aqui sera apresentada a estabilidade BIBO (Bounded Input
Bounded Output) que significa Entrada Limitada Sada Limitada.
Estabilidade BIBO: Um sistema e dito BIBO estavel se para qualquer entrada limitada
a sada tambem sera limitada.
A entrada do sistema e limitada se existir um n
umero M tal que
|x(t)| M
para todo t. Da mesma forma, a sada e limita se existir um n
umero R tal que a desigualdade
|y(t)| R
seja valida para todo t.
y(t) = K.
x( ).d
0
Se um sinal limitado for aplicado `a sua entrada, por exemplo um degrau unitario, a
sada sera uma rampa (a integral de uma constante e a propria variavel de integracao),
y(t) = K.t, a qual e ilimitada pois para t a sada y(t) . Logo, o integrador
nao e BIBO estavel.
1.6
Sinais Contnuos B
asicos
Nesta secao serao apresentados os os modelos de sinais contnuos mais comuns para o
estudo de sistemas.
Sistemas Lineares
1.6.1
20
Sinal exponencial
di(t)
R
di(t)
+ R.i(t) = 0
= .i(t)
dt
dt
L
Para este exemplo, s = R/L. A corrente e dada entao pela solucao da equacao diferencial acima, ou seja:
R
i(t) = i(0).e L .t
Este tipo de sinal exponencial e muito comum em engenharia eletrica e deve ser dada
muita importancia no seu correto entendimento.
No exemplo anterior, a constante s resultante era um n
umero real. No entanto, podem
haver situacoes onde, para facilitar a solucao dos sistemas, e mais interessante que s seja
um n
umero complexo do tipo s = + j, onde denota a parte real do n
umero e a
parte imaginaria. Entao, o sinal exponencial fica:
es0 .t = e(0 +j0 ).t = e0 .t .ej0 .t
Este sinal exponencial sera apresentado agora de acordo com a caracterstica de s0 , ou
seja, s0 sendo um n
umero real puro, imaginario puro e combinado.
s real puro:
Neste caso, x(t) = e0 .t . Os graficos para valores positivos e negativos e 0 sao vistos
na Figura 1.15.
s imagin
ario puro:
Neste caso, x(t) = ej0 .t . Utilizando-se a relacao de Euler,
ej0 .t = cos 0 .t + j. sin 0 .t
pode-se verificar que o sinal em questao e periodico, pois:
Sistemas Lineares
21
O perodo fundamental e
T0 =
2.
0
Este sinal representa um fasor de amplitude unitaria que gira no plano complexo,
representando um sinal senoidal, conforme a Figura 1.16.
s real e imagin
ario:
Quando s possui parte real e imaginaria nao nulas, o sinal torna-se: x(t) = e(0 +j0 ).t
Este sinal comporta-se como o sinal anterior, ou seja, representa um fasor girando
no plano complexo. No entanto, neste caso, o modulo nao sera constante. Havera um
amortecimento no modulo, dado por 0 . O amortecimento so e verificado se 0 < 0,
caso contrario havera um incremento no modulo. Na figura 1.17 pode-se ver uma representacao deste sinal (para o caso amortecido) em um grafico tridimensional onde, alem
das componentes real e imaginaria, foi representada a variavel tempo t.
Sistemas Lineares
22
Imag.
0.5
0.5
1
0.5
1
0
0
0.2
0.4
Tempo
0.6
0.5
0.8
Real
1.6.2
Degrau unit
ario
O degrau unitario tambem e um sinal muito utilizado na analise de sistemas. Normalmente, o degrau funciona como sinal de excitacao dos sistemas fsicos. Nos sistemas
realimentados, e comum a aplicacao de um degrau na entrada de referencia do sistema.
A notacao utilizada para representar um degrau unitario e u(t), que e definido da
seguinte forma:
0, se t < 0
1, se t > 0
u(t) =
u (t) =
0,
t
,
1,
se t < 0
se 0 t
se t >
Sistemas Lineares
23
1.6.3
Impulso unit
ario
(t) =
d
u (t)
dt
0,
se t < 0
= 1 , se 0 t
0, se t >
Observa-se que (t) tem area unitaria. Ao se aplicar o limite nesta funcao, chega-se
a definicao do impulso unitario:
du (t)
dt
0
(t) = dtd u(t)
d
lim u (t)
dt 0
(1.1)
(1.2)
Figura 1.18: Representacao do sinal impulso. (a) Unitario; (b) Multiplicado por k; (c)
Multiplicado e deslocado no tempo.
Propriedades da func
ao impulso
A seguir sao apresentadas algumas propriedades da funcao impulso.
1) x(t).(t) = x(0).(t) = Se um sinal x(t) e multiplicado pela funcao impulso,
entao o resultado sera uma amostra do sinal x(t) no instante zero, ja que o impulso
e definido neste instante.
Sistemas Lineares
24
Interpretacao:
t<0
4)
k.( ).d = k
(1.3)
1.7
f (t).(t t0 ).dt =
(t).dt = f (t0 )
Transformac
oes de Vari
avel Independente
0,
1,
f (t) =
1 0, 5.t,
0,
1.7.1
se
se
se
se
t<1
1 t < 0
0t<2
t2
Compress
ao no tempo:
Sistemas Lineares
25
0,
se
2t
<
1
0,
1,
se 1 2t < 0
1,
f (2t) =
f (2t) =
1 0, 5.(2t), se 0 2t < 2
1 t,
0,
se 2t 2
0,
se
se
se
se
t < 1/2
1/2 t < 0
0t<1
t1
1.7.2
Expans
ao no tempo:
0,
1,
f (1/3t) =
1 1/6t,
0,
t < 3
3 t < 0
0t<6
t6
1.7.3
Deslocamento no tempo:
O sinal e deslocado para a direita se um valor negativo e adicionado `a variavel independente. Reciprocamente, o sinal e deslocado para a esquerda se um valor positivo e
adicionado.
Sistemas Lineares
26
Exemplo, f (t + 4):
Susbstituindo t por t + 4 na definicao de f (t), chega-se ao seguinte sinal modificado:
0,
1,
f (t + 4) =
1 1/2t,
0,
se
se
se
se
t < 5
5 t < 4
4 t < 2
t 2
1.7.4
Invers
ao no tempo:
0,
1 + 1/2t,
f (t) =
1,
0,
se
se
se
se
t 2
2 < t 0
0<t1
t>1
Sistemas Lineares
27
Exerccios:
Calcule analiticamente e desenhe os seguintes sinais:
a) f (t 3);
b) f (1/2t);
c) f (3 t);
d) f (2t 3)
1.8
Soluc
ao de Sistemas Lineares Via Convoluc
ao
Para determinar a sada de um sistema linear para uma dada entrada x(t), pode-se utilizar duas tecnicas basicas. A solucao via equacoes diferenciais e a solucao via convolucao
da entrada com a resposta impulsiva do sistema. Esta segunda tecnica sera descrita a
seguir, porem, antes de apresenta-la, deve-se responder algumas questoes importantes,
como, por exemplo, o que e resposta impulsiva e o que e convolucao.
1.8.1
Representac
ao de sinais contnuos em termos de impulsos
(t) =
0,
1
,
0,
se t < 0
se 0 t
se t >
Sistemas Lineares
28
Entao, uma funcao contnua x(t) pode ser aproximada por uma somatorio de impulsos
deslocados e ponderados. Denotando x(t) como o sinal aproximado:
x(t) =
k=
x(k).. (t k)
(1.4)
Figura 1.23: Aproximacao de uma funcao contnua x(t) por impulsos de largura .
Tomando o limite para tendendo a zero, pode-se encontrar x(t) a partir de sua
aproximacao x(t).
Relembrando, quando 0:
ser escrito como:
x(t) =
x( ).(t ).d
(1.5)
A verificacao da equacao anterior se da atraves da utilizacao da propriedade de multiplicacao de um sinal pelo impulso (segunda propriedade do impulso unitario, da secao
1.6.3) e sabendo que a integral do impulso unitario e 1. Entao:
1.8.2
x( ).(t ).d =
(t ).d = x(t)
Nesta secao sera apresentada uma maneira de se calcular a sada y(t) de um sistema
linear e invariante no tempo (LIT) para uma dada entrada x(t) a partir do conhecimento
da resposta impulsiva do sistema. Para isto, serao utilizadas as propriedades da linearidade e invariancia no tempo. Portanto, todo o desenvolvimento apresentado sera valido
apenas para sistemas lineares e invariantes no tempo (LIT).
Sistemas Lineares
29
Resposta impulsiva
O que e a resposta impulsiva de um sistema LIT?
Considere a Figura 1.24 abaixo. Dado um sistema LIT onde a entrada aplicada e um
impulso unitario x(t) = (t), a sada para esta entrada sera definida como a resposta
impulsiva do sistema y(t) = h(t). A notacao h(.) sempre representa uma resposta ao
impulso unitario.
x( ).(t ).d
x( ).h(t ).d
Note que a representacao em blocos do sistema da Figura 1.24 fui substituda por uma seta apontando
para a direita. Esta notacao sera utilizada de agora em diante para representar o sistema.
Sistemas Lineares
30
x(t)
y(t) =
x( ).h(t ).d
(1.6)
y(t) =
x( ).h(t ).d =
(1.7)
Este resultado e de fundamental importancia no estudo de sistemas lineares e invariantes no tempo. A equacao (1.7) estabelece que a resposta do sistema para qualquer
entrada x(t) e expressa como uma integral envolvendo apenas a entrada e sua resposta
comum expressar que a sada e igual `
impulsiva h(t). E
a entrada convoluda com a resposta impulsiva. Tambem e comum denotar a integral de convolucao de uma maneira
mais simplificada:
y(t) = x(t) h(t) = h(t) x(t)
Conclus
ao
O calculo da resposta de estado zero (parcela da resposta total do sistema sem considerar as condicoes iniciais) de qualquer sistema linear e invariante no tempo a uma entrada
qualquer reduz-se a:
1) Determinacao da resposta impulsiva do sistema h(t);
2) Calculo da integral de convolucao:
y(t) =
x( ).h(t ).d =
t0
1.8.3
Representac
ao gr
afica da integral de convoluc
ao
A integral de convolucao da equacao (1.7) pode ser resolvida graficamente e por partes, conforme a forma do sinal de entrada e da resposta impulsiva do sistema. Este
procedimento sera apresentado atraves de um exemplo.
Sistemas Lineares
31
Considere um sistema com a resposta impulsiva h(t) que foi submetido a uma entrada
x(t). Ambos os sinais sao apresentados na Figura 1.25.
Sistemas Lineares
32
y0<t<2 (t) =
0
(2 ).d =
2.
2
2
t
0
= 2.t
t2
2
(1.8)
Para t 2:
Quando t 2, a area em comum e sempre a mesma conforme t e aumentado. Esta
situacao e mostrada na Figura 1.28.
yt2 (t) =
0
(2 ).d =
2.
2
2
2
0
= 2.2
4
=2
2
(1.9)
Sistemas Lineares
33
1.8.4
Propriedades da operac
ao de convoluc
ao
Sistemas Lineares
1.8.5
34
Sera derivada agora uma expressao geral para a sada de um sistema com resposta
impulsiva h(t) quando uma entrada degrau unitario e aplicada em sua entrada. Denotase s(t) como a sada do sistema quando x(t) = u(t).
Adequando o sinal de entrada assim como foi feito no exemplo da secao 1.8.3, chega-se
a seguinte expressao para a integral de convolucao:
s(t) =
h( ).u(t ).d
Atraves Figura 1.31 que mostra a entrada deslocada u(t ) utilizada na convolucao,
percebe-se que a sada de um sistema para esta entrada pode ser descrita como a integral
de ate t da propria resposta impulsiva do sistema. Ou seja:
t
1.h( ).d
s(t) =
Bibliografia
PHILLIPS, C. L.Signals, Systems, and Transforms, Prentice-Hall Inc.:New Jersey, USA,
1995, 707p.
GIROD, B.; RABENSTEIN, R.; STENGER, A. Sinais e Sistemas, LTC Editora:Rio de
Janeiro, RJ, 2003, 340p.
Sistemas Lineares
36
Captulo 2
S
erie de Fourier
2.1
Introduc
ao
t2
+ ...
2!
Sistemas Lineares
38
linear para uma entrada complexa somando-se as sadas para entradas mais simples, o
que e computacionalmente mais atraente.
Jean Baptiste Joseph Fourier (1768-1830) foi um matematico e fsico frances que estudou extensivamente a conducao do calor. Ele desenvolveu o que agora e chamado de
An
alise de Series de Fourier, para aplicar `a solucao de equacoes diferenciais parciais que
surgiram nos seu estudos de conducao de calor. A princpio ele nao foi capaz de provar
de uma forma geral que sua infinita serie de senos e cossenos realmente convergia para a
funcao que deveria representar. Esta prova so foi obtida aproximadamente 100 anos mais
tarde.
Antes de se apresentar a serie de Fourier propriamente dita, deve-se definir o que sao
funcoes periodicas e a resposta de sistemas LIT a funcoes exponenciais, uma vez que a
grande maioria dos sinais utilizados neste captulo se enquadram nestas categorias.
2.1.1
Sinais peri
odicos
Uma funcao x(t) e dita periodica se para um perodo T a relecao x(t) = x(t + T ) e
satisfeita para todo t. Por exemplo, a funcao cos .t e periodica ( = 2f = 2/T ) pois:
cos (t + T ) = cos(t + T ) = cos(t + 2) = cos t
As funcoes periodicas tem as seguintes propriedades:
1) Funcoes periodicas sao validas para todo tempo, ou seja, na equacao x(t) = x(t+T )
nao ha limite imposto em t;
2) A funcao periodica com perodo T e tambem periodica com perodo nT , onde n e
um inteiro. Entao, para uma funcao periodica:
x(t) = x(t + T ) = x(t + nT )
3) O perodo fundamental T0 e o mnimo valor do perodo T > 0 que satisfaz x(t) =
x(t + T ). A freq
uencia fundamental e definida como 0 = 2f0 = 2/T0 .
2.1.2
h( ).x(t ).d =
h( ).es(t ) .d
Sistemas Lineares
39
st
y(t) = e
h( ).es .d = H(s).est
2.2
As S
eries Trigonom
etricas de Fourier
(2.1)
onde T e o perodo.
O teorema de Fourier mostra que uma funcao periodica f (t) pode ser representada por
um somatorio ponderado de senos e cossenos. Para isso, f (t) deve satisfazer algumas
propriedades, que sao condicoes suficientes para atestar a existencia de uma serie de
Fourier. Estas condicoes, tambem chamadas de Condicoes de Dirichlet, sao descritas a
seguir. Qualquer sinal que apareca em um sistema fsico ira satisfazer estas condicoes.
1) Existe um u
nico valor de f (t) para qualquer t;
2) A integral
t0 +T
t0
|f (t)| .dt
(2.2)
Sistemas Lineares
40
f (t) = a0 +
(2.3)
k=1
A freq
uencia fundamental 0 e relacionada ao perodo T atraves de
0 =
2.
T
ak =
bk =
t0 +T
1
T
a0 =
f (t).dt
(2.4)
(2.5)
(2.6)
t0
t0 +T
2
T
2
T
t0
t0 +T
t0
Percebe-se que o termo a0 representa o valor medio da funcao f (t) no perodo, ou seja,
o termo com freq
uencia nula.
A Figura 2.1 apresenta um exemplo de um sinal periodico quadrado decomposto em
serie de Fourier. Sao apresentadas quatro curvas, representando o sinal formado por uma
serie de uma, tres, cinco e infinitas componentes.
Exemplo 2.1 Determine a serie de Fourier do sinal da Figura 2.2.
O perodo deste sinal e T = 2.. A freq
uencia angular e 0 = 2./T = 1 rad/s
C
alculo dos termos da serie de Fourier:
a0 =
ak =
bk =
2
2
2
2
1
2
f (t).dt = 0
0
1. cos kt.dt +
0
1. cos kt.dt = 0
1. sin kt.dt +
0
1. sin kt.dt =
4
,
k
0,
se k e mpar
se k e par
Sistemas Lineares
41
f (t)
1.5
k=1
k=3
k=5
1.1
0.7
0.3
-0.1
-0.5
-0.9
-1.3
0
A express
ao para a serie de Fourier deste sinal e ent
ao dada pela equacao abaixo, onde
o termo (2k 1) e utilizado para anular a express
ao quando k e par.
4
sin(2k 1)t
f (t) = .
k=1 (2k 1)
Expandindo em serie:
f (t) =
2.2.1
4
4
4
. sin t +
. sin 3t +
. sin 5t + . . .
3
5
Sistemas Lineares
42
classificacao, o calculo da serie de Fourier para cada uma pode ser simplificado em funcao
de sua simetria.
Fun
c
oes Pares: Funcao par e aquela em que f (t) = f (t). Esta funcao e simetrica em
relacao ao eixo vertical. Considere o sinal exemplo da figura abaixo:
f (t).dt
f (t).dt e
0
f (t).dt
f (t).dt = 2
0
f (t).dt = 0
a
Sistemas Lineares
43
1
[x(t) + x(t)]
2
xi (t) =
1
[x(t) x(t)]
2
Conseq
uentemente:
A partir destas propriedades pode-se perceber que o calculo dos coeficientes da serie
de Fourier para sinais com simetria pode ser simplificado fazendo as integrais com limites
de T /2 a T /2 ao inves de 0 a T . Desta forma:
a0 =
ak =
2
T
bk =
2
T
T /2
T /2
T /2
T /2
1
T
T /2
T /2
f (t).dt
(2.7)
(2.8)
(2.9)
Utilizando estes novos limites de integracao pode-se verificar que dependendo da funcao
f (t) ser par ou mpar alguns termos da serie de Fourier nao precisam ser calculados, pois
ja sabe-se de antemao que sao nulos.
Sistemas Lineares
44
Considere que f (t) e par. Com isso, f (t). cos k.0 t e tambem par, pois conforme visto
anteriormente, o produto de duas funcoes pares e tambem uma funcao par. Da mesma
forma, o produto f (t). sin k.0 t e mpar. Entao as equacoes (2.7), (2.8) e (2.9) tornam-se:
a0 = 2.
ak = 2.
2
T
1
T
T /2
f (t).dt
(2.10)
T /2
para k = 0
(2.11)
bk = 0
para k = 0
(2.12)
a0 = ak = 0
bk = 2.
2
T
(2.13)
T /2
para k = 0
(2.14)
Estes resultados sao coerentes com o raciocnio intuitivo: Se senos sao funcoes mpares,
estes nao devem estar presentes na representacao em serie de Fourier de funcoes pares.
Caso isto acontecesse, a serie de Fourier seria um somatorio de funcoes pares e mpares,
o que resultaria em uma funcao nem par nem mpar. Por isso, necessariamente os termos
bk tem que ser nulos para a representacao em serie de Fourier de funcoes pares. O mesmo
raciocnio e valido para representacao de funcoes mpares, onde os termos pares (ak ) e o
termo constante sao nulos. Alem disso, quando se representa funcoes simetricas atraves
de series de Fourier, os termos nao nulos podem ser obtidos atraves da integracao sobre
meio perodo do sinal. Multiplicando-se as integrais por 2 para corresponderem ao perodo
total T .
Em alguns casos e interessante deslocar provisoriamente o eixo vertical ou horizontal
da uma funcao, ou mesmo ambos, de modo a criar uma funcao par ou mpar e com isso
tirar proveito das simplificacoes possveis para funcoes simetricas. A serie de Fourier
encontrada para esta funcao modificada e ao final deslocada novamente para se obter a
representacao correta da funcao original. O exemplo a seguir ilustrara este procedimento.
Sistemas Lineares
45
Exemplo 2.2 Encotre as series de Fourier para f1 (t) e f2 (t), dadas pelas Figuras 2.5 e
2.6. Utilize as simplificacoes para funcoes simetricas.
1
2
2
2
=
= 1 rad/s
T
2
J
a que f1 (t) e uma funcao par:
bk = 0
e
ak = 2.
t
. cos kt.dt
ak =
2
.
1
cos k
2
2
k
k
4
,
2 k2
0,
se k e mpar
se k e par
f1 (t) =
1
4
1
1
+ 2 . cos t + . cos 3t + . cos 5t + . . .
2
9
15
Sistemas Lineares
46
Caso n
ao fossem utilizadas as simplificacoes da simetria, as seguintes express
oes deveriam ser calculadas:
ak =
bk =
t
. cos kt.dt +
t
. sin kt.dt +
cos kt.dt
sin kt.dt
bk = 2.
1
1
. sin kt.dt =
. (cos k. 1) =
2
k.
2
,
k
0,
se k e mpar
se k e par
f3 (t) =
1
1
2
. sin t + . sin 3t + . sin 5t + . . .
3
5
f2 (t) =
1
1
1 2
+ . sin t + . sin 3t + . sin 5t + . . .
2
3
5
Alternativamente, o eixo vertical poderia ter sido deslocado para criar uma funcao par,
f4 (t), vista na Figura 2.8.
Sistemas Lineares
47
1
2
Os termos ak s
ao dados pela seguinte integral:
ak = 2.
/2
1. cos kt.dt =
0
ak =
2
,
k
2
,
k
se k e mpar
se k e par
(k1)
2
.(1) 2
k
Expandindo em serie:
f4 (t) =
1
1
1 2
+ . cos t . cos 3t + . cos 5t + . . .
2
3
5
1 2
3
1
+ . cos t
. cos 3t
2
2
3
2
+ ...
Utilizando a relacao trigonometrica cos(x + y) = cos x. cos y sin x. sin y obtem-se uma
express
ao final simplificada para f2 (t):
f2 (t) =
2.2.2
1
1
1 2
+ . sin t + . sin 3t + . sin 5t + . . .
2
3
5
Observac
oes Adicionais
k=1
Sistemas Lineares
48
tanto os termos seno quanto cosseno estao presentes, algumas vezes e preferida a utilizacao
do tipo:
f (t) = a0 +
Ak . cos(k.0 t + k )
k=1
onde:
Ak =
a2k + b2k
bk
ak
pois desta forma se mostra claramente a amplitude Ak e o angulo de fase da kesima
harmonica.
k = tan1
2.3
Propriedades das S
eries de Fourier
+
x(t
a ) + x(ta )
a0 +
(ak . cos(k.0 t) + bk . sin(k.0 t)) =
2
k=1
onde x(t
e o valor limite de x(t) quando t se aproxima de ta pela esquerda e x(t+
a)
a)
e o valor limite pela direita. Percebe-se que a equacao acima tambem e satisfeita
para ta como ponto de continuidade;
3) Quase toda funcao contnua x(t) de perodo T0 pode ser uniformemente aproximada
por uma serie truncada de Fourier com qualquer grau pre-definido de exatidao. O
erro desta aproximacao e dado por:
N
e(t) = x(t) a0
Esta propriedade diz que este erro pode ser limitado por qualquer valor nao nulo
com a escolha de N suficientemente grande;
4) O erro medio quadratico (EM Q) e definido como:
EM Q =
1
T0
t0 +T
e2 (t).dt
t0
Sistemas Lineares
49
k=
1
k=
k=3
k=5
7
A terceira propriedade e bastante importante pois diz que quase qualquer funcao
contnua pode ser aproximada por uma serie de Fourier truncada, com qualquer grau
de exatidao. Considere agora a aproximacao de uma funcao nao contnua, como um sinal
periodico quadrado. A Figura 2.9 ilustra esta aproximacao. Em (a), esta o primeiro
harmonico, em (b) a soma do primeiro e terceiro harmonico e em (c) o somatorio ate o
nono harmonico. Conforme e aumentado o n
umero de componentes na serie de Fourier,
o erro de aproximacao e evidentemente diminudo.
3
4
2
4
2
4
3
4
2.4
S
erie Exponencial de Fourier
Nesta secao, sera derivada uma outra forma da serie de Fourier, que e reescrita nao
mais como um somatorio de senos e cossenos, mas sim como um somatorio de termos
Sistemas Lineares
50
f (t) = a0 +
ak . cos(k.0 t) +
k=1
bk . sin(k.0 t)
(2.15)
k=1
Para determinar a forma exponencial (ou complexa) da serie de Fourier sao utilizadas
as seguintes relacoes:
ej + ej
2
j
e ej
sin =
j2
j
e
= cos j sin
cos =
(2.16)
(2.17)
(2.18)
f (t) = a0 +
k=1
ak
bk
+
2
j2
jk0 t
.e
k=1
ak
bk
2
j2
.ejk0 t
f (t) = a0 +
k=1
1
1
(ak jbk ) .ejk0 t +
(ak + jbk ) .ejk0 t
2
2
k=1
(2.19)
Fazendo
1
Fk = .(ak jbk )
2
1
Fk = .(ak + jbk )
2
(2.20)
(2.21)
e verificando que Fk = Fk
, pode-se reescrever a equacao (2.19) como:
f (t) =
Fk .ejk0 t
(2.22)
k=
Para se chegar a expressao geral dos coeficientes da serie exponencial de Fourier, usanse as equacoes dos termos ak e bk , dados pelas equacoes (2.5) e (2.6), substituindo-as na
Sistemas Lineares
51
1
T
t0 +T
t0
entao:
Fk =
t0 +T
1
T
f (t).ejk0 t .dt
(2.23)
t0
ak jbk
=
2
a2k + b2k
tan1
2
Fk =
bk
ak
Ak
k
2
Fk = |Fk | k ; Fk = |Fk | k
1
|Fk | = |Fk | = . a2k + b2k
2
bk
k = arctan
ak
k = k
(2.24)
(2.25)
(2.26)
(2.27)
A serie exponencial de Fourier nao e uma serie diferente da trigonometrica, mas sim,
ambas sao formas diferentes de se expressar a mesma serie. Com isso, os coeficientes de
uma podem ser obtidos da outra. Em parte, isto pode ser ser verificado nas equacoes
(2.20) e (2.21) onde os coeficientes da serie exponencial sao determinados atraves dos
coeficientes da serie trigonometrica. Ja a relacao inversa, ou seja, a determinacao dos
coeficientes da serie trigonometrica em funcao daqueles da serie exponencial e dada pelas
equacoes a seguir.
a0 = F0
ak = Fk + Fk
bk = j(Fk Fk )
(2.28)
(2.29)
(2.30)
Sistemas Lineares
2.4.1
52
S
erie de Fourier exponencial para sinais sim
etricos
2
T
2
=
T
T /2
F0 =
Fk
f (t).dt
(2.31)
(2.32)
(2.33)
0
T /2
0
2
Fk = j
T
T /2
0
A mudanca de eixos feita para tornar o sinal par ou mpar tambem pode ser feita para
a serie exponencial, da mesma forma do que na serie trigonometrica.
Problema 2.1 Deduza as equacoes expostas acima dos coeficientes da serie exponencial
de Fourier para sinais simetricos.
Exemplo 2.3 Encontre a serie de Fourier exponencial da funcao f3 (t) do exemplo 2.2,
reapresentada na Figura 2.10 abaixo, com 0 = 1 rad/seg.
A express
ao geral dos coeficientes da serie, a partir da equacao (2.23), fica:
Sistemas Lineares
53
Fk =
1
2.
1
2
.ejkt .dt +
0
1 jkt
.e
.dt
2
Resolvendo a integral:
Fk =
1
. sin(k.t).dt
2
1
j
=
90
k.
k.
Entao, a serie exponencial de Fourier para esta funcao e escrita abaixo, apresentando
apenas os termos para 5 k 5.
j j3t
j j5t
j
j j3t
j j5t
j
.e
.e + .ejt +
.e
+
.e
f3 (t) = .ejt
3.
5.
3.
5.
Para relacionar este resultado com `
aquele obtido com a serie trigonometrica do exemplo
2.2, combinam-se os termos exponenciais da equacao anterior para que a relacao de Euler
possa ser utilizada. Com isso, o resultado em termos de seno e cosseno fica:
f3 (t) =
1
1
2
. sin t + . sin 3t + . sin 5t + . . .
3
5
Nota-se que este resultado, conforme esperado, e identico ao resultado obtido no exemplo 2.2.
2.5
O Espectro de Freq
u
encias de Sinais Peri
odicos
Sistemas Lineares
54
dos coeficientes de serie) e o espectro de fase (fase dos coeficientes). Na maioria dos casos
os coeficientes sao reais, sendo apenas necessario o espectro de amplitude das componentes de freq
uencia.
Exemplo 2.4 Faca a expans
ao em serie exponencial de Fourier da funcao porta periodica
mostrada na Figura 2.11. Plote o espectro de freq
uencia.
f (t) =
Fk =
1
T
T /2
f (t).ejk0 t .dt
/2
/2
A.ejk0 t .dt
=
/2
A jk0 t
.e
jk0 T
/2
(2.34)
/2
Sistemas Lineares
55
O Gr
afico de Sa(x) e apresentado na Figura 2.12. Esta funcao oscila com perodo 2,
com amplitude decrescente em ambas direcoes de x e tem zeros em x = , 2, . . . ,
etc.
Sa(x)
x
6 5 4 3 2
sin(x)
.
x
A.
.Sa(k0 /2)
T
Como
0 =
2.
T
k..
k.0
=
2
T
pode-se reescrever a equacao dos coeficientes da serie de Fourier como:
Fk =
A.
.Sa
T
k..
T
k=
Sa
k..
T
.ejk0 t
Sistemas Lineares
56
etc., com amplitude A./T, (A./T ).Sa(./T ), (A./T ).Sa(2../T ), . . . , etc., respectivamente.
O grafico do espectro de freq
uencia ser
a tracado para uma largura do pulso de 1/20
segundo. O perodo T ser
a tomado como 1/4, 1/2 e 1 segundo sucessivamente.
Para = 1/20 e T = 1/4 a freq
uencia fundamental e 0 = 8. e o espectro de freq
uencia
e dado por:
k.
A
Fk = .Sa
5
5
Nota-se que como 0 = 8., o espectro existe em = 0, 8., 16. 24., . . . ,
etc., conforme pode ser visto na Figura 2.13(a).
Para = 1/20 e T = 1/2 a freq
uencia fundamental e 0 = 4. e o espectro de freq
uencia
e dado por:
k.
A
Fk = .Sa
10
10
Este espectro existe apenas em = 0, 4., 8. 12., . . . , etc. e pode ser visto
na Figura 2.13(b).
Finalmente, para T = 1 a freq
uencia fundamental e 0 = 2. e o espectro existe apenas
em = 0, 2., 4. 6., . . . , etc. dado pela equacao abaixo e visto na Figura
2.13(c).
Sistemas Lineares
57
A
5
1
,
20 T
= 51 , T =
1
4
120 100 80 60 40 20
20
40
60
80
100
120
(a)
A
10
1
,
20 T
1
,
10
T =
1
2
120 100 80 60 40 20
20
40
60
80
100
120
(b)
A
20
1
,
20 T
1
,
20
T =1
120 100 80 60 40 20
20
40
60
80
(c)
100
120
Sistemas Lineares
58
Bibliografia
PHILLIPS, C. L. Signals, Systems, and Transforms, Prentice-Hall Inc.:New Jersey,
USA, 1995, 707p.
CLOSE, C. M. The Analysis of Linear Circuits, Harcourt, Brace & World, Inc.:New
York, USA, 1966, 716p.
LATHI, B. P. Sistemas de Comunicacao, Editora Guanabara, S. A.:Rio de Janeiro, RJ,
1979, 401p.
Sistemas Lineares
60
Captulo 3
Transformada de Fourier
Neste captulo estudaremos a transformada de Fourier e como essa transformada pode
nos auxiliar na analise de sinais e suas propriedades. A transformada de Fourier transforma um sinal f (t) numa funcao complexa F (), conhecida como espectro do sinal f (t).
As equacoes que definem a transformacao de variaveis (f (t) para F () e vice-versa) sao
dadas em (3.1). Estas operacoes, conhecidas como transformada de Fourier e sua inversa
sao ilustradas na figura 3.1.
F[f (t)]
F ()
f (t)
F 1 [F ()]
f (t) = F 1 [F ()] =
F () = F[f (t)] =
1
2
F ()ejt d
(3.1)
jt
f (t)e
dt
Existem sinais f (t) para os quais nao e possvel se calcular a Transformada de Fourier.
Uma condicao suficiente para a existencia da Transformada de Fourier e indicada a
seguir:
|F ()| = |
jt
f (t)e
dt|
jt
|f (t)e
| dt =
|f (t)| dt <
Logo, se o sinal f (t) e integravel em modulo a sua transformada vai seguramente existir.
No entanto o contrario nao e verdade em geral, pois sinais como seno, cosseno, e degrau
nao sao integraveis em modulo mas suas transformadas existem como casos limites e
podem ser expressas com o auxlio de funcoes impulsos.
Sistemas Lineares
3.1
62
Energia de sinais
E=
f (t)2 dt
(3.2)
Por exemplo, se f (t) representa a tensao ou corrente num resistor unitario, a energia do
sinal f (t) e dada pela integral acima. Os sinais que possuem energia limitada (E < )
sao portanto de grande interesse pratico.
Veremos a seguir que a energia de um sinal esta ligada ao modulo da Transformada de
Fourier do sinal em questao. Antes porem, devemos relembrar algumas propriedades da
funcao F () dadas no exemplo 3.1.
Exemplo 3.1 Mostre que |F ()| e uma funcao par de e F () e uma funcao mpar
de .
Solucao: Seja F () = M ()ej() onde M () e () denotam respectivamente o
modulo e a fase de F (). Note que f (t) e um sinal real logo f (t) e seu conjugado
ao iguais, isto e f (t) = f (t). Entao a conjugacao complexa de F ()
complexo f (t) s
resulta:
F () =
f (t)ejt dt =
f (t)ejt dt =
f (t)ejt dt = F ()
as conclus
oes desejadas:
Portanto F () = F () = M ()ej() que nos leva `
|F ()| = |F ()| = M () (funcao par)
F () = F () (funcao mpar)
Podemos agora mostrar que a energia de um sinal esta ligada ao modulo da Transformada de Fourier do sinal em questao. Com f (t) de (3.1) obtemos:
E =
1
=
2
f (t)2 dt =
f (t) f (t) dt =
f (t)
1
f (t)ejt dt d =
F ()
2
1
2
F ()ejt d dt
F ()F ()d
1
2
|F ()|2 d
(3.3)
Sistemas Lineares
63
A formula acima e conhecida como teorema de Parseval. Note ainda que como |F ()| e
uma funcao par temos:
1
E=
|F ()|2 d
(3.4)
0
A quantidade
|F ()|2
ET =
f (t) dt =
1
2
1
1+j
0
e2t dt =
1
2
|F ()|2 d =
0
0
1
1
d = tan1 (0 )
2
1+
4
Logo metade da energia do sinal est
a no intervalo de freq
uencia entre [1, 1].
3.2
3.2.1
C
alculo de algumas transformadas
Sinal Exponencial Unilateral
a>0
Entao:
F () = F[f (t)] =
eat u(t)ejt dt =
0
e(a+j)t dt =
1
a + j
Sistemas Lineares
64
G (t)
1
t
2
3.2.2
Sinal Porta
Usaremos a notacao G (t) para definir o sinal porta (gate) de largura como indicado
a seguir.
1, |t| < /2
G (t) =
0, |t| > /2
A transformada do sinal porta e calculada da seguinte forma.
/2
1 j /2
(e
ej /2 )
j
/2
sen( /2)
= Sa [ ]
=
( /2)
2
ejt dt =
F () = F[G (t)] =
A funcao Sa(x) =
na figura 3.3.
sen(x)
x
)
2
(3.5)
Note que F () nesse caso e real pois f (t) e par. Se f (t) for mpar entao F () e
imaginario puro.
3.2.3
Sinal Impulso
A notacao para a funcao impulso unitario que ocorre no instante zero e (t). Como
(t) = 0 para t = 0 temos:
(t)f (t)dt =
(t) dt = f (0)
Logo:
F () = F[(t)] =
(t)ejt dt = ej0 = 1
(t) 1
Sistemas Lineares
65
Sa(x)
1.1
1
0.9
0.7
0.5
0.3
0.1
-0.1
-0.3
-16
-12
-8
-4
12
16
sen(x)
x
f (t) =
K
Sa (Kt)
Podemos mostrar que a area dessa funcao e unitaria para qualquer valor do parametro
K, isto e,
K
Sa (Kt) dt = 1, K
Com esse resultado podemos ainda mostrar que quando K tende `a infinito a funcao f (t)
tende `a funcao impulso.
K
Sa (Kt)
K
(t) = lim
3.2.4
(3.6)
Func
oes Constante, Sinal e Degrau
A transformada do degrau nao pode ser facilmente obtida pela aplicacao da definicao.
A seguir veremos como obte-la com o auxlio das transformadas das funcoes sinal e
constante.
A funcao constante unitaria pode ser vista como o caso limite da funcao porta, isto e
Sistemas Lineares
66
lim Sa ( /2)
Sa ( /2)
= 2 lim
2
= 2()
(3.7)
sgn(t) =
sgn(t)
1
0
-1
2
2j
=
2
2
a0 a +
j
= lim
(3.8)
O resultado acima nos permite calcular agora a Transformada do degrau. Como u(t) =
1
(1 + sgn(t)) temos:
2
1
1
1
F[u(t)] = F[1] + F[sgn(t)] = () +
2
2
j
Resumindo:
1 2()
2
sgn(t)
j
1
u(t) () +
j
(3.9)
Sistemas Lineares
3.2.5
67
Sinais Senoidais
Nos ocuparemos agora das transformadas das funcoes senoidais. Pela definicao temos
F[cos(0 t)] =
ej0 t +ej0 t
2
como cos(0 t) =
T /2
jt
cos(0 t)e
cos(0 t)ejt dt
dt = lim
T /2
temos:
T
( 0 )
( + 0 )
{Sa [T
] + Sa [T
]}
2
2
2
(3.10)
(3.11)
Note que a Transformada de Fourier do sen(0 t) e cos(0 t) so nao e nula nas freq
uencias 0 . Isto mostra que esses sinais possuem energia concentrada nessas freq
uencias.
Isso nao ocorreria se as funcoes fossem sen(0 t)u(t) ou cos(0 t)u(t). Nesse caso obteramos:
j
F[cos(0 t)u(t)] = [( 0 ) + ( + 0 )] + 2
2
0 2
0
F[sen(0 t)u(t)] = [( 0 ) ( + 0 )] + 2
2j
0 2
3.2.6
Exponencial Eterna
3.2.7
(3.12)
Func
oes Peri
odicas
k=
Fk .ejk t
t0 < t < t0 + T
(3.13)
Sistemas Lineares
68
F0 =
1
Fn =
T
1
T
t0 +T
t0 +T
f (t)dt
(3.14)
f (t)ejn t dt
(3.15)
t0
t0
F[f (t)] =
k=
jk0 t
Fk F[e
] = 2
k=
Fk .( k.0 )
(3.16)
Problema 3.2 Mostre que se f (t) e uma funcao par, isto e f (t) = f (t), ent
ao Fn
e F[f (t)] s
ao ambos reais e se f (t) e mpar, isto e f (t) = f (t), Fn e F[f (t)] s
ao
jx
puramente imagin
arios. Sugest
ao: use a formula de Euler e = cos(x) + jsen(x).
Sistemas Lineares
69
f (t)
1
2
t
0
-1
T
Figura 3.5: Funcao onda quadrada de perodo 2.
1
T
j
=
2
Fk =
j
=
2
t0 +T
2
1
f (t)[cos(k t) jsen(k t)]dt
2 0
t0
2
2
j
f (t)sen(k t)dt =
sen(k0 t) dt
sen(k0 t)dt +
2 0
0
f (t) ejk t dt =
cos(k0 t)
k0
cos(k0 t)
+
k0
j
(cos(k) + 1 + cos(k2) cos(k))
2k
j
(2 2cos(k))
=
2k
Ou seja:
2j
Fk = k.
0
, se k e mpar
(3.17)
, se k e par
2j
4j
F[f (t)] = 2
( k) =
( k)
k
k
k=
k=
O resultado demonstra que a Transformada de Fourier de um sinal periodico existe
apenas nas freq
uencias das harmonicas que o compoem. Isto pode ser visto no grafico do
modulo da transformada do sinal, dado pela Figura 3.6.
Sistemas Lineares
70
|F (j)|
4
4
3
4
5
4
3
6 5 4 3 2 1
4
5
[rad/s]
T (t) =
k=
(t kT )
F[T (t)] = 2
k=
Fk .k.( k0 )
onde Fk , k = 1, 2, 3, ... s
ao os coeficiente harmonicos do sinal que s
ao dados por:
Fk =
1
T
T /2
T (t)ejk0 t dt =
T /2
Logo
F[T (t)] =
2
T
1
T
k=
T (t)
T /2
(t)ejk0 t dt =
T /2
( k0 ) = 0 0 ()
F [T (t)]
(t nT )
...
...
T
1
1 jk0 0
e
=
T
T
0 ( n0 )
...
...
T (t) 0 0 ()
3.3
Propriedades da transformada
Em primeiro lugar, vale a pena salientar que para funcoes integraveis em modulo podemos obter a transformada de Fourier diretamente da transformada de Laplace com
Sistemas Lineares
71
3.3.1
Linearidade
3.3.2
Simetria
t
)] = 2 G ()
2
3.3.3
G2 ()
(3.18)
Escalonamento
1
F( )
|a| a
Sistemas Lineares
3.3.4
72
Deslocamento em freq
u
encia e Modula
c
ao
Ao multiplicar um sinal f (t) pelo cos(0 t) estamos atenuando pela metade e deslocando
todo o espectro do sinal f (t) para as frequencias 0 . Este artifcio e conhecido como
modulacao em amplitude pois o sinal f (t), conhecido como sinal modulado, e a amplitude
do cosseno. A funcao cos(0 t) recebe o nome de portadora de freq
uencia 0 . A figura
3.8 mostra a transformada de Fourier do sinal porta de largura unitaria G1 (t). A figura
3.9 ilustra a o espectro do sinal modulado cos(100t)G1 (t).
1.1
F [G1 (t)]
0.9
0.7
0.5
0.3
0.1
-0.1
-0.3
-100
-80
-60
-40
-20
20
40
60
80
100
Sistemas Lineares
0.5
73
F [cos(100t)G1 (t)]
0.4
0.3
0.2
0.1
-0.1
-200
-160
-120
-80
-40
40
80
120
160
200
f(t)
MOD
cos(0 t)
f (t)cos(0 t)
MOD
f (t)cos2 (0 t)
FILTRO
IDEAL
f(t)
3.3.5
Deslocamento no Tempo
Sistemas Lineares
3.3.6
74
Diferenciac
ao e Integrac
ao no Tempo
F[
f ( )d ] =
df (t)
] = jF ()
dt
1
F () , se F () = 0 para = 0
j
A restricao F () = 0 para = 0 implica que o valor medio do sinal deve ser nulo, isto
e f (t)dt = 0. Essa restricao pode ser eliminada mas a expressao acima se torna mais
complicada.
Exemplo 3.6 Obtenha a Transformada de Fourier do sinal f (t) da figura 3.11.
f(t)
A
a
b
b
t
A
ba
Sistemas Lineares
75
A
ba
d2 f (t)
dt2
A
ba
(t + b)
A
ba
(t b)
(t + a)
A
ba
(t a)
A
2A
d2 f
jb
ja
ja
jb
]
=
[e
e
+
e
]
=
(cos(b) cos(a))
dt2
ba
ba
3.3.7
2A cos(a) cos(b)
ba
2
Diferenciac
ao em freq
u
encia
3.3.8
dF ()
d
Convoluc
ao
f1 ( )f2 (t )d
Analogamente `a transformada de Laplace podemos transformar a integral de convolucao em produto no domnio da freq
uencia. Seja F[f1 (t)] = F1 () e F[f2 (t)] = F2 ()
entao:
F[f1 (t) f2 (t)] = F1 ()F2 ()
A prova desse resultado e bastante simples.
F[f1 (t) f2 (t)] =
=
f1 ( )f2 (t )d ejt dt
f1 ( )
f2 (t )ejt dt d
Sistemas Lineares
como
76
3.4
1
F1 () F2 ()
2
(3.19)
Amostragem
A/D
amostragem e
codificac
ao
transmiss
ao
.....
D/A
sinal recebido
decodificac
ao e
reconstruc
ao
Sistemas Lineares
77
f (t)
t
(a) Sinal contnuo
T a (t)
Ta
t
(c) Sinal amostrado
k=
( k.a ) = a .a ()
1
F[f (t)] F[Ta (t)]
2.
1
Ta
k=
F () ( k.a ) =
( k.a )
k=
1
Ta
k=
F ( k.a )
Sistemas Lineares
78
F ()
A
B a B
a + B
segundos
Fs ()
F ()
A
Ta
A
...
...
0
a
a
Sistemas Lineares
79
G()
filtro ideal
Fs ()
Ta
F ()
G()
Fs ()G() = F ()
No domnio do tempo a filtragem acima e dada pela convolucao:
f (t) = fs (t) g(t)
onde g(t) = F 1 [G()] = Sa (
t). Este resultado pode ser obtido facilmente com (3.18) e
Ta = /
. Assim, ficamos com
f (t) =
=
=
k=
k=
k=
f (k.Ta )Sa (
(t k.Ta ))
A expressao acima mostra como se reconstroi exatamente o sinal f (t) a partir das amostras f (k.Ta ) coletadas no processo de amostragem. As amostras formam um conjunto de
amplitudes de funcoes sampling que quando somadas resultam no sinal original.
Note que a reconstrucao exata requer um filtro ideal o que nao pode ser implementado
na pratica. Apesar disso podemos obter boas aproximacoes do sinal a ser reconstrudo
substituindo-se o filtro ideal por um real que tenha uma funcao de transferencia cujo
espectro seja parecido com G().
Vejamos agora o que acontece quando a =
Caso a < 2
ou seja Ta >
2
Ta
< 2
.
segundos
Nesse caso o espectro do sinal antes e apos amostragem estao ambos ilustrados na
figura 3.20.
Note agora que o espectro do sinal amostrado Fs () contem o espectro do sinal original
F () porem distorcido com as superposicoes dos espectros deslocados. Essa distorcao
provocada pela superposicao dos espectros inviabiliza a reconstrucao do sinal e portanto
na escolha da freq
uencia de amostragem deve-se evitar o caso a < 2
.
Sistemas Lineares
80
Fs ()
F ()
A
Ta
A
...
...
a
0
a
Sistemas Lineares
81
F ()
j ( + 10) ( 100)
( + 100)
j ( 10)
1
R()
F ()
Ta =
Rs ()
y(t)
Ta F ()
1
80
70
70
1
G100 () = G100 ()
50
50
O sinal na sada do primeiro filtro e dado por R() = X()F () de onde conclumos
1
F ()
X()
1
50
X()
50
50
70 50
50
70
(a)
R() = X()F () = X()
(b)
que R() = X(). Veja figura 3.23. O espectro Rs () do sinal amostrado Rs (t) est
a
indicado na figura 3.24 e e dado por
1
Rs () =
Ta
onde a =
2
Ta
1
R( k.a ) =
Ta
k=
= 160 e a freq
uencia de amostragem.
k=
X( k.a )
Sistemas Lineares
82
O espectro do sinal ap
os o segundo filtro e dado por
Y () = Ta F ()Rs () =
Ta F ()
k=
1
X( k.a ) = X()
Ta
O produto Ta F ()Rs () pode ser facilmente obtido atraves da figura 3.24. Logo conclumos que y(t) = x(t).
Rs ()
Ta F ()
1
X(
Ta
+ 2a )
1
X(
Ta
+ a )
1
X()
Ta
1
X(
Ta
a )
1
X(
Ta
2a )
...
...
2a
a
70
50
a
50
2a
70
3.5
3.5.1
Aplicac
oes da Transformada de Fourier
Resposta em freq
u
encia de sistemas lineares
Sistemas Lineares
83
V2 (t) = L.
di(t)
dt
V2 ()
jL
Manipulando algebricamente estas duas equacoes pode-se chegar a relacao entradasada do sistema:
V2 ()
jL
H() =
=
V1 ()
R + jL
A funcao H() e tambem chamada de funcao resposta em freq
uencia do sistema, pois
informa a relacao da sada com a entrada em funcao da freq
uencia.
A resposta em freq
uencia pode ser levantada tambem experimentalmente. Sabendo
que a resposta em freq
uencia e um n
umero complexo:
H() = |H()|
()
cos(x .t)
A. cos(x .t + )
H()
3.5.2
H()|=x
Sistemas Lineares
84
x(t)
h(t)
1
X()
Y ()
+ j =
1
Y ()
=
H() =
X()
1 + j.
Se x(t) = u(t) ent
ao:
X() = .() +
1
j
Logo:
Y () =
1
1
+ .() .
j
1 + j.
j.( 1
K1
+
+ .() =
j
+ j)
K2
+ .()
+ j
Sistemas Lineares
85
Onde:
j. 1
j.( 1 + j)
K1 =
(1 + j. ). 1
j.( 1 + j)
K1 =
Assim,
Y () =
=1
j=0
= 1
j= 1
1
+ .()
j
1
+ j
3.5.3
Determinac
ao da resposta de sistemas LIT a sinais peri
odicos
x(t) =
Fk .ejk0 t
k=
jk0 t
Fk .e
k=
k=
Fk . F ejk0 t
k=
Fk .( k.0 )
k=
Fk .( k.0 )
k=
Sistemas Lineares
86
Como F 1 [2..( k.0 )] = ejk0 t , a sada do sistema em funcao do tempo para uma
entrada periodica fica:
y(t) =
k=
Com isso, conclui-se que o sinal de sada e tambem periodico, com mesmo perodo
T = 20 do sinal de entrada e os coeficientes da serie de Fourier do sinal de sada sao
iguais aos da entrada, multiplicados pela resposta em freq
uencia do sistema na freq
uencia
considerada.
3.5.4
Determinac
ao da resposta impulsiva de sistemas LIT
Quando um sistema e descrito por uma equacao diferencial, sua resposta impulsiva
tambem pode ser determinada com o uso da Transformada de Fourier . Considere o
seguinte sistema de primeira ordem:
a1 .
dy(t)
dx(t)
+ a0 .y(t) = b1 .
+ b0 .x(t)
dt
dt
A Transformada de Fourier pode ser utilizada para a determinacao da resposta impulsiva deste sistema aplicando-a na equacao diferencial acima. Desta forma:
(a1 .j + a0 ).Y () = (b1 .j + b0 ).X()
b1 .j + b0
Y ()
= H() =
X()
a1 .j + a0
Calculando-se a transformada inversa da resposta em freq
uencia H(), chega-se a
resposta impulsiva do sistema:
a0
b0 b1 . a 1
b1
b1 .j + b0
h(t) = F 1 [H()] = F 1
= F 1 +
a1 .j + a0
a1 a . j + a0
1
h(t) =
b1
.(t) +
a1
b0 b1 .a0
2
b1
a1
di(t)
+ Vo (t)
dt
a1
Sistemas Lineares
87
Vo
,
R
L dVo (t)
.
+ Vo (t)
R dt
L
.j.Vo () + Vo () =
R
j.
L
+ 1 .Vo ()
R
R/L
Vo ()
= H() = R
Vi ()
+ j
L
Pela transformada inversa chega-se a resposta impulsiva:
h(t) =
R (R/L)t
.e
.u(t)
L
Sistemas Lineares
88
Bibliografia
PHILLIPS, C. L. Signals, Systems, and Transforms, Prentice-Hall Inc.:New Jersey,
USA, 1995, 707p.
CLOSE, C. M. The Analysis of Linear Circuits, Harcourt, Brace & World, Inc.:New
York, USA, 1966, 716p.
LATHI, B. P. Sistemas de Comunicacao, Editora Guanabara, S. A.:Rio de Janeiro, RJ,
1979, 401p.
Sistemas Lineares
90
Captulo 4
Diagramas de Bode
Neste captulo sera apresentado o tipo de grafico mais utilizado para a representacao de
uma funcao complexa do tipo H(j) sao os diagramas de Bode. Estes graficos se consagraram com os trabalhos de Bode sobre amplificadores realimentados na decada de 1940 e
hoje sao muito utilizados na analise de sinais e sistemas. Nesses diagramas representa-se
o modulo em decibel e a fase em graus, ambos em funcao da freq
uencia (tipicamente em
Hertz) numa escala logartmica. Assim, os diagramas de Bode apresentam a resposta em
freq
uencia de um determinado sistema.
4.1
Introduc
ao
K (j + z1 ) (j + z2 ) . . . (j + zn )
(j + p1 ) (j + p2 ) . . . (j + pn )
(4.1)
onde os termos zn sao os chamdados zeros, que sao os valores que anulam o numerador
da equacao. Ja os termos pn sao os polos, termos que anulam o denominador.
Deseja-se estudar a resposta deste sistema a sinais senoidais. Considere entao entradas
senoidais do tipo
x(t) = A. cos(0 t)
Aplicando a transformada de Fourier, chega-se a uma expressao para X():
X() = A.. [( 0 ) + ( + 0 )]
Sistemas Lineares
92
A resposta em freq
uencia e entao dada por:
Y () = A.H()..( 0 ) + H()..( + 0 )]
Y () = A.H(0 ).. [( 0 ) + ( + 0 )]
(4.2)
A equacao (4.2) mostra que para uma entrada do tipo senoidal, a sada em regime
permanente de um sistema LIT tambem e senoidal, com uma amplitude que depende do
modulo da funcao resposta em freq
uencia e com um angulo de fase igual ao angulo de
H()|=0 .
Exemplo 4.1 Suponha um sistema com a funcao de resposta em freq
uencia dada abaixo:
1
H() =
1 + j
O modulo e a fase da funcao resposta em freq
uencia s
ao:
1
e () = tan1
|H()| =
1 + 22
Plotando o modulo e a fase em funcao da freq
uencia , encontra-se a resposta em
freq
uencia deste sistema, conforme a Figura 4.1.
Normalmente utiliza-se uma escala logartmica para representar |H(0 )|. Ou seja,
tracam-se graficos de log |H(0 )| ou 20 log |H(0 )| utilizando-se a unidade Decibel, abreviada por dB.
A vantagem em se utilizar uma escala logartmica e que tanto o modulo quanto a fase do
produto (ou divisao) de n
umeros complexos sao transformados em soma (ou subtracao)
dos modulos em dB e fases individuais de cada n
umero multiplicado (dividido). Assim,
pode-se analisar separadamente a contribuicao de cada um dos termos. Isso facilita
bastante a construcao manual dos graficos de modulo e fase.
Como exemplo desta propriedade, considere a seguinte funcao resposta em freq
uencia:
1
1
H() =
|H()| =
j.(1 + j)
1 + 22
A representacao de seu modulo em dB e:
1
20 log |H()|dB = log . log(1 + 2 2 )
2
Sistemas Lineares
93
4.2
Sistemas Lineares
94
i=1
log |j zi |
i=1
log |j pi |
(4.4)
[1 sgn(K)]
. +
2
i=1
(4.5)
i=1
4.3
Diagramas b
asicos
4.3.1
Constante
Sistemas Lineares
95
4.3.2
P
olo na origem
Quando a funcao de transferencia possui somente um polo na origem, ela tem a seguinte
forma:
K
H() =
j
Seu modulo em dB fica:
|H()|db = 20 log
K
= 20 log K 20 log
j
Para outros valores de K, a curva e tambem uma linha reta com a mesma inclinacao,
mas deslocada pela constante K, conforme a Figura 4.4.
A fase e constante () = 90 se K for positivo ou H() = 180 90 = 90 se K
for negativo.
Se o polo na origem for m
ultiplo, entao
H() =
K
j n
Sistemas Lineares
96
4.3.3
Zero na origem
Sistemas Lineares
4.3.4
97
1
1 + j
1
j+1
que
Percebe-se pela Figura 4.6 que as assntotas (linhas pontilhadas) possuem, no modulo,
inclinacoes de zero e -20 dB/decada para baixas e altas freq
uencias respectivamente. A
fase vale zero graus nas baixas freq
uencias, -90 graus nas altas freq
uencias e nas medias
freq
uencias pode ser aproximada por uma assntota de inclinacao -45 graus/decada. Aqui
Sistemas Lineares
98
Magnitude
db
-10
-20
-30
0
-10
-20
-30
-40
-50
-60
-70
-80
-90
Hz
1
T
0.1
T
10
T
Phase
degrees
Hz
1
T
0.1
T
10
T
1
T s+1
e assntotas
1
,
1j
Vo ()
1/jC
1
=
1 =
V1 ()
1 + jRC
R + jC
Sistemas Lineares
99
A freq
uencia de corte neste caso e dada por c =
exemplo est
a tracado na Figura 4.8.
1
.
RC
= 20 log
1+
dB
1
j 2RC
RC
1
2RC
= 20 log
= 1dB
1 + j/2
= 26, 6
Para = c :
H
1
RC
= 20 log
1+
dB
1
RC
j RC
1
RC
= 20 log
= 3dB
1+j
= 45
Para = 2c :
H
2
RC
= 20 log
dB
4.3.5
1
1+
2
j RC
RC
2
RC
= 20 log
= 7dB
1 + 2j
= 63, 4
1
1 + 2 jn +
j
n
(4.6)
Sistemas Lineares
100
2
n2
+ 4 2
2
n2
2/n
1 2 /n2
Diagramas Assint
oticos
Usando o mesmo raciocnio do caso de primeira ordem, pode-se analisar as assntotas
para valores extremos de freq
uencia. Assim, verifica-se que:
2
n
Em baixas freq
uencias, lim0 (j)2 +2
2 = 1, ou seja, GdB = 0. Logo nas
n (j)+n
baixas freq
uencias o termo se comporta como um fator constante unitario.
2
n
n
Em altas freq
uencias, lim (j)2 +2
2 = (j)2 ( ), ou seja, GdB =
n (j)+n
40 log /n . Logo nas altas freq
uencias o termo se comporta como um fator
2
n
do tipo (j)2 que possui fase -180 graus e modulo decrescendo na razao de -40
dB/decada.
Em medias freq
uencias, na freq
uencia de corte = n temos
n2
1
=
2
2
(j) + 2n (j) + n
j2
que possui modulo 20 log(2) dB e fase tan1 () = 90 graus.
Pela Figura 4.9 verifica-se que as assntotas (linhas pontilhadas) possuem, no modulo,
inclinacoes de zero e -40 dB/decada para baixas e altas freq
uencias respectivamente. A
fase vale zero graus nas baixas freq
uencias, -180 graus nas altas freq
uencias e nas medias
freq
uencias pode ser aproximada por uma assntota de inclinacao -90 graus/decada. Aqui
consideramos medias freq
uencias o intervalo entre uma decada abaixo e uma decada acima
da freq
uencia de quebra.
A frequencia e o pico de ressonancia sao calculados da seguinte forma:
d
d
|G(j)| = 0
d
d
1 ( n )2
+ 2 n
=0
2
2
2
r = n 1 2 , 0
2
(4.7)
Sistemas Lineares
20
101
Magnitude
db
= 0.1
10
0
assntotas
-10
=5
-20
-30
-40
Hz
-50
0.1 n
0
-20
-40
-60
-80
-100
-120
-140
-160
-180
10 n
Phase
degrees
assntotas
=5
= 0.1
Hz
10 n
0.1 n
Se >
`a zero.
2
n
2
(j)2 +2n j+n
e assntotas
Quando 0
1
2
1 2
(4.8)
Sistemas Lineares
102
Sistemas Lineares
103
x(t)
y(t)
I(t)
-
Entrada
Sada
SISTEMA
F (j) =
n2
1
=
LC(j)2 + RCj + 1
(j)2 + 2n j + n2
Logo:
n2 =
R
R
1
; 2n =
=
LC
L
2
C
L
10(j + 10)
j(j + 1)(j 2 + 100j + 104 )
Sistemas Lineares
104
102 (0, 1j + 1)
j(j + 1)(104 j 2 + 102 j + 1)
Magnitude
db
-20
-30
-40
-50
rad/s
-60
-1
10
0
-10
-20
-30
-40
-50
-60
-70
-80
-90
-1
10
10
10
10
10
Phase
degrees
rad/s
0
10
10
10
10
0.01(0.1j+1)
j
e assntotas
Sistemas Lineares
-20
-30
-40
-50
-60
-70
-80
-90
-100
-110
-120
105
Magnitude
db
rad/s
-1
10
-90
10
10
10
10
Phase
degrees
-100
-110
-120
-130
-140
rad/s
-150
-1
10
10
10
10
10
1
Figura 4.13: Diagrama de Bode do termo G2 (j) = G1 (j) j+1
e assntotas
-20
Magnitude
db
-40
-60
-80
-100
-120
-140
rad/s
-160
-1
10
-90
-110
-130
-150
-170
-190
-210
-230
-250
-270
-1
10
10
10
10
10
Phase
degrees
rad/s
0
10
10
10
10
1
Figura 4.14: Diagrama de Bode do termo G(j) = G2 (j) 104 j2 +10
ntotas
2 j+1 e ass
Sistemas Lineares
4.4
106
Nesta secao serao estudadas algumas propriedades associadas aos zeros da funcao de
transferencia.
Defini
c
ao 4.1 Um sistema e dito ser de Fase Mnima se todos os zeros da funcao de
transferencia desse sistema est
ao no semi-plano complexo esquerdo. Caso contr
ario, isto
e se existir algum zero no semi-plano direito ou sobre o eixo imagin
ario, o sistema e dito
ser de Fase Nao-mnima.
Para que um sistema de controle tenha algum interesse pratico ele deve ser estavel,
isto e todos os zeros da sua funcao de transferencia devem ter parte real estritamente
negativa. No entanto alguns sistemas fsicos estaveis podem possuir zeros no semi-plano
direito.
Exemplo 4.5 O circuito da figura 4.15 possui x(t) como tensao de entrada e y(t) como
tensao de sada.
r2
r1
r1
+
C
r2
1 r0 Cj
1 + rCj
(4.9)
1
Note que G(s) possui um polo em s = rC
e um zero em s = r01C . Portanto o sistema
e est
avel de fase n
ao mnima se escolhemos r2 > r1 , pois nesse caso o zero de G(s) est
a
no semi-plano direito. Se escolhemos r2 < r1 o sistema e est
avel de fase mnima pois
agora o zero de G(s) est
a no semi-plano esquerdo. O diagrama de Bode desse sistema e
indicado na figura 4.16 para caso (a): r1 = 10K, r2 = 20K, C = 1F e na figura
4.17 para caso (b): r2 = 10K, r1 = 20K, C = 1F . Veja que o diagrama de modulo
e igual para os dois casos mas o diagrama de fase e diferente.
Sistemas Lineares
107
Magnitude
db
0
-1
-2
-3
-4
-5
-6
-7
-8
-9
-10
Hz
-1
10
10
10
10
Phase
degrees
0
-20
-40
-60
-80
-100
-120
-140
-160
-180
10
Hz
-1
10
10
10
10
10
Figura 4.16: Caso (a): Sistema de fase nao mnima (r2 > r1 )
Magnitude
db
0
-1
-2
-3
-4
-5
-6
-7
-8
-9
-10
Hz
-1
10
10
10
10
10
Phase
degrees
-2
-6
-10
-14
-18
-22
-26
Hz
-30
-1
10
10
10
10
10
Sistemas Lineares
108
Sistemas de fase mnima possuem propriedades bastante interessantes. Sao mais simples de serem controlados e os seus diagramas de Bode (modulo e fase) sao assintoticos
nas altas e baixas freq
uencias e alem disso podemos relacionar a assntota de modulo com
a de fase atraves do grau relativo do sistema. Grau relativo de um sistema e a diferenca
de grau entre o denominador e o numerador da funcao de transferencia do mesmo.
Veja o que ocorre se for considerado um sistema que possui uma funcao de transferencia
do tipo:
K(am (j)m + + a1 j + 1)
(4.10)
G(s) =
bn (j)n + + b1 j + 1
com an , bn , K reais positivos.
O grau relativo desse sistema e n m. Note que para todo sistema de interesse pratico
o grau relativo e sempre positivo, ou seja, n m 0.
Nas baixas freq
uencias tem-se que:
lim G(j) = K
am
(j)mn
bn
Bibliografia
CLOSE, C. M., The Analysis of Linear Circuits, Harcourt, Brace & World, Inc.:New
York, USA, 1966, 716p.
DAZZO, J. J.; HOUPIS, C. H., An
alise e Projeto de Sistemas de Controle Lineares,
Editora Guanabara:Rio de Janeiro, RJ, 1984, 660p.
Sistemas Lineares
110
Captulo 5
Transformada de Laplace
5.1
Introduc
ao e No
c
oes de Func
oes Complexas
V(t)
Vc (t)
dv(t)
= v(t)
dt
Sistemas Lineares
112
LAPLACE
F(s)
Transf. Inversa
PROPRIEDADES DE FUNC
OES
COMPLEXAS:
Neste curso vamos nos restringir, com poucas excessoes, `a funcoes complexas racionais.
Defini
c
ao 5.1 (Fun
c
ao Racional) Uma funcao G(s) da vari
avel complexa s = +j e
racional se G(s) pode ser expressa como a divisao de dois polin
omios da vari
avel complexa
s.
A figura abaixo ilustra uma funcao complexa G(s) em termos de suas coordenadas
retangular e polar. onde |G(s)| = G2x + G2y e G(s) = tan1 Gy /Gx .
Im[G(s)]
Gy
Gx
Re[G(s)]
Sistemas Lineares
113
O n
umero complexo:
ej = cos + jsen
cos2 + sen2 = 1
ej = tan1
sen
=
cos
Defini
c
ao 5.3 (Fun
c
ao Analtica) Uma funcao G(s) e analtica numa regiao se G(s)
e todas as suas derivadas existem nessa regiao.
1
s+1
Sistemas Lineares
5.2
114
Definic
ao e Regi
ao de Converg
encia
Para uma funcao f (t) com t 0, define-se Transformada de Laplace de f (t) como
sendo a funcao complexa F (s) obtida atraves da integral:
f (t)est dt
(5.1)
f (t) = L1 [F (s)] =
1
2j
c+j
F (s)est ds
(5.2)
cj
onde t 0 e c e um n
umero real associado `a regiao do plano s = + j onde a funcao
F (s) esta definida. Esta regiao e chamada regiao de convergencia da Transformada de
Laplace . Dentro dessa regiao as funcoes f (t) para t 0 e F (s) estao ligadas de maneira
biunvoca, como ilustra a figura a seguir.
Trans. Direta
F (s)
Re[s] > c
f (t)
t0
Tranf. Inversa
1 (s2)t
e
|0
s2
0
1
1
1
=
[ lim e(s2)t lim e(s2)t ] =
lim e(s2)t
t0
s 2 t
s 2 s 2 t
e2t est dt =
Note que s = + j e
|ejt | = |cost + jsent| = 1.
Assim,
lim e(s2)t
0
para Re[s] = > 2.
Sistemas Lineares
115
1
s2
A regiao do plano complexo onde a Integral de Laplace esta definida e e finita recebe o nome de regiao de convergencia da Transformada de Laplace . Mostra-se que ao
escolhermos um contorno para a integral:
1
2j
c+j
F (s)est ds
cj
de tal forma que c > 2 (contorno dentro da regiao de convergencia) entao o resultado da
integral acima e e2t para t 0.
2
eat est dt =
1 (sa)t
1
e
|0 =
sa
sa
L[u(t)] =
0, t < 0
1, t 0
1 st 1
e |0 = .
s
s
1est dt =
L[f (t)] = A
( udv = uv
vdu)
test dt = At
est
|
s 0
0, t < 0
At, t 0, A constante
A
Aest
dt =
s
s
est dt =
A
.
s2
Sistemas Lineares
116
0,
t<0
sen(0 t), t 0, 0 cte
L[f (t)] =
sen(0 t)est dt =
1
1
1
=
2j s j0 s + j0
ej0 t ej0 t st
e dt
2j
0
= 2
s + 02
RESUMO
u(t) 1s : Polo simples na origem. Funcao Constante no tempo.
tu(t)
1
:
s2
1
et u(t) s+
: Polo em s = . Cresce exponencialmente no tempo se polo for positivo
( < 0). Decresce exponecialmente no tempo se polo for negativo ( > 0). Valor
constante no tempo se o polo for na origem.
0
sen(0 t)u(t) s2+
olos complexos conjugados sobre o eixo imaginario (s = j0 ).
2 : P
0
Funcao oscila no tempo sem amortecimento.
5.3
Propriedades
5.3.1
Operac
ao Linear
(1 f1 (t) + 2 f2 (t))est dt
= 1
f1 (t)est dt + 2
f2 (t)est dt
Sistemas Lineares
5.3.2
117
Func
ao Transladada em Atraso
Seja f (t) uma funcao, u(t) o degrau unitario e uma constante. Entao:
L[f (t )u(t )] = es L[f (t)]
f (t)
f (t )u(t )
f (t )u(t )est dt
f ( )u( )es( +) d
es
f ( )u( )es d
es
0
f ( )u( )es d
f ( )es d
es
es L[f (t)]
5.3.3
Func
oes Porta-deslocada e Impulso
Sistemas Lineares
118
fp (t)
1
t0
t0
0
1
u(t)
t0
1
u(t
t0
t0 ).
1
1
u(t) u(t t0 )
t0
t0
1
1
=
L[u(t)] L[u(t t0 )]
t0
t0
t0 s
1 e
1 1
=
t0 s t0 s
1
=
(1 et0 s )
t0 s
L[fp (t)] = L
Fun
c
ao Impulso: A Funcao Impulso Unitario que ocorre no instante t = t0 e representada por (t t0 ) e satisfaz as seguintes condicoes:
(t t0 ) =
0, t = t0
, t = t0
(t t0 )dt = 1
A Funcao Impulso e uma abstracao matematica e nao existe na pratica. Porem, variacoes bruscas de energia podem ser aproximadas pela funcao impulso. Alem disso, o
conceito da funcao impulso e bastante u
til na diferenciacao de funcoes descontnuas, como
veremos na sequencia.
Para calcular a transformada da funcao impulso devemos notar que o impulso na origem
e o caso limite da funcao porta quando t0 0, isto e:
(t) = lim
t0 0 t0
[u(t) u(t t0 )]
Sistemas Lineares
119
Assim temos:
1
(u(t) u(t t0 ))
t0 0 t0
1
= lim L
(u(t) u(t t0 ))
t0 0
t0
1
= lim
(1 et0 s )
t0 0 t0 s
d
(1 et0 s )
dt0
=
d
(t s)
dt0 0
= 1
L[(t)] = L lim
f(t)
A (t)
t0
t
t0
0
A (t t0 )
5.3.4
Multiplicac
ao por exponencial
L[e
f (t)] =
f (t)et est dt = F (s + )
Exemplo 5.10 J
a vimos que:
L[sen(0 t)u(t)] =
0
= F (s)
s2 + 02
Sistemas Lineares
120
Logo:
L[et sen(0 t)u(t)] =
0
= F (s + )
(s + )2 + 02
5.3.5
Mudan
ca na Escala de Tempo
Este resultado e u
til quando se deseja analisar sinais numa escala de tempo diferente
daquela em que ele ocorre na pratica. Pode ser o caso por exemplo de sinais muito lentos
ou muito rapidos.
Exemplo 5.11 Dado que L[et u(t)] =
5.3.6
1
s+1
5
.
5s+1
Teorema da Diferenciac
ao Real
df (t) def
= f (t)
dt
(5.3)
def
udv = uv|
0
vdu .
Quando uma funcao possui descontinuidade na origem, a sua derivada temporal ira
possuir um impulso na origem. Nesses casos precisamos tomar cuidado com o limite
inferior da transformada da derivada. Vamos entao definir:
Sistemas Lineares
121
L+ [f (t)] =
L [f (t)] =
0+
f (t)est dt
f (t)est dt
Note que se f (t) envolve um impulso na origem entao L+ [f (t)] = L [f (t)]. Quando
f (t) nao possui impulso na origem teremos L+ [f (t)] = L [f (t)] = L[f (t)].
Para o caso em que f(t) possui impulso na origem (f (t) possui descontinuidade na
origem) ficamos com:
L+ f(t) = sF (s) f (0+ )
L f(t) = sF (s) f (0 )
Note que na definicao L+ o tempo comeca em t = 0+ e portanto o impulso na origem
fica fora do intervalo considerado, oque nao nos interessa. Assim apenas a definicao L ,
por comecar a contagem dos tempos em t = 0 , nos sera u
til para tratar impulsos na
origem.
Exemplo 5.12 Seja f (t) = et , para t 0. Calcule L[f(t)].
Solucao:
f(t) = (t) et , t 0
s
L[f(t)] = 1
=
s+
s+
s
s
0=
s+
s+
OBSERVAC
OES:
Se a distincao entre L+ e L for necessaria basta substituir t = 0 por t = 0+ ou
t = 0 respectivamente.
Para que L[ n f (t)] exista e preciso que todas as derivadas de f (t) de ordem inferior
`a n existam e sejam transformaveis por Laplace.
Quando todas as condicoes iniciais forem nulas entao:
L [ n f (t)] = sn F (s)
Sistemas Lineares
122
0
s2 +02
podemos obter:
d sen(0 t)
u(t)
dt
0
1
d
(sen(0 t)u(t))
L
0
dt
1
(sF (s) f (0))
0
1
s 0
( 2
0)
0 s + 02
s
2
s + 02
L[cos(0 t)u(t)] = L
=
=
=
=
5.3.7
Quando uma funcao f (t) tende `a um valor constante em regime estacionario, isto e
quando t , este valor constante pode ser diretamente obtido atraves do limite:
lim f (t) = lim sF (s)
s0
onde L[f (t)] = F (s). Note que quando f (t) tende `a um valor constante em regime entao
f(t) tende a zero em regime. Como toda funcao que tende a zero em regime deve possuir
transformada com todos os polos no semi-plano esquerdo conclumos que todos os polos
de L[f(t)] = sF (s) devem estar no semi-plano esquerdo para que o limite acima possa ter
algum sentido. Caso contrario, se algum polo de sF (s) tem parte real nula ou positiva a
funcao f (t) nao tende a um valor constante em regime e portanto a igualdade acima nao
mais se verifica.
Exemplo 5.14 Qual e o valor de regime (se ele existe) da funcao f (t) cuja transformada
1
e F (s) = s(s+1)
?
Solucao: Como os polos de sF (s) n
ao possuem parte real nula nem positiva (os polos
s
ao s = 1) ent
ao f (t) tende `
a um valor constante em regime. E esse valor e dado por:
lim f (t) = lim sF (s) = 1
s0
1
.
s(s+1)
Problema 5.2 Calcule o valor de regime da funcao no tempo cuja transformada e F (s) =
1
. Diga se o teorema do valor final pode ser aplicado e qual e a funcao no tempo.
(s2)
5.3.8
Sistemas Lineares
123
5.3.9
2s+1
.
s2 +s+1
Teorema da Integrac
ao Real
f (t)dt =
F (s)
+
s
f (t)dt
|t=0
s
(5.4)
OBSERVAC
OES:
Se o valor inicial da integral for zero entao:
L
f (t)dt =
F (s)
s
f (t)dt =
f (t)dt
Sendo
f (t)dt|t=0 .
f (t)dt =
0
F (s)
s
Se f (t) possui impulso na origem entao deve-se especificar que a integral comeca
em t = 0 .
5.3.10
Teorema da Diferenciac
ao Complexa
Se f (t) e transformavel por Laplace, entao, exceto nos polos de F (s) vale a seguinte
relacao:
d
L[tf (t)] = F (s).
ds
No caso geral:
n
n
n d
F (s), n = 1, 2, . . . .
L[t f (t)] = (1)
dsn
Sistemas Lineares
5.3.11
124
Integral de Convoluc
ao
Sejam f1 (t) e f2 (t) duas funcoes nulas para t < 0. A Convolucao dessas duas funcoes
f1 (t) e f2 (t) sera representada pela notacao f1 (t) f2 (t) e e definida pela integral:
t
f1 (t) f2 (t) =
f1 (t )f2 ( )d
Propriedades:
f1 (t) f2 (t) = f2 (t) f1 (t)
f1 (t) (f2 (t) + f3 (t)) = f1 (t) f2 (t) + f1 (t) f3 (t)
L[f1 (t) f2 (t)] = L[f1 (t)]L[f2 (t)]
Au
ltima propriedade e muito importante e mostra que fazer a convolucao no tempo e
fazer o produto das transformadas na frequencia.
Prova:
L[f1 (t) f2 (t)] =
f1 (t )f2 ( )d
est dt
=
0
s(t )
f1 (t )e
dt =
f1 ()e
d =
0
Note ainda que L[f1 (t)] e uma funcao complexa da variavel s e nao depende de . Logo
obtemos:
L[f1 (t) f2 (t)] =
= L[f1 (t)]
f2 ( )es d
Sistemas Lineares
125
f (t)
f(t)
...
a
0
...
a(t 2)
t
0
a(t 1)
f (t) = au(t)
a(t n)
n=1
n=1
L[(t n)]
1
sF (s) = a a
ens L[(t)]
s
n=1
ens
a
F (s) = 2 a
s
s
n=1
f(t)
1
a2
0
-
2a
1
a2
Sistemas Lineares
126
Exemplo 5.17 Calcule a transformada de Laplace da funcao x(t) que resolve a seguinte
equacao diferencial a
x + bx + cx = 0, onde x = dx(t)
e x(0) = k1 , x(0)
= k2 .
dt
Solucao: Seja X(s) = L[x(t)]. Tomando a transformada dos dois lados da equacao
temos:
L[a
x + bx + cx] = L[0] = 0
aL[
x] + bL[x]
+ cL[x] = 0
L[x] = X(s)
L[x]
= sX(s) x(0)
L[
x] = s2 X(s) sx(0) x(0)
Exemplo 5.18 Calcule a transformada de Laplace do sinal f (t) = sen(0 t+)u(t), onde
e 0 s
ao constantes.
Solucao: Existem varias formas de se resolver o problema. A seguir apresenta-se uma
forma que explora as propriedades de funcoes senoidais e da funcao impulso.
f (t) = sen(0 t + )u(t)
s sen() + 0 cos()
s2 + 02
Sistemas Lineares
5.4
127
Transformada Inversa
F (s)
Re[s] > c
f (t)
t0
Tranf. Inversa
5.4.1
Frac
oes parciais para p
olos distintos
k(s + z1 )(s + z2 ) . . . (s + zm )
,
(s + p1 )(s + p2 ) . . . (s + pn )
n>m
a1
a2
an
+
+ +
s + p1 s + p2
s + pn
(5.5)
onde ai sao constantes conhecidas como resduos dos polos pi , respectivamente, e sao
calculados da seguinte forma:
ai = (s + pi )F (s)|s=pi
(5.6)
Isto pode ser facilmente verificado. Veja no caso do resduo do polo s = p1 . Multiplicando (5.5) por s + p1 temos:
(s + p1 )F (s) = a1 +
an
a2
(s + p1 ) + +
(s + p1 )
s + p2
s + pn
Sistemas Lineares
128
f (t) = L1 [F (s)] = L1
Note que a expansao por fracoes parciais (5.5) e valida para polos reais e complexos nao
repetidos. Para polos reais os resduos (5.6) sao reais e para polos complexos os resduos
sao complexos.
Exemplo 5.19 (P
olos Reais) Calcule a funcao no tempo cuja transformada e
F (s) =
s+3
(s + 1)(s + 2)
F (s) =
a1
a2
Assim,
t0
Exemplo 5.20 (N
ao Causal) Calcule a transformada inversa da funcao
s3 + 5s2 + 9s + 7
G(s) =
(s + 1)(s + 2)
Solucao: Como o grau do numerador e maior que o grau do denominador devemos
dividir um pelo outro ate que o resto da divisao seja uma funcao com grau do numerador
menor que o grau do denominador, como indicado a seguir.
s+3
(s + 1)(s + 2)
1
2
= s+2+
s+1 s+2
G(s) = s + 2 +
Sistemas Lineares
129
Logo:
g(t) = L1 [G(s)] = L1 [s] + L1 [2] + L1
=
d
(t) + 2(t) + 2et e2t ,
dt
1
2
+ L1
s+1
s+2
t 0
Exemplo 5.21 (P
olos Complexos) Calcule a transformada inversa da funcao
F (s) =
2s + 12
+ 2s + 5
s2
0
(s )2 + 02
L[et cos0 t] =
s
(s )2 + 02
2s + 12
0
s
=A
+B
2
2
2
2
(s + 1) + 2
(s ) + 0
(s )2 + 02
L1 [F (s)] = 5L1
5.4.2
Frac
oes Parciais para p
olos repetidos
Os metodos da secao anterior sao validos para polos distintos. Nesta secao estudaremos
o caso de polos repetidos baseado num exemplo que pode ser facilmente generalizado.
Sistemas Lineares
130
s2 + 2s + 3
(s + 1)3
d
[(s + 1)3 F (s)]s=1
ds
d2
[(s + 1)3 F (s)] = 2b1
ds2
1 d2
[(s + 1)3 F (s)]s=1
2! ds2
M
etodo 2 Atribuindo-se valores para s na igualdade
s = 0 3 = b3 + b2 + b1
s = 1 2 = b3
s = 1 6 = b3 + 2b2 + 4b1
Os dois metodos acima levam aos mesmos valores dos resduos: b3 = 2, b2 = 0, b1 = 1
e portanto:
2
0
1
+ L1
+ L1
3
2
(s + 1)
(s + 1)
s+1
2 t
t
= t e +0+e ,
t 0.
L1 [F (s)] = L1
Sistemas Lineares
5.4.3
131
Frac
oes Parciais para casos especiais
Exemplo 5.23 (P
olos distintos e repetidos) Calcule a transformada inversa da funcao
F (s) =
s2 + 2s + 3
(s + 1)2 (s + 2)
b2
b1
b0
+
+
2
(s + 1)
(s + 1) (s + 2)
Exemplo 5.24 (P
olos reais e complexos) Calcule a transformada inversa da funcao
F (s) =
2s + 12
(s2 + 2s + 5)(s + 1)
Solucao: A funcao possui dois polos complexos e um real. Para utilizarmos os resultados das secoes anteriores devemos primeiro separar os polos complexos dos reais da
seguinte forma:
F (s) =
b1 s + b0
b2
+
+ 2s + 5) (s + 1)
(s2
Sistemas Lineares
5.5
132
Rela
c
ao entre sinais e aloca
c
ao dos p
olos
Transformada F (s)
F (s) =
1
s
F (s) =
1
s2
f (t) = t.u(t)
F (s) =
1
s+a
s2
0
+ 02
Polos de F (s)
Sistemas Lineares
Funcao f (t)
f (t) = cos 0 t.u(t)
133
Transformada F (s)
F (s) =
s2
s
+ 02
F (s) =
0
(s + a)2 + 02
F (s) =
s+a
(s + a)2 + 02
Polos de F (s)
Sistemas Lineares
5.6
134
T
ecnicas de An
alise de Circuitos com Transformada de Laplace
Domnio s
v(t) = R.i(t)
V (s) = R.I(s)
Capacitor
Resistor
Elem.
v(t) =
1
C
t
0
i(t).dt + v(0)
1
.I(s)
s.C
+ 1s .v(0)
Indutor
i(t) = C. dv(t)
dt
V (s) =
i(t) =
1
L
t
0
v(t).dt + i(0)
v(t) = L. di(t)
dt
I(s) =
1
.V
s.L
(s) + 1s .i(0)
Sistemas Lineares
135
Exemplo 5.25 Considere o circuito da Figura 5.11, onde if (t) = 3.et .u(t) mA, R =
10k e C = 25F . Determine o circuito no domnio s e a resposta temporal.
3
s+1
1
1
C
.If (s)
If (s) =
1
s.C
s + RC
Vs (s) =
0, 003
s+1
120
(s + 4).(s + 1)
40
40
s+1 s+4
Sistemas Lineares
136
Sistemas Lineares
5.7
137
f (t)2 dt
(5.7)
Esta definicao de energia e uma generalizacao do conceito de energia dissipada em resistores. Por exemplo, se f (t) representa a tensao ou corrente num resistor unitario, a
energia dissipada no resistor e dada pela integral acima. Os sinais que possuem energia
limitada (E < ) sao portanto de grande interesse pratico.
Veremos a seguir que um sinal cuja transformada de Laplace e uma funcao racional
que possui todos os polos no semi-plano esquerdo estrito, isto e, polos com parte real
estritamente negativa, e um sinal de energia limitada. Para esses sinais a integral acima
existe e e finita. Seja o seguinte sinal:
x(t) = 1 + 2 e2t + 3 et + k1 et sen0 t + k2 et cos0 t,
t0
A transformada de x(t) e:
X(s) = L[x(t)] =
1
2
3
k1 0
k2 (s + 1)
+
+
+
+
s
s + 2 s + 1 (s + 1)2 + 02 (s + 1)2 + 02
Note que todos os polos possuem parte real negativa, exceto o polo na origem. Assim, os
polos reais de X(s) tornam-se expoentes de funcoes exponenciais decrescentes no tempo.
Os polos complexos estao associados `a sinais que causam oscilacoes amortecidas. O
amortecimento dessas oscilacoes e definido pela parte real dos polos (Re[polos] = 1 no
caso) e a frequencia de oscilacao e definida pela parte imaginaria do polo (Im[polo] =
0 ). O efeito temporal dos polos com parte real negativa diminui exponencialmente e
desaparece completamente em regime permanente, isto e, quando t .
Um sinal x(t) cuja transformada seja analtica no semi-plano direito 1 mas tenha um
polo simples na origem vai ter um nvel DC igual ao resduo desse polo (1 no caso acima).
O valor do sinal x(t) acima em regime permanente (t ) e constante e igual `a . Note
que nesse caso o sinal nao tem energia limitada pois a integral acima vai divergir dado
que o sinal nao converge para zero em regime.
Assim, um sinal qualquer x(t) vai ter um valor zero em regime (converge para zero
quando t ) apenas quando todos os polos da transformada possuem parte real
negativa. Se a transformada possui um polo na origem ( e os demais no semi-plano
esquerdo estrito) o sinal sera constante com um nvel DC nao nulo em regime. Em
qualquer outra situacao o sinal e divergente, isto e, nao tera um valor de regime finito.
A energia do sinal sera limitada apenas no primeiro caso, isto e, quando o sinal converge
para zero quando t .
1
Sistemas Lineares
5.8
138
Resolu
c
ao de Equac
oes Diferenciais
Atraves das leis da fsica podemos obter um modelo de comportamento para todos os
sistemas. Para sistemas dinamicos esse modelo e uma equacao diferencial. Este e o caso
por exemplo de motores, circuitos, turbinas e todos os outros dispositivos estudados na
engenharia. Saber como o sistema se comporta para dadas condicoes iniciais e uma dada
excitacao e equivalente a saber resolver a equacao diferencial.
A Transformada de Laplace pode ser utilizada para resolver equacoes diferenciais
lineares invariantes no tempo. Para isso basta transformar por Laplace cada um dos
termos da equacao dif. obtendo assim a transformada da funcao que resolve a equacao.
Em seguinda, utiliza-se a transformada inversa para encontrar a solucao no domnio do
tempo.
Exemplo 5.26 Resolva a seguinte equacao diferencial x +2x+5x
= b s
ao constantes dadas e g(t)=0.
Solucao: Note que
L[x] = X(s)
L[x]
= sX(s) x(0)
L[
x] = s2 X(s) sx(0) x(0)
1
s+2
x(0) + 2
x(o)
+ 2s + 5
s + 2s + 5
De forma similar ao exemplo 5.21 temos:
X(s) =
X(s) =
s2
s2
s+1
1
1
x(0) + 2
x(0) + 2
x(o)
+ 2s + 5
s + 2s + 5
s + 2s + 5
e consequentemente
x(t) = L1 [X(s)] = [et cos(2t) + 0.5et sen(2t)]x(0) + 0.5et sen(2t)x(0)
= 0. Encontre
a resposta desse sistema quando o mesmo e excitado por um degrau de amplitude 3, isto
e, g(t) = 3u(t).
Solucao: Note que
3
s
L[x] = X(s)
L[x]
= sX(s) x(0) = sX(s)
L[
x] = s2 X(s) sx(0) x(0)
= s2 X(s)
L[3u(t)] =
Sistemas Lineares
Logo:
139
3
s
3
3
3 s+2
X(s) =
=
2
2
s(s + 2s + 5)
5s 5 s + 2s + 5
s2 X(s) + 2sX(s) + 5X(s) =
2
2
5s 5 (s + 1) + 2
5 (s + 1)2 + 22
X(s) =
Logo:
3
1
s+1
3
1 3
L
L1
5s
5 (s + 1)2 + 22
5 (s + 1)2 + 22
3
3
3
=
et sen2t et cos2t, t 0.
5 10
5
L1 [X(s)] = x(t) = L1
=
0, x(0) = 0 ; (d) g(t) = 3u(t), x(0)
= 1, x(0) = 1
x(0)
g(t)
+
-
1
s
x(0)
x
1
s
2
5
Figura 5.14: Diagrama de simulacao analogica
5.9
A resposta de todo sistema linear invariante no tempo pode ser decomposta em duas
parcelas: uma que depende do sistema e do sinal de entrada e outra que depende do
sistema e das condicoes iniciais. A primeira parcela chamaremos de Resposta de Estado
Zero ja que esta parcela indica como um sistema, inicialmente em repouso (condicoes
iniciais nulas), responde a um dado sinal de entrada. A segunda parcela chamaremos de
Resposta de Entrada Zero pois ela indica como um sistema se comporta quando e deixado
para responder livremente `as suas condicoes inicias (sem excitacao externa). .
Sistemas Lineares
1.5
1.3
1.1
0.9
(a) 0.7
0.5
0.3
0.1
-0.1
-0.3
-0.5
1.5
1.3
1.1
0.9
0.7
(b) 0.5
0.3
0.1
-0.1
-0.3
-0.5
1.5
1.3
1.1
0.9
0.7
(c) 0.5
0.3
0.1
-0.1
-0.3
-0.5
1.5
1.3
1.1
0.9
0.7
(d) 0.5
0.3
0.1
-0.1
-0.3
-0.5
140
x(t)
0.0
0.7
1.4
2.1
2.8
3.5
4.2
4.9
5.6
6.3
7.0
1.4
2.1
2.8
3.5
4.2
4.9
5.6
6.3
7.0
1.4
2.1
2.8
3.5
4.2
4.9
5.6
6.3
7.0
2.1
2.8
3.5
4.2
4.9
5.6
6.3
7.0
x(t)
0.0
0.7
x(t)
0.0
0.7
x(t)
0.0
0.7
1.4
Sistemas Lineares
141
As respostas de Estado Zero e Entrada Zero de um sistema descrito por (5.11) podem
ser determinadas atraves da Transformada de Laplace .
Exemplo 5.28 Encontre as respostas de Estado Zero e Entrada Zero do circuito RLC
serie da figura 5.17.
Solucao: Do exemplo 5.29 temos que o comportamento dinamico entrada/sada do
circuito e dado por (5.12). Tomando a transformada dos dois lados da equacao temos:
L[a2 y + a1 y + a0 y] = L[b0 x]
(5.8)
Portanto:
a2 s + a1
a2
b0
X(s) +
y(0) +
y(0)
2
2
a2 + a1 s + a0
a2 s + a1 s + a0
a2 s + a1 s + a0
Y (s) = F (s)X(s) + F0 (s)y(0) + F1 (s)y(0)
(5.9)
Y (s) =
s2
onde
F (s) =
a2
s2
b0
a2 s + a1
a2
, F0 (s) =
, F1 (s) =
2
2
+ a1 s + a0
a2 s + a1 s + a0
a2 s + a1 s + a0
Considerando f (t) = L1 [F (s)], f0 (t) = L1 [F0 (s)] e f1 (t) = L1 [F1 (s)] podemos ent
ao
reescrever a express
ao acima com o auxlio da anti-transformada na forma:
1
y(t) = L1 [Y (s)] = L1 [F (s)X(s)] + y(0)L1 [F0 (s)] + y(0)L
[F1 (s)]
= f (t) x(t) + y(0)f0 (t) + y(0)f
1 (t)
x(t)
F (s)
y(0)
y(0)
y(t)
y(0)
F0 (s)
y(t)
x(t)
F (s)
y(0)
F1 (s)
(5.10)
Sistemas Lineares
142
nio do tempo.
1 (t) no dom
Podemos agora generalizar os resultados acima para sistemas de ordem mais elevada.
Considere um sistema descrito pela seguinte equacao diferencial:
an n y(t) + + a1 y(t) + a0 y(t) = bm m x(t) + + b1 x(t) + b0 x(t)
(5.11)
V(t)
-
Vc (t)
-
Exemplo 5.29 Considere o circuito RLC serie descrito na figura 5.17. A entrada do
sistema e a tensao V (t) e a sada e a tensao no capacitor Vc (t). Em termos da notacao
acima temos x(t) = V (t) e y(t) = Vc (t) e o comportamento dinamico entrada/sada e
regido pela seguinte equacao diferencial:
a2 y + a1 y + a0 y = b0 x
(5.12)
= V c (0).
Se ao inves do sistema de segunda ordem do exemplo acima, considerarmos um sistema
de ordem generica como em (5.11) obteramos:
n1
fi (t)ci
i=0
(5.13)
Sistemas Lineares
onde ci =
di y(t)
|
dti t=0
143
sao as condicoes iniciais.
(5.14)
3. A sada de um sistema depende das condicoes iniciais do mesmo. Este fato pode ser
verificado em (5.13) pela presenca das constantes ci que sao as condicoes iniciais.
Esta parcela da resposta recebe o nome de resposta de entrada zero do sistema.
Esta e a resposta do sistema quando a entrada e nula.
n1
yenz (t) =
n1
Fi (s)ci
(5.15)
i=0
5.10
Func
ao de Transfer
encia e Estabilidade
Veremos a seguir que a resposta de Estado Zero de um sistema esta associada `a duas
nocoes muito importantes: funcao de transferencia e estabilidade.
Sistemas Lineares
144
Defini
c
ao 5.4 (Fun
c
ao de Transfer
encia) Funcao de transferencia e uma funcao complexa que representa a relacao sada/entrada do sistema para condicoes iniciais nulas.
Pela definicao acima nota-se que a nocao de funcao de transferencia esta relacionada
com a resposta de Estado Zero do sistema. A relacao complexa sada/entrada de um
sistema com condicoes iniciais nulas pode ser obtida diretamente da resposta de Estado
Zero (5.14): Y (s)/X(s) = F (s). Assim, um sistema que possua a resposta de Estado
Zero (5.14) tera F (s) como funcao de transferencia. Quando se conhece a funcao de
transferencia F (s) de um sistema e a transformada do sinal de entrada X(s) se conhece
importante
tambem a resposta de Estado Zero do mesmo que e dada por (5.14). E
notar que a funcao de transferencia depende apenas dos parametros fsicos do sistema e
da estrutura entrada/sada do mesmo. Veja o exemplo 5.28. A entrada e as condicoes
inicias nao afetam a funcao de transferencia.
Quando as condicoes iniciais sao nulas resposta total do sistema e a propria resposta
de Estado Zero do mesmo, como pode ser visto nas equacoes (5.9) e (5.10).
Domnio do Tempo: y(t) = yesz (t) = f (t) x(t)
Domnio da Frequencia: Y (s) = Yesz (s) = F (s)X(s)
A funcao f (t) = L1 [F (s)] recebe o nome de Resposta Impulsional pois f (t) e a resposta do sistema quando as condicoes iniciais sao nulas e a entrada x(t) e um impulso no
instante t = 0 (X(s) = 1).
Defini
c
ao 5.5 (Sistemas Causais ou N
ao-Antecipativos) Um sistema dinamico e
dito ser Causal ou Nao-Antecipativo se sua Resposta Impulsional e nula para t < 0.
Pela definicao acima nota-se que a resposta y(t) de um sistema causal excitado com um
sinal x(t), apresenta a seguinte propriedade: o valor de y(t)|t=tf so depende da entrada
x(t) e da resposta impulsional f (t) para valores de tempo t tf . Em outras palavras,
a dinamica de um sistema causal em qualquer instante de tempo t = tf depende (nao
depende) da entrada e da resposta impulsional para valores de tempo menores (maiores)
que tf . Essa propriedade e mostrada a seguir.
y(t) = f (t) x(t) =
f (t )x( )d
para t = tf temos f (tf ) = 0 para > tf . Logo f (tf )x( ) = 0 para > tf e
portanto:
tf
y(tf ) =
0
f (tf )x( )f
Sistemas Lineares
145
Defini
c
ao 5.6 (Estabilidade de Sistemas) Um sistema e dito ser estavel se todos os
polos da sua funcao de transferencia est
ao localizados no semi-plano esquerdo estrito,
isto e, Re[polos] < 0. Caso contr
ario o sistema e dito ser instavel.
Pela definicao acima nota-se que a estabilidade e uma propriedade intrnseca do sistema.
Ela so depende da sua funcao de transferencia e portanto dos seus parametros fsicos e
da estrutura entrada/sada.
5.11
Diagrama de Blocos
Sistemas Lineares
146
x(t)
y(t)
I(t)
-
Entrada
Sada
SISTEMA
LCs2
1
X(s)
+ RCs + 1
Portanto:
F (s) =
LCs2
(5.16)
1
+ RCs + 1
X(s)
Y(s)
F(s)
Sistemas Lineares
147
entrada x(t) e dada por y(t) = x(t) f (t) onde f (t) = L1 [F (s)] e a resposta impulsional
do circuito.
Note que no diagrama de blocos acima foram eliminadas as informacoes sobre todas as
outras vari
aveis do circuito (corrente, etc). A Funcao de Transferencia da informacao
apenas sobre a relacao de causa-efeito entre as vari
aveis de entrada e de sada.
possvel, no entanto, explicitar a dependencia de outras vari
E
aveis no diagrama de
blocos atraves de simples manipulacao de equacoes. Por exemplo, para fazer aparecer a
vari
avel corrente no diagrama de blocos do circuito temos:
X(s) = RI(s) + LsI(s) + Y (s)
CsY (s) = I(s)
X(s) Y (s) = (R + Ls)I(s) I(s) =
1
Y (s) = Cs
I(s)
1
(X(s)
R+Ls
Y (s))
Agora essas equacoes podem ser transformadas em diagramas como mostra a figura
5.20.
X(s)
+
-
1
R+Ls
I(s)
1
Cs
Y(s)
Y(s)
5.12
Sistemas Realimentados
A presenca de uma malha fechada num diagrama de blocos caracteriza o que se chama
de sistema realimentado. De uma maneira geral um sistema realimentado pode ser caracterizado pelo diagrama da figura 5.21 onde
X(s) e a transformada do sinal de entrada.
Y (s) e a transformada do sinal de sada.
G(s) funcao de transferencia do sistema a ser controlado, incluindo acionadores, medidores e controladores (Filtros para fins de controle).
Sistemas Lineares
148
X(s)
E(s)
+
Y(s)
G(s)
H(s)
H(s) funcao de transferencia de realimentacao que inclui transdutores e eventuais controladores adicionais.
A Funcao de Transferencia entre X(s) e Y (s) no diagrama acima e conhecida como
F.T. de malha fechada e pode ser obtida atraves das equacoes inicadas no diagrama.
E(s) = X(s) H(s)Y (s)
Y (s) = G(s)E(s)
Para se obter a funcao de transferencia entre X(s) e Y (s) deve-se eliminar todas as
variaveis intermediarias, E(s) no caso acima. Com isso temos a seguinte relacao:
Y (s) =
X(s)
G(s)
X(s)
1 + G(s)H(s)
G(s)
1+G(s)H(s)
(5.17)
Y(s)
F.T.M.F.
X(s)
+
-
Y(s)
(R+Ls)(Cs)
Exemplo 5.32 Vimos que a F.T. entre X(s) e Y (s) no circuito da figura 5.19 e F (s) =
1/(LCs2 + RCs + 1). Vimos tambem que ao fazer aparecer a corrente no diagrama
Sistemas Lineares
149
F (s) =
5.12.1
1
G(s)
=
2
1 + G(s)
LCs + RCs + 1
Estabilidade de Conex
oes
X(s)
Y(s)
+
G2 (s)
G1 (s)
Y(s)
G2 (s)
Sistemas Lineares
150
Como as razes de D1 (s) e de D2 (s) possuem parte real negativa entao as razes de
D1 (s)D2 (s) possuem as mesma caractersticas. Logo a funcao de transferencia de X(s)
para Y (s) na coneccao da Figura 5.24 e estavel.
Ja no caso da conexao da Figura 5.25 temos:
Y (s) =
=
G1 (s)
X(s)
1 + G1 (s)G2 (s)
1+
N1 (s)
D1 (s)
N1 (s) N2 (s)
D1 (s) D2 (s)
N1 (s)D2 (s)
X(s)
D1 (s)D2 (s) + N1 (s)N2 (s)
Agora as razes do polinomio D1 (s)D2 (s) + N1 (s)N2 (s) podem ter parte real positiva
mesmo se as razes de D1 (s) e D2 (s) possuem parte real negativa. Esse e o caso, por
exemplo, se N1 (s) = 2, N2 (s) = 1 e D1 (s) = D2 (s) = s + 1.
5.12.2
Referencia
R(s)
+
-
G1 (s)
+
G2 (s)
C(s)
H(s)
G2 (s)
D(s)
1 + G2 (s)H(s)G1 (s)
Quando D(s) = 0 tem-se o diagrama da figura 5.28 e novamente com (5.17) temos:
CR (s) =
G1 (s)G2 (s)
R(s)
1 + G1 (s)G2 (s)H(s)
Sistemas Lineares
151
D(s)
G1 (s)
+
G2 (s)
CD (s)
H(s)
+
-
G2 (s)
CR (s)
H(s)
Logo:
Ctotal (s) = CD (s) + CR (s) =
G1 (s)G2 (s)
G2 (s)
D(s) +
R(s)
1 + G2 (s)H(s)G1 (s)
1 + G1 (s)G2 (s)H(s)
5.13
Problemas complementares
Sistemas Lineares
152
Problema 5.12 Sabendo que a resposta impulsional do sistema da figura 5.29(a) e w(t) =
2exp(t) , t 0 verifique se o sistema realimentado da figura 5.29(b) e est
avel. Justifique sua resposta.
u
G(s)
e
-
G(s)
1
s
10
(a)
(b)
Bibliografia
OGATA, K. Engenharia de Controle Moderno. 2.ed., Prentice Hall do Brasil:Rio de
Janeiro, RJ, 1993, 781p.
PHILLIPS, C. L. Signals, Systems, and Transforms. Prentice-Hall Inc.:New Jersey,
USA, 1995, 707p.
CLOSE, C. M., The Analysis of Linear Circuits, Harcourt, Brace & World, Inc.:New
York, USA, 1966, 716p.
Sistemas Lineares
154
Captulo 6
Resposta ao Degrau
6.1
Introduc
ao
Um grande n
umero de problemas de controle consiste em se manter constante a variavel
de sada. Veja por exemplo o problema de controle de posicionamento de uma antena
indicado na figura 1.5.A entrada do sistema, que representa o valor desejado da variavel
controlada (sada) e neste caso um degrau com amplitude igual ao valor desejado para
a sada. Quando se quer mudar a posicao da antena de uma posicao inicial, digamos
posicao zero, para uma nova posicao, digamos posicao um, o sinal de entrada deve ser
um degrau unitario. Ao se aplicar um degrau na entrada desse sistema de controle, a
posicao da antena vai evoluir da posicao zero para a posicao um segundo uma curva
que depende de como o sistema de controle foi projetado. Curvas tpicas dessa evolucao
podem ser encontradas na figura 6.1. Normalmente deseja-se um transitorio rapido, com
poucas oscilacoes e que a variavel controlada, posicao da antena no caso, va para o valor
desejado sem erro significativo de posicao em regime, isto e, erro de regime despresvel.
Para atender todos esses requisitos de performance, quando isso e possvel, o engenheiro
deve saber projetar adequadamente os filtros de controle do sistema. O primeiro passo,
no entanto, e saber especificar matematicamente os ndices de performance desejados
para a resposta. Veja na figura 6.1 que a resposta (a) e mais oscilatoria que as demais. A
resposta (c) atinge o valor de regime mais rapido que as demais e todas as tres possuem
erro de regime nulo (valor final da resposta e exatamente o valor desejado).
Neste captulo estudaremos alguns ndices de performance da resposta ao degrau que
nos permitira quantificar matematicamente o tamanho das oscilacoes da resposta, a rapidez da resposta e o erro de regime cometido.
Outros sinais de entrada como impulso e funcao rampa (x(t) = t) tambem sao de
interesse. No entanto, para condicoes iniciais nulas, a resposta de um sistema (linear
invariante) ao impulso, degrau, e rampa estao ligadas entre si. Para ilustrar este fato,
seja F (s) a F.T. de um sistema linear invariante indicado na figura 6.2.
f (t) = L1 [F (s)]
Resposta Impulsional: X(s) = 1 Y (s) = F (s) y(t) = f (t)
Sistemas Lineares
156
2.0
1.8
1.6
1.4
1.2
(a) 1.0
0.8
0.6
0.4
0.2 +
0.0
0
12
15
18
21
24
27
30
2.0
1.8
1.6
1.4
1.2
(b) 1.0
0.8
0.6
0.4
0.2 +
0.0
0
12
15
18
21
24
27
30
2.0
1.8
1.6
1.4
1.2
(c) 1.0
0.8
0.6
0.4
0.2 +
0.0
0
12
15
18
21
24
27
30
2.0
1.8
1.6
1.4
(d) 1.2
1.0
0.8
0.6
0.4
0.2 +
0.0
0
12
15
18
21
24
27
30
X(s)
Y(s)
F(s)
Sistemas Lineares
157
1
s
1
s2
t
0
1
s2
y(t) =
t
0
t
0
f (t)dt
f (t)dtdt
Note que, para condicoes iniciais nulas, a resposta ao impulso e a resposta `a rampa sao
respectivamente a derivada e a integral da resposta ao degrau. Por esse motivo vamos
nos concentrar na resposta ao degrau de agora em diante.
6.2
An
alise de Sistemas de Primeira Ordem
Sistemas cuja funcao de transferencia possui apenas um polo sao conhecidos como
sistemas de primeira ordem.
Exemplo 6.1 Verifique que o circuito da figura 6.3 e um sistema de primeira ordem.
Solucao: Para mostrar que o sistema e de primeira ordem precisamos encontrar a
funca
o de transferencia do mesmo e para isso se supoe que o circuito possui condicoes
iniciais nulas. As equacoes que regem o comportamento desse sistema s
ao indicadas
abaixo.
x + RI + y = 0
, condicao inicial nula (y(0) = 0)
I = C y
Aplicando Laplace temos:
R
+
x(t)
y(t)
x + RC y + y = 0
Y (s)
1
1
=
=
X(s)
RCs + 1
Ts + 1
1
Ts+1
Y(s)
Sistemas Lineares
158
1
T s+1
1
1
1
X(s) =
Ts + 1
Ts + 1 s
T
1
y(t) = 1 et/T ,
s Ts + 1
t0
entrada
y(t)
0
sada
t
1)
dy
|
dt t=0
1
T
Sistemas Lineares
6.3
159
An
alise de Sistemas de Segunda Ordem
Sistemas de segunda ordem sao aqueles cuja funcao de transferencia possui dois polos.
Nesta secao vamos estudar um tipo especial de sistemas de segundo ordem conhecido na
literatura de controle como sistema de segunda ordem padrao:
F (s) =
n2
s2 + 2n s + n2
(6.1)
onde e n recebem o nome de taxa de amortecimento e frequencia natural do sistema respectivamente. Os valores desses parametros dependem dos parametros fsicos do sistema
estudado, como ilustra o exemplo a seguir.
x(t)
y(t)
I(t)
-
Entrada
Sada
SISTEMA
LCs2
1
+ RCs + 1
1
n2
=
LCs2 + RCs + 1
s2 + 2n s + n2
Logo:
n2 =
R
R
1
; 2n =
=
LC
L
2
C
L
Sistemas Lineares
160
6.3.1
n2
(s2 + n2 )s
6.3.2
Quando 0 < < 1 os polos da funcao de transferencia indicada em (6.1) sao complexos
e do lado esquerdo do eixo imaginario. Isto pode ser verificado da seguinte forma. Os
polos sao dados pela equacao:
p1,2 =
2n
4 2 n2 4n2
= n n
2
2 1
Sistemas Lineares
161
n2
(s2 + 2n s + n2 )s
et
1 2
sen(d t + ),
= tan1
1 2
(6.2)
6.3.3
Y (s) =
n2
s2 + 2n s + n2
n2
s(s + n )2
6.3.4
n2
(s p1 )(s p2 )s
Sistemas Lineares
162
ep1 t ep2 t
+
y(t) = 1 +
p1
p2
n
2 2 1
Esta resposta pode ser vista na figura 6.1(d) para = 2 e n = 2. Note que se >> 1
entao, para o mesmo valor de n , temos |p1 | >> |p2 | e portanto o efeito do polo p1 sobre
a resposta desaparece bem mais rapido que o efeito do polo p2 que esta mais proximo do
eixo imaginario. Sendo assim para valores de >> 1 o sistema se torna extremamente
lento. Um sistema de primeira ordem com um polo p2 teria uma resposta muito parecida.
6.3.5
Caso inst
avel ( < 0)
6.4
Indices de desempenho
Nesta secao estudaremos formas de classificar quao boas sao as respostas da figura
6.1. Como a resposta transitoria `a um degrau normalmente apresenta oscilacoes antes
de atingir o regime permanente, torna-se imperativo a criacao de ndices de desempenho
que permitam quantificar tamanho de oscilacoes, tempo de duracao do transitorio, etc.
Sao comuns os seguintes ndices:
o tempo necessario para a resposta atingir o seu valor maximo.
tp (instante de pico) E
o valor relativo da diferenca entre o valor maximo da
Mp (sobressinal m
aximo) E
sada (ao longo do tempo) e o valor da sada em regime.
Mp =
y(tp ) y()
y()
ts (tempo de acomoda
c
ao) Tempo necessario para confinar a resposta numa faixa
em torno do seu valor de regime. Esta faixa caracteriza a tolerancia de erro, que
tipicamente vale 2 ou 5% do valor de regime).
Sistemas Lineares
163
2.0
1.8
1.6
Mp
faixa de erro toler
avel em regime
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 +
0
tp
ts
3
12
15
18
21
24
27
30
A figura 6.7 ilustra os ndices de desempenho descritos acima. Existem outros ndices de
performance que nao foram indicados acima e podem ser encontrados em qualquer livro
de controle de sistemas, por exemplo Ogata (1993).
Em geral nao e possvel se determinar expressoes analcas para os ndices de desempenho da resposta ao degrau indicados acima. No entanto, para sistemas de segunda ordem
do tipo (6.1) subamortecidos (1 > > 0) isto e possvel e essas expressoes sao obtidas a
seguir.
Instante de Pico (tp ):
O instante de pico pode ser caracterizado como sendo
o primeiro instante de tempo (exceto a origem) para o qual a derivada temporal da
resposta e nula. Tomando a derivada temporal da resposta y(t) em (6.2) e igualando `a
zero encontramos:
tp =
=
(6.3)
d
n 1 2
Sobressinal M
aximo (Mp ): Note que num sistema do tipo (6.1) o valor de regime da
resposta ao degrau unitario e um. Como y() = 1 vamos utilizar (6.2) para obter:
Mp =
y(tp ) y()
= y(tp ) 1 = e d = e 12
y()
(6.4)
Note que Mp depende somente de e quando 1 nao existe oscilacao e Mp nao tem
mais sentido.
Tempo de Acomoda
c
ao (ts ): Diferentemente do Sobressinal e do intante de pico,
nao existe uma expressao analtica exata para o tempo de acomodacao ts . Existem
abacos que permitem a determinacao exata de ts . Veja por exemplo Ogata (1993). A
seguir apresentamos duas possibilidades para se obter uma aproximacao de ts .
Sistemas Lineares
164
et
1 2
sen(d t + ) ,
= tan1
1 2
Impondo que a amplitude do seno esteja dentro da faixa de tolerancia que caracteriza
o tempo de acomodacao temos uma condicao suficiente para garantir que o tempo de
acomodacao foi atingido com a dada tolerancia. Note que o valor de regime da resposta
e y() = 1 e a amplitude do seno tende `a zero quando t . Seja a tolerancia de
erro que define o tempo de acomodacao. Impondo que a amplitude do seno esteja dentro
dessa tolerancia temos:
y(ts ) y()
=
y()
ets
1 2
sen(d + )
ts =
ln(
1 2)
(6.5)
(6.6)
1 2 = 4 e a parte imagin
aria dos polos.
ln(0,02 1 2 )
= 1, 38seg e o tempo de acomodacao com 2% de erro.
ts (2%) =
3
ln(0,05 1 2 )
ts (5%) =
= 1, 07seg e o tempo de acomodacao com 5% de erro.
3
Sistemas Lineares
165
1.20
1.08
0.96
0.84
0.72
0.60
0.48
0.36
0.24
0.12
0.00 +
0.0
0.3
0.6
0.9
1.2
1.5
1.8
2.1
2.4
2.7
3.0
6.5
tens
ao de erro e(t)
amplificador de potencia
posic
ao da antena
antena
torque do eixo da antena
engrenagens
torque do eixo do motor
tens
ao do motor Ea
motor
Sistemas Lineares
166
k0
k0
c(t)
e(t)
Amplificador de pot
encia: Esse dispositivo tem como funcao suprir com energia
o sistema de controle. Note que o sinal de entrada do amplificador e(t) e um sinal de
erro oriundo de medidores e portanto nao possui energia suficiente para acionar o motor.
Vamos considerar que o amplificador e ideal e possui um ganho de tensao k1 . Assim o sinal
de sada do amplificador Ea (t) e dado por Ea (t) = k1 e(t). Incorporando o amplificador
no diagrama de blocos 6.10 obtemos um novo diagrama indicado na figura 6.11.
r(t)
k0
k0
c(t)
e(t)
k1
Ea (t)
Motor DC: A funcao do motor DC e acionar a antena para que ela esteja sempre
apontada para a direcao desejada. Sao comums as palavras acionador e servomotor para
designar a funcao do motor nesse tipo de sistema de controle.
O servomotor pode ser operado de dois modos. Num modo a corrente de campo
(estator) e mantida constante e uma tensao ajustavel e aplicada `a armadura (rotor) e
no outro modo se faz o contrario. Esses modos de operacao possuem caractersticas
diferentes e apenas o primeiro sera considerado aqui.
Sistemas Lineares
167
(6.8)
circuito de campo
Vcc
La
+
+
Ia
Ea (t)
Tm
(6.9)
(6.10)
Podemos agora incluir o motor no diagrama de blocos da figura 6.11 para obter o diagrama
da figura 6.13.
Engrenagens: O sistema de engrenagens tem como funcao adequar a velocidade de
rotacao do eixo da antena ao eixo do rotor. Um sistema de engrenagens possui funcao
analoga do transformador em sistemas eletricos. Nos dois casos, a potencia do primario
deve ser igual `a do secundario: no caso do transformador a potencia e o produto da tensao
pela corrente V1 I1 = V2 I2 e no caso da engrenagem a potencia e o produto do torque
Sistemas Lineares
168
r(t)
k0
c(t)
k0
e(t)
k1
Ea (t)
La
T
k2 m
Ra
T
k2 m
+ k3 m = Ea
Tm , m
c(t)
k0
e(t)
k1
Ea (t)
La
T
k2 m
Ra
T
k2 m
+ k3 m = Ea
Tm , m
Tc , c
Jc
Tc c
bc
Sistemas Lineares
169
Torques = 0
Tc = Jc c + bc c
(6.11)
e com a expressao acima podemos incluir a antena no diagrama 6.14 para obter o diagrama
da figura 6.16. Note que as variaveis m , Tm do eixo do motor e as variaveis c , Tc do eixo
r(t)
k0
c(t)
k0
dt
c
e(t)
Tc = Jc c + bc c
k1
Tc , c
Ea (t)
La
T
k2 m
Ra
T
k2 m
+ k3 m = Ea
Tm , m
k2 k3
k2
Ea
m
Ra
Ra
k2 k3
k2
)c =
Ea
2
n Ra
n Ra
(6.12)
Jc c(t) + (bc +
k2 k3
k2
)c(t)
=
Ea
2
n Ra
n Ra
(6.13)
e como c = c(t)
temos
2
C(s)
n Ra
=
Ea (s)
Jc s2 + (bc +
k2 k3
)s
n2 Ra
(6.14)
Sistemas Lineares
r(t)
170
c(t)
k0
k2
n Ra
k k
Jc s2 +(bc + n22R3 ) s
e(t)
k1
Ea (t)
Com isto o diagrama 6.16 pode ser simplificado como indicado no diagrama 6.17. Por
conveniencia de notacao iremos definir a funcao G(s) indicada a seguir.
G(s) =
K
J s2 +B s
(6.15)
K=
K0 K1 K2
n Ra
, B = bc +
K2 K3
n2 Ra
, J = Jc
Com G(s) acima o diagrama 6.17 pode ser rescrito como indicado na figura 6.18 que e
uma forma mais conveniente para nossos propositos. Agora a Funcao de Transferencia
r(t)
+
c(t)
K
Js2 +B
C(s)
G(s)
K
=
=
2
R(s)
1 + G(s)
Js +Bs+K
B
K
, 2n =
J
J
(6.16)
Sistemas Lineares
171
podemos entao corrig-lo ajustando os parametros fsicos do sistema, tais como o ganho
do amplificador k1 , ou o ganho do potenciometro k0 . Esse ajuste deve ser tal que o novo
valor da taxa de amortecimento seja compatvel com as oscilacoes admissveis para o
sistema. Lembre que quanto menor o valor de maior as oscilacoes da resposta ao degrau.
Exemplo 6.5 Suponha que o sistema de controle da figura 6.18 tenha um momento de
inercia J = 1 e um coeficiente de atrito viscoso B = 1. Determine o valor do ganho K
de tal forma que o sobressinal Mp na resposta ao degrau seja de 20% . Verifique o tempo
de acomodacao obtido.
Solucao: Com os valores dados a funcao de transferencia de malha fechada e:
G(s)
K
C(s)
=
= 2
R(s)
1 + G(s)
s +s+K
Comparando com (6.1) obtemos os valores de , n seguintes:
n2 = K , 2n = 1
Para que o sobressinal seja de 20% devemos ter a seguinte igualdade satisfeita:
Mp = 0, 2
1 2
= ln(0, 2)
= 0, 456
Sistemas Lineares
172
2.0
1.8
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 +
0
12
15
18
21
24
27
30
r(t)
k0
c(t)
k0
dt
tac
ometro
e(t)
k4
vT (t)
Tc = Jc c + bc c
k1
Tc , c
Ea (t)
La
T
k2 m
Ra
T
k2 m
+ k3 m = Ea
Tm , m
R(s)
+
K
J s+B
c (s)
1
s
C(s)
Sistemas Lineares
173
6.18. Isto nos leva ao diagrama da figura 6.21. Agora, com os par
ametros K e Q para
serem ajustados temos mais graus de liberdade para ajustarmos o sobressinal e o tempo
de acomodacao simultaneamente. Para ver como fazer isso vamos deduzir da figura 6.21
a nova Funcao de Transferencia do sistema de controle indicada a seguir.
C(s)
K
=
R(s)
J s2 + (B + KQ)s + K
Comparando a funcao de transferencia acima com a forma padrao da equacao (6.1) podemos encontrar os novos valores da taxa de amortecimeto e da frequencia natural n
do novo sistema de controle.
2n =
B + KQ
J
n2 =
K
J
(6.17)
1 2
= ln(0, 2)
= 0, 456
Com o valor de = 0, 456 podemos agora encontrar o valor de n impondo que o tempo
de acomodacao dado por (6.5) seja de 20% :
ts (5%) =
ln(0, 05 1 2 )
= 2 n = 3, 41
n
2n = 1 + KQ
Q=
2 n 1
K
= 0, 1814
Sistemas Lineares
174
2.0
1.8
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 +
0
10
6.6
Problemas complementares
r(t)
10
kp
-
sistema
1
s
(t)
kv
kp
Sistemas Lineares
175
ao degrau unit
ario desse sistema e indicada na figura 6.25. Encontre os valores de k1 e
k2 sabendo que o motor e regido pela equacao diferencial + 10 = Va .
r
Va
-
1
s
motor
k1
k2
2.0
1.8
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 +
0
10
12
14
16
18
20
Sistemas Lineares
176
Bibliografia
DAZZO, J. J.; HOUPIS, C. H., An
alise e Projeto de Sistemas de Controle Lineares,
Editora Guanabara:Rio de Janeiro, RJ, 1984, 660p.
K. OGATA, Engenharia de Controle Moderno, Prentice Hall 1993.
K. OGATA, Discrete Time Control Systems, Prentice Hall, 1995.
Sistemas Lineares
178
Captulo 7
Sistemas Discretos e Amostrados
7.1
Introduc
ao
Os termos tempo contnuo e analogico sao identicos quando empregados para caracterizar sinais e sistemas. Sinais analogicos sao funcoes de uma variavel de tempo contnuo
e sistemas analogicos sao aqueles que manipulam sinais analogicos. De maneira analoga,
os termos tempo discreto e digital sao tambem identicos. Um sinal de tempo discreto
existe apenas em instantes especficos de tempo. Sistemas de tempo discreto sao aqueles
que manipulam sinais digitais.
Microcomputadores e microprocessadores digitais sao largamente utilizados na ind
ustria
atual, seja para fins de supervisao ou de controle dos processos. No entanto, um grande
n
umero de sistemas industriais sao de natureza analogica. Sempre que um microcomputador faz parte de um sistema analogico a presenca de conversores A/D e D/A se faz
necessaria.
Cada sinal analogico que sera processado por um computador digital deve primeiro ser
convertido de analogico para digital por um conversor A/D. Paralelamente, cada valor
digital que ira influenciar o sistema analogico devera primeiro ser convertido de digital
para analogico por um conversor D/A. Como a sada do computador digital nao muda ate
que os proximos calculos e conversoes D/A sejam completados, o sinal analogico gerado
por alguns conversores D/A sao mantidos constantes durante cada ciclo. Isto e feito
por um dispositivo chamado sample-and-hold (S/H). Veremos mais tarde que conversores
A/D tambem utilizam disposivos S/H.
7.1.1
Convers
ao A/D
A grande vantagem de se manipular variaveis discretas e que elas podem ser armazenadas e processadas em computadores digitais. Para isso basta transformar os valores
discretos em codigo binario.
A conversao para codigo binario nao e exata em geral. Sempre existe um erro entre o
Sistemas Lineares
180
valor discreto a ser codificado e codigo binario que representa o valor em questao. Por
exemplo, um sinal de tensao entre 0 e 10V pode ser representado em codigo binario de 4
bits de acordo com a tabela 7.1 e a figura 7.1.
tens
ao
analogica
0`
a 0.625
0.625 `a 1.25
1.25 a` 1.875
1.875 `a 2.5
2.5 a` 3.75
3.75 a` 4.375
4.375 `a 5
5`
a 5.625
5.625 `a 6.25
6.25 a` 6.875
6.875 `a 7.5
7.5 `
a 8.125
8.125 `a 8.75
8.75 a` 9.375
9.375 `a 10
representacao
bin
aria
0000
0001
0010
0011
0101
0110
0111
1000
1001
1010
1011
1100
1101
1110
1111
Tabela 7.1: Representacao de um sinal de tensao analogico nao negativo em codigo binario
de 4 bits
Cada incremento do codigo binario representa um salto de 24 = 6.25% em relacao ao
valor maximo do sinal analogico, isto e 6.25% de 10 volts no caso acima. Assim cada
codigo binario representa um intervalo de tensao analogica e portanto existe um erro de
quantizacao associado `a conversao. Num conversor de 4 bits o erro e de 6.25% ou seja,
uma relacao sinal/rudo de 20log(24 )dB. Para um conversor de 16 bits teramos um erro
de 0.0015% que corresponde `a uma relacao sinal/rudo de 20log(216 ) = 96.3dB. Bons
dispositivos de audio possuem relacao sinal/rudo entre 60 e 70 dB. Esta faixa e atingida
com conversores de 12 bits ou mais.
7.1.2
Convers
ao D/A e Sample-and-Hold
O dispositivo sample-and-hold (S/H) e normalmente utilizado na entrada de conversores A/D e na sada conversores D/A. A sua funcao basica e coletar amostras (sample) e
mante-la constante (hold) durante todo o intervalo de amostragem.
A figura 7.2 mostra o diagrama de blocos e um esquema eletronico simplificado do
dispositivo S/H. A chave logica s e controlada por um relogio. Com a chave na posicao 1
1
o dispositivo funciona como um circuito RC cuja funcao de transferencia e VVout
= RCs+1
.
in
A sada se torna praticamente igual `a entrada pois a frequencia de quebra do circuito
1
e escolhida grande em relacao `a maxima frequencia de quebra do espectro do
q = RC
sinal de entrada. Valores tpicos sao R = 1000 ohms e C = 30 1012 farads oque implica
1
= 5.3 MHz. Com a chave na posicao 1 a
numa frequencia de quebra de fq = 2RC
Sistemas Lineares
181
c
odigo bin
ario
1111
1110
1101
1100
1011
1010
1001
1000
0111
0110
0101
0100
0011
0010
0001
0000
-0
tens
ao anal
ogica
10
Figura 7.1: Representacao de um sinal de tensao analogico nao negativo em codigo binario
de 4 bits
1
2
sinal de controle
s
Vout
S/H
Vout
da chave
A seguir veremos de forma simplificada o funcionamento desse dispositivo quando acoplado a conversores A/D e D/A.
Quando um sinal analogico vai ser codificado, o primeiro passo e coletar as amostras
Sistemas Lineares
182
v(t)
S/H
A/D
(a)
w(t)
v(t)
v(t)
w(t)
v(t)
v(t)
w(t)
w(t)
T
h(t)
t
(b)
Figura 7.3: (a) Diagrama de blocos de um conversor A/D com sample-and-hold e (b)
funcionamento do sistema
A figura 7.3(b) mostra os sinais de entrada e sada nas duas fases de funcionamento
do dispositivo. Chave na posicao 1 corresponde `a fase segurar o sinal representada por
H (do ingles hold) e chave na posicao 2 `a fase de rastreamento T (do ingles tracking).
Cada amostra coletada do sinal de entrada e mantida constante no conversor durante
todo o intervalo de amostragem. O conjunto de amostras e atualizado apenas com a
chegada de uma nova amostra depois de cada intervalo de amostragem. Assim, o dispositivo S/H juntamente com o conversor A/D (codificador) possuem a funcao de amostrar
e segurar a amostra (ja codificada) durante todo o intervalo de amostragem.
Durante o processo de conversao D/A a sada do conversor pode flutuar muito. Para
evitar esse inconveniente utiliza-se um dispositivo S/H na sada do conversor. O S/H
mantem constante o valor da amostra precedente ate que uma nova amostra esteja disponvel (decodificada). O funcionamento desse dispositivo esta indicado na figura 7.4. O
sinal temporizador de controle indica as duas fases de funcionamento: manter o sinal de
entrada (H) e rastrear o sinal de entrada (T).
Note que o S/H executa duas operacoes: amostrar (sample) e segurar (hold). Afim de
obter um modelo matematico de funcionamento do S/H, estudaremos a seguir essas duas
Sistemas Lineares
183
v(t): sada anal
ogica do conversor
entrada digital
S/H
(a)
v(t)
y(t)
v(t)
y(t)
y(t)
v(t)
v(t)
H
y(t)
h(t)
t
(b)
Figura 7.4: (a) Conversor D/A com S/H e (b) Sinais de entrada e sada
operacoes separadamente.
A partir de agora assumiremos que a operacao amostrar do S/H pode ser representada
por um amostrador ideal, isto e ela pode ser representada pela multiplicacao do sinal a
ser amostrado por um trem de impulsos como ilustra a figura 7.5. A amostra do sinal
r(t) coletada no instante t = kT corresponde `a area do impulso que ocorre no instante
t = kT , isto e r(kT )(t kT ).
r(t)
r (t) =
k=
r(t)(t kT )
perodo de chaveamento: T
para kT t < kT + T
(7.1)
O bloco ZOH representa um sistema cuja funcao de transferencia pode ser obtida. Para
isso basta calcular a resposta impulsional desse sistema que chamaremos de zoh(t), isto
Sistemas Lineares
184
rh (t)
r (t)
ZOH
r(kT )
r(kT )(t kT )
kT
kT
kT+T t
e o valor especfico do sinal rh (t) obtido com r (t) = (t). Pela figura 7.6 deduzimos que
a resposta impulsional vale
rh (t) = zoh(t) = 1
para 0 t < T
(7.2)
que pode ser rescrita de uma forma mais conveniente com o auxlio da funcao degrau
unitario u(t) na forma
zoh(t) = u(t) u(t T )
(7.3)
A funcao de transferencia ZOH(s) do bloco ZOH pode entao ser calculada com o auxlio
da transformada de Laplace
ZOH(s) = L[zoh(t)] =
1 esT
s
s
(7.4)
Assim conclumos que o dispositivo S/H pode ser representado por um amostrador ideal
em cascata com um segurador de ordem zero como ilustra a figura 7.7.
r(t)
S/H
rh (t)
r (t)
r(t)
ZOH
rh (t)
Figura 7.7: Sample-and-Hold visto como um amostrador ideal em cascata com um segurador de ordem zero
importante notar que o amostrador ideal nao pode ser implementado na pratica
E
devido `a presenca de impulsos no sinal amostrado. No entanto a representacao do dispositivo S/H indicada na figura 7.7 possui duas propriedades interessantes: (i) existem
dispositivos S/H cujos comportamentos entrada/sada sao similares ao acima descrito
e (ii) o segurador de ordem zero pode ser tratado de maneira muito conveniente pela
Transformada Z, como veremos mais tarde.
Um outro ponto importante a ser notado e que no controle de sistemas normalmente se
assume a priori, por razoes de simplicidade, que o erro de conversao binaria e despresvel.
Isso implica que o conversor A/D pode ser representado por um amostrador ideal e o
conversor D/A por um segurador de ordem zero. Essas hipoteses sao comuns em todos os
livros classicos de controle e tambem serao assumidas nesse captulo sempre que houver
conversores A/D e D/A presentes na malha de controle.
Sistemas Lineares
7.2
185
Diferentemente dos sinais analogicos, que podem ser representados por funcoes do tipo
x(t) onde t e a variavel tempo contnuo , um sinal discreto e uma sequencia de valores
organizados no tempo e pode ser representado por funcao do tipo x(kT ) onde k e a
variavel tempo discreto (k = 0, 1, 2, ...) e T denota o intervalo de tempo entre dois
valores consecutivos de x(kT ). Neste captulo usaremos indistintamente os termos sinal
discreto ou sequencia.
Exemplo 7.1 A dinamica da vari
avel corrente no circuito da figura 7.8 e descrita por
uma equacao diferencial pois I(t) e uma vari
avel analogica (tempo contnio).
I(t)
vC (t) + RI(t) = 0
com vC (0) = v0 .
com I(0) = v0 /R. Logo:
+ RI(t) = 0
C I(t)
I(t) = I(0)et/RC ,
t 0.
Suponha agora que estamos interessados em saber os valores da corrente I(t) apenas
nos instantes t = kT , onde k = 0, 1, 2, . . . e T e um intervalo de tempo dado. Os valores
da corrente nesses instantes s
ao representados agora por uma sequencia I(kT ) e n
ao mais
por um sinal analogico como mostra a figura 7.9.
I(kT ) = I(0)ekT /RC
0 T 2T ...
t = kT
a = eT /RC
(7.5)
Sistemas Lineares
186
Sistemas contnuos s
ao aqueles que manipulam sinais analogicos e s
ao representados por
equacoes diferenciais, como e o caso do sistema na figura 7.8. Sistemas discretos s
ao
aqueles que manipulam sequencias e s
ao representados por equacoes recursivas, como e
o caso do sistema representado pela equacao recursiva 7.5. Note que o sistema discreto
n
ao e equivalente ao sistema contnuo que lhe deu origem pois apesar dos dois sistemas
representarem o comportamento exato da corrente no circuito nos instantes t = kT ,
apenas o sistema contnuo pode fornecer o valor exato da corrente para todo instante de
tempo t 0. Alem disso, o sistema discreto descrito pela equacao recursiva 7.5 pode ser
interpretado como um algortmo cuja evolucao define a dinamica da corrente do circuito
RC nos instantes t = kT .
Exemplo 7.2 Obtenha a equacao recursiva que rege o comportamento dinamico do circuito da figura 7.10 nos instantes t = kT sendo T um intervalo de tempo dado, k =
0, 1, 2, . . . uma vari
avel discreta e e(t) constante por trechos, isto e, e(t) = e(kT ) para
kT t < kT + T .
R
+
e(t)
+
C
x(t)
x(t0 ) = x(kT )
e(kT )
s
Logo:
e(kT )
e(kT ) + sRCx(kT )
1
+ RCx(kT )
=
s
RCs + 1
s(RCs + 1)
e(kT ) x(kT ) e(kT )
+
=
s
s + 1/RC
X(s) =
tkT
RC
kT t < kT + T
tkT +T
Sistemas Lineares
187
k = 0, 1, 2, . . .
onde a e b s
ao duas constantes dadas por:
a = eT /RC
b = 1 eT /RC
O sistema discreto dado pela equacao recursiva acima define o comportamento da corrente
I(t) (sada) em funcao da tensao e(t) (entrada) nos instantes t = kT como indicado na
figura 7.11. Mais tarde iremos calcular a funcao de transferencia discreta desse sistema
com o auxlio da transformada Z.
e(kT )
circuito
x(kT )
format(F16.8)
Y0=0.
Y1=0.
R1=0.
R0=0.
read(1,100)R2
Y2=3.*Y1 - 2.*Y0+2.*R2+5.*R0
Y0=Y1
Y1=Y2
R0=R1
R1=R2
write(2,100)Y2
go to 150
Sistemas Lineares
188
Controlador
r(t)
e(t)
+
-
A/D
e(kT)
Computador
u(kT)
D/A
u(t)
sistema
a ser
controlado
y(t)
Medidores
230 stop
240 end
O conversor A/D e tomado como periferico 1 com formato de leitura F16.8 e o pe assumido que o processador
riferico 2 e o conversor D/A com formato de escrita F16.8. E
espera no passo 150 ate que a vari
avel R2 esteja disponvel para leitura, da mesma forma
como ele esperaria caso a entrada de dados fosse via teclado. Tambem se assume que o
periferico 2 possui um buffer que armazena cada amostra da sada ate que a conversao
D/A seja efetuada. Identifique as condicoes iniciais, o sinal de entrada e o sinal de sada.
7.3
Transformada Z
Sistemas Lineares
7.3.1
189
Defini
c
ao e exemplos
x(kT )z k
k=0
x(kT ) = 0, k < 0
(7.6)
x(0)
1r
(7.7)
Para que o resultado da soma da serie seja dado pela formula acima e preciso que a serie
seja convergente, isto e a razao da serie deve possuir modulo menor que a unidade |r| < 1.
Veja o exemplo a seguir.
Exemplo 7.3 Calcule a transformada Z da sequencia degrau unit
ario (u(kT )) definida
como u(kT ) = 1 para k 0 e u(kT ) = 0 para k < 0.
Solucao: Pela definicao temos:
Z[u(kT )] =
k=0
u(kT )z
z k = 1 + z 1 + z 2 + . . .
k=0
A serie acima possui razao r = z 1 e a soma dos termos dessa serie e dada por 7.7
desde que a vari
avel complexa z esteja na regiao onde |r| = |z 1 | < 1. Nessas condicoes
temos:
1
z
U (z) = Z[u(t)] =
=
1
1z
z1
Analogamente `
a Transformada de Laplace e Fourier, a Transformada Z tambem possui
uma regiao de convergencia. Uma serie e convergente se em modulo a razao e menor que
Sistemas Lineares
190
Im[z]
Im[s]
crculo unitario
Re[z]
Plano z
Plano s
transformada de Laplace
transformada Z
U (z) =
Re[s]
U (s) =
z
z1
1
s
X(z)
z Rconv
Z 1 [X(z)]
Sistemas Lineares
7.3.2
191
Relac
ao com a transformada de Laplace
x (t) =
k=0
x(kT )(t kT )
seja a representacao do sinal x(t) amostrado com amostrador ideal (veja figura 7.5). Note
que a representacao do sinal amostrado x (t) e diferente da representacao da sequencia
x(kT ) obtida com os valores de x(t) nos instantes t = kT . Enquanto x (t) e um sinal
analogico com impulsos, a sequencia x(kT ) e um sinal discreto.
Tomemos agora a Transformada de Laplace da expressao acima:
x(kT )ekT s ,
x(t) = 0, t < 0
k=0
(7.8)
x(kT )z k
k=0
x(kT ) = 0, k < 0
(7.9)
Comparando (7.9) com (7.6) conclumos que a mudanca de variavel (7.8) define a relacao
entre a variavel s da transformada de Laplace do sinal amostrado x (t) e a variavel z da
transformada Z da sequencia x(kT ).
Veja por exemplo a relacao que existe entre os polos da transformada Z e Laplace do
degrau unitario indicadas na figura 7.13. O polo da transformada de Laplace esta na
origem s = 0. O polo da transformada Z esta em z = 1. Este mapeamento de s = 0 em
Laplace para z = 1 no plano Z e dado pela equacao (7.8).
Exemplo 7.4 (Fun
c
ao Pot
encia) Calcule a transformada Z da funcao potencia ak
onde a e uma constante e k 0 e uma vari
avel discreta.
Solucao: Com (7.7) temos:
Z[ak u(k)] =
k=0
ak z k =
z
1
=
1 a z 1
za
Sistemas Lineares
192
Z[sen(0 kT )u(kT )] =
sen(0 kT )z
k=0
1
2j
k=0
k=0
ej0 kT ej0 kT k
z
2j
ej0 kT z k ej0 kT z k
1
z
z
j
T
0
2j z e
z ej0 T
z sen(0 T )
= 2
z 2 z cos(0 T ) + 1
Z[(k)] =
(k)z k = 1
k=0
k=0
(k m)z k = z m
As figuras 7.15, 7.16 e 7.17 ilustram a relacao entre a localizacao dos polos da transformada Z do sinal e o seu comportamento temporal.
7.4
7.4.1
Propriedades da Transformada Z
Linearidade
Sistemas Lineares
193
polos de F (z)
crculo unitario
1.00 +
X
X
0.99
0
polo z = 1
10
f (k) = (0.7)k
1.0 +
0.9
crculo unitario
0.8
+
0.7
0.6
0.5
0.4
+
0.3
+
0.2
+
+
0.1
0
0
polo z = 0.7
10
f (k) = (0.7)k
1.1
+
0.9
crculo unitario
0.7
+
0.5
0.3
+
+
0.1
+
+
-0.1
-0.3
-0.5
-0.7
0
+
1
polo z = 0.7
Figura 7.15: Relacao entre localizacao polos e evolucao temporal
10
Sistemas Lineares
194
polos de F (z)
0.8
crculo unitario
0.6
0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1.0
0
polo z = 1
+
1
+
3
+
5
+
7
+
9
10
f (k) = (1.2)k
8
crculo unitario
6
+
4
+
+
+
0
+
-2
+
+
-4
+
-6
polo z = 1.2
10
f (k) = (1.2)k
7
crculo unitario
+
6
+
+
4
+
+
3
+
+
+
+
polo z = 1.2
1 +
0
+
1
10
Sistemas Lineares
195
polos de F (z)
1.0
crculo unitario
0.8
0.6
+
0.4
0.2
0 +
-0.2
-0.4
-0.6
-0.8
-1.0
polos z = ej 0.5
10
crculo unitario
0.5
0.4
0.3
+
0.2
0.1
+
0 +
+
-0.1
-0.2
0
10
crculo unitario
X
+
+
1
+
0 +
-1
-2
-3
-4
-5
-6
+
+
10
Sistemas Lineares
7.4.2
196
k=0
7.4.3
z1
x(k)z k
k=0
x(k + 1)z k
k=0
k=0
Supondo que a sequencia x(k) converge para um valor finito em regime, temos que X(z)
pode ter no maximo um polo sobre o crculo unitario e nenhum polo fora dele (veja figuras
7.15-7.17). Isso implica que a funcao auxiliar (z 1)X(z) nao pode ter polos sobre nem
fora do crculo unitario, ou seja devem estar dentro do crculo unitario. Logo para z 1
temos o seguinte resultado:
lim
z1
k=0
k=0
Sistemas Lineares
7.4.4
197
Obtenc
ao de F (z) a partir de F (s)
Z[f (kT )]
amostrador ideal
F (s)
L1 [F (s)]
f (t)
T
L[f (kT )]
f (kT ) ln(z)
s= T
sinal amostrado
F (z)
Z[F (s)]
Tabelas ou Teorema dos Resduos
F (z) =
R(Pi )
i=1
z
z esT
s=Pi
Sistemas Lineares
198
Polo m
ultiplo (multiplicidade m)
R(Pi ) =
z
dm1
1
(s Pi )m F (s)
m1
(m 1)! ds
z esT
s=Pi
1
.
(s+a)(s+b)
1
z
= ba
zeaT
1
z
= ab
zebT
s=a
s=b
Logo:
F (z) =
z
z
1
aT
ba ze
z ebT
7.4.5
Convoluc
ao Discreta
x1 ( )x2 (t )d =
x2 ( )x1 (t )d
Tempo discreto:
x1 (kT ) x2 (kT ) =
n=
n=
Normalmente temos x1 (t) = 0 e x2 (t) = 0 para t < 0 e nesses casos podemos tomar
t0 = 0 como limite inferior, tanto na integral como no somatorio.
Com a Transformada de Laplace vimos que convolucao no domnio do tempo se transforma em produto no domnio da frequencia. Mostraremos a seguir que isto tambem e
verdade em relacao `a convolucao discreta e a Transformada Z .
L[x1 (t) x2 (t)] = X1 (s)X2 (s)
Sistemas Lineares
199
y(kT )z
k=0
k=0
n=0
m=0 n=0
x2 (mT )z m
x1 (nT )z n
n=0
m=0
= X1 (z)X2 (z)
que e o resultado desejado.
Note que a convolucao de uma sequencia qualquer x1 (kT ) com um pulso unitario (kT )
resulta na propria sequencia x1 (kT ) pois, como ja vimos Z[(kT )] = 1.
(kT ) =
1 k=0
0 k=0
Z[(kT )] =
(kT )z k = 1
k=0
Logo:
Z[f (kT ) (kT )] = Z[f (kT )]Z[(kT )] = Z[f (kT )]
f (kT ) (kT ) = f (kT )
A funcao pulso unitario (kT ) tem (em relacao a Transformada Z ) as mesmas propriedades que a funcao impulso unitario (t) tem em relacao `a Transformada de Laplace
.
7.5
Transformada Z Inversa
Sistemas Lineares
7.5.1
200
M
etodo da divis
ao polinomial
k=0
Como normalmente X(z) e expressa em termos de uma fracao polinomial, isto e X(z) =
sendo N (z) e D(z) dois polinomios, temos:
N (z)
D(z)
N (z)
= x(0) + x(T )z 1 + x(2T )z 2 + . . .
D(z)
Para obter a igualdade acima atraves das regras usuais de divisao polinomial seguimos
o seguinte procedimento: Suponha que o grau de N (z) nao e superior ao grau de D(z) e
defina n=grau(D(z)). Construa dois polinomios auxiliares
(z 1 ) = z n N (z) ,
N
1 ) = z n D(z)
D(z
(z 1 ) por D(z
1 ) para encontrar os valores de x(0), x(T ),
Faca agora a divisao de N
x(2T ), . . . .
(z 1 )
N (z)
N
= x(0) + x(T )z 1 + x(2T )z 2 + . . .
(7.10)
=
1
D(z)
D(z )
Exemplo 7.8 Determine o valor numerico de x(4T ) dado que X(z) =
10z
(z1)(z2)
Solucao:
N (z)
10z
= 2
D(z)
z 3z + 2
1
2
D(z ) = z D(z) = 1 3z 1 + 2z 2
(z 1 ) = z 2 N (z) = 10z 1
N
X(z) =
D(z )
Logo, por igualdade polinomial com (7.10) conclumos que: x(0) = 0, x(T ) = 10, x(2T ) =
30, x(3T ) = 70, x(4T ) = 150.
Quando se deseja obter uma forma analtica para x(kT ) este metodo nao e adequado
e o metodo seguinte pode ser utilizado.
Sistemas Lineares
201
Tempo contnuo
x(t)
Tempo discreto
x(kT + T )
Comparacao
Transformada de Laplace
L[x(t)]
= sX(s) x(0)
Transformada Z
Z[x(kT + T )] = zX(z) zx(0)
e Z[x(kT + T )]
7.5.2
M
etodo das frac
oes parciais de X(z)/z
10z
.
(z1)(z2)
10
|
z2 z=1
= 10 e B =
10
|
z1 z=2
= 10.
Logo:
X(z) = 10
Lembrando que Z[ak ] =
z
za
temos:
z
z
+ 10
z1
z2
7.6
k = 0, 1, 2, . . .
Soluc
ao de Equac
oes recursivas
Sistemas Lineares
202
Z[x(kT )] =
Z[x(kT + T )] =
x(kT )z k ,
x(t) = 0 t < 0
k=0
x(kT + T )z
x(nT )z (n1)
(n = k + 1)
n=1
k=0
n=1
n=0
= zZ[x(nT )] zx(0)
que prova a propriedade desejada. Analogamente temos:
Z[x(kT + 2T )] = zZ[x(kT + T )] zx(T )
= z[zZ[x(kT )] zx(0)] x(T )
= z 2 X(z) z 2 x(0) zx(T )
Podemos enfim generalizar a propriedade do deslocamento no tempo aplicando sucessivamente os resultados acima e obtemos apos m sucessivos deslocamentos:
Z[x(k + m)] = z m X(z) z m x(0) z m1 x(1) zx(m 1)
(7.11)
(7.12)
x(0) = 0, x(1) = 1
Sistemas Lineares
203
z
za
z
z
z
z
=
=
z 2 + 3z + 2
(z + 2)(z + 1)
z+1 z+2
obtemos:
x(k) = (1)k (2)k ,
k = 0, 1, 2, . . .
Note que os polos da Transformada X(z) se tornam a base das exponenciais no tempo.
Logo uma sequencia x(kT ) e convergente, x(k) tende assintoticamente `a zero quando
k , se todos os polos da sua transformada X(z) sao em modulo inferiores `a unidade.
Confirme esse resultado na figuras 7.15-7.17.
Im[z]
sequencias convergentes: |polos| < 1
X
X
Re[z]
X
sequencias divergentes: |polos| > 1
PLANO z
1
z
z
j
T
2j z e 0
z ej0 T
possui dois polos (z = ej0 T e z = ej0 T ) que sao complexos conjugados (z = cos(0 T )
jsen(0 T )) e possuem modulo unitario indicando que a serie e oscilatoria sem amortecimento.
Exemplo 7.11 Determine a resposta do sistema descrito pela seguinte equacao recursiva:
x(k + 2) 3x(k + 1) + 2x(k) = u(k)
onde u(k) e o pulso unit
ario e x(k) = 0 para k 0.
Solucao: Para resolver a equacao acima precisamos das condicoes iniciais x(0) e x(1).
O valor de x(0) = 0 e dado e o valor de x(1) = 0 se obtem da pr
opria equacao recursiva
com k = 1. Alem disso, com a Transformada Z encontraremos:
z 2 X(z) 3zX(z) + 2X(z) = U (z)
Sistemas Lineares
204
z
za
z
z
+
z1 z2
obtemos finalmente:
k = 0, 1, 2, . . .
Exemplo 7.12 J
a vimos no exemplo 7.2 que no circuito RC da figura 7.10 onde e(t) e
constante por trechos (e(t) = e(kT ), kT t kT + T ) os sinais de entrada e sada nos
instantes t = kT s
ao dados pela equacao recursiva:
x(KT + T ) a x(kT ) = b e(kT )
a = eT /RC
b = 1 eT /RC
z
z1
temos:
X(z) =
z
za
z0.5
z
z
=
+
(z 0.5)(z 1)
z 0.5 z 1
encontramos:
k = 0, 1, 2, . . .
(Sistema Estavel)
k0
Sistemas Lineares
7.7
205
Func
ao de Transfer
encia Discreta e Estabilidade
Todo sistema discreto pode ser representado por um diagrama similar ao da figura 7.20
onde x(k) representa a sequencia de entrada dada, y(k) a sequencia de sada obtida, CI as
condicoes iniciais (que sao os n 1 valores iniciais da variavel de sada) e o bloco sistema
representa um sistema que sera descrito por uma equacao recursiva linear e invariante no
tempo (coeficientes constantes) do tipo
an y(k + n) + + a1 y(k + 1) + a0 y(k) = bm x(k + m) + + b1 x(k + 1) + b0 x(k) (7.13)
Por conveniencia de notacao estamos utilizando x(k), y(k) ao inves de x(kT ) e y(kT ).
Isto nao significa que estamos assumindo o intervalo T = 1 (tempo entre dois valores
consecutivos da sequencia). Esta notacao , muito utilizada em livros de controle, significa que estamos representando a sequencia numa escala de tempo normalizado k = t/T .
Note entretanto que os coeficientes da equacao recursiva dependem de T e nao podemos eliminar essa dependencia. Veja no caso do exemplo 7.12: poderamos rescrever a
equacao recursiva na forma x(k + 1) a x(k) = b e(k) mas os coeficientes a, b seriam os
mesmos anteriores que dependem de T e dos parametros fsicos do sistema (capacitancia
e resistencia nesse exemplo particular).
x(k)
SISTEMA
y(k)
C.I.
Figura 7.20: Sistema discreto generico
7.7.1
A resposta de todo sistema linear invariante no tempo pode ser decomposta em duas
parcelas: uma que depende do sistema e do sinal de entrada e outra que depende do
sistema e das condicoes iniciais. A primeira parcela chamaremos de Resposta de estado
zero ja que esta parcela indica como um sistema, inicialmente em repouso (condicoes
iniciais nulas), responde a um dado sinal de entrada. A segunda parcela chamaremos de
Resposta de Entrada Nula pois ela indica como um sistema se comporta quando e deixado
para responder livremente `as suas condicoes inicias (sem excitacao externa).
As respostas de Estado Zero e Entrada Zero de um sistema descrito por (7.13) pode
ser determinada atraves da Transformada Z .
Exemplo 7.13 Considere o circuito descrito no exemplo 7.12. Calcule a resposta de
Entrada Zero, para uma dada condicao inicial x(0) = x0 e a resposta de Estado Zero
para uma entrada generica e(k).
Sistemas Lineares
206
b
za
e F0 (z) =
fi (k)ci
(7.14)
i=0
onde ci = y(i) sao as condicoes iniciais do sistema, f (k) e fi (k) sao sequencias que
dependem apenas dos coeficientes da equacao recursiva (7.13), x(k) e a sequencia de
entrada e y(k) a de sada.
Da expressao acima observe que:
1. A sada de um sistema discreto depende dos parametros fsicos e do perodo de
amostragem que determinam os coeficientes da equacao recursiva e que por sua vez
determinam as funcoes f (k) e fi (k) em (7.14).
2. A sada de um sistema discreto depende da entrada que lhe e aplicada e essa dependencia se expressa atraves da convolucao discreta y(kT ) = f (kT ) x(kT ). Esta
parcela da resposta recebe o nome de resposta de Estado Zero.
yesz (kT ) = f (kT ) x(kT ).
3. A sada de um sistema depende das condicoes iniciais ci = y(iT ) (i = 0, . . . , n 1).
Esta parcela da resposta recebe o nome de resposta de Entrada Zero.
n1
yenz (kT ) =
fi (kT )ci
i=0
Sistemas Lineares
207
7.7.2
Quando as condicoes iniciais sao nulas a sada de um sistema discreto linear invariante
so depende da entrada e da Funcao de Transferencia Discreta, como pode ser visto em
(7.14).
Domnio do Tempo: y(k) = f (k) x(k).
Domnio da Frequencia: Y (z) = F (z)X(z).
A funcao f (k) = Z 1 [F (z)] recebe o nome de resposta ao pulso unitario pois f (k) e a
resposta do sistema quando as condicoes inciais sao nulas e a entrada e um pulso unitario
no instante k = 0. (X(z) = 1).
Defini
c
ao 7.1 (Sistemas Causais ou N
ao-Antecipativos) Um sistema discreto e dito
ser Causal (ou Nao-Antecipativo) se a resposta de Estado Zero para um pulso unit
ario e
nula para k < 0.
Num sistema causal o valor da resposta num dado instante de tempo y(kT ) nao depende
do sinal de entrada x(nT ) para valores de n > k. Caso contrario o valor da resposta no
instante t = kT passa a depender de valores futuros do sinal da entrada (x(t) para t > kT )
e que portanto ainda nao estao disponveis no instante t = kT .
Mostraremos a seguir que um sistema e causal quando o polinomio do numerador da
Funcao de Transferencia F (z) possui grau ao do denominador.
Com (7.13) temos que:
an y(k + n) + + a0 y(k) = bm x(k + m) + + b0 x(k)
Observe que y(k + n) depende de x(k + m) e portanto se m > n a sada no instante
k + n depende de valores da entrada em instantes de tempo futuros pois k + m > k + n
para todo k. Logo para que um sistema discreto seja causal devemos ter m n.
Defini
c
ao 7.2 (Estabilidade) Um sistema discreto linear invariante no tempo e exponencialmente est
avel se todos os polos da sua Funcao de Transferencia Pulsada possuem
modulo inferior `
a unidade. Caso contr
ario e dito ser instavel.
Sistemas Lineares
208
Pela definicao acima, note que a estabilidade e uma propriedade intrnseca do sistema.
Ela so depende dos parametros fsicos do mesmo. Nao depende da entrada nem das
condicoes iniciais. O nome exponencialmente estavel apenas enfatiza que os polos da
Funcao de Transferencia serao a base de exponenciais no domnio do tempo e portanto
a resposta ao pulso converge exponencialmente para zero.
Exemplo 7.14 Verifique a estabilidade do sistema cuja funcao de transferencia e
F (z) =
z
(z 1)(z + 2)
z
za
1 z
1 z
3z 1 3z +2
temos:
1
1
Z 1 = f (k) = (1)k (2)k
3
3
7.8
Sistemas Amostrados
Sistemas Lineares
209
sera discretizada para posterior tramento num computador digital, isto e apenas os valores y(kT ) serao considerados. Precisamos entao saber qual seria o sistema discreto
equivalente que tem as sequencias x(kT ) como entrada e y(kT ) como sada, como indica
a figura 7.21(b). O sistema discreto que procuramos possui as mesmas informacoes do
sistema amostrado e alem disso a transformada Z pode ser aplicada diretamente.
x (t)
x(t)
G(s)
y(t)|t=kT
(a)
x(kT )
G(z)
y(kT )
(b)
Note que a entrada do bloco analogico G(s) na figura 7.21(a) e um sinal amostrado
onde os impulsos possuem areas de valores x(kT ) e a entrada do bloco discreto G(z) na
figura 7.21(b) e uma sequencia de valores x(kT ). A sada y(t) figura 7.21(a) e analogica
mas apenas os valores y(kT ) medidos nos instantes t = kT sao de interesse. Ja a sada
na figura 7.21(b) e a propria sequencia y(kT ).
A seguir apresentamos um procedimento para, dado o sistema analogico G(s), encontrar
o bloco discreto equivalente G(z). E equivalencia `a qual nos referimos e no sentido de
que os valores do sinal de entrada x(kT ) e sada y(kT ) do sistema discreto sao os mesmos
do sistema contnuo x(t), y(t) nos instantes t = kT .
Seja entao g(t) = L1 [G(s)] a resposta impulsional do sistema representado por G(s).
Da, a resposta y(t) e dada pela convolucao contnua da entrada com a resposta impulsional.
ENTRADA: x (t) =
SAIDA: y(t) =
k
n=0
k=0
x(kT )(t kT ),
0 t kT .
Sistemas Lineares
210
y(t) =
n=0
k
=
n=0
0 t kT
x(nT ) g(t nT )
y(t)|t=kT = y(kT ) =
n=0
x(nT )g(kT nT )
(7.15)
(7.16)
Exemplo 7.15 Considere o sistema amostrado da figura 7.21(a) com G(s) = (s+a)(s+b)
.
Calcule a funcao de transferencia discreta G(z) entre a sequencia x(kT ) de entrada e
y(kT ) de sada indicada na figura 7.21(b).
z
=
sT
ze
b a z eaT
s=a
z
z
R(P2 ) = (s + b)G(s)
=
z esT s=b a b z ebT
z
z
G(z) =
aT
ba ze
z ebT
R(P1 ) =
(s + a)G(s)
que e a funcao de transferencia discreta desejada. Alem disso, com a notacao Y (z) =
Z[y(kT )] e X(z) = Z[x(kT )] podemos escrever a relacao de entrada/sada do sistema
Y (z) = G(z)X(z).
Sistemas Lineares
X(s)
211
X (s)
G(s)
Y (s)
Y (s)
Y (s)
G(s)
Y (s)
X1 (s)
Y (s)
Y (s)
X2 (s)
ou de forma equivalente
X1 (s)
X1 (s)
Y (s)
X2 (s)
X2 (s)
Sistemas Lineares
212
ba
1
1
s+a s+b
(eat ebt )
ba
(eakT ebkT )
g(t)|t=kT = g(kT ) =
ba
tomando a Transformada Z da sequencia g(kT ) obtemos:
G(z) = Z[g(kT )] =
ba
z
z
aT
ze
z ebT
e finalmente temos:
Y (z) = G(z)X(z) Frequencia
y(kT ) = g(kT ) x(kT ) Tempo
A maioria dos sistemas amostrados possuem algum dispositivo sample-and-hold no seu
interior. A figura 7.23 mostra um sistema desse tipo. A funcao ZOH(s) e a funcao de
transferencia do segurador de ordem zero como indicado em (7.4).
x (t)
x(kT )
y(t)|t=kT
ZOH(s)
G(s)
Segurador
de Ordem
Zero
Sistema
a ser
Controlado
T
Amostrador
Ideal
Da figura 7.23 vamos definir a funcao auxiliar H(s) = ZOH(s)G(s). O problema agora
e encontrar a funcao de transferencia discreta H(z) que corresponde `a funcao H(s). Para
isso vamos utilizar a notacao H(z) = Z[H(s)]. Lembrando que esT = z temos com (7.4)
H(z) = Z[ZOH(s)G(s)] = Z[
G(s)
1 eT s
G(s)] = (1 z 1 )Z
s
s
(7.17)
1
Exemplo 7.16 Considere o sistema da figura 7.23 onde G(s) = s(s+1)
. Calcule a funcao
de transferencia Pulsada entre a sequencia x(kT ) e a sada y(t) nos instantes t = kT
com T = 1seg.
1
1
1 eT s
= (1 z 1 )Z 2
s
s(s + 1)
s (s + 1)
Sistemas Lineares
213
Definindo F (s) =
1
s2 (s+1)
onde:
R(P1 ) = (s + 1)F (s)
e
d
ds
d
ds
z
z eT s
=
s=1
z
z e1
z
z esT
s=0
1
z
=
s + 1 z esT
s=0
1
z
1
z
=
+
(T esT )
2
sT
(s + 1) z e
s + 1 (z es T )2
Tz
z
+
=
z 1 (z 1)2
R(P2 ) =
Logo:
s2 F (s)
F (z) =
e portanto:
s=0
z
z
z
+
1
2
ze
(z 1)
z1
H(z) = (1 z 1 )F (z) =
que e a funcao de transferencia desejada.
1
z1
+
1
1
ze
z1
Exemplo 7.17 Considere o circuito RC do exemplo 7.2 onde o sinal de entrada e constante por trechos. Represente o sinal constante por trechos como sendo a sada de um
segurador de ordem zero como indicado na figura 7.24. Calcule a equacao recursiva que
rege o comportamento do sistema nos instantes t = kT . Calcule tambem a resposta do
circuito para um degrau unit
ario.
e (t)
e(kT )
ZOH(s)
e(t)
y(t)|t=kT
G(s)
T
Segurador
de Ordem
Zero
Amostrador
Ideal
Circuito RC
1 eT s
1
1
1
= (1 z )Z
H(z) = Z
s
RCs + 1
s(RCs + 1)
F (s)
Sistemas Lineares
214
s=0
z
z1
H(z) = (1 z 1 )
s=1/RC
z
zeT /RC
1 eT /RC
z
z
=
z 1 z eT /RC
z eT /RC
z
z1
e portanto Y (z) =
b
z
]
z 1z a
k = 0, 1, 2, . . .
y(kT ) = Z 1 [Y (z)] = Z 1 [
= 1 ekT /RC ,
x (t)
1
s+a
x(t)
y1 (t)
y1 (t)
1
s+b
y2 (t)|t=kT
(a)
x (t)
1
s+a
1
s+b
y(t)|t=kT
Figura 7.25: (a) Dois sistemas amostrados em cascata; (b) Dois sistemas contnuos em
cascata
1
1
] Z[x(kT )] Y1 (z) = Z[
] X(z)
s+a
s+a
(b)
Sistemas Lineares
215
1
1
] Z[y1 (kT )] Y2 (z) = Z[
] Y1 (z)
s+b
s+b
z
zeaT
1
e Z[ s+b
]=
1
1
] Z[
] X(z)
s+a
s+b
z
zebT
1
1
=
(s + a)(s + b)
ba
z
z
aT
ze
z ebT
e portanto
Z
1
1
1
Z
=Z
(s + a)
(s + b)
(s + a)(s + b)
G(s)
s
1
Sistemas Lineares
216
Note que G(s) e um sistema fsico e portanto causal. Logo G(s)/s tambem e causal e
portanto h0 (t) = 0 para t < 0. Assim:
Z[h0 (kT T )] =
k=0
h0 (kT T )z k = h0 (T ) + h0 (0)z 1 + h0 (T )z 2 + . . .
como h0 (T ) = 0 temos:
Z[h0 (kT T )] = z 1 [h0 (0) + h0 (T )z 1 + h0 (2T )z 2 + . . .
1
= z [
k=0
1 eT s
G(s)
s
G(s)
]
s
logo:
Z[x(kT )] = H(z) Z[y(kT )]
X(z)
Y (z)
7.9
Sistemas Realimentados
Sistemas Lineares
programa do
computador
e(kT )
e(t)
modelo do
conversor D/A
com S/H
modelo do
modelo do
conversor
A/D
sinal de
referencia +
r(t)
217
u(kT )
G1 (z)
ZOH(s)
modelo do
sistema a ser
controlado
vari
avel
controlada
u(t)
G2 (s)
y(t)
(a)
computador
algortmo
(eq. recursiva)
sinal de
referencia +
r(kT )
e(kT )
modelo do
conversor D/A
com S/H
u(kT )
ZOH(s)
G1 (z)
modelo do
sistema a ser
controlado
u(t)
modelo do
medidor com
conversor A/D
G2 (s)
y(kT )
vari
avel
controlada
(b)
sinal de
referencia +
r(kT )
e(kT )
algortmo
(eq. recursiva)
u(kT )
G1 (z)
vari
avel
controlada
G3 (z)
y(kT )
(c)
Sistemas Lineares
sistema, controlador e
medidor anal
ogico
S/H
referencia
e(t)
r(t)
218
+
-
e (t)
c (t)
G(s)
ZOH(s)
(a)
T
v(t)
variavel controlada
H(s)
filtros auxiliares
S/H
referencia
e(t)
r(t)
+
-
sistema, controlador e
medidor digital
S/H
e (t)
c (t)
G(s)
ZOH(s)
v(t)
(b)
variavel controlada
ZOH(s)
H(s)
filtros auxiliares
Figura 7.27: Sistema de controle digital com medidor analogico (a) e digital (b)
C(z)
G(z)
e portanto:
C(z) =
G(z)
R(z)
1 + GH(z)
que expressa a relacao entre a sequencia r(kT ) e c(kT ) no sistema da figura 7.27(a).
Passemos agora ao sistema da figura 7.27(b). Para encontrar a funcao de transferencia
entre as sequencias e(kT ) e c(kT ) temos
C(z) = G(z)E(z)
C(z)
C(z)
+ H(z)C(z) = (1 + G(z)H(z))
G(z)
G(z)
Sistemas Lineares
219
e portanto:
C(z) =
G(z)
R(z)
1 + G(z)H(z)
que expressa a relacao entre a sequencia r(kT ) e c(kT ) no sistema da figura 7.27(b).
e(t)
r(t)
+
e (t)
ZOH(s)
Va (t)
c (t)
1
s(s+1)
T=1seg
T=1seg
e1 z + 1 2e1
0, 368z + 0, 264
= 2
2
1
1
z (1 + e )z + e
z 1, 368z + 0, 368
Logo:
C(z)
0, 368z + 0, 264
= 2
R(z)
z z + 0, 632
z
z1
temos:
0, 368z + 0, 264 z
z 2 z + 0, 632 z 1
Sistemas Lineares
7.10
220
7.11
Resposta em Frequ
encia
Sistemas Lineares
221
z = eT s . Esta relacao mostra como os polos do plano s de Laplace sao mapeados para o
plano z.
Exemplo 7.21 O sinal y(t) = eat cos(bt), t 0 com aT = 0, 3567 e bT = /4 resulta
na seguinte sequencia para T = 1seg.
y(kT ) = (e3567 )k cos(k/4),
= 0, 7k cos(k/4)
k = 0, 1, 2, . . .
F(z)
Y(z)
Sistemas Lineares
222
1
2
z
zej0 T
Y (z) =
z
zej0 T
temos:
z
z
1
F (z)
+
j
T
0
2 ze
z ej0 T
Supondo F (z) estavel, isto e, que todos os seus polos estejam dentro do crculo unitario,
a resposta em regime permanente e dada pelos termos da expansao por fracoes parciais
de Y (z) correspondentes aos polos sobre o crculo unitario, pois todos os outros termos
irao desaparecer quando k . Dessa forma, expandindo Y (z)/z por fracoes parciais e
desprezando os termos associados aos polos dentro do crculo unitario ficamos com:
Yss (z) =
1
z
z
j0 T
F
(e
)
+
F (ej0 T )
j
T
2 ze 0
z ej0 T
z
za
= ak vem:
1
(ej0 kT )M ej + (ej0 kT )M ej
2
= M cos(0 T k + )
yss (kT ) =
7.12
Problemas Complementares
Problema 7.6 O sinal de entrada do circuito na figura 7.31 e constante por trechos, isto
e v(t) = v(kT ) para kT t < kT + T . Sendo T = 1 segundo pede-se: (i) a equacao
recursiva que define o comportamento entrada/sada nos instantes t = KT ; (ii) a funcao
de transferencia discreta; (iii) a resposta de estado zero ao degrau de 10 volts; (iv) a
resposta de entrada zero para vc (0) = 1V, vc (T ) = 0V ; (v) a resposta total para uma
entrada degrau unit
ario e condicoes iniciais vc (0) = 3V, vc (T ) = 0V .
Problema 7.7 O sistema da figura 7.32 mostra um esquema de controle de velocidade de um motor DC controlado pela armadura. O perodo de amostragem e de T =
Sistemas Lineares
223
R = 3
+
L = 1H
C = 1F
v(t)
vc (t)
e(kT)
+
computador
u(kT)
S/H
va (t)
motor
w(t)
Sistemas Lineares
(a)
224
x 1 = x2
x 2 = x1 x2 + e
(b)
x1 (k + 1) = x2 (k)
x2 (k + 1) = x1 (k) x2 (k) + e(k)
x1
sistema
+
R
x(t)
v(t)
-
Figura 7.34: Entrada: tensao x(t) ; sada: tensao v(t) ; R=1 , C=1 F
ZOH(s)
-
T=1seg
1
s+2
1
s+1
T=1seg
y(t)
Bibliografia
K. OGATA, Engenharia de Controle Moderno, Prentice Hall 1993.
K. OGATA, Discrete Time Control Systems, Prentice Hall, 1995.
G.F.Franklin, J.D.Powell, M.L.Workman, Digital control of Dinamic Systems,
Addison-Wesley, 1990.