Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

Lista de Exercícios
Estruturas para sistemas em tempo discreto

1. Problema 6.25 do livro texto: Um sistema LIT causal é definido pelo diagrama
de fluxo da figura a seguir. Esta figura representa o sistema como uma cascata
de um sistema de segunda ordem em série com um sistema de primeira ordem.

a. Qual é a função do sistema em cascata global?

b. O sistema global é estável?
c. O sistema global é de fase mínima?
d. Desenhe o diagrama de fluxo de sinais na implementação forma direta II
transposta.
e. Escreva a equação de diferenças que representa o sistema.

Resolução

a. Função do sistema:

1 + 0,81𝑧 −2 1 + 2𝑧 −1
𝐻 (𝑧 ) = ⋅
1 − 0,3𝑧 −1 − 0,4𝑧 −2 1 + 0,8𝑧 −1

A primeira parte do polinômio representa um sistema de segunda ordem que será dividido
em dois sistemas de primeira ordem para identificar os polos e zeros.

1 + 0,81𝑧 −2 1 + 2𝑧 −1 𝑁 1 + 2𝑧 −1
𝐻 (𝑧 ) = ⋅ = ⋅
1 − 0,3𝑧 −1 − 0,4𝑧 −2 1 + 0,8𝑧 −1 𝐷 1 + 0,8𝑧 −1

Numerador 𝑁:

𝑁 = 1 + 0,81𝑧 −2 = (𝑎 + 𝑏𝑧 −1 )(𝑎 − 𝑏𝑧 −1 )

𝑁 = 𝑎2 + 𝑎𝑏𝑧 −1 − 𝑎𝑏𝑧 −1 − 𝑏2 𝑧 −2 = 𝑎2 − 𝑏2 𝑧 −2

𝑁 = 𝑎2 − 𝑏2 𝑧 −2 = 1 + 0,81𝑧 −2

𝑎2 = 1 ⟹ 𝑎 = 1 −𝑏2 = 0,81 ⟹ 𝑏 = √−0,81 = 𝑗0,9

Estruturas 1
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

𝑁 = (1 + 𝑗0,9𝑧 −1 )(1 − 𝑗0,9𝑧 −1 )

Denominador 𝐷:

𝐷 = 1 − 0,3𝑧 −1 − 0,4𝑧 −2 = (𝑐 + 𝑑𝑧 −1 )(𝑐 + 𝑒𝑧 −1 )

𝐷 = 𝑐 2 + 𝑐𝑒𝑧 −1 + 𝑐 𝑑𝑧 −1 + 𝑑𝑒𝑧 −2 = 𝑐 2 + 𝑐(𝑒 + 𝑑 )𝑧 −1 + 𝑑𝑒𝑧 −2

𝐷 = 𝑐 2 + 𝑐(𝑒 + 𝑑 )𝑧 −1 + 𝑑𝑒𝑧 −2 = 1 − 0,3𝑧 −1 − 0,4𝑧 −2

𝑐2 = 1 ⟹ 𝑐 = 1 (𝑒 + 𝑑 )𝑧 −1 = −0,3𝑧 −1 ⟹ 𝑒 + 𝑑 = −0,3 𝑑𝑒𝑧 −2 = −0,4𝑧 −2 ⟹ 𝑑𝑒 = −0,4

𝑒 + 𝑑 = −0,3
(1)
𝑑𝑒 = −0,4

Resolvendo o sistema de equações (1):

𝑑 = −0,8 𝑒 = 0,5

𝐷 = (1 − 0,8𝑧 −1 )(1 + 0,5𝑧 −1 )

Portanto a função do sistema com os polos e zeros identificados fica:

𝟏 + 𝟎, 𝟖𝟏𝒛−𝟐 𝟏 + 𝟐𝒛−𝟏 (𝟏 + 𝒋𝟎, 𝟗𝒛−𝟏 )(𝟏 − 𝒋𝟎, 𝟗𝒛−𝟏 )(𝟏 + 𝟐𝒛−𝟏 )

