function rr_int = ApplyResamplingToRR(t_win,rr_0,HRVparams)
%
% rr_int = ApplyResamplingToRR(t_win,rr_0,HRVparams)
%
% OVERVIEW: Resampling of RR intervals before to apply methods that
% require resampled data for frequency domain analysis
%
% INPUT:
% t_win - Vector containing start time of each window
% or ECG/PPG waveform
% rr_0 - vector of rr intervals to be resampled
% HRVparams - struct of settings for hrv_toolbox analysis that can
% be obtained using InitializeHRVparams.m function
% HRVparams = InitializeHRVparams();
% OUTPUT
% rr_int - resampled rr intervals
%
% REPO:
% https://github.com/cliffordlab/PhysioNet-Cardiovascular-Signal-Toolbox
% ORIGINAL SOURCE AND AUTHORS:
% This script written by Giulia Da Poian
% Dependent scripts written by various authors
% (see functions for details)
% COPYRIGHT (C) 2018
% LICENSE:
% This software is offered freely and without warranty under
% the GNU (v3 or later) public license. See license file for
% more information
sf = HRVparams.freq.resampling_freq; % (Hz) resampling frequency
ti = t_win(1):1/sf:t_win(end); % time values for interp.
interp_method = HRVparams.freq.resample_interp_method;
% Chose the resampling method to use
switch interp_method
case 'cub'
rr_int = interp1(t_win,rr_0,ti','spline')'; % cubic spline interpolation
case 'lin'
rr_int = interp1(t_win,rr_0,ti','linear')'; %linear interpolation
otherwise
warning('using cubic spline method')
rr_int = interp1(t_win,rr_0,ti','spline')'; % cubic spline interpolation
end