Xử lý tín hiệu ECG bằng Wavelet phần 8

LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm 2007


78

Phụ lục

1. Các họ wavelet được sử dụng
• Haar

• Daubechies Wavelets: dbN


• Biorthogonal:

PDF created with pdfFactory Pro trial version www.pdffactory.com
www.bme.vn
LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm 2007


79


• Coiflets



• Symlets

PDF created with pdfFactory Pro trial version www.pdffactory.com
LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm 2007

80


2. Các lệnh sử dụng:
§ Awgn
Mục đích: Cộng nhiễu trắng Gaussian vào tín hiệu
Cú pháp: y =awgn(x,snr) : cộng nhiễu trắng Gaussian vào vector tín hiệu x
y=awgn(x,snr,sigpower) : sigpower là công suất của x tính theo
dbW y=awgn(x,snr,’measured’): đo công suất trước khi cộng
nhiễu
y=awgn(x,snr,state): khởi động lại trạng thái của bộ phát ngẫu
nhiên rand về trạng thái số nguyên
y=awgn(…,powertype): chỉ ra đơn vị của snr và sigpower
§ Dwt Biến đổi wavelet rời rạc đơn mức
Mục đích: biến đổi wavelet 1 chiều rời rạc đơn mức
Cú pháp: [cA,cD]=dwt(X,’wname’): tính toán vecto hệ số xấp xỉ cA và
chi tiết cD đạt từ phân tách wavelet của vecto X
[cA,cD]=dwt(X,’wname’,’mode’,MODE)
[cA,cD]=dwt(X,Lo_D,Hi_D): tính phân tách wavelet với các
bộ lọc phân tách Lo_D thông thấp, bộ lọc phân tách Hi_D thông cao
[cA,cD]=dwt(X,Lo_D,Hi_D,’mode’,MODE)
Mô tả: Thực hiện phân tách wavelet một chiều đơn mức với một wavelet
nhất định
MODE chuỗi chứa chế độ mở rộng yêu cao, có thể tự chọn
§ Wavedec
PDF created with pdfFactory Pro trial version www.pdffactory.com
LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm 2007

81
Mục đích: Phân tách wavelet một chiều đa mức
Cú pháp: [C,L]=wavedec(X,N,’wname’)
[C,L]=wavedec(X,N,Lo_D,Hi_D)
§ Wmaxlev
Mục đích: Mức phân tách tối đa
Cú pháp: L=wxmaxlev(S,’wname’)
Mô tả:wmaxlev cho phân tách mức cho phép tối đa, nhưng khái quát giá
trị nhỏ hơn được nhân Giá trị thông thường là 5 cho trường hợp 1 chiều
Các hàm tổng hợp và tái tạo
§ Idwt
Mục đích:Biến đổi wavelet 1 chiều rời rạc ngược đơn mức
Cú pháp: X=idwt(cA,cD,’wname’)
X=idwt(cA,cD,Lo_R,Hi_R)
X=idwt(cA,cD,’wname’,L)
X=idwt(cA,cD,Lo_R,Hi_R,L)
X=idwt(..,’mode’,MODE)
§ Waverec
Mục đích: Tái tạo wavelet 1 chiều đa mức
Cú pháp: X=waverec(C,L,’wname’)
X=waverec(C,L,Lo_R,Hi_R)
§ Wrcoef
Mục đích : Tái tạo nhánh đơn từ các hệ số wavelet một chiều
Cú pháp: X=wrcoef(‘type’,C,L,’wname’,N): tái tạo mức N
X= wrcoef(‘type’,C,L,Lo_R,Hi_R,N):bộ lọc tái tạo mức N
X= wrcoef(‘type’,C,L,’wname’): tái tạo mức tối đa
N=length(L)-2
X= wrcoef(‘type’,C,L,Lo_R,Hi_R)
Mô tả: tái tạo các hệ số của tín hiệu một chiều

Type: xác định chỉ số tái tạo chi tiết hoặc xấp xỉ
§ Upcoef
PDF created with pdfFactory Pro trial version www.pdffactory.com
LUẬN VĂN TỐT NGHIỆP TRƯỜNG ĐHBK TP.HCM-năm 2007


82
Mục đích:Tái tạo trực tiếp từ các hệ số wavelet 1 chiều
Cú pháp: Y=upcoef(O,X,’wname’,N)
Y=upcoef(O,X,’wname’,N,L)
Y=upcoef(O,X,Lo_R,Hi_R,N)
Y=upcoef(O,X,Lo_R,Hi_R,N,L)
Y=upcoef(O,X,’wname’)
Y=upcoef(O,X,’wname’,N)
Mô tả: O xác định loại chi tiết ‘d’ hay xấp xỉ ‘a’
Các tiện ích cấu trúc phân tách
§ Detcoef
Mục đích:Tách ra các hệ số chi tiết một chiều
Cú pháp D=detcoef(C,L,N)
D=detcoef(C,L)
Mức N cần 1≤N≤NMAX NMAX=length(L)-2
§ Appcoef
Mục đích: Tách ra các hệ số chi tiết một chiều
Cú pháp: A=appcoef(C,L,N)
A=appcoef(C,L)
§ Upwlev
Mục đích: tái tạo đơn mức của phân tách wavelet 1 chiều
Cú pháp: [NC,NL,cA]=upwlev(C,L,’wname’)
[NC,NL,cA]=upwlev(C,L,Lo_R,Hi_R)
Mô tả: thực hiện sự tái tạo đơn mức của cấu trúc phân tách[C,L] cho giá

trị mới [NC,NL] và trích ra vector các hệ số của xấp xỉ cuối cA
Mức phân tách N=length(L)-2
Khử nhiễu
§ Ddencmp
Mục đích:Các giá trị mặc định để khử nhiễu hoặc nén

PDF created with pdfFactory Pro trial version www.pdffactory.com