JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 623945 10.1155/2013/623945 623945 Research Article Pricing Currency Option Based on the Extension Principle and Defuzzification via Weighting Parameter Identification Xu Jixiang 1,2 Tan Yanhua 3 Gao Jinggui 4 Feng Enmin 1 Palhares Reinaldo Martinez 1 School of Mathematical Sciences Dalian University of Technology Dalian Liaoning 116024 China dlut.edu.cn 2 School of Sciences Liaoning Shihua University Fushun Liaoning 113001 China lnpu.edu.cn 3 School of Sciences Hebei University of Technology Tianjin 300130 China hebut.edu.cn 4 College of Information Science and Engineering Shandong University of Science and Technology Qingdao Shandong 266510 China sdust.edu.cn 2013 14 03 2013 2013 24 11 2012 20 01 2013 2013 Copyright © 2013 Jixiang Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We present a fuzzy version of the Garman-Kohlhagen (FG-K) formula for pricing European currency option based on the extension principle. In order to keep consistent with the real market, we assume that the interest rate, the spot exchange rate, and the volatility are fuzzy numbers in the FG-K formula. The conditions of a basic proposition about the fuzzy-valued functions of fuzzy subsets are modified. Based on the modified conditions and the extension principle, we prove that the fuzzy price obtained from the FG-K formula for European currency option is a fuzzy number. To simplify the trade, the weighted possibilistic mean (WPM) value with a weighting function is adopted to defuzzify the fuzzy price to a crisp price. The numerical example shows our method makes the α-level set of fuzzy price smaller, which decreases the fuzziness. The example also indicates that the WPM value has different approximation effects to real market price by taking different values of weighting parameter in the weighting function. Inspired by this example, we provide a method, which can identify the optimal parameter.

1. Introduction

With the fast growing of the trading volume in the foreign exchange market, the trading of currency option has increased. It is well know that currency option manages the risk of the foreign exchange market. Hence, an appropriate formula for pricing currency option is becoming extremely significant. Garman and Kohlhagen (1983)  derived the closed-form formula (G-K formula) for pricing European currency option by the method of Black and Scholes (1973) . There are six variables in the G-K formula, that is, the spot exchange rate, the volatility of spot exchange rate, the domestic and foreign risk-free interest rate, the strike price, and the time to maturity date. However, some variables in the G-K formula are assumed as constants, which is inconsistent with the empirical phenomena, such as the volatility smile. According to these empirical phenomena, many researchers have been devoted to modify the G-K model, such as Chesney and Scott (1989) , Amin and Jarrow (1991) , Heston (1993) , Bates (1996) , Sarwar and Krehbiel (2000) , and Carr and Wu (2007) .

In general, the data, for instance, the spot exchange rate, the domestic or foreign risk-free interest rate, cannot be recorded or collected precisely. Usually, these financial tools may have different values in different commercial banks and financial institutions. Meanwhile, there are differences between the buy (bid) and the sale (ask) price for these financial tools, which are similar to bid-ask spreads in dealer markets for the stock. The existence of bid-ask spreads implies that we cannot precisely get the true market price. The bid-ask spreads are considered to be natural bands to represent the uncertainty around the market price, which is similar to fuzzy number. Furthermore, expert opinions or statistical estimators of market parameters may be expressed as fuzzy numbers. Therefore, the fuzzy set theory proposed by Zadeh (1965)  is useful to model such uncertainty. In a parallel development, the fuzzy set theory is applied to option pricing. Yoshida (2003) [10, 11] discussed the valuation of European call and put options in fuzzy environment. Muzzioli and Torricelli (2004)  discussed the pricing of European options in a multiperiod binomial model. Wu (2007)  obtained the fuzzy pattern of Black-Scholes formula using the extension principle. Xu et al. (2009)  and Zhang et al. (2012)  investigated a jump-diffusion option pricing model under fuzzy environment. Liu (2009)  applied fuzzy approach in  to the G-K formula of currency option. Xu et al. (2010)  and Guerra et al. (2011)  calculated the Greek letters of currency option and stock option under uncertainty environment, respectively, where the Greek letters are useful tools for managing option risk for an option writer. But the division operation of fuzzy numbers in Liu , Xu et al.  is invalid occasionally. In this paper, we adopt the extension principle to obtain the fuzzy version of G-K formula for European currency option.

