function [signal,check] = filtersignal(signal,fs,type,order,hpf,lpf,stopf)
% INPUT: signal - Signal
% fs - Sampling frequency
% type - This is a vector of ones and zeros.
% order - This is a cell with the high pass and low pass filter orders
%
% hpf - cutoff frequency for highpass filter
% lpf - cutoff frequency for lowpass filter
% stopf - cutoff frequency for band stop filter
% other - 1 or 0, define a new filter when other ==1
% OUTPUT: signal - Filtered signal
%
%
% Author(s): Jonathan Moeyersons (Jonathan.Moeyersons@esat.kuleuven.be)
% Sabine Van Huffel (Sabine.Vanhuffel@esat.kuleuven.be)
% Carolina Varon (Carolina.Varon@esat.kuleuven.be)
%
% Version History:
% - 06/05/2019 JM Initial version
%
% Copyright (c) 2019, Jonathan Moeyersons, KULeuven-ESAT-STADIUS
%
% This software is made available for non commercial research purposes only
% under the GNU General Public License. However, notwithstanding any
% provision of the GNU General Public License, this software may not be
% used for commercial purposes without explicit written permission after
% contacting jonathan.moeyersons@esat.kuleuven.be
%
% 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 <https://www.gnu.org/licenses/>.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Pre-define check
check = 1;
% Set the signal in the correct form (collumn vectors)
[R_temp,K_temp] = size(signal);
if R_temp < K_temp
signal = signal';
end
% Get the Nyquist frequency
Nf = fs/2;
% High pass filter
if type(1)
try
Cf = hpf/Nf; % Cutoff freq. Passband corner freq. 0.5Hz
[bh,ah] = butter(order{1},Cf,'high'); % High pass filter
signal = filtfilt(bh,ah,signal);
catch ME
% Show an error dialog
errordlg(ME.message)
% Change check
check=0;
end
end
% Low pass filter
if type(2)
try
Cf = lpf/Nf; % Cutoff freq. Passband corner freq. 100Hz
[bl,al] = butter(order{2},Cf,'low'); % Low pass filter
signal = filtfilt(bl,al,signal);
catch ME
% Show an error dialog
errordlg(ME.message)
% Change check
check=0;
end
end
% Notch filter
if type(3)
try
% Define zeros
omega = stopf/Nf*pi;
zero(1) = cos(omega)+1i*sin(omega);
zero(2) = cos(omega)-1i*sin(omega);
% Define poles
pole = 0.99*zero;
% Define coefficients
bs = [1 -zero(1)-zero(2) zero(1)*zero(2)] ;
as = [1 -pole(1)-pole(2) pole(1)*pole(2)] ;
signal = filtfilt(bs,as,signal);
catch ME
% Show an error dialog
errordlg(ME.message)
% Change check
check = 0;
end
end
% Self designed filter
if type(4)
try
% Design filter
Filt = designfilt;
if ~isempty(Filt)
if strcmp(Filt.FrequencyResponse,'differentiator') || strcmp(Filt.FrequencyResponse,'hilbert')
% Filter the signal
signal = filter(Filt,signal);
else
% Filter the signal
signal = filtfilt(Filt,signal);
end
else
% Change check
check = 0;
end
catch ME
% Show an error dialog
errordlg(ME.message)
% Change check
check = 0;
end
end