paradigma 0.1.5__py3-none-any.whl → 0.2.0__py3-none-any.whl

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/fold_vector_full.m

@@ -0,0 +1,40 @@
+%-------------------------------------------------------------------------------
+% fold_vector_full: fold the vector for frequency decimation:
+%
+%    y[n] = y[n] + sum_{p=0}^{a-l} x[pJ+n]
+%
+%    for n=0,1,...,J-1
+%
+%    which returns DFT{y[n]} = Y[ak]
+%
+% Syntax: y=fold_vector_full(x,J,a)
+%
+% Inputs: 
+%     x,J,a - 
+%
+% Outputs: 
+%     y - 
+%
+% Example:
+%     
+%
+
+% John M. O' Toole, University College Cork
+% Started: 23-04-2014
+%
+% last update: Time-stamp: <2014-04-23 17:11:54 (otoolej)>
+%-------------------------------------------------------------------------------
+function y=fold_vector_full(x,J,a)
+% Assume that all parameters (J,a) are sane.
+
+y=zeros(J,1);
+x=x(:);
+
+n=0:J-1;
+y(n+1)=x(n+1);    
+for p=1:a-1
+  y(n+1)=y(n+1)+x(p*J+n+1);    
+end
+
+% scale by decimation factor:
+y=y./a;

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/fold_vector_half.m

@@ -0,0 +1,34 @@
+%-------------------------------------------------------------------------------
+% fold_vector_half: Fold vector
+%
+%    y[n] = y[n] + sum_{p=0}^{a-l} x[pJ+n]
+%
+%    for n=0,1,...,Jh
+%
+% Syntax: y=fold_vector_half(x,J,Jh,a)
+%
+% Inputs: 
+%     x,J,Jh,a - 
+%
+% Outputs: 
+%     y - 
+%
+% Example:
+%     
+%
+
+% John M. O' Toole, University College Cork
+% Started: 23-04-2014
+%
+% last update: Time-stamp: <2014-04-23 14:05:04 (otoolej)>
+%-------------------------------------------------------------------------------
+function y=fold_vector_half(x,J,Jh,a)
+% Assume that all parameters (J,Jh,a) are sane.
+
+y=zeros(Jh+1,1);
+x=x(:);
+
+n=0:(Jh);
+for p=0:a-1
+  y(n+1)=y(n+1)+x(p*J+n+1);    
+end

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/gen_LFM.m

@@ -0,0 +1,29 @@
+%-------------------------------------------------------------------------------
+% Generate a linear frequency modulated (LFM) signal
+%
+%  USE: sig=gen_LFM(N,fstart,fstop,phase)
+%
+%  INPUT:
+%       N      = length of LFM signal
+%       fstart = starting frequency
+%       fend   = ending frequency
+%       phase  = (optional) phase offset
+%
+%  OUTPUT:
+%       sig = LFM signal of length-N
+%
+%  EXAMPLE:
+%       N=256; fstart=0.1; fend=0.4;
+%       x=gen_LFM(N,fstart,fend);
+%       plot(x); 
+%-------------------------------------------------------------------------------
+function sig=gen_LFM(N,fstart,fstop,phase)
+if(nargin<1 || isempty(N)) N=128; end
+if(nargin<2 || isempty(fstart)) fstart=0.1; end
+if(nargin<3 || isempty(fstop)) fstop=0.4; end
+if(nargin<4 || isempty(phase)) phase=0; end
+
+
+
+n=0:N-1;
+sig=cos( 2*pi.*(fstart.*n + ((fstop-fstart)/(2*N)).*(n.^2)) + phase );

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/get_analytic_signal.m

@@ -0,0 +1,76 @@
+%-------------------------------------------------------------------------------
+% get_analytic_signal: analytic signal as per [1]
+%
+% Syntax: [z,N2,N,Nh]=get_analytic_signal(z)
+%
+% INPUT:
+%      s = real-valued signal
+%
+% OUTPUT:
+%      z = complex-valued analytic signal
+%
+% EXAMPLE:
+%      N=64;  f=0.1; N2=2*N;
+%      s=cos( 2*pi*f.*(0:N-1));     
+%      z=get_analytic_signal(s);
+%      plot(1:N2,real(z),1:N2,imag(z));
+%      legend('real','imaginary');
+%     
+%
+% [1] J. M. O' Toole, M. Mesbah, and B. Boashash, "A New Discrete Analytic Signal for
+%     Reducing Aliasing in the Discrete Wigner-Ville Distribution", IEEE Trans.  on Signal
+%     Processing, vol. 56, no. 11, pp. 5427-5434, Nov. 2008.
+
+
+% John M. O' Toole, University College Cork
+% Started: 14-04-2014
+%
+% last update: Time-stamp: <2016-08-11 17:18:41 (otoolej)>
+%-------------------------------------------------------------------------------
+function [z,N2,N,Nh]=get_analytic_signal(z)
+if(rem(length(z),2))
+    warn_str=sprintf(['odd-length signal.\nCurrently, code works for even-length signal ' ...
+                      'only.\n\nRemoving last sample of the signal to force to ' ...
+                     'even-length.\n']);
+    warning(warn_str);
+    z=z(1:(end-1));
+end
+
+
+if(isreal(z)) 
+  z=gen_analytic(z); 
+else
+    warn_str=sprintf(['using complex-valued signal; assuming this is an analytic signal ' ...
+                      'zero-padded to length-2N?\n\n If unsure, input the real part ' ...
+                     'of signal only, i.e. real(x)\n'] );
+    warning(warn_str);
+end
+N2=length(z); N=N2/2; Nh=ceil(N/2);
+
+if(N~=fix(N))
+  error('Analytic signal must be length 2N.');
+end
+
+
+
+function z=gen_analytic(s1)
+%---------------------------------------------------------------------
+% generate the analytic signal (with zero-padding in the time-direction)
+%---------------------------------------------------------------------
+s1=real(s1);
+N=length(s1); N2=2*N;
+
+% 1. zero-pad N-point signal to 2N
+s1=[s1(:); zeros(N,1)];
+S1=fft(s1);
+
+% 2. Get analytic signal of 2N-point signal
+H=zeros(N2,1); 
+H([1 N+1])=1;
+H(2:N)=2;
+z_cb=ifft(S1.*H);
+
+
+% 3. Force the second half of the time-domain signal to zero.
+z=[z_cb(1:N); zeros(N,1)];
+

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/get_window.m

