Sinais de Tempo Discreto e Sistemas
Sinais de Tempo Discreto e Sistemas
Digital de Sinais
Módulo II
Marco Antonio Grivet Mattoso Maia
CETUC-PUC/Rio
Sinal de Tempo Discreto
Do Inglês: Discrete-Time Signals
x[0] x[7]
x[-4] x[10]
x[1]
x[-5] x[6]
x[2] x[3]
x[4]
n
x[5]
n
0
n
0
x[-2]
x[-2] x[0] x[2]
x[-1] x[0]
x[1] x[2] x[3]
x[3]
x[1]
x[-1]
x[9]
x[4]
x[-1]
x[-4]
x[5] x[8]
x[-2]
x[-3] x[6]
y[n] = T { x[n]}
n n
x[5] y[7]
nd = 2
n
y[n] =T { x[n]} =... + x [ n − 3] + x [ n − 2] + x [ n − 1] + x [ n ] =∑ x[k ]
k = −∞
M2
1
y[n] T { x=
= [n]} . ∑ x[n = − k]
M 1 + M 2 + 1 k = − M1
x[n + M 1 ] + x[n + M 1 − 1] + x[n + M 1 − 2] + ..... + x[n − M 2 ]
=
M1 + M 2 + 1
a [a=
1 a2 a3 ....ana +1 ] b [b1 b2 b3 ....bnb +1 ]
Um sistema é dito sem memória quando a saída y[n] para cada valor de n
depende somente de x[n] para o mesmo valor de n.
Exemplos:
y[n] = (x[n])2
y[n] = sen(x[n])
Obs:
O sistema de médias móveis não é um sistema sem memória.
a1 a1
x1[n] x1[n]
x T(.) x
x2[n]
a2 Σ T(.) = x2[n]
a2 Σ
x T(.) x
EXEMPLO:
Um sinal s[n] é dito limitado quando existe um valor fixo positivo e finito Bs
tal que |s[n]|≤Bs< ∞
+∞
y[= ]} T ∑ x[k ].δ [n − k=
n] T { x[n= ]
k = −∞
+∞
∑ x[k ].T {δ [n − k ]}
k = −∞
Logo: +∞
=
y[n] ∑ x[k ].=
k = −∞
h[n − k ] x[n]* h[n]
=h[n] 0
Se o sistema for causal: para n < 0
Logo:
n
=y[n] ∑ x[k ].h[n − k ]
k = −∞
y[0] = x[0].h[0]
=y[1] x[0].h[1] + x[1].h[0]
y[2] = x[0].h[2] + x[1].h[1] + x[2].h[0]
.................
14/07/2023 Curso CIAW-PUC/Rio 24
Tipos de Resposta Impulsional
FIR (Finite Impulse Response)
Se h[n] só é não nulo numa coleção finita
de valores de n
h1[n] h2[n]
h2[n] h1[n]
h1[n]*h2[n]
n
y[n] = ∑ x[k ]
k = −∞
n n −1
z[n=
] y[n] − y[n − 1]= ∑ x[k ] − ∑ x[k =]
k = −∞ k = −∞
x[n]
h1[n]
Σ
h2[n]
h1[n]+h2[n]
x[n] 1 1 y[n] 1
1/2
0 1 0 1 2 3 4 5
-1/2
y[0] 2
y[0] = x[0].h[0] ⇔ h[0] = = = 2
x[0] 1
y[1] − x[1].h[0] 1 − (1).(2)
y[1] = x[0].h[1] + x[1].h[0] ⇔ h[1] = = =
−1
x[0] 1
y[2] = x[0].h[2] + x[1].h[1] + x[2].h[0]
y[2] − x[1].h[1] −0.5 − (1).(−1)
= h[2] = = 0.5
x[0] 1
y[3] = x[0].h[3] + x[1].h[2] + x[2].h[1] + x[3].h[0]
y[3] − x[1].h[2] 0.5 − (1).(0.5)
= h[3] = = 0
x[0] 1
y[4] = x[0].h[4] + x[1].h[3] + x[2].h[2] + x[3].h[1] + x[4].h[0] = 0
14/07/2023 Curso CIAW-PUC/Rio 34
Solução 2 (Inteligente)
x[n]= δ [n] + δ [n − 1] ⇔ y[n]= h[n] + h[n − 1]
y[0]= h[0] + h[−1]= h[0]= 2
y[1]= h[1] + h[0]= h[1] + 2= 1 ⇒ h[1]= −1
= h[2] + h[1]
y[2] = h[2] −= 1 −0.5 ⇒ h[2] = 0.5
y[3] = h[3] + h[2] = h[3] + 0.5 = 0.5 ⇒ h[3] = 0
y[4] = h[4] + h[3] = h[4] = 0
=
y[n] 0 para n ≥ 3
+∞
n
ms ( t ) =m ( t ) . T .δ T ( t ) ⇔ M s ( f ) =M ( f ) ∗ δ1/ T ( f ) =∑ M f −
n = −∞ T
M(f)
f
-2/T -1/T W 1/T 2/T
14/07/2023 Curso CIAW-PUC/Rio 36
14/07/2023 Curso CIAW-PUC/Rio 37
14/07/2023 Curso CIAW-PUC/Rio 38
14/07/2023 Curso CIAW-PUC/Rio 39
Acumulador
n
y[n] = T { x[n]} = ∑ x[k ] ⇔ y[n] =
k = −∞
x[n] + y[n − 1]
=a [1=
1] b [1]
acumulador
x[n] y[n]
Σ
atraso
unitário
1
y[n] − y[n − 1]
= x1[n=
] .{ x[n] − x[n − M 2 + 1]}
M 2 +1
1
=a [1=
1] b .[1 0 ... −1]
M 2 +1
média móvel
1
x[n] M2+1 y[n]
+
x Σ acumulador
-
atraso
nd=M2+1
=
Suponha que x[ n] e j .2.π . f .n para -∞ < n < +∞
Então:
+∞ +∞
= =
y[n] x[n]* h[n] ∑
h = −∞
h[k ].x[=
n − k] ∑
k = −∞
.2. . f .( n − k )
h[k ].e j π=
+∞
. ∑ h[k ].e − j .2.π . f .k e j .2.π . f .n .H ( f )
j .2.π . f .n
e=
k = −∞
H(f)
resposta em freqüência
(função complexa)
arg ( H ( f ) )
H ( f ) = H ( f ) .e
A resposta em freqüência de um slit T pode ser conhecida, observando
a sua saída quando a entrada é da forma x[n]=exp(j.2.π.f.n)
− j .2.π . f .nd
H ( f ) = 1
=
Logo H( f ) e ⇒
arg [ H ( f ) ] = −2.π . f .nd
Se x[n] = e j .2.π . f .n
n n
=
Então y[n] ∑
= e e
k = −∞
j .2.π . f .k j .2.π . f .n
.∑ =
e j .2.π . f .( k − n )
k = −∞
∞
1
.∑ e − j .2.π . f .k e j .2.π . f .n .
e=j .2.π . f .n
k =0 1 − e − j .2.π . f
1
Logo H ( f ) =
1 − e − j .2.π . f
+∞ +∞
H( f ) ∑
= h[k ].e
k = −∞
e − j .2.π . f .k − j .2.π .k
∑
= h[ k
k = −∞
].e − j .2.π . f .k
+∞
= ∑
k = −∞
h[k ].e − j .2.π .( f=
+1).k
H ( f + 1)
+∞
X ( f ) = ∑ x[k ].e − j .2.π . f .k
k = −∞
+∞ ⇔ x[n] ↔ X ( f )
x[n] = ∫ X ( f ).e j .2.π . f .n
.df
−∞
x[−n] X (− f )
n.x[n] j dX ( f )
.
2.π df
T T
C/D
-W W f T.δT(t)
f
0
2.W
T.δT(t)
-W W f
<1/T filtro passa-baixa
1/T
f
0
xa(t)
-3T -2T -T 0 T 2T 3T 4T t
D/C
xa[0] xa[7]
xa[-4] xa[10]
xa[1]
xa[-5] xa[6]
xa[2] xa[3]
xa[4]
n
xa[5]
xa[0] xa[7]
xa[-4] xa[10]
xa[1]
xa[-5] xa[6]
xa[2] xa[3]
xa[4]
n
xa[5]
ℑ u ( t − n.=
T ) e − j .2π . f .( n.T ) .ℑ u =
( t ) e− j.2π . f .( n.T ) .U ( f )
{ }
+∞
X DA ( f ) = ∑ xa [ n ]. e − j .2π . f .( n.T ) − e − j .2π . f .( ( n +1).T ) .U ( f ) =
n = −∞
+∞
= (1 − e − j .2π . f .T
) .U ( f ) . ∑ x [ n].e
a
− j .2π . f .( n.T )
n = −∞
n = −∞
= (1 − e − j .2π . f
) .U ( f ) .X ( f )
Logo
1
(1 − e
H(f )= − j .2π . f
) .U ( f ) =
(1 − e − j .2π . f
) . 0.5δ ( f ) + j.2.π . f =
1 − e − j .2π . f
=
j.2.π . f
zero-order holder
Esta soma não converge de forma geral. É comum que ela só convirja
para determinados valores de z. Os valores para os quais a soma
converge formam o que é denominado de Região de Convergência (ROC)
x[n] = δ [n] Im
+∞
( z)
X= ∑
n = −∞
z − n 1.=
x[n].= z −0 1
+∞ n
1 1 1
( z ) ∑ =
Re
X= = −1
n =0 z 1−
1 1 − z
z
ROC= { z ∈ | z > 1}
x[n] = δ [n] − δ [n − 1] Im
+∞
∑
X ( z) = x[ n ].
n = −∞
z −n
=
1. z −0
− 1. z −1
=
1 − z −1
X ( z) ∑
= = x[n].z ∑= A.α .z
−n
A.∑
n −n
n = −∞ n= 0 n= 0 z
α |α|
Esta soma só converge sse < 1 e neste caso
z Re
1 z 1
X ( z ) A=
= . A=
. A.
α z − α 1 − α . z −1
1−
z
para z > α ROC= { z ∈ | z > α }
α β
+∞ +∞ n n
A.∑ +B.∑
z
n 0= n 0 z
A 1a. soma só converge sse z > α e a 2a. soma só converge sse z > β
Então as duas somas só existem sse z > max { α , β }
Neste caso:
X ( z ) A.
=
1
1 − α .z −1
+ B.
1
1 − β .z −1
{ }
ROC= z ∈ | z > max { α , β }
Método de Inspeção
Expansão por Frações Parciais
Expansão por Série de Potencias
Integral de Contorno
x[n − nd ] X ( z ).z − nd
z0n .x[n] X ( z / z0 )
x[−n] X (1/ z )
n.x[n] dX ( z )
− z.
dz
Seqüência Transformada Z
x[n], y[n] X(z),Y(z)
=
k
k 0= k 0
∑ a . y[n=
− k ] ∑ b .x[n − k ] k
na nb
=
k
k 0= k 0
∑ a .z −k
Y ( z ) = ∑ bk .z − k X ( z )
na −k nb −k
∑ k a . z .Y ( z ) = ∑ k b . z .X ( z)
= k 0= k 0
nb
Y ( z) ∑ b .z k
−k
( z)
H= = k =0
na
X ( z)
∑ a .z
k =0
k
−k
Y ( z) ∑k
b . z −k
( z)
H= = k =0
na
X ( z)
∑ k
a
k =0
. z −k
n = −∞ n= 0 1 − α .z −1
ROC = α
H ( z ) = 0 ⇒ 1 − α .z −1 = ∞ ⇒ z = 0 ⇒ não há zeros devido ao ROC
H ( z ) = ∞ ⇒ 1 − α .z −1 = 0 ⇒ z = α
x[n] x[n-1]
z-1 atraso unitário
z-1
a1 y[n-1]
+
z-1
a2
y[n-2]
y[=
n] a1. y[n − 1] + a2 . y[n − 2] + b0 .x[n] + b1.x[n − 1]
) a1.z −1.Y ( z ) + a2 .z −2 .Y ( z ) + b0 . X ( z ) + b1.z −1. X ( z )
Y ( z=
Y ( z) b0 + b1.z −1
=
H ( z) =
X ( z ) 1 − a1.z −1 − a2 .z −2
x[n] y[n]
b0 w[n]
+ +
z-1 z-1
b1 a1
+ y[n-1]
x[n-1]
z-1
a2
y[n-2]
a0 . y=
[n] b0 .x[n] + b1.x[n − 1] + b2 .x[n − 2] + ..... + bN .x[n − N ] +
- a1. y[n − 1] - a2 . y[n − 2] -.....- aM . y[n − M ]
N
Y ( z) ∑
−k
b . z −1 −2 −3 −N
k
b + b . z + b . z + b . z .... + b . z
( z ) = kM=0 = 0 1 −1 2 −2 3 −3
H= N
−M
+ + + +
∑ ak .z − k 0 1
X ( z) a a . z a2 . z a3 . z .... a M . z
k =0
z-1 z-1
a1
x[n-1] b1 + y[n-1]
+
z-1 z-1
a2
b2 + y[n-2]
x[n-2] +
bN-1 aM-1
x[n-N-1] + + y[n-M-1]
z-1 z-1
bN aM
x[n-N] y[n-M]