𝑯(𝒛) = ⋅ =
𝟏 − 𝟎, 𝟑𝒛−𝟏 − 𝟎, 𝟒𝒛−𝟐 𝟏 + 𝟎, 𝟖𝒛−𝟏 (𝟏 − 𝟎, 𝟖𝒛−𝟏 )(𝟏 + 𝟎, 𝟓𝒛−𝟏 )(𝟏 + 𝟎, 𝟖𝒛−𝟏 )

O sistema tem zeros em 𝑧 = −𝑗0,9, 𝑧 = 𝑗0,9 e 𝑧 = −2, e polos em 𝑧 = 0,8, 𝑧 = −0,5 e 𝑧 =

−0,8.

b. Verificando o ponto (a), o sistema é estável porque todos os polos estão dentro do
círculo de raio unitário. Isto garante a estabilidade de um sistema causal.

c. O sistema, de acordo com o ponto (a) não é de fase mínima devido a que tem um zero
fora do círculo de raio unitário em 𝑧 = −2.

Estruturas 2
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

d. Para implementar o sistema na forma direta II transposta deveremos reescrever a

função do sistema.

(1 + 𝑗0,9𝑧 −1 )(1 − 𝑗0,9𝑧 −1 )(1 + 2𝑧 −1 ) 1 + 0,81𝑧 −2 1 + 2𝑧 −1 𝑁

𝐻 (𝑧 ) = −1 −1 −1
= −1 −2
⋅ −1
=
(1 − 0,8𝑧 )(1 + 0,5𝑧 )(1 + 0,8𝑧 ) 1 − 0,3𝑧 − 0,4𝑧 1 + 0,8𝑧 𝐷

Numerador 𝑁:

𝑁 = (1 + 0,81𝑧 −2 )(1 + 2𝑧 −1 )

𝑁 = 1 + 2𝑧 −1 + 0,81𝑧 −2 + 1,62𝑧 −3

Denominador 𝐷:

𝐷 = (1 − 0,3𝑧 −1 − 0,4𝑧 −2 )(1 + 0,8𝑧 −1 )

𝐷 = 1 + 0,8𝑧 −1 − 0,3𝑧 −1 − 0,24𝑧 −2 − 0,4𝑧 −2 − 0,32𝑧 −3

𝐷 = 1 + 0,5𝑧 −1 − 0,64𝑧 −2 − 0,32𝑧 −3

1 + 2𝑧 −1 + 0,81𝑧 −2 + 1,62𝑧 −3
𝐻 (𝑧 ) =
1 + 0,5𝑧 −1 − 0,64𝑧 −2 − 0,32𝑧 −3

e. Equação de diferenças:

1 + 2𝑧 −1 + 0,81𝑧 −2 + 1,62𝑧 −3 𝑌 (𝑧 )
𝐻 (𝑧 ) = −1 −2 −3
=
1 + 0,5𝑧 − 0,64𝑧 − 0,32𝑧 𝑋 (𝑧 )

𝑌 (𝑧)(1 + 0,5𝑧 −1 − 0,64𝑧 −2 − 0,32𝑧 −3 ) = 𝑋(𝑧)(1 + 2𝑧 −1 + 0,81𝑧 −2 + 1,62𝑧 −3 )

Estruturas 3
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

𝑌(𝑧) + 0,5𝑌 (𝑧)𝑧 −1 − 0,64𝑌 (𝑧)𝑧 −2 − 0,32𝑌 (𝑧)𝑧 −3 =

= 𝑋(𝑧) + 2𝑋(𝑧)𝑧 −1 + 0,81𝑋(𝑧)𝑧 −2 + 1,62𝑋 (𝑧)𝑧 −3

Aplicando transformada z inversa:

𝑦[𝑛] + 0,5𝑦[𝑛 − 1] − 0,64𝑦[𝑛 − 2] − 0,32𝑦[𝑛 − 3]

= 𝑥 [𝑛] + 2𝑥 [𝑛 − 1] + 0,81𝑥 [𝑛 − 2] + 1,62𝑥[𝑛 − 3]

Reescrevendo a equação:

