hi
you could use
np = 2; % number of poles
nz = 1; % number of zeros
exp_data = readmatrix('data.csv'); %read data
u = exp_data(:,1); % input data
y = exp_data(:,2); % output data
Ts = 0.01; % sampling time in seconds
data = iddata(y,u,Ts);
sys = tfest(data,np,nz)
in MATLAB for estimating transfer function from data
more details can be found here