我們需要改變/重新實現在GSL標準DFT實現,這是 int
FUNCTION(gsl_dft_complex,transform) (const BASE data[],
const size_t stride, const size_t n,
BASE result[],
const gsl_fft_directio
我創建了一個簡單的積分函數和DFT函數,我可以使用我寫的其他代碼。 from math import sin,pi
from time import time
def aintegral(d,step):
return sum(d)*step
def partialdft(d,step,w):
i = 0
x = d
while i/step < l
我試圖實現3D DFT,但我遇到了一些麻煩。我相信我應該做的只是連續3次一維DFT,每個方向一個。假設一維DFT是正確的,你可以看到什麼是錯,此代碼: def dft3d(self, real3d, img3d, nx, ny, nz, dir):
#Transform depth
for i in range(nx):
for j in range(ny):