軟閾值與硬閾值

　　cxdsoft=c;

　　cxdhard=c;

　　for j=1:jN j是分解尺度

　　flk = first（j）：last（j）;flk是di在c中的索引

　　thr（j）=sqrt（log（length（flk）））/log（j+1）;

　　for k=0：（length（flk）-1）k是位移尺度

　　djk=c（first（j）+k）;為了簡化程序

　　djk2=c（first（j）+k）;

　　absdjk=abs（djk）;

　　thr1=thr（j）;

　　閾值處理

　　if absdjkthr1

　　djk=djk-thr1/exp（（djk-thr1）/n）;

　　else

　　djk=djk+thr1/exp（（-djk-thr1）/n）;

　　end

　　cxd51（first（j）+k）=djk;

　　end

　　end

　　新方法二

　　cxd52=c;

　　for j=1:jN

　　flk = first（j）：last（j）;

　　for k=0：（length（flk）-1）

　　djk=c（first（j）+k）;為了簡化程序

　　absdjk=abs（djk）;

　　thr1=thr（j）;

　　if absdjkthr1

　　djk=djk-thr1+thr1/（2*b+1）;

　　else

　　djk=djk+thr1-thr1/（2*b+1）;

　　end

　　cxd53（first（j）+k）=djk;

　　end

　　end

　　新方法四

　　cxd54=c;

　　for j=1:jN

　　flk = first（j）：last（j）;

　　for k=0：（length（flk）-1）

　　djk=c（first（j）+k）;為了簡化程序

　　absdjk=abs（djk）;

　　thr1=thr（j）;

　　bb=（2*aa+1+m*sqrt（（2*aa+1）.^2-pi.^2））/（pi*thr1.^（2*n+1））;

　　kk=（1+（bb*thr1.^（2*n+1））.^2）/（（2*a+1）*bb*thr1.^（2*n））;

　　if absdjkthr1

　　djk=djk-thr1+kk*atan（bb*djk.^（2*aa+1））;

　　else

　　djk=djk+thr1-kk*atan（bb*djk.^（2*aa+1））;

　　end

　　cxd54（first（j）+k）=djk;

　　end

　　end

　　信號重構

　　snew51=waverec（cxd51，l，wname）;

　　snew52=waverec（cxd52，l，wname）;

　　snew53=waverec（cxd53，l，wname）;

　　snew54=waverec（cxd54，l，wname）;

　　信噪比及均方差分析

　　snry51=20*log10（norm（x）/norm（snew51-x））;

　　snry52=20*log10（norm（x）/norm（snew52-x））;

　　snry53=20*log10（norm（x）/norm（snew53-x））;

　　snry54=20*log10（norm（x）/norm（snew54-x））;

　　mmsey51=mmse（snew51-x）;

　　mmsey52=mmse（snew52-x）;

　　mmsey53=mmse（snew53-x）;

　　mmsey54=mmse（snew54-x）;

　　閾值處方法六

　　m=40;

　　cxd6=c;

　　for j=1:jN

　　flk = first（j）：last（j）;

　　thr1=thselect（s，’heursure‘）;

　　for k=0：（length（flk）-1

　　）

　　djk=c（first（j）+k）;為了簡化程序

　　absdjk=abs（djk）;

　　thr2=thr1/3;

　　if absdjk》thr1

　　djk=djk’*（absdjk.^m-thr1.^m）/absdjk.^m;

　　elseif djk