我已經開始（一個小項目）來計算頻域圖像的功率譜。 所以，我有什麼到現在是這樣的： %// close all; clear all; %// not generally appreciated
img = imread('ajw_pic.jpg','jpg'); % it is a color image
img = rgb2gray(img); %// change to gray
psd
根據卷積定理，時域卷積是fft域中的乘積。有了正確的補零，它的工作原理： % convolution in time domain
a = [1 2 3];
b = [4 5 6];
c = conv(a,b);
a_padded=[a 0 0]; b_padded=[b 0 0];
c_bis=ifft(fft(a_padded).*fft(b_padded));
% we do
我具有周期性術語 v(x) = sum over K of [exp(iKx) V(K) ]
其中K = 2 * PI * N/A，其中A是術語，且n = 0,1,2的週期性， 3 .... 現在我想找到傅立葉係數V（K）對應於特定K.假設我有（x）的具有10000點v的矢量線爲 x = 0,0.01a,0.02a,...a,1.01a,....2a....100a
使得我的格子的大小是1