BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ………………………. Đồ án Nghiên cứu mã hoá band con Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 1 MỤC LỤC LỜI MỞ ĐẦU 3 Chương 1 5 LÝ THUYẾT CHUNG VỀ XỬ LÝ TÍN HIỆU SỐ 5 1.1. TÍN HIỆU VÀ HỆ THỐNG RỜI RẠC THEO THỜI GIAN 5 1.2. BIỂU DIỄN SỰ BIẾN ĐỔI CỦA TÍN HIỆU VÀ HỆ THỐNG 6 1.2.1. Biến đổi sang miền Z 7 1.2.2. Biến đổi Fourier 8 1.3. BỘ LỌC SỐ 9 1.3.1 Hệ thống FIR 11 1.3.2. Hệ thống IIR 13 1.4. LẤY MẪU 17 1.5. DFT VÀ FFT 19 1.5.1. DFT 19 1.5.2. FFT 21 Chương 2 28 Mà HÓA BAND CON 28 2.1. CÁC HỆ THỐNG LỌC SỐ NHIỀU NHỊP 28 2.1.1. Bộ lọc phân chia 28 2.1.2. Bộ lọc nội suy 33 2.1.3. Bé läc biÕn ®æi nhÞp lÊy mÉu víi hÖ sè h÷u tØ 36 2.2. BANK LỌC SỐ QMF 43 2.2.1. Bank lọc số phân tích 44 2.2.2. Bank lọc số tổng hợp 45 2.2.3. Bank lọc hai kênh QMF 45 2.3. Mà HÓA BAND CON CỦA TÍN HIỆU TIẾNG NÓI 51 2.3.1. Cấu trúc dạng cây phân giải đều 52 2.3.2. Cấu trúc dạng cây đa phân giải (Multiresolution) 55 2.4. MỘT SỐ LOẠI Mà 57 Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 2 2.4.1. Lượng tử hóa (Quantizing) 57 2.4.2. Mã hóa đều theo phương pháp so sánh 59 2.4.3. Mã hóa theo phương pháp phản hồi phi tuyến 64 2.5. GIẢI Mà 66 Chương 3 67 MÔ PHỎNG Mà HÓA BAND CON 67 3.1. GIỚI THIỆU VỀ MATLAB 67 3.2. CÁC KHỐI TRONG SIMULINK…………………………………… 68 3.2.1. Bộ lọc số (Digital Filter) 68 3.2.2. Bộ nội suy và bộ phân chia 68 3.2.3. Bộ mã hóa và giải mã 69 3.2.4. Bộ khuếch đại (Gain) 70 3.3. MÔ PHỎNG Mà HOÁ BAND CON 71 KẾT LUẬN 73 CÁC THUẬT NGỮ VÀ BẢNG CHỮ VIẾT TẮT DÙNG TRONG ĐỒ ÁN 74 TÀI LIỆU THAM KHẢO 75 Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 3 LỜI MỞ ĐẦU Trước kia, thông tin được xử lý hoàn toàn bằng tín hiệu tương tự hay khi tín hiệu số với các linh kiện điện tử các mạch logic phức tạp và cồng kềnh, giá thành lại cao. Ngày nay, đi liền với sự phát triển của khoa học kỹ thuật và công nghệ là sự phát triển vượt bậc của máy tính đã làm gia tăng một cách mạnh mẽ các ứng dụng của XỬ LÝ TÍN HIỆU SỐ (Digital Signal Proccesing). Trong xử lý tín hiệu, nhờ các linh kiên điện tử đã được tích hợp sẵn cùng những chiếc máy tính hiện đại gọn nhẹ, dễ sử dụng thì tin tức được số hóa và xử lý bằng các thuật toán đã được lập trình với tốc độ ngày càng cao. Do đó, xử lý tín hiệu số đã được ứng dụng rộng rãi trong nhiều lĩnh vực khác nhau như: thông tin liên lạc, phát thanh truyền hình, trong đo lường và tự động điều khiển và các nghành công nghệ khác. Trong xử lý tín hiệu do dải tần số đưa vào xử lý rất rộng, các thành phần tần số không mong muốn sẽ gây nhiễu cho tín hiệu sau xử lý. Đặc biệt là trong các lĩnh vực như xử lý tín hiệu hình ảnh, tín hiệu âm thanh trong khi dải tần phải xử lý là rất rộng, các thành phần tần số cao sẽ gây ra nhiễu tín hiệu khi xử lý, vậy vấn đề đặt ra ở đây là làm thế nào để có thể nén được tín hiệu hay thu hẹp dải tần tín hiệu xử lý mà vẫn không làm mất thông tin. Hiện nay, có rất nhiều phương pháp đã và đang được nghiên cứu rộng rãi để để xử lý tín hiệu âm thanh. Tất cả đều với mục đích chung là làm thế nào để biểu diễn tín hiệu âm thanh với ít bít nhất để giảm bề rộng của dải tần xử lý và loại bỏ các thành phần không mong muốn ở dải tần cao trong khi vẫn giữ được âm thanh trung thực. Do tín hiệu âm thanh (tiếng nói) thì năng lượng của phổ tiếng nói tập trung ở miền tần số thấp, ở miền tần số cao thì năng lượng của phổ âm thanh rất nhỏ. Các phương pháp nén tín hiệu trước đây, tiếng nói được mã hóa trong toàn bộ dải tần của tín hiệu, như vậy gây ra sự dư thừa thông tin khi mã hóa trong miền tần số cao. Ý tưởng của đề tài Mà HÓA BAND CON là chia dải tần Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 4 của tín hiệu âm thanh thành nhiều dải con và mã hóa ở mỗi dải tần một số lượng bít khác nhau, ở dải tần cao thì mã hóa với số bít thấp hơn ở dải tần số thấp sẽ làm giảm một cách đáng kể không gian lưu trữ trong truyền phát, điều này làm cho việc mã hóa hay nén tín hiệu âm thanh tối ưu hơn, và nó cũng làm giảm bớt các thành phần tín hiệu không mong muốn. Nội dung của đề tài được chia ra làm ba phần: Chương 1. Lý thuyết chung về xử lý tín hiệu số. Chương 2. Mã hóa band con Chương 3. Mô phỏng hệ thống mã hóa band con bằng Matlab-simulink. Trong quá trình làm đồ án được sự giúp đỡ rất nhiệt tình của các thầy, các cô cùng các bạn trong lớp. Đặc biệt là thạc sĩ Nguyễn Văn Dương người đã trực tiếp hướng dẫn tôi hoàn thành đồ án này. Tôi xin chân thành cảm ơn thạc sĩ Nguyễn Văn Dương, các thầy cô trong tổ bộ môn điện tử viễn thông đồng các thầy cô trong trường ĐHDL Hải Phòng và các bạn trong lớp ĐT901 đã giúp tôi hoàn thành tốt nhiệm vụ đồ án mà nhà trường và tổ bộ môn giao cho. Hải phòng, ngày 10 tháng 7 năm 2009 Sinh viên thực hiện: Nguyễn Thị Tuyến Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 5 Chương 1 LÝ THUYẾT CHUNG VỀ XỬ LÝ TÍN HIỆU SỐ 1.1. TÍN HIỆU VÀ HỆ THỐNG RỜI RẠC THEO THỜI GIAN Trong hầu hết các lĩnh vực có liên quan đến xử lý tin tức hoặc thông tin đều bắt nguồn với việc biểu diễn tín hiệu như một dạng mẫu thay đổi liên tục. Sóng âm tạo ra tiếng nói của con người cũng tuân theo nguyên tắc này. Từ các mẫu tín hiệu, để thuận tiện thì người ta biểu diễn tín hiệu bằng các hàm toán học, như các hàm biến đổi theo thời gian. Ví dụ, chúng ta có thể biểu diễn x a (t) là dạng sóng liên tục thay đổi theo thời gian (tín hiệu analog ). Ngoài ra, cũng có thể biểu diễn tín hiệu bằng các hàm toán học của các biến rời rạc, chúng ta biểu diễn dãy tín hiệu này là x(n). Nếu tín hiệu được lấy mẫu từ tín hiệu tương tự với chu kỳ lấy mẫu T, khi đó ta biểu diễn dạng của tín hiện là x a (nT). Trong các hệ thống xử lý số tín hiệu, chúng ta thường dùng đến các dãy đặc biệt như: Dãy nhảy đơn vị: Dãy nhảy bậc đơn vị: Dãy hàm mũ thực: Dãy mũ ảo: )).sin(.) (cos( . . ) ( )(      njn n e jn e n e jn en r e    (1.1.4) Dãy sin: Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 6 )sin()( 0  nnS  (1.1.5) Dãy chữ nhật: Xử lý tín hiệu, tức chúng ta phải chuyển đổi tín hiệu về các dạng mẫu của tín hiệu mà ta muốn. Do đó chúng ta phải quan tâm đến ở đây là các hệ thống rời rạc, hoặc khi chúng ta biến đổi một dãy tín hiệu vào để được một dãy tín hiệu ra. Sự biến đổi tín hiệu này có thể được miêu tả như hình 1.1. Hình 1.1. Mô phỏng hệ thống Những hệ thống như trên có thể hoàn toàn xác định được bằng đáp ứng xung của nó đối với mẫu xung đơn vị đưa vào. Đối với những hệ thống này, đầu ra có thể được tính khi đưa vào dãy x(n) và đáp ứng xung đơn vị h(n) nhờ công thức tổng chập:     k nhnxknhkxny )(*)()().()( (1.1.7) Dấu * ở đây dùng cho tổng chập. Ta có công thức tổng chập tương đương là:     k nxnhknxkhny )(*)()().()( (1.1.8) 1.2. BIỂU DIỄN SỰ BIẾN ĐỔI CỦA TÍN HIỆU VÀ HỆ THỐNG Phân tích và thiết kế của các hệ thống tuyến tính sẽ rất đơn giản nếu chúng ta sử dụng trong miền Z và miền tần số cho cả hệ thống và tín hiệu, khi đó chúng ta cần thiết phải xét đến sự biểu diễn Fourier, miền Z của hệ thống và tín hiệu rời rạc theo thời gian. T[] x(n) y(n)=T[x(n)] Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 7 1.2.1. Biến đổi sang miền Z Sự biến đổi sang miền Z của một dãy được định nghĩa bằng hai phương trình sau:      n n ZnxZX ).()( (1.2.1) Trong miền thời gian: dZZZX j nx C n    1 ).( 2 1 )(  (1.2.2) Từ công thức (1.2.1) ta thấy dãy tín hiệu x(n) được biến đổi sang miền Z (biến đổi Z thuận). Sau khi tín hiệu được biến đổi sang miền Z thì tín hiệu là một dãy lũy thừa đối với biến 1 X , giá trị của dãy x(n) biểu diễn bộ các hệ số trong dãy lũy thừa. Một cách chung nhất, điều kiện đủ để biến đổi sang miền Z là dãy lũy thừa phải hội tụ tại một giá trị giới hạn.      n n Znx )( (1.2.3) Một bộ các giá trị cho các dãy hội tụ đựợc định nghĩa bằng một vùng trong mặt phẳng Z. Nói chung miền này có dạng: R 1 < Z < R 2 (1.2.4) Bảng 1.1. là các tính chất của phép biến đổi Z ngược. Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 8 Các tính chất Dãy miền n Biến đổi Z 1. Tính tuyến tính a.x 1 (n) + b.x 2 (n) a.X 1 (Z) +b.X 2 (Z) 2.Tính dịch chuyển theo thời gian x(n+n 0 ) )(. 0 ZXZ n 3. Thay đổi thang tỉ lệ a n .x(n) X(a -1 .Z) 4. Vi phân của X(Z) theo Z n.x(n) -Z )( )( Zd ZdX 5. Đảo trục thời gian X(-n) )( 1 ZX 6. Tích chập của hai dãy x(n)*h(n) X(Z).H(Z) 7. Tích của hai dãy x(n).w(n)   C dVVVZWVX j .)./().( 2 1 1  Bảng 1.1. Các tính chất của phép biến đổi Z ngược Biến đổi Z ngược được đưa ra bởi tích phân đường trong phương trình (1.2.2), trong đó C là đường cong kín bao quanh gốc tọa độ trong mặt phẳng Z, nằm trong miền hội tụ của X(Z). Trong những trường hợp đặc biệt của phép biến đổi Z ngược, như xử dụng tính chất của phép biến đổi Z ngược. 1.2.2. Biến đổi Fourier Phép biến đổi Fourier của tín hiệu rời rạc theo thời gian được biểu diễn bằng công thức:      n njj enxeX . ).()(  (1.2.5)         deeXnx njj . ).( 2 1 )( (1.2.6) Những phương trình trên có thể nhận ra dễ dàng, là trường hợp đặc biệt của phương trình (1.2.1) và (1.2.2). Ngoài ra biểu diễn Fourier có thể đạt đựợc bằng cách giới hạn phép biến đổi Z vào vòng tròn đơn vị của mặt phẳng Z, thay Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P M M à à H H Ó Ó A A B B A A N N D D C C O O N N Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 9  j eZ  biểu diễn như hình (1.2), biến số  có thể biểu diễn bằng góc trong mặt phẳng Z. Điều kiện đủ để tồn tại biến đổi Fourier có thể tính bằng cách gán 1Z trong phương trình (1.2.3), ta có:     n nx )( (1.2.7) Hình 1.2. Vòng tròn đơn vị trong mặt phẳng Z Một đặc điểm quan trọng của biến đổi Fourier )(  j eX là một hàm tuần hoàn của  , tuần hoàn với chu kỳ là  2 , điều này có thể dễ nhận thấy bằng cách thay thế  2 vào phương trình (1.2.5). Một cách khác, bởi vì )(  j eX được tính bằng X(Z) trên vòng tròn đơn vị, nên chúng ta có thể thấy rằng )(  j eX phải lặp lại mỗi lần khi  quay hết một vòng quanh vòng tròn đơn vị (tương ứng với một góc là  2 Radian). Bằng cách thay  j eZ  vào mỗi công thức trong bảng (1.1), chúng ta có thể có được các công thức cho biến đổi Fourier. Tất nhiên kết quả này chỉ đúng với biến đổi Fourier khi phép biến đổi đã tồn tại. 1.3. BỘ LỌC SỐ Bộ lọc số là hệ thống tuyến tính bất biến theo thời gian. Thông số vào ra của hệ thống quan hệ với nhau bằng tổng chập trong phuơng trình (1.1.7) và (1.1.8), quan hệ trong miền Z được đưa ra trong bảng (1.1). Re[Z] Im[Z]  [...]... Do tính chất ưu việt của bộ lọc số nhiều nhịp này mà nó đã được nghiên cứu và ứng dụng nhiều trong kỹ thuật viễn thông, đặc biệt là trong xử lý tín hiệu số: xử lý tiếng nói, xử lý hình ảnh, các hệ thống antenna, kỹ thuật audio số Đặt biệt hơn là ứng dụng chính của nó là mã hóa band con (subband coding) trong xử lý tiếng nói, ta sẽ nghiên cứu ở phần sau Hệ thống xử lý số nhiều nhịp là hệ thống xử lý số... nhân tổng cộng là N log 2 N , bằng với phép nhân trong cách tính theo phương pháp phân chia 2 theo thời gian, số phép cộng cũng như vậy 27 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP Mà HÓA BAND CON Chương 2 Mà HÓA BAND CON 2.1 CÁC HỆ THỐNG LỌC SỐ NHIỀU NHỊP Kĩ thuật lọc số nhiều nhịp ngày càng được ứng dụng nhiều trong lĩnh vực xử lý số tín hiệu, như là nó có thể dùng để tăng tốc độ tính toán... cộng, nhân các giá trị của dãy với hằng số (các hằng số trên nhánh hàm ý phép nhân), và chứa các giá trị của dãy trước đưa vào Vì vậy biểu đồ khối đưa ra chỉ dẫn rõ rang về tính phức tạp của hệ thống 12 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP Mà HÓA BAND CON 1.3.2 Hệ thống IIR Nếu hàm hệ thống của phương trình (1.3.7) có các điểm cực cũng như điểm không, thì phương trình sai phân (1.3.5)... phần liên quan đến cực liên hợp phức, H(Z) có thể viết như sau: c0 k  c1k Z K H Z    k 1 1  a1k Z 1 1  a 2 k Z 2 (1.3.17) 15 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP Mà HÓA BAND CON Điều này gợi ý một dạng sơ đồ song song biểu diễn như hình (1.5b) cho N= 4 x(n) b10 + + b20 + Z-1 + + + Z-1 a11 a12 y(n) b11 Z-1 b12 + + + + a21 a22 b21 Z-1 b22 + + Hình 1.5a Cấu trúc dạng tầng c10... (Low pass filter) có C  M Mạch lọc thông thấp này có nhiệm vụ loại bỏ các thành phần tần số 30 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP    M Mà HÓA BAND CON , chỉ giữa lại thành phần tần số    M , như vậy sẽ tránh được hiện tượng chồng phổ Sơ đồ tổng quát của mạch lọc phân chia: FS FS h(n) x(n) FS M M y H (n) y H  M (n) Bộ lọc phân chia Hình 2.1.3 Mạch lọc phân chia Trong đó... có C   L Trong miền biến số n mạch lọc này làm nhiệm vụ nội suy ra các mẫu biên độ 0 Còn trong miền tần số nó loại bỏ các ảnh phụ cơ bản 34 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP Mà HÓA BAND CON Sơ đồ tổng quát của mạch lọc nội suy được biểu diễn trên hình 2.1.6 L x(n) y  L (n) h(n) y  LH (n) Bộ lọc thông thấp có  C   L h(n): đáp ứng xung của bộ lọc Hình 2.1.6 Bộ lọc nội suy... k 1 1  a1k Z   (1.3.16) K là phần nguyên của (N+1)/2 Hệ thống cấp hai này được biểu diễn như trong hình (1.5a) cho trường hợp N = M = 4 14 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP x(n) b0 Mà HÓA BAND CON y(n) + + Z-1 Z-1 b1 a1 + + Z-1 Z-1 b2 a2 + + Z-1 Z-1 b3 + + a3 Hình 1.4a Cấu trúc dạng trực tiếp x(n) w(n) b0 + y(n) + Z-1 a1 + b1 Z-1 a2 + b2 + + Z-1 a3 + b3 + Hình 1.4b Cấu trúc... 1 Hàn hệ thống H(Z) là một hàm hữu tỉ của Z Nó có thể được biểu diễn bằng dạng điểm cực và điểm không trong mặt phẳng Z Như vậy H(Z) có thể viết 10 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP Mà HÓA BAND CON dưới dạng: M H Z    A 1  c r Z r 1  1  d k Z   1 1 k 1 (1.3.7) Như chúng ta đã xét trong miền Z, hệ thống nhân quả sẽ có miền hội tụ dạng Z  R1 Nếu hệ thống cũng là... những ứng dụng lọc tuyến tính, dạng song song đưa ra những đặc tính cao hơn về phương diện làm tròn giảm tiếng ồn, các sai số hệ số, và tính ổn định 16 Sinh viên: NGUYỄN THỊ TUYẾN LỚP ĐT901 ĐỒ ÁN TỐT NGHIỆP Mà HÓA BAND CON 1.4 LẤY MẪU Để sử dụng phương pháp xử lý số tín hiệu đối với tín hiệu tương tự, chúng ta cần biểu diễn tín hiệu như một dãy các giá trị Để thực hiện các biến đổi, thông thường người... hợp 1/T2FN) . BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ………………………. Đồ án Nghiên cứu mã hoá band con Đ Đ Ồ Ồ Á Á N N T T Ố Ố T T N N G G H H I I Ệ Ệ P P . Bộ mã hóa và giải mã 69 3.2.4. Bộ khuếch đại (Gain) 70 3.3. MÔ PHỎNG Mà HOÁ BAND CON 71 KẾT LUẬN 73 CÁC THUẬT NGỮ VÀ BẢNG CHỮ VIẾT TẮT DÙNG TRONG ĐỒ