cycle_conv.m

来自「是自己自创的好程序」· M 代码 · 共 15 行

clear all;
close all;
tic;
a=(randn(1,16)+sqrt(-1)*randn(1,16));
b=[a(13:16),a(1:12)];
% F=ifft(eye(16))*sqrt(16);
% F=ifftshift(F,2);
% a_fft=F'*a.';
% b_fft=F'*b.';
a_fft=fft(a,16)/sqrt(16);
b_fft=fft(b,16)/sqrt(16);
test=(a_fft).*sqrt(-1).*exp(-sqrt(-1)*(2*pi/16)*4*(1:16));
test-b_fft
toc;