非抽取小波图像去噪

稀疏图像与信号处理的一个数值实验，利用非抽取小波进行图像去噪，对比正交小波和冗余小波的去噪效果。里面是完整的代码，用MATLAB可以运行。

function y = aconv（fx） 
% aconv -- Convolution Tool for Two-Scale Transform 
%  Usage 
%    y = aconv（fx） 
%  Inputs 
%    f    filter 
%    x    1-d signal 
%  Outputs 
%    y    filtered result 
% 
%  Description 
%    Filtering by periodic convolution of x with the 
%    time-reverse of f. 
% 
%  See Also 
%    iconv UpDyadHi UpDyadLo DownDyadHi DownDyadLo 
% 
 
	n = length（x）; 
	p = length（f）; 
	if p < n 
	   xpadded = [x x（1:p）]; 
	else 
	   z = zeros（1p）; 
	   for i=1:p 
		   imod = 1 + rem（i-1n）; 
		   z（i） = x（imod）; 
	   end 
	   xpadded = [x z]; 
	end 
	fflip = reverse（f）; 
	ypadded = filter（fflip1xpadded）; 
	y = ypadded（p:（n+p-1））; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .CA....       697  2018-12-05 12:36  My Project\aconv.m

    .CA....       352  2018-12-06 11:25  My Project\CircularShift.m

    .CA....       405  2018-12-05 12:38  My Project\DownDyadHi.m

    .CA....       384  2018-12-05 12:35  My Project\DownDyadLo.m

    .CA....    263222  2018-12-05 10:53  My Project\einstein.bmp

    .CA....       884  2018-12-05 12:34  My Project\FWT2_PO.m

    .CA....      2385  2018-12-05 12:33  My Project\FWT2_TI.m

    .CA....       240  2018-12-06 12:30  My Project\HardThresh.m

    .CA....       651  2018-12-05 12:40  My Project\iconv.m

    .CA....       989  2018-12-06 12:33  My Project\IWT2_PO.m

    .CA....      1640  2018-12-06 12:32  My Project\IWT2_TI.m

    .CA....       257  2018-12-05 12:39  My Project\lshift.m

    .CA....      8720  2018-12-05 12:23  My Project\MakeONFilter.m

    .CA....       397  2018-12-05 12:38  My Project\MirrorFilt.m

    .CA....       402  2018-12-05 12:33  My Project\packet.m

    .CA....       860  2018-12-05 22:43  My Project\ParseInputs.m

    .CA....       629  2018-12-05 12:19  My Project\quadlength.m

    .CA....      1704  2018-12-06 12:29  My Project\RedundantWavelet.m

    .CA....       271  2018-12-05 12:37  My Project\reverse.m

    .CA....       280  2018-12-05 12:59  My Project\rshift.m

    .CA....       382  2018-12-06 13:57  My Project\UpDyadHi.m

    .CA....       342  2018-12-06 13:52  My Project\UpDyadLo.m

    .C.D...         0  2018-12-15 15:21  My Project

----------- ---------  ---------- -----  ----

               286093                    23