Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Find Transfer Function from Experimental Data

  1. Dec 16, 2012 #1
    I would like to find a useful function to predict the output based on the input. I'm trying to figure out how to get the transfer function given actual input and output data.

    What I've tried is getting a transfer function using fourier transforms. I'm attemping to do a fft on both the input and the output and dividing the output over the input and using fft of the absolute value to get the transfer function. The problem is it doesn't allow me to predict the output very well, though it does give me something that looks like a sum of exponential which is what I expect. Do I need to subtract an offset or apply gain and an offset? (y = mx + b sort of thing). I think I'm missing some key concept that I don't know how to describe.

    Here's a simple sketch of what I'm trying to do.


    Here's some pseudo code of what I've tried. I don't have matlab or anything fancy to do this for me, I'm trying to write it for an arduino project.

    input_fft = fft(input);
    output_fft = fft(output);
    tf_fd = output_fft / input_fft;
    tf = fft(abs(tf_fd));

    Just a point in the right direction would be great! Maybe I'm going about this all wrong.
  2. jcsd
  3. Dec 17, 2012 #2
    You can download opensource Octave. It has system identification functions, arx for example.
  4. Dec 17, 2012 #3

    I like Serena

    User Avatar
    Homework Helper

    Welcome to PF, arcticMouse! :smile:

    Are you aware that you should calculate a newoutput for some newinput with:
    newoutput = fftinverse(fft(newinput) * tf_fd)
    Last edited: Dec 17, 2012
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook