View License

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

» Watch video

Highlights from
Fresnel Integral

Join the 15-year community celebration.

Play games and win prizes!

» Learn more

Be the first to rate this file! 8 Downloads (last 30 days) File Size: 2.34 KB File ID: #34134 Version: 1.1

Fresnel Integral

by

Drew Compston (view profile)

 

08 Dec 2011 (Updated )

Compute Fresnel integral.

| Watch this File

File Information
Description

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.

Acknowledgements

Fresnel Integrals inspired this file.

This file inspired Fresnel S And Fresnel C.

MATLAB release MATLAB 7.5 (R2007b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Updates
06 Sep 2012 1.1

Fixed problem for inputs less than -1.6.

Contact us