偏微分方程和数值方法

第11章偏微分方程和数值方法

本章的学习目标

了解布莱克-休尔斯方程所属的二阶偏微分方程的类型；

了解热传导方程的推导过程和它代表的物理意义；

理解偏微分方程所附带的边界条件和初始条件的具体形式和现实意义；

熟悉傅里叶变换方法及其主要性质；

熟悉用傅里叶积分变换来求热传导方程；

掌握通过变量代换来求解布莱克-休尔斯方程；

了解求解偏微分方程的主要几种有限差分格式以及各自的优缺点；

掌握柯尔莫格罗夫方程的推导过程，并理解扩散过程和数学期望之间的联系；

掌握费曼-卡茨定理并了解它同风险中性定价之间的关系；

了解产生随机数、随机分布和随机过程的技术方法；

掌握使用蒙特卡罗模拟技术计算期权的实际操作方法和两种算法优化技术。

完整学习有着数百年知识积淀的偏微分方程理论本身是一项艰巨而耗费时日的工作。幸运的是：在我们所关注的微观金融领域，几乎所有被用到的偏微分方程都只属于其中一个很小的类别——二阶线性偏微分方程。因此这一章的设计思想（用软件行业的术语来说）是完全面向任务的，目标很明确：求解布莱克-休尔斯方程。

我们这样安排本章的结构：首先了解一下偏微分方程的数学表达形式，并在直观的热物理背景下，讨论偏微分方程的定解问题。然后使用经典的傅里叶变换（Fourier transform）技术来求出一般热传导方程的解析解，这就使我们可以遵循布莱克-休尔斯（1973）的方法求解布莱克-休尔斯偏微分方程获得期权价格。但是由于获得解析解的机会并不多，因此我

第11章 偏微分方程和数值方法

·533·

们要学习偏微分方程的数值解法（numerical method ）——有限差分方法（finite difference method ）。此外，根据风险中性定价原理和费曼-卡茨（Feynman-Kac ）理论，布莱克-休尔斯方程还有一种概率解。但是获得这种概率解同样需要一种被称为蒙特卡罗的数值方法，因此本章最后我们要考察这日益重要的种计算机模拟（simulation ）技术。

11.1 介 绍

11.1.1 基本概念

我们在第4、9、10章中都见到过著名的布莱克-休尔斯方程：

0212

2 22=−∂∂+∂∂+∂∂rf S f rS S f S t f σ 其中S 是某种基础产品的价格，r 是无风险收益率，t 是时间，σ是S 的波动率，f 则是基于S 的衍生产品的价格。这个偏微分方程包含了衍生产品价格运动信息，是对所有基础产品和基于它的衍生产品之间相对价格运动形式的高度概括。它在金融理论中的重要性，怎样强调都不过分。

就其本身而言，这是一个有着两个自变量的偏微分方程，这种偏微分方程的更一般形式为：

G Ff Ef Df Cf Bf Af t S tt St SS =+++++ （11-1）

该方程有这样几个特征：

（1）它（被称为）是二阶的，因为（如果）它的最高阶是22 S

