fft8_s1.m

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

function y=fft8_s1(x)
phi_fft(1:4)=0;w=j*2*pi/4;


for k1=0:1
    for k0=0:1
        index=2*k1+k0+1;
        phi_fft(index)=floor(exp((-1)*w*k1*k0)*16384*sqrt(2));
      
    end
end

y=x;

for k1=0:1
    k0=1:4;index=4*k1+k0;
    x(index)=floor(x(index).*phi_fft(k0)/32768);
end

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



y=x;