随机振动讲义全文

随机振动讲义全文(总49页)

-CA^^M^HAT>(YICAI)-CompanyOnel

校作为文酗面，使用请直接删除

目录

第一章绪论.....................................错误!未指定书签。

1.1随机振动的基本概念和特征....................错误!未指定书签。

1.2随机振动研究的内容和意义...................................3

第二章随机振动的数学描述.........................错误!未指定书签。

2.1随机过程的基本概念和特征..................错误!未指定书签。

2.2随机过程的数学描述........................错误!未指定书签。

2.2.1随机变量定义.........................错误!未指定书签。

2.2.2一维随机变量的概率分布函数与概率密度函数错误！未指定书

签。

2.2.3多维随机变量.........................错误!未指定书签。

2.2.4随机变量的数字特征...................错误!未指定书签。

2.2.5随机变量的分布以及运算...............错误!未指定书签。

2.3随机过程的幅域描述........................错误!未指定书签。

2.3.1随机过程概率统计特征量...............错误!未指定书签。

2.3.2平稳随机过程.........................错误!未指定书签。

2.4随机过程的时域描述........................错误!未指定书签。

2.4.1各态历经随机过程.....................错误!未指定书签。

2.4.2平稳随机过程的自相关函数.............错误!未指定书签。

2.4.3互相关函数...........................错误!未指定书签。

2.5随机过程的频域描述：.......................................17

2.5.1典型函数的傅里叶变换.................错误!未指定书签。

2.5.2功率谱密度函数.......................错误!未指定书签。

2.5.3平稳随机过程的谱分类：...............错误!未指定书签。

2.5.4随机过程的分布.......................错误!未指定书签。

2.6随机过程的运算.............................错误!未指定书签。

2.6.1微分运算.............................错误!未指定书签。

2.6.2积分运算.............................错误!未指定书签。

2.6.3随机振动位移、速度和加速度的相关函数和谱密度函数关系.错

误!未指定书签。

第三章SDOF系统的随机响应......................错误!未指定书签。

3.1系统的脉冲响应函数和频率响应函数描述.....错误!未指定书签。

3.2单自由度系统随机响应分析.................错误!未指定书签。

第四章多自由度系统的随机响应分析..................错误!未指定书签。

4.1多自由度系统的脉冲响应函数、频率响应函数...错误!未指定书签。

4.2单输入问题的MDOF系统的随机响应...........错误!未指定书签。

4.3多输入问题的MDOF系统的随机响应...........错误!未指定书签。

4.4MDOF系统随机响应分析的模态方法...........错误!未指定书签。

4.5随机响应分析的虚拟激励方法...............错误!未指定书签。

第五章连续系统的随机响应分析......................错误!未指定书签。

参考文献...........................................错误!未指定书签。

第一章绪论

1.1随机振动的基本概念和特征

前面研究的振动问题都属于确定性振动(deterministicvibration),所谓

的确定性就是指振动是有一定规律的，或者可以用一个确定的函数来描述，或

者可以用若干离散的值来描述，而且这个规律是可以重复的，可以预先估计

的。例如，无阻尼自由振动问题：

inx+kx-O(1-1)

在确定的初始条件作用下，系统的振动响应规律为：

x«)=Asin(d+cr)(1-2)

其中，而，是由表征系统特性的物理参数确定的，A和a由初始条件确

定。只要已知初始时刻的振动值%,就可以预知之后任意时刻的振动值。

该系统在另外一次相同的初始激励下，系统振动规律理论上会得到完全的重

复。再看一个有外激励力作用的系统的振动规律：

nvc+kx=f(x)(1-3)

这个系统的振动规律为：

x(f)=/(1-4)

其中，/为任意的外激励，〃为系统的脉冲响应函数。这个杜哈梅积分如果可

以精确积分，振动规律可以表示成一个确定的函数表达式，如果不能，需要利

用数值积分，得到的振动规律是一组给定的离散时刻的确定的数值。同样，在

下一次相同的外激励作用下，振动规律还可以得到完全的重复。

在自然界和工程实际中还存在另外一种截然不同的现象，其变化是高度

不规则，无规律的，不可预估也不可重复，物理现象的这种变化规律称为随机

的。例如，海浪，地震，阵风(湍流)，火箭的喷气噪声以及不平路面。在随

机现象作用下，系统产生的振动规律也同样有随机的特征，振动过程是不确定

的，这样振动称为随机振动。工程中有很多这样的实际例子：

在海浪作用下，海洋平台结构、水面舰船、出入水的导弹的振动

在湍流作用下，飞行器结构的振动

