fileID = fopen('D:\\Andy\\Logs\\MATLAB\\td_20.txt','r');
formatSpec = '%f';
td = fscanf(fileID,formatSpec);
td = transpose(td)
[meas,timept]=size(td); %%% epo_tem = signal --> 21*15000
disp('meas');
disp(meas);
disp('timept');
disp(timept);
	
%disp('td')
%for i=1:length(td)
%    x = sprintf("%f", td(i));
%    disp(x)
%end

sr = 300
fn=[0.5 6]
[a,b]=butter(6,fn(1)/(sr/2),'high')

[a,b]=butter(6,fn(1)/(sr/2),'high');  % create parameters of butterworth filter
x = sprintf('fn(1)=%d', fn(1));
disp(x);
disp('a')
disp(a)
disp('b')
disp(b)
		
for i=1:meas
    Y(i,:)=filtfilt(a,b,td(i,:));  % zero-phase shift filtering
end
		
disp('td Y')
for i=1:length(td)
    x = sprintf("%f %f", td(i), Y(i));
    disp(x)
end