thtab64.m

％ Atomizer Main Directory, Version .802 里面信号含有分解去噪合成过程的代码 ％---------------------------------------

% thtab64: BP Thesis Table 6.4 -- Barrier Iteraion Log of BPDN_Interior
% Signal: Noisy Bumps, SNR = 7
% Signal Length: 1024
% Dictionary: Wavelet Packets with D=log2(n), qmf = (Symmlet, 8)
% Observation:
%	Table 6.3 records the monitoring output of BPDN_Interior
%

help thtab64

n = 1024; D = log2(n);
m = 2*n * (D+1);
qmf = MakeONFilter('Symmlet', 8);

disp('Bumps');
x = MakeSignal('Bumps', n)';
SNR = 7;
randn('seed',931316785);
[x y] = NoiseMaker(x, SNR);

[xrec c] = BPDN_Interior(y, 'WP', D, qmf, 0, 1e-5, 1e-5, 1e-3, 1e-4);