𝒚[𝒏] = 𝒙[𝒏] + 𝟐𝒙[𝒏 − 𝟏] + 𝟎, 𝟖𝟏𝒙[𝒏 − 𝟐] + 𝟏, 𝟔𝟐𝒙[𝒏 − 𝟑] − 𝟎, 𝟓𝒚[𝒏 − 𝟏] + 𝟎, 𝟔𝟒𝒚[𝒏 − 𝟐]

+ 𝟎, 𝟑𝟐𝒚[𝒏 − 𝟑]

--------------------------------------------------------------------------------------------------------------------------

2. Problema 6.12 do livro texto: O diagrama de fluxo da figura descreve um sistema LIT.

a. Determine a função de transferência.

b. Determine a equação de diferenças que relaciona a entrada 𝑥[𝑛] com a saída
𝑦[𝑛].
c. Desenhe o diagrama de fluxo de acordo com a função de transferência obtida no
ponto a.
d. Calcule a resposta ao impulso ℎ[𝑛].
e. Verifique se o sistema é estável.
f. Verifique se o sistema tem fase mínima.

Resolução

a. Função do sistema:

Para facilitar a resolução definiremos as variáveis 𝑤1 [𝑛], 𝑤2 [𝑛] e 𝑤3 [𝑛]:

Estruturas 4
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

𝑤1 [𝑛] = −𝑥 [𝑛] + 𝑤2 [𝑛] + 𝑤3 [𝑛]

𝑤2 [𝑛] = 𝑥 [𝑛 − 1] + 2𝑤3 [𝑛]

𝑤3 [𝑛] = 𝑤2 [𝑛 − 1] + 𝑦[𝑛 − 1]

𝑦[𝑛] = 2𝑤1 [𝑛]

Aplicando transformada z às equações:

𝑊1 (𝑧) = −𝑋(𝑧) + 𝑊2 (𝑧) + 𝑊3 (𝑧) (2)

𝑊2 (𝑧) = 𝑧 −1 𝑋(𝑧) + 2𝑊3 (𝑧) (3)

𝑊3 (𝑧) = 𝑧 −1 𝑊2 (𝑧) + 𝑧 −1 𝑌(𝑧) (4)

Substituindo a equação (4) na (2):

𝑊1 (𝑧) = −𝑋(𝑧) + 𝑊2 (𝑧) + 𝑧 −1 𝑊2 (𝑧) + 𝑧 −1 𝑌(𝑧)

𝑊1 (𝑧) = −𝑋(𝑧) + 𝑊2 (𝑧)(1 + 𝑧 −1 ) + 𝑧 −1 𝑌(𝑧) (5)

Substituindo a equação (4) na (3):

𝑊2 (𝑧) = 𝑧 −1 𝑋 (𝑧) + 2(𝑧 −1 𝑊2 (𝑧) + 𝑧 −1 𝑌(𝑧))

𝑊2 (𝑧) = 𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑊2 (𝑧) + 2𝑧 −1 𝑌(𝑧)

𝑊2 (𝑧) − 2𝑧 −1 𝑊2 (𝑧) + 𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧)

𝑊2 (𝑧)(1 − 2𝑧 −1 ) = 𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧)

𝑧 −1 𝑋 (𝑧) + 2𝑧 −1 𝑌(𝑧)
𝑊2 (𝑧) = (6)
1 − 2𝑧 −1
Equação (6) em (5):

𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧)
𝑊1 (𝑧) = −𝑋(𝑧) + (1 + 𝑧 −1 ) + 𝑧 −1 𝑌(𝑧) (7)
1 − 2𝑧 −1
Resolvendo o termo:

Estruturas 5
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧)
(1 + 𝑧 −1 )
1 − 2𝑧 −1

Numerador:

(𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧))(1 + 𝑧 −1 ) = 𝑧 −1 𝑋 (𝑧) + 2𝑧 −1 𝑌 (𝑧) + 𝑧 −2 𝑋(𝑧) + 2𝑧 −2 𝑌(𝑧)

Portanto o termo fica:

𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧) + 𝑧 −2 𝑋(𝑧) + 2𝑧 −2 𝑌(𝑧)

