## Diff of /gabor/zak.m[5f391e] .. [84b1c9]  Maximize  Restore

### Switch to unified view

a/gabor/zak.m b/gabor/zak.m
1
function c=zak(f,a);

1
function c=zak(f,a);

2
%ZAK  Zak transform

2
%ZAK  Zak transform

3
%   Usage:  c=zak(f,a);

3
%   Usage:  c=zak(f,a);

4
%

4
%

5
%   ZAK(f,a) computes the Zak transform of f with parameter _a.

5
%   zak(f,a) computes the Zak transform of *f* with parameter *a*.  The

6
%   The coefficients are arranged in an _a x _L/a matrix, where _L is the

6
%   coefficients are arranged in an $a \times L/a$ matrix, where *L* is the

7
%   length of f.

7
%   length of *f*.

8
%

8
%

9
%   If f is a matrix, then the transformation is applied to each column.

9
%   If *f* is a matrix then the transformation is applied to each column.

10
%   This is then indexed by the third dimension of the output.

10
%   This is then indexed by the third dimension of the output.

11
%

11
%

12
%   Assume that $c=zak(f,a)$, where *f* is a column vector of length *L* and

13
%   $N=L/a$. Then the following holds for $m=0,\ldots,a-1$ and $n=0,\ldots,N-1$

12
%

14
%

13
%   Assume that c=ZAK(f,a), where f is a column vector of length L and

15
%   ..                     N-1

14
%   N=L/a. Then the following holds for $m=0,...,a-1$ and $n=0,...,N-1$

16
%     c(m+1,n+1)=1/sqrt(N)*sum f(m-k*a+1)*exp(2*pi*i*n*k/N)

17
%                          k=0

15
%

18
%

16
%M                         N-1

17
%M    c(m+1,n+1)=1/sqrt(N)*sum f(m-k*a+1)*exp(2*pi*i*n*k/N)

18
%M                         k=0

19
%F  \begin{eqnarray*}

20
%F  c(m+1,n+1) & = & \frac{1}{\sqrt{N}}\sum_{k=0}^{N-1}f(m-ka+1)e^{2\pi ink/M}

19
%   .. math:: c(m+1,n+1) = \frac{1}{\sqrt{N}}\sum_{k=0}^{N-1}f(m-ka+1)e^{2\pi ink/M}

21
%F  \end{eqnarray*}

22
%

20
%

23
%   See also:  izak

21
%   See also:  izak

24
%

22
%

25
%   References: ja94-4 bohl97-1

23
%   References: ja94-4 bohl97-1

26
24
...
...
60
  % of the Zak transform.

58
  % of the Zak transform.

61
  c(:,:,ii)=dft(reshape(f(:,ii),a,N),[],2);

59
  c(:,:,ii)=dft(reshape(f(:,ii),a,N),[],2);

62
end;

60
end;

63
61
64
62
65
%OLDFORMAT