_{1}

Many attempts have been made to derive a simple expression for the characteristic function of the lognormal distribution. This would be useful for computing the sum of lognormal variables, either with each other, or with other statistical variables. In this paper , we provide a simple formula for the characteristic function, which is exact, closed and computable. An extension to the sum of <i>correlated</i> lognormals, used in the pricing of Asian options, is a consequence of this approach.

The usefulness of finding a simple expression for the characteristic function (CF, a.k.a Fourier transform FT) of the lognormal distribution has been made clear in a variety of contexts: from insurance and finance to engineering to biology. The difficulty is that the CF has an essential singularity at 0, and therefore a branch cut must be made [

Thus, we have the conclusions in [

“A closed form expression of the Laplace transform of the lognormal distribution does not exist.”

and further in [

“Recently, Asmussen et al. (2016a) reported a closed form approximation for the Laplace transform … that works quite well for small values of the σ parameter and is asymptotically equivalent to the approximated Laplace transform.”

and from [

“The usual method of evaluating the distribution function requires the computation of many convolutions of the conditional distribution of the amount of a claim given that a claim has occurred. When the expected number of claims is large, the computation can become unwieldy even with modern large scale electronic computers.”

The last quotation relates to convolutions, which are obviated if the CF of the claims amount distribution can be found.

In the context of Asian options, we find a similar view from [

“Contrary to the case of geometric averaging, the distribution for the arithmetic average is not known. This implies that no exact closed form solution is available for European-style Asian options based on the arithmetic average. Because no analytical formula is available for European-style Asian options, it does not come as a surprise that an analytical solution to the similar American-style problem is not feasible either.”

Nonetheless, when viewed from complex function theory, the lognormal is a well behaved function; in fact it is entire (i.e. integral, with no finite singularities) in the complex plane. [

In this paper, we provide and prove an exact expression for the lognormal CF, which can be easily verified numerically by the inverse Fourier transform. It is not the purpose of this paper to assess how it relates to the approximations, asymptotic or otherwise, in the literature. In fact, that is a task for the authors of those approximations.

For a variance parameter σ 2 and mean μ the pdf of the lognormal distribution may be written as

p ( x ) = 1 2 π σ 2 x e − ( ln x + μ ) 2 2 σ 2

(where μ = σ 2 2 , for x to have unit mean, which is assumed without loss of generality). Its FT is using x = e ζ

p ^ ( z ) = 1 2 π σ 2 ∫ ζ ∈ R e i z e ζ e − ( ζ + μ ) 2 2 σ 2 d ζ (FT)

which is a difficult integral to evaluate numerically because of the presence of the oscillatory term e i z e ζ , which grows exponentially for real z.

Remark 1 There is considerable ambiguity as to the meaning of “closed form” for a mathematical expression. In its strictest sense, it admits only arithmetic expressions. In its most general sense, it admits the operations of integration and differentiation, as many functions (e.g. the gamma function) need to be defined in terms of these operations. However the term precludes algorithms and numerical solutions of a system of equations.

A simple approach, introduced by [

· the horizontal segments x = ± X , and

· the vertical segments y = 0 and y = π 2 .

Alternatively we can use the contour in the lower half plane bounded by

· the horizontal segments x = ± X , and

· the vertical segments y = 0 and − π 2

This is depicted in

This is justified since the lognormal is an entire function of ζ , provided we show that the integrals on the vertical segments tend to zero as | X | → ∞ . We distinguish between contours in the upper and lower half planes, as is seen below.

For the upper half plane, the integral along the vertical segment for 0 ≤ y ≤ π 2 , given X and z is:

1 2 π σ 2 ∫ 0 ≤ y ≤ π 2 e i z e X + i y ︸ oscillatory e − ( X + i y + μ ) 2 2 σ 2 ︸ dampening d y (1)

Now the dampening term in 1 has the real part:

ℜ ( X + μ + i y ) 2 = 1 2 [ ( X + μ + i y ) 2 + ( X + μ − i y ) 2 ] = ( X + μ ) 2 − y 2 .

Thus the dampening term has no definitive influence on the convergence of the vertical integral as X → ± ∞ . The oscillatory term can be dealt with as follows.

The integral is over a finite range for y in the upper half plane and tends to 0 if ℜ ( i z e i y ) < 0 , as with z = a + i b :

ℜ ( i z e i y ) = ℜ ( a + i b ) ( i cos y − sin y ) = − a sin y − b cos y .

If 0 ≤ y ≤ π 2 , so that s i n y , c o s y ≥ 0 , then ℜ ( i z e i y ) ≤ 0 for a = ℜ ( z ) ≥ 0 , b = ℑ ( z ) ≥ 0 .

The integral is over a finite range for y in the lower half plane and again tends to 0 if ℜ ( i z e i y ) < 0 . For − π 2 ≤ y ≤ 0 , s i n y ≤ 0 and c o s y ≥ 0 , this applies to ℜ ( z ) ≤ 0 , ℑ ( z ) ≥ 0 .

Thus for ℑ ( z ) ≥ 0 we have in the limit as | X | → ∞ , we derive an alternative expression for the lognormal CF:

p ^ ( z ) = 1 2 π σ 2 { ∫ ζ ∈ R e − z e ζ e − ( ζ + μ + π 2 i ) 2 2 σ 2 d ζ if ℜ ( z ) ≥ 0 ∫ ζ ∈ R e z e ζ e − ( ζ + μ − π 2 i ) 2 2 σ 2 d ζ if ℜ ( z ) ≤ 0 (2)

But the above expressions for p ^ ( z ) may be analytically continued to ℑ ( z ) < 0 within the indicated ranges of ℜ ( z ) , as the integrals converge, regardless of the sign of ℑ ( z ) ^{1}. It is also evident that p ^ ( z ) then obeys the general equality for FTs of a real distribution:

p ^ ( z ) ¯ = p ^ ( − z ¯ ) .

For convenience, denote z ˜ = − sign [ ℜ ( z ) ] z , so that 2 may be written succinctly as:

p ^ ( z ) = 1 2 π σ 2 ∫ ζ ∈ R e z ˜ e ζ e − ( ζ + μ + π 2 i ) 2 2 σ 2 d ζ (3)

and so

p ( x ) = 1 ( 2 π ) 3 / 2 σ ∬ e − i z x e z ˜ e ζ e − ( ζ + μ + π 2 i ) 2 2 σ 2 d ζ d z

Remark 2 It has been suggested (anonymously) that Equation (2) is equivalent to the formula in [

p ^ ( z ) = 1 2 π i ∫ k − i ∞ k + i ∞ Γ ( s ) e − μ + ln ( z ) − sign ( z ) i π 2 + σ 2 s 2 2 d s for z ∈ R \ { 0 } (4)

where k is any number k > 0 and Γ ( s ) is presumably the gamma function. This alleged equivalence is clearly false:

· The term for Γ ( s ) does not appear in 2;

· The above expression is an integral over a vertical strip, unlike the horizontal strip in 2;

· Moreover the author does not exploit or attempt to compute the expression in 4; instead he resorts to asymptotic approximations.

Remark 3 The approach of [

g ( x ) = λ e − λ p ( x ) + λ x ∫ 0 x y p ( y ) g ( x − y ) d y (5)

where g ( x ) is the density of total claims x over a unit period and λ a Poisson parameter. The equation was developed to avoid convolutions in computing g ( x ) = e − λ ∑ n = 0 λ n n ! p ∗ n . In fact the FT has a simple form g ^ ( z ) = e λ ( p ^ − 1 ) , if we can derive the FT p ^ ( z ) , and use the inverse FT. However, given the form of p ( x ) , 5 has a very complicated kernel. It is claimed (anonymously) that this is “state of the art” for insurance purposes. However:

· the equation is an inhomogeneous Volterra integral equation of the 2nd kind, the most intractable possible. Since it applies to any claim amount density p ( x ) , it cannot be used to solve for p ^ ( x ) unless other conditions are imposed.

· Generally it needs to be solved by discretization of the variable x and numerical quadrature of a system of equations, which renders the approach only slightly less problematic than the convolution approach;

· No solution, nor an approach to one, is offered, approximate or otherwise;

· It does not and cannot deal in any way with the CF of the lognormal, the theme of this paper, which is useful and important in fields other than insurance, and for a multitude of different purposes (e.g. electrokinetics).

There is clearly a branch cut along the imaginary axis (in fact, two branch cuts originating from 0 along the imaginary axis). By taking the inverse transform 1 2 π ∫ e − i z x p ^ ( z ) d z , we should recover p ( x ) . This is confirmed in

We calculate the FT of the lognormal average, where each u i is lognormally distributed, such as:

v = u 1 + u 2 + ⋯ + u n n = U n . (6)

Taking averages allows us to compare the effect of convolution for various levels of n. The lognormal sum U has the CF p ^ ( z ) n and thus the average lognormal v has CF p ^ ( z / n ) n . This simplifies greatly the procedures discussed in [

The FTs of v ( x ) and that of U ( x ) may be written and are related by

v ( x ) = 1 2 π ∫ − ∞ ∞ e − i x z v ^ ( z ) d z (7)

= 1 2 π ∫ − ∞ ∞ e − i x z p ^ ( z / n ) n d z = n 2 π ∫ − ∞ ∞ e − i n x z p ^ ( z ) n d z (8)

= n ( 2 π ) n / 2 + 1 σ n { ∫ 0 ∞ [ e − i x z ∫ e − z e ζ e − ( ζ + μ + π 2 i ) 2 2 σ 2 d ζ ] n d z if ℜ ( z ) ≥ 0 ∫ − ∞ 0 [ e − i x z ∫ e z e ζ e − ( ζ + μ − π 2 i ) 2 2 σ 2 d ζ ] n d z if ℜ ( z ) ≤ 0 (9)

so that

v ( x ) = n 2 π ∫ − ∞ ∞ e − i n x z p ^ ( z ) n d z (10)

= 2 n ( 2 π ) n / 2 + 1 σ n ℜ ∫ 0 ∞ [ e − i x z ∫ e − z e ζ e − ( ζ + μ + π 2 i ) 2 2 σ 2 d ζ ] n d z (11)

so that

v ( x ) = n 2 π ∫ − ∞ ∞ e − i n x z p ^ ( z ) n d z (12)

= 2 n ( 2 π ) n / 2 + 1 σ n ℜ ∫ 0 ∞ [ e − i x z ∫ e − z e ζ e − ( ζ + μ + π 2 i ) 2 2 σ 2 d ζ ] n d z (13)

The charts of v for various values of n are shown in

Remark 4 The integrations used for

Remark 5 It is easy to see that the FT of a mixture of lognormals and other distributions, such as the normal, may be derived in the same way. Clearly the lognormals do not need to be identically distributed. Where the lognormals are correlated, it is possible to extend the one dimensional results to n dimensions, as is shown in the next section.

Remark 6 Where the sum for U in 6 has negative components, it is straightforward to allow for them. In this event U can assume negative values, so that the expression in 10 does not apply. To deal with this, consider the term

p ^ ( − z ) = 1 2 π σ 2 ∫ ζ ∈ R e − i z e ζ e − ( ζ + μ ) 2 2 σ 2 d ζ .

Then 10 still applies, with the roles of ℜ ( z ) ≥ 0 and ℜ ( z ) ≥ 0 reversed. This may be regarded as a special case of the following section. This, and the next section, generalize the results in [

Remark 7 [

p ^ ( z ) = ∫ e − z 2 π x e − i z x ( 2 σ ) 2 e − 1 2 x 2 d x .

· [

· The integral has a quadratic exponent, and is thus solvable analytically. This contradicts all the literature.

For a vector x , the multivariate lognormal is the exponential of a multivariate normal of x , and thus has the form

p ( x ) = e − 1 2 ( l n x ) T Σ − 1 ( l n x ) ( 2 π ) n | Σ | ∏ x 1 x 2 ⋯ x n

where Σ is the covariance matrix of the elements of x = [ x 1 , x 2 , ⋯ , x n ] , Since Σ − 1 is positive definite, it can be written as

Σ − 1 = A T A

for some symmetric matrix A , so that

p ( x ) = e − 1 2 ‖ A ( l n x ) ‖ 2 ( 2 π ) n | Σ | ∏ x 1 x 2 ⋯ x n

and

p ^ ( z ) = | Σ | ( 2 π ) n ∫ 1 ∏ x 1 x 2 ⋯ x n e i z ⋅ x e 1 2 ‖ A ( l n x ) ‖ 2 d x

and letting x = e ζ we get the simplification for z ∈ C n :

p ^ ( z ) = | A | ( 2 π ) n ∫ e i z ⋅ e ζ e − 1 2 ‖ A ζ ‖ 2 d ζ .

here we have the same problem as in FT: where any component of z = [ z 1 , z 2 , ⋯ , z n ] has a positive real part, there is an uncomfortable oscillation. For convenience, denote z ˜ = − sign [ ℜ ( z ) ] z for the multidimensional case, as suggested by 3.

Fortunately we can deal with this problem by applying contour integration, component by component, according to the sign of ℜ ( z ) . There are 2^{n} possible signs, so that there are a similar number of branches around z = 0 , and we have the following 3:

p ^ ( z ) = | A | ( 2 π ) n ∫ e − i z ˜ ⋅ e ζ e − 1 2 ‖ A ζ ‖ 2 d ζ .

A further simplification can be made by writing r = A ζ so that

p ^ ( z ) = 1 ( 2 π ) n ∫ e i z ˜ ⋅ e A − 1 r e − 1 2 ‖ r ‖ 2 d r .

Letting A − 1 = [ a 1 , a 2 , ⋯ , a n ] , this becomes

p ^ ( z ) = 1 ( 2 π ) n ∫ e i z ˜ ⋅ e A − 1 r e − 1 2 ‖ r ‖ 2 d r = 1 ( 2 π ) n ∫ e i z ˜ ⋅ e a 1 r 1 e i z ˜ ⋅ e a 2 r 2 ⋯ e i z ˜ ⋅ e a n r n e − 1 2 ‖ r ‖ 2 d r 1 d r 2 ⋯ d r n = 1 ( 2 π ) n ∏ j = 1 n [ ∫ e i z ˜ ⋅ e a j r e − 1 2 r 2 d r ] . (14)

since 14 is a separable function of r . For computation this requires n integrations for each value of z .

Once p ^ ( z ) is calculated, the inverse FT may be used to compute p ( x ) in the usual way. From there, the distribution of any function of x , including a weighted sum of the form y = b T x may be found. In particular, this may be found as follows.

If q ( y ) denotes the density of the weighted sum of y, (which must be positive) then Q ( y ) = ∫ y ∞ q ( s ) d s is its cdf, whose extremes are 0 and 1. Thus:

Q ( y ) = ∫ 0 ≤ y ≤ b T x p ( x ) d x

and since q ( y ) = − ∂ Q ∂ y :

q ^ ( z ) = − ∫ e i y z ∂ Q ∂ y d y = − e i y z Q | + i z ∫ e i y z Q d y = 1 + i z Q ^ ( z )

q ^ ( z ) = 1 ( 2 π ) n ∂ ∂ z ∏ j = 1 n [ ∫ e i z ˜ b ⋅ e a j r e − 1 2 r 2 d r ] .

and the density of y is q = − ∂ Q ∂ y , and its FT is

q ^ = − ∫ e i y z ∂ Q ∂ y d y = − e i y z Q | 0 ∞ + ∫ i z e i y z Q d y = 1 + ∫ i z e i y z Q d y = 1 − ∫ [ e i ( b T x ) z − 1 ] p ( x ) d x = ∫ e i ( b T x ) z p ( x ) d x = p ( z b ) = 1 ( 2 π ) n ∫ e i z b ˜ ⋅ e A − 1 r e − 1 2 ‖ r ‖ 2 d r .

Taking the inverse FT, we get finally:

p ^ ( z ) = 1 ( 2 π ) n ∫ e i z ˜ ⋅ e A − 1 r e − 1 2 ‖ r ‖ 2 d r .

p ( y ) = 1 ( 2 π ) n ∫ e i y z e i z b ˜ ⋅ e A − 1 r e − 1 2 ‖ r ‖ 2 d r d z . (15)

and

Q ^ ( z ) = ∫ 0 ≤ y ≤ b T x e i y z p ( x ) d x d y = 1 i z ∫ [ e i ( b T x ) z − 1 ] p ( x ) d x = 1 ( 2 π ) n 1 i z ∫ [ e i ( b T x ) z ] p ( x ) d x − 1 i z

and so using the inverse FT:

Q ( y ) = 1 ( 2 π ) n + 2 ∫ e − i y z i z [ e i ( b T x ) z ] p ( x ) d x d z − ∫ e − i y z i z d z = 1 ( 2 π ) n + 2 [ ∫ 1 + e i ( b T x − y ) z i z p ( x ) d x d z ] (16)

since the last improper integral becomes:

− ∫ e − i s z i z d z = − 2 π [ ∫ 0 ∞ e − i s z − e + i s z ] 1 z d z = 2 π ∫ 0 ∞ sin s z z d z = 1

so finally

Q ( y ) = 1 ( 2 π ) n [ ∫ e i ( b T x ) z y p ( x ) d x − 1 2 π ] .

Asian options may be priced according to the distribution of the multivariate lognormal. The Asian call option is payable as the excess (over a given strike K) of the (arithmetic) average of prices at predetermined times T = [ t 1 , t 2 , ⋯ , t n ] . Assuming that prices P t follow a standard unit Wiener process W t (and assuming for simplicity that the stock has zero drift, no dividends and unit variance),

we have d P t P t = d W t , and thus

P t = e W t − 1 2 t .

which is lognormal. Thus the covariance of the exponent at two times, s and t, is m i n ( s , t ) , and the covariance matrix is given by:

Σ i j = min ( t i , t j ) .

(This does not depend on any drift that might be introduced into the stock price process). The distribution of the average price at the relevant times may then be determined from the above procedure. This contrasts with the approximations based on two moments [

We take T = [ 1,2,3,4 ] / 4 over a unit period, for which

Σ = [ 0.25 0.25 0.25 0.25 0.25 0.50 0.50 0.50 0.25 0.50 0.75 0.75 0.25 0.50 0.75 1.00 ]

This contradicts the assertion in [

A = [ + 0.11401 − 0.17467 + 0.15360 − 0.06066 − 0.21426 + 0.07441 + 0.18842 − 0.13985 + 0.28868 + 0.28868 − 0.00000 − 0.28868 + 0.32827 + 0.61694 + 0.83121 + 0.94521 ]

A − 1 = [ + 1.61073 − 2.01175 + 1.15470 + 0.15838 − 2.46778 + 0.69867 + 1.15470 + 0.29765 + 2.17013 + 1.76910 − 0.00000 + 0.40102 − 0.85705 − 1.31308 − 1.15470 + 0.45603 ]

Since the price P t = e W t − 1 2 t has a weighting of e − 1 2 t when applied to the lognormal, we take b = e − 1 2 T / n . The distribution and cumulative distribution of the average stock price Q = ∑ i = 1 n P t i / n are depicted in

With a strike of K = 0.4 , the Asian or European options can be priced as

C = ∫ x ≥ K ( x − K ) p ( x ) d x

where p ( x ) is the density for either the Asian price 14, or it can be taken as the European price at expiration:

p ( x ) = e − ( x + 1 ) 2 2 2 π

In numerical terms these are easily computed:

Asian = 0.8912

European = 0.4485

Evidently the Asian is the more valuable as the average price over a period of times has a greater likelihood of exceeding the strike than at any particular given time.

Contrary to much of what is asserted in the literature, the lognormal distribution can be comfortably accommodated within complex function theory, as given in 3. No approximations are required, except to the extent that numerical integrations over infinite domains may not be exactly precise. This facilitates the application of the lognormal to all its areas of relevance, including the pricing of Asian options.

The author declares no conflicts of interest regarding the publication of this paper.

Leung, A.P. (2020) The Lognormal Characteristic Function in Several Dimensions, with Application to Asian Options. Journal of Mathematical Finance, 10, 399-411. https://doi.org/10.4236/jmf.2020.103024