f ∂∂； （2）它（被称为）是齐次的，因为（如果）0),(=t S G ；

（3）它（被称为）是线性的，因为（如果）A 、B 、C 、D 这些系数只是自变量S 和t 的函数。

二阶线性偏微分方程有不少类别，数学上对二阶线性偏微分方程的最重要的分类方式受到解析几何中对二次曲线：

022=+++++F Ey Dx Cy Bxy Ax

的分类方法的启发。

我们知道根据上面曲线方程的系数可以判断二次曲线的形状究竟是双曲线、抛物线或者椭圆。类似的在这里，方程（11-1）中：

（1）如果042>−AC B ，就称之为双曲型（hyperbolic ）偏微分方程；

（2）如果042=−AC B ，则称之为抛物型（parabolic ）偏微分方程；

微观金融学及其数学基础

·534·

（3）如果042<−AC B ，则称之为椭圆型（elliptic ）偏微分方程。

对于布莱克-休尔斯方程，因为有：

222/1S A σ=

0===G C B

rS D =

r F E −==,1

而且

042=−AC B

所以现在可以知道，我们面对的是一个二阶齐次线性抛物型偏微分方程。听上去相当复杂，我们接下来的任务就是一步一步揭开这种方程所蕴涵的实际物理意义，并得到求解它的一般方法。

11.1.2 物理意义

一个二阶齐次线性抛物型偏微分方程究竟描述了怎样的一种现象呢？鉴于它最初的来源和最广泛的应用都发生在热物理领域，我们不妨来看一下它的实验背景，这对于我们理解这个方程将提供足够的洞察力。

根据日常生活的经验，我们知道当物体内部各处的温度不一致时，热量就会从高温处向低温处传递，这被称为“热传导”现象。现在假定存在一导热物体，它在3维空间占据

的区域为G ，

边界面为G ∂，我们怎样才能知道它其中的某一部分的温度变化情况呢？用温度函数),,,(t z y x u 表示该物体在t 时刻和),,(z y x 位置的温度，我们来建立该温度函数需要满足的关系式①。

设想从物体G 内任意割取一个由光滑曲面L 所围成的区域D （见图11-1）

。根据热量守恒定律，D 内各点的温度由任一时刻1t 的),,,(1t z y x u 改变为2t 时刻的),,,(2t z y x u 所吸收（或释放）的热量Q ，应当等于从1t 到2t 时间内通过L 进入（或流出）D 内的热量1Q 和D 内热源提供的热量2Q 的总和。

图11-1 热量在物体内部的传导

① 以下分析均假定u 对x ，y ，x 具有二阶连续偏导数，对t 具有一阶连续偏导数。

第11章 偏微分方程和数值方法

·535·

下面我们分别决定这些热量，首先是:

（1）D 内温度改变所需要的热量Q 。假定物体的比热（使单位质量的物体温度改变1摄氏度所需要的热量）为),,(z y x c ，密度为),,(z y x ρ。那么根据物理中的实验规律，无穷小体积

z y x V d d d d =

的温度由),,,(1t z y x u 升高到),,,(2t z y x u 所需要的热量Q d 为：

V t z y x u t z y x u c Q d )],,,(),,,([d 21−=ρ （11-2）

整个D 由于温度改变需要的热量是：

∫∫∫−=D

V t z y x u t z y x u c Q d )],,,(),,,([21ρ （11-3）

根据牛顿-莱布尼兹公式

∫∂∂=−2

1 12d ),,,(),,,(t t t t

u t z y x u t z y x u （11-4） 前式可以改写为

∫∫∫∫∫∫∫∫

∂∂=∂∂=D t t D t t t V t u c V t t u c Q 2121 d d d d ρρ （11-5） （2）通过L 进入D 的热量1Q 。这里要使用热物理中的傅里叶（Fourier ）热传导定律。该定律证明了：在无穷小时间间隔d t 内通过一个法矢量为n 的无穷小曲面d S ，流向n 所指那一侧的热量为：

t S n

u z y x k Q d d ),,(d ∂∂−= （11-6） 其中),,(z y x k 是该物体在点),,(z y x 处的热传导系数。它恒为正，数值大小取决于组成物体的材料的性质；n 是曲面的外法线，n u ∂∂/是温度函数在),,(z y x 处沿外法线n 的方向导数。我们规定n 所指的那一侧为S d 的正侧，因此该式表示了在t d 时间内从S d 的负侧流向正侧的热量。之所以用负号表示热流的方向与温度梯度方向相反，因为热量总是从温度高的一侧流向温度低的一侧。

现考虑光滑封闭曲面L ，

设在其上确定了一连续变动的单位外法线n ，则在两个时刻1t 和2t 内，经由该物体内任意封闭曲面L 进入D 的热量为：

t S n u k Q t t L

d d 2

1 1∫∫∫ ∂∂−= （11-7）