test.m

来自「这是学习分数阶傅里叶变换」· M 代码 · 共 24 行

N=128;
T=1;
t=linspace(0,T,N);
%x=12*sin(20*pi*t+pi/4)+5*cos(80*pi*t);
x=cos(100*t);
dt=t(2)-t(1);
x=x(:);
f=1/dt;
X=frft(x,1);
%X=X';
y=ifft(x);
z=fft(x);
p1=real(X);
p2=real(y);
p3=real(z);
F=X(1:N);
Y=y(1:N);
f=f*(0:N-1)/N;
Z=z(1:N);
figure(1)
   subplot(311);plot(f,F);
   subplot(312);plot(f,Y);
   subplot(313);plot(f,Z);
    norm(z-X,2)