testfrft.asv

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

N=128;
T=1;
t=linspace(0,T,N);
y =cos(100*t);
dt=t(2)-t(1);
f=1/dt;
clear p_saved hn_saved E_saved
a=3;
%for a=0:0.05:2
    fy=Disfrft(y,a);
    %fys=cdpei(y,a);
    fyss=frft(y,a);
    fysa=fft(y);
    fyre=frft(fyss,-3);
    fyif=ifft(y);
    % blue,green,red,cyan
    figure(1)
    subplot(611);plot(t,real([fyss]));
    title(['a = ',num2str(a)]);
    subplot(612);plot(t,real([fysa]));
    subplot(613);plot(t,real([y]));
    subplot(614);plot(t,real([fyre]));
    subplot(615);plot(t,real([fyif]));
    %subplot(312);plot(x,imag([fyss]));
    %subplot(313);plot(x,abs([fyss]));
    %pause(0.7);
  %end