matlab wden nan,Automatic 1-D de-noising

Examples

% The current extension mode is zero-padding (see dwtmode).

% Set signal to noise ratio and set rand seed.

snr = 3; init = 2055615866;

% Generate original signal and a noisy version adding

% a standard Gaussian white noise.

[xref,x] = wnoise(3,11,snr,init);

% De-noise noisy signal using soft heuristic SURE thresholding

% and scaled noise option, on detail coefficients obtained

% from the decomposition of x, at level 5 by sym8 wavelet.

lev = 5;

xd = wden(x,'heursure','s','one',lev,'sym8');

% Plot signals.

subplot(611), plot(xref), axis([1 2048 -10 10]);

title('Original signal');

subplot(612), plot(x), axis([1 2048 -10 10]);

title(['Noisy signal - Signal to noise ratio = ',...

num2str(fix(snr))]);

subplot(613), plot(xd), axis([1 2048 -10 10]);

title('De-noised signal - heuristic SURE');

% De-noise noisy signal using soft SURE thresholding

xd = wden(x,'heursure','s','one',lev,'sym8');

% Plot signal.

subplot(614), plot(xd), axis([1 2048 -10 10]);

title('De-noised signal - SURE');

% De-noise noisy signal using fixed form threshold with

% a single level estimation of noise standard deviation.

xd = wden(x,'sqtwolog','s','sln',lev,'sym8');

% Plot signal.

subplot(615), plot(xd), axis([1 2048 -10 10]);

title('De-noised signal - Fixed form threshold');

% De-noise noisy signal using minimax threshold with

% a multiple level estimation of noise standard deviation.

xd = wden(x,'minimaxi','s','sln',lev,'sym8');

% Plot signal.

subplot(616), plot(xd), axis([1 2048 -10 10]);

title('De-noised signal - Minimax');

% If many trials are necessary, it is better to perform

% decomposition once and threshold it many times:

% decomposition.

[c,l] = wavedec(x,lev,'sym8');

% threshold the decomposition structure [c,l].

xd = wden(c,l,'minimaxi','s','sln',lev,'sym8');

% Editing some graphical properties,

% the following figure is generated.

Denoise a signal consisting of a 2-Hz sine wave with transients

at 0.3 and 0.72 seconds. Use Donoho and Johnstone's universal threshold

with level-dependent estimation of the noise. Obtain denoised versions

using the DWT and MODWT. Compare the results.

N = 1000;

t = linspace(0,1,N);

x = 4*sin(4*pi*t);

x = x - sign(t-.3)-sign(.72 - t);

y = x+0.15*randn(size(t));

xdDWT = wden(y,'sqtwolog','s','mln',3,'db2');

xdMODWT = wden(y,'modwtsqtwolog','s','mln',3,'db2');

subplot(2,1,1)

plot(xdDWT), title('DWT Denoising');

axis tight;

subplot(2,1,2)

plot(xdMODWT), title('MODWT Denoising');

axis tight;

Denoise a blocky signal using the Haar wavelet with MODWT and

DWT denoising. Compare the L2 and L-infty norms of the difference

between the original signal and the denoised versions.

[x,xn] = wnoise('blocks',10,3);

xdMODWT = wden(xn,'modwtsqtwolog','s','mln',6,'haar');

xd = wden(xn,'sqtwolog','s','mln',6,'haar');

plot(x)

hold on

plot(xd,'r--')

plot(xdMODWT,'k-.')

legend('Original','DWT','MODWT')

hold off

norm(abs(x-xd),2),

norm(abs(x-xd),Inf)

norm(abs(x-xdMODWT),2),

norm(abs(x-xdMODWT),Inf)