2.3-最常见的随机过程或随机模型
最常见的随机过程或随机模型
Brown运动或Wiener过程 二项过程 Poission过程 白噪声过程
自回归过程
移动平均过程 混合自回归移动平均过程 利率期限结构或均值回复模型 ARCH类模型与GARCH类模型
1
Brown运动或Wiener过程
引言
Brown运动是1827年英国生物学家Brown在研究 花粉运动时被发现的。
10
1965年,法玛 (Fama)提出了著名 的效率市场假说。该 假说认为,证券价格 对新的市场信息的反 应是迅速而准确的, 证券价格能完全反应 全部信息。
有效 市场 三个 层次
1、弱式效率市场假说 2、半强式效率市场假说 3、强式效率市场假说
根据众多学者的 实证研究,发达国 家的证券市场大体 符合弱式效率市场 假说。一般认为, 弱式效率市场假说 与马尔可夫随机过 程(Markov Stochastic Process)是内在一 致的。因此我们可 以用数学来刻画股 票的这种特征。
维纳过程的性质
z (T ) z (0) i t
i 1
n
[z (T ) – z (0)]也是正态分布 均值等于 0 方差等于T 标准差等于 T 方差可加性
9
为何使用布朗运动?
正态分布的使用:经验事实证明,股票价格 的连续复利收益率近似地服从正态分布 数学上可以证明,具备特征1 和特征2的维 纳过程是一个马尔可夫随机过程 维纳过程在数学上对时间处处不可导和二次 变分(Quadratic Variation)不为零的性质, 与股票收益率在时间上存在转折尖点等性质 也是相符的
特别地,当a(S,t)=Su, b(S,t)=S 时, (3.1)式变为
dS udt dw S
随机过程理论及其应用研究
随机过程理论及其应用研究一、前言随机过程理论是概率论重要分支之一,涉及到各种随机模型和随机变量的演化问题。
它在现代数学、物理学、工程学、生物学、金融学等领域有广泛的应用。
本文将简要介绍随机过程的定义、分类、性质及其应用研究。
二、随机过程的基本概念随机过程是一种数学模型,用来描述随机事件随时间的演化规律。
它是一族随机变量{X(t), t∈T}的集合,其中T是一个表示时间的指标集。
通常,T是时间轴上的一个连续实数集,或者离散的有限集或无限集。
简单地说,随机过程X(t)是在时间t上的一种不确定性量化,X(t)可能随时间t逐渐变化或保持不变。
设X(t)为随机过程的第t个时刻的观测值,通常称X(t)为该随机过程在t时刻的状态。
如果T是一个有限集,那么对应的随机过程称为离散时间随机过程;如果T是一个几何无限集,那么对应的随机过程是连续时间随机过程;如果T是一个更一般的无限集合,那么这个随机过程就是一种更一般的随机过程。
三、随机过程的分类根据时间指标集T的性质,随机过程可以分为离散时间随机过程和连续时间随机过程两种。
1、离散时间随机过程离散时间随机过程定义在离散时间集合上,通常表示为{Xn , n∈ N},其中N是自然数集合, Xn是该过程在第n次观察时的状态。
离散时间随机过程通常被用于表示计数过程、排队过程、随机游走等。
2、连续时间随机过程连续时间随机过程定义在连续时间上,通常表示为{X(t), t >0},其中t∈[0,∞)。
连续时间随机过程通常被用于描述信号传输、通信系统、金融市场等。
四、随机过程的性质随机过程的性质包括时域分布、均值函数和自相关函数。
1、时域分布随机过程在任意时刻t的状态随机变量X(t)的概率分布称为该随机过程在时域上的分布。
时域分布可通过概率密度函数(PDF)、累积分布函数(CDF)和概率质量函数(PMF)来描述。
2、均值函数随机过程的均值函数描述了其期望值随时间的变化规律,通常表示为E[X(t)],它是随机过程X(t)在时间t上的平均值。
随机过程的自回归模型
随机过程的自回归模型随机过程是描述随机事件随时间变化的数学模型。
自回归模型是一种常用的随机过程模型,它假设当前时刻的随机变量值与前一时刻以及过去的随机变量值有关。
一、引言随机过程在众多领域中都有广泛的应用,如金融领域的股票价格变动、通信领域的信号传输、天气预测等。
为了更好地描述随机过程中的随机性和变化规律,研究者提出了各种各样的统计模型。
其中,自回归模型是一种重要的方法。
二、自回归模型的基本概念自回归模型是指当前时刻的随机变量值与前一时刻以及过去的随机变量值之间存在一定的关系。
自回归模型可以用数学表达式表示为:X(t) = c + Σ(ai * X(t-i)) + ε(t)其中,X(t)表示当前时刻的随机变量值,c为常数项,ai为系数,X(t-i)表示过去时刻的随机变量值,ε(t)为噪声项。
三、自回归模型的特点1. 随机性:自回归模型中的噪声项ε(t)具有随机性,能够很好地描述随机过程中的不确定性。
2. 滞后效应:自回归模型中的系数ai表示随机变量值与过去时刻的关系,不同的系数对应不同的滞后效应。
3. 参数估计:自回归模型中的系数ai可以通过最小二乘法等统计方法进行估计,得到模型的参数。
四、自回归模型的应用1. 金融领域:自回归模型可以用于股票价格预测、汇率波动预测等金融领域的分析和建模。
2. 信号处理:自回归模型可以用于信号压缩、降噪等信号处理的应用中。
3. 时序数据分析:自回归模型可以用于时序数据的分析和预测,如天气预测、销售预测等。
五、自回归模型的改进和扩展1. 非线性自回归模型:在自回归模型的基础上引入非线性关系,提高模型的拟合能力。
2. 高阶自回归模型:考虑更多过去时刻的随机变量值,提高模型的时序预测能力。
3. 多变量自回归模型:考虑多个随机变量之间的关系,更好地描述多维随机过程。
六、总结自回归模型是一种常用的随机过程模型,能够很好地描述随机性和变化规律。
它在金融、信号处理、时序数据分析等领域有广泛的应用。
研究生数学建模e题常用的模型
研究生数学建模e题常用的模型
研究生数学建模中常用的模型包括:
1.线性模型:线性回归、线性规划等模型,适用于描述一些简单的线性关系。
2.非线性模型:非线性回归、非线性规划等模型,适用于描述一些复杂的非线性关系。
3.随机模型:包括随机过程、马尔可夫链、随机优化模型等,适用于描述具有随机性或不确定性的问题。
4.动态模型:包括差分方程、微分方程等模型,适用于描述随时间变化的问题。
5.优化模型:包括线性规划、整数规划、多目标规划等模型,适用于求解最优化问题。
6.网络流模型:包括最小生成树、最短路径、最大流等模型,适用于描述网络中的最优路径或流量问题。
7.图论模型:包括图的匹配、图的着色、图的遍历等模型,适用于描述图论问题。
8.排队论模型:包括排队系统、服务系统等模型,适用于描述排队等待问题。
9.时间序列模型:包括ARIMA模型、ARCH模型等,适用于描述时间序列数据的变化规律。
10.复杂系统模型:包括Agent-Based模型、神经网络模型等,适用于描述复杂系统内部的交互和演化过程。
以上模型只是研究生数学建模中常用的一部分,具体的模型选择要根据问题的特点和要求进行决定。
随机过程和随机序列
3. 自相关函数和协方差函数 设X(t1)和X(t2)是随机过程X(t)在t1和t2 二个任意时刻的状态,fX(x1,x2;t1,t2)是 相应的二维概率密度,称它们的二阶 联合原点矩为X(t)的自相关函数,简称 相关函数
RX (t1, t2 ) E[ X (t1 ) X (t2 )]
n
为随机过程X(t)的n维概率密度。
随机过程X(t)和Y(t)的四维联合概率密度
f XY ( x1, x2 , y1, y2 ; t1, t2 , t1 ' , t2 ' ) FXY ( x1, x2 , y1, y2 ; t1, t2 , t1 ' , t2 ' ) x1x2y1y2
FX ( x1, x2 ,, xn ; t1, t2 ,, tn ) P{X (t1 ) x1, X (t2 ) x2 ,, X (tn ) xn }
为随机过程X(t)的n维概率分布函数。
FX ( x1, x2 ,, xn ; t1, t2 ,, tn ) f X ( x1, x2 ,, xn ; t1, t2 ,, tn ) x1x2 xn
则有
RXY (t1, t2 ) mX (t1 )mY (t2 )
当X(t)和Y(t)互相独立时, X(t)与Y(t) 之间一定不相关;反之则不成立。
研究随机过程有两条途经:
侧重于研究概率结构
侧重于统计平均性质的研究
4.2.3 随机过程的特征函数 对于某一固定时刻t,随机变量X(t)的 特征函数就定义为随机过程的一维特 征函数
f X ( x1, x2 ; t1, t2 ) f X ( x1, t1 ) f X ( x2 , t2 )
4.2.2 随机过程的数字特征 1. 数学期望 对于任意的时刻t,X(t)是一个随机变 量,将这个随机变量的数学期望定义 为随机过程的数学期望,记为mx(t), 即
随机模型
与之对应,则称为 (t , w);w );t T 随机过程,一 般简化为 (t ) 。
(t , w), w
定义2 (马尔可夫过程) 设随机过程 (t ) ,如果在已知时 间t系统处于状态x的条件下,在时刻 ( >t)系统所处状态 (t 和时刻t以前所处的状态无关,则称 ) 为马尔可夫过程。 从定义2可知马氏过程只与t时刻有关,与t时刻以前无关。 定义3 (马尔可夫链) 设随机过程 (t ) 只能取 可列个值 r1 , r2 ,rn ,, 把 (t ) rn 称为在时刻 t 系统处于 t E 状态 n (n 1,2,) 若在已知时刻 系统处于 状态的条件下 En ,在时刻 t ( ) 系统所处的状态情况与t时刻以前所处状 (t 态无关,则称 ) 为时间连续,状态离散的马氏过程。而 状态的转移只能在 1,2,) 发生的马氏过程称为马尔 t t n (n 可夫链。 从定义3可知,马氏链是状态离散,时间离散的马尔可夫 过程。
i
i k 2
i1
P ( k ) 有: 定义5 若(2)式中 i j
Pij( k ) Pij
k 1,2,
(6.3)
则称为均匀马氏链 (与第几次转移无关) Pij P( E j / Ei ) P( A(j k ) / Ai( k 1) 即
定义6 转移概率与转移矩阵 令转移概率 Pi j (i 1,2, , j 1,2, ) 为矩阵M 1 的第 i 行 ,第j列元素则有
P11 M 1 P21 P31
P12 P22 P32
P13 P23 P33
(6.4)
M 1 称为马氏链的转移矩阵,其中
Pij 0 Pij 1 j 1
概率论中的随机过程和时间序列
概率论中的随机过程和时间序列随机过程和时间序列是概率论中重要的两个概念,它们在许多领域中有广泛的应用,如统计学、物理学、工程学、经济学等。
随机过程是一个随时间变化的概率分布的集合,而时间序列是一组随时间变化的相关观测值。
一、随机过程随机过程是一组随时间变化的概率分布的集合。
即,对于一个随机过程,每个时间点的随机变量都服从某种概率分布。
随机过程可以看作是一个在时间和状态空间中变化的随机变量。
随机过程可以用数学形式表示为:$$ X(t,\omega) $$其中,t表示时间,ω表示一个样本点或一个事件,X(t,ω)表示在时间点t和样本点ω下的随机变量。
随机过程可以是离散的,也可以是连续的。
根据t的取值范围,随机过程可以分为时域随机过程和频域随机过程。
时域随机过程指的是随机过程在时间上的变化情况,而频域随机过程指的是将随机过程变换到频域中的变化情况。
随机过程的常见模型有马尔可夫过程、布朗运动等。
马尔可夫过程是指在任何时刻t,未来状态的概率分布只与当前状态有关,并且与过去状态无关。
布朗运动是一种连续时间随机过程,它的变化是随机的,但是具有连续性和平稳性。
二、时间序列时间序列是一组随时间变化的相关观测值。
时间序列的分析要求观察数据的时间趋势、季节性、周期性和随机性等方面的规律。
因此,时间序列是一种用来研究随时间变化的数据的分析方法。
时间序列的建模一般有两种方式:统计模型和机器学习模型。
统计模型常用的包括平稳时间序列模型(ARMA、ARIMA、ARCH等)和非平稳时间序列模型(趋势模型、季节模型、协整模型等)。
机器学习模型主要包括回归模型、神经网络模型和支持向量机模型等。
时间序列分析方法中,最常用的是平稳时间序列模型。
平稳时间序列模型是指时间序列具有稳定的均值和方差。
ARIMA模型是一种经典的平稳时间序列模型,主要用于描述时间序列的自相关和移动平均性质。
ARIMA模型具有较好的预测性能和可解释性。
三、应用随机过程和时间序列在许多领域中有广泛的应用,如金融、经济、信号处理、控制系统等。
数学建模中的随机过程与随机优化理论研究
数学建模中的随机过程与随机优化理论研究随机过程是一类重要的数学模型,广泛应用于自然、社会、经济等各领域的研究中。
在数学建模中,随机过程能够对问题进行精确的表述,并且通过对其进行优化能够最优地解决问题。
随机优化理论是基于随机过程的优化理论,通过对随机过程进行分析和改进来提高问题的优化效果。
一、随机过程随机过程是描述随机事件在时间或空间上的演化过程的数学模型。
通俗地讲,就是在一个长时间内,随机事件会发生一些令人难以预料的变化,但是这些变化仍然具有一定的规律性。
随机过程可以用数学语言来描述这种变化的规律性,从而帮助我们更好地理解和应对这种随机性。
随机过程中的随机性可以是在时间上的随机,例如某个事件的发生概率可能在某个时间点会突然增大,也可以是在空间上的随机,例如在一张土地利用图中,某个区域的耕地数量可能会因为自然灾害等原因发生变化。
常见的随机过程有马尔科夫链、布朗运动、泊松过程等等。
二、随机优化理论随机优化理论是在随机过程的基础上发展而来的,旨在通过对随机过程的优化来解决实际问题。
在随机过程中,我们可以使用各种方法来分析变化的规律性,包括概率论、统计学、微积分等等。
而在随机优化理论中,我们需要对这种规律性进行探究和改进,以实现更加准确和有效的优化。
一个典型的随机优化问题是参数优化问题。
在参数优化问题中,我们需要找到一个最好的参数值,以使得某个目标函数达到最优状态。
一般来说,目标函数可能会受到各种随机性的影响,因此需要使用随机优化理论来解决。
三、应用实例随机过程与随机优化理论广泛应用于物理学、统计学、经济学、天文学、信息学、信号处理、控制论等多个领域。
以下列举几个实例:1. 声波传递模型声波传递模型是一种描述声波在空间传递的数学模型。
声波在传递过程中可能受到各种干扰和随机性的影响,因此需要使用随机过程来描述其变化规律,并使用随机优化理论来优化传递过程中的参数,以实现最佳效果。
2. 股市预测分析股市行情的变化受到众多因素的影响,包括政治、经济等多种因素。
几种常用的随机过程
第十讲 几种常用的随机过程10.1 马尔可夫过程 10.1.1马尔可夫序列马尔可夫序列是指时间参数离散,状态连续的马尔可夫过程。
一个随机变量序列x n (n=1,2,…),若对于任意的n 有)|(),...,,|(1121x x F x xx x F n n X n n nX---= (10.1)或)|(),...,,|(1121xx f x xx x f n nXn n nX---=(10.2)则称x n 为马尔可夫序列。
x n 的联合概率密度为)()|( )|()|(),...,,(11221121x f x x f xx f x x f x x x f XXn n Xn nXnX⋅⋅---=(10.3)马尔可夫序列有如下性质:(1) 一个马尔可夫序列的子序列仍为马尔可夫序列。
(2) )|(),...,,|(121xx f x x x x f n nXk n n n n X -+++=(10.4)(3) )|(),...,|(111xX x x X n n n n E E --=(10.5)(4) 在一个马尔可夫序列中,若已知现在,则未来与过去相互独立。
即)|()|()|,(1x x f xx f x x x f r sXn nXrsnX-=,n>r>s (10.6)(5) 若条件概率密度)|(1x x f n n X -与n 无关,则称马尔可夫序列是齐次的。
(6) 若一个马尔可夫序列是齐次的,且所有的随机变量X n 具有同样的概率密度,则称该马尔可夫序列为平稳的。
(7) 马尔可夫序列的转移概率满足切普曼—柯尔莫哥洛夫方程,即)|()|()|(x x fx x fx x fsr Xrn Xsn X⎰∞∞-=,n>r>s (10.7)10.1.2马尔可夫链马尔可夫链是指时间参数,状态方程皆为离散的马尔可夫过程。
1 马尔可夫链的定义 设),2,1( =n X n 为离散时间随机过程,其状态空间},,,{21a a a NI =。
随机过程知识点
随机过程知识点随机过程是现代概率论的重要分支之一,它描述的是一个或多个随机变量随时间的变化规律。
在实际应用中,随机过程经常被用来建立模型,进行仿真以及预测未来的变化趋势等。
随机过程知识点众多,本文将从概念、分类、建模等方面进行探讨。
一、概念随机过程指的是一个定义在时间集合T上的随机变量的集合{Xt:t∈T}。
其中,T表示时间的取值范围,Xt是一个随机变量。
每个时刻t对应一个随机变量Xt,称为随机过程在时刻t的取值。
二、分类根据随机变量的值域,随机过程可以分为离散随机过程和连续随机过程两类。
1. 离散随机过程离散随机过程的取值集合为有限或可数集合。
在离散随机过程中,随时间变化的变量通常被称为时间序列。
离散随机过程可以进一步分为如下几类:(1)马尔可夫链马尔可夫链是最简单的离散随机过程模型,假设当前时刻状态只与前一时刻状态有关。
马尔可夫链的基本性质是:状态转移概率只与当前状态有关,而与历史状态无关。
(2)泊松过程泊松过程是一种间断性随机过程,它描述了单位时间或者单位面积内,某事件发生次数的概率分布。
泊松过程的关键特征是时间和事件之间的指数分布关系,即事件之间的时间间隔是独立且指数分布的。
2. 连续随机过程连续随机过程是取值集合为实数(或实数集合的子集)的随机过程。
在连续随机过程中,随时间变化的变量通常被称为随机过程信号。
连续随机过程可以进一步分为如下几类:(1)布朗运动布朗运动是最基本的连续随机过程,描述了物体在连续介质中的随机运动。
其轨迹连续但不光滑,呈现出瞬时变化的特点。
(2)随机游走随机游走是一种简单的随机过程模型,它描述了物体在一组不断变化的环境下进行的随机运动。
其主要特征是不规则的移动和不可预测性。
三、建模在实际应用中,随机过程的建模是非常重要的。
通过从数学模型中提取重要的特征和参数,可以更好地理解随机过程的行为,从而更好地预测未来的变化。
1. 马尔可夫模型马尔可夫模型是一种广泛使用的随机过程模型,其基本假设是状态的未来只与当前状态有关。
第二章 随机过程的概念与基本类型
2 sin( t ) sin( cos[( s s ) ] t ) DZ 2 cos( s ) cos( t ) sin( s ) sin( t )
2
2
cos[( s t ) ]
例2.2.2 设X(t)=Y+Zt, t>0,Y, Z~N(0, 1),求{X(t), t>0}的一、二维概率密度族。
ft ( x) x 2 exp 2 2 2 1 1 2 1 t
2
2 x exp 2 2 1 t
B X s , t E X s X t m X s m X t E Y Zs Y Zt E Y
举例
例2.2.1 设X(t)=Ycos(t)+Zsin(t), t>0,Y, Z相互 独立,EY=EZ=0,DY=DZ=2。求{X(t), t>0} 的均值函数和协方差函数。 解 :m ( t ) EX ( t ) E [Y cos( t ) Z sin( t )] X
cos( t ) EY sin( t ) EZ 0
1
1 t 2 B 1 st
1 st 2 1 s
2
x 2
T
1
1 s 1 t 1
2 2
x2 1 2 2 1 s
1 s 1 t
2 2
x1 x 2
2 1 t x2
B X ( s , t ) E [( X ( s ) EX ( s ))( X ( t ) EX ( t ))] E [ X ( s ) X ( t )] EX ( s ) EX ( t ) E [ X ( s ) X ( t )]
第二章 随机过程的概念和类型
第二章 随机过程的概念和基本类型2.1 随机过程的基本概念随机过程是随机数学一个十分广泛的分支,它研究的是客观世界中随机现象演变过程的统计规律性.随机过程理论不仅广泛应用于自然科学的各个领域(例如物理学、生物学、电子技术等),而且在社会科学的许多领域也日益受到重视.我们都知道,初等概率论的主要研究对象是随机现象,可以用一个或有限个随机变量来描述随机试验所产生的随机现象.但是,随着科学技术的不断发展,我们必须对一些随机现象的过程进行研究,也就是要考虑无穷多个随机变量,而且解决问题的出发点不是随机变量的独立样本,而是无穷个随机变量的一次具体观测.这时,必须用一簇随机变量才能刻画这种随机现象的全部统计规律,这种随机变量簇就是随机过程.下面先考察几个例子.例 2.1 某人不断地掷一颗骰子,设()X n 表示第n 次掷骰子时出现的点数,1,2,n =⋅⋅⋅,对于任意一个n ,在第n 次掷骰子前不知道试验的结果会出现几点,因此,()X n 是一个随机变量.这样,随机现象可以用一簇随机变量{(),1}X n n ≥来描述.例2.2 设()X t 表示某流水线从开工(0t =)到时刻t 为止的累计次品数,在开工前不知道时刻t 的累计次品数将有多少,因此,()X t 是一个随机变量,假设流水线不断工作,随机现象可以用一簇随机变量{(),0}X t t ≥来描述.例2.3 在天气预报中,若以()X t 表示某地区第t 次统计所得到的该天最高气温,则()X t 是一个随机变量,为了预报未来该地区的气温,我们必须用一簇随机变量{(),0}X t t ≥来描述它的统计规律性.例2.4 在海浪分析中,需要观测某固定点海平面的垂直振动,设()X t 表示在时刻t 该点海平面相对于平均海平面的高度,则()X t 是一个随机变量,我们可以用一簇随机变量{(),0}X t t ≥来描述它的统计规律性.上述例子的共同点是,不是静止地研究某种随机现象,从而研究个别随机变量,而是动态地关心某种随机现象如何随时间变化而发展的,也就是说,需要研究许多随机变量组成的一簇随机变量.一般地,这簇随机变量包含无限多个随机变量,如果这簇随机变量包含有限多个随机变量(例如例 2.1),那么,这类问题用初等概率论中多维随机变量来解决.一簇随机变量描述了随机现象的变化发展过程.为了更深入地研究随机过程的相关性质,我们先给出随机过程的一般定义.定义2.1 设(,ΩF ,)P 是一概率空间,T 是给定的参数,若对于任意t T ∈,有一个随机变量(,)X t ω与之对应,则称随机变量簇{(,),}X t t T ω∈是(,ΩF ,)P 上的随机过程(stochastic process ),简记为随机过程{(),}X t t T ∈,在不致引起混淆的情况下,也可记为()X t .T 为参数集(或指标集),通常表示时间,t 为参数(或指标).需要说明的是:上述定义中的参数集T 可以是时间集,也可以是长度、重量、速度等物理量的集合,随机过程本来通称随机函数,当参数集T 是时间集时称为随机过程,但现在将参数集不是时间集的随机函数也称随机过程,对参数集T 不再有时间限制.在例2.1中,{1,2,}T =⋅⋅⋅,在例2.2, 例2.3和例2.4中[0,)T =+∞,一般地,如果T 由有限多个或可列无限个元素组成的集合,则称{(),}X t t T ∈为离散时间(或离散参数)的随机过程,例2.1是离散时间的随机过程,当T 为有限集时,{(),}X t t T ∈就是概率论中多维随机变量;如果T 是一区间,则称{(),}X t t T ∈为连续时间(或连续参数)的随机过程,例2.2, 例2.3 和例2.4都是连续时间的随机过程.从数学的角度看,随机过程{(),}X t t T ∈是定义在T R ⨯上的二元函数,对固定的t ,(,)X t ω是(,ΩF ,P )上的随机变量,随机变量()X t 所取的值称为随机过程在时刻t 所处的状态(state ),随机过程{(),}X t t T ∈所有随机变量的全体称为随机过程的状态空间(state space ),记为I ;对固定ω,(,)X t ω是定义在T 上的函数,称为随机过程{(),}X t t T ∈的一个样本函数(sample function )或轨道(orbit ),样本函数的全体称为样本函数空间.在例2.1中,{1,2,3,4,5,6}I =;在例2.2中,{0,1,2,}I = ;在例2.3中,(,)I =-∞+∞,在例2.4中[0,)I =+∞.不难看出,在上述例子中,把状态空间作适当扩大,仅仅是为了数学上处理的方便,如果I 是由有限个或可列无限个元素组成的集合,则称{(),}X t t T ∈为离散状态的随机过程,例2.1和例2.2都是离散状态的随机过程;如果I 是一个区间,则称{(),}X t t T ∈为连续状态的随机过程,例2.3和例2.4都是连续状态的随机过程.现将这一分类列表如下:表2-1随机过程的分类随机过程的分类,除了按照参数集和状态集是否可列外,还可以进一步根据过程之间的概率关系进行分类,如独立增量过程、Poisson 过程、Markov 过程、平稳过程、鞅过程等.2.2 随机过程的分布概率论基本内容之一是研究随机变量的分布,随机变量的分布刻画了随机变量的统计规律,分布的表现形式是分布函数(或离散型随机变量的概率函数,或连续型随机变量的概率密度).我们知道,随机过程{(),}X t t T ∈由一簇随机变量组成,当参数集T 为有限集时,随机过程{(),}X t t T ∈由有限个随机变量组成,它本质上与概率论中的多维随机变量相同,可以用多维随机变量的分布函数(或概率函数,或密度函数)来表示随机过程{(),}X t t T ∈的分布;当T 为无限集时,也可以借助有限个随机变量的联合分布来刻画随机过程{(),}X t t T ∈的分布.对于任意一个t T ∈, ()X t 是一维随机变量,其分布函数为(;){()},F x t P X t x x R =≤∈称(;)F x t 为随机过程{(),}X t t T ∈的一维分布函数,显然,对于不同的t ,()X t 是不同的随机变量,因此,(;)F x t 一般也不同,全体一维分布函数组成的集合{(;),:}F x t x R t T ∈∈ F 1称为随机过程{(),}X t t T ∈的一维分布函数簇.对于任意两个12,t t T ∈, ()12(),()X t X t 是二维随机变量,其分布函数为{}21212112212(,;,)(),(),(,)F x x t t P X t x X t x x x R ≤≤∈称1212(,;,)F x x t t 为随机过程{(),}X t t T ∈的二维分布函数,显然,对于不同的12,t t ,()12(),()X t X t 是不同的随机变量,因此,1212(,;,)F x x t t 一般也不同,全体二维分布函数组成的集合212121212{(,;,),(,):,}F x x t t x x R t t T ∈∈ F 2称为随机过程{(),}X t t T ∈的二维分布函数簇.一般地,对于任意n 个12,,,n t t t T ∈ , ()12(),(),,()n X t X t X t ⋅⋅⋅是n 维随机变量,其分布函数为{}121211(,...;,,,)(),,(),n n n n F x x x t t t P X t x X t x ≤⋅⋅⋅≤ 1(,,)n n x x R ⋅⋅⋅∈称11(,,;,,)n n F x x t t 为随机过程{(),}X t t T ∈的n 维分布函数,显然,对于不同的12,,,n t t t ,()1(),,()n X t X t ⋅⋅⋅是不同的随机变量,因此,11(,,;,,)n n F x x t t 一般也不同,全体n 维分布函数组成的集合1111{(,,;,,),(,,):,,}n n n n n F x x t t x x R t t T ⋅⋅⋅⋅⋅⋅⋅⋅⋅∈⋅⋅⋅∈ F n 称为随机过程{(),}X t t T ∈的n 维分布函数簇.定义2.2 {(),}X t t T ∈全体一维分布函数簇F 1、二维分布函数簇F 2⋅⋅⋅的并集F 1n n F ∞== 11111{(,,;,,),(,,):,,,1}n n n n n n F x x t t x x R t t T n ∞=∈⋅⋅⋅∈≥称为随机过程{(),}X t t T ∈的有限维分布函数簇.如果随机过程{(),}X t t T ∈是一个连续状态的随机过程,对于任意,()t T X t ∈通常是连续型随机变量,其密度函数为(;)f x t .称(;)f x t 为随机过程{(),}X t t T ∈的一维密度函数,全体一维密度函数组成的集合称为随机过程{(),}X t t T ∈的一维密度函数簇;一般地,称()1(),,()n X t X t 的密度函数11(,,;,,)n n f x x t t 为随机过程{(),}X t t T ∈的n 维密度函数,全体n 维密度函数组成的集合称为随机过程{(),}X t t T ∈的n 维密度函数簇. 随机过程{(),}X t t T ∈一维密度函数簇、二维密度函数簇 的并集111{(,,;,,:,,,1)}n n n f x x t t t t T n ∈≥ 称为随机过程{(),}X t t T ∈的有限维密度函数簇.类似可以得到离散状态随机过程{(),}X t t T ∈的有限维概率函数簇.随机过程{(),}X t t T ∈有限维分布函数簇、有限维密度函数簇、有限维概率函数簇统称为随机过程{(),}X t t T ∈的有限维分布簇.随机过程{(),}X t t T ∈有限维分布函数簇满足如下两条性质:(1)(对称性) 设12,,,n i i i 为1,2,,n 的任意排列,12,,,n t t t T ∀∈ ,则1111(,,;,,)(,,;,,)n n n n i i i i F x x t t F x x t t =(2)(相容性 consistent )设121,,,,,,,m m n m n t t t t t T +<∀∈ ,则1111(,,,,;,,)(,,;,,)m n m m F x x t t F x x t t ∞∞=反之,对于给定的满足对称性和相容性的分布函数簇,是否存在一个以它作为其有限维分布函数簇随机过程?Kolmogorov 在1931年证明了下述定理肯定地回答了.定理2.1 (Kolmogorov 存在定理)设已知参数集T 满足对称性和相容性的分布函数簇F ,则必存在一概率空间(,ΩF ,P )及定义在上的随机过程{(),}X t t T ∈,它的有限维分布函数簇是F .下面举例说明求随机过程的一维、二维分布.例2.4 设随机过程(),0X t tV t =≥,V 为随机变量,概率函数为{1}0.4,P V =-= {1}0.6P V == 求随机过程()X t 的一维分布函数();12F x 与(;2)F x 及二维分布函数()12,;12,2F x x解 当12t =时,(12)2X V =是离散型随机变量;当2t =时,(2)2X V =是离散型随机变量,它们的概率函数分别为分布函数分别为 0,121;0.4,121221,12x F x x x <-⎧⎪⎛⎫=-≤<⎨ ⎪⎝⎭⎪≥⎩ 和 0,2(;2)0.4,221,2x F x x x <-⎧⎪=-≤<⎨⎪≥⎩当1212,2t t ==时,()()(12),(2)2,2X X V V =是二维离散型随机变量,它的概率函数为因此,()(12),(2)X X 分布函数为 (){}1212,;12,2(12),(2)F x x P X x X x =≤≤=121212121220,0.4,12122,12221,122x x x x x x x x <-<-⎧⎪-≤<≥-≥--≤<⎨⎪≥≥⎩或且且且2.3 随机过程的数字特征定义2.3 随机过程{(),}X t t T ∈,如果对于任意,()t T EX t ∈存在,称()(),X m t EX t t T =∈ (2.1)为随机过程{(),}X t t T ∈的均值函数(expectation function ),简记()m t .定义 2.4 随机过程{(),X t t T ∈,如果对于任意,,s t T ∈ [()()][()()]E X s m s X t m t --存在,称(,)[()()][()()],X C s t E X s m s X t m t -- ,s t T ∈ (2.2)为{(),}X t t T ∈的自协方差函数(self covariance - function ),简称协方差函数,简记(,)C s t ;称 (,)[()()],X R s t E X s X t ,s t T ∈ (2.3) 为随机过程{(),}X t t T ∈的自相关函数(self correlation - function ),简称相关函数,简记为(,)R s t .自协方差函数(,)C s t 是随机过程{(),}X t t T ∈本身在不同时刻状态之间线性关系程度的一种描述,特别地,当s t =时,称为随机过程{(),}X t t T ∈的方差函数(variance function ).2()(,)[()()],X X D t C t t E X t m t t T =-∈ (2.4)由Schwarz 不等式知,随机过程{(),}X t t T ∈的协方差函数和相关函数一定存在,且有下面的关系式(,)(,)()()X X X X C s t R s t m s m t =-.特别地,当均值函数()0X m t ≡时,(,)(,)X X C s t R s t =.从定义可以知道,均值函数()m t 是反映随机过程{(),}X t t T ∈在时刻t 的平均值; 方差函数()X D t 是反映随机过程{(),}X t t T ∈在时刻t 对均值函数()m t 的偏离程度,而协方差函数(,)C s t 和相关函数(,)R s t 反映的是随机过程{(),}X t t T ∈在时刻s 和t 的线性相关程度.例2.5 设随机过程()cos()sin(),0X t Y t Z t t θθ=+>,其中,Y Z 是相互独立的随机变量,且20,EY EZ DY DZ σ====,求{(),0}X t t >的均值函数()m t 和协方差函数(,)C s t .解 由数学期望的性质()[cos()sin()]cos()sin()0EX t E Y t Z t t EY t EZ θθθθ=+=+=又由,Y Z 的相互独立,因此(,)(,)[()()]X X C s t R s t E X s X t ==[cos()sin()][cos()sin()]E Y s Z s Y t Z t θθθθ=++222cos()cos()sin()sin()cos[()]s t EY s t EZ t s θθθθσθ=+=-类似可以定义两个随机过程的互协方差函数和互相关函数.定义2.5 设随机过程{(),}X t t T ∈,{(),}Y t t T ∈,称(,)[()()][()()],,XY X Y C s t E X s m s Y t m t s t T --∈ (2.5)为{(),}X t t T ∈与{(),}Y t t T ∈的互协方差函数(mutual covariance function ),称(,)[()()],,XY R s t E X s Y t s t T ∈ (2.6)为{(),}X t t T ∈与{(),}Y t t T ∈的互相关函数(mutual correlation function ).如果对任意,s t T ∈,有(,)0XY C s t =,则称{(),}X t t T ∈与{(),}Y t t T ∈互不相关.显然有(,)(,)()()XY XY X Y C s t R s t m s m t =- (2.7)例 2.6 设(),Z t X Yt t R =+∈,若已知二维随机变量(,)X Y 的协方差矩阵为2122σρρσ⎛⎫ ⎪⎝⎭,求()Z t 的协方差函数. 解 由数学期望的性质121122(,){[()()][()()]}Z X Y X Y C t t E X Yt m m t X Yt m m t =+-++-+1122{[()()][(()]}X Y X Y E X m Yt m t X m Yt m t =-+--+-2[()()][()()]X X X Y E X m X m E X m t Y m =--+--112[()()][()()]Y X Y Y E t Y m X m E t t Y m Y m +--+--222112112122()XY XY XX YY C t C t C t t C t t t t σρσ=+++=+++例 2.7 设两个随机过程()sin()X t A t ωθ=+与()sin()Y t A t ωθϕ=+-,其中,,,A B ωϕ为常量ϕ为[0,2]π上的均匀分布的随机变量,求12(,)XY R t t .解 设12t t <,则212121201(,)[()()]sin()sin()2XY R t t E X t Y t A t B t d πωθωθϕθπ==++-⎰ 211211210sin(){sin()cos[()]cos()sin[()]}2AB t t t t t t t d πωθωθωϕωθωϕθπ=++--++--⎰222110{cos[()]()2AB t t sin t d πωϕωθθπ=--+⎰ 221110sin[()]sin()cos()}t t t t d πωϕωθωθθ+--++⎰ 21cos[()]2AB t t ωϕ=-- 例 2.8 设()X t 为信号过程,()Y t 为噪音过程,令()()()W t X t Y t =+,则()W t 的均值函数为()()()w X Y m t m t m t =+其相关函数为(,)[()()][()()]w R s t E X s Y s X t Y t =++[()()][()()]E X s X t E X s Y t =+[()()][()()]E Y s X t E Y s Y t ++(,)(,)(,)(,)X XY YX Y R s t R s t R s t R s t =+++上式表明两个随机过程之和的相关函数可以表示为各个随机过程的相关函数之和.特别地,若两个随机过程的均值函数恒为0且互不相关时,有(,)(,)(,)W X Y R s t R s t R s t =+2.4 复值随机过程在工程技术上,常把随机过程表示成复数的形式进行研究更为方便.例如,在许多有关谱函数的运算要用到Fourier 变换,就需要复数形式.定义2.6 设{(),}X t t T ∈,{(),}Y t t T ∈是取值实数的两个随机过程,若对于任意t T ∈, ()()()Z t X t iY t =+其中i =,则称{(),}Z t t T ∈为复随机过程.类似可以定义复随机过程的均值函数、协方差函数、相关函数、方差函数如下: 均值函数: ()[()]()(),Z X Y m t E Z t m t im t t T ==+∈相关函数: 121212(,)()(),,Z R t t E Z t Z t t t T ⎡⎤=∈⎣⎦协方差函数:{}121122(,)[()()][()()]Z Z Z C t t E Z t m t Z t m t =--=121212(,)()(),,Z Z Z R t t m t m t t t T +∈ 方差函数:2()[|()()|](()())(()())(,)Z Z Z Z Z D t E Z t m t E Z t m t Z t m t C t t ⎡⎤=-=--=⎣⎦对于两个随机过程可以定义互相关函数和互协相关函数.互相关函数:12121122(,)[()()]Z Z R t t E Z t Z t =互协相关函数:(){}1212121122111222(,)(),()[()()][()()]Z Z Z Z C t t Cov Z t Z t E Z t m t Z t m t ==--2.5 随机过程的主要类型随机过程可以根据状态空间和参数集离散或连续进行分类,现在我们将根据随机过程的统计特征进一步将随机过程分类,这些常见的随机过程在以后的章节中将作进一步说明,这里只作简单介绍如下:2.5.1 二阶矩过程(two order - moment process )定义2.7 设{(),}X t t T ∈是(取值实数或复值)的随机过程,若对于任意t T ∈,都有2[|()|]E X t <∞(二阶矩存在),则称{(),}X t t T ∈是二阶矩过程二阶矩过程{(),}X t t T ∈的均值函数()()X m t EX t =一定存在,一般假定()0X m t =,这时,协方差函数化为(,)[()()],,X C s t E X s X t s t T =∈.二阶矩过程的协方差函数具有以下性质:(1)(Hermite 性)(,)(,),X X C s t C t s = ,s t T ∈(2)(非负定性)对任意i t T ∈及复数,1,2,,,1i i n n α=≥ 有11(,)0n n X i j i j i j Ct t αα==≥∑∑2.5.2正交增量过程(orthogonal incremental process )定义2.8 设{(),}X t t T ∈是零均值的二阶矩过程,若对于任意1234t t t t T <≤<∈,有 2143[()()][()()]0E X t X t X t X t ⎡⎤--=⎣⎦(2.8) 则称{(),}X t t T ∈为正交增量过程.从定义可以看出,正交增量过程的协方差函数可由其方差确定,且()2(,)(,)min(,)X X X C s t R s t s t σ== (2.9)事实上,不妨设[,]T a b =为有限区间,且规定()0X a =,取12340,,t t t s t b ====,则当a s t b <<<时,有()(()()E X s X t X s ⎡⎤-⎣⎦()()()()(()()E X s X a X t X s ⎡⎤=--⎣⎦0= 因此,(,)(,)()()(,)X X X X X C s t R s t m s m t R s t =-= =()()()()()()()E X s X t E X s X t X s X s ⎡⎤⎡⎤=-+⎣⎦⎣⎦ ()2()(()()()()()X E X s X t X s E X s X s s σ⎡⎤⎡⎤=-+=⎣⎦⎣⎦ 同理,当b s t a >>>时,2(,)(,)()X X X C s t R s t t σ==于是 ()2(,)(,)min(,)X X X C s t R s t s t σ== 2.5.3 独立平稳增量过程(independent stationary incremental process ) 定义2.9 给定随机序列{,1}n X n ≥,如果随机变量12,,X X 相互独立,那么随机序列{,1}n X n ≥为独立过程(或独立随机序列).在例2.1中,如果骰子每次出现的点数是相互独立的,那么得到一个独立随机过程.值得注意的是,就物理意义来说,连续参数独立过程是不存在的,因为,当1t 和2t 很接近时,我们完全有理由说1()X t 和2()X t 有一定的依赖关系,因此,连续参数独立过程只是理想化的随机过程.定义2.10 设随机过程{(),}X t t T ∈,若对任意正整数n 和12n t t t T <<<∈ ,随机变量21321()(),()(),,()()n n X t X t X t X t X t X t ----相互独立,则称随机过程{(),}X t t T ∈为独立增量过程.同独立过程一样,独立增量过程中的参数集T 可以是离散的,也可以是连续的.独立增量过程的直观含义是:随机过程{(),}X t t T ∈在各个不相重叠的时间间隔上状态的增量是相互独立的.在实际应用中,某服务系统在某时间间隔的“顾客”数,电话传呼站的“电话”次数等都可用这种过程来描述.正交增量过程与独立增量过程都是根据不相重叠的时间间隔上增量的统计相依性来定义的,前者增量是不相关,后者增量是独立的.显然,正交增量过程不一定是独立增量过程;而独立增量过程只有在二阶矩存在,且均值为零的条件下才是正交增量过程.定理2.2 设二阶矩过程{(),}X t t T ∈是独立增量过程,若[,),()0T a X a =+∞=,则{(),}X t t T ∈的协方差函数为()2(,)min{,},,X X C s t s t s t a σ=≥. 证明 假设s t <,由()()(),()()()X s X s X a X t X t X a =-=-相互独立性, ()()(,)(),()(),[()()()]X C s t Cov X s X t Cov X s X t X s X s ==-+()()(),()()(),()Cov X s X t X s Cov X s X s =-+2()()X DX s s σ==定义2.11 设随机过程{(),}X t t T ∈,对于任意,,,s t T s t T ττ∈++∈,增量()()X s X s τ+-与()()X t X t τ+-服从相同的分布,则称{(),}X t t T ∈为平稳增量过程.平稳增量过程的直观含义是:随机过程{(),}X t t T ∈在时间间隔(,]t t τ+上状态的增量()()X t X t τ+-仅仅依赖终点和起点的时间差τ,与时间起点无关.如果一个独立增量过程同时又是平稳增量过程,则称它为平稳独立增量过程.平稳独立增量过程是一种很重要的随机过程,后面将反复提到.定理2.3 设随机序列{,0}n X n ≥,且00X =(1){,0}n X n ≥是独立增量过程的充要条件是n X 可以表示为独立随机变量序列的部分和(1)n ≥;(2){,0}n X n ≥是平稳独立增量过程的充要条件是n X 可以表示为独立同分布随机变量序列的部分和(1)n ≥.证明 充分性由定义直接得到,下面证明必要性.令随机变量 1,1n n n U X X n -=-≥,则1,1nn ii X U n ==≥ (1){,0}n X n ≥是独立增量随机过程,对任意n ,增量12,,,n U U U 相互独立,因此,12,,U U 是独立随机变量序列;(2){,0}n X n ≥是平稳独立增量过程时,对任意,m n ,增量,m n U U 同分布,因此,12,,U U 是独立同分布随机变量序列.2.5.4 维纳过程(W i e n e r process )在概率论中我们都知道,正态分布是一种十分重要的分布,正态过程在随机过程中的地位类似于正态随机变量在概率论中的地位,尤其在电讯技术中,正态过程有着十分广泛的应用.定义2.12 设随机过程{(),}X t t T ∈,对任意正整数n 和12,,,,n t t t T ∈()12(),(),,()n X t X t X t 是n 维正态分布,即有密度函数2121211()exp ()()(2)||2n f x x B x B μμπ-⎧⎫=---⎨⎬⎩⎭其中()1212(,,,),(),(),,()T T n n x x x x EX t EX t EX t μ== ,()ij n n B b ⨯=为正定矩阵,{[()()][()()]}ij i i j j b E X t EX t X t EX t =--.则称{(),}X t t T ∈为正态过程或Gauss 过程.19世纪英国植物学家布朗(Brown )发现,浸在水中的微小花粉粒子,受到作不规则运动的水分子的随机碰撞在水面上做不规则的运动,后来,人们把这种运动称为布朗运动.爱因斯坦(Einstein )于1905年第一次给出它的物理解释.1918年,控制论创始人维纳(Wiener )首先对这个随机过程进行了严格的数学论证,奠定了研究这类随机过程的基础.定义2.13 设随机过程{(),}X t t T ∈满足下列条件:(1)(0)0X =;(2)()X t 是独立增量过程;(3)对任意0s t ≤<,增量()2()()0,()X t X s N t s σ-⋅- ,其中,常数20σ>,则称随机过程{(),}X t t T ∈为参数为2σ的Wiener 过程.从定义可以看出,Wiener 过程的参数集[0,)T =∞,状态空间(,)I =-∞+∞,而且Wiener 过程也是平稳增量过程,因此,Wiener 过程是平稳独立增量过程,另外,当s t ≥时,2()()(0,||)X t X s N t s σ-- 依然成立,特别地,当21σ=时,随机过程{(),}X t t T ∈为标准Wiener 过程.定理2.4设随机过程{(),}X t t T ∈为参数为2σ的Wiener 过程.则(1) Wiener 过程是一个正态过程;(2) 22()0,()X X m t t t σσ==; 0t > 且 2121212(,)(,)min(,),X X R t t C t t t t σ==⋅ 12,0t t ≥ (2.10)证明:(2)()()[()(0)]0X m t EX t E X t X ==-=当12t t <时,1212(,)[()()]X R t t E X t X t =1211[()(0)][()()()(0)]E X t X X t X t X t X =--+-21211[()(0)][()()][()(0)]E X t X X t X t E X t X =--+-21t σ=当12t t >时,同样可以得到2122(,)X R t t t σ=因此 21212(,)min(,)X R t t t t σ=例 2.9 设随机过程{(),}X t t T ∈为参数为4的Wiener 过程,定义随机过程()2(3),0Y t X t t =>,则有()Y t 的均值函数为: ()()2(3)0Y m t EY t EX t ===;()Y t 的相关函数为:121212(,)()()4(3)(Y R t t EY t Y t EX t X t ==12121644min(3,min(,)3t t t t =⨯= 2.5.5 泊松过程(Poisson process )在现实世界中有很多例子,例如:盖格记数器上的粒子数,二次大战时,伦敦空袭的弹着点,电话总机所接听的呼唤次数,交通流中事故数,某地区地震发生次数等.这类过程有如下两个性质:一是时间和空间上的均匀性,二是未来的变化与过去的变化没有关系,为了描述这类过程的特性,我们来建立Poisson 过程的模型.定义2.14 给定随机过程{(),0}N t t ≥,如果()N t 表示时间段[0,]t 出现的质点数,状态空间{0,1,2,}I = ,且满足(1)(0)0N =;(2)当s t <时,()()N s N t ≤, 则称{(),0}N t t ≥为记数过程(counting process ).记数过程的样本函数是单调不减的右连续函数(阶梯函数),当跳跃度为1时,称为简单记数过程.简单记数过程表示同一时刻至多出现一个的记数过程.记数的对象不仅仅是电话呼叫次数、来到商店的顾客数,也可表示质点流. 记数过程是时间连续状态离散的随机过程.定义2.15 设随机过程{(),0}N t t ≥是记数过程,如果()N t 满足条件:(1)(0)0N =;(2)()N t 是独立增量过程;(3)对任意0a ≥,0t >,区间(,]a a t +(0a =是应理解为[0,]t )上的增量()()N a t N a +-服从参数为t λ的Poisson 分布,即(){()()},0,1,2,!kt t P N a t N a k e k k λλ-+-=== (2.11) 则,称{(),0}N t t ≥为参数为λ的泊松过程(Poisson process ).0λ>条件(3)表明,()()N a t N a +-的分布只依赖时间t 而与时间起点a 无关,因此,Poisson 过程具有平稳增量性,当0a =时,(){()},0,1,2,,0!kt t P N t k e k k λλλ-===> 因此,Poisson 过程的均值函数为()()N m t EN t t λ==,它表明在时间段[0,]t 出现的平均次数为t λ,λ称为Poisson 过程的强度. 因此,Poisson 过程表明前后时间的独立性和时间上的均匀性,强度λ描述了随机时间发生的频率.有关Poisson 过程的更多结果,后面将进一步论述.2.5.6 马尔可夫过程(Markov process )定义 2.16 设随机过程{(),}X t t T ∈,对于任意正整数n 及12,n t t t <<< 1111{(),,()}0n n P X t x X t x --==> ,且条件分布1111{()|(),,()}n n n n P X t x X t x X t x --≤== 11{()|()}0n n n n P X t x X t x --=≤=> 则称{(),}X t t T ∈为马尔可夫过程(Markov process ).定义中给出的性质称为马尔可夫性,或称无后效性,它表明若已知系统“现在”的状态,则系统“未来”所处状态的概率规律性就已确定,而不管系统“过去”的状态如何.也就是说,系统在现在所处状态的条件下,它将来的状态与过去的状态无关.Markov 过程{(),}X t t T ∈的状态空间和参数集可以是连续的,也可以是离散的.有关Markov 过程的进一步讨论,我们将在第四章进行.2. 5.7 鞅过程(martingale process )最近几十年才迅速发展起来的现代鞅(过程)论是概率论的一个重要分支,它给随机过程论、随机微分方程等提供了基本工具.定义2.17 设参数集{0,1,2,}T = ,如果随机序列{(),0}X n n ≥对任意0,n ≥,且|()|E X n <∞,若[](1)|(1),(2),,()()E X n X X X n X n += (2.12)则称{(),0}X n n ≥为离散参数鞅(discrete parameter martingale ).定义 2.15 设参数集[0,)T =∞,如果随机过程{(),}X t t T ∈对任意|()|,E X t t T <∞∈,若[()|(),](),,..E X s X u u t X t s t a s ≤=> (2.13)则称{(),}X t t T ∈为连续参数鞅(continuous parameter martingale ).上式中,如果将“=”换成“≤”或“≥”,则分别称为离散参数(连续参数)上(或下)鞅.鞅是用条件期望来定义的,关于离散时间鞅,我们可以作下面的直观解释:设()X n 表示赌徒在第n 次赌博时的资本,(1)X 表示最初赌本(这是一常数)而()X n (2)n ≥由于赌博的输和赢是一个随机变量,如果赌博是公平的,那么每次他的资本增益的期望为零,在以后的赌博中,他资本的期望值还是他最近一次赌完的资本数()X n ,用数学模型表示,就是定义中的等式,因此,鞅表示一种“公平”的赌博,上鞅和下鞅表示一方赢利的赌博.例 2.8 设{(),0}Y n n ≥相互独立的随机变量序列,(0)0,Y = 且|()|,E Y n <∞ ()0,0EY n n =≥,令1(0)0,()(),1ni X X n Y i n ===≥∑,则{(),0}X n n ≥是鞅.证明 因为11|()||()||()|n ni i E X n E Y i Y i ===≤<∞∑∑,且[(1)|(0),(1),,()][()(1)|(0),(1),,()]E X n X X X n E X n Y n X X X n +=++[()|(0),(1),,()][(1)|(0),(1),,()]E X n X X X n E Y n X X X n =++()[(1)]()X n E Y n X n =++=定理2.5 设{(),0}X t t ≥是Wiener 过程,则它是鞅.证明:对于任意0s t <<,由独立增量性得 [()()|()][()()]0E X t X s X s E X t X s -=-=因此,对于任意参数01,,,,n t t t t ,01(0)n t t t t =<<<< 有[()|(),0][()()()|(),0]i n n i E X t X t i n E X t X t X t X t i n ≤≤=-+≤≤[()()]()()n n n E X t X t X t X t =-+=习 题 二2.1 设随机变量Y 具有概率密度()f y ,令(),(0,0)Yt X t e t Y -=>>,求随机过程()X t 的一维概率密度及12(),(,)X EX t R t t .2.2 设随机过程()cos()sin()X t A t B t ωω=+,其中ω为常数,,A B 是相互独立且服从正态2(0,)N σ的随机变量,求随机过程的均值和相关函数.2.3 随机过程()X t 的均值函数()X m t 和协方差函数12(,),()X C t t t ϕ为普通函数,令()()()Y t X t t ϕ=+,求随机过程()Y t 的均值和相关函数.2.4 设随机过程2()X t X Yt Zt =++,其中,,X Y Z 是相互独立的随机变量,且均值为0,方差为1,求随机过程()X t 的协方差函数.2.5 设()f t 是一个周期为T 的周期函数,随机变量Y 在(0,)T 上均匀分布,令()()X t f t Y =-,证明:随机过程()X t 满足 01[()()]()()T E X t X t f t f t dt T ττ+=+⎰ 2.6 设随机过程()X t 和()Y t 的互协方差函数为12(,)XY C t t ,证明1212|(,)|()()XY X Y C t t t t σσ≤2.7 设{(),0}X t t ≥是实正交增量过程,(0)0,X V =是标准正态随机变量,对任意的0t ≥,()X t 与V 相互独立,令()()Y t X t V =+,求随机过程{(),0}Y t t ≥的协方差函数.2.8 设,Y Z 是独立同分布随机变量,12{1}{1}P Y P Y ===-=, ()cos()sin(),X t Y t Z t θθ=+t -∞<<∞,其中θ为常数,证明:随机过程()X t 是广义平稳过程,但不是严平稳过程.。
随机过程的基本概念与分类
随机过程的基本概念与分类随机过程是概率论的一个重要分支,在不同领域如金融、通信、生物学等都有广泛的应用。
它描述的是一组随机变量的演化规律,具有许多重要的特性和分类方式。
本文将介绍随机过程的基本概念和分类方法。
一、基本概念随机过程由一个或多个随机变量组成,这些随机变量的取值取决于一个或多个参数,如时间。
随机过程可以定义为函数的族,其中函数的输入参数是随机变量,输出是实数或向量。
常用的随机过程有离散时间和连续时间两种。
在离散时间随机过程中,随机变量类似于离散的时间点,通常用n表示。
每个时间点上都有一个随机变量X(n)与之相关。
连续时间随机过程则对应于时间变量连续变化的情况,通常用t表示。
每个时间点上都有一个随机变量X(t)与之相关。
随机过程的演化可以通过转移概率描述。
转移概率表示从一个时间点到另一个时间点的跳转概率,常用P(i,j)表示从状态i到状态j的概率。
二、分类方法1. 马尔可夫链马尔可夫链是一个简单的、具有重要应用的随机过程。
它具有马尔可夫性质,即未来状态只与当前状态有关,与历史状态无关。
马尔可夫链有着平稳分布,并且可以通过转移概率矩阵进行描述。
2. 马尔可夫过程马尔可夫过程是一种时间连续的随机过程。
它的转移概率与时间无关,但与前一状态有关。
常见的马尔可夫过程有泊松过程、连续时间马尔可夫链等。
3. 马尔可夫决策过程马尔可夫决策过程是一种在马尔可夫过程基础上引入决策的模型。
它包括状态空间、决策空间、转移概率、奖励函数等要素。
马尔可夫决策过程在决策分析、控制理论等领域有广泛应用。
4. 平稳随机过程平稳随机过程是指在统计特性上不随时间改变的过程。
平稳随机过程具有恒定的概率分布和自相关函数。
常见的平稳随机过程有白噪声、自回归过程等。
5. 随机游走随机游走是一种具有随机性的移动方式。
它可以用来模拟股市价格、随机漫步等现象。
随机游走中的步长和方向通常是随机变量,可以是离散的或连续的。
6. 马尔可夫随机场马尔可夫随机场是一种描述多变量间关系的图模型。
随机微分方程(stochastic differential equation,sde)
随机微分方程(stochastic differential equation,sde) 1. 引言1.1 概述随机微分方程(Stochastic Differential Equation,SDE)是一类描述随机现象的微分方程。
相比于传统的确定性微分方程,SDE中包含了一个或多个随机项,能够更准确地描述现实世界中的不确定性和变动性。
SDE在各个领域中广泛应用,特别是金融学、物理学和生物学等领域。
1.2 文章结构本文将从以下几个方面介绍随机微分方程及其应用:定义与基本概念、解随机微分方程的方法与技巧,以及在实际问题中的应用。
具体可以分为三个主要部分:引言、主体内容和结论展望。
1.3 目的本文旨在介绍随机微分方程的基本概念、解法和应用,并探讨其在金融学、物理学和生物学等领域中的实际应用。
通过对随机微分方程的深入了解,读者可以更好地理解和利用该方法来解决实际问题,并对未来研究提出展望。
以上为“1. 引言”部分的内容。
2. 随机微分方程的定义与基本概念2.1 随机过程简介随机过程是一类描述随着时间推移而随机变化的数学模型。
它可以看作是时间参数上的一族随机变量的集合。
随机过程常用于描述具有随机性质的现象,如金融市场中的股票价格、天气预报中的温度变化等。
2.2 随机微分方程的定义随机微分方程是一类描述含有随机项(通常为噪声)的微分方程。
它通常采用以下形式表示:dX(t) = a(X(t), t)dt + b(X(t), t)dW(t)其中,X(t)是未知函数,a(X(t), t)和b(X(t), t)是已知函数,dW(t)表示Wiener 过程(也称为布朗运动或白噪声)。
这个方程表示了X在无穷小时间段dt内发生微小变化dX(t),其中包含一个确定性项a(X(t), t)dt和一个随机项b(X(t), t)dW(t)。
2.3 常见的随机微分方程模型在实际应用中,有许多不同类型的随机微分方程模型被广泛使用。
- Ornstein-Uhlenbeck 过程:该模型描述了维持平衡状态的粒子在受到随机扰动时的演化过程。
随机模型
随机模型在现实世界中, 不确定现象是普遍存在的. 例如, 漂浮在液面上的微小粒子不断地进行着杂乱无章运动, 粒子在任一时刻的位置是不确定的; 又如公共汽车站等车的人数在任一时刻也是不确定的, 因为随时都可能有乘客的到来和离去. 这类不确定现象, 表面看来无法把握, 其实, 在其不确定的背后, 往往隐藏着某种确定的概率规律, 因此, 以概率和数理统计为基础的随机模型就成为解决此类问题最有效的工具之一.依随机规律是否随时间的变化而变化, 随机时模型可分为静态和动态两类, 前者只涉及到随机变量(向量)的概率分布及其数字特征, 后者则要处理随机过程和随机微分方程, 本讲章主要讨论前者.§1 电梯问题有r 个人在一楼进入电梯, 楼上有n 层, 设每个乘客在任何一层出电梯的概率相同, 试建立一个概率模型, 求直到电梯中的乘客下完时, 电梯需停次数的数学期望.分析: 对于此问题, 容易想到该问题为离散随机变量求数学期望的问题, 既然求电梯需停次数的数学期望, 那么电梯需停次数便为要定义的随机变量, 用X 表示, 显然X 的取值范围为}},min{,,2,1{n r ", 然后需计算}{i X P =, 那么所求的答案即为}{},min{1i X p i n r i =⋅∑=, 用上述思路求解电梯问题理论上完全正确, 然而我们却很难得到问题的一个简洁的表达式结果, 原因是古典概率}{i X P =的计算相当复杂, 而且要依据r , n 的不同情况具体来求. 下面我们看一些具体的例子,例1: 取2=r , 3=n 时, 此时X 的取值范围为}2,1{13211{1}33C P X ⋅=== 23222{2}33C P X ⋅=== 从而所求的解为35322311=⋅+⋅ 例2: 取3=r , 2=n 时, 此时X 的取值范围为}2,1{12311{1}24C P X ⋅=== 43222}2{33=−==X P 从而所求的解为47432411=⋅+⋅ 例3: 取3=r , 4=n 时, 此时X 的取值范围为}3,2,1{14311{1}416C P X ⋅=== 3433!6{3}416C P X ⋅=== 313443413!9{2}416C C P X −⋅−⋅=== 从而所求的解为1637166316921611=⋅+⋅+⋅ 上述三例是在r , n 给出具体数据的情况下的计算, 可以看出, 随着r , n 数据的增大, 计算变的愈加复杂, 且没有明显的规律可言. 而当r , n 未给出具体数据时, 用上述思路求解问题, 想要得到具体的表达式就更为困难.下面我们换个角度考虑该问题, 我们将电梯在第i 层是否停下来这一事件作为随机变量i Y , 1=i Y 表示停下来, 0=i Y 表示电梯未停, 其中i 取值为},,2,1{n ", 这样问题便转化为求i Y 的期望之和, 由题意容易得知电梯在任何楼层上是否停留这一概率完全相同即⎪⎪⎩⎪⎪⎨⎧−−−=r r i n n n n Y )1(0)1(11 从而i Y 的期望即为r i nn Y E )1(1)(−−= 那么原问题的解即为]1(1[)(1r ni i n n n Y E −−⋅=∑= 显然要比刚才的方法来得简单, 而且得到了统一的表达式结果, 避免了用第一种方法在计算古典概率时对r , n 大小的具体讨论.我们用前面的三个例子验证上面的结论:例1中, 取2=r , 3=n 时, 得到的结果为35, 而 35])313(1[32=−−⋅ 例2中, 取3=r , 2=n 时, 得到的结果为47, 而 47])212(1[23=−−⋅例3中, 取3=r , 4=n 时, 得到的结果为1637, 而 1637])414(1[43=−−⋅ 通过电梯问题求解的讨论, 可以看出在解决带随机性现象的问题中, 方法的选取是非常重要的, 只有采取合适的方法才会事半功倍.§2钓鱼问题为了估计湖中鱼的数量, 先从湖中钓出r 条鱼做上记号后又放回湖中, 然后再从湖中钓出S 条鱼, 结果发现S 条鱼中有x 条鱼标有记号. 问应该如何估计湖中鱼的数量N ?分析与求解该问题就是要从第二次钓出的标有记号的鱼所占的比例估计出湖中鱼的数量. 首先我们假设放回湖中的鱼在湖中的分布是均匀的. 则第二次钓出的标有记号的鱼数X 是一个随机变量, X 服从超几何分布{}x s x r N r s NC C P X x C −−⋅== (*) 其中x 为整数, 且],min[)](,0max[s r x r N s ≤≤−−. 用),(N x L 表示(*)式的右端, 则取使),(N x L 达到极大值的N 作为N 的估计量. 直接对N 求导考察极值比较困难, 我们用比值法来研究),(N x L 的变化(,)(,)(,1)()()L x N N r N s A x N L x N N N r s x −−==⋅−−−− 22()()N r s N rs N r s N Nx−++=−++ (**) 从(**)式看出, 当且仅当xrs N <时, )1,(),(−>N x L N x L , 而当且仅当xrs N >时, )1,(),(−<N x L N x L . 因此),(N x L 在x rs 附近取得极大值, 于是N 的估计值N ˆ为][x rs 或[x rs +1, 取使得),(N x L 的值更大的一个即可.上面的求解方法实际上是运用了概率统计中的极大似然原理, 即现在这个事件发生了, 那么客观情况使得它最有可能发生.下面我们换个角度考虑上述问题, 既然假设放回湖中的鱼在湖中的分布是均匀的, 我们可以认为湖中整个鱼群中含带有记号的鱼群比例与湖中任意一部分鱼群中含带有记号的鱼群比例完全相同, 即sx N r = 从而srx N =, 取整即与上述分析所得的结果完全相同. 从这一问题, 我们可以学到估计类似问题的一种实际操作方法.§3 广告中的数学在我们的现实生活中, 广告无所不在.广告给商家带来了丰厚的利润, 广告中蕴藏着诸多学问.以房产销售广告为例, 房产开发商为了扩大销售, 提高销售量, 通常会印制精美的广告分发给大家.虽然买房人的买房行为是随机的, 他可能买房, 也可能暂时不买, 可能买这家开发商的房子, 也可能买另一家开发商的房子, 但与各开发商的广告投入有一定的关联.一般地, 随着广告费用的增加, 潜在的购买量会增加, 但市场的购买力是有一定限度的.表3.1给出了某开发商以往9次广告投入及预测的潜在购买力.表3.1 广告投入与潜在购买力统计( 单位: 百万元)广告投入0.2 0.4 0.5 0.52 0.56 0.65 0.67 0.69 1购买力10340 10580 10670 10690 10720 10780 10800 10810 10950下面从数学角度, 通过合理的假设为开发商制定合理的广告策略, 并给出单位面积成本700元, 售价为4000元条件下的广告方案.模型假设(1) 假设单位面积成本为1p 元, 售价为2p 元, 忽略其他费用, 需求量r 是随机变量, 其概率密度为()p r .(2) 假设广告投入为p 百万元, 潜在购买力是p 的函数记作()s p ,实际供应量为y .模型建立开发商制定策略的好坏主要由利润来确定, 好的策略应该获得好的利润(平均意义下), 为此, 必须计算平均销售量()E x .0()()()y y E x rp r dr yp r dr +∞=+∫∫上面右边第二项表示当需求量大于等于供应量时, 取需求量等于供应量.因此, 利润函数为21()()R y p E x p y p =−−利用0()1p r dr +∞=∫得到2120()()()()yR y p p y p p y r p r dr =−−−−∫ ( 3.1)上式中, 第一项表示已售房毛利润, 第二项为广告成本, 第三项为未售出房的损失.模型求解为了获得最大利润, 只需对(3.1)式求导并令其为零, 设()R y 获得最大值时y 的最优值为*y , 则2120()()()0y dR y p p p p r dr dy=−−=∫因此, *y 满足关系式*2102()y p p p r dr p −=∫ ( 3.2) 通过(3.2)式知道, 在广告投入一定的情况下, 可以求出最优的供应量, 但依赖于需求量的概率分布.为使问题更加明确, 增加如下假设:(3) 假设需求量r 服从[0,()]U s p 分布, 即10()()()0r s p s p p r ⎧≤≤⎪=⎨⎪⎩其他 (3.3) 将(3.3)代人(3.2)得到212*()p p y s p p −= ( 3.4) 即最优的供应量等于毛利率与由广告费确定的潜在购买力的乘积.将( 3.4)式代入( 3.1)式, 得到最大利润为2212()(*)()2p p R y s p p p −=− ( 3.5) 对(3.5)式关于p 求导, 得驻点*p 满足的方程为22212'(*)()p s p p p =− ( 3.6) 因此, 只要知道了潜在购买力函数, 就可以给出最优的广告投入.下面根据开发商获得的相关数据, 来确定潜在购买力函数. 通过对表3.1数据分析, 得知其符合log istic 型曲线增长率, 经拟合得到52()10/(9)p s p e −=+ ( 3.7)记52221210()p l p p −=×− 将(3.7)式代入(3.6)式, 当1180l −>时, 求得*11ln(19ln 22p l l =−−−+ (3.8) 将120.0007,0.004p p ==代入(3.8)式得到*0.49p =(百万元).§4 报童的策略背景简介报童问题即单周期库存问题(Single-Period Problem), 是供应链管理中最重要的模型之一, 其历史可以追溯到1888年著名经济学家Edgeworth 应用它解决银行的现金流(cash-flow)问题, 1955年, Whitin 首次建立了受价格影响的报童问题模型. 目前, 报童问题在生产、服务、管理和金融等领域成功地取得了广泛的应用. 现在报童问题衍生出许多扩展问题.与众多扩展模型相比, 经典报童问题模型是最简单最基本的问题, 它可以描述为: 报童每天早晨以单位批发价b 从报社买进报纸, 然后以单位零售价a 出售, 晚上将没有卖掉的报纸当作废品以价格c (c b a >>)处理掉. 同时假设:(1) 报童拥有购买足够多报纸的资金;(2) 报纸过剩只能以低于零售价的价格v 处理;(3) 报纸供应不足, 会遭受缺货惩罚;(4) 不计其他费用(如交通费、摊位费等).报童应该如何确定订购量而获得最高的利润呢? 显然, 报童应该根据市场需求量来确定订购量, 而市场需求量是随机的. 假设报童通过经验已经掌握了市场需求量的随机规律, 我们就可以建立随机优化模型来求解报童问题了.报童每天清晨从报社购进报纸零售, 晚上将没有买掉的报纸退回. 每份报纸的购进价为元b, 零售价为a元, 退回价a−元, 退回一份报纸赔为c元,. 报童售出一份报纸赚bcb−元. 报童每天如果购进的报纸太少, 不够卖时会少赚钱, 如果购进的报纸太多, 卖不完时会赔钱. 试为报童筹划每天应如何确定购进的报纸数使得收益最大?模型一问题分析报童应该根据需求量确定购进量, 而需求量是随机的,所以这是一个风险决策问题. 假定报童已经通过自己每天的卖报经验或其它渠道掌握了需求量的分布规律, 需求量r 为一连续型随机变量, 密度函数为)(r f ,假设每天的购进量为n , 由于需求量r 是随机的, 可以小于n 、等于n 或大于n , 这就导致报童每天的收入也是随机的, 所以作为优化模型的目标函数, 不能是报童每天的收入函数, 而应该是他长期卖报的日平均收入. 从概率论大数定律的观点看这相当于报童每天收入的期望值, 以下简称它为平均利润.模型建立显然, 若n r >, 则以a 价售出r 份报纸, 以c 价售出n r −份报纸, 若n r ≤, 则全部n 以a 价售出, 故平均利润为 0()[()()()]()n F n a b r b c n r f r dr =−−−−∫()()n a b nf r dr +∞+−∫ 0()()()()na b n a c n r f r dr =−−−−∫0()()()()n dF n a b a c f r dr dn=−−−∫ 由一元函数极值存在的必要条件可知()0dF n dn=, 得 0()n a b f r dr a c−=−∫ 而22()()()0d F n a c f n dn=−−<. 从而知满足上式的n 可以使平均利润达到最大. 易知上式等价于()()n nf r dr a b b cf r dr +∞−=−∫∫ 上式左边是报童订购n 份报纸时, 不能将它卖完的概率与能将它卖完的概率之比, 右边则表示卖出一份报纸的盈利与退回一份报纸亏损之比, 该式表明最优购进量是使这两个比相等的购进量.实例: 如3.0,6.0,1===c b a , 需求量r 服从正态分布)10,100(2N , 则当报童的订报量n 满足0()43()n nf r dr a b b c f r dr +∞−==−∫∫或等价于 04()7n a b f r dr a c −==−∫ 不难计算得到102n =时长期平均收益最大. 若3.0,6.0,2.1===c b a 此时 02()3na b f r dr a c −==−∫, 104n =时长期平均收益最大. 若2.0,6.0,8.0===c b a 此时01()3n a b f r dr a c −==−∫,96n =时长期平均收益最大.再如, a b −=0.3, b c −=0.1, 需求量r 为服从[2000,4000]上均匀分布的连续型随机变量, 其密度函数为1200020004000()0x f r ≤≤⎧=⎨⎩其它 由0()n a b f r dr a c−=−∫知, 2000()n f r dr ∫a b a c −=−=0.30.30.1+=0.75, 于是, 有 n -2000=1500, 因此, 报童的最优策略是订购3500份.模型二(报童问题)报童每天要到邮局去订报, 出售一份报纸可获得利润()A a b =−元, 但如卖不出退回邮局, 每份报纸要损失()B b c =−元. 根据以往经验, 得知每天需求量为k 份的概率为k p . 问报童每天应订购多少份报纸, 才能使它获利的期望值最大.如果卖报童问题中的顾客每天需求量X 是一个离散型随机变量, 设报童每天订购的份数为n 份, 于是有()k P X k p ==, 记出售一份报纸可获得利润()A a b =−(元)退回邮局一份报纸要损失()B b c =−(元)则报童每天的利润()f X 可用下列公式来表示:()()An X n f X AX n X B X n≥⎧=⎨−−<⎩ 因此报童获利的期望值为:()C n =0(())()()k E f X f k P X k ∞===∑=10[()]n k k k k n Ak n k B p Anp −∞==−−+∑∑ (4.1)报童需要做出的决策: 确定一个订购数n , 使得(())E f X 最大. 于是所求模型为10max ()[()]n k k k k n C n Ak n k B p Anp −∞===−−+∑∑我们采用边际分析法来求解, 也即利用价格结构来检验和判断在什么情况下, 再多订一份报纸是合算的.假设报纸订购数取n 份是合算的, 现考察再多订一份报纸是否合算, 也就是考察第n +1件报纸的利润期望值. 第n +1份报纸售出时, 获利为A 元, 售不出去时获利为B −元. 因此, 此时多订一份报纸的利润期望值为:()(1)()Ap B p A B p B +−−=+−其中(1)p P X n =≥+. 所谓合算, 就是利润期望值大于零.故由()0A B p B +−>, 可解得售出概率p 应满足下述不等式B p A B>+ (4.2) 其中100(1)()1()1n n i i n i i p P X n P X i P X i p ∞=+===≥+===−==−∑∑∑.即001nn i i i i B A p p A B A B ==−>⇒<++∑∑ 于是, 报纸的最佳订购量*n 应满足:*()B P X n A B≥>+, (4.3) 其中B A B+称为临界值. 那么, 最好的n 应满足(1)(),(1)()c n c n c n c n −≤+≤即10n n i i i i A p p A B −==≤≤+∑∑ 例如, 已知某种报纸每天需求量N j 的概率分布如下: 表1 需求量N j 0 1 2 3 4 5 6 7 8 P (N j ) 0.05 0.10 0.100.250.200.150.050.05 0.05 每出售一份报纸, 可获利4角; 若当天卖不掉, 每份报纸将损失3角. 试问每日应进多少份报纸?因为我们引进了需求量随机变量X , 所以我们将表1作一点修改并就在表中进行计算, 如表2:表2n0 1 2 3 4 5 6 7 8P(X =n) 0.05 0.10 0.100.250.200.150.050.05 0.05 P(X≤n) 0.05 0.15 0.250.500.700.850.900.95 0.10 P(X≥n+1) 0.95 0.85 0.750.500.300.150.100.05 0.00现在本问题中, A=4, B=3,所以BA B+≈0.43, 在n=3时,P(X≥3+1)= 0.50满足式(3), 所以最佳订购量n*=4份.该模型也可以理解为单周期随机库存问题, 即假定在一个周期末库存的货物对下一个周期没有任何价值, 即模型适用于仅有一次机会存贮以供需求的产品如时装、新鲜食品、月饼等.§5 最佳订票问题一.问题提出在激烈的市场竞争中, 航空公司为争取更多的客源而开展的一个优质服务项目是预订票业务. 公司承诺, 预先订购机票的乘客如果未能按时前来登机, 可以乘坐下一班飞机或退票, 无需附加任何费用. 当然也可以订票时只订座, 登机时才付款, 这两种办法对于下面的讨论是等价的.设某种型号的飞机容量为n, 若公司限制预定n张机票, 那么, 由于总会有一些订了机票的乘客不按时来登机, 致使飞机因不满员飞行而利润降低, 甚至亏本, 如果不限制订票数量呢, 那么当持票按时前来登机的乘客超过飞机容量时, 必然会引起那些不能登机飞走的乘客(以下称被挤掉者)的抱怨. 公司不管以什么方式予以补救, 都会导致受到一定的经济损失, 如客源减少, 或挤到以后班机的乘客, 公司要无偿供应食突宿或者付给一定的赔偿金等. 这样, 综合考虑公司的经济利益, 必然存在一个恰当的订票数量和限额.假设飞机容量为300, 乘客准时到达机场而未乘上飞机的赔偿费是机票价格的%10, 飞行费用与飞机容量、机票价格成正比(由统计资料知, 比例系数为0.6, 乘客不按时前来登机的概率为0.03), 请你:(1)建立一个数学模型, 给出衡量公司经济利益和社会声誉的指标, 对上述预定票业务确定最佳的预定票数量.(2)考虑不同客源的不同需要, 如商人喜欢上述这种无约束的预定票业务, 他们宁愿接受较高的票价; 而按时上下班的雇员或游客, 愿意以若不能按时前来登机, 则机票失效为代价, 换取较低额的票价. 公司为降低风险, 可以把后者作为基本客源. 根据这种实际情况, 制定更好的预订票策略.二.模型的假设及符号说明1、模型的假设①假设预订票的乘客是否按时前来登机是随机的.②假设已预订票的乘客不能前来登机的乘客数是一个随机变量.③假设飞机的飞行费用与乘客的多少无关.2、符号说明n: 飞机的座位数, 即飞机的容量;g: 机票的价格;f: 飞行的费用;b: 乘客准时到达机场而未乘上飞机的赔偿费;m: 售出的机票数;k: 已预订票的乘客不能前来登机乘客数, 即迟到的乘客数, 它是一个随机变量;p: 已预订票的m个乘客中有k个乘客不能按时前k来登机的概率;p: 每位乘客迟到的概率;()m P j : 已预订票前来登机的乘客中至少挤掉j 人的概率, 即社会声誉指标;S : 公司的利润;ES : 公司的平均利润.三.问题的分析及数学模型 1、问题的分析通过上面引进的符号易知, 赔偿费g b 1.0=, 习行费用ng f 6.0=, 每位乘客迟到概率03.0=p , 已预订票的m 个乘客中, 恰有k 个乘客不能按时前来登机, 即迟到的乘客数k 服从二项分布()p m B ,, 此时,(1)(0,1,2,,)k km k k m p C p p k m −=−="当n k m ≤−时, 说明k m −个乘客全部登机, 此时利润f g k m S −−=)(当n k m >−时, 说明有n 个乘客登机, 有n k m −−个乘客没有登上飞机, 即被挤掉了, 此时利润b n k m f ng S )(−−−−=根据以上的分析, 利润S 可表示为:⎩⎨⎧−<>−−−−−−≥≤−−−=)()()()(n m k n k m bn k m f ng n m k n k m fg k m S 迟到的乘客数1,,2,1,0−−=n m k "时, 说明有n k m −−个乘客被挤掉了;迟到的乘客数m n m n m k ,,1,"+−−=时, 说明已来的k m −个乘客全部登机了.于是平均利润∑∑−=−−=−−+−−−−=mnm k kn m k kpf g k m pb n k m f ng ES ])[())((1因为∑∑∑∑∑∑−−=−−=−−==−−=−=+−−−−=−−−−=−−101110)()()()()1)((])[(n m k kn m k k n m k km k kn m k kmnm k kkpgk gE p f mg f mg kp kp g p f mg pf g k m所以111(())()()()m n kk m n m n k kk k ES ng f m k n b pmg f mg f p gE k gkp−−=−−−−===−−−−+−−−−+∑∑∑1010(())(()())(())()()m n kk m n kk m E k g f ng f m k n b mg f gk pm E k g f b g m k n p −−=−−==−−+−−−−−−+=−−−+−−∑∑由于k m k kmk p p C p p n B k −−=)1(),,(~, 可知, 随机变量k的数学期望mp k E =)(, 此时,∑−−=−−−−+−−−=10)1()()()1(n m k k m k k mp p Cn k m g b f mg p ES2、数学模型通过以上对问题的分析, 可以在一定的社会声誉指标()m P j 范围内, 寻求合适的m , 根据Ng f 6.0=的关系, 使得目标函数f ES 达到最大, 即 1011max [(1)(1)()(1)]10.61[0.97 1.1(300)(1)]1180m n k km k m k m n k km k m k ES b p m m k n C p p f N g m m k C p p −−−=−−−==−−+−−−−=−−−−−∑∑下面考虑社会声誉指标.由于j k n m ++=, 所以j n m k −−=, 即当被挤掉的乘客数为j 时, 等价的说法是恰有j n m −−个迟到的乘客.公司希望被挤掉的乘客人数不要太多, 被挤掉的概率不要太大, 可用至少挤掉j 人的概率作为声誉指标, 相应地k 的取值范围为j n m k −−=,,2,1,0", 社会声誉指标()(1)m n jk k m k j m k P m C p p −−−==−∑四、模型求解为了对模型进行求解, 可以分别给定m , 比如,350,,306,305"=m , 计算f ES /, 同时, 给定j , 比如取5=j , 计算社会声誉指标()m P j , 从中选取使f ES /最大, 且社会声誉指标()m P j 小于等于某个α(比如取05.0=α)最佳订票数m .下面给出MATLAB 计算程序.%飞机最佳订票策略ch43 %文件名: ch43.m%m表示售出的票数; Es表示平均利润; p表示声誉指标; for m=305: 325sm=0;p=0;for k=0;m-305pp=(prod(m-k+1:m)/prod(1:k))*0.03^k*0.97^(m-k);p=p+pp;sm=sm+(m-k-300)*pp/prod(1:k);endEs=(1/180)*[0.97*m-1.1*sm]-1;mEspend执行后可输出以下结果:m ES P305 0.6436 9.2338e-005306 0.6490 9.3723e-004307 0.6543 0.0048308 0.6596 0.0167309 0.6649 0.0442310 0.6703 0.0952311 0.6756 0.1742312 0.6810 0.2796313 0.6864 0.4028314 0.6917 0.5314315 0.6971 0.6525316 0.7024 0.7566317 0.7078 0.8388318 0.7132 0.8890319 0.7185 0.9399320 0.7239 0.9661321 0.7293 0.9818322 0.7347 0.9907323 0.7400 0.9954324 0.7454 0.9979325 0.7508 0.9990从计算结果易见, 当m=309时, 社会声誉指标()53090.04420.05,p =<当m =310时, 社会声誉指标()53100.09520.05,p =>所以为了使尽/ES f 量大, 且要满足社会声誉指标()50.05,p m <则最佳订票数可取m =309.。
2.3-最常见的随机过程或随机模型
3
二项分布是指随机变量满足概率分布
P(ξ = k ) = C p (1- q)
k n
k
n- k
其中,k=1,2, …,0<p<1,q=p-1。 二项过程实质上是将二项分布作为一个过程来描 述金融资产价格变化的。
4
假设股票价格在t时刻为S(t),当时间变化到 t+t时,价格要么以概率p从S上涨到uS(u >1), 要么以概率q下降到dS(d<1);时间为t+2t时有 三种可能:u2S、udS、d2S,以此类推,见树型 结构
15
利率期限结构或均值回复模型
在金融市场中,许多情况下的金融资产价格的变化,随着时 间的推移常常趋于某个长期平均水平,称为均值回复现象, 例如利率的变化就常常如此。具体的利率期限结构或均值回 复模型定义为
dS (u S )dt S dt
其中λ>0,ε服从标准正态分布。当股票价格S低于均值μ时, μ-S取正值,即S具有正的漂移率,dS将会变为正值。反之, 当股票价格S高于均值μ时,μ-S取负值,即S具有负的漂移 率,dS将会变为负值。尽管变化过程中价格可能会偏离均 值μ ,但长期来看S都会向均值μ靠近。过程中偏离的程度 由参数λ>0决定的。注意:资产价格表现出来的某种长期可 预测性,与市场有效性的假定是不符合的。
14
混合自回归—移动平均过程 若时间序列{t }在t时刻,不仅与其以前的自身值 有关,而且与以前时刻的冲击或扰动存在着一定的 依存关系,则称为混合自回归—移动平均过程,其 一般形式(记作ARMA(m,n))为 t =a1t -1+ a2t -2+…+ amt -m+t +b1t -1+ b2t -2+…+ bn t –n
典型的随机过程
2 X
0 0 rX ( )d rX (0 ) 0.05
随机过程的遍历性
平稳随机过程满足:mX
P
mX
P
RX ( ) RX ( )
rX ( )
1
0
0
相关时间示意图
mˆ X
1 2T
T
x(t )dt
T
Rˆ X
(
)
1 2T
T
x(t ) x(t)dt
n
fx xi过程 X (t) Acos0t B sin 0t ,
其中A、B是两个独立的正态随机变量,且
有 E[A] E[B] 0 ,E[ A2 ] E[B2 ] ,2 0为常数,
求此过程的一维概率密度。
例2、一零均值高斯过程X(t),其协方差函数为:
求在时刻t21随机过程的基本概念及定义22随机过程的统计描述23平稳随机过程24随机过程的联合分布和互相关函数25随机过程的功率谱密度26典型的随机过程随机过程随机过程随机过程的统计特性随机过程的统计特性随机过程是样本函数的集合随机过程是随机变量的集合随机过程的概率分布随机过程的数字特征随机过程的统计特性随机过程的统计特性随机过程的特征函数离散形式
( )
E[lim T
1 2T
XT () 2 ]
平稳随机过程:维纳-辛钦定理 RX ( ) GX ()
性质: 实平稳随机序列的功率谱是实的、非负的偶函数。
RX
(0)
E[ X
2
(n)]
1 2
GX ()d
互功率谱定义:GXY ()
E{lim T
条件过程模型
条件过程模型条件过程模型,也被称为随机过程模型,是一种用于描述随时间变化的随机量的数学模型。
它是现代概率论和统计学中一种重要的工具,可以广泛应用于信号处理、时间序列分析、风险管理等领域。
本文将从以下几个方面探讨条件过程模型。
一、基本概念条件过程模型是一类时间序列模型,其描述的是随着时间推移变化的随机量。
随机过程可用于描述众多自然现象,如热力学系统、金融市场、天气预报等。
在条件过程模型中,时间通常被视为连续的。
简单来说,条件过程模型描述的是在给定一些已知历史信息的情况下,未来某些事件发生的概率。
二、建模方法条件过程模型的建模方法通常涉及以下步骤:1.选择过程类型:根据研究对象、数据特征等信息,选择合适的过程类型,如布朗运动、随机游走、自回归过程等。
2.确定过程参数:过程的参数通常通过最大似然估计等方法确定。
3.进行预测:在确定了过程参数后,可以利用已知历史信息进行预测,推断未来的概率。
三、常见的条件过程模型1.布朗运动模型:布朗运动是一种最简单、最基本的随机过程模型。
它描述的是一个粒子在无规则的热运动下的轨迹。
布朗运动模型被广泛应用于金融市场中的股票价格、汇率等的预测。
2.马尔可夫过程模型:马尔可夫过程模型是指一类具有状态转移概率的随机过程。
在这种模型中,当前状态的概率只与前一时刻状态有关,与之前状态无关。
马尔可夫过程模型被广泛应用于天气预报、环境污染等领域。
3.随机游走模型:随机游走模型描述的是一种类似于赌徒随机行进的过程。
在随机游走模型中,每步的移动是由随机因素和先前位置决定的。
该模型被广泛应用于股票价格和商品价格的预测。
四、应用领域条件过程模型在各个领域都有很广泛的应用。
以下列举几个典型的应用领域:1.信号处理:条件过程模型可以应用于信号处理中的噪声抑制、谱估计、信源分离等方面。
2.时间序列分析:条件过程模型是时间序列分析的基础,可以用于预测股票价格、汇率、气象数据等。
3.风险管理:条件过程模型可以用于风险管理中市场风险、信用风险、操作风险等方面的预测和计量。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ξt=aξt-1+ εt,t=1,2, …, ξ ,
其中,a为常数,εt为白噪声过程,称为扰动项。当|a|<1 其中, 为常数, 为白噪声过程,称为扰动项。 为常数 时为平稳过程; 时称为随机游走过程; 时为平稳过程;a=1时称为随机游走过程;|a|>1为非平 时称为随机游走过程 为非平 稳过程。 稳过程。
i = ∑c2 pi (1 − p) 2−i u i d 2−i S i =0 2
在t + n∆t 时刻,股票的期望价格为: ∆ 时刻,股票的期望价格为:
n
E ( S t + n∆t ) =
å
i= 0
C p (1- p) u d S
6
i n
i
n- i
i
n- i
Poission过程 过程
引言: 引言: Brown运动是用以描述连续时间下金融资产价格 运动是用以描述连续时间下金融资产价格 运动的, 运动的,但金融资产价格并不都是随时间而连续变 化的,有时会出现跳跃, 化的,有时会出现跳跃,Poission过程就是经常 过程就是经常 用以模拟跳跃的一类随机过程。 用以模拟跳跃的一类随机过程。
11
自回归过程
按时间次序排列的随机过程{ξt}( t=1,2,…)称为时间序 按时间次序排列的随机过程 ξ , , 称为时间序 列。 若时间序列是相互独立的, 若时间序列是相互独立的,则说明事件后一刻的行为与前一 刻毫无关系,即系统无记忆性。 刻毫无关系,即系统无记忆性。 若情况相反,则前后时刻事件之间就有一定的依存性。其中 若情况相反,则前后时刻事件之间就有一定的依存性。 最简单的关系就是事件后一刻的行为只与前一刻的行为有关, 最简单的关系就是事件后一刻的行为只与前一刻的行为有关, 而与其前一刻以前的行为无直接联系, 主要与ξ 相关。 而与其前一刻以前的行为无直接联系,即ξt主要与ξt -1相关。 从记忆的角度理解,是最短的记忆,即一期记忆, 从记忆的角度理解,是最短的记忆,即一期记忆,描述这种 关系的模型称为一阶自回归过程,记为AR( ), ),即 关系的模型称为一阶自回归过程,记为 (1),即
13
移动平均过程 自回归过程表示在t时刻的事件ξt 只与其以前的响 自回归过程表示在 时刻的事件ξ 时刻的事件 有关, 应ξt -1,ξt -2,…,ξt -m 有关,而与以前时刻的扰 , 动无关。若时间序列{ξ 与其以前的冲击或扰动 动无关。若时间序列 ξt }与其以前的冲击或扰动 有关, εt -1,εt -2,…,εt -n有关,而与以前时刻的响应 , 无关,那就是n阶移动平均过程 记为MA(n),即 阶移动平均过程, 无关,那就是 阶移动平均过程,记为 即 ξt = b0+εt +b1εt -1+ b2εt -2+…+ bnεt –n ε t=1,2,… , , 当|bj|<1时,表示冲击在一段时间内会消失; 时 表示冲击在一段时间内会消失; 表示冲击永远保持下去; |bj|=1表示冲击永远保持下去;|bj|>1表示冲 表示冲击永远保持下去 表示冲 击将放大,其中i=1,2,…,n。 击将放大,其中 , , 。
dS = λ (u − S ) dt + σSε dt
其中λ>0,ε服从标准正态分布。当股票价格 低于均值 时, , 服从标准正态分布 当股票价格S低于均值 服从标准正态分布。 低于均值µ时 其中 µ-S取正值,即S具有正的漂移率,dS将会变为正值。反之, 取正值, 具有正的漂移率, 将会变为正值 反之, 将会变为正值。 取正值 具有正的漂移率 当股票价格S高于均值 高于均值µ时 取负值, 当股票价格 高于均值 时,µ-S取负值,即S具有负的漂移 取负值 具有负的漂移 率,dS将会变为负值。尽管变化过程中价格可能会偏离均 将会变为负值。 将会变为负值 但长期来看S都会向均值 靠近。 都会向均值µ靠近 值µ ,但长期来看 都会向均值 靠近。过程中偏离的程度 由参数λ>0决定的。注意:资产价格表现出来的某种长期可 决定的。 由参数 决定的 注意: 预测性,与市场有效性的假定是不符合的。 预测性,与市场有效性的假定是不符合的。
10
P(ξt −ξs = k) =
λk (t − s)k e−λ(t−s)
, k = 0,1,2L, λ > 0
白噪声过程 随机过程{ξt}t≥0称为白噪声过程,若Eξt=0,且 随机过程 ⋅εt−j ) = 0, j ≠ 0
显然,白噪声过程一个平稳的纯粹随机过程, 显然,白噪声过程一个平稳的纯粹随机过程,在金 融研究中主要用于模型无法解释的波动。 融研究中主要用于模型无法解释的波动。
3
二项分布是指随机变量满足概率分布
P (ξ = k ) = C p (1- q )
k n
k
n- k
其中, 其中,k=1,2, …,0<p<1,q=p-1。 , 。 二项过程实质上是将二项分布作为一个过程来描 述金融资产价格变化的。
4
假设股票价格在t时刻为 假设股票价格在 时刻为S(t),当时间变化到 时刻为 , t+∆t时,价格要么以概率 从S上涨到 上涨到uS(u >1), ∆ 时 价格要么以概率p从 上涨到 , 要么以概率q下降到 下降到dS(d<1);时间为 要么以概率 下降到 ;时间为t+2∆t时有 ∆ 时有 三种可能: 三种可能:u2S、udS、d2S,以此类推,见树型 、 、 ,以此类推, 结构
k! 则称{ξ 为参数为λ − 的 过程。 则称 ξt }t≥0为参数为λ(t−s)的Poission过程。 过程
直接计算可知, ξ 所以λ 直接计算可知,Eξt =Vξt =λt,即,所以λ表示单 ξ λ , 位时间内事件出现的平均次数,因而λ 位时间内事件出现的平均次数,因而λ也常被称为 发生率或强度。 发生率或强度。
最常见的随机过程或随机模型
1
主要内容
Brown运动或 运动或Wiener过程 运动或 过程 二项过程 Poission过程 过程 白噪声过程 自回归过程 移动平均过程 混合自回归移动平均过程 利率期限结构或均值回复模型 ARCH类模型 类模型
2
二项过程
1979年Cox、Ross和Rubinstein利用二项过程 年 、 和 利用二项过程 提出了二叉树期权定价模型, 提出了二叉树期权定价模型,用以构造股票价格运 动过程,进行股票期权定价分析。 动过程,进行股票期权定价分析。 目前, 目前,二叉树模型已被广泛应用于金融资产定价 领域, 领域,并为直观理解金融资产价格的复杂随机行为 提供了最佳认识工具, 提供了最佳认识工具,为金融计算提供了可行的数 值方法。 值方法。
⋅
⋅ ⋅
⋅
⋅ ⋅
5
显然, 显然,在t +∆t 时刻,股票的期望价格为 ∆ 时刻, E(St+∆t)=puS+(1-p)dS, ∆ , 在t +2∆t 时刻,股票的期望价格为: ∆ 时刻,股票的期望价格为:
,
E(St+2∆t ) = p u S + 2p(1− p)udS+ (1− p) d S
2 2 2 2
14
混合自回归—移动平均过程 混合自回归 移动平均过程 若时间序列{ξt }在t时刻,不仅与其以前的自身值 若时间序列 ξ 在 时刻, 时刻 有关, 有关,而且与以前时刻的冲击或扰动存在着一定的 依存关系,则称为混合自回归—移动平均过程 移动平均过程, 依存关系,则称为混合自回归 移动平均过程,其 一般形式(记作ARMA(m,n))为 一般形式(记作 , ) ξt =a1ξt -1+ a2ξt -2+…+ amξt -m+εt +b1εt -1+ ε b2εt -2+…+ bn εt –n
12
更一般地, 阶自回归过程 阶自回归过程{ξ 更一般地,m阶自回归过程 ξt }( t=1,2,…), , , 记为AR(m), 满足: 记为 ( ) 满足: ξt =a1ξt -1+ a2ξt -2+…+amξt -m+εt ε t=1,2,… , , m阶自回归过程具有 期记忆或者说 阶动态性。 阶自回归过程具有m期记忆或者说 阶动态性。 阶自回归过程具有 期记忆或者说m阶动态性 若滞后算子多项式1− 若滞后算子多项式 −a1z−…-amzm=0的根在单位 − 的根在单位 圆之外时,为平稳过程。否则,就是非平稳的。 圆之外时,为平稳过程。否则,就是非平稳的。
7
计数过程: 计数过程: 如果用ξ 表示[0,t]内随机事件发生的总数,则随机 内随机事件发生的总数, 如果用ξt表示 内随机事件发生的总数 过程{ξ 称为计数过程,且满足: 过程 ξt }t≥0称为计数过程,且满足: (a) ξt ≥ 0; ; 是整数值; (b) ξt是整数值; (c) 对于任意两个时刻 ≤ s<t,有ξs<ξt; 对于任意两个时刻0≤ 有 ξ (d) 对于任意两个时刻 ≤ s<t, ξt -ξs等于在区间 ( s, t ] 对于任意两个时刻0≤ ξ 中发生的事件的个数。 中发生的事件的个数。
15
利率期限结构或均值回复模型
在金融市场中,许多情况下的金融资产价格的变化, 在金融市场中,许多情况下的金融资产价格的变化,随着时 间的推移常常趋于某个长期平均水平,称为均值回复现象, 间的推移常常趋于某个长期平均水平,称为均值回复现象, 例如利率的变化就常常如此。 例如利率的变化就常常如此。具体的利率期限结构或均值回 复模型定义为
17
传统的随机过程和模型对金融资产收益的模拟和描述主要 是线性的,不能很好处理上述特征,因而也常常无法准确 是线性的,不能很好处理上述特征, 估计和预测金融资产的收益及其波动性。 估计和预测金融资产的收益及其波动性。 ARCH类模型一般由条件均值方程和条件方差方程两个方 类模型一般由条件均值方程和条件方差方程两个方 程组成。但由于此类方程主要用于估计波动性和相关性, 程组成。但由于此类方程主要用于估计波动性和相关性, 所以重点在条件方差方程,而条件均值方程常常比较简单. 所以重点在条件方差方程,而条件均值方程常常比较简单