1 − 2𝑧 −1

Voltando para a equação (6):

𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌 (𝑧) + 𝑧 −2 𝑋(𝑧) + 2𝑧 −2 𝑌(𝑧)

𝑊1 (𝑧) = −𝑋(𝑧) + + 𝑧 −1 𝑌(𝑧)
1 − 2𝑧 −1

Aplicando denominador comum:

−𝑋(𝑧)(1 − 2𝑧 −1 ) + 𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌 (𝑧) + 𝑧 −2 𝑋(𝑧) + 2𝑧 −2 𝑌(𝑧) + 𝑧 −1 𝑌(𝑧)(1 − 2𝑧 −1 )

𝑊1 (𝑧) =
1 − 2𝑧 −1

𝑊1 (𝑧)
−𝑋(𝑧) + 2𝑧 −1 𝑋(𝑧) + 𝑧 −1 𝑋(𝑧) + 2𝑧 −1 𝑌(𝑧) + 𝑧 −2 𝑋 (𝑧) + 2𝑧 −2 𝑌(𝑧) + 𝑧 −1 𝑌(𝑧) − 2𝑧 −1 𝑌(𝑧)
=
1 − 2𝑧 −1

−𝑋(𝑧) + 3𝑧 −1 𝑋 (𝑧) + 𝑧 −2 𝑋(𝑧) + 3𝑧 −1 𝑌 (𝑧)

𝑊1 (𝑧) =
1 − 2𝑧 −1

−𝑋(𝑧) + 3𝑧 −1 𝑋 (𝑧) + 𝑧 −2 𝑋(𝑧) + 3𝑧 −1 𝑌(𝑧)

𝑌 (𝑧 ) = 2
1 − 2𝑧 −1

𝑌 (𝑧)(1 − 2𝑧 −1 ) = 2(−𝑋(𝑧) + 3𝑧 −1 𝑋(𝑧) + 𝑧 −2 𝑋(𝑧) + 3𝑧 −1 𝑌 (𝑧))

𝑌 (𝑧) − 2𝑧 −1 𝑌(𝑧) = −2𝑋(𝑧) + 6𝑧 −1 𝑋(𝑧) + 2𝑧 −2 𝑋(𝑧) + 6𝑧 −1 𝑌 (𝑧)

𝑌 (𝑧) − 8𝑧 −1 𝑌(𝑧) = −2𝑋(𝑧) + 6𝑧 −1 𝑋(𝑧) + 2𝑧 −2 𝑋(𝑧)

𝒀(𝒛) −𝟐 + 𝟔𝒛−𝟏 + 𝟐𝒛−𝟐

𝑯(𝒛) = =
𝑿(𝒛) 𝟏 − 𝟖𝒛−𝟏

b. Equação de diferenças:

𝑌 (𝑧) − 8𝑧 −1 𝑌(𝑧) = −2𝑋(𝑧) + 6𝑧 −1 𝑋(𝑧) + 2𝑧 −2 𝑋(𝑧)

Aplicando transformada z inversa:

𝒚[𝒏] − 𝟖𝒚[𝒏 − 𝟏] = −𝟐𝒙[𝒏] + 𝟔𝒙[𝒏 − 𝟏] + 𝟐𝒙[𝒏 − 𝟐]

Estruturas 6
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

c. Diagrama de fluxo:

d. Resposta ao impulso ℎ[𝑛]:

Por divisão longa

1,2187
𝐻 (𝑧) = −0,25𝑧 −1 − 0,7812 −
1 − 8𝑧 −1

𝒉[𝒏] = −𝟎, 𝟐𝟓𝜹[𝒏 − 𝟏] − 𝟎, 𝟕𝟖𝟏𝟐𝜹[𝒏] − 𝟏, 𝟐𝟏𝟖𝟕(𝟖)𝒏 𝒖[𝒏]

e. Estabilidade

𝑌 (𝑧) −2 + 6𝑧 −1 + 2𝑧 −2
𝐻 (𝑧 ) = =
𝑋 (𝑧 ) 1 − 8𝑧 −1

