分数Black-Scholes模型下美式期权定价的积分方程式

为得到分数Black-Scholes模型下美式期权价格的公式,文章以看涨期权为例,应用偏微分方程法,推导期权价格的积分方程式.由于美式期权的价格可分解为欧式期权的价格和由于提前实施需要增付的期权金,而提前实施期权金与最佳实施边界的位置有关,所以为导出最佳实施边界所满足的方程,文章首先研究分数Black-Scholes方程的基本解,然后建立美式看涨期权的分解公式,推导最佳实施边界适合的非线性积分方程,从而得到美式看涨期权价格的积分方程式.美式看跌期权价格的积分方程式类似得到.     

随机波动率与双指数跳扩散组合模型的美式期权定价

在股价满足Cox-Ingersoll-Ross(CIR)随机波动率与Kou的双指数跳扩散组合模型下,利用随机分析方法讨论了美式看跌期权函数及最佳实施边界的性质.应用一阶线性近似实施边界获得了期权价格的拟解析式和实施边界满足的非线性方程.进一步,应用梯形法离散处理方程式内积分表达式,建立了期权最佳实施边界和价格的数值算法.最后分别给出了常数波动率或CIR随机波动率的数值实例.     

非线性Black-Scholes模型下障碍期权定

研究了原生资产价格遵循非线性Black-Scholes模型时障碍期权的定价问题.首先,根据混合分数布朗运动的Ito公式和金融市场的复制策略,得到了障碍期权适合的抛物初边值问题.其次,利用扰动理论中单参数摄动展开方法,给出了障碍期权的近似定价公式.最后,利用Feyman-Kac公式分析了近似定价公式的误差估计问题,结果表明近似解一致收敛于相应期权价格的精确解.     

跳扩散模型下美式期权的对称性

利用分析方法得到了跳扩散模型下美式看涨、看跌期权的价格和最佳实施边界间的对称性公式．美式看涨和看跌期权价格问的对称关系通常是利用概率理论得到,这里给出了这些结果在跳扩散模型下的另一种证明．此外,由本文所得结果和偏微分方程理论,可以得到跳扩散模型下美式看涨期权的最佳实施边界以及永久美式期权的若干性质．     

美式看跌期权的两点Geske-Johnson近似定价法

在分数Black-Scholes模型下,应用两点Geske-Johnson定价法推导连续支付红利为常数的美式看跌期权的近似公式.首先假定期权没有提前实施,其价格为对应欧式看跌期权的价格;再将期权的实施时刻指定为两个时刻,通过中性风险定价法推导价格公式,然后利用两点Geske-Johnson定价法得到美式看跌期权价格的近似公式.最后给出一个数值算例,结果显示Hurst参数和到期日对价格的影响.     

跳跃扩散型离散算术平均亚式期权的近似价格公式

在标的资产价格遵循跳跃扩散过程条件下 ,研究没有封闭形式解的离散算术平均亚式期权 ,运用二阶 Edgeworth逼近得到离散算术平均亚式期权的近似价格公式 .     

永久经理期权的最优实施策略

利用偏微分方程数值方法研究金融市场上永久经理期权(ESOs)的最优实施策略问题.讨论了两种实施情况,即通常的整体实施情况以及非限制实施情况.在非限制实施情况下,持有者在任意可实施时刻可以实施其持有的任意份ESOs.两种实施情况下的最优实施策略分别对应着一个抛物型变分不等式定解问题的自由边界(最佳实施边界).通过数值分析的方法分别研究了自由边界的性质,比较了两种情况下自由边界的异同及其所对应的金融意义.     

跳跃-扩散过程下百慕大交换期权定价

在等价鞅测度下,利用条件期望等知识导出在风险中性定价模型中,标的资产服从跳跃-扩散过程时百慕大交换期权的解析定价公式,依此结合Richardson两点外推加速法得到美式交换期权近似解.提出的数值算例阐明提前执行特征具有重要经济价值.定价结果可以评估场外交易的金融期权价格尤其是实物期权定价.     

具有分数O-U过程的永久美式看跌期权的定价

该文研究具有分数Ornstein-Uhlenbeck过程的永久美式看跌期权的定价问题.首先, 利用分析金融衍生品定价的delta对冲方法和无套利原理, 遵循标准的讨论步骤, 建立了标的资产价格服从分数Ornstein-Uhlenbeck过程的欧式看涨期权和看跌期权的定价公式.然后, 通过求解一个自由边界问题, 对标的资产价格服从分数Ornstein-Uhlenbeck过程的永久美式看跌期权的定价以及实施该期权时的临界标的资产价格给出了显式解.     

求解美式回望期权的有限元方法

期权作为一种金融衍生产品,在欧美国家一直很受欢迎.由于其规避风险的特性,期权也吸引了中国投资者的兴趣.基于市场的需求,2015年初,上海证券交易所推出了中国首批期权产品,期权定价问题的研究热潮正席卷全球.本文研究的美式回望期权,是一种路径相关的期权,其支付函数不仅依赖于标的资产的现值,也依赖其历史最值.分析回望期权的特点,不难发现:1)这类期权空间变量的变化范围为二维无界不规则区域,难以应用数值方法直接求解;2)最佳实施边界未知,使得该问题变得高度非线性.本文的主要工作就是解决这两个困难,得到回望期权和最佳实施边界的数值逼近结果.现有的处理问题1)的有效方法是采用标准变量替换、计价单位变换以及Landau变换将定价模型化为一个[0,1]区间上的非线性抛物问题,本文也将沿用这些技巧处理问题1).进一步,采用有限元方法离散简化后的定价模型,并论证了数值解的非负性,提出了利用Newton法求解离散化的非线性系统.最后,通过数值模拟,验证了本文所提算法的高效性和准确性.     

