Đang tải... (xem toàn văn)
Biến đổi wavelet và ứng dụng trong các thiết bi ghi sóng tim công nghệ mới
BIếN Đổi WAVELET Và ứNG DụNG TRONG CáC THIếT Bi GHI SóNG TIM CÔNG NGHệ Mới (Xử lý tín hiệu trong đo lường) Tóm tắt: Bài viết này giới thiệu phép biến đổi Wavelet và ứng dụng của nó trong các thiết bị y tế hiện đại, như máy đo nhịp tim, điện tim thông minh, máy theo dõi loạn nhịp và máy đặt nhịp có thể cấy ghép. Chất lượng của chúng phụ thuộc chính vào hệ phân tích. Hệ này chịu ảnh hưởng rất nhiều vào độ chính xác và độ tin cậy của phép tách phức hợp QRS, sóng P và T. Biến đổi Wavelet có thể mô tả đặc trưng cục bộ của tín hiệu. Tính chất đa thang của biến đổi Wavelet có thể nhận biết được phức hợp QRS từ các sóng P, T có biên độ cao, từ các hiện tượng giả, trôi đường đẳng trị và nhiễu. Do hiệu quả của kỹ thuật này trong phép xử lý các tín hiệu không dừng và do đặc tính bền vững với nhiễu, phép biến đổi Wavelet đã nổi lên như là một công cụ rất mạnh trong kỹ thuật phân tích tín hiệu điện tim. I. ĐặT VấN Đề : Điện tim (Electrocardiograph - ECG) là một thiết bị đo thuộc nhóm thiết bị ghi. Vì vậy, cũng như các thiết bị đo lường khác, hệ xử lý tín hiệu là cốt lõi và quyết định độ tin cậy, độ chính xác của thiết bị. Với sự ra đời các loại vi mạch cỡ lớn (LSI) và cực lớn (VLSI) như m P, DSP hay vi mạch hệ những Neural ., cùng công nghệ phần mềm có giải bước đi nhảy vọt, chúng ta đã có khả nǎng được nhiều bài toán mà trước đây mới chỉ là những ý tưởng. ở các thiết bị ghi sóng tim công nghệ mới, vấn đề cơ bản là tách chính xác các điểm đặc trưng như phức hợp QRS, sóng P và T. Một đặc thù của tín hiệu y sinh là chúng thường chìm trong môi trường nhiễu và các hiện tượng giả (artifact). Thông thường các bộ tách phức hợp QRS được thiết kế ở môi trường nhiễu vừa phải và nó phải chứa ít nhất một bộ lọc thông dải để loại bỏ các sóng P, T và nhiễu. Sau đó tín hiệu được đẩy qua phép biến đổi không tuyến tính như đạo hàm, bình phương . Cuối cùng để làm rõ phức hợp QRS người ta phải dùng luật quyết định để xem phức hợp này có trong tín hiệu hay không. Kỹ thuật này bắt gặp một số trở ngại: - Dải thông tín hiệu của phức hợp QRS là khác nhau đối với các đối tượng khác nhau, thậm chí đối với các nhịp khác nhau của cùng một đối tượng. - Dải thông của phức hợp QRS và nhiễu phủ đè lên nhau. - Hiện tượng giả và trôi đường đẳng trị là thường trực và đôi khi khá trầm trọng trong các tín hiệu điện tim. Nhờ tính đa thang (multiscale) của phép biến đổi Wavelet mà ta có thể khắc phục được những trở ngại trên. Do hiệu quả của kỹ thuật này đối với phép xử lý tín hiệu không dừng và do tính bền vững với nhiễu mà phép biến đổi Wavelet đang nổi lên như là những công cụ rất mạnh trong kỹ thuật xử lý tín hiệu và giải các bài toán toán tử vi phân không gian. Không chỉ ở các nước phát triển mà ở một số nước đang phát triển gần chúng ta như Trung Quốc, Thái Lan, ấn Độ . phép biến đổi Wavelet đang được nghiên cứu rất nghiêm túc. Trong khuôn khổ bài viết này, chúng tôi xin phép được giới thiệu những vấn đề cơ bản của phép biến đổi Wavelet và ứng dụng điển hình trong kỹ thuật tách các điểm đặc trưng của sóng điện tim. II. Phép biến đổi wavelet : Như ta đã biết, phổ nǎng lượng là một đặc trưng cho tín hiệu trong miền tần số. Muốn đo hay tính toán nó ta thường dùng: - Phép biến đổi Fourier - Qua hàm tương quan - và các bộ lọc Mỗi phương pháp đều có ưu nhược điểm riêng. Để có hệ thống và dễ dàng tiếp cận phép biến đổi Wavelet, chúng ta hãy đi từ phép biến đổi Fourier 2.1. Phép biến đổi FOURIER : Phép biên đổi Fourier (Fourier Transform - FT) được thể hiện bởi: X (f) = e -2p jft .dt ( 1 ) Định nghĩa này bao trùm toàn miền không gian. Tần số f ở đây mang tính ý nghĩa toàn cục của tín hiệu, vì vậy nó được gọi là tần số toàn cục (global frequency). Bất kỳ một đột biến nào trong miền thời gian của tín hiệu x(t) được kéo theo và lan ra toàn trục tần số. Từ phép biến đổi Fourier của hàm x(t), chúng ta phải đánh giá các đột biến. Nhưng những thông tin này lại không được sắp xếp để định vị trong trường không gian, vì vậy chúng ta gặp khó khǎn trong việc tìm vị trí của chúng. Để giải quyết vấn đề này, người ta đưa khái niệm tần cố cục bộ (Local Frequency) và dẫn đến phép biến đổi Fourier thời gian ngắn (Short Time Fourier Transform - STFT). 2.2. Phép biến đổi Fourier thời gian ngắn : 2.2.1. Định nghĩa: Để vẫn dùng được phép biến đổi Fourier cho việc định vị thông tin, tác giả GABOR xây dựng một phương pháp mới. Ông dùng hàm cửa sổ không gian g(t) trong tích phân Fourier. Hàm cửa sổ này được truyền tịnh tiến dọc theo trục không gian nhằm phủ toàn bộ tín hiệu. ở vị trí t và tần số f, phép biến đổi Fourier cửa sổ (hay còn gọi là phép biến đổi Fourier thời gian ngắn) được xác định như sau: STFT (t ,f) = . G(t-t ). e -2p jft .dt (2) Trong đó x(t) là tín hiệu được quan sát qua cửa sổ g(t). STFT đã đưa tín hiệu x(t) vào hàm 2 chiều trong mặt phẳng thời gian và tần số (t, f). Như vậy, tần số ở đây là tần số cục bộ. 2.2.2. Phân tích tần số với Q thay đổi Phép phân tích dùng STFT phụ thuộc trước hết vào việc chọn cửa sổ g(t). Dạng tổng quát của nó thường là hàm chẵn, thực. Nó có thể được coi như là đáp ứng xung của bộ lọc thông thấp, và để chuẩn hoá, nǎng lượng của g(t) được coi là bằng 1 . || g || 2 = | g(t) | 2 . dt = 1 (3) Phổ nǎng lượng của FT g(t) tập trung ở tần số thấp (xem hình 1). Đường cong (a), (b) và (c) có các tần số khác nhau trong khoảng thời gian riêng. Đường cong (a'), (b') và (c') có độ rộng dải giống nhau, nhưng khác nhau về vị trí trên trục tần số vì vậy phương pháp này mang tên "phân tích tần số với Q thay đổi ". Như vậy để làm rõ biến thời gian của tần số, người ta đã đưa ra khái niệm cửa sổ. Mục đích đầu tiên của cửa sổ là giới hạn sự lan truyền của tín hiệu cần khảo sát sao cho các đặc tính phổ là dừng trong toàn vùng cửa sổ. Độ rộng dải D f lúc đó sẽ là: 2 | G(f) | 2 .df | G(f) | 2 .df Mẫu số của biểu thức chính là nǎng lượng của g(t). Về mặt tần số, ta chỉ phân biệt được hai tín hiệu khi chúng cách xa nhau một lượng Hình 1 : Phân tích tần số bằng STFT (a): Hàm cửa sổ g(t) (b):Đồ thị g(t).cos (w o t) (c): Đồ thị g(t).cos (2w o t) (a’), (b’), (c’) biến đổi Fourier của (a), (b), (c) D f. Vậy độ phân giải tần số của STFT được đặc trưng bằng ngay D f. (Kỳ sau đǎng tiếp) Tương tự ta cũng có chiều rộng thời gian D t : 2 .| g(t) | 2 .dt | g(t) | 2 .dt Và ta cũng chỉ phân biệt được 2 xung trong miền thời gian khi chúng cách nhau một khoảng lớn hơn D t. Như vậy độ phận giải thời gian và độ phân giải tần số không thể nhỏ tùy ý. Theo hai tác giả O. Rionl và M. Vetterli thì: D t. D f ? p (6) Nghĩa là ta có thể hy sinh độ phân giải thời gian vì độ phân giải tần số hoặc ngược lại, chứ không thể có được sự hoàn hảo cho cả hai. Đây là một nhược điểm cơ bản của STFT. Để khắc phục nhược điểm này, tác giả Mallat đã đưa ra phương pháp phân tích trên cơ sở phép giãn (dilation) và đưa đến một phép biến đổi mới đó là WAVELET TRANSFORM (WT). Phép biến đổi này sẽ sử dụng Q không đổi, hay độ rộng dải tương đối không đổi. 2.3. Biến đổi WAVELET (WAVELET TRANSFORM) : 2.3.1. Định nghĩa : Mallat thực hiện WT bằng phép khai triển tín hiệu thành các hàm họ, chúng là các hàm tịnh tiến và giãn của hàm đơn trị y (t). WT của hàm x(t) ở vị trí t và thang s được xác định bằng công thức sau: Wx(s, t ) = x(t) * y s (t) = .Y s (t-t ) .dt (7) [* là ký hiệu của tích chập] Trong đó y s (t) là phép giãn của y (t) với hệ số s và được tính: y s (t) = y ( ) (8) Khi tín hiệu chứa các cấu trúc quan trọng, các cấu trúc này thuộc về các thang khác nhau, thì thông tin tín hiệu phải được tổ chức lại thành các tập thành phần chi tiết của kích cỡ thay đổi. Bằng các phép giãn (dilations) và co (contractions) ta có thể nhận được những thành phần này từ hàm mẫu y (t). Hàm y (t) có thể được coi như bộ lọc thông dải (band pass filter - BDF) và tiếp theo là đặc tính Q không đổi của các bộ BPF khác. Bởi lẽ, trong WT, khái niệm thang được đưa vào như là phép luân phiên tần số dẫn đến việc biểu diễn thang thời gian. Điều này có nghĩa là tín hiệu đã được biểu diễn trong mặt phẳng thời gian - thang (t - s). 2.3.2. Phân tích tần số với Q không đổi Trái ngược với STFT có độ phân giải cố định trong miền tần số và không gian, còn độ phân giải của WT lại thay đổi theo thông số thang s. WT tách tín hiệu thành tập dải tần có kích thước không đổi trên thang logarith (xem hình 2). Các đường cong a’, b’, c’ có độ rộng dải giống nhau trên thang logarith. Hình 2 : Phân tích tần số bằng WT Điểm ngược nhau thứ 2 là STFT sử dụng cửa sổ phân tích đơn, còn WT sử dụng các cửa sổ bội theo những mục đích thích hợp: dùng cửa sổ hẹp ở tần số cao, cửa sổ rộng ở tần số thấp, nhằm đảm bảo là Q không đổi, hay độ rộng tương đối của bǎng không đổi. Để đảm bảo được điều kiện này, độ phân giải D f và D t (công thức 4 và 5) phải thay đổi trong mặt phẳng tần số - thời gian để nhận được phép phân tích đa phân giải (Multiresolution). Độ phân giải thời gian trở nên tùy ý ở tần số cao trong khi đó độ phân giải tần số cũng tùy ý ở tần thấp. Hay nói cách khác D f tỷ lệ với f : = C = constant (hằng số) Điều này rất có ý nghĩa, bởi vì trong thực tế phần lớn tín hiệu là tổng hợp của các thành phần tần cao trong khoảng thời gian ngắn và các thành phần tần thấp trong khoảng thời gian dài. Sóng điện tim là thí dụ điển hình về ý nghĩa này. III. WT và kỹ thuật tách đỉnh sóng điện tim : 3.1. Cơ sở của kỹ thuật tách điểm đặc trưng : Sóng điện tim bao gồm tập hợp các đỉnh đặc trưng, đó là : P, Q, R, S, T. Việc xác định vị trí của các đỉnh này là vô cùng quan trọng. WT được đánh giá là kỹ thuật thích hợp nhất, lý tưởng nhất trong việc phân tích các điểm đặc trưng của sóng ECG. Giá trị cực đại modul cục bộ của WT cung cấp thông tin để phân tích các đặc trưng quan trọng của tín hiệu ECG. Từ đó ta có thể phân biệt được các sóng ECG chìm trong nhiễu, tín hiệu giả và trôi đường đẳng trị. Hình 2 : Phân tích tần số bằng WT (a): Đồ thị wavelet y (t) (b): Đồ thị của y s1 (t) với s 1 >1 (c): Đồ thị của y s2 (t) với s 2 >1 (a’), (b’), (c’) biến đổi Fourier của (a), (b), (c) Tín hiệu ECG trước hết được làm trơn ở các thang khác nhau, sau đó các điểm đột biến được tách ra từ đường đạo hàm bậc 1 của chúng. Giả sử ta có q (x) là hàm làm trơn. Hàm này phải có tích phân bằng 1 và hội tụ tại O ở vô cùng. Ta ký hiệu y (a) (x) và y (b) (x) là đạo hàm bậc 1 và bậc 2 của hàm q (x) : y (a) (x) = (9) và y (a) (x) = (10) Theo định nghĩa, y (a) (x) và y (b) (x) có thể được coi như là Wavelet vì tích phân của chúng bằng O. y (a) (x) . dx = 0 (11) và y (b) (x) . dx = 0 (12) Phép giãn hàm y (x) sẽ là : y s (x) = y ( ) (13) WT được tính bởi phép nhân chập giữa tín hiệu và Wavelet đã được làm giãn. WT của f(x) ở thang s và vị trí x được tính với quan hệ Wavelet y (a) (x) và y (b) (x) được xác định bởi : W s a . f(x) = f(x) * y s (a) (x) (14) và W s b . f(x) = f(x) * y s (b) (x) (15) Hai hàm trên chính là đạo hàm bậc 1 và bậc 2 của tín hiệu đã được làm trơn ở thang s. Vị trí cực trị của W s a . f(x) tương ứng với các điểm cắt O của hàm y s (b) .f(x) và các điểm uốn của f*q s (x), nghĩa là chúng tương ứng với các điểm đột biến của sóng điện tim (xem hình 3). Hình 3 : Quan hệ giữa các biến thiên tín hiệu và WT của nó Điểm uốn của f * q (x) có thể là giá trị lớn nhất hoặc nhỏ nhất của giá trị tuyệt đối đạo hàm bậc 1. Giá trị cực đại là các điểm đột biến của f * q (x), còn giá trị cực tiểu tương ứng với những biến đổi chậm. Với toán tử đạo hàm bậc 2 sẽ khó phân biệt được hai loại cắt O này. Cần nói thêm rằng, các điểm cắt O cho chúng ta về thông tin vị trí, nhưng không phân biệt được các dao động biên độ nhỏ từ các điểm gián đoạn quan trọng. Khi tách giá trị cực đại cục bộ, các giá trị của | W s a f(x) | được ghi ở các vị trí lớn nhất, chúng thể hiện đạo hàm ở các điểm uốn. Các giá trị này rất cần thiết cho các luật ra quyết định các thông số quan trọng của sóng điện tim. 3.2. Phát triển thuật toán : 3.2.1. WT của tín hiệu điện tim : Hàm Gaussian được chọn làm hàm làm trơn q (x). Đạo hàm bậc 1 của nó chính là Wavelet. WT được tính với thang s = 2 j , jẻ z, z là tập nguyên. Vì vậy WT này có tên là WT nhị nguyên và được định nghĩa bằng tích chập sau : W 2 j f(x) = f(x) * Y 2 j (x) (16) Trong đó : Y 2 j = y ( ) Biến đổi Fourier của Wf (2 j , x) theo biến x sẽ là : f(2 j , w ) = (2 j ,w ) * (w ) (17) và ù (2 j ,w )ù 2 = 1 (18) ټ /font> là ký hiệu của phép biến đổi Fourier rời rạc. Toàn bộ trục tần số được phủ bởi phép dãn hàm Y (w ) với hệ số thang 2 j , j ẻ z. Trong thực tế chúng ta không thể tính WT cho toàn thang 2 j khi j biến thiên từ - Ơ đến + Ơ , mà phải giới hạn j hữu hạn và có độ mịn lớn hơn 0. Để chuẩn hoá ta chọn thang mịn bằng 1 và 2 j là thang lớn nhất. Vậy hàm f (x) sẽ được xác định bởi : ù (w )ù 2 = ù (2 j ,w )ù 2 (19) Từ (18) ta có : ù (w )ù 2 = 1 - ù (2 j ,w )ù 2 (20) ù (w )ù 2 = 1 (21) Bởi lẽ nǎng lượng của biến đổi Fourer của (w ) tập trung ở tần số thấp và đó là hàm trơn. Toán tử làm trơn S 2 j được định nghĩa bằng công thức sau: S 2 j f(x) = f(x) * F 2 j (x) (22) Trong đó: F 2 j (x) = . F 2 j (x/2 j ) Biến đổi Fourier của S 1 f(x), S 2 j f(x) và W s j f(x) sẽ là: 1 f(w ) = (w ) . (w ) (23) 2 j f(w ) = ( 2 j w ) . (w ) (24) 2 j f(w ) = ( 2 j w ) . (w ) (25) ứng dụng định lý Parseval ta có: ù ù 1 f(x) ù ù 2 = ù W 2 j f(x)ù 2 +ù S 2 j f(x)ù 2 (26) Phương trình này chứng tỏ rằng các thành phần tần cao của S 1 f(x) đã bị mất trong S 2 j f(x) có thể thu lại từ WT nhị nguyên giữa s = 1 và s = 2 j . Trong thực tế, tín hiệu cần xử lý thường đã được rời rạc hoá. Bất kỳ tín hiệu rời rạc có nǎng lượng hữu hạn nào đều có thể biểu diễn bằng phép lấy mẫu đơn trị của cùng một hàm đã được làm trơn với thang s = 1. Giả sử D = d n , nẻ z là tín hiệu rời rạc có nǎng lượng hữu hạn : ù d nù 2 < + (27) (Chú ý : Đây chính là điều kiện tồn tại phép biến đổi Fourier của tín hiệu rời rạc. Có thể tham khảo thêm ở 3 tập (1), (2), (3). Lúc này tồn tại một hàm f(x) với bất kỳ n là nguyên để : S 1 f(n) = d n Và tín hiệu rời rạc D của d n sẽ là : [...]... môđun này được tính trong cùng một cửa sổ thời gian c Tách sóng T và P 4 5 Nǎng lượng sóng T và P tồn tại phần lớn ở thang s = 2 và 2 Nhưng thông thường trôi đường 5 4 đẳng trị lại rất trầm trọng ở s = 2 Do vậy, thang s = 2 được chọn để tách sóng T và P Phép tách này cũng được thực hiện trong cùng một cửa sổ thời gian III KếT LUậN : Phép bi n đổi Wavelet đã và đang là công cụ rất mạnh trong kỹ thuật... nghiên cứu các phương pháp thử nghiệm chất lượng điện tim thế hệ mới, đó là một trong những nhu cầu cấp thiết trong lĩnh vực quản lý chất lượng thiết bị y tế WT là một lĩnh vực khá mới mẻ ở Việt Nam, vì vậy bài viết không tránh khỏi các sai sót Chúng tôi rất mong được sự góp ý và thông cảm của các bạn đồng nghiệp Chúng tôi xin chân thành cảm ơn PGS.TS Phạm Thượng Hàn - Chủ nhiệm Bộ môn kỹ thuật đo và. .. (29) j-1 f(n-2 k) (30) Trong đó: là toán tử đơn j 1 f(n) = dn là tín hiệu cần xử lý ở đây là tín hiệu tim j W2 f(x) là bi n đổi Walelet của f(x) hk với k ẻ z là các hệ số của bộ lọc thông thấp H(w ) gk với k ẻ z là các hệ số của bộ lọc thông cao G(w ) Nghĩa là: -ikw (31) -ikw (32) H(w ) = hk e G(w ) = gk e Khi Wavelet là dạng lưỡi toàn phương có trụ bền vứng và 1 mômen tắt, thì bi n đổi Fourier rời rạc... (hay đáp ứng xung hữu hạn FIR = Finite Impulse Response), có pha tuyến tính tổng quát Nó cũng là bộ lọc phản đối xứng và có độ trễ thời gian từ điểm trung tâm j-1 (trễ được xét bằng 2 -1 điểm trong thuật toán) Hình 4 : Quan hệ giữa các điểm đặc trưng của các sóng mô phỏng và WT của chúng ở các thang khác nhau Hình 4 cho chúng ta thấy rằng : tuỳ theo dạng lưỡi toàn phương của Wavelet, các sóng đơn... f(n) ở các thang khác nhau Sườn lên của sóng tương ứng với cực tiểu âm, và sườn xuống tương ứng với cực đại dương Các mođun của các cực đại, cực tiểu được gọi là đường cực đại môđun Nếu sóng đơn pha đối xứng qua đỉnh của nó (Hình 1(a)), thì đỉnh của nó tương ứng với điểm cắt j-1 O của cặp đại dương- tiểu âm với độ trễ chính xác là 2 điểm Nếu là phản đối xứng (hình 1b) thì j-1 đỉnh của nó tương ứng với... PGS.TS Phạm Thượng Hàn - Chủ nhiệm Bộ môn kỹ thuật đo và tin học công nghiệp - Đại học Bách khoa Hà Nội, PGS.TS Phạm Thượng Cát, Phân viện trưởng Phân viện Tự động hoá - Viện Công nghệ thông tin - Trung tâm Khoa học tự nhiên và công nghệ Quốc gia, đã cung cấp nhiều tài liệu quý giá và đóng góp nhiều ý kiến bổ ích trong quá trình nghiên cứu và bi n soạn tài liệu này TàI LIệU THAM KHảO 1 Phạm Thường Hàn:... Thang càng lớn thì sai số giữa trễ và 2 - 1 cũng càng lớn j Dựa vào các cơ sở nêu trên, ta có đáp ứng bi n tần của các bộ lọc tương đương Q (w ) với j = 1,2.; 5, theo tốc độ lấy mẫu 250/s như hình 5 j Hình 5 : Đáp ứng bi n tần của bộ lọc tương đương Q (w ) ở các thang khác nhau Đến đây, ta đã có đủ những cơ sở để xây dựng phương pháp tách điểm đặc trưng cho sóng điện tim 3.2.3 Phương pháp tách điểm đặc... thuật này mà độ chính xác, độ tin cậy của các hệ chẩn đoán ứng dụng trí tuệ nhân tạo đang được nâng cao Sau một thời gian nghiên cứu WT, chúng tôi đang chuẩn bị cơ sở kỹ thuật để tiến hành các thực nghiệm tại Bộ môn kỹ thuật đo và tin học công nghiệp - Trường Đại học Bách khoa Hà Nội Hy vọng rằng trong một thời gian gần đây chúng tôi sẽ có được những kết quả nghiên cứu Chúng tôi cũng tin tưởng rằng... loại bỏ nó * Định vị đỉnh R Từ WT của các tín hiệu khác nhau, ta dễ dàng nhận ra rằng bi n độ của WT phụ thuộc vào độ dốc của các thành phần tín hiệu đầu vào Phức hợp QRS gần như đối xứng, do đó bi n độ của cực tiểu âm và cực đại dương là gần như nhau Ngược lại nhiễu chứa các thành phần phi đối xứng Trong hầu hết các trường hợp, thuật toán đều tách chính xác đỉnh R Tuy nhiên ta có thể nâng cao chất... của sóng S (hoặc R khi S vắng) Sóng Q và S là các sóng tần cao, bi n độ thấp Nǎng lượng của chúng phân bố tập trung ở các 1 thang thấp Vì vậy điểm đầu và cuối của QRS được tách ở thang s =2 Điểm đầu QRS tương ứng với khởi điểm của cực đại môđun thứ nhất trước khi sóng R tạo ra cặp cực đại môđun Điểm cuối QRS tương ứng với điểm kết thúc của cực đại môđun thứ nhất sau khi sóng R tạo ra cặp cực đại môđun . BI N Đổi WAVELET Và ứNG DụNG TRONG CáC THIếT Bi GHI SóNG TIM CÔNG NGHệ Mới (Xử lý tín hiệu trong đo lường) Tóm tắt: Bài viết này giới thiệu phép bi n. (dilation) và đưa đến một phép bi n đổi mới đó là WAVELET TRANSFORM (WT). Phép bi n đổi này sẽ sử dụng Q không đổi, hay độ rộng dải tương đối không đổi. 2.3. Bi n