Đang tải... (xem toàn văn)
Xử lý tín hiệu ECG bằng Wavelet
LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 26 2.4. CÁC PHƯƠNG PHÁP XỬ LÝ TÍN HIỆU TRƯỚC ĐÂY: Xử lý tín hiệu : - Trong miền thời gian: dùng bộ lọc. - Trong miền tần số: phương pháp Fourier Transform. - Trong miền thời gian và tần số: STFT( Short Time Fourier Transform). 2.4.1. Phương pháp Fourier: 2.4.1.1. Biến đổi Fourier: Cho một hàm f(t) khả tích tuyệt đối, biến đổi Fourier của nó được định nghĩa: () () ()tfedtetfFtitj,ωωω ==∫∞∞−− (2.14) Biến đổi Fourier ngược là: () ()∫∞∞−= ωωπωdeFtftj21 (2.15) Giả sử rằng biến đổi Fourier và Fourier ngược của nó tồn tại, ta kí hiệu : () ( )ωFtf = là một cặp biến đổi Fourier. ü Biến đổi liên tục: >=<⋅=⋅=∫∫∞∞−−∞∞−tjtjtjetxdetxjXdejXtxωωωωωωωπ),()()()(21)( (2.16) ü Rời rạc: njnjnjjenxeXeeXnxωωπωωπ−∞−∞=∑∫⋅=⋅=][)()(21][2 (2.17) PDF created with pdfFactory Pro trial version www.pdffactory.comwww.bme.vn LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 27 2.4.1.2. Tính chất: Dịch :Nếu trong miền thời gian nếu f(t) bị dịch đi một đoạn là t0 thì trong miền tần số biến đổi Fourier của nó được nhân với một hệ số pha: ( )00tjettfω−↔− Ngược lại, nếu trong miền tần số sẽ nhân một hệ số pha trong miền thời gian: ( )0)(0ωωω−↔ Ftfetj Tỷ lệ :Một tỷ lệ trong miền thời gian sẽ sinh ra một tỷ lệ nghịch trong miền tần số: ( )↔aFaatfω1 Moment : Gọi mn là moment thứ n của F(t): ()∫∞∞−= dttftmnn, n = 1,2 … 2.4.1.3. Chuỗi Fourier: Cho hàm tuần hoàn F(t) với chu kỳ T: F(t+T) = F(t) Nó có thể được biểu diễn là một tổ hợp tuyến tính của các số mũ phức với tần số nω0, trong đó ω0 = 2π/T. Khai triển chuỗi Fourier của F(t) () ()∑∞−∞==ktjkekFtf0ω (2.18) với () ()∫−−=2201TTtjkdtetfTkFω • Phân tích liên tục: PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 28 ><=⋅=⋅===∫∫∑∑−−∞−∞=∞−∞=tjkTtTjkTtjkkktTjkkktjkketxTdtetxTdtetxTaeaeatx000),(1)(1)(1)(22ωπωπω (2.19) • Phân tích rời rạc: ><=⋅=⋅===−>=<−>=<−>=<>=<∑∑∑∑njkNnnNjkNnnjkkNknNjkkNknjkkenxNenxNenxNaeaeanx000],[1][1][1][22ωπωπω (2.20) Nhược điểm : Khi biến đổi sang miền tần số thông tin thời gian bị mất. Nếu thuộc tính tín hiệu không thay đổi nhiều theo thời gian thì nhược điểm là không đáng kể. Tuy nhiều tín hiệu có chứa các thông số động: trôi, nghiêng, biến đổi đột ngột, khởi đầu và kết thúc của sự kiện thì Fourier không phát hiện được. 2.4.2. Phương pháp STFT: Để đạt được một biến đổi Fourier cục bộ, chúng ta có thể định nghĩa một biến đổi Fourier cửa sổ. Tín hiệu đầu vào được nhân với một hàm cửa sổ W(t - τ) và sau đó lấy biến đổi Fourier của nó. Kết quả là một biến đổi hai chỉ số STFTf(ω,τ) được cho bởi: [ ]∫−⋅−⋅=ttjxdtetWtxSTFTωωτωτ )()(),( (2.21) )(τ−tW : hàm cửa sổ dùng để phân tích (tập trung tại τ=t) tjeω−: hàm cơ bản (trọng số của FT) Khi đó ta thể biểu diễn dưới dạng tích trong: ( ) ( ) ( )tftgSTFTf τωτω,, = với ( ) ( )tjettgωτωτω −=, (2.22) PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 29 Mỗi hàm thành phần trong khai triển có cùng độ phân giải thời gian và tần số, một cách đơn giản là một sự định vị khác nhau ở miền thời gian – tần số. f(t) có thể được khôi phục lại trong L2(R) bằng tích phân kép () ( ) ()∫∫∞∞−∞∞−= τωτωπτωddtgSTFTtff ,,21 (2.23) STFT cũng có tính chất bảo toàn năng lượng: () ( )∫∫∞∞−∞∞−= τωτωπddSTFTtff22,21 (2.24) Cung cấp phân bố năng lượng cả t1in hiệu trong mặt phẳng thời gian tần số 2202)()()()(),(),(∫−∗⋅−⋅===tftjxxSdtetttxftSTFTftSPECPπγγγ (2.25) Phân tích rời rạc : dtenTttxkFnTSTFTkFtjtxπγγ2*)()()(),(−∫−= (2.26) ∑∑−⋅=nkkFtjxenTtgkFnTSTFTtxπγ 2)()(),()( (2.27) Cửa sổ hẹp, phân giải thời gian tốt Cửa sổ rộng phân giải tần số tốt Nhược điểm: Độ chính xác giới hạn phụ thuộc vào kích thước hàm cửa sổ. Chọn kích thước cụ thể cho cửa sổ thời gian, bằng nhau cho mọi tần số. Đối với một số tín hiệu cần sự mềm dẻo hơn thì STFT không đáp ứng được PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 30 2.5. BIẾN ĐỔI WAVELET: 2.5.1. Phân tích đa phân giải: Phân tích tín hiệu ra các thành phần dựa trên dãy xấp xỉ. Một tín hiệu sẽ được xấp xỉ thô cộng với các chi tiết. Trong đó không gian con thô và không gian con chi tiết là trực giao với nhau. Nói cách khác, tín hiệu chi tiết là hiệu của phiên bản thô và phiên bản tinh của tín hiệu . Bằng cách áp dụng một cách đệ quy các dãy xấp xỉ, chúng ta sẽ thấy không gian các tín hiệu đầu vào L2(R) có thể được sinh ra bởi các không gian của dãy xấp xỉ tại tất cả các độ phân giải. Khi độ phân giải chi tiết tiến đến vô cùng thì sai số xấp xỉ tiến đến 0. Độ phân giải được đưa ra bởi Mallat và Meyer, nó không chỉ là cơ sở cho Wavelet mà còn là công cụ toán học rất mạnh để liên kết Wavelet và phân tích băng con tín hiệu. Định nghĩa : Một phân tích đa phân tích đa phân giải bao gồm một chuỗi của khối các không gian con đóng. 2112 −−⊂⊂⊂⊂ VVVVVo Sao cho : - Đầy đủ hợp : U−−∈=ZmmRLV )(2 - Đầy đủ giao : { }IZmmV∈= 0 - Bất biến tỷ lệ : ommVtfVtf ∈⇔∈ )2()( - Bất biến dịch : ZnVntfVtfoo∈∀∈−⇒∈ ,)()( - Tồn tại oV∈∞ϕ sao cho { }znnt ∈− )(ϕ là một cơ sở trực chuẩn Vo. 2.5.2. Chuỗi wavelet: Một hàm f ∈ L2(R) có thể được biểu diễn: PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 31 Công thức phân tích ( ) ( ) ( )∑∈=ZnmnmtnmFtf,,, ψ (2.28) Công thức tổng hợp ( ) () () ()()∫∞∞−== dttfttftnmFnmnm ,,,, ψψ (2.29) • Tính chất: - Tuyến tính: Giả sử rằng toán tử T được định nghĩa: ( ){ } ( ) ( ) ( )tftnmFtfTnm,,,ψ== Với a, b ∈ R thì: ( ) ( )[ ]( )[ ]( )[ ]tgbTtfaTtbgtafT +=+ - Dịch : Nếu một tín hiệu có một khai triển tỷ lệ hữu hạn () ( ) ()∑∑∈ −∞==ZnMmnmtnmFtf2,, ψ Thì tín hiệu này sẽ có tính dịch yếu tương ứng được dịch đi kM22, đó là : ( ) ( )knmFktfmMM −−↔−222,2 với - ∞ ≤ m ≤ M2. - Tỷ lệ: Nếu tín hiệu f(t) có hệ số biến đổi Wavelet là F(m,n) thì: PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 32 ()( )ZknkmFtfkk∈−↔−,,222 - Đẳng thức Pareval. Họ Wavelet trực chuẩn { })(,tnmψ thỏa mãn ()2,2,ftZnmnm=∑∈ψ f ∈ L2(R). Bảo toàn năng lượng - Lấy mẫu đôi và Tiling thời gian - tần số: Quá trình lấy mẫu ở miền thời gian, tại tỷ lệ m, được thực hiện với chu kỳ 2m, lúc đó ψm,n(t) = ψm,0(t –2mn). Ở dạng tỷ lệ, số mũ của 2 thường được đề cập. Khi tần số là đảo của tỷ lệ thì ta thấy nếu Wavelet tập trung quanh ω0 thì Ψm,n(ω) tập trung quanh m20ω. Điều này sinh ra một quá trình lấy mẫu đôi của miền thời gian - tần số Lấy mẫu đôi của miền thời gian Tiling của miền thời gian tần số. tần số của khai triển chuỗi Wavelet. Hình 2.10: Lấy mẫu đôi và tiling thời gian. Shift nScale mm = 2m = 1m = 0m = -2m = -2tfPDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 33 2.5.3. Các tính chất của hàm tỷ lệ: 2.5.3.1. Phương trình hai tỷ lệ: Hàm tỷ lệ có thể được xây dựng từ chính nó. Phương trình hai tỷ lệ: () ()( )∑−=nntngt 220ϕϕ (2.30) Tương tự đối với Wavelet ()10VWt ⊂∈ψ thì: () ()( )∑−=nntngt 221ϕψ (2.31) 2.5.3.2. Tính chất moment: Bộ lọc thông thấp g0(n) trong phần dãi lọc, nó có ít nhất một zero tại ω = π, và vì vậy g1(n) có ít nhất một zero tại ω = 0. Khi Φ(ω) = 1, thì nó kéo theo ψ(ω) có ít nhất một zero tại ω = 0, do đó: () ()00=Ψ=∫∞∞−dttψ Một cách tổng quát: ()0=∫∞∞−dtttnψ n = 0, …, N –1. Điều đó có nghĩa N moment đầu tiên của Wavelet là zero. 2.5.4. Biến đổi wavelet : Một đặc điểm nỗi bậc của khai triển Wavelet là dựa trên nền tảng cấu trúc đa phân giải để đưa ra một giải thuật rời rạc thời gian hiệu quả, bằng cách thực hiện một dãi lọc, giải thuật này được đưa ra bởi Mallat và được gọi là giải thuật Mallat. 2.5.4.1.Biến đổi wavelet liên tục: § Wavelet mẹ: )(1)(,abtatba−=ψψ (2.32) PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 34 a: thông số dịch chuyển b: thông số tỉ lệ a, b ∈ R (a ≠ 0) và chuẩn ( ) ( )ttbaψψ =, a1:yếu tố bình thường hóa đảm bảo cho các wavelet có cùng mức năng lượng § Khai triển: Cho trước hàm wavelet mẹ )(tψ có giá trị thực, một hàm bất kì )()(2RLtf ∈ : - Công thức phân tích: )(],[)(,,tnmFtfnmZnmψ∑∈= (2.33) - Công thức tổng hợp:dttfttftnmFnmnm)()()(),(],[,,∫∞∞−== ψψ (2.34) ],[ nmF được gọi là các hệ số khai triển chuỗi wavelet 2.5.4.2 . Biến đổi rời rạc: § Chọn các giá trị cố định: 0maa = và 00 manbb = , .2,1,0, ±±=nm dtnbtatxadtttxnmWxmmnm)()()()(),(002/0,−==∫∫+∞∞−−−+∞∞−ψψ (2.35) Sử dụng phân tích đa phân giải x(t) được phân tách thành nhiều mức khác nhau: x(t)= ∑=Kj 1∑∞−∞=kdi(k)ψj,k(t)+ ∑∞−∞=kaK(k)φK,k(t) (2.36) ψj,k(t): wavelet rời rạc dung phân tích φK,k(t) là hàm tỉ lệ rời rạc dj(k): tín hiệu chi tiết (hệ số wavelet) tại mức 2j PDF created with pdfFactory Pro trial version www.pdffactory.com LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm2007 35 aK(k) tín hiệu xấp xỉ ( hệ số tỉ lệ) tại mức 2K § Cơ sở wavelet rời rạc gồm hai hàm cơ bản: -Trường hợp cơ sở trực chuẩn: Cơ sở phân tích tín hiệu trùng với cơ sở tổng hợp tín hiệu. Cơ sở gồm hai hàm cơ bản: ][0nφ và ][1nφ - Trường hợp trường hợp cặp cơ sở trực giao: Cơ sở phân tích tín hiệu khác với cơ sở tổng hợp tín hiệu. Cơ sở phân tích gồm hai hàm cơ bản: ][0nφ và ][1nφ Cơ sở tổng hợp gồm hai hàm cơ bản: ][~0nφ và ][~1nφ § Dãy bộ lọc hai kênh: Hai hàm cơ sở của wavelet rời rạc chính là những hệ số của hai bộ lọc thông thấp và thông cao h[n]=)2(),(21ntt −φφ (2.37) g[n]=)1()1()2(),(21nhntt −−=−φψ (2.38) 2.5.5. Biến đổi chuỗi wavelet của tín hiệu : Wavelet có ba hàm cơ sở. Hiện nay người ta quan tâm nhiều đến nghiên cứu hai loại là wavelet trực chuẩn (như wavelet Haar, wavelet Daubechies…)và wavelet cặp trực giao ( wavelet B-spline). Biến đổi chuỗi wavelet gồm hai qua trình: khai triển tín hiệu thành những chuỗi wavelet và tổng hợp chúng thành tí hiệu ban đầu. Biến đổi chuỗi wavelet thực hiện bằng cơ sở trực chuẩn hay cặp trực giao đều có sơ đồ chung chỉ khác ở hàm cơ sở tổng hợp tín hiệu. Thực hiện dựa trên biến đổi wavelet rời rạc do phân tích rời rạc đảm bảo mã hóa tiế kiệm không gian và là vừa đủ cho sự tái tạo chính xác. Tín hiệu vào x[n] được phân tích bởi dãy bộ lọc băng hai kênh thành hai dãy hệ số co kí hiệu là cD và cA. Dãy hệ số của nhánh lọc thông thấp cA, lại tiếp tục được phân tích PDF created with pdfFactory Pro trial version www.pdffactory.com [...]...TRƯỜNG ĐHBK TP.HCM-năm2007 LUẬN VĂN TỐT NGHIỆP bởi cùng một dãy bộ lọc hai kênh, và sự phân tích lại tiếp tục trên nhánh thông thấp Mỗi lần khai triển sẽ tạo ra hai dãy hệ số wavelet : - Bộ lọc thông cao tạo ra những hệ số chi tiết - Bộ lọc thông thấp tạo ra các hệ số xấp xỉ Sơ đồ rút gọn của phân tích các tầng: Hình 2.1: : Sơ đồ phân tích DWT Sơ đồ tái tạo các tầng: Hình . TRƯỜNG ĐHBK TP.HCM-năm2007 26 2.4. CÁC PHƯƠNG PHÁP XỬ LÝ TÍN HIỆU TRƯỚC ĐÂY: Xử lý tín hiệu : - Trong miền thời gian: dùng bộ lọc. - Trong miền. TP.HCM-năm2007 30 2 .5. BIẾN ĐỔI WAVELET: 2 .5. 1. Phân tích đa phân giải: Phân tích tín hiệu ra các thành phần dựa trên dãy xấp xỉ. Một tín hiệu sẽ được xấp