thtab63.m

来自「％ Atomizer Main Directory, Version .802 」· M 代码 · 共 23 行

% thtab63: BP Thesis Table 6.3 -- Barrier Iteraion Log of BPDN_Primal
% 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_Primal
%

help thtab63

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_Primal(y, 'WP', D, qmf, 0, 1e-10, 1e-3);