bpfig23.m

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

% bpfig23: BP Figure 2.3 -- MOF representations are not sparse
%----------------------------------------------------------------------
%
%    MOF tends to use all the basis functions nonorthogonal to
%    signals.  Hence ends up with very non-sparse representations.
%
% Signal:        Hydrogen
% Signal Length: 256
% Dictionary:    WP with D, qmf = (Symmlet 8)
% Observations:
%    (a) The signal: Hydrogen.
%    (b) The ideal phase plane indicates one signal atoms.
%    (c) The brightened phase plane from the frame solution is spread
%        out, exaggerating greatly the intrinsic complexity of the
%        signal. The phase plane shown is brightened; the actual phase 
%        plane is much weaker.
%
% Use:
%    bpfig23                   uses the current solver.
%    ATOMIZER_ENGINE = 1998;   selects original solver.
%    ATOMIZER_ENGINE = 2001;   selects later    solver.
%    Default is most recent solver.
%----------------------------------------------------------------------

%----------------------------------------------------------------------
%        1998: (S. Chen) Original script for BP paper.
% 09 Apr 2001: (M. Saunders) Choice of solvers implemented.
%----------------------------------------------------------------------

help bpfig23
bpengine;

n    = 256;   D = log2(n);   t = (1:n)' / n;
qmf  = MakeONFilter('Symmlet', 8);
x    = InputSignal('Hydrogen', n);
cMOF = MOF(x, 'WP', D, qmf, 0);

switch ATOMIZER_ENGINE
  case 1998
    cBP = BP_Interior (x, 'WP', D, qmf, 0);
  case 2001
    cBP = BP_Interior2(x, 'WP', D, qmf, 0);
end

%--------------------
% Plots
%--------------------
fprintf('\nConstructing figure(%1g) ...\n', FIGURE)
figure(FIGURE);   clf reset;

colormap(1-gray);
subplot(2,2,1);   plot(t, x);
                  title('(a)  Signal: Hydrogen');
subplot(2,2,3);   PhasePlane(cBP, 'WP', n, qmf);
                  title('(b) Ideal Phase Plane')
subplot(2,2,4);   PhasePlane(cMOF, 'WP', n, qmf);
                  title('(c) Phase Plane by MOF')
brighten(FIGURE,-.5);
brighten(FIGURE,-.5);