平稳过程,基本理论是在20世纪30~40年代建立和发展起来的,并已相当完善。其后的研究主要是向某些特殊类型以及多维平稳过程、平稳广义过程和齐次随机场等方面发展。平稳过程理论在无线电技术和自动控制等领域有着广泛的应用,并且是诸如时间序列分析、信号分析、滤波、预测理论以及控制理论等应用学科的重要工具。
正文
设 X =( X ( t ), t ∈ T )是一个取 复数 值的随机过程,其中指标集 T 为整数或 实数 全体(分别称为离散指标和连续指标)。如果对任意的自然数 n 及任意的 t 1, t 2,… tn , 的概率分布与( X ( t 1), X ( t 2),…, X ( t n))的概率分布相同,则称 X 为严平稳过程。如果二阶绝对矩 ,而且对任意的 t ,τ∈ T ,均值(见数学期望)E X ( t ) m( 常数 ),协方差 (与τ无关),则称 X 为宽平稳过程。Г( t )称为 X 的协方差函数。一个严平稳过程,如果它的二阶矩有穷,则一定也是宽平稳的(见矩)。
宽平稳过程的谱分解 将傅里叶分析方法应用于宽平稳过程,可以把过程表成有不相关随机振幅的简谐振动的叠加,这就产生了过程按频率的谱分解,它是宽平稳过程理论的一个基本结果。
宽平稳过程的协方差函数是非负定的,即对任意 t 1, t 2,… t n∈ T ,任意 复数 z 1, z 2,… zn ,都有 。根据这一性质,对于离散指标,Г( t )和过程本身分别有如下的谱分解式:
。
对于连续指标,如果还假设过程为均方连续的,即对任意 (它的一个等价条件是Г( t )为连续函数),则有谱分解式:
,
上面的 F ( λ )是区间【- , 】或(-∞,∞)上的有界非降函数,并可取为右连续的,称为过程 X 的谱分布函数。η( λ )是相应区间上的正交增量过程,即满足Eη( λ )呏0,且当 时,有
及
。
对η 的积分是随机积分中一种常见的对正交增量过程的积分。如果 F ( λ )绝对连续,则称它的导函数 ƒ ( λ )= F ┡( λ )为过程 X 的谱密度(或功率谱密度)。最重要的特例是有理谱密度,对于离散指标,它的一般形式是
,
式中 多项式 与 当| z |≤1时不为0。这时过程 X 称为自回归滑动平均序列,通常简记成ARMA( p , q )序列,它满足如下的随机差分方程:
式中( ξ ( t ), t ∈ T )为正交随机序列,即满足E ξ ( t )呏0,当 t ≠ s 时, ,且 。对于连续指标,有理谱密度的一般形式是
,
式中 q < p ,且多项式 与 当 z 有非负实部时不为0。从物理上看,上述性质表明,具有有理谱密度的平稳过程可以看成白噪声(即具有常值谱密度的平稳过程)输入一个有限阶非时变线性系统所得的输出。这一特性使得它们在实际应用中占有重要的地位。
从平稳过程的谱分解可以推出一些重要的结论。例如均方大数律 (见大数律):如果宽平稳过程 X 的谱分布函数 F ( λ )在 λ =0处连续,则在均方收敛的意义下,成立
,
。
后一积分是均方随机积分。
宽平稳过程的线性预测 这是由Α.Η.柯尔莫哥洛夫和N.维纳在1940年左右提出并解决的问题。对于一个均值为0的宽平稳过程 X =( X ( t ), t ∈ T ),随机变量集{ X ( s ), β ∈ T , s ≤ t }表示到时刻 t 为止所能观测到的过程的历史,用M t ( X )表示由这些随机变量的一切有限线性组合及其均方极限所构成的希尔伯特空间。设τ>0,所谓"τ步"线性预测,就是要用M t ( X )中的随机变量来估计还未曾观测的 X ( t +τ)( t +τ∈ T );而线性最优预测,就是要在M t ( X )中选择 ( t ,τ),使预测误差的方差
达到最小。可以证明 ( t ,τ)惟一存在,它就是 X ( t +τ)在空间 上的投影。如果对任意的 t , β ∈ T ,都有 =M s ( X ),则称过程为决定的或奇异的,这时 ,即可以无差误地进行预测;如果所有M t ( X )的公共部分仅包含零,即 ,则称过程为纯非决定的或正则的。一般的宽平稳过程 X 可以分解成相互正交的两部分之和: ,其中 X r是纯非决定的宽平稳过程, X s是决定的宽平稳过程。这时, X r可以按离散指标或连续指标而分别表为一个正交随机序列( ξ ( t ), t ∈ T )的向后的滑动和:
,
或一个正交增量过程 ( ξ ( t ), t ∈ T ) 的向后的滑动积分(对正交增量过程的积分):
,
而且对任何 t ∈ T 有 。宽平稳过程的这一分解称为沃尔德分解。由此可得,对于离散指标,
,
;
对于连续指标,
,
。
此外,宽平稳过程 X 本身为纯非决定的充分必要条件是,它有谱密度 ƒ ( λ ),而且满足
,
。
具有有理谱密度的过程是纯非决定过程的重要特例。
多维宽平稳过程 设 X ( t )=( X 1( t ), X 2( t ),…, Xk ( t ))┡是由 k 个分量组成的,而且均值E X ( t )呏m为常值向量,协方差阵E( X ( t +τ)-m)( X (τ)-m)*=Г( t )与τ无关,则称 X =( X ( t ), t ∈ T )为 k 维宽平稳过程,其中记号“┡”与“*”分别表示向量或矩阵的转置与共轭转置。这时,Г( t )与 X ( t )在形式上有和一维情形一样的谱分解,只是 F ( λ )(或 ƒ ( λ ))变为 k × k 的函数矩阵,它的对角线元就是各分量本身的谱分布(或密度)函数,而它的非对角线元称为相应分量的互谱分布(或密度)函数,而且 ƒ ( λ )= ƒ *( λ )是非负定矩阵。关于多维宽平稳过程的线性预测问题,也有类似于一维的结果。
齐次随机场 如果随机过程 X =( X ( t ), t ∈ T )的指标集 T 是 k 维整值向量或实值向量的全体,且其均值与协方差函数满足与宽平稳过程的定义相同的条件,则称 X 为齐次随机场或 k 指标平稳过程。这时 X ( t )与Г( t )也有相应的谱分解。如果进一步,Г( t )只与指标向量 t 的长度 ( tj 为 t 的分量)有关,则称 X 为迷向场。齐次场在力学的湍流理论中很有用。齐次场的线性预测问题比宽平稳过程的情形要复杂得多,中国学者 江泽培 开始了这方面的工作。
严平稳过程的遍历 定理 关于严平稳过程,最重要的性质是以概率1成立的遍历定理,或称为各态历经定理。这一名词来自物理学。任何严平稳过程 X =( X ( t ), t ∈ T )都可看作是由某个概率空间 ( ,F, P )上的保测变换群{ S t, t ∈ T }作用于随机变量 X (0)而产生的,其中对任意的 , 实数 x 及随机变量 ξ ,满足条件 ,而 X ( t )= S t X (0)。用 表示使 St A = A 对任何 t ∈ T 成立的事件 A 的全体,则 是F的子 σ 域。如果均值E X (0)存在,则 以概率1收敛于条件期望E{ X (0)| }。如果 中只含概率为1或为0的事件,则称过程 X 为遍历的,这时 以概率1收敛于E X (0)=E X ( t )。后一结果也称为严平稳过程的强大数律,它表明,过程几乎所有的样本对时间的平均都趋近于每一时刻的过程值对概率分布的平均。
参考书目
J.L.Doob, sto chastic Processes, John Wiley & Sons,New York, 1953.
E.J.Hannan, Mul tip l e 时间 Series,John Wiley &Sons,New York, 1970.