I would use this syntax (from the documentation):
- data = iddata(y,u,Ts) creates an iddata object containing a time-domain output signal y and input signal u, respectively. Ts specifies the sample time of the experimental data.
Specifically with respect to your data, I would do something like this:
time = data(:,1);
Ts = mean(diff(time));
time_data = iddata(data(:,2),data(:,3),Ts);
freq_data = fft(time_data);
freq_units = freq_data.units;
Fr = freq_data.frequency;
sys = tfest(freq_data,2)
NOTE — I do not have your data so this is UNTESTED CODE. I cannot be certain it will work without modification.
Best Answer