paradigma 0.3.2__py3-none-any.whl → 0.4.1__py3-none-any.whl

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

@@ -1,40 +0,0 @@
-%-------------------------------------------------------------------------------
-% 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

@@ -1,34 +0,0 @@
-%-------------------------------------------------------------------------------
-% 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

@@ -1,29 +0,0 @@
-%-------------------------------------------------------------------------------
-% 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

@@ -1,76 +0,0 @@
-%-------------------------------------------------------------------------------
-% 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

@@ -1,176 +0,0 @@
-%--------------------------------------------------------------------------------
-% 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

@@ -1,11 +0,0 @@
-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

@@ -1,97 +0,0 @@
-%--------------------------------------------------------------------------------
-%
-% 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

@@ -1,149 +0,0 @@
-%-------------------------------------------------------------------------------
-% 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

@@ -1,15 +0,0 @@
-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

@@ -1,13 +0,0 @@
-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