%裝載原始信號
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('修改后的小波分解系數');