小波变换 阈值去噪

clear all
clc
load leleccum;%装在采集的信号leleccum.mat
s=leleccum（1:3920）;%将信号中第1到第3920个采样点赋给s
ls=length（s）;

%画出原始信号
subplot（221）;
plot（s）;
title（‘原始信号‘）;grid;

%用db1小波对原始信号进行3层分解并提取系数
[cL]=wavedec（s3‘db1‘）;
ca3=appcoef（cL‘db1‘3）;
cd3=detcoef（cL3）;
cd2=detcoef（cL2）;
cd1=detcoef（cL1）;

%对信号进行强制性去噪处理并图示结果
cdd3=zeros（1length（cd3））;
cdd2=zeros（1length（cd2））;
cdd1=zeros（1length（cd1））;
c1=[ca3 cdd3 cdd2 cdd1];
s1=waverec（c1L‘db1‘）;
subplot（222）;