function x = lpr(y, pfilt)% LPR   Multi-level Laplacian pyramid reconstruction%%	x = lpr(y, pfilt)%% Input:%   y:      output of a Laplacian pyramid in a cell vector%   pfilt:  pyramid filter name (see PFILTERS)%% Output:%   x:      reconstructed signal%% See also: LPD% Get the pyramidal filters from the filter name[h, g] = pfilters(pfilt);% Decide extension modeswitch pfilt    case {'9-7', '9/7', '5-3', '5/3', 'Burt'}        extmod = 'sym';            otherwise        extmod = 'per';        endx = y{1};for n = 2:length(y)    x = lprec1(x, y{n}, h, g, extmod);end