clc;
load?‘ECG_100.MAT‘?;
E=sig（:2）;
E=E‘;
n=size（E）;
s=E（1:2000）;%取部分數據
%wavelet
[cl]=wavedec（E3‘db5‘）;
ca3=appcoef（cl‘db5‘3）;
cd1=detcoef（cl1）;
cd2=detcoef（cl2）;
cd3=detcoef（cl3）;

%使用stein的無偏似然估計原理進行選擇各層的閾值
%cD1cD2cD3為各層小波系數，?
%‘rigrsure’為無偏似然估計閾值類型
thr1=thselect（cd1‘rigrsure‘）;
thr2=thselect（cd2‘rigrsure‘）;
thr3=thselect（cd3‘rigrsure‘）;
%各層的閾值
TR=[thr1thr2thr3];
%‘s‘為軟閾值;‘h‘硬閾值。
SORH=‘s‘;
%---------去噪----------------
%XC為去噪后信號
%[CXCLXC]為的小波分解結構
%PERF0和PERF2是恢復和壓縮的范數百分比。
%‘lvd‘為允許設置各層的閾值
%‘gbl‘為固定閾值。
%3為閾值的長度
[XCCXCLXCPERF0PERF2]=wdencmp（‘lvd‘E‘db5‘3TRSORH）;
%---------去噪效果衡量（SNR越大效果越好MSE越小越好）
%選取信號的長度
N=n（2）;
x=E;
y=XC;
F=0;
M=0;
for?ii=1:N
???m（ii）=（x（ii）-y（ii））^2;
???t（ii）=y（ii）^2;
???f（ii）=t（ii）/m（ii）;
????F=F+

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????240184??2001-05-14?16:43??ECG?denoising?based?on?wavelet?transform\ECG_100.MAT

?????文件???????1187??2019-03-14?18:21??ECG?denoising?based?on?wavelet?transform\ECG_wavelet.m

?????文件???????9119??2019-03-12?19:59??ECG?denoising?based?on?wavelet?transform\locationP_QRS_T.m

?????目錄??????????0??2019-03-12?19:55??ECG?denoising?based?on?wavelet?transform

-----------?---------??----------?-----??----

???????????????250490????????????????????4