testfrft.asv

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

x=-pi:0.02:pi;
y =cos(x);
clear p_saved hn_saved E_saved
for a=0:0.05:2
    fy=Disfrft(y,a);
    fys=cdpei(y,a);
    fyss=frft(y,a);
    fysa=fft(y);
    % blue,green,red,cyan
    figure(1)
    subplot(511);plot(x,real([fyss]));
    title(['a = ',num2str(a)]);
    subplot(512);plot(x,real([fysa]));
    subplot(514);plot(x,real([y]));
    subplot(512);plot(x,imag([fyss]));
    subplot(513);plot(x,abs([fyss]));
    pause(0.7);
end