Choosing FFT points

5 次查看（过去 30 天）
joeDiHare
joeDiHare 2011-10-28
Hi All. It is not clear to me how to choose points in a FFT, and matlab Help, doesn't help on this. E.g., I have an L long array and I want an fft window of 32ms (which is 704 samples with a fs of 22kHz); fft(X(1:704)) gives me 704 frequency bins by default, but I want only a subset of frequency "lines", and specifically I need M=128 frequency bins. I tought I would just use the syntax fft(x(1:704),128), however:
x1=ifft( fft( x(1:704),128 ), 704)
is not equal to
x1=ifft( fft( x(1:704) ) )
Why?? (thanks)
  1 个评论
Wayne King
Wayne King 2011-10-28
By the way:
x1=ifft( fft( x(1:704),128 ), 704)
does not take the DFT of x(1:704) as you think, it takes the DFT of x(1:128)
isequal(fft(x(1:128)),fft(x(1:704),128))

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Wayne King
Wayne King 2011-10-28
The bins in the DFT are Fs/N where N is the length of the input vector and Fs is the sampling frequency.
If you only want to evaluate the DFT at a subset of 128 DFT bins, then you can select those from the output of fft()
Also, see the help for goertzel(), which evaluates the DFT at specific bins.
  1 个评论
joeDiHare
joeDiHare 2011-10-29
Thanks for your answer.
Goertzel algorithm looks very good, but there is no IDFT.
Is it correct to use ifft(X), where X is the output of the goertzel "zeropadded" in the frequency bins discarded by the DFT, so that I have back a time domain signal 704 samples long?.
something like this:
ind_freq=1:2:128*2;
G=goertzel(x(1:704),ind_freq);
magn=abs(G); phas=phase(G);
X=zeros(704,2);
X(ind_freq,:)=[magn phas];
x=ifft( X(:,1).*exp(1j*X(:,2) )

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Transforms 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by