bobdn.m

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

function [xrec, coef] = BOBDN(y, NameOfDict, par1, par2, par3)
% BOBDN -- Best Orthogonal Basis DeNoising Method for WP/CP
%  Usage
%	[xrec, coef] = BOBDeNoise(y, NameOfDict, par1, par2, par3)
%  Inputs
%       y               1-d signal to be de-noised; column vector
%       NameOfDict      name of the dictionary, 'WP' or 'CP'
%       par1		degree of finest frequency partition, in case of 'WP'
%			depth of finest time splitting, in case of 'CP'
% 	par2		orthonormal quadrature mirror filter, in case of 'WP'
%			0, in case of 'CP'
%	par3		0
%  Outputs
%       xrec            cleaned 1-d signal; column vector
%       coef            coef of cleaned signal; column vector
%  Description
%	Assumes noise level == 1
%  References
%	D.L. Donoho and I.M. Johnstone
%	``Ideal de-noising in an orthonormal basis chosen from a
%	libaray of bases''
%	Comptes Rendus de l'Academie des Sciences, To apper.
%  See also
%	WPDeNose, CPDeNoise
%

cmmdstr = [NameOfDict 'DeNoise(y, par1, par2)'];
[clean bb st] = eval(cmmdstr);
cmmdstr = [NameOfDict 'Analysis(clean, par1, par2)'];
pkt = eval(cmmdstr);
coef = UnpackBasisCoeff(bb, pkt);
pkt = zeros(size(pkt));
temp = PackBasisCoeff(bb, pkt, coef);
coef = temp(:);
xrec = clean';