1function c=ref_dfft(f,order)
2%-*- texinfo -*-
3%@deftypefn {Function} ref_dfft
4%@verbatim
5%REF_DFFT  Reference Discrete Fractional Fourier Transform
6%   Usage:  c=ref_dfft(f,order);
7%
8%@end verbatim
9%@strong{Url}: @url{http://ltfat.github.io/doc/reference/ref_dfft.html}
10%@end deftypefn
11
12% Copyright (C) 2005-2016 Peter L. Soendergaard <peter@sonderport.dk>.
13% This file is part of LTFAT version 2.3.1
14%
15% This program is free software: you can redistribute it and/or modify
16% it under the terms of the GNU General Public License as published by
17% the Free Software Foundation, either version 3 of the License, or
18% (at your option) any later version.
19%
20% This program is distributed in the hope that it will be useful,
21% but WITHOUT ANY WARRANTY; without even the implied warranty of
22% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
23% GNU General Public License for more details.
24%
25% You should have received a copy of the GNU General Public License
26% along with this program.  If not, see <http://www.gnu.org/licenses/>.
27
28L=size(f,1);
29
30% Create matrix representation of the DFT
31F=idft(eye(L));
32
33c=(F^order)*f;
34
35
36
37