File: ifft2-1.m

package info (click to toggle)
octave2.1 1%3A2.1.73-19
  • links: PTS, VCS
  • area: main
  • in suites: lenny
  • size: 37,108 kB
  • ctags: 20,884
  • sloc: cpp: 106,508; fortran: 46,978; ansic: 5,720; sh: 4,991; makefile: 3,230; yacc: 3,132; lex: 2,892; lisp: 1,715; perl: 778; awk: 174; exp: 134
file content (22 lines) | stat: -rw-r--r-- 349 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
## ifft2-1.m
##
## Author: David Billinghurst (David.Billinghurst@riotinto.com.au)
##         Comalco Research and Technology
##         02 May 2000
M=12;
N=7;

m=3;
n=2;

x = 2*pi*(0:1:M-1)/M;
y = 2*pi*(0:1:N-1)/N;

sx = cos(m*x);
sy = cos(n*y);

S = kron(fft(sx)',fft(sy));
answer=kron(sx',sy);
s = ifft2(S);

all( all( abs(s-answer) < 30*eps ) )