92 lines
2.5 KiB
Matlab
92 lines
2.5 KiB
Matlab
|
function indx = residual_resampling(weights)
|
||
|
|
||
|
%@info:
|
||
|
%! @deftypefn {Function File} {@var{indx} =} resample (@var{weights},@var{noise})
|
||
|
%! @anchor{particle/traditional_resampling}
|
||
|
%! @sp 1
|
||
|
%! Resamples particles (Resampling à la Kitagawa or stratified resampling).
|
||
|
%!
|
||
|
%! @sp 2
|
||
|
%! @strong{Inputs}
|
||
|
%! @sp 1
|
||
|
%! @table @ @var
|
||
|
%! @item weights
|
||
|
%! n*1 vector of doubles, particles' weights.
|
||
|
%! @item noise
|
||
|
%! n*1 vector of doubles sampled from a [0,1] uniform distribution (stratified resampling) or scalar double
|
||
|
%! sampled from a [0,1] uniform distribution (Kitagawa resampling).
|
||
|
%! @end table
|
||
|
%! @sp 2
|
||
|
%! @strong{Outputs}
|
||
|
%! @sp 1
|
||
|
%! @table @ @var
|
||
|
%! @item indx
|
||
|
%! n*1 vector of intergers, indices.
|
||
|
%! @end table
|
||
|
%! @sp 2
|
||
|
%! @strong{This function is called by:}
|
||
|
%! @ref{particle/resample}
|
||
|
%! @sp 2
|
||
|
%! @strong{This function calls:}
|
||
|
%!
|
||
|
%! @end deftypefn
|
||
|
%@eod:
|
||
|
|
||
|
% Copyright (C) 2011 Dynare Team
|
||
|
%
|
||
|
% This file is part of Dynare.
|
||
|
%
|
||
|
% Dynare is free software: you can redistribute it and/or modify
|
||
|
% it under the terms of the GNU General Public License as published by
|
||
|
% the Free Software Foundation, either version 3 of the License, or
|
||
|
% (at your option) any later version.
|
||
|
%
|
||
|
% Dynare is distributed in the hope that it will be useful,
|
||
|
% but WITHOUT ANY WARRANTY; without even the implied warranty of
|
||
|
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
||
|
% GNU General Public License for more details.
|
||
|
%
|
||
|
% You should have received a copy of the GNU General Public License
|
||
|
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
||
|
|
||
|
% AUTHOR(S) frederic DOT karame AT univ DASH evry DOT fr
|
||
|
% stephane DOT adjemian AT univ DASH lemans DOT fr
|
||
|
|
||
|
% What is the number of particles?
|
||
|
number_of_particles = length(weights);
|
||
|
|
||
|
% Set vectors of indices.
|
||
|
jndx = 1:number_of_particles;
|
||
|
indx = zeros(1,number_of_particles);
|
||
|
|
||
|
% Multiply the weights by the number of particles.
|
||
|
WEIGHTS = number_of_particles*weights;
|
||
|
|
||
|
% Compute the integer part of the normalized weights.
|
||
|
iWEIGHTS = fix(WEIGHTS);
|
||
|
|
||
|
% Compute the number of resample
|
||
|
number_of_trials = number_of_particles-sum(iWEIGHTS);
|
||
|
|
||
|
if number_of_trials
|
||
|
WEIGHTS = (WEIGHTS-iWEIGHTS)/number_of_trials;
|
||
|
EmpiricalCDF = cumsum(WEIGHTS);
|
||
|
u = fliplr(cumprod(rand(1,number_of_trials).^(1./(number_of_trials:-1:1))));
|
||
|
j=1;
|
||
|
for i=1:number_of_trials
|
||
|
while (u(i)>EmpiricalCDF(j))
|
||
|
j=j+1;
|
||
|
end
|
||
|
iWEIGHTS(j)=iWEIGHTS(j)+1;
|
||
|
end
|
||
|
end
|
||
|
|
||
|
k=1;
|
||
|
for i=1:number_of_particles
|
||
|
if (iWEIGHTS(i)>0)
|
||
|
for j=k:k+iWEIGHTS(i)-1
|
||
|
indx(j) = jndx(i);
|
||
|
end
|
||
|
end
|
||
|
k = k + iWEIGHTS(i);
|
||
|
end
|