Owing to the vague of the currency market data as we just described above, some input variables, especially the interest rate, the spot exchange rate, and the volatility in the G-K formula, cannot always be expected in a precise sense. Therefore, it is reasonable to assume these input variables in the G-K formula are fuzzy numbers. Then, the fuzzy version of the G-K formula for pricing European currency option is obtained via the extension principle. It is found that one of the conditions of the basic proposition in  is not verified. Actually, it is difficult to be verified. So we modify the conditions such that they can be satisfied easily. Based on the modified conditions and the extension principle, we prove that the fuzzy price obtained from the FG-K formula for European currency option is a fuzzy number. Meanwhile, we obtain the explicit expressions of the endpoints of every α-level set (the closed interval) of the fuzzy price for European currency option. In so doing, the investor can pick any value from some closed interval (α-level set) of the fuzzy price with an acceptable belief degree α as European currency option trading price. To simplify the trade convenient, the weighted possibilistic mean (WPM) value with weighting function  is adopted to defuzzify the fuzzy price. Usually, the process of defuzzification is to find a crisp number that synthesizes the fuzzy price. Compared with the methods in [16, 18], the numerical example shows our method makes the α-level set of fuzzy price smaller, which decreases the fuzziness. The reason of the differences is also discussed. The example also indicates that the WPM value has different approximation effects to real market price by taking different values of weighting parameter in the weighting function. Inspired by this example, we provide a method, which can identify the the optimal parameter.

The rest of the paper is organized as follows. In Section 2, some terminologies, notations of fuzzy number, and the main results of [16, 18] are introduced. The fundamental theories about the fuzzy-valued functions of fuzzy subsets are proposed in Section 3. In Section 4, the fuzzy version of the G-K formula for European currency option is obtained by the extension principle. About the fuzzy price for European currency option, defuzzification via weighting parameter identification and numerical analysis are discussed in Section 5. Finally, the conclusions are stated in Section 6.

2. Preliminaries

In this section, we will introduce some definitions and propositions that will be used in the sequel. In order to weaken the conditions of the basic proposition in , we modify the definition of fuzzy number by another definitions of the 0-level set and the support set of fuzzy subset. Let X be the universe of discourse. Throughout this paper, the universe set X is assumed to be the set of all real number () endowed with a usual topology.

Definition 1 (see [<xref ref-type="bibr" rid="B23">9</xref>]).

A fuzzy subset A~ of X is a set of ordered pairs A~={(x,μA~(x))xX}, where μA~:X[0,1] is called membership function of A~.

The concept of fuzzy subset was first introduced by . The value μA~(x) can be interpreted as the membership degree of a point x to the set A~.

Definition 2 (see [<xref ref-type="bibr" rid="B23">9</xref>]).

A fuzzy subset A~ of X is called normal if there exists at least one element xX such that μA~(x)=1.

Definition 3 (see [<xref ref-type="bibr" rid="B7">21</xref>]).

Let A~ be a fuzzy subset of X. The support of A~, denoted by Supp(A~), is the crisp set cl{xXμA~(x)>0}.

Definition 4 (see [<xref ref-type="bibr" rid="B23">9</xref>]).

Let A~ be a fuzzy subset of X. The α-level set of A~, denoted by A~α, is the crisp set {xXμA~(x)α}, where α[0,1].

Definition 5 (see [<xref ref-type="bibr" rid="B23">9</xref>]).

A fuzzy subset A~ of X is called convex if for all α[0,1] its α-level sets A~α are (crisp) convex sets.

Alternatively, a fuzzy subset A~ of X is a convex fuzzy subset if and only if μA~(λx+(1-λ)y)min{μA~(x),μA~(y)} for all λ[0,1]; that is, μA~(x) is a quasi-concave function.

Definition 6 (see [<xref ref-type="bibr" rid="B4">22</xref>]).

A real-valued function f: is called upper semicontinuous at a point x if f(x)limsupkf(xk) for every sequence {xk} that converges to x. If f: is upper semicontinuous at every point in , we say that f(x) is upper semicontinuous.

A real-valued function f(x) is upper semicontinuous if and only if {xf(x)α} is closed for all α .

Definition 7 (fuzzy number).

Let a~ be a fuzzy subset of . Then a~ is called a fuzzy number if the following conditions are satisfied:

a~ is a normal and convex fuzzy subset,

its membership function μa~(x) is upper semicontinuous,

the Supp(a~) is bounded.

Note that if a~ is a fuzzy number, then the α-level set a~α, for all α(0,1], is a compact (closed and bounded in ) and convex set; that is, a~α is a bounded and closed interval for all α(0,1]. Then, the α-level set of a~ is denoted as a~α=[a~αL,a~αU] for all α(0,1], and a~0= especially. For convenience, we consider the fuzzy input variables as triangular fuzzy numbers; that is, the graph of the membership function looks like triangles.

Definition 8.

