frft_test.asv

分数阶Fourier变换程序

%% FrFT test program.
clc;
close all;
clear all;
k = 100;
fc = 30;
fs = 10000;
ts = 1/fs;
n = 0:ts:3;
s = exp(j*(2*pi*fc*n + pi*k*n.^2));
real_s = real(s);
imag_s = imag(s);
%%%%% plot %%%%%%%
figure;
subplot(2,1,1);
plot(n,real_s);
subplot(2,1,2);
plot(n,imag_s);

F_s = fft(s);
figure;
plot(abs(F_s));

F