Research article Special Issues

Contribution of a power multivector to distorting load identification

  • Received: 21 September 2022 Revised: 08 March 2023 Accepted: 14 March 2023 Published: 17 March 2023
  • The identification of harmonic generating loads and the assignation of responsibility for harmonic pollution is an important first step for harmonic control in modern power systems. In this paper, a previously introduced power multivector is examined as a possible tool for the identification of such loads. This representation of power is based on the mathematical framework of Geometric Algebra (GA). Components of the power multivector derived at the point of connection of a load are grouped into a single quantity, which is a bivector in GA and is characterized by a magnitude, direction and sense. The magnitude of this bivector can serve as an indicator of the distortion at the terminals of the load. Furthermore, in contrast to indices based solely on magnitude, such as components derived from any apparent power equation, the proposed bivectorial representation can differentiate between loads that enhance distortion and those with a mitigating effect. Its conservative nature permits an association between the distortion at specific load terminals and the common point of connection. When several loads connected along a distribution line are considered, then an evaluation of the impact of each one of these loads on the distortion at a specific point is possible. Simulation results confirm that information included in the proposed bivector can provide helpful guidance when quantities derived from apparent power equations deliver ambiguous results.

    Citation: Anthoula Menti, Dimitrios Barkas, Pavlos Pachos, Constantinos S. Psomopoulos. Contribution of a power multivector to distorting load identification[J]. AIMS Energy, 2023, 11(2): 271-292. doi: 10.3934/energy.2023015

    Related Papers:

    [1] Shancy Petsel Jacob, Chikkamenahalli Lakshminarayana Lakshmikanth, Thomas M. McIntyre, Gopal Kedihitlu Marathe . Platelet-activating factor and oxidized phosphatidylcholines do not suppress endotoxin-induced pro-inflammatory signaling among human myeloid and endothelial cells. AIMS Allergy and Immunology, 2017, 1(3): 108-123. doi: 10.3934/Allergy.2017.3.108
    [2] Ryuta Muromoto, Kenji Oritani, Tadashi Matsuda . IL-17 signaling is regulated through intrinsic stability control of mRNA during inflammation. AIMS Allergy and Immunology, 2022, 6(3): 188-199. doi: 10.3934/Allergy.2022014
    [3] Manal Alkan, Fadel Sayes, Abdulraouf Ramadan, Francois Machavoine, Michel Dy, Elke Schneider, Nathalie Thieblemont . Basophil activation through TLR2 and TLR4 signaling pathways. AIMS Allergy and Immunology, 2018, 2(3): 126-140. doi: 10.3934/Allergy.2018.3.126
    [4] Fernanda Costa Brandão Berti, Karen Brajão de Oliveira . IL-10 in cancer: Just a classical immunosuppressive factor or also an immunostimulating one?. AIMS Allergy and Immunology, 2018, 2(2): 88-97. doi: 10.3934/Allergy.2018.2.88
    [5] Hiroto Matsuse, Tohru Yamagishi, Norio Kodaka, Chihiro Nakano, Chizu Fukushima, Yasushi Obase, Hiroshi Mukae . Therapeutic modality of plasmacytoid dendritic cells in a murine model of Aspergillus fumigatus sensitized and infected asthma. AIMS Allergy and Immunology, 2017, 1(4): 232-241. doi: 10.3934/Allergy.2017.4.232
    [6] Sabrine Louhaichi, Mariem Salhi, Anissa Berraïes, Besma Hamdi, Jamel Ammar, Kamel Hamzaoui, Agnès Hamzaoui . Co-inhibitory receptors in female asthmatic patients: Correlation with IL-17 and IL-26. AIMS Allergy and Immunology, 2018, 2(1): 10-23. doi: 10.3934/Allergy.2018.1.10
    [7] Stephen C Allen . Cytokine signaling in the modulation of post-acute and chronic systemic inflammation: a review of the influence of exercise and certain drugs. AIMS Allergy and Immunology, 2020, 4(4): 100-116. doi: 10.3934/Allergy.2020009
    [8] Izabela Krzemień, Angelika Szatan, Paulina Skalska, Martyna Cieślik, Angelika Domagała, Magdalena Gębicka, Krzysztof Bryniarski, Katarzyna Nazimek . Serum cytokine concentrations correlate with the severity of imiquimod and IL-23-induced psoriatic inflammation in mice. AIMS Allergy and Immunology, 2024, 8(4): 265-282. doi: 10.3934/Allergy.2024016
    [9] Mansur Aliyu, Sayed-Hamidreza Mozhgani, Omid Kohandel Gargari, Mustapha Ahmed Yusuf, Ali Akbar Saboor-Yaraghi . The host immune responses to SARS-CoV-2 and therapeutic strategies in the treatment of COVID-19 cytokine storm. AIMS Allergy and Immunology, 2021, 5(4): 240-258. doi: 10.3934/Allergy.2021018
    [10] Kamyar Nasiri, Zohreh Mortezania, Sanaz Oftadehbalani, Mohammad Khosousi Sani, Amin Daemi, Seyyedeh Touran Hosseini, Yusuf Döğüş, Zafer Yönden . Parasitology and Covid-19: epidemiological, immunological, pathological, and therapeutic aspects. AIMS Allergy and Immunology, 2023, 7(1): 92-103. doi: 10.3934/Allergy.2023007
  • The identification of harmonic generating loads and the assignation of responsibility for harmonic pollution is an important first step for harmonic control in modern power systems. In this paper, a previously introduced power multivector is examined as a possible tool for the identification of such loads. This representation of power is based on the mathematical framework of Geometric Algebra (GA). Components of the power multivector derived at the point of connection of a load are grouped into a single quantity, which is a bivector in GA and is characterized by a magnitude, direction and sense. The magnitude of this bivector can serve as an indicator of the distortion at the terminals of the load. Furthermore, in contrast to indices based solely on magnitude, such as components derived from any apparent power equation, the proposed bivectorial representation can differentiate between loads that enhance distortion and those with a mitigating effect. Its conservative nature permits an association between the distortion at specific load terminals and the common point of connection. When several loads connected along a distribution line are considered, then an evaluation of the impact of each one of these loads on the distortion at a specific point is possible. Simulation results confirm that information included in the proposed bivector can provide helpful guidance when quantities derived from apparent power equations deliver ambiguous results.



    The idea of an option is not new. Options were traded in Ancient Romans against outdoing cargoes from their seaport. The options are the main dynamic segment of the security market, since the origin of the Chicago Board option exchange (CBOF) in April 1997. It is the largest options exchange in the world with more than one million contracts per day (Boyle, 1977). The option is a type of a derivative. It is defined as "It gives the holder the rights but not obligation to buy in case of a call option or sell in the case put option of an underlying asset at a fixed price (exercise or strike price) on the maturity (expiry date)". Options are used for hedging and speculation. Hedging is nothing it is simply a process to reduce the future loss and speculation is to make a position in the market (Hull, 2006). Option theory has a long history, but it was not until Black and Scholes presented the first complete equilibrium option pricing model in the year 1973. In the same year, Robert Merton explains the BSM in different ways. From 1973 when BSM came in the market it has been widely used by traders to determine the price of the option (Oduro, 2012). Another option model existed in the market namely BM developed by Cox-Ross-Rubinstein in 1979 (Cox et al., 1979). BM is a simple and easy to understand. Jump diffusion model suggested by Merton in 1973. The Monte Carlo simulation method was developed by Boyle in the year 1977.

    Pricing of options (call/put) is one of the most important aspects of trading the derivatives. The two models BSM and BM are famous in the financial market. In 1973 Fisher Black and Myron Scholes, they develop an option pricing model called BSM for European style options without dividends which allowed investors to approximate a price. Previous models by Sprenkle in 1961, Ayes in 1963, Boness in 1964 and others had determined the value of the option in terms of warrants (Black and Scholes, 1973). A warrant is a contract that confirms the right, but not obligation to buy or sell. It is issued by a firm, not a stock exchange. The assumption of the BSM is that the underlying assets follow a lognormal distribution before BSM researchers such as Mandelbrot in 1963, Fama in 1965 and Blattberg and Gonedes in 1974 suggesting that it is following student t distribution. In 1973, Robert Merton modified the BSM to account for dividends and variables interest rate. In 1976, Jonathan Ingersoll gave an assumption for BSM that is there is no tax or transaction cost. John Cox and Stephen Ross in 1976 showed that underlying assets do not have to move continuously, but the price of the underlying goes from one price to another.

    Feng and Kwan (2012), investigated that the BM which eventually converges to well-known BSM as the number of Binomial period increase. In the year, 2012 the Oduro and Dedu showed the relationship between the BM and the BSM with the help of MS VBA. In the year 2014, Lazarova et al. investigated that the two modes BM and the BSM do not affect the crucial idea that underlies the option of each model. In the year 2017, Dar and Anuradha showed that the BM and the BSM have the negligible difference between them with the help of MATLAB by simulations.

    Different parameters of data are estimated by using MS-Excel. Initially, the basic terminologies of options are studied like some definitions of derivatives, types of options, option styles are mentioned. In BM we use lattices to estimate the value of the option (call/put) at the initial time (t = 0) that provides a payoff for each call and put options at different strike prices that we generated by using excel function at an expiry/mature date based on the non-dividend paying shares and without premium. Also, we calculated the price of each option by applying a direct formula known as Black-Sholes formula. Finally, the results we got that BSM is approximately equal to BM at one period with the help of statistical t-test and Tukey method. MINITAB and SPSS have been used as a necessary tool to calculate all the important results.

    To understand how BM and BSM are approximately equal at one period by using t-test and Turkey method or in other words, there is no mean difference between the BM and BSM at one period.

    Data have taken from State Bank of India on monthly from 1st January 2010 to 31st December 2010 monthly basis from www.finance.yahoo.com. First of all, we calculated the price of a European option at one period by using the two famous pricing models BM and BSM. Later we compare both the models in order to check whether there are any mean differences between the BM and BSM at one period with the help of statistical tools -t-test and Turkey method.

    A derivative is a contract/security which gives us promise to make payment on some future date. It is a financial tool whose value depends on some underlying assets. The underlying assets are shares, bonds, index, interest rate, currency, commodity (gold, wheat etc).

    It is a contract between individuals or firms in which one party is ready to buy and another party is ready to sell. "An option is a contract that gives the owner the right to buy or sell any underlying asset at a fixed price (the predetermined price of an underlying asset) on the future date".

    There are two types of options:

    "It gives the rights to the holder of an option but not obligation to buy an underlying asset at a strike price on the future date".

    "It gives the rights to the holder of an option but not obligation to sell an underlying asset at a strike price on the future date".

    There are two styles of the options which are defined below:

    "American option is a contract that gives the holder right to buy or sell any underlying asset at a strike price on or before the expiry date". The owner can exercise this option anytime before the expiry date.

    "European Option is a contract that gives the right to buy or sell any underlying asset at a fixed strike price on the exercise date". The European option only exercises at maturity date not before.

    Below are some notations that are using in this paper:

    t -The current time;

    St -The underlying share price at time t;

    K -The strike/exercise price;

    T -The option exercise date; ­

    r is the interest rate;

    The money realised by the holder of an option at the end of its life (Hull, 2006). The profit made by the holder at the maturity date is known as a payoff.

    Let us consider a call option with St is the price of the underlying at time t and fixed strike price K.

    At expiry time T we have two different cases:

    (a) "At expiry time T, if the price of strike price K is less than the price of underlying asset ST, then the call option is exercised i.e. the holder buys the underlying asset at price K and immediately sell it in the market at priceST, the holder realize the profit".

    Payoff=STk

    (b) "At expiry time T, if the price of strike price K is greater than the price of underlying asset ST, then the call option is not exercised". The option expires worthless with payoff=0.

    Combine the above two cases, at time T. The value of a call option is given by a function

    Cu=max(STk,0)

    Figure 1 shows payoff of a European call option without premium.

    Figure 1.  Payoffs of a European call option.

    Let us consider a put option with St is the price of the underlying at time t and fixed strike price K.

    At expiry time T we have two different cases:

    (a) "At the expiry date T, if the strike price K is greater than the underlying asset ST, then the put option is exercised, i.e. the holder buys the underlying asset from the market price ST and sell it to the writer at price K". The holder realizes the profit of K -ST.

    (b) "At expiry time T, if the strike price K is less than the price of underlying asset ST, than the put option is not exercised". The option expires worthless with payoff = 0.

    Combine the above two cases, at time T. The value of a put option is given by a function

    cd=max(KST,0)

    Figure 2 shows payoff of a European put option.

    Figure 2.  Payoff of a European put option.

    "It reduces the future risk in the underlying by entering into a derivative contract whose value moves in the opposite direction to the underlying position and cancels part or all of it out (Hull, 2006)".

    "It is a sort of winning or losing the amount of money I.e; gambling, in this a person even gains a big margin of money or loses a lot of amounts (Hull, 2006)".

    "Arbitrage means risk-free trading profit"; it is something described as a free breakfast.

    Arbitrage opportunity exists if

    (ⅰ) If we make an immediate profit with a probability of loss is zero.

    (ⅱ) With a zero cost initial investment if we can get money in future.

    "No-arbitrage means when arbitrage opportunity does not exist in the market".

    The BM is also known as Cox-Ross-Rubinstein (CRR) model because J.C Cox, S.A Ross and M. Rubinstein developed BM in the year 1979. Options play a necessary role in the financial market as widely applied in financial derivatives (Liu and Chen, 2009).BM is a model to estimate the value of an option at time t = 0 (at the beginning) that provides a payoff at a future date based on the value of a non-dividend paying shares at a future date. The BM is based on the assumption that there is no arbitrage opportunity exist in the market. The assumption is that the stock price follows a random walk. In every step, it has a certain probability of moving up or down.

    We consider a one-step BM. In this model, we start at time t = 0 when the underlying asset price is equal to St. Over a time period (say one) the underlying asset prices do one of the two things (Dar and Anuradha, 2017). It will be either

    Jump upward with some probability say q

    Jump downward with some probability say 1 -q

    At time t = 1 we have two possibilities

    S1={S0u,iftheunderlyingassetpricegoesupS0d,iftheunderlyingassetpricegoesdown

    Note: u is always greater than 1 and d is less than 1.

    Figure 3 represents the one step BM

    Figure 3.  One step BM.

    Figure 3 is one period model because we consider only time interval from 0 to 1, it is Binomial because there are only two possibilities either goes up S0u with probability q or down S0d=S1 with probability 1q. The above model is referred as Binomial tree and each path from S0 to S0u and from S0 to S0d as a branch. The points at the end of the branch are known as nodes (Dar and Anuradha, 2017).

    We can now see the implication of there being no arbitrage in this model. In order to avoid arbitrage, it is necessary that we must have d<er<u.

    The derivative price (one period) at time t=0 is

    V0=er(Cuq+Cd(1q))(1)

    Where q=erdud and 1q=uerud

    Note that the no-arbitrage conditions d<er<u ensures that 0<q<1.

    Proof of equation 1:

    We consider a one-period binomial model as shown in figure 3; we have two possibilities at time t = 1 i.e.

    S1={SouiftheunderlyingassetpricegoesupS0diftheunderlyingassetpricegoesdown

    Here u>1 and d<1

    In order to divert arbitrage, we must have d<er<u

    Where,

    ● "r is the risk-free rate of interest"

    Let at t = 0, we hold a portfolio that consists of α unit of the underlying asset and β unit of cash. So the value of the portfolio at time t = 0 is:

    V0=αS0+β[2]

    The value of the portfolio at time t = 1 will be

    V1={αS0u+βeriftheunderlyingassetpricegoesupαS0d+βeriftheunderlyingassetpricegoesdown[3]

    Let the derivatives pay cu if the prices of the underlying asset go up and cd if the prices of the underlying asset go down.

    Let us choose α and β so that

    V1={cuiftheunderlyingassetpricegoesupcdiftheunderlyingassetpricegoesdown[4]

    So, from equation (3) and (4) we get

    αS0u+βer=cu[5]
    αS0d+βer=cd[6]

    Subtract equation (5) and (6) we get

    αS0uαS0d=cucd
    αS0(ud)=cucd
    α=cucdS0(ud)

    Put value of α in equation (6) we get

    β=er(cducudud)

    Substitute the value of α and β in equation (2) we get

    V0=er{cu(erdud)+cd(uerud)}
    V0=er{cu(erdud)+cd(1eruud)}
    V0=er(cuq+cd(1q))

    Where Synthetic probability (q) is simple number defined as:

    q=erdud

    If we denote the payoff of the options at time t = 1 by a random variable byct, we can write equation (1) as:

    V0=erEQ(ct)

    Where Q is a probability measure which gives probability q for upward and 1-q for downward and q is simple a real number and it is known as synthetic probability.

    In a CBM we have mean and variance implies by the BM correspond to the mean and variance of the lognormal distribution. For recombining BM, an additional condition that leads to a unique solution is:

    u=1d

    It means that an up step and a down step would mean the stock price, after two steps, is the same as it is at time t = 0

    Here we are using an assumption that the underlying stock follows a Geometric Brownian Motion (GBM). I.e. the underlying stock price changes continuously through time according to the stochastic differential equation (SDE):

    dSt=St(μdt+σdZt)

    Recall that this is same as the lognormal model. Here we showed that if dSt=St(μdt+σdZt), then StS0 has a lognormal distribution with parameters (α0.5σ2)t and σ2t.

    Under the risk-neutral law if we parameterise the lognormal distribution so that:

    ln{St/St0}N[(rσ22)(t0),σ2(tt0)]

    Then the conditions that must be met are:

    E[S(t+δt)/S(t)]=erδt(a)

    And

    var[ln{S(t+δt)/S(t)}]=σ2δt(b)

    Where:

    δt is the time interval of each step in the BM

    S(t) or St denotes the price of the underlying asset at time t

    Noting that:

    E[S(t+δt)/S(t)]=qu+(1q)d

    It follows from the equation (a) that:

    qu+(1q)d=erδt

    Solve the above equation

    q=erδtdud

    Using the equation b and assume that u=1d

    var[ln{S(t+δt)/S(t)}]=q(ln(u))2+(1q)(ln(u))2E[ln{S(t+δt)/S(t)}]2
    =(ln(u))2E[ln{S(t+δt)/S(t)}]2(c)

    The last term involves terms of higher order than δt

    i.e.

    E[ln{S(t+δt))/S(t)}]2=f{(δt)2}(d)

    When δt0 the equation d tends to 0

    So, if we ignore the above equation and equate the equation c to σ2t, we get

    var[ln{S(t+δt)/S(t)}]=(ln(u))2=σ2t(e)

    If we solve the equation e, we get the value of u

    u=eσδt

    And w know than d=1u, so

    d=eσδt

    The Black Scholes Merton model (sometimes also known as Black Scholes model) was developed in the year 1973 by the famous economists namely Fisher Black, Robert Merton and Myron Scholes. Because of this model R. Merton and M. Scholes awarded by the Nobel Prize in 1977. Nowadays it is a famous model for estimating the values of the options. It has some assumptions that are given are -the price of an underlying asset follow a GBM, Arbitrage opportunity does not exist in the market, there are no taxes, the risk-free rate of interest is constant and also the volatility (Hull and White, 1987). Essentially, BSM formula shows us how to find the price of an option contract (call and put option) can be determined by using a simple formula. The formula for European call "c" and put option "p" is:

    c=StN(d1)KerTN(d2)[7]
    p=KerTN(d2)StN(d1)[8]

    Where,

    d1=ln(StK)+(r+σ22)TσT
    d2=ln(StK)+(rσ22)TσT=d1σT

    "St is the present price of the price of an underlying asset, K is the strike price, r is the free risk rate interest, σ is the volatility of the stock and N(*) is the cumulative distribution function for a standard normal distribution, StN(d1) is the present value of the underlying asset if the option is exercised and Ker is the present value of the strike price if the option is exercised (Lazarova et al., 2014)".

    We have taken the historical data of State Bank of India-NSE of one year starting from 1st Jan 2010 to 31st Dec 2010 to estimate the parameters. The risk-free rates of interest, volatility have been estimated from historical data.

    The underlying asset price St is observed only at fixed intervals. When we speak about the historical volatility of an underlying asset price, it means we actually mean historical volatility of return. It is a measure of variation of underlying asset price over a time period (Baxter and Rennie, 1996). Mathematically historical volatility is an approximate same as the standard deviation of the return.

    σ=nt=1(vt¯vt)2t1

    Where, vt=ln(StSt1)

    St is the stock price, ¯St is the mean of the stock price and t is a number of observations.

    After calculation we have

    Standard deviation/Volatility σ=4.1317%

    A non-dividend paying underlying assets has volatility σ=4.1317%. We can easily calculate the value of the parameters u, d and q by using the following formulae

    u=eσδt

    d=eσδt

    q=erδtuddud

    Where σ(volatility) as is computed above, r=0.45513% is risk-free interest or expected return, δt=1 and er = 0.995459041

    After calculation we have

    u=1.042182

    d=0.959525

    q=0.434735

    1q=0.565265

    We can use above information to check our assumption i.e. no arbitrage opportunity exist in the market or not. In order to avoid arbitrage we must have d<er<u.

    The above information satisfies our inequality d<er<u. Hence there is no arbitrage opportunity exist in the market.

    In BM we use lattices to estimate the value of the option (call/put) at an initial time (t = 0) that provides a payoff for each call and put options at different strike prices that we generated by using excel function at a expiry/mature date based on the non-dividend paying shares and without premium. The below table 1 shows the value of a call and put option using the equation (1) at different strikes prices and we used ln(closeprice) of SBI as share price at time t = 0.

    Table 1.  Value of call and put option in BM.
    Date S0 Max(S1K,0) Max(KS1,0) Max(S2K,0) Max(KS2,0) call option/ V0 put option/ V0
    01-01-2010 7.62881 0.00000 3.04939 0.00000 3.67997 0.00000 3.42137
    01-02-2010 7.58797 1.90805 0.00000 1.28085 0.00000 1.56060 0.00000
    02-03-2010 7.63926 0.00000 7.03850 0.00000 7.66994 0.00000 7.42917
    01-04-2010 7.74097 4.06750 0.00000 3.42765 0.00000 3.72272 0.00000
    03-05-2010 7.72701 0.00000 0.94705 0.00000 1.58575 0.00000 1.31405
    01-06-2010 7.74153 4.06809 0.00000 3.42819 0.00000 3.72329 0.00000
    01-07-2010 7.82521 0.15529 0.00000 0.00000 0.49152 0.06782 0.27911
    02-08-2010 7.92530 0.00000 3.74039 0.00000 4.39548 0.00000 4.12944
    01-09-2010 8.08347 3.42445 0.00000 2.75629 0.00000 3.06066 0.00000
    01-10-2010 8.05535 0.39514 0.00000 0.00000 0.27069 0.17257 0.15371
    01-11-2010 8.00376 4.34138 0.00000 3.67981 0.00000 3.98552 0.00000
    01-12-2010 7.94162 0.00000 5.72339 0.00000 6.37982 0.00000 6.12225

     | Show Table
    DownLoad: CSV

    The table 1 simply shows the value of European call and put option at one period using BM.

    In BSM estimate the value of the option (call/put) at an initial time (t = 0) at different strike prices that we generated by using Excel function at an expiry/mature date based on the non-dividend paying shares. The below table 2 shows the value of a call and put option using BSM (equation 7, 8) at different strike prices and we used ln(closeprice) of SBI as share price at time t = 0.

    Table 2.  Value of call and put option in BSM.
    Date S0 N(d1) Equation 7 N(d2) Equation 8 Value of call option (c) N(d1) N(d2) Value of put option p
    01-01-2010 7.62881 1.83E-19 1.26E-19 6.26E-21 1 1 3.421369
    01-02-2010 7.58797 1 1 1.560599 1.11E-08 1.41E-08 5.93E-10
    02-03-2010 7.63926 6.82E-61 3.45E-61 1.3E-62 1 1 7.429168
    01-04-2010 7.74097 1 1 3.722722 3.71E-57 7.16E-57 7.42E-59
    03-05-2010 7.72701 7.82E-05 6.62E-05 5.85E-06 0.999922 0.999934 1.314055
    01-06-2010 7.74153 1 1 3.723287 3.6E-57 6.96E-57 7.21E-59
    01-07-2010 7.82521 0.266254 0.252866 0.051339 0.733746 0.747134 0.262627
    02-08-2010 7.92530 2.03E-24 1.33E-24 6.42E-26 1 1 4.129438
    01-09-2010 8.08347 1 1 3.060659 4.28E-31 6.91E-31 1.23E-32
    01-10-2010 8.05535 0.530838 0.514382 0.142253 0.469162 0.485618 0.123398
    01-11-2010 8.00376 1 1 3.985517 6.75E-63 1.35E-62 1.33E-64
    01-12-2010 7.94162 1.09E-43 6.14E-44 2.56E-45 1 1 6.122248

     | Show Table
    DownLoad: CSV

    The table 2 simply shows the value of European call and put option at one period using BSM.

    Suppose that x1, x2, x3......xn and y1, y2, y3......yn are two independent random samples. We set a null hypothesis

    H0 = there are no differences between the means

    H1 = there are differences between the means of the two samples.

    In this model, it is assumed that the variances between the two sample data are equal but unknown.

    The t-test statistic under H0 is:

    t=ˉxˉyS1n1+1n2tn1+n22

    Where ˉx is the mean of the sample x1, x2, x3... xn and ˉy is the mean of the sample y1, y2, y3... yn.

    S is the standard deviation, n1 + n2 -2 is degree of freedom and n1 and n2 is the total elements in respective groups.

    The hypothesis is:

    H0 = there are no mean difference between BM and BSM at one period

    H1 = There are mean difference between BM and BSM at one period

    Table 3.  The Basic statistics of options using BM and BSM.
    Paired Samples Statistics
    Mean N Std. Deviation Std. Error Mean
    Pair 1 Put option (BM) 1.90409167 12 2.692650519 .777301251
    Put option (BSM) 1.90019192 12 2.695360837 .778083653
    Pair 2 Call option (BM) 1.35776500 12 1.740456698 .502426572
    Call option (BSM) 1.35386515 12 1.743466636 .503295466

     | Show Table
    DownLoad: CSV

    The standard deviation or variance is almost equal in call and put options using BM and BSM.

    Table 4.  The correlation between the options using BM and BSM.
    Paired Samples Correlations
    N Correlation Sig.
    Pair 1 Put option (BM) & Put option (BSM) 12 1.000 .000
    Pair 2 Call option (BM) & Call option (BSM) 12 1.000 .000

     | Show Table
    DownLoad: CSV
    Table 5.  t-Test -Paired Two Sample for Means at 95% confidence interval between the call options.
    Call option (BM) Call option (BSM)
    Mean 1.357765 1.353865167
    Variance 3.029189518 3.039675874
    Observations 12 12
    Pooled Variance 3.034432696
    Hypothesized Mean Difference 0
    Df 22
    t Stat 0.005483817
    P(T < =t) one-tail 0.497836996
    t Critical one-tail 1.717144335
    P(T < =t) two-tail 0.995673993
    t Critical two-tail 2.073873058

     | Show Table
    DownLoad: CSV
    Table 6.  t-Test -Paired Two Sample for Means at 95% confidence interval between the put options.
    Put option (BM) Put option (BSM)
    Mean 1.904091667 1.900191917
    Variance 7.250366818 7.264970044
    Observations 12 12
    Pooled Variance 7.257668431
    Hypothesized Mean Difference 0
    Df 22
    t Stat 0.003545797
    P(T < =t) one-tail 0.498601413
    t Critical one-tail 1.717144335
    P(T < =t) two-tail 0.997202826
    t Critical two-tail 2.073873058

     | Show Table
    DownLoad: CSV

    The correlation between the call option using BM and BSM is 1 and in case of the put option, it is also 1. It indicates that it is perfect correlated.

    The 95%, confidence interval of the difference provides an estimate of the boundaries between which the true mean difference lies in 95% of all the possible random samples of 12 values. At two-tailed, we are accepting the null hypothesis because t-stat value is lying between -2.0738 and 2.0738 at 95% confidence interval in case of call option and 2.0738 and 2.0738 in case of put option and we may conclude that the mean of option values using BM and BSM at one period do not differ significantly. We can conclude that we don't have enough evidence to reject the null hypothesis.

    It is used in ANOVA, in order to create the confidence interval (C.I) of all pair-wise differences between factor level means. It applies simultaneously to the set of all pair-wise comparisons.

    The authors used the Turkey pair-wise comparisons in order to check whether there is any difference between a pair of groups or not. In other words, the pair of groups are statistically significant or not (the call option value using BM and BSM, put option using BM and BSM).

    In these results, the table 7 shows that group contains letter A only in all factors. That letter A in all the groups indicates that there are no significant differences between the factors i.e. call option using BM and BSM, put option value using BM and BSM.

    Table 7.  Grouping Information Using the Tukey Method and 95% Confidence.
    Factor N Mean Grouping
    call option value (BM) 12 1.358 A
    call option value (BSM) 12 1.354 A
    put option value (BM) 12 1.904 A
    put option value (BSM) 12 1.900 A
    Note: Means that do not share a letter are significantly different.

     | Show Table
    DownLoad: CSV

    Table 8, shows the Tukey simultaneous C.I. The pairs, the call option value using BM and BSM and the put option value using BM and BSM have C.I. in between positive and negative which means that all the intervals include zero. The confidence interval for the means of call option (BM) -call option (BSM) and put option value (BM) – put option value (BSM) are (-1.479, 1.471) and (-2.28, 2.28) respectively. Both the pairs include zeros. If zero is not included in intervals then the means are statistically significant. But in our result both the pairs include zeros that indicate there is no significant difference between the means of the groups. The figure (4) and (5) also shows the bolded dot, which also indicates that the intervals include zero. The individual confidence level is 95%, it means that you can be 95% confident that each individual interval contains the true difference between the specific pair of groups.

    Table 8.  Tukey Simultaneous Tests for Differences of Means.
    Difference of Levels Difference of Means SE of Difference 95% CI T-Value Adjusted P-Value
    call option value (BM) -call option value (BSM) -0.004 0.711 (-1.479, 1.471) -0.01 0.996
    put option value (BM) -put option value (BSM) -0.00 1.10 (-2.28, 2.28) -0.00 0.997
    Note: Individual confidence level = 95.00%.

     | Show Table
    DownLoad: CSV
    Figure 4.  Tukey simultaneous 95% C.Is for call option using BM and BSM.
    Figure 5.  Tukey simultaneous 95% C.Is for put option using BM and BSM.

    The BM and BSM is an important model in a financial market in order to calculate the option pricing. In this paper, we studied how two models (BM and BSM) are approximately equal to each other at one period. The procedures we did in this paper in order to understand the difference between the two option pricing models at one period. The correlation between the call and put option using both the pricing methods at one period is perfect positive correlated. It indicated that in both the models the price of the call and put option either goes up or down. The statistical t-test showed that there are no significant differences between the call and put option using BM and BSM at 95% confidence interval. In case of Tukey method, all the factors contain only same letter 'A' which indicates there is no significant difference between the BSM and BM. Also at 95% C.I., the intervals are involving the zero also. As per Tukey method, if zero value includes in the C.I. then there is no evidence to reject the null hypothesis. In other words, we can say that there is no statistically significant difference between the BM and BSM.

    This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

    All authors declare no conflicts of interest in this paper.



    [1] Das SR, Mishra DP, Ray PK, et al. (2021) Power quality improvement using fuzzy logic-based compensation in a hybrid power system, Int J Power Electron Drive Syst 12: 576–584. https://doi.org/10.11591/ijpeds.v12.i1.pp576-584 doi: 10.11591/ijpeds.v12.i1.pp576-584
    [2] Das SR, Mishra AK, Ray PK, et al. (2022) Application of artificial intelligent techniques for power quality improvement in hybrid microgrid system. Electronics 11: 3826: 1–19. https://doi.org/10.3390/electronics11223826 doi: 10.3390/electronics11223826
    [3] Das B, Panigrahi PK, Das SR, et al. (2021) Power quality improvement in a photovoltaic based microgrid integrated network using multilevel inverter. Int J Emerg Electr Power Syst 23: 197–209. https://doi.org/10.1515/ijeeps-2021-0040 doi: 10.1515/ijeeps-2021-0040
    [4] Sinvula R, Abo-Al-Ez KM, Kahn MT (2019) Harmonic source detection methods: A systematic literature review. IEEE Access 7: 74283–74299. https://doi.org/10.1109/ACCESS.2019.2921149 doi: 10.1109/ACCESS.2019.2921149
    [5] Saxena D, Bhaumik S, Singh SN (2014) Identification of multiple harmonic sources in power system using optimally placed voltage measurement devices. IEEE Trans Ind Electron 61: 2483–2492. https://doi.org/10.1109/TIE.2013.2270218 doi: 10.1109/TIE.2013.2270218
    [6] Carta D, Muscas C, Pegoraro PA, et al. (2019) Identification and estimation of harmonic sources based on compressive sensing. IEEE Trans Instrum Meas 68: 95–104. https://doi.org/10.1109/TIM.2018.2838738 doi: 10.1109/TIM.2018.2838738
    [7] Li C, Xu W, Tayjasanant T (2004) A "critical impedance"—based method for identifying harmonic sources. IEEE Trans Power Deliv 19: 671–678. https://doi.org/10.1109/TPWRD.2004.825302 doi: 10.1109/TPWRD.2004.825302
    [8] Stevanović D, Petković P (2015) A single-point method for identification sources of harmonic pollution applicable to standard power meters. Electr Eng 97: 165–174. https://doi.org/10.1007/s00202-014-0324-z doi: 10.1007/s00202-014-0324-z
    [9] Li C-S, Bai Z-X, Xiao X-Y, et al. (2016) Research of harmonic distortion power for harmonic source detection. International Conference on Harmonics and Quality of Power. IEEE, Belo Horizonte, Brazil, 1–4. https://doi.org/10.1109/ICHQP.2016.7783437
    [10] Anu G, Fernandez FM (2020) Identification of harmonic injection and distortion power at customer location. 19th International Conference on Harmonics and Quality of Power (ICHQP), Dubai, United Arab Emirates, 1–5. https://doi.org/10.1109/ICHQP46026.2020.9177869
    [11] Cataliotti A, Cosentino V, Nuccio S (2008) Comparison of nonactive powers for the detection of dominant harmonic sources in power systems. IEEE Trans Instrum Meas 57: 1554–1561. https://doi.org/10.1109/TIM.2008.925338 doi: 10.1109/TIM.2008.925338
    [12] Xu W, Liu X, Liu Y (2003) An investigation on the validity of power direction method for harmonic source determination. IEEE Trans Power Deliv 18: 214–219. https://doi.org/10.1109/TPWRD.2002.803842 doi: 10.1109/TPWRD.2002.803842
    [13] IEEE Std 1459-2010: IEEE Standard definitions for the measurement of electric power quantities under sinusoidal, nonsinusoidal, balanced, or unbalanced conditions, 2010.
    [14] Menti A, Zacharias T, Milias-Argitis J (2007) Geometric algebra: a powerful tool for representing power under nonsinusoidal conditions. IEEE Trans Circuits Syst I, Fundam Theory Appl 54: 601–609. https://doi.org/10.1109/TCSI.2006.887608 doi: 10.1109/TCSI.2006.887608
    [15] Menti A, Zacharias Th, Milias-Argitis J (2010) Power components under nonsinusoidal conditions using a power multivector. Proc 10th Conference-Seminar International School on Nonsinusoidal Currents and Compensation (ISNCC). Lagow, Poland, 174–179. https://doi.org/10.1109/ISNCC.2010.5524495
    [16] Doran C, Lasenby A (2003) Geometric Algebra for Physicists. Cambridge University Press. https://doi.org/10.1017/CBO9780511807497
    [17] Bollen M, Olofsson M, Larsson A, et al. (2014) Standards for supraharmonics (2 to 150 kHz). IEEE Electromag Compatibil Mag 3: 114–119. https://doi.org/10.1109/MEMC.2014.6798813 doi: 10.1109/MEMC.2014.6798813
    [18] Barkas D, Ioannidis G, Kaminaris S, et al. (2022) Design of a supraharmonic monitoring system based on an FPGA device. Sensors 22: 1–20. https://doi.org/10.3390/s22052027 doi: 10.3390/s22052027
    [19] Menti A, Barkas D, Kaminaris S, et al. (2012) Supraharmonic emission from a three-phase PV system connected to the LV grid. Energy Rep 7: 527–542. https://doi.org/10.1016/j.egyr.2021.07.100 doi: 10.1016/j.egyr.2021.07.100
  • This article has been cited by:

    1. Yue Liu, Dehao Ruan, Stability of a Class of Impulsive Neutral Stochastic Functional Partial Differential Equations, 2020, 2020, 1026-0226, 1, 10.1155/2020/9051396
    2. Dehao Ruan, Yue Liu, Generalized Halanay Inequalities with Applications to Generalized Exponential Stability and Boundedness of Time-Delay Systems, 2019, 2019, 1024-123X, 1, 10.1155/2019/6072481
    3. Amir Ahmad Dar, N. Anuradha, Ziadi Nihel, 2021, chapter 8, 9781799830535, 149, 10.4018/978-1-7998-3053-5.ch008
    4. Ramishka Nardhamuni, Kayleigh Greenslade, Wayne van Zijl, An analysis of the use of share-based payments by the JSE Top 100 companies, 2022, 1029-1954, 1, 10.1080/10291954.2021.2008610
    5. Sojin Kim, Jimin Kim, Jongwoo Song, Option pricing and profitability: A comprehensive examination of machine learning, Black-Scholes, and Monte Carlo method, 2024, 31, 2287-7843, 585, 10.29220/CSAM.2024.31.5.585
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1587) PDF downloads(87) Cited by(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog