• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제19권1호
• /
• pp.1-21
• /
• 2015

This paper presents two algorithms based on the Jamshidian equation which is from the Black-Scholes partial differential equation. The first algorithm is for American call options and the second one is for American put options. They compute numerically free boundary and then option price, iteratively, because the free boundary and the option price are coupled implicitly. By the upwind finite-difference scheme, we discretize the Jamshidian equation with respect to asset variable s and set up a linear system whose solution is an approximation to the option value. Using the property that the coefficient matrix of this linear system is an M-matrix, we prove several theorems in order to formulate a bisection method, which generates a sequence of intervals converging to the fixed interval containing the free boundary value with error bound h. These algorithms have the accuracy of O(k + h), where k and h are step sizes of variables t and s, respectively. We prove that they are unconditionally stable. We applied our algorithms for a series of numerical experiments and compared them with other algorithms. Our algorithms are efficient and applicable to options with such constraints as r > d, $r{\leq}d$, long-time or short-time maturity T.

• 대한조선학회지
• /
• 제10권1호
• /
• pp.27-32
• /
• 1973

This paper shows the method for obtaining the body flutter velocity and frequency for flight body which consists of low aspect ratio wing and body combination by assuming slender body of cylinderical shell structure. The stiffness matrix of the cylinderical shell is represented from Donnel eq. by the finite difference method, and also unsteady aerodynamic influence matrix is represented by the Doublet Lattice Method of Albano & Rodden. The flutter matrix can be obtained from those matrices.

• 전기학회논문지P
• /
• 제58권1호
• /
• pp.85-89
• /
• 2009

We investigated the quench characteristics in accordance with increase of turns number of trigger coil and shunt resistance of matrix-type superconducting fault current limiter (SFCL) with $2{\times}3$ array. The matrix-type SFCL consists of the trigger part to apply magnetic field and the current-limiting part to limit fault current. The fault current limiting characteristics according to the increase of magnetic field and applied voltage were nearly same. This is because the application of magnetic field hasn't an affect on total impedance of the SFCL. When turns number of a reactor increased, the voltage difference between two superconducting units in the current-limiting part according was decreased. The resistance difference generated in two superconducting units was also decreased. Therefore, we confirmed that the differences of the critical behaviors between superconducting units were reduced by application of magnetic field. By this results, we could decide the optimum turns number of reactor to apply magnetic field.

• Journal of Electrical Engineering and Technology
• /
• 제9권6호
• /
• pp.2049-2057
• /
• 2014

In this paper, a numerical method for analyzing coupling between high-altitude electromagnetic pulse (HEMP) as external field and a shielded twisted pair (STP) cable is proposed, which is based on an expanded chain matrix. Load responses of electromagnetic (EM) field excitation in uniform transmission line (TL) are solved by Baum-Liu-Tesche (BLT) equations in frequency domain, however, it is difficult to apply BLT equations to solve load responses of STP cable because the iteratively changing configuration of each twisted pairs are involved in cable. To avoid this problem and decrease memory and CPU time, we proposed the expanded chain matrix modeling method that is calculated using ABCD parameters, and applied multi-conductor transmission line (MTL) theory to consider the EMP coupling effectiveness of each twisted pairs. The results implemented by the proposed method are presented and compared with those obtained by the finite-difference time domain (FDTD) method as a kind of 3D full wave analysis.

• Journal of the Optical Society of Korea
• /
• 제6권3호
• /
• pp.121-127
• /
• 2002

The operating principal of a ring laser gyroscope depends on the phase difference for the counter-propagating waves within a closed path. The reflecting mirrors mounted on the monoblock form the traveling waves. The manufacturing accuracy of the monoblock influences the traveling path of ray, the sensitivity of laser resonator for misalignments, and diffraction losses. A 3 $\times$ 3 ray transfer matrix was derived for optical components with centering and squaring errors in a ring resonator. The matrix can be utilized to predict the optical ray paths on the basis of the manufacturing errors of the monoblock as well as the misalignment of mirrors. Then the distance and orientation (o. slope) at the arbitrary plane inside the resonator along the ideal optical path can be calculated from the chain multiplication of the ray transfer matrix for each optical component in one round trip. We also show that the counter-propagating rays In a ring resonator with errors does not coincide in each round trip, which results in gain difference between two beams, and how these errors can be adjusted through the alignment procedure. Finally this 3 $\times$ 3 ray matrix formalism can be used to calculate the beam size and its displacement from the optical axis and the deviation at the diaphragm.

• Kyungpook Mathematical Journal
• /
• 제55권4호
• /
• pp.909-931
• /
• 2015