A fuzzy number a~ is a triangular fuzzy number if its membership function μa~(x) is characterized as follows: (1)μa~(x)={x-a1a2-a1a1xa2,x-a3a2-a3a2xa3,0otherwise. Here, [a1,a3] is the support set and the membership function has a peak at a2. The triangular fuzzy number a~ usually is denoted as a~=(a1,a2,a3).

Definition 9 (see [<xref ref-type="bibr" rid="B9">20</xref>]).

Let a~ be a fuzzy number with a~α=[a~αL,a~αU], α[0,1]. A function f:[0,1] is said to be a weighting function if f(x) is nonnegative, monotone increasing and satisfies the normalization condition: (2)01f(α)dα=1. We define the f-weighted possibilistic mean (WPM) value of the fuzzy number a~ as (3)M-f(a~)=01a~αL+a~αU2f(α)dα.

Next, the decomposition theory is introduced which establishes an important connection between fuzzy sets and crisp sets. That is, we can get the expression of membership function of a fuzzy set if all its α-level sets are available.

Proposition 10 (decomposition theorem [<xref ref-type="bibr" rid="B24a">23</xref>–<xref ref-type="bibr" rid="B24c">25</xref>]).

Let A~ be a fuzzy subset with membership function μA~(x) and α-level set A~α={xμA~(x)α}. Then (4)μA~(x)=sup0α1α·1A~α(x), where 1A~α(x) is an indicator function of set A~α, that is, 1A~α(x)=1 if xA~α and 1A~α(x)=0 if xA~α.

Let denote all the fuzzy subsets of . Assume A~i, i=1,2,,n. Let A~A~1×A~2××A~n and Supp(A~)Supp(A~1)×Supp(A~2)××Supp(A~n), where × is Cartesian product. A fruitful and powerful tool of fuzzy set theory for calculating the fuzzy-valued functions on fuzzy sets is the extension principle.

Proposition 11 (extension principle [<xref ref-type="bibr" rid="B24a">23</xref>–<xref ref-type="bibr" rid="B24c">25</xref>]).

Let f:n be a real-valued function and let A~1,A~2,,A~n be fuzzy subsets of . Then one can induce a fuzzy-valued function f~:n according to the real-valued function f(x1,x2,,xn). In fact, f~(A~1,A~2,,A~n) is a fuzzy subset of . Then the membership function of f~(A~1,A~2,,A~n) is characterized as follows: (5)μf~(A~1,A~2,,A~n)(r)=sup{(x1,,xn)r=f(x1,,xn)} Supp (A~)min1in{μA~i(xi)}, where supxg(x)=0.

Remark 12.

It is necessary to point out that there is a term Supp(A~) in (5), which does not exist in . Since supxg(x)=0, the term Supp(A~) does not affect the proposition. However, it plays an important role in our results.

Proposition 13 (see [<xref ref-type="bibr" rid="B4">22</xref>]).

Let S be a nonempty compact set in n. If f is upper semicontinuous on S, then f attains maximum over S and if f is lower semicontinuous on S, then f attains minimum over S.

Proposition 14.

(i) [26, Theorem 4.25] Let f:ST be a function from one metric space (S,dS) to another (T,dT). If f is continuous on a compact subset X of S, then the image f(X) is a compact subset of T; in particular, f(X) is a closed and bounded in T.

(ii) [26, Theorem 4.37] Let f:ST be a function from one metric space (S,dS) to another (T,dT). Let X be a connected subset of S. If f is continuous on X, then the f(X) is a connected subset of T.

We also review the results in [16, 18] about the formula for currency option under uncertain environment as follows: (6)(C~t)αL=(Q~t)αL·e-(r~f)αU×(T-t)·N  ((d~1)αL)-K·e-(r~d)αL×(T-t)·N((d~2)αU),(C~t)αU=(Q~t)αU·e-(r~f)αL×(T-t)·N((d~1)αU)-K·e-(r~d)αU×(T-t)·N((d~2)αL), where (7)(d~1)αL=(ln((Q~t)αLK)+[(r~d)αL-(r~f)αU+12((σ~Q)αL)2]×(T-t)((Q~t)αLK))((σ~Q)αUT-t)-1,(d~1)αU=(ln((Q~t)αUK)  +[(r~d)αU-(r~f)αL+12((σ~Q)αU)2]×(T-t)((Q~t)αUK))((σ~Q)αLT-t)-1,(d~2)αL=(d~1)αL-(σ~Q)αUT-t,(d~2)αU=(d~1)αU-(σ~Q)αLT-t.

3. Expressing of <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M151"><mml:mrow><mml:mi>α</mml:mi></mml:mrow></mml:math></inline-formula>-Level Set of the Fuzzy-Valued Functions of Fuzzy Numbers

In this section, we will establish the fundamental proposition for this paper, which based on the decomposition theorem, the extension principle, and some properties about continuous function. The following proposition is a result about the fuzzy-valued functions of fuzzy subsets.

Proposition 15.

Let f:n be a continuous, surjective, real-valued function and let A~1,A~2,,A~n be fuzzy subsets of . Let f~:n be a fuzzy-valued function induced by f(x1,x2,,xn) via the extension principle defined in (5). Suppose that

μA~i(xi)  is upper semicontinuous on , i=1,2,,n,

Supp(A~)  is a bounded set of n,

then (8)(a)(f~(A~1,A~2,,A~n))α={f(x1,x2,,xn)x1(A~1)α,x2(A~2)α,,xn(A~n)α},α[0,1],(9)(b) Supp (f~(A~1,A~2,,A~n))={f(x1,x2,,xn)x1 Supp (A~1),x2 Supp (A~2),,xn Supp (A~n)}.

Proof.

(a) For α=0, (8) holds from Definition 4 and the property of surjection of f(x1,x2,,xn).

For α(0,1], let r{f(x1,x2,,xn)x1(A~1)α,x2(A~2)α,,xn(A~n)α}. Then there exists (x1*,x2*,,xn*)(A~1)α×(A~2)α××(A~n)α such that r=f(x1*,x2*,,xn*). Thus min1in{μA~i(xi*)}α. It follows that (10)μf~(A~1,A~2,,A~n)(r)=sup{(x1,,xn)r=f(x1,,xn)}Supp(A~)min1in{μA~i(xi)}α, which means that (11){f(x1,x2,,xn)x1(A~1)α,x2(A~2)α,,xn(A~n)α}(f~(A~1,A~2,,A~n))α. On the other hand, μA~i(xi) is upper semicontinuous on , then Ui={xiμA~i(xi)α} is a closed subset of , i=1,2,,n. Therefore, Ci=××i-1×Ui×××n-i is a closed subset of ni=1,2,,n. Since (12){(x1,x2,,xn)min1in{μA~i(xi)}α}={(x1,x2,,xn)μA~i(xi)α,i=1,2,,n}=i=1nCi is a closed subset of n, we have min1in{μA~i(xi)} is upper semicontinuous on n. Let r(f~(A~1,A~2,,A~n))α, that is, (13)sup{(x1,,xn)r=f(x1,,xn)}Supp(A~)min1in{μA~i(xi)}α. Since {(x1,,xn)r=f(x1,,xn)} is a closed subset of nand Supp(A~) is a compact set of n and min1in{μA~i(xi)} is upper semicontinuous on n, we conclude from Proposition 13 that there exists (x-1,x-2,,x-n) such that r=f(x-1,x-2,,x-n) and min1in{μA~i(x-i)}α. Therefore μA~i(x-i)α, that is, x-i(A~i)α,  i=1,2,,n. This means (14)r{f(x1,x2,,xn)x1(A~1)α,x2(A~2)α,,xn(A~n)α}. Therefore, (15)(f~(A~1,A~2,,A~n))α{f(x1,x2,,xn)x1(A~1)α,x2(A~2)α,,xn(A~n)α} which completes the proof of (a).

(b) Let rSupp(f~(A~1,A~2,,A~n)), we conclude from Definition 3 that there exists a sequence {rm}m=1 which satisfies rmr and μf~(A~1,A~2,,A~n)(rm)>0, that is, (16)sup{(x1,,xn)rm=f(x1,,xn)}Supp(A~)min1in{μA~i(xi)}>0. Since {(x1,,xn)rm=f(x1,,xn)} is closed and Supp(A~) is bounded and closed, then {(x1,,xn)rm=f(x1,,xn)}Supp(A~) is compact. Moreover, we deduce that min1in{μA~i(xi)} is upper semicontinuous on n according to the proof of (a). Consequently, we conclude from Proposition 13 that there exist {x-1m}m=1,{x-2m}m=1,,{x-nm}m=1 such that (17)μf~(A~1,A~2,,A~n)(rm)=min1in{μA~i(x-im)}>0formN,i=1,2,,n, where rm=f(x-1m,,x-nm) and x-imSupp(A~i). Since Supp(A~){(x1,,xn)rm=f(x1,,xn)} is compact, for each i{1,2,,n}, there exists a subsequence {x-imk}k=1 of {x-im}m=1 such that x-imkx-iSupp(A~i) as k. Notice that (18)rmk=f(x-1mk,,x-nmk)r,f(x-1mk,,x-nmk)f(x-1,,x-n)ask. We have r=f(x-1,,x-n), where x-iSupp(A~i). Then (19)Supp(f~(A~1,A~2,,A~n)){f(x1,x2,,xn)x1Supp(A~1),,xnSupp(A~n)}. For the converse, let z{f(x1,x2,,xn)x1Supp(A~1),,xnSupp(A~n)}. Then there exist sequences {xim}m=1, i=1,2,,n, which satisfy μA~i(xim)>0 and ximxi as m, where xiSupp(A~i) and z=f(x1,x2,,xn). Recall the continuity of f(x1,x2,,xn), we have (20){f(x1,x2,,xn)x1Supp(A~1),,xnSupp(A~n)}Supp(f~(A~1,A~2,,A~n)). This completes the proof.

Remark 16.

One of the conditions of Proposition in  is that {(x1,,xn)r=f(x1,,xn)} is compact forallrRange(f). This one is not verified there. Actually, it is difficult to be verified. Meanwhile, (8) only holds for α(0,1]. We change the above conditions to (i) μA~i(xi) is upper semicontinuous on , i=1,2,,n; (ii) Supp(A~) is bounded of n and; (iii) f:n is a continuous, surjective, real-valued function. Conditions (i)–(iii) can be easily satisfied. At the same time, (8) holds for α[0,1] and (9) is also obtained.

The following is the fundamental proposition for discussing the Garman-Kohlhagen formula for European currency option via the extension principle.

Proposition 17.

Let f:n be a continuous, surjective, real-valued function and let a~1,a~2,,a~n be fuzzy numbers. Let f~:n be a fuzzy-valued function induced by f(x1,x2,,xn) via the extension principle defined in (5). Then f~(a~1,a~2,,a~n) is a fuzzy number and its α-level set is (21)(f~(a~1,a~2,,a~n))α={f(x1,x2,,xn)x1(a~1)α,x2(a~2)α,,xn(a~n)α}={f(x1,x2,,xn)(a~1)αLx1(a~1)αU,,(a~n)αLxn(a~n)αU}.

Proof.

Let Xα{(x1,x2,,xn)x1(a~1)α,x2(a~2)α,,xn(a~n)α} for all α(0,1]. Then it follows from Definition 7 that Xα is an n-dimensional bounded interval; that is, Xα is a compact and connected subset of n. Next, we conclude from Propositions 15 and 14 that

{μf~(a~1,a~2,,a~n)α}=(f~(a~1,a~2,,a~n))α={f(x1,x2,,xn)(x1,x2,,xn)Xα} is a compact and connected subset of for all α(0,1];

(f~(a~1,a~2,,a~n))0= is closed and convex set;

Supp(f~(a~1,a~2,,a~n))={f(x1,x2,,xn)x1Supp(a~1),x2Supp(a~2),,xnSupp(a~n)}.

Furthermore,

the membership function μf~(a~1,a~2,,a~n) of f~(a~1,a~2,,a~n) is upper semicontinuous;

f~(a~1,a~2,,a~n) is a convex fuzzy subset;

Supp(f~(a~1,a~2,,a~n)) is bounded.

Obviously, f~(a~1,a~2,,a~n) is a normal fuzzy subset, as X1 is nonempty. Therefore, from Definition 7 we conclude that f~(a~1,a~2,,a~n) is a fuzzy number and its α-level set is (21). This completes the proof.

4. Fuzzy Currency Options Pricing Formula via the Extension Principle 4.1. The G-K Formula for European Currency Option

The closed-form solution for European currency option pricing formula has been derived by Garman and Kohlhagen . The G-K formula for European call currency option with the expiry date T and the strike price K is described as follows. Let Qt, Ct denote the spot exchange rate, and the price of European call currency option at time t[0,T], respectively, then (22)Ct=Qt·e-rfτ·N(d1)-K·e-rdτ·N(d2), where τ=T-t, d1=[ln(Qt/K)+(rd-rf+(1/2)σQ2)τ]/(σQτ), d2=d1-σQτ, rf, rd, σQ, N(·) denote the foreign interest rate, the domestic interest rate, the volatility of spot exchange rate and the standard Gaussian cumulative distribution function, respectively. Let (23)f(q,t,K,rd,rf,σ)=q·e-rft·N(d1)-K·e-rdt·N(d2), then (24)Ct=f(Qt,T-t,K,rd,rf,σQ), for all t[0,T]. Furthermore, the price Pt of European put option at time t[0,T] with same expiry date T and strike price K can be obtained by the following put-call parity relationship (Musiela and Rutkowski ): (25)Ct-Pt=e-rf(T-t)·Qt-e-rd(T-t)·K. Let (26)g(q,t,K,rd,rf,σ)=f(q,t,K,rd,rf,σ)-e-rft·q+e-rdt·K for all t[0,T], then (27)Pt=g(Qt,T-t,K,rd,rf,σQ).

4.2. The Fuzzy Version of the G-K Formula for European Currency Option

As the discussion in Section 1, it is suitable to represent the four input variables in the G-K formula by fuzzy numbers. Thus, we give the following assumptions.

Assumption 1.

The spot exchange rate at time t[0,T] is fuzzy number Q~t;

the domestic interest rate is fuzzy number r~d;

the foreign interest rate is fuzzy number r~f;

the volatility of spot exchange rate is fuzzy number σ~Q.

Based on Assumption 1, (24) and the extension principle (Proposition 11), the original price for European call currency option Ct at time t[0,T] turns into fuzzy subset C~t with the membership function (28)μC~t(c)=sup{(q,rd,rf,σQ)c=f(q,T-t,K,rd,rf,σQ)}Supp(A~)×min{μQ~t(q),μr~d(rd),μr~f(rf),μσ~Q(σQ)}, where Supp(A~)=Supp(Q~t)×Supp(r~d)×Supp(r~f)×Supp(σ~Q). Since the function f(Qt,T-t,K,rd,rf,σQ) satisfies conditions of Proposition 17, we conclude that the fuzzy price C~t at time t[0,T] is a fuzzy number. Furthermore, its α-level set (C~t)α is (29)(C~t)α={f(q,T-t,K,rd,rf,σQ)q(Q~t)α,rd(r~d)α,rf(r~f)α,σQ(σ~Q)α}, for all α[0,1]. On the other hand, from Definition 7, the α-level set (C~t)α of C~t is a closed interval (30)(C~t)α=[(C~t)αL,(C~t)αU]. Notice that the α-level set of Q~t, r~d, r~f, and σ~Q are (Q~t)α=[(Q~t)αL,(Q~t)αU], (r~d)α=[(r~d)αL,(r~d)αU], (r~f)α=[(r~f)αL,(r~f)αU], and (σ~Q)α=[(σ~Q)αL,(σ~Q)αU], respectively. Therefore, from (29) and (30), the left endpoint (C~t)αL and the right-end point (C~t)αU can be displayed as (31)(C~t)αL=min(Q~t)αLq(Q~t)αU,(r~d)αLrd(r~d)αU,(r~f)αLrf(r~f)αU,(σ~Q)αLσQ(σ~Q)αUf(q,T-t,K,rd,rf,σQ),(C~t)αU=max(Q~t)αLq(Q~t)αU,(r~d)αLrd(r~d)αU,(r~f)αLrf(r~f)αU,(σ~Q)αLσQ(σ~Q)αUf(q,T-t,K,rd,rf,σQ). Moreover, according to the sensitivity analysis (Greek letters) of the G-K formula , we have (32)fq=e-rftN(d1)>0,frd=t·K·e-rdtN(d2)>0,frf=-t·q·e-rftN(d1)<0,fσQ=q·e-rftN(d1)·t>0. This mean f is increasing with respect to q, rd, and σQ and is decreasing with respect to rf. Therefore, from (31), we conclude that (33)(C~t)αL=f((Q~t)αL,T-t,K,(r~d)αL,(r~f)αU,(σ~Q)αL),(C~t)αU=f((Q~t)αU,T-t,K,(r~d)αU,(r~f)αL,(σ~Q)αU). On the other hand, from Proposition 10, the membership function of C~t can also be written as (34)μC~t(c)=sup0α1α·1(C~t)α(c). Finally, (28), (33), and (34) constitute the fuzzy version of the G-K formula for European call currency option.

Similarly, the conclusions about European put currency option are as follows, through the above discussions.

The original price of European put currency option Pt at time t[0,T] turns into fuzzy subset P~t with the membership function under Assumption 1(35)μP~t(p)=sup{(q,rd,rf,σQ)p=g(q,T-t,K,rd,rf,σQ)}Supp(A~)×min{μQ~t(q),μr~d(rd),μr~f(rf),μσ~Q(σQ)}, where Supp(A~) is as above. Furthermore, the fuzzy price P~t at time t[0,T] is a fuzzy number and the α-level set (P~t)α of P~t is a closed interval (36)(P~t)α=[(P~t)αL,(P~t)αU]. Moreover, the left endpoint (P~t)αL and the right endpoint (P~t)αU can be displayed as (37)(P~t)αL=g((Q~t)αU,T-t,K,(r~d)αU,(r~f)αL,(σ~Q)αL),(P~t)αU=g((Q~t)αL,T-t,K,(r~d)αL,(r~f)αU,(σ~Q)αU). On the other hand, the membership function of P~t can also be written as (38)μP~t(p)=sup0α1α·1(P~t)α(p).

Finally, (35), (37), and (38) constitute the fuzzy version of the G-K formula for European put currency option.

5. Defuzzification via Weighting Parameter Identification and Numerical Analysis

In this section, we defuzzify the fuzzy price for European currency option C~t(P~t). Then an numerical example is presented. Up to now, it is proved that the fuzzy price for European currency option is a fuzzy number under Assumption 1. For practical purposes, it may be convenient to find a crisp number that synthesizes the fuzzy number C~t(P~t). This type of problem is known as the “defuzzification procedure" in the literature. Defuzzification is the process to select an appropriate crisp value based on a fuzzy set such that the selected crisp value may represent the fuzzy set in some sense.

Of all the methods of defuzzification, two most popular methods are the center of area (COA) and the mean of maxima (MOM) which require explicit membership function. Now, the membership function of C~t(P~t) cannot be expressed explicitly; it can be deduced by the decomposition theorem (see (34) or (38)) or the extension principle (see (28) or (35)) only. Since its α-level sets of C~t(P~t) are available (see (33) or (37)), the weighted possibilistic mean value of Fullér and Majlender  about fuzzy number is fited to be the crisp number that we needed (see Definition 9) in the following.

Next, the proposed fuzzy version of the G-K formula (28) is tested with the daily market price data of EUR/USD currency option. To compare with the results in [16, 18] conveniently, the market data from  are used in Example 18. In the following example, we assumed that the weighting function is fn(x)(1+n)αn, where 0α1 is the belief degree and n is the belief degree of weighting parameter.

Example 18.

A European EUR/USD currency call option is studied on march 16, 2006. For the EUR/USD currency option, the spot exchange rate is 1.215, the 3-month volatility is around 9.0%, the domestic 3-month interest rate is around 4.93%, the foreign 3-month interest rate is around 2.71%, and the strike price is 1.21 with 3-month to expiry. The market price of this currency option is 0.0274 USD. The fuzzy foreign and domestic interest rates, the fuzzy spot exchange rate, and the fuzzy volatility are assumed to triangular fuzzy numbers: Q~t=(1.2138,1.2150,1.2162), σ~Q=(7.2%,9.0%,10.8%), r~d=(4.91%,4.93%,4.95%), and r~f=(2.69%,2.71%,2.72%), respectively.

The fuzzy price C~t of European call currency option can be obtained by (28). Table 1 gives the α-level set of the fuzzy price C~t for different belief degrees α and the f-weighted possibilistic mean value M-fn(C~t) for different belief degrees of weighting parameters n, respectively. In fact, if Qt=1.215, rd=4.93%, rf=2.71%, and σQ=9.0%, then European call currency option price is 0.027872\$ by the G-K formula (24). This situation matches the observation Ct=0.027872\$ with belief degree 1.

Fuzzy EUR/USD currency option price.

α -level set of this paper α -level set of paper [16, 18] WPM Market price
α [ ( C ~ t ) α L , ( C ~ t ) α U ] [ ( C ~ t ) α L , ( C ~ t ) α U ] n M - f n ( C ~ t )
0.0 [ 0.022898,0.032859 ] [ - 0.05526,0.11108 ] 0 0.02787
0.1 [ 0.023394,0.032360 ] [ - 0.04644,0.10225 ] 1 0.027589
0.2 [ 0.023890,0.031860 ] [ - 0.037785,0.093576 ] 2 0.02745
0.3 [ 0.024386,0.031360 ] [ - 0.029274,0.085052 ] 3 0.027312
0.4 [ 0.024883,0.030860 ] [ - 0.020886,0.076652 ] 4 0.027174
0.5 [ 0.025379,0.030361 ] [ - 0.012604,0.068359 ] 5 0.027036 0.0274
0.6 [ 0.025876,0.029861 ] [ - 0.0044076,0.060155 ] 6 0.0269
0.7 [ 0.026373,0.029362 ] [ 0.0037199,0.052020 ] 7 0.026763
0.8 [ 0.026870,0.028862 ] [ 0.011796,0.043938 ] 8 0.026628
0.9 [ 0.027367,0.028363 ] [ 0.019839,0.035892 ] 9 0.026492
1.0 [ 0.027864,0.027864 ] [ 0.027864,0.027864 ] 10 0.026358
Remark 19.

There are some differences between the α-level sets of fuzzy call option price in this paper and those in [16, 18], which are shown in Table 1. The reason is that the division operation of fuzzy numbers in [16, 18] is invalid occasionally. Table 1 also shows that M-fn(C~t) is the best approximation to the real price when n=2.

Next, assuming the currency option prices Ct1,Ct2,,CtN at t1,t2,,tN dates are available, how to acquire a more accurate CtN+1 at tN+1 date is an important question. Based on minimizing the mean square error between the market prices and WPMs, a method of weighting parameter identification is presented as follows.

Frame:

collect the market price data: Ct1,Ct2,,CtN;

calculate the performance index J(n)=i=1N(01(((C~ti)αL+(C~ti)αU)/2)(1+n)αndα-Cti)2;

solve min0n1J(n) and optimal solution n* is obtained and;

calculate M-fn(C~tN+1)=01(((C~tN+1)αL+(C~tN+1)αU)/2)(1+n*)αn*dα, and let CtN+1=M-fn(C~tN+1).

6. Conclusions

In this paper, the Fuzzy version of the Garman-Kohlhagen (FG-K) formula for pricing European currency option via the extension principle is proposed. In the FG-K formula, the interest rate, the spot exchange rate, and the volatility in conventional G-K formula are replaced by fuzzy numbers. So the FG-K formula is able to keep consistent with the real market. We prove that the fuzzy price for European currency option is a fuzzy number. And then the fuzzy price C~t(P~t) is defuzzified by the weighted possibilistic mean value. The example shows that the fuzzy price in this paper is more accurate than the ones in [16, 18] and the fuzzy prices in [16, 18] are just approximations to our fuzzy price for European currency option. In addition, the example indicates that the WPM value has different approximation effects to real market price by taking different values of weighting parameter in the weighting function. So an identification method is presented. And the optimal parameter is obtained if we have already known some market prices for European currency option. The method in this paper can be extended to other option pricing. How to find out different kinds of weighting function in defuzzification requires a further study.

Acknowledgments

This work was supported by 863 Program (Grant no. 2007AA02Z208), 973 Program (Grant no. 2007CB714304), and the National Natural Science Foundation of China (Grants nos. 10871033 and 11171050).

Garman M. B. Kohlhagen S. W. Foreign currency option values Journal of International Money and Finance 1983 2 3 231 237 10.1016/S0261-5606(83)80001-1 Black F. Scholes M. The pricing of options and corporate liabilities Journal of Political Economics 1973 81 3 637 659 10.1086/260062 Chesney M. Scott L. Pricing european currency options: a comparision of the modified black-scholes model and a random variance model Journal of Financial and Quantitative Analysis 1989 24 3 267 284 10.2307/2330812 Amin K. Jarrow R. Pricing foreign currency options under stochastic interest rates Journal of International Money and Finance 1991 10 3 310 329 10.1016/0261-5606(91)90013-A Heston S. A closed-form solution for options with stochastic volatility with applications to bond and currency options Review of Financial Studies 1993 6 2 327 344 10.1093/rfs/6.2.327 Bates D. Jumps and stochastic volatility: exchange rate processes implicit in deutsche mark options Review of Financial Studies 1996 9 1 69 107 10.1093/rfs/9.1.69 Sarwar G. Krehbiel T. Empirical performance of alternative pricing models of currency options The Journal of Futures Markets 2000 20 3 265 291 10.1002/(SICI)1096-9934(200003)20:3&lt;265::AID-FUT4&gt;3.0.CO;2-4 Carr P. Wu L. R. Stochastic skew in currency options Journal of Financial Economics 2007 86 1 213 247 10.1016/j.jfineco.2006.03.010 Zadeh L. A. Fuzzy sets Information and Computation 1965 8 338 353 MR0219427 ZBL0139.24606 Yoshida Y. The valuation of European options in uncertain environment European Journal of Operational Research 2003 145 1 221 229 10.1016/S0377-2217(02)00209-6 MR1947161 ZBL1011.91045 Yoshida Y. A discrete-time model of American put option in an uncertain environment European Journal of Operational Research 2003 151 1 153 166 10.1016/S0377-2217(02)00591-X MR1997013 ZBL1112.91328 Muzzioli S. Torricelli C. A multiperiod binomial model for pricing options in a vague world Journal of Economic Dynamics & Control 2004 28 5 861 887 10.1016/S0165-1889(03)00060-5 MR2030666 Wu H. C. Using fuzzy sets theory and Black-Scholes formula to generate pricing boundaries of European options Applied Mathematics and Computation 2007 185 1 136 146 10.1016/j.amc.2006.07.015 MR2298436 Xu W. D. Wu C. F. Xu W. J. Li H. Y. A jump-diffusion model for option pricing under fuzzy environments Insurance: Mathematics & Economics 2009 44 3 337 344 10.1016/j.insmatheco.2008.09.003 MR2519077 ZBL1162.91388 Zhang L. H. Zhang W. G. Xu W. J. Xiao W. L. The double exponential jump diffusion model for pricing European options under fuzzy environments Economic Modelling 2012 29 3 780 786 10.1016/j.econmod.2012.02.005 Liu F. Y. Pricing currency options based on fuzzy techniques European Journal of Operational Research 2009 193 2 530 540 10.1016/j.ejor.2007.10.059 MR2463282 ZBL1153.91544 Wu H. C. European option pricing under fuzzy environments International Journal of Intelligent Systems 2005 20 1 89 102 10.1002/int.20055 Xu W. J. Xu W. D. Li H. Y. Zhang W. G. A study of Greek letters of currency option under uncertainty environments Mathematical and Computer Modelling 2010 51 5-6 670 681 10.1016/j.mcm.2009.10.041 MR2594717 ZBL1190.91144 Guerra M. L. Sorini L. Stefanini L. Option price sensitivities through fuzzy numbers Computers & Mathematics with Applications 2011 61 3 515 526 10.1016/j.camwa.2010.11.024 MR2764045 ZBL1217.91219 Fullér R. Majlender P. On weighted possibilistic mean and variance of fuzzy numbers Fuzzy Sets and Systems 2003 136 3 363 374 10.1016/S0165-0114(02)00216-6 MR1984582 ZBL1022.94032 Chalco-Cano Y. Román-Flores H. Rojas-Medar M. Saavedra O. R. Jiménez-Gamero M. D. The extension principle and a decomposition of fuzzy sets Information Sciences 2007 177 23 5394 5403 10.1016/j.ins.2007.07.007 MR2357263 ZBL1124.03326 Bertsekas D. P. Nedić A. Ozdaglar A. E. Convex Analysis and Optimization 2003 Belmont, Mass, USA Athena Scientific MR2184037 ZBL1140.90001 Zadeh L. A. The concept of a linguistic variable and its application to approximate reasoning—I Information Sciences 1975 8 3 199 249 10.1016/0020-0255(75)90036-5 Zadeh L. A. The concept of a linguistic variable and its application to approximate reasoning—II Information Sciences 1975 8 3 301 357 Zadeh L. A. The concept of a linguistic variableand its application to approximate reasoning—III Information Sciences 1975 9 1 43 80 MR0386371 ZBL0404.68075 Apostol T. M. Mathematical Analysis 1974 2nd Reading, Mass, USA Addison-Wesley MR0344384 ZBL0365.60096 Musiela M. Rutkowski M. Martingale Methods in Financial Modelling 2005 36 2nd Berlin, Germany Springer Stochastic Modelling and Applied Probability MR2107822