在阵风作用下，高耸建筑物、桥梁的振动

在地震作用下，所有地面建筑结构的振动

在发动机喷气噪声以及大气气动噪声的作用下，火箭、导弹等飞行器结构

的振动

在不平路面的作用下，各种车辆的振动。

这些振动都是确定的工程结构在随机的外激励力或运动激励作用下产生

的，都是随机振动。上述例子共同的特征是：

激励和响应都不能用时间的确定函数来描述；

对于某一特定时刻取值不确定；

对于单个试验记录，从当前时刻的值无法预估之后时刻的值；

两次相同条件的试验结果不可能重复，但多次的试验结果放在一起却可以

发现现象的某些统计规律。

就是说振动运动是随机的，所以在任一给定时刻f=时X的精确值不可能

精确预计，我们最多只能求出在时刻"，X取值于某一区间的可能性或概率，

给出在某一时刻的统计规律，而且统计规律也可能是随时间变化的。

1.2随机振动研究的内容和意义

随机问题，主要分为两大类：

I）系统是确定性的，激励是随机的

前面所列举的例子都属于这一类。确定性的系统在随机的激励作用下，系

统的响应也是随机的。在这类问题中，主要研究激励以及由其引起的随机振动

响应的统计规律，研究这些规律与系统特性之间的关系。通常的随机振动研究

主要属于这一类。

2）系统是随机的，激励或确定，或随机

自然界和工程中也有这样的问题，例如，雨天，输电线的振动问题，这

里，输电线的质量是随机变化的，也就是系统的特性是随机的。这类问题，同

样也是研究随机现象的统计规律以及它们之间的相互关系。

当然，随机振动也有其它的分类，

按系统自由度可分为：单自由度随机振动；多自由度随机振动；无限多

自由度随机振动。

按振动微分方程的特点可分为：线性随机振动；非线性随机振动。

按随机振动频带宽窄可分为：宽带随机振动，窄带随机振动。

按振动的特性随时间变化情况可分为：平稳随机振动；非平稳随机振

动。

我们主要研究线性单、多自由度、连续体系统在单个和多个平稳随机激励作用

下的响应分析。

实际工程中，随机振动现象是十分普遍的，严格地说，一切实际系统的

振动都是随机的，只不过有些振动随机的成分很小，可以忽略，当作确定性系

统来研究。但是对于象湍流引起的飞机、火箭的振动、海浪导致出入水的导弹

的振动，以及前面介绍的其它例子，都必须考虑振动的随机性，用随机振动的

研究方法进行研究，才能得出更符合实际情况的结论。

第二章随机振动的数学描述

由于确定性的结构系统在随机变化的激励力作用下，系统的振动响应也

是随机变化的，所以随机振动主要研究激励以及由其引起的随机振动响应的统

计规律，以及这些规律与系统特性之间的关系。对这些规律我们可以利用概率

论的知识对他们进行定量或定性的研究，所以，首先我们要对随机激励或者随

机响应进行赋值，也就是用一个变量来表示，也就是要对随机振动的各个量进

行数学描述。

2.1随机过程的基本概念和特征

随机过程是对在空间和时间上高度不规则，事先无法预估，其变化也无

法重复，其统计规律随时间演化的物理现象的一种数学描述。工程中存在着很

多这种物理现象，如在第一章所举的例子，这些物理现象无法用确定性的理论

来描述，但可以用随机过程来描述。随机振动的数学抽象即为随机过程。

随机过程的每一次测量所得结果可看作一次实现，或叫样本函数。所有

可能的样本函数的集合构成一个随机过程。因此，随机过程是由时间上无限

长、样本的无限多个的样本函数构成的，可以写为：

X（D/GT,j=l,2,...}（2-1）

图2-1:随机过程示意图

随机过程的每次实现是一个确定的非随机函数，但各个实现各不相同，

因此为了得到随机过程的统计特性也必须做大量的独立测量。例如在同一条件

的海域内，布置n个同一类型的波高仪，可同时测得n个记录，得到n个实

现，不（。,9（。,•，演«）。在某一固定时刻4可得各样本瞬时波面高度

占⑥,％^）,…,怎,），它们构成了通常的随机变量x（4），在另一时刻弓又构成

另一个随机变量无&）。因此随机过程也可以是样本空间上的随机变量x（f）的集

合。下文就将X(r)表示为随机过程。随机过程是随机变量进一步发展得到的，

是随机变量随时间的变化，是随机变量的推广。

可以看出随机过程是对随机现象的完全描述，严格的随机过程应包含随机现