This paper presents some new paranormed sequence spaces $X(r,s,t,p;{\Delta})$ where $X{\in}\{l_{\infty}(p),c(p),c_0(p),l(p)\}$ defined by using generalized means and difference operator. It is shown that these are complete linear metric spaces under suitable paranorms. Furthermore, the ${\alpha}$-, ${\beta}$-, ${\gamma}$-duals of these sequence spaces are computed and also obtained necessary and sufficient conditions for some matrix transformations from $X(r,s,t,p;{\Delta})$ to X. Finally, it is proved that the sequence space $l(r,s,t,p;{\Delta})$ is rotund when $p_n$ > 1 for all n and has the Kadec-Klee property.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 추계 학술대회논문집
• /
• pp.48-58
• /
• 1999

The Characteristic Flux Difference Splitting (CFDS) scheme designed to adapt the characteristic boundary conditions at the wall and inflow/outflow boundary planes satisfies Roe's property U, although the CFDS Jacobian matrix is decomposed by a product of elaborate transformation matrices and explicit eigenvalue matrix. When the CFDS algorithm, thus a variant of Roe's scheme, is applied straightforwardly to hypersonic flows over a blunt body, the strong bow shock gradually breaks down near the stagnation point. This numerical instability is widely observed by many researchers employing flux-difference method, known in the literature as the carbuncle phenomenon. Many remedies have been proposed and resulted in partial cures. When the idea of Sanders et al. which identifies the minimum eigenvalues near the discontinuity present is applied to CFDS method, it is shown that the instability problem can be controlled successfully. A few flux splitting methods have also been tested and results are compared against the Nakamori's Mach 8 blunt body flow.

• ETRI Journal
• /
• 제45권4호
• /
• pp.704-712
• /
• 2023

This paper presents a flexible finite-difference technique for analyzing the nonuniform guiding structures. Because the voltage and current variations along the nonuniform structure differ for each segment, this work considers the adaptable discretization steps. This technique increases the accuracy of the final response. Moreover, by applying the singular value decomposition and discarding the nonprincipal singular values, an optimal lower rank approximation of the discretization matrix is obtained. The computational cost of the introduced method is significantly reduced using the optimal discretization matrix. Also, the proposed method can be extended to the nonuniform waveguides. The technique is verified by analyzing several practical transmission lines and waveguides with nonuniform profiles.

• 한국수자원학회논문집
• /
• 제42권4호
• /
• pp.331-341
• /
• 2009

본 연구에서는 generalized logistic(GL) 분포의 최우도 추정량(maximum likelihood estimate)에 대한 불확실성 추정을 위하여 사용되는 관측정보행렬(observed information matrix)과 Fisher 정보행렬(Fisher information matrix)의 정확도를 비교해 보고자 하였다. 타 분포형에 대한 기존의 연구결과에서 표본의 크기가 클 경우 매개변수 추정시 관측정보행렬이 동시에 추정되어 계산시간도 단축되고 Fisher 정보행렬의 정확도와도 차이도 거의 없어 관측정보행렬의 사용이 추천된 바 있으나, 최근 사용이 증가되고 있는 GL 분포에 대한 연구결과는 아직 전무한 실정이며 기존 연구문헌의 결과를 토대로 구체적인 연구 없이 관측정보행렬을 사용하고 있는 상황이다. 따라서 본 연구에서는 이를 위해 모의실험을 수행하였으며, 모의 결과 최우도법에 의한 매개변수의 분산 및 공분산은 기존의 연구 결과와 비슷한 결과를 보이나, quantile에 대한 불확실성 추정에는 관측정보행렬보다 Fisher 정보행렬의 사용이 더 적절할 것으로 판단되었다.

• 전산구조공학
• /
• 제10권1호
• /
• pp.193-200
• /
• 1997

탄성지반상의 원형탱크해석에는 여러방법이 있지만 최근에 널리 사용되는 방법은 유한요소법이다. 그러나 이 방법은 탄성지반상의 탱크해석시 많은 절점수가 필요하게 된다. 이것은 곧 많은 계산기 기억용량 및 계산시간 뿐만 아니라 노력이 필요하게 된다. 본 연구에서는 유사탄성지반보(Analogy of Beam on Elastic Foundation) 및 지반강성행렬(Foundation Stiffness Matrix)을 이용하여 축대칭하중을 받는 축대칭탱크를 뼈대 구조화 할 수 있었다. 또한 이 뼈대 구조를 유한요소로 분할하고, 각 요소 강성행렬(Stiffness Matrix)을 전달행렬(Transfer Matrix)로 전환하여 전달행렬법으로 원형탱크를 해석 할 수 있었다. 유한요소법과 전달행렬법을 탄성지반상의 원형탱크 해석에 적용한 결과 두 해석결과의 차이는 없고, 전달행렬법을 적용한 경우 최종 연립방정식수가 4개로 간략화 되었다.