I have a system with measured input (u
) and output (y
).
I assume that this is an linear time-invariant (LTI) system and I want to find its impulse response (ir
).
This is the Matlab code I use to find ir
:
ir = real(ifft(fft(y)./fft(u)));
And this is how I reconstruct the system output:
Y = conv(ir, u);
However, as can be seen below, the reconstructed output differs greatly from the true output. What am I missing here?
The full matlab code including data can be found here