象的无穷多个独立测量样本，而且每个样本应该在时间上是无限长。实际分析

中，我们只能用样本长度有限，样本数目有限的样本集合来代替随机过程。所

得结果仅是随机现象统计特征的一个估计，一个近似。

2.2随机过程的数学描述

随机过程的概念一方面定义为无穷多个样本函数的集合，另一方面可以

看作无穷多个随机变量的集合

X&),i=l,2,...8(2-2)

其中X&)是由随机过程X在乙时刻所有可能的取值为«)构成的随机变量，j

是样本函数的编号，j=l,2,...ooo正因为它可以认为是由无穷多个随机变量构

成的，所以我们首先从随机变量的概率描述角度，来对随机过程进行描述。

2.2.1随机变量定义

对所研究的随机现象赋值便得到了一个随机变量，例如，哈尔滨地区每

年冬天的最低气温。在同一海域内布置n个同一类型的波高仪，在某一时刻所

测得的n个波高值，就构成一个描述波高可能取值的随机变量。在相同随机激

励的多次作用下，结构系统在某一固定时刻振动响应可能的取值，都属于随机

变量。

许多随机现象的试验结果表现为数量，用来表示随机试验各种结果的变

量叫做随机变量。随机试验的一种结果也就是随机变量的一个可能取值，这些

所有可能的取值的集合就是一个随机变量，用集合符号表示就是：

X={xj,/=1,2,3,...〃(2-3)

式中看为随机变量X的一种可能取值。〃取有限值就是离散随机变量，〃取无

穷大就是连续随机变量。

研究一个随机变量，不但要知道它在每次试验时的取值，更重要的是要

知道它取这个数值的概率。综上所述随机变量的基本特征，用数学的语言来描

述给出的定义为：定义于某样本空间C上的实变量如果对于每一

个实数X,X(〃)〈x的概率Prob{X(〃)《x}都存在，那么就称x(〃)为随机变

量。通常主要考虑随机变量X(〃)的值取在整个实数轴(70,”)上的问题。以下

为行文方便X(〃)简写为X。

对一个随机变量作完整的概率描述就是给出它的概率分布，也就是给出X取

值小于每一个xe(-oo,8)的概率，就是给出函数：

F(x)=Prob{X<x,x&(-00,00)}=P(X<x)(2-4)

F(x)称为X的概率分布函数。概率分布函数的性质：

1)F(+oo)=1(2-5)

由定义可知实变量X取值小于+8的概率是100%,或说X〈”是肯定的

2)F(-oo)=0(2-6)

X取值小于-8的概率是0,或说X<-oo是不可能的

3)/(x)是单调增函数

由定义可知,若巧>X”则尸小)〉尸(%)

4)尸(%)NO,恒非负

5)对任意元素有X取值在区间(巧,々)内的概率为：

F(x2)-F(Xj)=Prob{x]<x<x2)(2-7)

6)Pro仪X>x)=1-/⑺(2-8)

注意：对连续型随机变量，取值为一个特定值的概率为零，

ProZ?(X=x)=0o

当F(x)连续可导时，可以得到其导数函数

/、dF(x)F(x+dx)-F(x)

〃(x)--------lim--------------(2-8)

dxdiodx

其意义可解释为随机变量X取值在x附近的单位区间的概率大小，因为：

因此，p(x)大表示网光)在该点的变化较大，也就是在这个区间概率分布密度也

大，所以也称p(x)为概率分布密度函数，简称概率密度函数。概率密度函数表

示X取值在x点附近的单位区间内的概率大小。

概率密度函数的性质：

1)jp(u)du=F(x)(2-9)

—00

oc

2)jp(u)du—F(-i-oo)-F(-oo)=F(+oo)=1(2-10)

—00

x2x2

3)jp(u)du-dF(u)-F(x2)-F(xt)(2-11)

w再

4)p(x)NO(2-12)

单调增函数的导函数恒非负。

5)〃(+8)="(-oc)=0(2-13)

有些问题需要考虑两个或两个以上的随机现象同时发生的概率，例如打靶，

就需要考虑在刘y两个方向同时射中区间的概率，这就是

二维联合概率问题，还有更多维，仅以二维为例。

对于二维的随机变量Z=[X,Y],它的联合概率分布函数定义为：

F(x,y)=Prob(X<x,Y<y)=P(X<x,Y<y)(2-14)

即F(x,y)为随机变量X取小于x同时丫小于y的概率，性质：

1)F(x,y)>0,x,y&R(2-15)

2)F(-oo,y)-F(x,-oo)=F(-oo,-oo)=0(2-16)

3)尸(g”)=l(2-17)

