Mercurial > octave-nkf

view test/octave.test/signal/ifft2-1.m @ 4960:ce01dbd7e026 ss-2-1-58

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2004-09-02 03:47:49 by jwe]
author jwe
date Thu, 02 Sep 2004 03:47:49 +0000
parents 53acb394c728
children
line wrap: on
line source

## 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 ) )