fourierin_1d {fourierin} | R Documentation |
It computes Fourier integrals of functions of one and two variables on a regular grid.
fourierin_1d(f, lower_int, upper_int, lower_eval = NULL, upper_eval = NULL, const_adj, freq_adj, resolution = NULL, eval_grid = NULL, use_fft = TRUE)
f |
function or a vector of size m. If a function is provided, it must be able to be evaluated at vectors. If a vector of values is provided, such evaluations must have been obtained on a regular grid and the Fourier integral is faster is m is a power of 2. |
lower_int |
Lower integration limit(s). |
upper_int |
Upper integration limit(s). |
lower_eval |
Lower evaluation limit(s). It can be NULL if an evaluation grid is provided. |
upper_eval |
Upper evaluation limit(s). It can be NULL if an evaluation grid is provided. |
const_adj |
Factor related to adjust definition of Fourier transform. It is usually equal to 0, -1 or 1. |
freq_adj |
Constant to adjust the exponent on the definition of the Fourier transform. It is usually equal to 1, -1, 2pi or -2pi. |
resolution |
A vector of integers (faster if powers of two) determining the resolution of the evaluation grid. Not required if f is a vector. |
eval_grid |
Optional matrix with d columns with the points where the Fourier integral will be evaluated. If it is provided, the FFT will not be used. |
use_fft |
Logical value specifying whether the FFT will be used. |
See vignette for more detailed examples.
If w is given, only the values of the Fourier integral are returned, otherwise, a list with the elements
w |
A vector of size m where the integral was computed. |
values |
A complex vector of size m with the values of the integral |
##--- Example 1 --------------------------------------------------- ##--- Recovering std. normal from its characteristic function ----- library(fourierin) #' Function to to be used in integrand myfun <- function(t) exp(-t^2/2) # Compute Foueien integral out <- fourierin_1d(f = myfun, lower_int = -5, upper_int = 5, lower_eval = -3, upper_eval = 3, const_adj = -1, freq_adj = -1, resolution = 64) ## Extract grid and values grid <- out$w values <- Re(out$values) plot(grid, values, type = "l", col = 3) lines(grid, dnorm(grid), col = 4) ##--- Example 2 ----------------------------------------------- ##--- Computing characteristic function of a gamma r. v. ------ library(fourierin) ## Function to to be used in integrand myfun <- function(t) dgamma(t, shape, rate) ## Compute integral shape <- 5 rate <- 3 out <- fourierin_1d(f = myfun, lower_int = 0, upper_int = 6, lower_eval = -4, upper_eval = 4, const_adj = 1, freq_adj = 1, resolution = 64) grid <- out$w # Extract grid re_values <- Re(out$values) # Real values im_values <- Im(out$values) # Imag values # Now compute the real and # imaginary true values of the # characteric function. true_cf <- function(t, shape, rate) (1 - 1i*t/rate)^-shape true_re <- Re(true_cf(grid, shape, rate)) true_im <- Im(true_cf(grid, shape, rate)) # Compare them. We can see a # slight discrepancy on the # tails, but that is fixed # when resulution is # increased. plot(grid, re_values, type = "l", col = 3) lines(grid, true_re, col = 4) # Same here plot(grid, im_values, type = "l", col = 3) lines(grid, true_im, col = 4)