%應(yīng)用db5作為小波函數(shù)進(jìn)行3層分解
%利用無(wú)偏似然估計(jì)閾值
%對(duì)100.dat?from?MIT-BIH-DB的單導(dǎo)聯(lián)數(shù)據(jù)進(jìn)行去噪處理
clear
clc
load（‘fangzhen.txt‘）;
fangzhen=fangzhen（:1）;
A=0.04:0.04:10.68;
fangzhen=fangzhen‘;
n=size（fangzhen）;
s=fangzhen（1:267）;
%小波分解
[C?L]=wavedec（fangzhen4‘db5‘）;%?從c中提取尺度3下的近似小波系數(shù)?第一部分C爲(wèi)小波分解係數(shù)（是一個(gè)矩陣形式），而L則是表示相應(yīng)小波係數(shù)的個(gè)數(shù)（矩陣內(nèi)係數(shù)的個(gè)數(shù)）。
cA3=appcoef（CL‘db5‘4）;%從信號(hào)c中提取尺度CA3近似部分，低頻系數(shù)
cD1=detcoef（CL1）;
cD2=detcoef（CL2）;
cD3=detcoef（CL3）;
cD4=detcoef（CL4）;
%從信號(hào)c中提取尺度123下的細(xì)節(jié)小波系數(shù)

%使用stein的無(wú)偏似然估計(jì)原理進(jìn)行選擇各層的閾值
%cD1cD2cD3為各層小波系數(shù)，?
%‘rigrsure’為無(wú)偏似然估計(jì)閾值類型
thr1=thselect（cD1‘rigrsure‘）;
thr2=thselect（cD2‘rigrsure‘）;
thr3=thselect（cD3‘rigrsure‘）;
thr4=thselect（cD4‘rigrsure‘）;
%各層的閾值
TR=[thr1thr2thr3thr4];
%‘s‘為軟閾值;‘h‘硬閾值。
SORH=‘s‘;
%---------去噪----------------
%XC為去噪后信號(hào)
%[CXCLXC]為的小波分解結(jié)構(gòu)
%PERF0和PERF2是恢復(fù)和壓縮的范數(shù)百分比。
%‘lvd‘為允許設(shè)置各層的閾值
%‘gbl‘為固定閾值。
%3為閾值的長(zhǎng)度
[XCCXCLXCPERF0PERF2]=wdencmp（‘lvd‘fangzhen‘db5‘4TRSORH）;
%使用