HỆ RỜI RẠC TUYẾN TÍNH BẤT BIẾN TRONG MIỀN TẦN SỐ... Nếu h[n] là đáp ứng xung của hệ rời rạc LTI trong miền thời gian thì Hej là đáp ứng tần số bằng cách lấy biến đổi Furie rời rạc của h
Trang 1HỆ RỜI RẠC TUYẾN TÍNH BẤT BIẾN TRONG MIỀN TẦN SỐ
Trang 2Nếu h[n] là đáp ứng xung của hệ rời rạc LTI trong miền thời gian thì H(ej) là đáp ứng tần số bằng cách lấy biến đổi Furie rời rạc của h[n]:
H
R.4.1
Trang 3R.4.2 H(ej) là hàm phức của với chu kỳ 2 và có thể được biểu diễn bằng phần thực và phần ảo hoặc theo độ lớn và pha
Do đó:
) (
im
re ( ) H ( ) H ( ) e
H ) ( ej ej j e j e j j
r
)
) (e j
H : đáp ứng biên độ của hệ LTI
Trang 4
R.4.3 Hàm tăng ích G() của hệ LTI được định nghĩa như sau:
A() = -G() : hệ số suy hao của hệ LTI
) (
log 20
)
Trang 5R.4.4 Đối với hệ thống rời rạc được đặc trưng bởi đáp ứng xung thực h[n], hàm biên độ là một hàm chẵn và pha là hàm
lẻ theo , nghĩa là:
) (
) (
) (
) (
Tương tự, ( j ), H im ( j )
re e e
theo
Trang 6R.4.6 Thời gian truyền nhóm của hệ rời rạc tuyến tính bất biến theo thời gian được định nghĩa như sau:
Trong đó: c() là hàm trãi pha Nếu pha được tính bằng
radian thì trễ nhóm được tính bằng s
biến hệ số thực có đáp ứng tần số H(e j) và tín hiệu vào
x[n] = A cos (0n + ) , với A là hệ số thực được cho như sau:
) ) ( ( cos ) (
] [ n A H ej n
y
Trang 7R.4.8 Từ công thức tổng chập biểu diễn cho hệ rời rạc tuyến tính bất biến theo thời gian, ta suy ra đáp ứng tần số của hệ rời rạc LTI được tính bằng tỉ số giữa Y(ej) và X(ej)
) (
)
( )
e X
e
Y e
Trong đó : Y(ej) là biến đổi Furie của tín hiệu ra y[n]
X(ej) là biến đổi Furie của tín hiệu vào x[n]
Trang 8R.4.9 Đối với hệ thống tuyến tính bất biến theo thời gian được đặc trưng bởi phương trình sai phân hệ số hằng thì đáp ứng tần số H(ej) được biểu diễn dưới dạng:
M
k
k
j k
j
e d
e
p e
Trang 9R.4.10.Biến đổi z, H(z) của chuỗi đáp ứng xung {h[n] }được gọi là hàm truyền đạt hoặc hàm hệ thống Từ công thức tổng chập biểu diễn cho hệ thống tuyến tính bất biến theo thời gian ta suy ra hàm truyền H(z) chính là tỉ số giữa Y(z) (biến đổi z của đầu ra y[n] ) và X(z), (biến đổi z của tín hiệu vào x[n])
) (
)
( )
(
z X
z
Y z
Trang 10R.4.11 Nếu miền hội tụ (ROC) của H(z) có chứa đường tròn đơn vị thì đáp ứng tần số của hệ rời rạc LTI được tính như sau:
j
e z
j j
j j
( ) ( )
( ) (
) (
Trang 11R.4.13 Một hệ rời rạc LTI được đặc trưng bởi phương trình sai phân tuyến tính hệ số hằng thì hàm truyền H(z) được biểu diễn:
N N
M M
z d z
d d
z p z
p
p z
X
z
Y z
)
( )
1 0
1 1 0
) 1
(
) 1
( )
(
1 1
0
1 1
z
p z
H
k
N k
k
M k
Trong đó: 1, 2, …., M là các điểm không
1, 2, …., N là các điểm cực
Nếu N > M thì thêm (N-M) điểm không tại z = 0
Nếu N < M thì thêm (M-N) điểm cực tại z = 0
Trang 12
R.4.15 Tất cả các điểm cực của hàm truyền nhân quả ổn định đều phải nằm hoàn toàn bên trong đường tròn đơn vị.
trong dải thông , bằng 0 trong dải chắn và có pha 0 bất kỳ Đáp ứng tần số của bốn bộ lọc số pha không lý tưởng với đáp ứng xung hệ số thực được biểu diễn như hình vẽ
Trang 13HBS(e j )
1
Trang 14R.4.17 Đáp ứng xung hLP[n] của bộ lọc thông thấp lý tưởng có dạng:
[
n
n n
h c LP
1)
H LP
Trong đó: <1 thì hệ ổn định Tần số c gọi là tần số cắt 3dB và liên hệ với tham số qua công thức:
c
Trang 15Hàm truyền của bộ lọc thông cao bậc 1 có đáp ứng xung dài
vô hạn có dạng:
1
11
1 2
1 )
Trang 16Đáp ứng biên độ sẽ bằng 0 tại = 0 và = Giá trị cực đại bằng 1 ứng với = 0, 0 được gọi là tần số giữa của bộ lọc thông dải Khi đó:
) (
Băng thông 3-dB (3dB) là độ sai khác giữa hai tần số cắt 3dB
và được tính như sau:
2
dài vô hạn có dạng:
2 - 1
-2
z
z ) (1
1
1
2
1 )
HBP
Trang 17R.4.21 Hàm truyền của bộ lọc chắn dải bậc hai có đáp ứng xung dài vô hạn được tính bằng công thức:
2 1
2 -1
)1
(1
z2
1
2
1)
z z
Trang 18R.4.22.Với K bộ lọc thông thấp bậc 1 mắc nối tiếp thì hàm truyền của hệ thống sẽ là:
K LP
z z
1
z1
1.2
1)
c
cos1
2sin
cos)1
Trong đó B 2K1/ K
Trang 19x[n] H(z) v[n] u[n] H(z)
w[n]
u[n] = v[-n], y[n] = w[-n]
với chiều dài hữu hạn không theo thời gian thực một cách dễ dàng nếu điều kiện nhân quả được relaxed Trong sơ đồ sau, tín hiệu vào hữu hạn được đưa qua bộ lọc nhân quả hệ số thực có hàm truyền là H(z), tín hiệu ra sau đó được lấy nghịch đảo và được đưa qua bộ lọc thêm một lần nữa Quá trình này được mô
tả như sau:
Trang 20R.4.24 Chúng ta có thể thiết kế hàm truyền với đáp ứng xung
hữu hạn (FIR) và đáp ứng pha tuyến tính Nếu hàm truyền có
đáp ứng xung đối xứng thì:
N n
0 ,
] [
] [ n h N n
h
hoặc đáp ứng xung phản đối xứng thì:
N n 0 , ] [
] [ n h N n
h
N : bậc của hàm truyền và h[n] có chiều dài là N + 1 Có 4 loại hàm truyền:
Loại 1: Đáp ứng xung đối xứng với chiều dài lẻ
Loại 2: Đáp ứng xung đối xứng với chiều dài chẵn
Loại 3: Đáp ứng xung phản đối xứng với chiều dài lẻ
Trang 21R.4.25 Hàm truyền FIR loại 2 có một không tại z = -1, do đó không thể dùng để thiết kế bộ lọc thông cao Hàm truyền FIR loại 3 có không tại z = 1 và z = -1 nên cũng không thể dùng để thiết kế bộ lọc thông thấp, thông cao và chắn dải Bộ lọc FIR loại 4 cũng không thể dùng để thiết kế lọc thông thấp do xuất hiện không tại z = 1 Chỉ có bộ lọc FIR loại 1 mới được sử dụng để thiết kế các loại bộ lọc.
Trang 22R.4.26. Hàm truyền H(z) hệ số thực gọi là nhân quả và ổn định nếu:
) 1, (e j
(e j 2A
thì A(z) được gọi là hàm truyền thông tất Khi đó hàm truyền thông tất IIR hệ số thực bậc M nhân quả có dạng:
) (
)
( )
(
1
z D
z D
z z
Trang 23R.4.28 Hàm truyền nhân quả và ổn định với tất cả các điểm không nằm bên trong hoặc trên đường tròn đơn vị gọi là hàm truyền pha cực tiểu Ngược lại hàm truyền nhân quả và ổn định với tất cả các điểm không nằm ngoài đường tròn đơn vị gọi là hàm truyền pha cực đại.
được gọi là bù trễ lẫn nhau nếu tổng các hàm truyền của chúng bằng một bội số nguyên lần của trễ đơn vị, nghĩa là:
0
, )
( 1
Trong đó n0 là số nguyên không âm
Trang 24Hàm truyền bù trễ H1(z) và hàm truyền FIR pha tuyến tính loại 1 H0(z) có chiều dài L lẻ liên hệ với nhau qua biểu thức:
)()
( 1
1
z A z
M
j
e H
Trang 25R.4.32 Để Am(z) là hàm thông tất hệ số thực bậc m thì:
m m
m m
m m
m m
m m
z d z
d z
d z
d
z z
d z
d z
d
d z
1 1
) 1 ( 1
2 2
1 1
1
) (
Suy ra hàm thông tất hệ số thực bậc m-1 có dạng:
) 1 ( '
1
) 2 ( '
2
1 ' 1
) 1 ( )
2 ( ' 1
1 '
2
' 1 - m 1
1
d
) ( 1
) (
m m
m m
m
m m
m
m m
z d
z d
z d
z z
d z
d
z A
k
k z
A z
z A
Trong đó:
1-m1,2, ,i
,
i i
d
d d d
d
Trang 26Đặt km = Am() Để Am(z) ổn định thì km2 < 1 và Am-1(z) là hàm thông tất ổn định Để kiểm tra tính ổn định của Am-1(z) ,
ta tạo hàm thông tất bậc m-2 và cứ tiếp tục ta sẽ có một tập hàm thông tất bậc giảm dần:
Am(z), Am-1(z),…, A2(z), A1(z), A0(z) = 1
Với các hệ số : km, km-1,……, k2 , k1
Hàm thông tất Am(z) ổn định khi và chỉ khi kl2 < 1,
với l = m, m-l,… ,l