Bài giảng Xử lý tín hiệu số: Chương 2 - TS. Đặng Quang Hiếu

Chia sẻ: _ _ | Ngày: | Loại File: PDF | Số trang:15

Thêm vào BST

Báo xấu

29
lượt xem 4
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài giảng "Xử lý tín hiệu số - Chương 2: Các phép biến đổi Fourier" cung cấp cho người học các kiến thức: Biến đổi Fourier, chuỗi Fourier rời rạc cho dãy tuần hoàn, biến đổi Fourier rời rạc. Mời các bạn cùng tham khảo nội dung chi tiết.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Bài giảng Xử lý tín hiệu số: Chương 2 - TS. Đặng Quang Hiếu

ET4020 - Xử lý tín hiệu số Chương 2: Các phép biến đổi Fourier TS. Đặng Quang Hiếu http://dsp.edabk.org Trường Đại học Bách Khoa Hà Nội Viện Điện tử - Viễn thông Năm học 2012 - 2013 Outline Biến đổi Fourier Chuỗi Fourier rời rạc cho dãy tuần hoàn Biến đổi Fourier rời rạc
Biến đổi Fourier FT n ω IFT ∞ X FT jω x(n) −−→ X (e ) = FT{x(n)} = x(n)e −jωn x=−∞ ◮ Tuần hoàn với chu kỳ 2π ◮ Phổ biên độ: |X (e jω )|, và phổ pha: arg{X (e jω )}. ◮ Biến đổi ngược: Z π IFT 1 X (e ) −−−→ x(n) = IFT{X (e jω )} = jω X (e jω )e jωn dω 2π −π Các ví dụ về FT 1. Tìm X (e jω ), |X (e jω )| và arg{X (e jω )} của các dãy sau đây: (a) x(n) = δ(n) (b) x(n) = δ(n − 2) (c) x(n) = δ(n − 2) − δ(n) (d) x(n) = rectN (n) (e) x(n) = (0.5)n u(n) (f) x(n) = u(n) 2. Xét bộ lọc thông thấp lý tưởng có đáp ứng tần số (trong một chu kỳ) như sau: 1, |ω| ≤ ωc Hlp (e jω ) = 0, ωc < |ω| ≤ π (a) Hãy tìm đáp ứng xung hlp (n) của bộ lọc này. (b) Giải bài toán cho trường hợp bộ lọc thông cao
Phổ biên độ và phổ pha của rect10 (n) 10 8 6 |X(jω)| 4 2 0 −8 −6 −4 −2 0 2 4 6 8 ω 4 2 arg{X(jω)} 0 −2 −4 −8 −6 −4 −2 0 2 4 6 8 ω Các tính chất ◮ Quan hệ với biến đổi z: X (e jω ) = X (z)|z=e jω ◮ Điều kiện hội tụ: ∞ X |x(n)| < ∞ n=−∞ Một hệ thống LTI có đáp ứng tần số khi và chỉ khi nó ổn định. ◮ Tuyến tính, dịch thời gian, dịch tần số, chập, v.v. ◮ Các tính chất đối xứng ◮ Quan hệ Parseval X∞ Z π 2 1 |x(n)| = |X (e jω )|2 dω n=−∞ 2π −π ◮ Định lý Wiener - Khintchine: Nếu x(n) ∈ R thì FT{rxx (n)} = SXX (e jω ) := |X (e jω )|2 trong đó SXX (e jω ) là phổ mật độ năng lượng của x(n).
Outline Biến đổi Fourier Chuỗi Fourier rời rạc cho dãy tuần hoàn Biến đổi Fourier rời rạc Khái niệm dãy tuần hoàn x˜(n) = x˜(n − N), ∀n ◮ Chu kỳ N ∈ Z → ký hiệu x˜(n)N . ◮ Tồn tại khai triển Fourier ◮ Khác hệ số N so với khái niệm chuỗi Fourier cho tín hiệu tuần hoàn trong môn Tín hiệu và hệ thống!
Định nghĩa cặp chuỗi Fourier rời rạc cho dãy tuần hoàn N−1 X 2π ˜ (k) = X x˜(n)e −j N kn n=0 N−1 1 X˜ 2π x˜(n) = X (k)e j N kn N k=0 2π ◮ WN = e −j N . ◮ ˜ (k)|, arg{X Biên độ và pha: |X ˜ (k)}. Ví dụ: Cho tín hiệu tuần hoàn x˜(n) với chu kỳ N: 1, ℓN ≤ n ≤ ℓN + M − 1, ∀n ∈ Z, M < N x˜(n) = 0, n còn lại ˜ (k), |X Hãy tìm X ˜ (k)|, arg{X ˜ (k)}. Khi N = 100, M = 10 0.1 0.09 0.08 0.07 0.06 |X[k]| 0.05 0.04 0.03 0.02 0.01 0 −100 −80 −60 −40 −20 0 20 40 60 80 100 k 3 2 1 arg{X[k]} 0 −1 −2 −3 −100 −80 −60 −40 −20 0 20 40 60 80 100 k
Các tính chất ◮ Tuyến tính, dịch thời gian, dịch tần số ◮ Đối ngẫu: Nếu DFS ˜ x˜(n) ←−−→ X (k) thì DFS ˜ (n) ←− X −→ N x˜(−k) ◮ Các tính chất đối xứng Chập tuần hoàn DFS ˜ x˜1 (n) ←−−→ X 1 (k) DFS ˜ x˜2 (n) ←−−→ X 2 (k) ˜3 (k) = X Nếu X ˜1 (k)X ˜2 (k) −→ Chập tuần hoàn: N−1 X ∗)N x˜2 (n) = x˜3 (n)N = x˜1 (n)(˜ x2 (n − m) x˜1 (m)˜ m=0
Các bước tính chập tuần hoàn Tìm x˜3 (n0 ), ∀n0 ∈ [0, (N − 1)] (1) Lấy đối xứng x˜2 (m) → x˜2 (−m) (2) Dịch theo trục thời gian đi n0 mẫu x2 (n0 − m) trong đoạn [0, (N − 1)] (3) Nhân: v˜n0 (m) = x˜1 (m)˜ (4) Tính tổng: Cộng tất cả thành phần khác không của v˜n0 (m) trong đoạn [0, (N − 1)] → x˜3 (n0 ) (5) Kết quả là một dãy tuần hoàn với chu kỳ N: x˜3 (n0 ) = x˜3 (n0 + rn), ∀r ∈ Z. Minh họa các bước tính phép chập tuần hoàn x˜1 (m) x˜2 (m) b b b b b b b b b b b b b b b b b b b b b b m m -4 -3 -2 -1 0 1 2 3 4 5 6 -4 -3 -2 -1 0 1 2 3 4 5 6 x˜2 (−m) v˜0 (m) x˜3 (0) = 1.75 b b b b b b b b b b b b b b b b b b b b m b b m -4 -3 -2 -1 0 1 2 3 4 5 6 -4 -3 -2 -1 0 1 2 3 4 5 6 x˜2 (1 − m) v˜1 (m) x˜3 (1) = 2 b b b b b b b b b b b b b b b b b b b b m b b m -4 -3 -2 -1 0 1 2 3 4 5 6 -4 -3 -2 -1 0 1 2 3 4 5 6 x˜2 (2 − m) v˜2 (m) x˜3 (2) = 2.25 b b b b b b b b b b b b b b b b b b b b m b b m -4 -3 -2 -1 0 1 2 3 4 5 6 -4 -3 -2 -1 0 1 2 3 4 5 6
Kết quả phép chập tuần hoàn x˜1 (n) b b b b b b b b b b n b b b -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 x˜2 (n) b b b b b b b b b b b b b n -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 x˜3 (n) b b b b b b b b b b b b b n -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 Bài tập 1. Viết chương trình Matlab để vẽ phổ biên độ và phổ pha của một dãy có chiều dài hữu hạn bất kỳ 2. Sử dụng hàm freqz trong Matlab để vẽ đáp ứng tần số của một hệ thống LTI từ phương trình sai phân tuyến tính hệ số hằng. 3. Lấy mẫu tần số. Cho dãy x(n) có chiều dài hữu hạn L với phổ X (e jω ) (chu kỳ 2π). Để biểu diễn phổ tín hiệu, người ta lấy jk 2π các mẫu tại tần số ω = k 2π N để thu được X (e N ) với chu kỳ lấy mẫu 2π N . Với những giá trị nào của N thì ta có thể tái tạo 2π lại hoàn toàn x(n) từ các mẫu X (e jk N )?
Outline Biến đổi Fourier Chuỗi Fourier rời rạc cho dãy tuần hoàn Biến đổi Fourier rời rạc Khái niệm Xét tín hiệu x(n) có chiều dài hữu hạn N, nếu lấy đủ mẫu (tối thiểu N / một chu kỳ) của phổ X (e jω ), thì có thể khôi phục lại được x(n). −→ Biến đổi Fourier rời rạc DFT cho dãy có chiều dài hữu hạn! Cho x(n) với chiều dài hữu hạn N: x(n) = 0, ∀n < 0, n > N − 1, ta có dãy tuần hoàn x˜(n): x˜(n) = x(n mod N) Lấy một chu kỳ từ DFS{X˜ (k)}: X˜ (k), 0 ≤ k ≤ (N − 1) X (k) = 0, k còn lại
Định nghĩa cặp biến đổi Fourier rời rạc N−1 X 2π X (k) = DFT{x(n)} = x(n)e −j N kn , ∀k ∈ [0, N − 1] n=0 N−1 1 X 2π x(n) = IDFT{X (k)} = X (k)e j N kn , ∀n ∈ [0, N − 1] N k=0 x(0) X (0) x(1) DFT X (1) .. .. . N - điểm . x(N − 1) X (N − 1) Ví dụ: Tìm DFT N-điểm của x(n) = rectM (n) cho ba trường hợp: M = 1, M = N và 1 < M < N. Dạng ma trận Xét ma trận WN×N trong đó Wkn = WNkn   1 1 1 ··· 1  1 WN1 WN2 · · · WN (N−1)     2(N−1)  W= 1 WN2 WN4 · · · WN   . .. .. ..   . ..   . . . . .  (N−1) (N−1)2 (N−1)2 1 WN WN · · · WN và X = [X (0), X (1), · · · , X (N − 1)]T x = [x(0), x(1), · · · , x(N − 1)]T DFT và IDFT có thể được biểu diễn dưới dạng: X = Wx 1 H x = W X N
Dịch vòng: Một chu kỳ của tín hiệu tuần hoàn sau dịch x(n) b b b b n -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 x˜(n)N b b b b b b b b b b b b b n -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 x˜(n − n0 )N b b b b b b b b b b b b b n -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 x(n − n0 )N b b b b n -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 Dịch vòng: Đặt lên một vòng tròn và quay quanh tâm b a c x(n)N = {a, b, c, d} a b x(n − 1)N = {d, a, b, c} d d c
Tính chất dịch ◮ Dịch thời gian DFT{x(n − n0 )N } = e −j(2π/N)kn0 X (k) ◮ Dịch tần số DFT{e j(2π/N)k0 n x(n)} = X (k − k0 )N Đối ngẫu Nếu DFT{x(n)} = X (k) thì DFT{X (n)} = Nx(−k)N Lưu ý: x(−k)N =?
Đảo trục thời gian Nếu DFT{x(n)} = X (k) thì DFT{x(−n)N } = X (−k)N Các tính chất đối xứng (a) DFT{x ∗ (n)} = X ∗ (−k)N (b) DFT{x ∗ (−n)N } = X ∗ (k) (c) DFT{Re[x(n)]} = 21 [X (k) + X ∗ (−k)N ] (d) DFT{ 21 [x(n) + x ∗ (−n)N ]} = Re[X (k)] (e) Nếu x(n) ∈ R ◮ X (k) = X ∗ (−k)N = X ∗ (N − k) ◮ Re[X (k)] = Re[X (N − k)] ◮ Im[X (k)] = −Im[X (N − k)] ◮ |X (k)| = |X (N − k)| ◮ arg{X (k)} = − arg{X (N − k)}
Chập vòng Định nghĩa chập vòng: N−1 X x3 (n)N = x1 (n)(∗)N x2 (n) = x1 (m)x2 (n−m)N , ∀n ∈ [0, N −1] m=0 Áp dụng DFT ta có: DFT{x1 (n)(∗)N x2 (n)} = X1 (k)X2 (k) Cách tính chập vòng: ◮ Miền thời gian ◮ Miền tần số Ví dụ: Tính chập vòng 5-điểm (N = 5) của hai dãy sau: x1 (n) = rect4 (n) + 0.5δ(n − 4) 1 − n4 , 0 ≤ n ≤ 4 x2 (n) = 0, n còn lại Dạng ma trận của chập vòng x3 = X2 · x1 trong đó x3 = [x3 (0), x3 (1), · · · , x3 (N − 1)]T , x1 = [x1 (0), x1 (1), · · · , x1 (N − 1)]T và X2 là (circulant matrix):   x2 (0) x2 (N − 1) · · · x2 (1)  x2 (1) x2 (0) · · · x2 (2)    X2 =  .. .. . . ..   . . . .  x2 (N − 1) x2 (N − 2) · · · x2 (0) ◮ Dạng ma trận của chập tuyến tính → ma trận Toeplitz! ◮ Làm thế nào để tính chập vòng bằng Matlab?
Mối quan hệ giữa chập vòng và chập tuyến tính Cho hai dãy có chiều dài hữu hạn, x(n): [0 · · · (N − 1)] và h(n): [0 · · · (M − 1)]. Nếu y1 (n) = x(n) ∗ h(n) và y2 (n) = x(n)(∗)L h(n) (a) Với những giá trị nào của L thì y1 (n) = y2 (n), ∀n? (b) Nếu L = N thì tại những thời điểm n nào ta có y1 (n) = y2 (n)? Quan hệ Parseval N−1 X N−1 ∗ 1 X x(n)y (n) = X (k)Y ∗ (k) N n=0 k=0 Nếu x(n) = y (n): N−1 X N−1 2 1 X |x(n)| = |X (k)|2 N n=0 k=0