fwppbssynthesis.m

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

function sig = FWPPBSSynthesis(wp,qmf,dqmf)
% FWPPBSSynthesis -- Synthesize signal from full wavelet packet table 
%                    (periodized, biorthogonal,symmetric)
%  Usage
%    sig = FWPSynthesis(wp,qmf,dqmf)
%  Inputs
%    wp       Wavelet Packet Table 
%    qmf      Wavelet orthonormal qmf
%    dqmf    dual quadrature mirror filter (symmetric, dual of qmf)
%  Outputs
%    sig      the synthesized signal
%  See Also
%    FWPSynthesis
% 
	[n,L] = size(wp);
	D = L-1;
%
	sig = wp(:,L)';
	for d= (D-1):-1:0,
		for b=0:(2^d-1)
			 odd = rem(b,2);
			 rson = 2*b + 1 - odd;
			 lson = 2*b + odd;
			 sig(packet(d,b,n)) = ...
				 UpDyadLo_PBS(sig(packet(d+1,lson,n)),dqmf) + ...
				 UpDyadHi_PBS(sig(packet(d+1,rson,n)),qmf) + ...
				 wp(packet(d,b,n),d+1)';
		end
	end
	sig = sig';