matlab fractional fft

little matlab demo for fractional fft...enjoy

N=64;
t=2*pi/N*[-N/2:N/2];c=1;
for k=-N/2:N/2 B(c,:)=1/sqrt(N+1)*exp(-j*k*t);c=c+1;end;
f=exp(-t.^2);
f=f-min(f);
f=f';
figure(1);clf;set(gcf,'Renderer','zbuffer');
for n=0:0.1:2.0
Bn=B^n;fn=abs(Bn*f);plot(fn,'b');hold on;plot(f,'r--');hold off;
text(10,1.7,sprintf('%0.2f',n));
grid on;axis([1 N+1 0 2]);axis square;
drawnow;
end;

creatively acquired from :

http://www.cs.dartmouth.edu/~farid/r...fracderiv.html