I need to apply the below filter code to the ECG as contained in the file

% DEFINE FILTER PARAMETERS 
Fs  = 2E+6;
Fn  = Fs/2;
Fc1 =  630;
Fc2 =  640;
Fs1 = Fc1*0.8;
Fs2 = Fc2/0.8;
Rp  =    1;
Rs  =   10;
% DESIGN FILTER
[n,Wp]  = ellipord([Fc1 Fc2]/Fn, [Fs1 Fs2]/Fn, Rp, Rs);
[z,p,k] = ellip(n,Rp,Rs,Wp);
[sos,g] = zp2sos(z,p,k);
figure
freqz(sos, 2^14, Fs) 

How to apply the above filter to the file (each row is ecg for 1 person): https://drive.google.com/file/d/1vXb3zY1hlKdXbn6cSvrm7O7OfL1Nx_qe/view?usp=share_link

(I got the code from https://in.mathworks.com/matlabcentral/answers/183972-implementation-of-narrow-band-pass-filter-butterworth#comment_1337312)

Answer 1

[Walter Roberson]

FilteredRowsOfECG = sosfilt(sos, YourRowsOfECG, 2);

[Paul]

I believe that should be:

FilteredRowsOfECG = g*sosfilt(sos, YourRowsOfECG, 2);

[Walter Roberson]

You might well be right, @Paul .

[Neha Sinha]

Thank you