O sistema tem um polo em 𝑧 = 8, portanto o sistema, sendo causal, é instável devido

a que tem um polo fora do círculo de raio unitário e a RDC não contém o mesmo

f. Fase mínima

Considerando o numerador da função de transferência:

𝑁 = −2 + 6𝑧 −1 + 2𝑧 −2 = 0

Dividindo tudo por -2:

𝑁 = 1 − 3𝑧 −1 − 𝑧 −2 = 0

(1 + 𝑎𝑧 −1 )(1 + 𝑏𝑧 −1 )

Estruturas 7
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

1 + 𝑏𝑧 −1 + 𝑎𝑧 −1 + 𝑎𝑏𝑧 −2

1 + (𝑏 + 𝑎)𝑧 −1 + 𝑎𝑏𝑧 −2

𝑏 + 𝑎 = −3 𝑎𝑏 = −1 (8)

𝑎 ≅ −3,3 𝑏 ≅ 0,3

𝑁 = (1 − 3,3𝑧 −1 )(1 − 0,3𝑧 −1 )

O sistema tem um zero em 𝑧 ≅ −3,3, fora do círculo de raio unitário, portanto não tem
fase mínima

--------------------------------------------------------------------------------------------------------------------------

3. Problema 6.10 do livro texto: O diagrama de fluxo da figura representa um sistema LIT:

a. Escreva a equação de diferenças representada pelo diagrama de fluxo.

b. Desenhe o diagrama de fluxo em cascata com dois sistemas de primeira ordem.
c. Calcular h[n].
d. O sistema é estável?

Resolução

a. Equação de diferenças

𝑤[𝑛] = 𝑥 [𝑛] + 0,5𝑦[𝑛] (9)

𝑣[𝑛] = 2𝑥 [𝑛] + 0,5𝑦[𝑛] + 𝑤[𝑛 − 1] (10)

𝑦[𝑛] = 𝑥 [𝑛] + 𝑣[𝑛 − 1] (11)

Substituindo a equação (9) na equação (10):

𝑣[𝑛] = 2𝑥 [𝑛] + 0,5𝑦[𝑛] + 𝑥 [𝑛 − 1] + 0,5𝑦[𝑛 − 1] (12)

Substituindo a equação (12) na equação (11):

𝑦[𝑛] = 𝑥 [𝑛] + 2𝑥 [𝑛 − 1] + 0,5𝑦[𝑛 − 1] + 𝑥 [𝑛 − 2] + 0,5𝑦[𝑛 − 2] (13)

Da equação (13):

Estruturas 8
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

𝒚[𝒏] − 𝟎, 𝟓𝒚[𝒏 − 𝟏] − 𝟎, 𝟓𝒚[𝒏 − 𝟐] = 𝒙[𝒏] + 𝟐𝒙[𝒏 − 𝟏] + 𝒙[𝒏 − 𝟐]

b. Diagrama de fluxo em cascata:

A função de transferência pode ser obtida do gráfico ou da equação de diferenças.

Aplicando transformada z à equação de diferenças:

𝑌 (𝑧)(1 − 0,5𝑧 −1 − 0,5𝑧 −2 ) = 𝑋(𝑧)(1 + 2𝑧 −1 + 𝑧 −2 )

𝑌(𝑧) 1 + 2𝑧 −1 + 𝑧 −2
𝐻 (𝑧 ) = =
𝑋(𝑧) 1 − 0,5𝑧 −1 − 0,5𝑧 −2

Trabalhando com o numerador:

𝑎2 + 2𝑎𝑏 + 𝑏2 = (1 + 𝑎)(1 + 𝑏)

1 + 2𝑧 −1 + 𝑧 −2 = (1 + 𝑧 −1 )2 = (1 + 𝑧 −1 )(1 + 𝑧 −1 )

Trabalhando com o denominador:

1 − 0,5𝑧 −1 − 0,5𝑧 −2 = (1 + 𝑎𝑧 −1 )(1 − 𝑏𝑧 −1 )

