Skip to main content

Browse the archive

https://github.com/cran/pracma
21 November 2023, 19:34:12 UTC
Revision b5e4bf28fcba9f5eaffbeecfb0bc307452d074ee authored by Hans W. Borchers on 01 November 2014, 00:00:00 UTC, committed by Gabor Csardi on 01 November 2014, 00:00:00 UTC 
version 1.7.7
1 parent d57c14d
Raw File
Tip revision: b5e4bf28fcba9f5eaffbeecfb0bc307452d074ee authored by Hans W. Borchers on 01 November 2014, 00:00:00 UTC
version 1.7.7
Tip revision: b5e4bf2
fft.R 
##
##  f f t . R  Fourier Transform
##


ifft <- function(x) {
    if (length(x) == 0)
        return(c())
    if ( (!is.vector(x, mode="numeric") && !is.vector(x, mode="complex")))
        stop("Argument 'x' must be real or complex vector.")
    
    fft(x, inverse = TRUE) / length(x)
}


fftshift <- function(x) {
    stopifnot(is.double(x) || is.complex(x) || is.integer(x))
    if (!is.vector(x))
        stop("Argument 'x' must be a real or complex vector.")

    m <- length(x)
    p <- ceiling(m/2)
    idx <- c((p+1):m, 1:p)
    x[idx]
}


ifftshift <- function(x) {
    stopifnot(is.double(x) || is.complex(x) || is.integer(x))
    if (!is.vector(x))
        stop("Argument 'x' must be a real or complex vector.")

    m <- length(x)
    p <- floor(m/2)
    idx <- c((p+1):m, 1:p)
    x[idx]
}
back to top