Lấy mẫu tín hiệu ( chương 7)

Chia sẻ: Dang Anh Tu | Ngày: | Loại File: PDF | Số trang:37

Thêm vào BST

Báo xấu

347
lượt xem 45
download

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

Lấy mẩu Tín hiệu liên tục có thể được xử lý bằng cách xử lý các mẩu của tín hiệu qua hệ thống rời rạc. Điều cần thiết là phải duy trì tốc độ lấy mẩu tín hiệu đủ lớn để khôi phục tốt tín hiệu (không có sai số hay sai số với độ dung sai chấp nhận được). Điều này có thể thực hiện được dùng định lý lấy mẩu. 5.1 Định lý lấy mẫu Ta sẽ chứng minh là tín hiệu thực có phổ băng thông giới hạn B Hz [F(w)= 0 với w 2pB ]...

Chủ đề:

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

Lưu

Nội dung Text: Lấy mẫu tín hiệu ( chương 7)

CHƯƠNG 7: LẤY MẨU TÍN HIỆU Nội dung Định lý lấy mẩu 5.1 Tính biến đổi Fourier: Biến đổi Fourier rời rạc (DFT) 5.2 Biến đổi Fourier nhanh (FFT) 5.3 Phụ chương 5.1 5.4 Tóm tắt 5.5 Tài liệu tham khảo: B.P. Lathi, Signal Processing and Linear Systems, Berkeley-Cambridge Press, 1998 Lấy mẩu Tín hiệu liên tục có thể được xử lý bằng cách xử lý các mẩu của tín hiệu qua hệ thống rời rạc. Điều cần thiết là phải duy trì tốc độ lấy mẩu tín hiệu đủ lớn để khôi phục tốt tín hiệu (không có sai số hay sai số với độ dung sai chấp nhận được). Điều này có thể thực hiện được dùng định lý lấy mẩu. 5.1 Định lý lấy mẫu Ta sẽ chứng minh là tín hiệu thực có phổ băng thông giới hạn B Hz [F(w)= 0 với w > 2pB ] có thể được khôi phục chính xác (không có sai số nào) từ các tốc độ lấy mẩu đồng đều với tốc độ Fs > 2B mẩu/giây. Nói cách khác, tốc dộ lấy mẩu tối thiểu là Fs = 2B Hz. Để chứng minh định lý lấy mẩu, xét tín hiệu f(t) (hình 5.1a) có phổ giới hạn B Hz (hình 5.1b). Để thuận tiện, ta vẽ phổ là hàm theo w cũng như theo F (Hz). Lấy mẩu f(t) với tốc độ Fs Hz (Fs mẩu/giây) có thể thực hiện bằng cách nhân f(t) với chuỗi xung dT(t)
(hình 5.1c), gồm các xung đơn vị lặp lại theo chu kỳ T giây, với T =1/Fs. Kết quả là tín hiệu được lấy mẫu f (t ) vẽ trong hình 5.1d, là tín hiệu gồm các xung cách nhau từng T giây (thời gian lấy mẫu). Xung thứ n, nằm tại t = nT, có cường độ f(nT), là giá trị của f(t) tại t = nT. f (t ) = f (t )d T (t ) = å f (nT )d (t - nT ) (5.1) n Để tìm F (w ) , biến đổi Fourier của f (t ) , ta lấy biến đổi Fourier của vế phải phương trình (5.3) từng thừa số một. Biến đổi của thừa số thứ nhất trong ngoặc là F(w). Biến đổi của thừa số thứ hai 2 f (t ) cos w s t là F(w –ws ) + F(w +ws ) (xem phương trình (4.41), cho thấy phổ F(w) dời ws và –ws. Tương tự, biến đổi Fourier của thừa số thứ ba 2 f (t ) cos 2w s t là F(w –2ws ) + F(w +2ws ), cho thấy phổ F(w) dời 2ws và –2ws, và tiếp tục cho tới vô hạn. Điều này tức là phổ F (w ) gồm F(w) lặp lại theo chu kỳ ws = 2p/T rad/s, hay Fs = 1/T Hz, như vẽ trong hình 5.1e. Ngoài ra còn có thêm hằng số nhân 1/T trong phương trình (5.3). Do đó 1¥ F (w ) = å F (w - nw s ) (5.4) T n=-¥ Nếu muốn khôi phục f(t) từ f (t ) , ta phải khôi phục được F(w) từ F (w ) . Có thể khôi phục được nếu không có trùng lắp giữa các chu kỳ liên tiếp của F (w ) . Hình 5.1e cho thấy cần có Fs > 2B (5.5) Đồng thời, thời gian lấy mẩu T =1/Fs. Do đó 1 T£ (5.6) 2B ¢ Thí dụ 5.1 Trong thí dụ này, ta xét ảnh hưởng của tín hiệu khi lấy mẩu theo tốc độ Nyquist, thấp hơn tốc độ Nyquist (lấy mẩu thiếu), cao hơn tốc độ Nyquist (lấy mẩu lố). Xét tín hiệu f (t ) = sin c 2 (5pt ) (hình 5.2a) có phổ F (w ) = 0,2D( 20p ) (hình 5.2b). Băng thông của w tín hiệu là 5 Hz (10p rad/s). Như thế, tốc độ Nyquist là 10Hz; tức là, ta phải lấy mẩu tín hiệu với tốc độ không nhỏ hơn 10 mẩu/s. Khoảng Nyquist là T = 1/2B = 0,1 giây. Nhắc lại là phổ của tín hiệu đã lấy mẩu gồm (1 / T )F (w ) = T D( 20p ) lặp lại theo chu w 2 kỳ bằng với tần số lấy mẩu Fs Hz. Ta trình bày thông tin này trong bảng sau với tốc độ lấy mẩu Fs = 5Hz (lấy mẩu thiếu). 10Hz (tốc độ Nyquist) và 20Hz (lấy mẩu lố). (1/T)F(w) Thời gian lấy mẩu T Nhận xét Tần số lấy mẩu Fs D( 20p ) w Lấy mẩu thiếu 5 Hz 0,2 2D( 20p ) w Tốc độ Nyquist 10 Hz 0,1 4D ( ) w Lấy mẩu lố 20 Hz 0,05 20p Trong trường hợp đầu (lấy mẩu thiếu), tốc độ lấy mẩu là 5Hz (5 mẩu/giây) , và phổ F (w ) lặp lại sau mỗi 5Hz (10p rad/s). Các phổ liên tiếp trùng lắp, như vẽ trong hình 1 T
5.2d, và phổ F(w) không thể được khôi phục từ F (w ) ; tức là f(t) không thể được khôi phục từ các mẩu f (t ) trong hình 5.2c. Trường hợp thứ hai, khi dùng tốc độ lấy mẫu Nyquist 10Hz (hình 5.2e). Phổ F (w ) gồm các thành phần phổ T F (w ) không trùng lắp lặp lại 1 từng 10 Hz. Do đó, phổ F(w) có thể được khôi phục từ F (w ) dùng mạch lọc thông thấp lý tưởng có băng thông 5 Hz (hình 5.2f). Sau cùng, trường hợp cuối là lấy mẩu lố (tốc độ lấy mẩu 20Hz); phổ F(w) gồm các gồm các thành phần phổ T F (w ) không trùng lắp (lặp 1
lại từng 20 Hz) với các băng tần trống giữa các chu kỳ liên tiếp. Do đó, phổ F(w) có thể được khôi phục từ F (w ) dùng mạch lọc thông thấp lý tưởng hay mạch lọc thực tế cũng được (tín hiệu vẽ chấm trong hình 5.2h). ¢ r Bài tập E5.1 Tìm tốc độ Nyquist là khoảng Nyquist cho tín hiệu: (a) sinc(100pt) và sinc(100pt) + sinc(50pt) . Đáp số: khoảng Nyquist là 0,01 giây và tốc độ Nyquist là 100Hz cho cả hai tín hiệu. s 5.1-1 Khôi phục tín hiệu: Công thức nội suy Quá trình khôi phục tín hiệu liên tục f(t) từ các mẩu còn được gọi là phép nội suy. Trong phần 5.1, ta thấy tín hiệu f(t) có băng thông giới hạn B Hz có thể được khôi phục (nội suy) chính xác từ các mẩu của mình. Quá trình khôi phục bằng cách cho tín hiệu đã lấy mẩu qua mạch lọc thông thấp lý tưởng có băng thông B Hz. Từ phương trình (5.3), tín hiệu đã lấy mẩu chứa các thành phần (1/T)f(t) và để khôi phục được f(t) (hay F(w)), thì tín hiệu đã lấy mẩu phải qua mạch lọc lý tưởng có băng thông B Hz và độ lợi T. Do đó, hàm truyền của bộ lọc khôi phục (hay nội suy) là: æwö H (w ) = Trect ç ÷ (5.7) è 4pB ø Quá trình nội suy được biểu diễn trong miền tần số như là tác động lọc. Ta xem xét tiếp quá trình này từ quan điểm khác, tức là trong miền thời gian. Để bắt đầu, ta hảy xét mạch lọc cực kỳ đơn giản có đáp ứng xung là rect ( T ) , vẽ t trong hình 5.3a. Đây là xung cổng (gate pulse) có tâm tại gốc, có chiều cao đơn vị, và độ rộng T (thời gian lấy mẩu). Ta tìm ngõ ra của bộ lọc khi ngõ vào là tín hiệu đã lấy mẩu f (t ) , Từng mẩu của f (t ) , đã được biến thành xung, tạo ngõ ra là xung cổng với chiều cao bằng với cường độ của mẩu. Thí dụ, mẩu thứ k là xung có cường độ f(kT) nằm tại t - kT và có thể viết thành f (kT )d (t - kT ) . Khi xung qua bộ lọc, tạo ra ngõ ra là xung
cổng có cao độ f(kT), tâm tại t = kT (phần nét gián đoạn trong hình 5.3b). Từng mẩu trong f (t ) sẽ tạo ra xung cổng tương ứng tại ngõ ra bộ lọc dạng xấp xỉ hình bậc thang của f(t), vẽ trong phần chấm trong hình 5.3b. Bộ lọc cho ta dạng thô của phép nội suy. Hàm truyền H(w) của bộ lọc là biến đổi Fourier của đáp ứng xung rect ( T ) . Giả t sử, ta dùng tốc độ Nyquist; tức là T = 1/2B. ætö h(t ) = rect ç ÷ = rect (2 Bt ) èT ø Và æ wT ö 1 æ w ö H (w ) = T sin cç ÷= ç ÷ (5.8) è 2 ø 2B è 4B ø Đáp ứng biên độ H (w ) của bộ lọc này, như vẽ trong hình 5.3c, giải thích lý do về tính thô của phép nội suy. Bộ lọc này, còn được gọi là mạch lọc giữ bậc zêrô (zero hold filter), dạng xấu nhất của mạch lọc thông thấp lý tưởng (phần tô bóng trong hình 5.3c) cần có cho phép nội suy chính xác. Ta có thể cải thiện dạng mạch lọc bậc zêrô bằng cách dùng mạch lọc giử bậc một (first order hold filter), cho ta phép nội suy tuyến tính thay vì dạng nội suy theo bậc thang. Bộ nội suy tuyến tính này, với đáp ứng xung là xung tam giác D( 2tT ) , cho phép xấp xỉ với đỉnh các mẩu được nối nhau bằng đường thẳng (xem bài tập 5.1-5). Hàm truyền bộ lọc nội suy lý tưởng lấy từ phương trình (5.7) và vẽ trong hình 5.4a. Đáp ứng xung của bộ lọa này là biến đổi Fourier nghịch của H(w) là h(t ) = 2 BT sin c (2pBT ) (5.9a) Giả sử lấy mẩu với tốc độ lấy mẩu Nyquist; tức là 2BT =1, thì h(t ) = sin c(2pBT ) (5.9b) f (t ) = å f (kT )h(t - kT ) k
f (t ) = å f (kT ) sin c[2pB(t - kT )] (5.10a) k f (t ) = å f (kT ) sin c(2pBt - kp ) (5.10b) k Phương trình (5.10) là công thức nội suy, tìm các giá trị của f(t) giữa các mẩu là phép cộng dồn (weighted sum) mọi giá trị mẩu. Thí dụ 5.2 ¢ Tìm tín hiệu f(t) với băng thông giới hạn B Hz, và có các mẩu là f(0) = 1 và f(±T) = f(±2T) = f(±3T) = . . . = 0 trong đó khoảng lấy mẩu là khoảng Nyquist của f(t); tức là T = 1/2B. Dùng công thức nội suy (5.10b) để khôi phục f(t) từ các mẩu. Do tất cả các mẩu Nyquist đều là zêrô, trừ một mẩu (tương ứng với k = 0) trong tổng bên vế phải của phương trình (5.10b). Do đó f(t) = sinc (2pBt) Tín hiệu này được vẽ trong 5.4b. Quan sát thấy là chỉ tín hiệu có băng thông B Hz và giá trị mẩu là f(0) = 1 và f(nT) = 0 (n ¹ 0). Các tín hiệu khác không thỏa được các điều kiện ¢ này. 5.1-2 Khó khăn thực tế khi khôi phục tín hiệu Nếu tín hiệu được lấy mẩu với tốc độ Nyquist Fs = 2B Hz, phổ F (w ) bao gồm F(w) được lặp lại nhiều lần mà không có khoảng hở giữa các chu kỳ kế tiếp, như vẽ trong hình 5.5a. Như đã thấy trong phần 4.5, thì mạch lọc này là dạng không thực hiện được; mà chỉ có thể xấp xỉ gần đúng dùng vô số các khâu trễ trong đáp ứng. Nói cách khác, ta có thể khôi phục lại tín hiệu f(t) từ các mẩu dùng dùng vô số các khâu trễ. Một giải pháp thực tế là lấy mẩu tín hiệu với tốc độ cao hơn tốc độ Nyquist (Fs > 2B hay w s > 4pB) . Kết quả là F (w ) , bao gồm F(w) được lặp lại nhiều lần với một số hữu hạn các khoảng hở giữa các chu kỳ kế tiếp, như vẽ trong hình 5.5b. Có thể khôi phục F(w) từ F (w ) dùng mạch lọc thông thấp có đặc tính ngắt chậm (từ từ), vẽ đường chấm trong hình 5.5b. Nhưng ngay cả trong trường hợp này, độ lợi mạch lọc phải là zêrô sau chu kỳ thức nhất
của F(w) (xem hình 5.5b). Theo tiêu chuẩn Paley-Wiener, thì không thể thực hiện được mạch lọc này. Ưu điểm duy nhất trong trường hợp này là ta có thể xấp xỉ mạch lọc cần có dùng số khâu trễ ít nhất. Điều này cho thấy là trong thực tế không thể khôi phục chính xác được tín hiệu f(t) có băng thông giới hạn từ các mẩu của tín hiệu, ngay cả khi tốc độ lấy mẩu có cao hơn tốc độ Nyquist đi nữa. Tuy nhiên, khi tốc độ lấy mẩu càng tăng thì tín hiệu khôi phục càng gần với tín hiệu mong muốn. Sự bội bạc của trùm phổ (The treachery of aliasing) Trong thực tế, khi khôi phục tín hiệu từ các mẩu còn một khó khăn cơ bản khác. Định lý lấy mẩu vừa chứng minh đã dựa trên giả thiết là tín hiệu có băng thông giới hạn. Các tín hiệu thực tế đều có thời gian giới hạn; tức là có độ rộng xung hữu hạn. Ta có thể chứng minh (xem bài tập 5.1-10) là tín hiệu không thể đồng thời vừa có thời gian giới hạn và băng thông giới hạn. Nếu tín hiệu có thời gian giới hạn, thì không thể có băng thông giới hạn và ngược lại (nhưng tín hiệu có thể có thời gian không giới hạn và băng thông không giới hạn). Rõ ràng, mọi tín hiệu thực tế, cần có thời gian giới hạn, thì băng thông không hạn chế; chúng có băng thông vô hạn, và phổ F (w ) gồm các chu kỳ chồng lắp của F(w) lặp lại mỗi Fs Hz (tần số lấy mẩu) như vẽ trong hình 5.6. Do có băng thông vô hạn trong trường hợp này, nên chồng lắp phổ là đặc trưng không thay đổi, bất chấp tốc độ lấy mẩu. Do có phần đuôi chồng lắp, nên F (w ) không còn đủ thông tin về F(w), và không còn khả năng, ngay cả trong lý thuyết, để khôi phục f(t) từ các tín hiệu đã lấy mẩu f (t ) . Nếu tín hiệu lấy mẩu đi qua mạch lọc thông thấp lý tưởng, thì ngõ ra không phải là F(w) mà là một phiên bản F(w) bị méo dạng do hai nguyên nhân riêng biệt sau: 1. Phần đuôi của F(w) bị mất khi |F| > Fs/2 Hz; 2. Phần đuôi này tái xuất hiện với dạng gấp lại trong phổ. Chú ý là phổ ngang qua Fs/2 = 1/2T Hz. Tần số này được gọi là tần số gấp lại. Do đó, phổ tự gấp lại tại tần số gấp lại. Thí dụ, thành phần tần số Fs/2+ Fx + Fx lộ diện để “đóng vai” thành phần tần số thấp hơn (Fs/2) – Fx trong tín hiệu khôi phục được. Do đó, thành phần tần số cao hơn (Fs/2) lại tái xuất hiện như là thành phần tần số thấp hơn (Fs/2). Yếu tố đảo phần đuôi được gọi là gấp phổ (spectral folding) hay trùm phổ (aliasing), vẽ phần tô bóng trong hình 5.6. Trong qua trình trùm phổ này, ta không chỉ mất mọi thành phần tần số cao hơn (Fs/2) Hz, nhưng các thành phần này lại tái hiện (trùm phổ) như thành phần tần số thấp. Sự xuất hiện lại này còn phá hủy tính toàn vẹn của các thành phần tần số thấp, như vẽ trong hình 5.6.
Giải pháp: Bộ lọc chống trùm phổ Thực hiện như sau, ta biết nguy cơ tiềm tàng nằm ở các thành phần tần số lớn hơn (Fs/2) = (1/2T) Hz. Ta cần loại bớt (triệt) các thành phần này khỏi f(t) trước khi lấy mẩu f(t). Theo phương pháp này, ta chỉ mất các thành phần lớn hơn tần số gấp (Fs/2)Hz; do đó các thành phần này không thể tái hiện để làm hỏng các thành phần có tần số thấp hơn tần số gấp. Các thành phần tần số cao được lọc nhờ mạch lọc thông thấp lý tưởng có băng thông là (Fs/2)Hz. Bộ lọc này được gọi là bộ lọc trùm phổ. Chú ý là phải thực hiện việc chống trùm phổ trước khi lấy mẩu tín hiệu. Bộ lọc chống trùm phổ là mạch lọc lý tưởng nên không thực hiện được. Trong thực tế, ta dùng mạch lọc có tần số cắt dạng dốc đứng, nhằm làm suy giảm sắc nét các phổ còn sót lại của tần số gấp (Fs/2). Lấy mẩu thực tế Khi chứng minh định lý lấy mẩu, ta giả sử các mẩu lý tưởng có được bằng cách nhân tín hiệu f(t) với chuỗi xung có độ rộng hữu hạn, vẽ trong hình 5.7b. Tín hiệu lấy mẩu được vẽ trong hình 5.7c. Điều kinh ngạc là ta có thể khôi phục hay tái tạo tín hiệu f(t) từ tín hiệu lấy mẩu f (t ) trong hình 5.7c. Đáng ngạc nghiên hơn là tốc độ lấy mầu lại không thấp hơn tốc độ Nyquist. Tín hiệu f(t) có thể khôi phục bằng các cho f (t ) qua lọc thông thấp dù nó đã được lấy mẩu dùng chuỗi xung. Kết quả này có vẽ đáng tin cậy khi ta xét thực tế là việc khôi phục f(t) đòi hỏi kiến thức về giá trị của mẩu Nyquist. Thông tin này có sẳn hay nằm trong tín hiệu lấy mẩu f (t ) trong hình 5.7c do cường độ của mẩu thứ k là f(kT). Để chứng minh điều này một cách giải tích, ta thấy là chuỗi xung lấy mẩu pT(t) vẽ trong hình 5.7b, là tín hiệu tuần hoàn, có thể biểu diễn thành chuỗi Fourier 2p ¥ pT (t ) = C0 + å C n cos(nw s t + q n ) ws = T n =1
Và é ù ¥ f (t ) = f (t ) pT (t ) = f (t ) êC0 + å C n cos(nw s t + q n )ú ë û n=1 ¥ = C0 f (t ) + å C n f (t ) cos(nw s t + q n ) (5.11) n =1 Tín hiệu đã lấy mẩu f (t ) gồm có C0 f (t ) , C1 f (t ) cos(w s t + q1 ), C2 f (t ) cos(2w s t + q 2 ), …. Chú ý là thừa số thứ nhất C0 f (t ) là tín hiệu mong muốn và tất cả các thừa số khác đều là tín hiệu điều chế với phổ có tâm tại ±ws, ±2ws, ±3ws, . . . , như vẽ trong hình 5.7e. Rõ ràng tín hiệu f(t) có thể được khôi phục bằng cách cho f (t ) qua mạch lọc thông thấp, để có ws > 4pB ( hay Fs > 2B). ¢ Thí dụ 5.3 Để minh họa việc lấy mẩu thực tế, xét tín hiệu f (t ) = sin c 2 (5pt ) được lấy mẩu dùng chuỗi xung vuông pT (t ) vẽ ở hình 5.8c. Chu kỳ của pT (t ) là 0,1 giây, để có tần số cơ bản là 10Hz. Do đó w s = 20p . Chuỗi Fourier của pT (t ) có thể biểu diễn thành ¥ pT (t ) = C0 + å Cn cos nw s t n =1 sin ( n4 ) ; tức là Dùng phương trình (3.66) ta có C0 = và C n = p 1 2 np 4 2 2 2 1 1 ,K , C4 = 0, C5 = - C1 = , C3 = C0 = C2 = , , p 3p 5p p 4
Do đó 1 2 1 2 f (t ) cos 60pt + L f (t ) cos 20pt + f (t ) cos 40pt + f (t ) = f (t ) pT (t ) = f (t ) + p p 3p 4 Và 1 1 1 F (w ) = F (w ) + [ F (w - 20p ) + F (w + 20p )] + [ F (w - 4 0p ) + F (w + 40p )] 2p p2 4 1 [ F (w - 6 0p ) + F (w + 60p )] + L + 3p 2 Trong trường hợp này F (w ) = 0,2D( 20p ) . Phổ F (w ) vẽ trong hình 5.8e. Quan sát thấy w phổ gồm F (w ) được lặp lại theo chu kỳ 20p rad/s (10Hz). Do đó, không có trùng lắp giữa các chu kỳ, và F (w ) có thể được khôi phục dùng bộ lọc thông thấp lý tưởng với băng thông 5 Hz. Bộ lọc thông thấp lý tưởng có độ lợi đơn vị (và băng thông 5 Hz) sẽ cho phép thừa số thứ nhất bên vế phải của phương trình trên đi qua đầy đủ và triệt mọi thừa số khác. Do đó, ngõ ra y(t) là 1 y (t ) = f (t ) ¢ 4 5.3-1 Một số ứng dụng của định lý lấy mẩu Định lý lấy mẩu rất quan trọng trong phân tích, xử lý và truyền dẫn tín hiệu do cho phép ta thay thế tín hiệu liên tục theo thời gian bằng chuỗi rời rạc các số. Do đó, việc xử lý tín hiệu liên tục tương đương với việc xử lý chuỗi rời rạc các số. Phép xử lý này đưa ta đến lĩnh vực lọc số. Trong lĩnh vực thông tin, việc truyền các bản tin liên tục theo thời gian rút lại thành việc truyền các chuỗi số dùng chuỗi xung. Tín hiệu liên tục theo thời
gian f(t) được lấy mẩu và các giá trị mẩu được dùng thay đổi một số tham số của chuỗi xung tuần hoàn. Ta có thể thay đổi biên độ (hình 5.9b), độ rộng (hình 5.9c), vị trí (hình 5.9d) của xung tỉ lệ với các giá trị mẩu của tín hiệu f(t). Từ đó, ta có phương pháp điều chế biên độ xung (PAM: pulse-amplitude modulation), điều chế độ rộng xung (PWM: pulse-width modulation), hay điều chế vị trí xung (PPM: pulse-position modulation). Dạng điều chế xung quan trọng nhất hiện nay là điều chế xung mã (PCM: pulse-code modulation). Trong mọi trường hợp, thay vì truyền tín hiệu f(t), ta truyền tín hiệu điều chế xung tương ứng. Tại máy thu, ta độc thông tin từ tín hiệu hiệu điều chế xung và khôi phục tín hiệu analog f(t). Một ưu điểm của điều chế xung là cho phép truyền đồng thời nhiều tín hiệu trên cơ sở chia sẻ thời gian (ghép kênh bằng cách phân chia theo thời gian TDM: time-division multiplexing). Do tín hiệu điều chế xung chỉ chiếm một phần thời gian của kênh truyền, nên ta có thể truyền nhiều tín hiệu điều chế xung trên cùng một kênh bằng cách chuyển vị (inter-weaving) chúng. Hình 5.10 vẽ phương pháp ghép kênh hai tín hiệu PAM dùng phương pháp TDM. Dùng phương pháp này, ta có thể ghép kênh nhiều tín hiệu trên cùng một kênh truyền bằng cách giảm độ rộng xung. Một phương pháp khác để truyền đồng thời nhiều tín hiệu analog (baseband) là phương pháp ghép kênh bằng cách phân chia theo tần số (FDM: frequency - division multiplexing), thảo luận trong phần 4.8-4. Trong phương pháp FDM, nhiều tín hiệu được truyền đồng thời dùng phương pháp chia sẻ khổ sóng của kênh truyền. Phổ của từng bản tin được dời đến dải tần đặc thù, chưa bị các tín hiệu khác chiếm. Thông tin của nhiều tín hiệu được đặt trong các dải tần số không trùng lắp của kênh (hình 4.45). Có thể xem TDM và FDM là dạng đối ngẫu của nhau. Điều chế xung mã (PCM: Pulse Code Modulation) PCM rất hữu ích và được dùng rộng rãi nhất trong các phương thức điều chế xung vừa nêu. Về cơ bản PCM là phương pháp chuyển đổi tín hiệu analog thành tín hiệu số (chuyển đổi A/D). Một tín hiệu analog có đặc trưng là biên độ có thể có giá trị bất kỳ trong một tầm liên tục. Do đó, tín hiệu analog có thể có vô hạn các giá trị. Ngược lại, biên độ tín hiệu số chỉ có thể có số hữu hạn các giá trị. Một tín hiệu analog có thể được chuyển đổi thành tín hiệu số bằng cách lấy mẩu và lượng tử hóa (làm tròn). Nếu chỉ lấy mẩu tín hiệu analog không, thì cũng chưa có tín hiệu số do tín hiệu đã lấy mẩu vẩn còn có nhiều giá trị bất kỳ trong một tầm liên tục. Để số hóa tín hiệu này thì cần có quá trình lượng tử hóa, để biên độ này nằm trong khoảng lượng tử như vẽ trong hình 5.11a.
Các biên độ của tín hiệu analog f(t) nằm trong tầm (– V, V). Tầm này là được chia thành L khoảng con, có biên độ là rv = 2V/L. Tiếp đến, mỗi mẩu biên độ được xấp xỉ dùng giá trị trung điểm của khoảng con của mẩu đang giảm (xem hình 5.11a với L = 16). Rõ ràng là mỗi mẩu được xấp xỉ với một trong L giá trị. Như thế, tín hiệu được số hóa với các mẩu lượng tử là một trong L giá trị. Các tín hiệu này là tín hiệu số L – phân (L – ary ). Theo quan điểm thực tế, tín hiệu nhị phân (là tín hiệu chỉ có hai giá trị) là tín hiệu rất hấp dẫn do tính đơn giản, kinh tế và dễ sử dụng trong kỹ thuật. Ta có thể chuyển đổi tín hiệu L – phân thành tín hiệu nhị phân dùng phương pháp mã hóa xung. Hình 5.11b vẽ mã trong trường hợp L = 16. Mã này, dùng biểu diễn nhị phân của hệ octal từ 0 đến 15, được gọi là mã nhị phân tự nhiên (NBC: natural binary code). Còn có nhiều phương pháp khác để mã hóa nhị phân. Mỗi mức trong 16 mức được truyền đi với một mã nhị phân 4 bit. Để truyền dữ liệu nhị phân này, ta cần định nghĩa các dạng xung riêng biệt để chỉ thị các trạng thái nhị phân. Một khả năng là dùng dạng xung âm cho bit 0 và xung dương cho bit 1 để các mẩu có thể được truyền đi thành nhóm 4 bit nhị phân, như vẽ trong hình 5.11b. Tín hiệu có được là tín hiệu PCM. Tín hiệu analog f(t) được chuyển đỗi thành tín hiệu số (nhị phân). Tín hiệu nhị phân được gọi là bit, hiện đang là dạng viết tắt trong chuẩn công nghiệp. Băng thông của tín hiệu âm thanh vào khoảng 15kHz, nhưng nhiều thử nghiệm cho thấy tín hiệu phát âm rõ (nghe được rõ) không bị ảnh hưởng nếu các thành phần tần số cao hơn 3400Hz bị loại bỏ. Do mục tiêu của thông tin điện thoại là nghe được rõ chứ không trung thực (hi-fi), nên các thành phần có tần số cao hơn 3400Hz được lọc bỏ dùng bộ lọc thông thấp. Tín hiệu sau đó được lấy mẩu với tốc độ 800 mẩu trong một giây (8 kHz). Tốc độ này cố tình cao hơn tốc độ Nyquist là 6,8 kHz để tránh phải có các mạch lọc không thực hiện được khi khôi phục tín hiệu. Sau cùng mỗi mẩu được lượng tử thành 256 mức (L = 256), cần có nhóm 8 xung nhị phân để mã hóa các mẩu (28=256). Do đó, tín hiệu điện thoai số hóa gồm 8 x 8000 = 64000 hay 64 kbits/s dữ liệu, cần có 64.000 xung nhị phân trong một giây khi truyền. Đĩa CD là ứng dụng của PCM. Đây là tình huống dùng cho tín hiệu âm thanh hi-fi có băng thông 15 kHz. Dù tốc độ lấy mẩu Nyquist chỉ là 30 kHz, hiện người ta dùng tốc độ lây mẩu 44,1 kHz như lý do đã trình bày trước đây. Tín hiệu được lượng tử với số lượng lớn các mẩu (L = 65.536) nhằm giảm sai số do lượng tử. Các mẩu được mã hóa nhị phân rồi ghi vào CD. Ưu điểm của tín hiệu số Một số ưu điển của tín hiệu số so với tín hiệu analog được liệt kê dưới đây: 1. Truyền tín hiệu số ưu việt hơn so với tín hiệu analog do tín hiệu số có tính chống nhiễu kênh truyền và méo tốt hơn nhiều (khi nhiễu và méo nằm trong giới hạn). Bản tin tín hiệu số (nhị phân) trong hình 5.12a bị méo dao kênh truyền, như vẽ trong hình 5.12b. Nếu méo nằm trong giới hạn, ta có thể khôi phục dữ liệu mà không bị lỗi do chỉ cần có quyết định nhị phân khi nhận được xung dương hay xung âm. Hình 5.12c vẽ cùng dữ liệu với méo do kênh truyền và nhiễu. Trường hợp này dữ liệu có thể được khôi phục chính xác bao lâu mà méo và nhiễu kênh truyền cón nằm trong giới hạn. Điều này khác với trường hợp tín hiệu analog, khi với nhiễu hay méo dạng bất kỳ, sẽ làm sái dạng tín hiệu thu được.
2. Tuy nhiên, ưu điểm lớn nhất của truyền dẫn tín hiệu số so với truyền dẫn analog là khả năng tái tạo lại tín hiệu. Khi truyền dẫn analog, một tín hiệu bản tin khi truyền qua kênh (đường truyền), tín hiệu ngày càng yếu đi, trong khi nhiễu kênh truyền và méo dạng tín hiệu, ngày càng tích lủy và tăng mạnh dần. Tín hiệu sau cùng bị áp đảo do nhiễu và méo dạng càng bị tổn thất lớn. Khuếch đại chỉ giúp rất ít do đồng thời tăng cường tín hiệu và nhiễu theo cùng tỉ lệ. Do đó, cự ly truyền bản tin analog bị giới hạn bởi công suất truyền. Khi truyền đủ dài, méo kênh truyền và nhiễu tích lủy đủ để áp đảo tín hiệu kể cả tín hiệu số. Có thể khắc phục bằng cách thiết lập các trạm lặp (repeater) dọc theo đường truyền với cự ly đủ ngắn để có thể phát hiện xung tín hiệu trước khi nhiễu và méo dạng có thể tích lủy đủ lớn. Mỗi trạm xung phát hiện xung là, sạch xung và truyền tiếp đến trạm lặp kế tiếp, để thực hiện lại các bước vừa nêu. Nếu nhiễu và méo vẫn còn nằm trong giới hạn (điều này thực hiện được do các trạm được đặc tương đối gần nhau), xung sẽ được phát hiện chính xác. Theo phương cách này, bản tin tín hiệu số có thể được truyền trên cự ly dài với độ tin cậy cao. Ngược lại, bản tin analog không thể làm sạch tín hiệu một cách tuần hoàn, nên truyền dẫn có độ tin cậy thấp hơn. Trong PCM thì sai số có ý nghĩa nhất là sai số lượng tử. Sai số này có thể được giảm thiểu bằng cách gia tăng mức lượng tử, với giá phải trả là làm băng thông lớn hơn trong môi trường truyền dẫn (kênh truyền). 3. Thiết lập phần cứng số rât mềm dẽo và cho phép dùng vi xử lý, các máy tính mini, chuyển mạch số, và mạch tích hợp cở lớn LSI. 4. Tín hiệu số có thể được mã hóa để có mức sai số thấp nhất và độ trung thực cao, cũng như tinh bảo mật cao. Đồng thời ta cũng có thể dùng những thuật toán tinh vi để xử lý tín hiệu số. 5. Phương pháp ghép kênh (multiplex) nhiều tín hiệu số thường dễ dàng và hiệu quả. 6. Phương thức lưu trữ tín hiệu số thường tương đối dễ và ít tốn kém. Ngoài ra còn có khả năng tìm kiếm và chọn lọc thông tín từ xa. 7. Việc tái tạo bản tin số thường có độ tin cậy cực lớn mà không bị xấu đi. Bản tin analog như bản photocopy hay phim, suy giảm chất lượng sau mỗi bước của nhiều bước khôi phục, và quá trình di chuyển từ nơi nay đến nơi khác thường tốn kém.
8. Giá các phần cứng số giảm nửa sau hai hay ba năm, trong khi tính năng và dung lượng tăng đôi với cùng thời gian này. Hiện nay các bước phát triển của công nghệ số ngày càng nhanh và càng tinh vi. Trong những năm gần đây, ta đã thấy đĩa CD, là thiết bị số đã làm cáo chung kỹ thuật ghi âm analog; báo chí truyền hình ảnh được scan theo dạng số; và ngày nay đã làm nước Mỷ phải hướng về các tiêu chuẩn mới cho truyền hình độ nét cao đối chọi lại với các chuẩn truyền hình analog từ Nhật và Châu Âu (?!!). Ngược lại thì công nghệ analog như giấy in, viđêo, âm thanh, và phim vẫn chưa có bước giảm giá nhanh. Dù gì đi nữa thì hiện nay công nghệ số đã và đang thống trị các lĩnh vực như công nghệ thông tin và công nghệ lưu trữ. Như đã nói trên, tín hiệu số có thể được tạo ra từ nhiều nguồn khác nhau. Một số nguồn tín hiệu như máy tính đã là số. Một số nguồn là analog, nhưng đã được chuyển đổi thành dạng số theo nhiều kỹ thuật khác nhau thí dụ như PCM. r Bài tập E5.2 Mã ASCII (American Standard Code for Information Interchange) có 128 ký tự dạng nhị phân. Một máy tính tạo ra 100.000 ký tự/giây. Chứng tõ là (a) Cần có 7 bit đễ mã hóa từng ký tự s (b) Cần có 700.000 bit/giây để truyền ngõ ra của máy tính. 5.1-3 Đối ngẫu của phép lấy mẩu theo thời gian: Định lý lấy mẩu phổ Tương tự như các trường hợp khác, định lý lấy mẩu cũng có phần đối ngẫu. Trong phần 5.1, ta đã thảo luận về định lý lấy mẩu theo thời gian, và đã chứng tõ được là tín hiệu với băng thông giới hạn B Hz có thể được khôi phục từ tín hiệu được lấy mẩu với tốc độ Fs > 2B mẩu/giây. Chú ý là phổ tín hiệu tồn tại trong tầm tần số từ – B đến B Hz. Do đó, 2B là độ rộng phổ (không phải là băng thông, là B) của tín hiệu. Điều này có nghĩa là tín hiệu f(t) có thể được khôi phục từ các mẩu được lấy với tốc độ Fs lớn hơn độ rộng phổ (tính bằng Hz) của tín hiệu. Đối ngẫu của định lý lấy mẩu theo thời gian là định lý lấy mẩu theo tần số. Định lý này áp dụng cho các tín hiệu có thời gian giới hạn, là đối ngẫu của tín hiệu có băng thông giới hạn. Ta tiếp tục chứng tõ là phổ F(w) của tín hiệu có thời gian giới hạn là t giây có thể được khôi phục từ mẩu của F(w) được lấy với tốc độ R > t (độ rộng của tín hiệu) mẩu/Hz. Hình 5.13a vẽ tín hiệu f(t) có thời gian giới hạn t giây cùng biến đổi Fourier F(w) của nó. Mặc dù F(w) thường có dạng phức, nhưng ta chỉ chứng minh cho trường hợp F(w) là hàm thực là đủ.
t ¥ F (w ) = ò f (t )e - jwt dt = ò f (t )e - jwt dt (5.12) -¥ 0 Tạo fTo(t) là tín hiệu tuần hoàn tạo ra bằng cách lặp lại f(t) mỗi T0 giây (T0 > t), như trong hình 5.13b. Tín hiệu tuần hoàn này có thể được biểu diễn theo chuỗi Fourier mủ 2p ¥ fT0 (t ) = å Dn e jnw0t w0 = T0 n = -¥ Trong đó (giả sử t < T0) 1t 1 T0 Dn = ò f (t )e - jnw0t dt = ò f (t )e - jnw0t dt T0 0 T0 0 Phương trình (5.12) cho 1 Dn = F (nw0 ) T0 Kết quả này cho thấy là các hệ số của chuỗi Fourier của fTo(t) là (1/T0) lần các giá trị mẩu của phổ F(w) lấy tại các khoảng w0. Điều này tức là phổ của tín hiệu tuần hoàn fTo(t) là phổ F(w) được lấy mẩu, như vẽ trong hình 5.13b. Bao lâu mà t còn < T0, chu kỳ kế tiếp của f(t) xuất hiện trong fTo(t) không bị trùng lắp, nên có thể khôi phục f(t) từ fTo(t). Khôi phục này gián tiếp hàm ý là F(w) có thể được khôi phục từ các mẩu của mình. Các mẩu này phân cách với tần số cơ bản F0 =1/T0 Hz của tín hiệu tuần hoàn fTo(t). Điều kiện để khôi phục là T0 ³ t; tức là F0 £ (1/t) Hz Do đó, để có thể khôi phục phổ F(w) từ các mẩu của F(w), thì các mẩu cần được lấy với khoảng tần số không lớn hơn F0 = (1/t) Hz. Nếu R là tốc độ lấy mẩu (mẩu/Hz) thì R = (1/F0) ³ t mẩu/Hz (5.13) Nội suy phổ Phổ F(w) của tín hiệu f(t) giới hạn theo thời gian t giây có thể được khôi phục từ các mẩu của F(w), Trong trường hợp này, dùng tính đối ngẫu của công thức nội suy tín hiệu trong phương trình (5.10). ta có công thức nội suy phổ 2p æ wt ö F (w ) = å F (nw0 ) sin cç - np ÷ w0 = (5.14) t è2 ø n ¢ Thí dụ 5.4 Phổ F(w) của tín hiệu có độ rộng đơn vị f(t) được lấy mẩu tại các khoảng 1 Hz hay 2p rad/s (tốc độ Nyquist). Các mẩu là: F(0) = 1 và F(±2pn) = 0 (n = 1, 2, 3, . . . ) Dùng công thức nội suy (5.14) để xây dựng F(w) từ các mẩu của mình. Do tất cả trừ một mẩu Nyquist đều là zêrô, chỉ còn một thừa số (tương ứng với n = 0) trong tổng bên vế phải của phương trình (5.14). Do đó, với F(0) = 1 và t = 1, ta có æw ö F (w ) = sin cç ÷ (5.15) è2ø Với tín hiệu có độ rộng đơn vị thì chỉ có một phổ với giá trị mẩu là F(0) = 1 và F(±2pn) = 0 (n ¹ 0 ). Các phổ khác không thỏa các điều kiện này. ¢
Tìm biến đổi Fourier bằng phương pháp số - Biến đổi Fourier rời rạc DFT 5.2 Tính toán số biến đổi Fourier của f(t) cần giá trị các mẩu của f(t) do máy tính số chỉ có thể xử lý các dữ liệu dạng rời rạc (chuỗi các số hạng). Do đó, máy tính chỉ có thể tính F(w) tại các giá trị rời rạc của w [mẩu của F(w)]. Do đó cần có các mẩu của F(w) để lấy mẩu f(t). Công việc này được thực hiện dùng kết quả của hai định lý đã phát triển trong phần 5.1. Bắt đầu với tín hiệu hữu hạn theo thời gian f(t) (hình 5.14a) và các phổ F(w) (hình 5.14b). Do f(t) bị giới hạn theo thời gian, nên F(w) không có giới hạn về băng thông. Để thuận tiện, ta cho các phổ đều là hàm của biến tần số F (tính theo Hz) thay vì theo w. Theo định lý lấy mẩu, thì phổ F (w ) của tín hiệu đã lấy mẩu f (t ) gồm có F(w) được lặp lại theo Fs (Hz) với Fs = 1/T. Điều này được vẽ trong hình 5.14c và 5.14d. Bước tiếp theo, tín hiệu đã lấy mẩu trong hình 5.14c được lặp lại tuần hoàn mỗi T0 giây, như vẽ trong hình 5.14e. Định lý lấy mẩu phổ đòi hỏi phải lấy mẩu phổ với tốc độ T0 mẩu/Hz. Tốc độ này tức là các mẩu cách nhau F0 = 1/T0 Hz , như vẽ trong hình 5.14f. Phần thảo luận trên cho thấy là khi tín hiệu f(t) được lấy mẩu và lặp lại theo chu kỳ thì phổ tương ứng cũng được lấy mẫu và lặp lại theo chu kỳ. Ta cần tìm quan hệ giữa tốc độ lấy mẩu của f(t) với tốc độ lấy mẩu của F(w). Số lượng mẩu Hình 5.14e và 5.14f cho thấy điều thú vị là số mẩu N0 của tín hiệu trong một chu kỳ T0 tại hình 5.14e giống hệt số mẩu N 0/ của phổ trong một chu kỳ Fs tại hình 5.14f. Lý do là: N0 = (T0/T) và N 0/ = (Fs/ F0) (5.16a) Nhưng, do Fs= 1/T và F0= 1/T0 (5.16b) / N0 = (T0/T) = (Fs/ F0) = N 0 (5.16c) Trùm phổ và rò phổ khi tính toán số Hình 5.14f cho thấy sự hiện diện của trùm phổ trong các mẩu của phổ F(w). Sai số trùm phổ có thể giảm thiểu được bằng cách tăng tần số lấy mẩu Fs (giảm khoảng lấy mẩu T = 1/ Fs). Tuy nhiên, không thể loại trừ hoàn toàn trùm phổ của tín hiệu có thời gian giới hạn f(t), do có phổ F(w) không giới hạn băng thông. Nếu ta đã bắt đầu với tín hiệu có phổ bị giới hạn về băng thông F(w), thì chúng không bị trùm phổ trong phổ ở hình 5.14f. Điều không may là tín hiệu không bị giới hạn về thời gian và các phần lặp lại (trong hình 5.14e) sẽ làm tín hiệu bị chồng khớp (trùm phổ trong miền thời gian). Trong trường hợp này, ta sẽ đấu tranh với sai số khi lấy mẩu tín hiệu. Nói cách khác, khi tính toán biến đổi Fourier trực tiếp hay biến đổi nghịch bằng phương pháp số, ta có thể giảm sai số tủy ý, nhưng không bao giờ có thể loại trừ được sai số. Đây là điều thực tế khi tính toán số biến đồi Fourier thuận và nghịch, bất chấp phương pháp nào được dùng. Thí dụ, nếu xác định biến đổi Fourier bằng cách lấy trực tiếp tích phân bằng phương pháp số, dùng phương trình (4.8), thì sai số do khoảng lấy tích phân rt không bao giờ có thể là zêrô. Chú ý tương tự cho trường hợp tìm biến đổi nghịch với phương pháp số. Do đó, ta cần luôn nhớ về bản chất của sai số này trong kết quả của mình. Trong phần thảo luận (hình 5.14) của mình, ta đã giả sử f(t) là tín hiệu có giới hạn về thời gian. Nếu f(t) không bị giới hạn về
thời gian, ta cần phải giới hạn thời gian do tính toán số chỉ hoạt động được với dữ liệu hữu hạn mà thôi. Hơn nữa, việc cắt cụt dữ liệu tạo sai số do việc trãi phổ (nhòe phổ) và rò phổ, như trình bày trong phần 4.9. Rò phổ cũng tạo ra trùm phổ. Có thể giảm thiểu rò phổ bằng cách cắt cụt tín hiệu bằng cửa sổ hình nón. Nhưng lựa chọn này làm tăng yếu tố trãi phổ hay nhòe phổ. Trãi phổ có thể được giảm thiểu bằng cách tăng độ rộng cửa sổ (tức là cho nhiều dữ liệu hơn), điều này làm tăng T0, và làm giảm F0 (tăng độ phân giải tần số hay phổ). Ảnh hưởng hàng rào (picket fence effect) Phương pháp tính toán số chỉ cho các giá trị mẩu đồng đều của F(w) [hay f(t)]. Dùng phương pháp này giống như việc nhìn tín hiệu các phổ của nó qua một ‘hàng rào”. Các đỉnh chủ yếu có thể nằm giữa hai mẩu và vẫn có thể không xuất hiện, tình hình này tạo một hình ảnh sai về thực tế. Các kết quả sai lệch này có thể tránh được khi dùng số lượng mầu N0 đủ lớn, làm tăng độ phân giải. Ta còn có thể dùng công thức nội suy phổ (phương trình 5.14) để xác định giá trị của F(w) giữa hai mẫu. Tìm biến đổi Fourier rời rạc (DFT: Discrete Fourier Transform) Nếu f(kT) và F(rw0) lần lượt là các mẩu thứ k và thứ r, ta định nghĩa các biến mới fk và Fr như sau T f k = Tf (kT ) = 0 f (kT ) (5.17a) N0 Fr = F (rw0 ) Và (5.17b) Trong đó w0 = 2pF0 = (2p/T0) (5.17c) Ta sẽ chứng minh fk và Fr quan hệ với nhau theo các phương trình sau: N 0 -1 åfe - jrW0 k Fr = (5.18a) k k =0 2p 1 N 0 -1 jrW0k å Fr e W 0 = w0T = fk = (5.18b) N 0 k =0 N0 Các phương trình này định nghĩa biến đổi Fourier rời rạc thuận và nghịch (DFT), với Fr là biến đổi Fourier rời rạc trực tiếp (DFT) của fk, và fk là biến đổi Fourier rời rac nghịch (IDFT: Inverse Fourier Transform) của Fr. Ý niệm
fk Û Fr còn được dùng để cho thấy fk và Fr là cặp DFT. Xin nhớ rằng fk là (T0/N0) nhân với mẩu thứ k của Fr là mẩu thứ r của F(w). Biết giá trị mẩu của f(t), ta có thể tính được giá trị mẩu của F(w) – và ngược lại – dùng DFT. Tuy nhiên, chú ý là fk là hàm của k (k = 0, 1, 2, . . . , N0 – 1) chứ không phải của t và Fr là hàm của r (r = 0, 1, 2, . . . , N0 – 1) thay vì w. Hơn nữa, cả fk và Fr là chuỗi tuần hoàn với chu kỳ N0 (Các hình 4.14e và 5.14f). Các chuỗi này được gọi là chuỗi tuần hoàn N0. Việc chứng minh quan hệ DFT trong phương trình (5.18) lấy trực tiếp từ kết quả của định lý lấy mẩu. Tín hiệu được lấy mẩu f (t ) (hình 5.14c) có thể viết thành N 0 -1 å f (kT )d (t - kT ) f (t ) = (5.19) k =0 Do d(t – kT) ó e– jkwT , biến đổi Fourier của phương trình (5.19) cho N 0 -1 å f (kT )e - jkwT F (w ) = (5.20) k =0 ws Nhưng từ hình 5.1e (hay phương trình 5.4), rõ ràng là trong khoảng w £ , với F (w ) là 2 F (w ) biến đổi Fourier cỉa f (t ) là , giả sử trùm phổ có thể bỏ qua, thì T N 0 -1 ws F (w ) = TF (w ) = T å f (kT )e - jkwT w£ 2 k =0 Và N 0 -1 Fr = F (rw0 ) = T å f (kT )e - jkrw0T (5.21) k =0 Nếu ta đặt w0T = W 0 , thì từ phương trình (5.16a) và (5.16b) W 0 = w0T = 2pF0T = (2p/N0) (5.22) Đồng thời, từ phương trình (5.17a), Tf(kT) = fr Do đó, phương trình (5.21) trở thành 2p N 0 -1 Fr = å f r e - jrW0k W0 = (5.23) N0 k =0 Quan hệ của biến đổi nghịch (5.18b) có thể được tìm dùng phương pháp tương tự với vai trò của t và w được đảo ngược lại, nhưng ở đây là dùng cách chứng minh trực tiếp hơn. Để chứng minh quan hệ nghịch trong phương trình (5.18b), ta nhân hai vế của phương trình (5.23) với ejmW0r , rồi lấy tổng trong khoảng r theo é ù N 0 -1 N 0 -1 N 0 -1 å Fr e jmW0r = å ê å f r e - jrW0k ú e jmW0r r =0 ë k =0 û r =0 Bằng cách thay đổi thứ tự lấy tổng bên vế phải é N0 -1 ù N 0 -1 N 0 -1 å Fr e jmW0r = å f k ê å f r e j ( m-r )W0k ú ë k =0 û r =0 r =0 Phụ chương 5.1 cho thấy tổng bên trong của vế phải là zêrô với k ¹ n và tổng là N0 khi k = m. Do đó, tổng bên ngoài sẽ chỉ có một thừa số khác zerô khi k = m, và là
N 0 f k = N 0 f m . Do đó thường ta xác định Fr trong tầm (0, N0 – 1) hơn là trong tầm (- , N0 - 1 ). N0 2 2 Chọn T, T0 và N0. Khi tính DFT, đầu tiên cần chọn giá trị thích hợp cho T, T0 và N0. Nhằm mục đích này, đầu tiên phải tìm băng thông cơ bản B (Hz) của tín hiệu. Tần số lấy mẩu Fs ít nhất phải bằng 2B; tức là, (Fs/2) ³ B (5.25a) Hơn nữa, khoảng lấy mẩu T = 1/ Fs , (phương trình 5.16b), và T £ 1/2B (5.25b) Khi đã chọn xong B, ta có thể chọn T từ phương trình (5.25b). Đồng thời F0 =1/T0 (5.26) Với F0 là độ phân giải tần số [mức phân biệt giữa các mẩu của F(w)]. Do đó, nếu đã có F0 , ta có thể chọn T0 theo phương trình (5.26). Biết được T0 và T, ta xác định N0 từ: N0 = T0/T (5.27) Các điểm gián đoạn (discontinuity) Nếu f(t) có bước nhảy gián đoạn tại điểm lấy mẩu, thì nên lấy giá trị mẩu là trung bình các giá trị hai bên của điểm gián đoạn do biểu diễn Fourier tại điểm gián đoạn hội tụ tại trị trung bình. Đệm zêrô (Zero padding) Nhắc lại khi quan sát Fr thì giống như quan sát phổ F(w) qua “ hàng rào”. Nếu khoảng tần số lấy mẩu F0 không đủ nhỏ, ta có thể bị mất một số chi tiết quan trọng và có được hình ảnh bị nhòe. Để có số mẩu cao hơn, ta cần giảm F0 . Do F0 = 1/T0, là số mẩu cao nhất cần có để tăng giá trị của T0, chu kỳ lặp lại của f(t). Chọn lựa này làm tăng N0, số mẩu của f(t), bằng cách thêm một số mẩu đệm có giá trị zêrô. Việc thêm các mẩu đệm này được gọi là đệm zêrô. Do đó, các zêrô đệm này làm tăng số mẩu và có thể trợ giúp để có thêm ý tưởng tốt hơn về phổ F(w) từ các mẩu Fr. Các zêrô đệm không cải thiện tính chính xác hay độ phân giải. Có một điểm cần được hiểu rõ là đệm zêrô thường chỉ cho ta thêm mẩu mà không cải thiện được tính chính xác của các giá trị mẩu này. Đệm zêrô sẽ chỉ hữu ích hơn nếu thời khoảng lấy mẩu đủ nhỏ để có thể bỏ qua sai số trùm phổ. Đệm zêrô không bao giờ có thể cải thiện tính chính xác hay độ phân giải tần số theo nghĩa thực. Chỉ có thể tăng tính chính xác bằng cách giảm trùm phổ (aliasing), điều này đòi hỏi giảm khoảng lấy mẩu tín hiệu T (T < 1/2B, với B là băng thông hiệu quả của tín hiệu). ¢ Thí dụ 5.5 Tín hiệu f(t) có độ rộng 2 ms và băng thông cơ bản là 10 kHz. Mong muốn có độ phân giải tần số là 100 Hz khi dùng DFT (F0 = 100). Tìm N0. Độ rộng hiệu quả của tín hiệu T0 là T0 = (1/ F0) = 1/100 = 10ms Do độ rộng tín hiệu chỉ có 2ms, ta cần đệm zêrô trong khoảng 8 ms. Đồng thời, với B = 10.000, nên Fs = 2B = 20.000 và T = 1/ Fs = 50ms. Vậy: