fft8_s2.m

来自「快速离散傅立叶变换」· M 代码 · 共 29 行

function y=fft8_s2(x)
phi_fft(1:8)=0;w=2*pi/8;
for k2=0:1
    for k1=0:1
        for k0=0:1
            index=4*k2+2*k1+k0+1;
            phi_fft(index)=floor(exp((-1)*j*w*k2*(2*k1+k0))*16384*sqrt(2));
            
        end
    end
end
y=x;


k0=1:8;index=k0;
x(index)=floor(x(index).*phi_fft(k0)/32768);
   

for k2=0:1
    for k0=0:1
        k1=[0 1];index=4*k2+2*k1+k0+1;
        x(index)=wfft2(x(index));
    end
end

 

y=x;