wavelet_constrain.asv

小波变换的条件限制的实现的编程

%利用小波对信号的高频进行抑制
%
%利用coif3小波进行4层分解
load sumsin;
s=sumsin;

w='coif3';
maxlev=4;
[c,l]=wavedec(s,maxlev,w);
newc=c;
%将分解后的第三层第四层细节系数设置为0
newc=wthcoef('d',c,l,[3,4]);
%在原始信号的时间区间【400,600】内将第一层细节系数设置为0
%并且将其他系数进行衰减，求出系数的起始点和中止点的索引值
k=maxlev+1;
first=sum(l(1:k-1))+1;
last=first+l(k-1);
indd1=first:last;
%将系数除以3，进行信号衰减
newc(indd1)=c(indd1)/3;

%在区间【400,600】上求出第一层的系数索引
indd1=(first+400/2):(first+600/2);

%将改索引值设置为0
newc(indd1)=zeros(size(indd1));
%将第二层的相应与原始信t=500的时间点处的位置为4

k=maxlev;
first=sum((1:k-1))+1;
newc(first+500/2^2)=4;
%
%综合修改后的分解结构
synth=waverec(newc,l,w);
%用图示所示
subplot(2,2,1);
plot(s);
title('原始信号');
subplot(2,2,2);
plot(c);
title('coif3小波分解后的系数')
subplot(2,2,3);
plot(synth);
title('小波抑制后的信号');
subplot(2,2,4);
plot(newc);
title('修改后的小波系数');