File Exchange

image thumbnail

Fresnel Integral

version 1.1 (2.34 KB) by

Compute Fresnel integral.

5 Downloads

Updated

View License

Usage: F = FCS(X, TERMS) or [C, S] = FCS(X, TERMS)

Returns the Fresnel cosine integral C and sine integral S for each element in X. If only one output is requested, F = C + j*S is output. The convention of this function is to use a pi/2 normalization in the argument of the trignometric function in the definition. To get the Fresnel integral without this normalization, input sqrt(2/pi)*x and multiply the the output by sqrt(pi/2).

For inputs with magnitude less than or equal to 1.6, the Taylor series expansion using TERMS (default 12) number of terms is used. Otherwise, an auxillary function is used. The algorithm is described in the paper:

Klaus D. Mielenz, Computation of Fresnel Integrals. II J. Res. Natl. Inst. Stand. Technol. 105, 589 (2000), pp 589-590 (accessible online: http://nvlpubs.nist.gov/nistpubs/jres/105/4/j54mie2.pdf)

The paper claims the error should be less than 1e-9.

Inputs:
  -X: values to compute Fresnel integral of.
  -TERMS: number of terms to use in Taylor series expansion (optional, default is 12).

Outputs:
  -C: Fresnel cosine integral (with a pi/2 normalization).
  -S: Fresnel sine integral (with a pi/2 normalization).

This file was inspired by the file Fresnel integrals by Peter Volegov. That file uses a MEX file to compute the Fresnel integral using the same method that this file uses. I implemented the same technique directly in MATLAB.

Comments and Ratings (0)

Updates

1.1

Fixed problem for inputs less than -1.6.

MATLAB Release
MATLAB 7.5 (R2007b)
Acknowledgements

Inspired by: Fresnel integrals

Inspired: FresnelS and FresnelC

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video