clear;?
%讀取ECG信號?
fid=fopen（‘ecg.txt‘‘r‘）;?
x=fscanf（fid‘%f?%f?%f‘[600?3]）;?
fclose（fid）;?
N=length（x）;?
for?i=0:N/3-1?
????t（i+1）=x（3*i+1）;?
????sig1（i+1）=x（3*i+2）;?
end?
s=sig1+0.6*randn（1length（sig1））;?
?
figure（1）;
subplot（211）;
plot（ts‘LineWidth‘2）;
xlabel（‘時間‘）;
axis（[0?1.2?-inf?inf]）;
ylabel（‘幅值‘）;?
title（‘原始ECG信號‘）;?
grid?on

%選用基波函數?
wavelet=‘db5‘;?
%分解級數?
level=3;
%壓縮參數（去噪時alpha=3）
alpha=1.5;?
%閾值類型（硬閥值）
sorh=‘h‘;?
%小波分解?
[cl]=wavedec（slevelwavelet）;?
%birge-massart算法計算閾值?
[thrnkeep]=wdcbm（clalpha）;?
%使用硬閾值壓縮信號?
%‘lvd’為允許設置各層閥值，‘gbl閥值’為固定閥值
[xccxclxcperf0perfl2]=wdencmp（‘lvd‘clwaveletlevelthrsorh）;?
subplot（212）;
plot（txc‘LineWidth‘2）;
xlabel（‘時間‘）;
axis（[0?1.2?-inf?inf]）;
ylabel（‘幅值‘）;
title（‘去噪ECG信號‘）;
grid?on

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

?????文件????????893??2011-01-05?17:34??ecg1.m

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

??????????????????893????????????????????1