@@ -0,0 +1,176 @@
+%--------------------------------------------------------------------------------
+% General function to calculate a window function (mostly a wrapper for the existing
+% window functions).
+% 
+% USE: win = get_window( win_length, win_type, [win_param], [DFT_WINDOW], [Npad] )
+%
+% INPUT:
+%       win_length = length of window
+%       win_type   = type of window: { 'delta' | 'rect' | 'bart' | 'hamm'
+%                    | 'hann' | 'tukey' | 'gauss' | 'cosh' | 'blackmanharris' }
+%       win_param  = (optional) window parameter.
+%       DFT_WINDOW = (optional) parameter { 0 | 1 }. 
+%                     If 1 returns DFT of window. 
+%       Npad       = (optional) zero-pad window to length Npad.
+%
+% OUTPUT:
+%       win = window 
+%
+% See also: BARTLETT, HANNING, HAMMING, GAUSSWIN, TUKEYWIN, CIRCSHIFT, FFT
+%
+% EXAMPLE
+%      N=64; 
+%      win=get_window(N,'tukey',0.4,0,2*N); 
+%      plot(win);
+
+
+% Copyright (C) 2007,2008 John M. O' Toole, The University of Queensland
+% 
+% This program 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.
+% 
+% This program 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 this program.  If not, see <http://www.gnu.org/licenses/>.
+%--------------------------------------------------------------------------------
+function win = get_window( win_length, win_type, win_param, DFT_WINDOW, Npad )
+if( nargin<3 ) win_param=[]; end
+if( nargin<4 ) DFT_WINDOW=0; end
+if( nargin<5 ) Npad=0; end
+
+
+win=get_win(win_length,win_type,win_param,DFT_WINDOW);
+win=shift_win(win);
+
+if(Npad>0)
+  win=pad_win(win,Npad);
+end
+
+
+
+function win=get_win(win_length, win_type, win_param, DFT_WINDOW )
+%--------------------------------------------------------------------------------
+% Get the window. Negative indices are first.
+%--------------------------------------------------------------------------------
+
+switch win_type
+ case 'delta'
+  win=zeros(win_length,1);
+  wh = floor(win_length/2);
+  win(wh+1)=1;
+ case 'rect'
+  win(1:win_length) = 1;
+ case {'bart', 'bartlett'}
+  win = bartlett( win_length );
+ case {'hamm', 'hamming'}
+  win = hamming( win_length );
+ case {'hann', 'hanning'}
+  win = hanning( win_length );
+ case 'tukey'
+  % NOTE: seems to be problem with Octave's (v2.9.12) tukeywin.m for N odd.
+  if(isempty(win_param))
+    win = tukeywin( win_length );
+  else
+    win = tukeywin( win_length, win_param );
+  end
+ case 'gauss'
+  if(isempty(win_param))
+    win = gausswin( win_length );
+  else
+    win = gausswin( win_length, win_param );
+  end
+ case 'cosh'
+  win_hlf = fix( win_length / 2);
+
+  if(isempty(win_param))
+      win_param = 0.01;
+  end
+  for m = -win_hlf:win_hlf
+      win(mod(m,win_length)+1) = cosh( m ).^( -2 * win_param );
+  end
+  win = fftshift(win);
+  case {'blackmanharris', 'bmh'}  
+    win = blackmanharris(win_length, 'symmetric');
+  case {'nuttall'}  
+    win = nuttallwin(win_length, 'symmetric');
+  case {'dolph', 'chebwin'}
+    win = chebwin(win_length, win_param);
+  
+  otherwise
+  error(['Unknown window type ' win_type]);
+end
+  
+
+%---------------------------------------------------------------------
+% If want the DFT of win
+%---------------------------------------------------------------------
+if(DFT_WINDOW)
+  win=circshift(win(:),ceil(win_length/2));
+  win=fft(win);
+  win=circshift(win(:),floor(win_length/2));
+end
+
+
+function w=shift_win(w)  
+%--------------------------------------------------------------------------------
+% Shift the window so that positive indices are first.
+%--------------------------------------------------------------------------------
+N=length(w);
+w=circshift(w(:),ceil(N/2));
+
+
+
+function w_pad=pad_win(w,Npad)
+%--------------------------------------------------------------------------------
+%
+% Pad window to Npad. 
+%
+% Presume that positive window indices are first. 
+%
+% When N is even use method described in [1]
+%
+%   References:
+%     [1] S. Lawrence Marple, Jr., Computing the discrete-time analytic
+%     signal via FFT, IEEE Transactions on Signal Processing, Vol. 47,
+%     No. 9, September 1999, pp.2600--2603.
+%
+%--------------------------------------------------------------------------------
+w=w(:);
+w_pad=zeros(Npad,1);
+N=length(w);
+Nh=floor(N/2);
+if(Npad<N) error('Npad is less than N'); end
+
+
+% Trival case:
+if(N==Npad) 
+  w_pad=w; 
+  return; 
+end
+
+
+ % For N odd:
+if( rem(N,2)==1 )
+  n=0:Nh;
+  w_pad(n+1)=w(n+1); 
+  n=1:Nh;
+  w_pad(Npad-n+1)=w(N-n+1);
+ 
+  % For N even:
+  % split the Nyquist frequency in two and distribute over positive
+  % and negative indices.
+else
+  n=0:(Nh-1);
+  w_pad(n+1)=w(n+1); 
+  w_pad(Nh+1)=w(Nh+1)/2; 
+
+  n=1:Nh-1;
+  w_pad(Npad-n+1)=w(N-n+1);
+  w_pad(Npad-Nh+1)=w(Nh+1)/2; 
+end

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/isreal_fn.m

@@ -0,0 +1,11 @@
+function flag = isreal_fn(fn,E)
+% ISREAL_FN - is 'fn' real valued
+
+if(nargin<2)
+  E=sum(abs(fn(:)).^2) * 1e-12;
+end
+flag=1;
+
+if( max( abs(imag(fn(:))) ) > E ) flag=0; end
+
+

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/padWin.m

