signal processing - Frequency response using FFT in MATLAB - Stack Overflow

Looks OK at first glance, but a couple of things you're missing: you should apply a window function to the time domain data before the FFT, see e.g. http://en.wikipedia.org/wiki/Window_function for windowing in general and http://en.wikipedia.org/wiki/Han

stackoverflow.com