4)F(x,+oo)=Proh{X<x,Y<+oo)=F(x)-Prob(X<x)(2-18)

5)F(+oo,y)-Prob(X<+oo,y<y)-F(y)-Prob(X<y)(2-19)

6)F(x,y)单独对是单调增函数

x2

:

7)F(x2,+oo)-F(x,,+oo)=Prob(x]<X<x2,Y<-H»)=F(X2)-/'(X1)=jp(x)dx(2-20)

当尸(x,y)有二阶偏导数时，有

夕尸(x,y)

p(x,y)=(2-21)

dxdy

这个二阶偏导函数定义了二维联合概率密度函数。由定义及F(x,y)的性质可

知，

xy

尸(x,y)=JJp&Mdn(2-22)

-O0-00

二维联合概率密度函数性质：

1)p(x,y)>0(2-23)

-KC+sC

2)JJp(x9y)dxdy=1(2-24)

―30—00

X+00X400x

3)尸(x,”)=jj4d77=J(Jp(4,〃W"W&=F(x)=]p记)

-00-CO-00-00-00

所以有

+50

Jp(x,yWy=p(x)(2-25)

—00

同理，由于

有

-KC

Jp(x,y)dx=p(y)(2-26)

-oc

这就给出了二维联合概率密度函数与一维的关系。

对于二维随机变量，还定义有条件概率密度函数为：

<、p(x,y)/、p(x,y)

p(x:y)='~,p(y:x)='-,

P(y)，(x)

其中p(x:y)表示在y条件下，x发生的概率，且有

p(x,y)=p(x:y)-p(y)=p(y:x)•p(x)(2-27)

若X,丫统计独立，则

p(x:y)=p(x),p(y:x)=p(y)(2-28)

且有

p(x,y)=p(x)-p(y)(2-29)

随机变量的统计特征可以用概率分布函数，或概率密度函数作完整描述，

但要确定这些函数一般不大容易，通常也不是总有这个必要，实际问题是只需

主要的统计特征即可，这些主要的数字特征称为随机变量的矩。

原点矩：实随机变量X的n阶矩定义为X"的集合平均，也称n阶原点矩，

即有

8

£[%"]=Jx"p(x)dx(2-30)

—00

其中最常用的是一阶原点矩和二阶原点矩。

一阶原点矩定义为

00

E[X]=Jxp(x)dx(2-31)

—oo

也就是随机变量的均值，也称数学期望，常记为从。(对离散随机变量有

〃,=仇幻=才王〃(刈，如果随机试验得到一系列独立的观测值演

/=1

(7=1,2,3…〃)，那么其样本均值为：

n,=i

一阶原点矩性质：

l.E(a)=a”是常数(2-31)

2.E[aX]^aE[X](2-32)

3.E[a+X]^a+E[X](2-33)

4.E[X+Y]=E[X]+E[Y],或者典(2-33)

证明：

E[X+Y]^[f(x+y)dxdy=[fxp(x,y)dxdy+[fyp(x,y)dxdy

J-00J—ooJ-OOJ-00J-00J-a>

p+<»p-Hx>/»+00「+00

=1x\p(x,y)dy-dx+yp(x,y)dxdy

J-COJ-COJ-COJ—00

广+8,+00

=J-coxp(x)dxJ+-x>yp(yWy

=E[X]+E[y]

5.若二者相互统计独立

E[xr]=E[xiE[r]或者E[nx』=naxj(2-34)

ii

证明：

二阶原点矩定义为：

E[X2]=x2p(x)-dx(2-35)

J-00

也称为随机变量的均方值，常记为上;，通常表示随机变量的能量水平。

上面讨论的都是随机变量相对于坐标原点的矩，也称为原点矩，还有一种常

见的矩，是相对于均值的，称为中心矩。〃阶中心矩定义为：

EKX-4力=公(2-36)

一阶中心矩为：

+CO

<•+□0广广+oo

I(x-jLix)p(x)dx=Ixp(^x)dx-/Lt\p{x}dx=)nx-/nx=0(2-37)

二阶中心矩为：

22

E[(X-E[X])]=r(X-AV)P(X)^=D[x](2-38)

J-oo

也称为X的方差，常记为b；,其平方根4称为标准差。

对离散随机变量有

2

D[X]=戊=£[(X-//.V)]=之(七—4)2p(x,)(2-39)

i=l

样本方差(Samplevariance)

D[X]=a-=-X(^-^)2(2-40)

方差表明随机变量偏离均值的程度。方差性质：

LO[a]=0,a是常数(2-41)

2.D[aX]=a2D[X](2-42)

3.D[a+X]^D[X](2-43)

4.D[aX+bY]=a2D[X]+b2D[Y],若x,y统计独立(2-44)

证明:

均值，均方值(均方根值)，方差4(标准差)是随机变量最重要的三个数字

特征量，它们之间有如下关系:

22

<r；=E[(x-4力=E[X-2X^X+^]=E[X]-2^XE[X]+^=/一4(2-45)

联合矩：多个随机变量的矩的关系是联合矩，以两个随机变量X,Y为例，其

(九,加)阶的联合原点矩定义为:

E[X'P"]=「'「"x"y'np(x,y)dxdy(2-46)

J-00J—co

〃=m=1时有

E[XY]=rrxyp(x,y)dxdy,也称为相关矩(2-47)

J-ocJ-00

当X,Y统计独立时有，E[XY]=E[X]E[Y]0

同理有(〃,加)阶的联合中心矩定义为：

E[(X————渥y(2-48)

〃="=1时有

仇(X-4V)]=J：J：(%_〃,.)(y_〃v)p(x,y)dxdy

=「JI(孙一切"一)%+N*、)P(x,y)dxdy(2-49)

^E[XY]-E[X]E[Y]

也称为随机变量X,Y的协方差(Covariance)，两个随机变量之间的协方差表征

了它们之间的相关性，通常用J,=Cov[X,y]表示，即

J,=Cov[X,Y]=E[(X-〃,)(y-〃v)]=E[XY]-E[X][Y](2-50)

当两个随机变量相互统计独立则有

E[XW=E[X][W,q.=0(2-51)

当c»,不等于o时，说明x,y之间具有相关性，但是相关程度的大小，通常用

。灯的无量纲化的系数来表征

内晤卜1(2-52)

称为相关系数。其绝对值小于一，为了证明这一点，利用如下着名的Schwarz

不等式

2

E[|xy|]<(E[X2]E[Y2]y(2-53)

特别地，当y=i时,有印印<(比闩)；

当q=°，即x,y统计独立时有%=o,所以

0<|A?.|<l,-l<Ay<l,(2-54)

当|Q』=I时，称为随机变量X,Y完全相关。

2.2.5随机变量的分布以及运算

随机变量的特定概率密度函数对应着特定的取值分布，常见的分布有均匀

分布，高斯分布(正态分布)等。

均匀分布的概率密度函数为

'1,

----,a<x<b

p(x)={b—a(2-55)

.0,其它

高斯分布的概率密度函数为

[(A"*)2

=--e2。；(2-56)

2和1

随机变量的初等函数仍然是随机变量，后者的分布由前者确定，且若已知

x的p(x),y=g(x),则有

E[Y]=ryp(y)dy=「'p(x)•gMdx(2-57)

J-30J-00

2.3随机过程的幅域描述

2.3.1随机过程概率统计特征量

上述对随机变量的成熟的概率描述手段，可以直接用于描述随机过程，只不

过为了表示随机过程是一个动态的，随时间变化的过程，需要加一个时间变

量，如p(x，4)表示随机过程在小寸刻的随机变量x&)的概率密度函数，一维概

率分布函数定义为：

F(x,r)=Prob[X(/)<x,x€/?,rGT]=fp(x,t)cbc(2-58)

J-00

对应的数字统计特征为：

4(f)=E[X⑺]=匚xp(x,t}dx(2-59)

/⑺=E[X2(t)]=匚x2p(x,t)dx(2-60)

b；(/)=E[(X⑺-4⑺力=村⑺-%⑺(2-61)

表明随机过程在每一时间截口的分布中心，能量水平和偏离分布中心的程度。

这些一维的概率分布只能描述各个独立时刻单个随机变量的概率特性，无法揭

示随机过程不同时刻之间的相互关系，为此必须使用二维以上的概率分布描

述。随机过程的二维概率分布函数定义为：

/(x”fI,%2,)=Pr帅[X(G<玉，X«2)<%2]=/J：P(x"，尤2，f2al公2Q-62)

其性质也和前述二维概率分布函数和二维概率密度函数性质类似。回忆前述描

述不同随机变量之间相关程度的数学特征量是协方差，对随机过程不同时刻之

间的相关性也可以用该量来描述，同样定义：

CovtX(Q,X&)]=£[(X&)-4&))(X«2)—〃M))](2-63)

为随机变量的自协方差，通常用6(44)表示。

4小幻=CoWXQJ.XG)]=a(X6)—4(G)(X«2)—〃&))]

(2-64)

=E[(X(GX«2

上式右侧第一项是X(G，X«2)的相关矩，一阶联合原点矩也称随机过程X")的

自相关函数，通常记为：

&.((I,弓)=E[X&)X(t2)]=II玉々P(X1，4，々，，2)的公2(2-65)

上述(2-64)公式表明若随机过程的均值4⑺=0,那么有

。«为)=R«4)(2-66)

R,(小功也就表示了随机过程不同时刻的随机变量之间相关程度。由于多数随机

过程，例如，海浪符合这个条件，所以，将二者统称为相关函数。用凡代替

Cto很显然有

-=E[X(t)X(t)]=^(t)(2-67)

Cxg)=£[(X(Z)-4⑺力=戊⑺(2-68)

以上考虑的是单一随机过程的概率描述。对不同的随机过程X⑺,丫⑺可分别派

生出两族随机变量x(G，y&)(i#=i,2,3...)。因而，有是需要考虑它们之间的联

合概率分布或联合矩。此时联合概率密度函数可以写为PG/,%,%)。他们之

间的二阶联合原点矩和中心矩分别为

r

(i,t2)=E[X(zt)r(r2)]=xyp(xt,ti，y2,t2)dxdy(2-69)

q,(22)=E[(x&)—〃M))(y«2)—〃>«2))]=仆(22)-〃«)巴,伉)(2-70)

Rn.,C.分别是称为互相关函数和互协方差函数，表示他们是来自于不同的随

机过程，对应的来自于同一随机过程都冠以“自”。

均方差，方差，自相关，协方差，统称为二阶矩。

若E[X2(r)]<oo,则均方差存在，由Schwarz不等式：

可以推知自相关函数必定存在。即可认为随机过程的二阶矩函数存在，X(f)表

示二阶矩过程。

与相关系数对应规范化的互协方差函数为:

,.(?],t)

2(2-71)

2.3.2平稳随机过程

在实际中经常遇到这样一类随机过程，他们随时间变化是在一平均值周围连

续地随机波动，其统计特征都基本上不随时间变化，称该过程为平稳随机过程

(Stationaryrandomprocess)。

平稳随机过程一般定义:若一个随机过程的概率特征量在时间参数做任意平

移时保持不变，则称此过程是平稳的。

严格平稳随机过程定义:若随机过程的〃维联合概率密度函数对任意实数7

都有

p{xx,t\,X2,t2,--Xn,tn)=p{Xx,t{+T,X2,t2+T,Xn,tn+T)(2-71)

则称此过程是〃阶平稳的，且低于〃的各阶也都是平稳的，如

这个定义是严格平稳的条件。严格平稳的条件工程上很难满足。因此引入了广

义平稳(弱平稳或者宽平稳)的概念：若一个随机过程均值和自相关函数或者

协方差不随时间变化，即满足

1.4(4)=4=const(2-72)

2.QCr,+r)=CV(/2,Z2+r)=...=Cx(r)-const(2-73)

两个条件，即均值不随时间变化，协方差也不与计时起点或时间原点有关，只

与时差7有关。这样的随机过程称为广义平稳随机过程。工程中的平稳的含义

通常是指广义平稳。平稳随机过程的协方差

协方差的一个重要性质是:在随机过程上增加一个确定性函数并不改变协方差函

数。例如：x（f）的均值为从⑺和协方差GG/）,〃⑺是一个确定函数，则

yQ）=x（t）+〃⑺的协方差不变。

显然：当〃⑺=-〃,“）时有，y（t）=x（f）-z（f）,且丫⑺的均值为零，

。,、（4,，2）=Cr（，l12）0

所以对协方差的要求就和对自相关函数的要求一样。此外，对平稳随机过

程而言，有时为了简化运算而假设平稳随机过程均值为零，工程中有许多过程

为零。

注意：由上述平稳随机过程定义可知，满足这个定义的随机过程的样本函

数无限长，而且在整个（0,+00）上统计特性对时间参数原点的选取有一定的均匀

性，即与参数f的初始时刻选取无关，而实际的随机过程通常很难满足这个条

件，因此在实际工程问题处理中，只要一个随机过程在一个较长的区间上呈现

上述均匀性，就可以近似看作平稳随机过程。例如，火车在启动和停止阶段，

就不满足均匀性的假设，但在中间较长一段时间内是基本匀速行驶的，因此可

看作广义平稳过程。

2.4随机过程的时域描述

随机振动的时域描述主要指时差域描述，用随机过程不同时刻之间的相关

情况来描述随机振动。这里主要指平稳随机过程，而且通常还假设均值为零。

2.4.1各态历经随机过程

平稳随机过程的均值和方差不依赖于时间，均值可由任意时刻的多个样本

的集合平均求得，协方差也仅取决于作相关的时差「，但仍需对随机过程进行

大量观测，取得足够多的样本函数，尽管样本函数可能不需要很长，但工作量

仍然是很大的。因此就猜想能否用仅用一个足够长的样本来代替大量样本构成

的总体，用该样本的时间平均特性代替样本空间的集合平均特性呢？为此引入

样本函数时间平均概念。

设平稳随机过程X⑺任一样本函数为不⑺，下文为书写简便用X⑺代替任一

无限长样本函数，其时间均值定义为：

T

4=/〈X（f）\）=..1C，90（2-74）

时间平均意义下的自相关函数定义为:

1y

??A.(r)=(X(Z)X(r+r))=lim—£x(r)x(r+r)dt(2-75)

时间平均意义下的均方值当r=0时有，

"(0)=(X«)X(f»=M=limi[，⑺力(2-76)

时间平均意义下的方差定义为

珑—4)2"㈣"：(x⑺—〃,.)2山(2-77)

各态历经随机过程：对一个平稳随机过程，若有

E[X(/)]=(X(f))=4(2-78)

则称该平稳随机过程关于均值遍历。若有

(X(f)X(f+7))=E[X(r)X(r+r)]=/?v(r)(2-79)

则称过程关于相关函数具有遍历性。具有一定遍历性的随机过程称为遍历过

程，或称各态历经随机过程。也可以写成如下形式：

(Xj(/»=£[X(。)](均值遍历)(2-80)

(X«)X,C+T))=E[X(O)X/+T)](相关函数遍历)(2-81)

其中，，为样本函数编号，/为时间采样点编号。

平稳随机过程遍历的基本含义就是样本函数的总体统计特征等于单个样本在

较长时间段内的时间统计特征。

2.4.2平稳随机过程的自相关函数

根据前述的集合平均意义以及时间平均意义上的自相关函数定义，可以得

到其性质如下：

1T

1.R、(0)=&Hx(r)M)•力=际20(2-82)

2.|/?v(r)|</?v(0)(2-83)

证明：

由于

所以

由此有

说明随机变量与自身的相关性最好。

3.&(7)=&(-T)(2-84)

证明：&(r)=limyJ；x(t)x(t+T)dt

令，t—+di'

所以有

由平稳性定义也可以直接得到&(T)是偶函数这个性质。

4.lim/?v(r)=O(2-85)

r—>oo

通常实际的物理系统总是有一点耗散的，随着时差的增大，一般来说随机

过程的相关性有所减弱，而且当「到8时有，尺(7)趋向于0。

2.4.3互相关函数

在随机振动分析中，通常要用到来自两个不同随机过程的相关，例如随机激

励力与随机振动响应的相关情况，还有两个以上不同的随机激励力作用在同一

结构上等情况。对各态历经的随机过程XQ),y(r),互相关函数定义为：

=lim—1[T+(2-86)

,X->00TJ。

性质：

一般不对称(2-87)

2.勺⑺=5”(—7)(2-88)

f

证明：7?vv(r)=limy£x(f)y{t+r)dt

令t+T=t',dt=df,t=t'-r

例2-1：%(。与马⑺为两个平稳随机过程，求=。内”)+”2工2⑺自相关函数

解：

对均值为零的平稳随机过程，若相互独立则有q,,(7)=0,即&V(T)=0

2.5随机过程的频域描述：

2.5.1典型函数的傅里叶变换

x(f)的连续傅里叶定义为：

X(0)=「％(/)•⑪.出(2-89)

J-00J-O0

x(t)=——「'X⑼•e^dco(2-90)

2万卜

线性性质：

以⑺+by(t)=aX(m+Z?Y(co)(2-91)

对称性质：

xQ)0X(co)

(2-92)

x(-t)。X(一⑼

平移性质：

必±幻0乂(0)".

(2-93)

X(<y±«y0)<=>尤(f)e'

变标尺性

八、1一一CD、

x(kt)<=>—X(—)

kk(2-94)

X(kty)o：x(()

共辗性

X*(o)=X(-0)(2-95)

微分特性

北。)=icoX(CO)

.(2-96)

无⑺o-"X(0)

乘积与卷积特性

%⑺•W⑺。X|(⑼*x2(<y)

%(①)•X2(0)oXj⑺*x2(t)

典型函数的傅里叶变换

1.脉冲函数

定义：

若有「'b«).左=1,称为单位脉冲函数，其性质为

J-00

傅里叶变换为

/⑼=♦[♦(/)]=r必)•e-jMdt=「'8(t-0)-eiMdt=e"°=1(2-98)

J-30J-CO

J(r)=Fl[<y(«)]=—f，Z3(co)-eiMd(o=eiMd(o(2-99)

2兀-

可以得出如下结论：

+00

JeiMd(o=1n8(t}(2-100)

-00

2.正余弦函数sin69(j,cosG(J的傅里叶变换

j(i(v

F[cosa)Qt]=gF[e^+e~]=g[2而侬-0)+2砌-例)一①)]=乃[3®+跳)+3(①-①0)](2-101)

1.1

+J(l)bt

F[sin(Dot]=—:F[e-—；[2•(g-幻)一2乃5(一/一⑼]

2j2/(2-102)

=-5(°-跳)一53+%)]=%胆(。+%)*3一⑻]

j

3.单位指数函数—"〃⑺«>0

F[e-a,u(t)]=「"e-a'u(t)-e-iMdt=fe"e"dr=f"叱如你=―-—(2-103)

J-()oa+jco

4.矩形脉冲函数

X(<y)=「"加)9的力=/a-e-j,0,dt=^-(e-jaT-ejaT)=-^(-2;sin(oT)=—sin(oT(2-104)

JFJ-R(0(00)

若TfO,则矩形函数相应于矩形脉冲，则有

X(a))=—coT=2aT=a-2T,为矩形脉冲T的面积。

co

252功率谱密度函数

相关函数的傅里叶变换称为功率谱密度函数(Powerspectraldensity

function),自相关函数的傅里叶变换称为自功率谱密度函数，互相关函数的傅

里叶变换称为互功率谱密度函数。分别叙述：

自功率谱密度函数

定义1：自相关函数的傅里叶变换

S、.(3)=[Rv⑺6-师公或(2-105)

(2-106)

也可以说自相关函数是自功率谱密度函数的逆傅里叶变换，即

&e)=(1'S,kco)ei<aTd(o或

(2-107)

R,(r)=J：S,(/iW(2-108)

由于6(0)表示均方值，因此上式当7=0时有

6(0)=/=J：S")4(2-109)

所以S,(7)在整个频带上的积分等于它的均方值，可以说，S,(7)表示x⑺在单

位带宽内具有的能量，具有能量(或功率)的密度的概念，所以称为功率谱密

度。

所以也有如下的定义：

1

92

Sv(/)=lim-|x(/)|,(2-110)

r->8]11

可以证明这两个定义是等价的。

证明：

同时有

由于对任意的随机函数X"),上两式均成立，因此有：

自功率谱性质：

表示振动功率按频率的分布

1.5,3)20(2-111)

2.4(0)=⑼4(2-112)

所以S,表示单位频带上信号的能量

3.自功率谱是偶函数

f+OO

=SK(a>)=2£7?v(r)cos<yrJr(2-113)

互功率谱密度函数：

对应的互功率谱也有两个等价定义

1.，(⑼=匚RqW4或Sxy(f)=匚R,,(7)e-J2"r.八(2-114)

2.S")=lim"/W)(2-115)

rT-X»/

互功率谱密度函数性质：

互功率谱密度函数一般是复数，不对称，且有

"3—)=黑3)(2-116)

证明：

对于实际的信号，一般没有负频率的概念，前述的意义是在(F,”)上，这仅

仅是理论上的定义，因此工程上便于应用，把负频率的谱密度折算到正的频率

上去，由S,(⑼是偶函数，所以定义:

2S(co')0<«<+00

(2-117)

八［0(y<0

移为单边自谱密度函数，对应S*⑼称双边自谱密度函数。

广8

GV(69)=2SA.(69)=4£Rx(T)coscordr(co>0)(2-118)

单边自谱下的面积同样等于均方值，因为：

+

=——jSx(a))dco=——£2Sx{co)dco=——jGx(co)dco(2-119)

类似的定义单边上的互谱密函数：

2S,,(⑼69>0

Gy3)=<(2-120)

0CO<0

对应的称为双边互谱密度函数。

相干函数：

在时域内用相关系数表示两个随机变量的相关程度，同样在频域内也定义一个

类似的无量纲数来表示随机函数的相关程度。

2,、1%(叫

(2-121)

S、.(0)Sv(。)

可以证明:

片，3)wi(2-122)

相干函数可以用来检查系统是否有随机干扰和非线性干扰，即如果片接近于

1,表示所经过的系统非线性程度很小，噪声干扰也很小，反之干扰比较大，得

到的谱密度函数不可信，因为输出y不完全是由输入x引起的。一