What should be FFT of a constant function?

Question

Richa Sharma 2012-12-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56277-what-should-be-fft-of-a-constant-function

As per theory,FFT of a constant fn. is a DC value. But when I take an array 'A' of ones of size 172, i.e. A[1 1 1 1 1.....172 times]

FFT(A)

gives DC as well as AC cofts. Why?

2 个评论
显示无隐藏无

Azzi Abdelmalek 2012-12-11

What do you mean by AC coefts?

Image Analyst 2012-12-11

The FFT breaks up the signal into a weighted sum of sinusoidal signals. Any point in an FFT not at the center is the weight (coefficient) of one of the sinusoidal signals that goes into making up the final signal. The center frequency is flat (no sinusoid) and is often called the DC component.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2012-12-11

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56277-what-should-be-fft-of-a-constant-function#answer_68141

Round-off error in the calculations. Look at the magnitudes: everything is down near 10E-15

3 个评论
显示 1更早的评论隐藏 1更早的评论

Muthu Annamalai 2012-12-11

You have the details right, Image Analyst; while the rest of the comments are only partly true.

Image Analyst 2012-12-11

If you're interested in seeing the sinc effect, I posted some nice demo code here: http://www.mathworks.com/matlabcentral/answers/56139#comment_116309

请先登录，再进行评论。

Answer 2

Azzi Abdelmalek 2012-12-11

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56277-what-should-be-fft-of-a-constant-function#answer_68158

编辑：Azzi Abdelmalek 2012-12-11

在 MATLAB Online 中打开

If you mean by AC, sinusoidal signal, when you calculate its FFT, it's important to define the interval. To give a sens to your FFT, you have to calculate it in one period. For the sinusoidal signal, the period is 2*pi, then theoretically, the result will be one value at k=1 and not at k=0 like in a constant signal.

t=0:.1:2*pi-0.1;
g=fft(sin(t));
stem(abs(g))

For a constant the result will be a constant for k=0; and 0 elsewhere

h=ones(1,length(t))
figure;
stem(abs(fft(h)))

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Richa Sharma 2012-12-12

Thank You very much for your response. Actually, my inputs are digital images.I was trying to check the fft response for constant background in images. I am going for 1D FFT in matlab. Say, my image is I and I write J=FFT(I) Then,for constant background in I, first coft. in columns of J should have high coft. value, the rest should be zeros is what I was expecting.But I was getting nonzero cofts for other frequency bins also.I think it is due to Round-off error in the calculations as Roberson said.Kindly correct me if I am wrong. Regards

请先登录，再进行评论。