fft8_s3.m

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

function y=fft8_s3(x)
y=x;

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