example6_8.m

关于小波分析的 具体是什么 我也不太清楚

%装载原始信号
load sumsin;
s=sumsin;
%==============================
%设置小波名并利用coif3小波进行4层分解
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(l(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('修改后的小波分解系数');