Frequency scale in fft after zero filling (odd/even number of data)

15 views (last 30 days)
FW
FW on 21 Apr 2019
Commented: Star Strider on 21 Apr 2019
Product: Matlab 2017b
If we have data set "A" which consists of time versus light absorbance.There is another set of data "B". Both data are sampled at the same rate Fs.
The two sets are to be divided in the frequency domain. C= fft(A)/fft(B). Prior to division, in order to make the length of A and B equal, zero filling in done or one can use the fft(X, n). Sometimes there are odd number of data points and sometimes, even in A and B after zero filling.
In such a case, how do we make a frequency axis for C in MATLAB? The number of original data points have changed after zero filling in either A or B.
Example: The key question is, assume we had initially 1000 data points in A, and 800 data points in B. Let us use C= fft (A, 1024) ./ fft(B, 1024). We get a new vector C. How should we make a frequency scale for C now given that the number of data points were changed.
Thanks.

Accepted Answer

Star Strider
Star Strider on 21 Apr 2019
I am not certain what you are doing. One way to make the frequencies equal would be to use the second argument, ‘n’, in the fft function, to be the same for all of them (and preferably relatively large, so that all your data have the same frequency vector). This has the effect of ‘zero-padding’ the argument vector, increasing the frequency resolution and equalising the frequencies of the various signals.
The only way to determine if this approach is appropriate to your data is to try it with them.
  6 Comments
Star Strider
Star Strider on 21 Apr 2019
As always, my pleasure.
... whenever we do zero padding of the original data by using fft( X, n) option, our frequency axis should be made as per our chosen "n" rather than original length of t.
Yes.
The ‘Fv’ vector should be:
Fv = linspace(0, 1, fix(numel(FTy1)/2)+1)*Fn;
My apologies for the error. (Since ‘FTy1’ and ‘FTy2’ are the same lengths here, there is no need to compute the ‘Fv’ vector separately for each one.) This also assumes the two are vectors. If not, use the length function instead of numel to calculate the correct length of ‘Fv’.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by