分数次布朗运动的欧式障碍期权定价

假定标的股票服从分数次布朗运动,应用偏微分方程的方法求出下降敲出欧式看涨障碍期权价格显示解,以及看涨-看跌的平价关系式.最后,通过有限差分法比较了显示解的准确性,分析了Hurst参数对期权价格和风险特征参数的影响.     

A numerical technique for solving fractional variational problems by Müntz–Legendre polynomials

This paper presents an efficient approximate method for solving a class of fractional variational problems (FVPs). The fractional derivatives are described in the Caputo sense. In the our method, we use the Müntz–Legendre orthonormal basis and the properties of Rayleigh–Ritz method for fractional calculus to reduce FVPs to solve a system of algebraic equations which solved using Newton’s iterative method. Also, an estimation of the error is given in the sense of Sobolev norms. The illustrative examples are provided to demonstrate the applicability and simplicity of the new technique.     

分数布朗运动驱动的幂型期权定价模型研究

股价运动分形特征的发现,说明布朗运动作为期权定价模型的初始假定存在缺陷.本文假定标的资产价格服从几何分数布朗运动,利用分数风险中性测度下的拟鞅(quasi-martingale)定价方法重新求解分数Black-Scholes模型,进而对幂型期权进行定价.结果表明,幂型期权结果包含了Black-Scholes公式和平方期权结果,且相比标准期权价格,分数期权价格要同时取决于到期日和Hurst参数H.     

Inverse problem for a space‐time fractional diffusion equation: Application of fractional Sturm‐Liouville operator

An inverse problem of determining a time‐dependent source term from the total energy measurement of the system (the over‐specified condition) for a space‐time fractional diffusion equation is considered. The space‐time fractional diffusion equation is obtained from classical diffusion equation by replacing time derivative with fractional‐order time derivative and Sturm‐Liouville operator by fractional‐order Sturm‐Liouville operator. The existence and uniqueness results are proved by using eigenfunction expansion method. Several special cases are discussed, and particular examples are provided.     

一类时间分数阶偏微分方程的解

考虑一类时间分数阶偏微分方程,该方程包含几种特殊情况:时间分数阶扩散方程、时间分数阶反应-扩散方程、时间分数阶对流-扩散方程以及它们各自相对应的整数阶偏微分方程． 通过Laplace-Fourier变换及其逆变换,该方程在空间全平面和半平面内的基本解可以求出,但其表达式则是通过适当的变形来求．另外,对于有限域上的初边值问题,则可由Sine(Cosine)-Laplace变换导出该方程的一种级数形式的解,并通过两个数值例子来说明该方法的有效性．     

A space-time Legendre spectral tau method for the two-sided space-time Caputo fractional diffusion-wave equation

The space-time fractional diffusion-wave equation (FDWE) is a generalization of classical diffusion and wave equations which is used in modeling practical phenomena of diffusion and wave in fluid flow, oil strata and others. This paper reports an accurate spectral tau method for solving the two-sided space and time Caputo FDWE with various types of nonhomogeneous boundary conditions. The proposed method is based on shifted Legendre tau (SLT) procedure in conjunction with the shifted Legendre operational matrices of Riemann-Liouville fractional integral, left-sided and right-sided fractional derivatives. We focus primarily on implementing this algorithm in both temporal and spatial discretizations. In addition, convergence analysis is provided theoretically for the Dirichlet boundary conditions, along with graphical analysis for several special cases using other conditions. These suggest that the Legendre Tau method converges exponentially provided that the data in the given FDWE are smooth. Finally, several numerical examples are given to demonstrate the high accuracy of the proposed method.     

Vieta–Fibonacci wavelets: Application in solving fractional pantograph equations

In this paper, the Vieta–Fibonacci wavelets as a new family of orthonormal wavelets are generated. An operational matrix concerning fractional integration of these wavelets is extracted. A numerical scheme is established based on these wavelets and their fractional integral matrix together with the collocation technique to solve fractional pantograph equations. The presented method reduces solving the problem under study into solving a system of algebraic equations. Several examples are provided to show the accuracy of the method.     

Tikhonov regularization for infeasible absolute value equations

We investigate the optimum correction of an absolute value equation by minimally changing the coefficient matrix and right-hand side vector using Tikhonov regularization. Solving this problem is equivalent to minimizing the sum of fractional quadratic and quadratic functions. The primary difficulty with this problem is its nonconvexity. Nonetheless, we show that a global optimal solution to this problem can be found by solving an equation on a closed interval using the subgradient method. Some examples are provided to illustrate the efficiency and validity of the proposed method.     

Iterative methods for solving fourth- and sixth-order time-fractional Cahn-Hillard equation

This paper presents analytical-approximate solutions of the time-fractional Cahn-Hilliard (TFCH) equations of fourth and sixth order using the new iterative method (NIM) and q-homotopy analysis method (q-HAM). We obtained convergent series solutions using these two iterative methods. The simplicity and accuracy of these methods in solving strongly nonlinear fractional differential equations is displayed through the examples provided. In the case where exact solution exists, error estimates are also investigated.     

A novel variable‐order fractional nonlinear Klein Gordon model: A numerical approach

In this article, a novel variable order fractional nonlinear Klein Gordon model is presented where the variable‐order fractional derivative is defined in the Caputo sense. The merit of nonstandard numerical techniques is extended here and we present the weighted average nonstandard finite difference method to study numerically the proposed model. Special attention is paid to study the convergence and to the stability analysis of the numerical technique. Moreover, the truncation error is analyzed. Three test examples are provided. Comparative studies are done between the used numerical technique and the weighted average finite difference method. It is found that the stability regions are larger by using the weighted average nonstandard finite difference method.