@@ -0,0 +1,97 @@
+%--------------------------------------------------------------------------------
+%
+% Pad window to Npad. 
+%
+% Presume that positive window indices are first. 
+%
+%  USE: w_pad=padWin(w,Npad)
+%
+%  INPUT: 
+%        w    = window (vector) of length N
+%        Npad = pad window to length N (Npad>N)
+%
+%  OUTPUT: 
+%        w_pad = window of length N zeropadded to length Npad
+%
+%
+% When N is even use method described in [1]
+%
+%   References:
+%     [1] S. Lawrence Marple, Jr., Computing the discrete-time analytic
+%     signal via FFT, IEEE Transactions on Signal Processing, Vol. 47,
+%     No. 9, September 1999, pp.2600--2603.
+%
+%
+% STARTED: 04-01-2008
+
+%   Copyright (c) 2010, John M. O' Toole, The University of Queensland
+%   All rights reserved.
+%
+%  Redistribution and use in source and binary forms, with or without
+%  modification, are permitted provided that the following
+%  conditions are met:
+%      * Redistributions of source code must retain the above
+%        copyright notice, this list of conditions and the following
+%        disclaimer.
+%      * Redistributions in binary form must reproduce the above
+%        copyright notice, this list of conditions and the following
+%        disclaimer in the documentation and/or other materials
+%        provided with the distribution.
+%      * Neither the name of the The University of Queensland nor the 
+%        names of its contributors may be used to endorse or promote 
+%        products derived from this software without specific prior 
+%        written permission.
+%  
+%  THIS SOFTWARE IS PROVIDED BY JOHN M. O' TOOLE ''AS IS'' AND ANY
+%  EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+%  IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+%  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL JOHN M. O' TOOLE BE
+%  LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+%  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
+%  OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR
+%  BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
+%  LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+%  (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE
+%  USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
+%  DAMAGE.
+%--------------------------------------------------------------------------------
+function w_pad=padWin(w,Npad)
+DB=0;
+w=w(:);
+w_pad=zeros(Npad,1);
+N=length(w);
+Nh=floor(N/2);
+if(DB) dispVars(N,Npad); end
+if(Npad<N) error('Npad is less than N'); end
+
+
+% Trival case:
+if(N==Npad) 
+  w_pad=w; 
+  return; 
+end
+
+
+ % For N odd:
+if( rem(N,2)==1 )
+  if(DB) dispDB('odd'); end
+ n=0:Nh;
+ w_pad(n+1)=w(n+1); 
+ n=1:Nh;
+ w_pad(Npad-n+1)=w(N-n+1);
+ 
+ % For N even:
+ % split the Nyquist frequency in two and distribute over positive
+ % and negative indices.
+else
+  if(DB) dispDB('even'); end
+ n=0:(Nh-1);
+ w_pad(n+1)=w(n+1); 
+ w_pad(Nh+1)=w(Nh+1)/2; 
+
+ n=1:Nh-1;
+ w_pad(Npad-n+1)=w(N-n+1);
+ w_pad(Npad-Nh+1)=w(Nh+1)/2; 
+end
+
+

paradigma/ppg/hr_functions/TFD toolbox JOT/utils/vtfd.m

@@ -0,0 +1,149 @@
+%-------------------------------------------------------------------------------
+% Plot TFD and time/frequency plots
+%
+% USE:  vtfd(tfd,sig,Fs)
+%
+%  INPUT:
+%        tfd    - the time-frequency distribution 
+%        sig    - time-domain signal
+%        Fs     - sampling frequency
+
+%
+%   Copyright (c) 2010, John M. O' Toole, The University of Queensland
+%   All rights reserved.
+%
+%  Redistribution and use in source and binary forms, with or without
+%  modification, are permitted provided that the following
+%  conditions are met:
+%      * Redistributions of source code must retain the above
+%        copyright notice, this list of conditions and the following
+%        disclaimer.
+%      * Redistributions in binary form must reproduce the above
+%        copyright notice, this list of conditions and the following
+%        disclaimer in the documentation and/or other materials
+%        provided with the distribution.
+%      * Neither the name of the The University of Queensland nor the 
+%        names of its contributors may be used to endorse or promote 
+%        products derived from this software without specific prior 
+%        written permission.
+%  
+%  THIS SOFTWARE IS PROVIDED BY JOHN M. O' TOOLE ''AS IS'' AND ANY
+%  EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+%  IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+%  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL JOHN M. O' TOOLE BE
+%  LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+%  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
+%  OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR
+%  BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
+%  LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+%  (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE
+%  USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
+%  DAMAGE.
+%-------------------------------------------------------------------------------
+function [h_image,h_time,h_freq]=vtfd(tfd,s1,FS,n,k)
+if(nargin<3 || isempty(FS)) 
+  FS=1; 
+end
+if(nargin<2 || isempty(s1)) 
+  TFD_ONLY=1; 
+  Ntime=size(tfd,1);  
+else 
+  TFD_ONLY=0; 
+  Ntime=length(s1);
+end
+if(nargin<4 || isempty(n)), n=[]; end
+if(nargin<5 || isempty(k)), k=[]; end
+
+
+[N,M]=size(tfd);
+
+
+% $$$ clf;
+
+
+%---------------------------------------------------------------------
+% Frame sizes:
+%---------------------------------------------------------------------
+X_AXIS_WIDTH=0.75;
+Y_AXIS_HEIGHT=0.75;
+X_AXIS_START=0.23;
+Y_AXIS_START=0.23;
+TIME_FREQ_PLOTS_WIDTH=0.15;
+TIME_FREQ_PLOTS_GAP=0.07;
+
+
+%---------------------------------------------------------------------
+% X/Y tick labels
+%---------------------------------------------------------------------
+if(isempty(n))
+    ntime=1:Ntime; ntime=ntime./FS;
+    n=linspace(ntime(1),ntime(end),N);
+else
+    ntime=n;
+end
+Mh_time=ceil(Ntime/2);  Mh=ceil(M/2);
+if(isempty(k))
+
+    k=linspace(0,0.5,Mh);
+    k_time=linspace(0,0.5,Mh_time);
+    k=k.*FS;
+else
+    k_time=k;
+end
+% $$$ keyboard;
+
+if(TFD_ONLY)
+    imagesc(k,n,tfd); axis('xy');
+    
+    % seems to bug with Octave and gnuplot, so have to force square:
+    if(is_octave), axis('square');  end
+  
+else
+  %---------------------------------------------------------------------
+  % 1. time plot
+  %---------------------------------------------------------------------
+  s1=real(s1); 
+  h_time=subplot(2,2,1);
+  set(h_time,'position',[(X_AXIS_START-TIME_FREQ_PLOTS_WIDTH-TIME_FREQ_PLOTS_GAP) ...
+                      Y_AXIS_START TIME_FREQ_PLOTS_WIDTH ...
+                      Y_AXIS_HEIGHT]);
+  plot(s1,1:length(s1)); 
+  axis('tight');
+  grid('on');
+  set(h_time,'xticklabel',[]); set(h_time,'yticklabel',[]);
+
+
+  %---------------------------------------------------------------------
+  % 2. freq plot
+  %---------------------------------------------------------------------
+  S1=fft(s1);
+  h_freq=subplot(2,2,4);
+  set(h_freq,'position',[X_AXIS_START (Y_AXIS_START- ...
+                                       TIME_FREQ_PLOTS_WIDTH-TIME_FREQ_PLOTS_GAP) ...
+                      X_AXIS_WIDTH TIME_FREQ_PLOTS_WIDTH]);
+  plot(1:Mh_time,abs(S1(1:Mh_time)).^2);
+  axis('tight');
+  grid('on');
+  set(h_freq,'xticklabel',[]); set(h_freq,'yticklabel',[]);
+
+
+  %---------------------------------------------------------------------
+  % 2. time-freq plot
+  %---------------------------------------------------------------------
+  h_image=subplot(2,2,2);
+  set(h_image,'position',[X_AXIS_START Y_AXIS_START X_AXIS_WIDTH ...
+                      Y_AXIS_HEIGHT]);
+
+  imagesc(k,n,tfd); axis('xy');
+  
+  % seems to bug with Octave and gnuplot, so have to force square:
+  if(is_octave), axis('square');  end
+end
+
+
+function r = is_octave ()
+  persistent x;
+  if (isempty (x))
+    x = exist ('OCTAVE_VERSION', 'builtin');
+  end
+  r = x;

paradigma/ppg/preprocessing/preprocessing_imu.m

@@ -0,0 +1,15 @@
+function [IMU_filt, t_tar] = preprocessing_imu(t_orig, IMU_sample, fs_imu) 
+%%----Preprocessing of the IMU pipeline----%%                        
+    % Resampling to a uniform sampling rate at 30 Hz                                  
+    t_tar = t_orig(1):1/fs_imu:t_orig(end);         % Correct target for resampling (by usage of t_start = t_orig(1) and t_end = t_orig(end));
+    signal_resampled = spline(t_orig, IMU_sample', t_tar);
+    
+    % High-pass filter for detrending
+    [b,a]=butter(4,0.2/(fs_imu/2),'high');  % high-pass filter for gravity removal
+    IMU_filt = filtfilt(b,a,double(signal_resampled'));
+   
+   % [zhi,phi,khi] = butter(4,0.2,'high');
+   % soshi = zp2sos(zhi,phi,khi);
+   % IMU_filt2 = filtfilt(soshi,double(signal_resampled'));
+end
+

paradigma/ppg/preprocessing/preprocessing_ppg.m

@@ -0,0 +1,13 @@
+function [PPG_data, t_tar] = preprocessing_ppg(t_orig, PPG_sample, fs_ppg)         
+              
+    % Resampling to a uniform sampling rate at 30 Hz                  
+    t_tar = t_orig(1):1/fs_ppg:t_orig(end);         % Correct target for resampling (by usage of t_start and t_end);
+    signal_resampled = spline(t_orig, PPG_sample, t_tar); % Spline interpolation used for resampling
+
+    % Band-pass filter for detrending
+
+    [b,a]=butter(4,[0.4, 3.5]/(fs_ppg/2),'bandpass');  % band-pass filter for detrending
+    PPG_filt = filtfilt(b,a,double(signal_resampled)); % Apply filter
+    PPG_data = PPG_filt';   % In case of a sample package error where PPG_filt consists of multiple "parts"
+
+end