1 − 0,5𝑧 −1 − 0,5𝑧 −2 = 1 − 𝑏𝑧 −1 + 𝑎𝑧 −1 − 𝑎𝑏𝑧 −2

Portanto:

𝑎 − 𝑏 = −0,5 𝑎𝑏 = 0,5
𝑎 = 0,5 𝑏=1

1 − 0,5𝑧 −1 − 0,5𝑧 −2 = (1 + 0,5𝑧 −1 )(1 − 1𝑧 −1 ) (14)

Para o filtro em cascata podemos escrever:

𝑌(𝑧) (1 + 𝑧 −1 ) (1 + 𝑧 −1 )
𝐻 (𝑧 ) = = .
𝑋(𝑧) (1 + 0,5𝑧 −1 ) (1 − 1𝑧 −1 )

O diagrama de fluxo do sistema em cascata será:

Estruturas 9
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

c. Para calcular h[n] é conveniente usar divisão longa devido a que o numerador e o
denominador têm o mesmo grau.

1 + 2𝑧 −1 + 𝑧 −2
𝐻 (𝑧 ) =
1 − 0,5𝑧 −1 − 0,5𝑧 −2

A função de transferência fica:

3 + 𝑧 −1
𝐻(𝑧) = −2 +
1 − 0,5𝑧 −1 − 0,5𝑧 −2
Definindo o polinômio para facilitar o cálculo e aplicando frações parciais (usando os polos
calculados na equação (14):

𝑁 3 + 𝑧 −1 3 + 𝑧 −1
𝑃 (𝑧 ) = = =
𝐷 1 − 0,5𝑧 −1 − 0,5𝑧 −2 (1 + 0,5𝑧 −1 )(1 − 1𝑧 −1 )

𝑁 𝐴 𝐵 𝐴(1 − 1𝑧 −1 ) + 𝐵(1 + 0,5𝑧 −1 )

𝑃 (𝑧 ) = = + =
𝐷 1 + 0,5𝑧 −1 1 − 1𝑧 −1 1 − 0,5𝑧 −1 − 0,5𝑧 −2

𝑁 = 𝐴(1 − 1𝑧 −1 ) + 𝐵(1 + 0,5𝑧 −1 ) = 3 + 𝑧 −1

𝐴 − 𝐴𝑧 −1 + 𝐵 + 0,5𝐵𝑧 −1 = 3 + 𝑧 −1

𝐴 + 𝐵 = 3 −𝐴𝑧 −1 + 0,5𝐵𝑧 −1 = 𝑧 −1 ⟹ −𝐴 + 0,5𝐵 = 1 (15)

Do sistema de equações (15):

𝐴 = 0,333 𝐵 = 2,666

𝐻 (𝑧) = −2 + 𝑃(𝑧)
0,333 2,666
𝐻 (𝑧) = −2 + −1
+
1 + 0,5𝑧 1 − 1𝑧 −1
Aplicando transformada z inversa:

𝒉[𝒏] = −𝟐𝜹[𝒏] + 𝟎, 𝟑𝟑𝟑. (−𝟎, 𝟓)𝒏 𝒖[𝒏] + 𝟐, 𝟔𝟔𝟔. 𝟏𝒏 𝒖[𝒏]

d. O sistema tem dois polos, um em 𝑧 = −0,5 e outro em 𝑧 = 1 Devido a este último polo,
o sistema não é estável (polo no círculo de raio unitário).

--------------------------------------------------------------------------------------------------------------------------

Referências
OPPENHEIM, A. V.; SCHAFER, R. W. Digital Signal Processing. New Jersey: Prentice-
Hall, 1975.

Estruturas 10
 Processamento Digital de Sinais

Prof. Eng. Viviana R. Zurro MSc.

OPPENHEIM, A. V.; SCHAFER, R. W. Processamento em Tempo Discreto de Sinais. 3.

ed. São Paulo: Pearson Education do Brasil, 2012.

OPPENHEIM, A. V.; WILLSKY, A. S. Sinais e Sistemas. 2a. ed. São Paulo: Pearson
Education do Brasil, 2010.

Estruturas 11