thfig24.m

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

% thfig24: BP Thesis Figure 2.4 -- 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)
% Observation:
%	(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 plase plane shown is brightened; the actual phase 
%	    plane is much weaker.
%

help thfig24

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

figure(1); clf
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')