fastmdcsynthesis.m

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

function x = FastMDCSynthesis(c, D, fineness, par3)
% FastMDCSynthesis -- Synthesis operator for Multi-Duration-Cosine dictionary
%  Usage 
%    	x = FastMDCSynthesis(c, D, fineness)
%  Inputs:
%	x	the signal, column vector
%	D	the depth, 0<=D<=log2(n)
%   fineness	how many times finer in frequency than the standard DCT
%  Output:
%	c	the coefs, column vector
%  See Also:
%	FastMDCAnalysis, FastDCTAnalysis, FastDCTSynthesis
%

m = length(c);
n = m / (fineness * (D+1));
t = (0:(n-1))' + .5;
x = zeros(n,1);
for d=0:D,
	s = 2^d;
	cc = c((1:(fineness*n)) + d*n*fineness);
	x = x + exp(-t/s) .* FastDCTSynthesis(cc,fineness) ...
		* (2*n/s*(1-exp(-2*n/s)))^.5;
end