随机振动完整课件

1、第五章 随机振动,5-1 引言,5-2 随机过程,5-3 随机过程的数字特征,5-4 相关函数,5-5 功率谱密度函数,5-6 线性系统在随机激励下的响应,1,学习交流PPT,5-1 引言,前面三章所考察的振动都是确定性振动振动系统的规律可以用时间的确定性函数来描述 振动系统的物理量可以用随时间变化的确定函数来描述，因此，确定性振动中的物理量在将来某一时刻的值是可以预测的。 举例：单自由度系统的简谐强迫振动,2,学习交流PPT,自然界和工程上还存在着另一类振动，它们的规律不能用时间的确定函数来描述，但又具有一定的统计规律性。 在数学上这类振动可以用随机过程来加以描述。 这类振动称为随机振动非确

2、定而又具有统计规律。 举例：车辆行驶时由于路面不平引起的振动,3,学习交流PPT,确定性系统+确定性激励 确定性响应 确定性系统+随机激励 随机响应 随机系统+任何激励 随机响应,4,学习交流PPT,随机振动与确定性振动的本质区别在于它一般指的不是单个现象，而是一个包含着大量现象的集合；从集合中的单个现象来看似乎是杂乱的，但从总体来看却存在着一定的统计规律性。因此，它虽然不能用时间的确定函数来描述，但能用统计特性来描述。 在确定性振动中，系统的激励与响应之间有着确定的函数关系，而在随机振动中，只能满足于确定它们的统计特性之间的关系。,5,学习交流PPT, 介绍随机振动中物理量的描述方法（相关函

3、数、功率谱密度）。 讨论受随机激励的振动系统的激励、系统特性、响应三者统计规律性之间的关系。,本章内容：,6,学习交流PPT,过程：物理量随时间变化的情况 随机过程：无法准确预知物理量随时间的变化情 况，但其变化规律服从统计规律 随机过程是大量现象的一个数学抽象，理论上是由无限多个无限长的样本组成的集合。 ：样本函数 对于随机现象，我们感兴趣的往往不是各个样本本身，而是力图从这些样本得出总体的统计特性。,5-2 随机过程,7,学习交流PPT, 随机过程 的所有样本函数 在时刻 的值 构成一个随机变量 对随机变量求集合平均 称为随机过程 在时刻 的集合平均值 。 一般情况下 依赖于采样时刻 。,

4、一、集合平均 . 平稳过程,5-3 随机过程的数字特征,8,学习交流PPT,、 构成两个随机变量 对它们的乘积 求集合平均 称为随机过程 于时刻 与 的自相关函数。它是时差 的函数，在一般情况下，它也依赖于采样时刻 ，反映这两个时刻的随机变量 的统计联系。,自相关函数,9,学习交流PPT,随机过程可以根据其统计特性是否随采样时刻而变化来进行分类。 统计特性依赖于采样时刻的过程非平稳过程 统计特性不依赖于采样时刻的过程平稳过程,平稳过程,10,学习交流PPT,平稳过程的特点,集合平均值为常数 相关函数仅仅依赖于时差,11,学习交流PPT,二、时间平均.各态历经过程,随机过程的每一个样本函数可以在

5、时域内求得: 时间平均值 对第k个样本求平均 补充：求函数的平均值,12,学习交流PPT, 自相关函数,在一般情况下，对于不同的样本将得到不同的 和,13,学习交流PPT,各态历经过程,定义：各个样本的时域统计值都是等同的，而且任一个样本函数在时域的统计值与任一时刻的随机变量的统计值相等。,14,学习交流PPT,结论：,集合平均与时刻t1无关，而时间平均与样本标号k无关，是过程为各态历经的充要条件。 从任何一个样本得出的时间特性就等于集合平均特性，将使数据处理容易的多。 举例：半主动悬挂特性评价,各态历经过程一定是平稳的，反之不一定。 平稳+各个样本的统计特性相同 各态历经,15,学习交流PP

6、T,例5-1：求正弦函数的相关函数,物理意义：表示样本函数 与其延时 时刻得到的 之间波形的相似程度。,，相似程度最高。 ，相似程度最低。,16,学习交流PPT,各个样本函数在时刻 t1的值 构成一个随机变量 ，考察 不大于某个特定的值 x这一随机事件，可以得出发生这一事件的概率 ，它是 x的确定函数，在一般情形下也依赖于采样时刻 t1。,三、概率分布、概率密度,（1）对于单个随机变量，最完整的统计描述是给出它的概率分布或概率密度。,17,学习交流PPT,由概率论公理有：,18,学习交流PPT,随机过程 在时刻 的概率分布函数,19,学习交流PPT,概率密度函数,显然有：,取值于区间 的概率为

7、：,20,学习交流PPT,此外，概率密度函数有下列性质：,，极小值,，极大值,对于平稳过程来说，其概率分布函数和概率密度函数也不依赖于采样时刻。,21,学习交流PPT,（2）多个随机变量的联合概率分布,设：采样时刻t1与t2的两个随机变量,定义：随机过程 于时刻 t1与t2的二维联合概率分布函数： 不大于x1 ，同时 不大于x2的联合概率,22,学习交流PPT,性质：,23,学习交流PPT,定义：二维概率密度函数,类似地，可以定义随机过程X(t)的n维概率分布与 n维概率密度函数。,对平稳过程，其二维概率密度只是时差 的函数,24,学习交流PPT,四、矩,从随机变量的概率分布出发，可以确定其它

8、一系列 统计特性。,定义：随机变量X(t1)的n次矩,25,学习交流PPT,一次矩：,二次矩：,26,学习交流PPT,均值x可视为信号的静态部分 x(t)-x则视为信号围绕其均值波动的动态成分 此动态成分的均方值即为方差。,二次中心矩：,27,学习交流PPT,对于各态历经过程，可以直接从时间平均求得各次矩。,28,学习交流PPT,除某些特殊情况外，确定随机变量的概率分布函 数和概率密度函数都比较困难。 在随机振动中经常遇见的正态分布过程和某些 各态历经过程，却可以用一定的程序来计算。,五、确定随机变量的概率分布函数和概率密度 函数,29,学习交流PPT,各态历经过程的集合概率与任何样本的时间概

9、率相同。,设各态历经过程的一个样本函数如图所示，T表示 样本总长， 幅值小于 x所对应的各个时间区间为 。,概率密度函数为：,单个样本的概率分布函数按如下公式计算：,样本的概率分布函数：,30,学习交流PPT,确定各态历经过程分布函数和密度函数的步骤,在样本曲线上，划一根平行于x轴的水平线， 其幅值为x。,用几何关系求出x(t)的幅值在此水平线下的时 间区段 。,31,学习交流PPT,例5-2：各态历经过程的一个样本是图示三角波，求其： 概率分布函数、概率密度函数及均值、均方值和方差。,解：,32,学习交流PPT,水平线x,几何关系：,每个周期T：,取n个周期，,则：,33,学习交流PPT,综

10、上,34,学习交流PPT,计算均值、均方值和方差：,或：,35,学习交流PPT,作业：,2. 设 为常数，X为随机变量，求 的平均值与方差。,1. 用时间平均法求例5-2的均值、均方值和方差,36,学习交流PPT,5-4 相关函数,引言： 均值和方差只是描述了随机过程单一时刻 （随机变量）的数学特征，要描述两个不同时刻的 随机变量之间的联系则要引入相关函数。,一、自相关函数 二、互相关函数,37,学习交流PPT,设随机过程 在两个任意时刻 的随机变 量为 是这两个随机变量 的二维概率密度函数。,的自相关函数，是随机变量,乘积的集合平均。 描述随机过程在两个不同时刻之间的线性依赖关系。,一、自相

11、关函数,定义：,38,学习交流PPT,对于平稳过程,39,学习交流PPT,随机过程在时刻t和t+形成的状态（随机变量）X(t)和X(t+)之间的相关程度。,物理意义：,规范化自协方差（自相关系数）：,40,学习交流PPT,自协方差,41,学习交流PPT,自相关函数的性质：, 自相关函数是偶函数, 周期平稳过程的自相关函数也是周期函数， 其周期与过程的周期相同。,42,学习交流PPT,自相关函数的性质：,=0时的自相关函数就是均方值,43,学习交流PPT, 如果随机过程不是周期过程，则：,44,学习交流PPT, 自相关函数是一个有界函数,一般越大，则两时刻的随机变量X(t1)和X(t1+)之 间

12、的相关性愈差。,，Rx()。,45,学习交流PPT,例5-3 ：初相位是随机的正弦函数随机过程,x=Asin(t+)，其中是随机变量，其均值为零。, 原随机过程是周期各态历经的，其周期为,。而Rx()亦是周期函数,46,学习交流PPT,例5-4：交流电电路中的电压和电流分别为 v=Vsint，i=Isin(t+) 瞬时功率P=vi=VIsintsin(t+) 相应的平均功率为：,而自相关函数,比较P和Rx()，可见Rx()体现了随机过程的平均功率随时差的变化。,47,学习交流PPT,二、互相关函数,两个随机过程：,48,学习交流PPT,定义：,的互相关函数。,描述：两个随机过程之间的线性依赖关

13、系。,一般：,49,学习交流PPT,对于平稳过程：,对于各态历经过程：,其中x(t)、y(t+)分别是随机过程 的 代表性样本函数。,50,学习交流PPT,（1）,互相关函数一般不是的偶函数，一般也不在 0时取极大值。,Rxy()与Ryx() = Rxy(-)之间一般并没有关系。,互相关函数的性质：,51,学习交流PPT,（2） Rxy()是有界函数,证明：,52,学习交流PPT,互协方差：,结论：,53,学习交流PPT,互相关系数：,当X，Y的均值都等于零时，互相关系数就等于 规范化互协方差。,定义：,规范化互协方差,54,学习交流PPT,（3）Rxy()与Rx(0)，Ry(0)之间的关系,

14、55,学习交流PPT,（4）两个统计独立的随机变量一定是不相关的， 但两个不相关的随机变量不一定是统计独立的。,证明：,56,学习交流PPT,5-5 功率谱密度函数,一、自功率谱密度函数 二、互功率谱密度函数,自相关函数Rx()描述“平均功率”随时差的变化“平均功率”的时间结构。 功率谱密度S x(f)：描述“平均功率”在频域（谱域）的分布频率结构。 二者在不同的域（时域或频域）反映着同一个统计特性。在不同的场合，各有所长，相辅相成。,57,学习交流PPT,一、自功率谱密度函数,定义：,圆频率（角频率）,存在上述傅立叶变换的条件：,一般地，Rx() , Rx()的傅立叶变换一般是存在的。,58

15、,学习交流PPT,为什么称为“功率谱” ？,设 是作用在R1上的电压信号，则 是瞬时功率信号，而平均功率,一方面，此式表示平均功率 的时间结构，即各个瞬时的功 率 对于平均功率的贡献。 另一方面，又表示了平均功率的频率结构，即各种频率的功 率成分Sx(f)df对于平均功率的贡献，因此称为功率谱。,59,学习交流PPT,为什么称为功率谱“密度”,量纲：,60,学习交流PPT,自谱密度Sx(f)的性质：,（1）Sx(f)0,（2）,61,学习交流PPT,相应地，,62,学习交流PPT,（3）随机过程的自谱在整个频域上的积分等于 随机过程的均方值。,（4） 双边谱,工程上，把自谱定义在正半轴上，称为

16、单边谱。,63,学习交流PPT,（5） 导数过程的自谱,64,学习交流PPT,从Parseval 定理角度来定义功率谱密度 信号在时域的总能量等与它在频域的总能量,65,学习交流PPT,设 是平稳随机过程的一个样本函数，一般情况下它不一定能满足绝对可积的条件，为此引入辅助函数：,根据Parseval 定理,66,学习交流PPT,对于各态历经过程：,67,学习交流PPT,二、互功率谱密度函数,定义：,对于各态历经过程：,对于平稳过程：,68,学习交流PPT,互功率谱密度函数性质（自学）,（1）互谱一般是复函数,（2）,（3）互谱是有界函数,（4）两个不相关且均值为零的随机过程,69,学习交流PP

17、T,作业：,5.1 5.3 5.5 5.9,70,学习交流PPT,5-6 线性系统在随机激励下的响应,设系统受到的平稳随机激励为 ，它的一个 样本函数 ，引起的响应可根据Duhamel积分 得到：,讨论线性振动系统只受一个平稳随机激励时的情况。,所有的样本函数引起的响应的全体 为一个随机过程，讨论： 统计性质以及它与 激励 的统计性质的关系。,71,学习交流PPT,1 响应的均值,对于平稳过程：,因此,72,学习交流PPT,2. 响应的自相关函数,根据Duhamel积分，有,73,学习交流PPT,因此,线性振动系统受到平稳过程激励后的响应也是平稳过程。 如果激励是各态遍历过程，则响应也是各态遍

18、历过程。,74,学习交流PPT,3. 响应的自谱以及均方值,在上面的证明过程中用到了傅里叶变换的性质：,75,学习交流PPT,响应的均方值,76,学习交流PPT,4. 激励与响应的互相关函数和互谱,激励与响应的互相关函数为,77,学习交流PPT,例5-5 如果激励为白噪声，求激励与响应的互相关 函数。,激励的自相关函数为,78,学习交流PPT,激励为白噪声时激励与响应的互相关函数,在白噪声激励下激励与响应的互相关函数与系统的脉冲响应成正比。 在实际中可以用这个性质来测定系统的脉冲响应。,79,学习交流PPT,激励与响应的互谱,80,学习交流PPT,例5-6：求单自由度系统的基础以白噪声运动时 响应的自谱和均方值,解：基础运动的自谱为,基础运动时系统的频响函数为,因此响应的自谱为,81,学习交流PPT,响应的均方值为：,82,学习交流PPT,大作业1： 汽车的