随机振动课件(全)

5-1 引言w前面三章所考察的振动都是确定性振动振动系统的规律可以用时间的确定性函数来描述w振动系统的物理量可以用随时间变化的确定函数来描述，因此，确定性振动中的物理量在将来某一时刻的值是可以预测的。w举例：单自由度系统的简谐强迫振动 1自然界和工程上还存在着另一类振动，它们的规律不能用时间的确定函数来描述，但又具有一定的统计规律性。在数学上这类振动可以用随机过程来加以描述。这类振动称为随机振动非确定非确定非确定非确定而又具有统计规统计规统计规统计规律律律律。举例：车辆行驶时由于路面不平引起的振动2确定性系统+确定性激励 确定性响应 确定性系统+随机激励 随机响应 随机系统+任何激励 随机响应3随机振动与确定性振动的本质区别在于它一般指的不是单个现象，而是一个包含着大量现象的集合；从集合中的单个现象来看似乎是杂乱的，但从总体来看却存在着一定的统计规律性。因此，它虽然不能用时间的确定函数来描述，但能用统计特性来描述。在确定性振动中，系统的激励与响应之间有着确定的函数关系，而在随机振动中，只能满足于确定它们的统计特性之间的关系。4 介绍随机振动中物理量的描述方法（相关函数、功率谱密度）。讨论受随机激励的振动系统的激励、系统特性、响应三者统计规律性之间的关系。本章内容：5过程：物理量随时间变化的情况随机过程：无法准确预知物理量随时间的变化情 况，但其变化规律服从统计规律随机过程是大量现象的一个数学抽象，理论上是由无限多个无限长的样本组成的集合。：样本函数对于随机现象，我们感兴趣的往往不是各个样本本身，而是力图从这些样本得出总体的统计特性。5-2 随机过程6 随机过程 的所有样本函数 在时刻 的值 构成一个随机变量 对随机变量求集合平均称为随机过程 在时刻 的集合平均值。一般情况下 依赖于采样时刻 。一、集合平均.平稳过程5-3 随机过程的数字特征7 、构成两个随机变量对它们的乘积 求集合平均称为随机过程 于时刻 与 的自相关函数。它是时差 的函数，在一般情况下，它也依赖于采样时刻 ，反映这两个时刻的随机变量 的统计联系。自相关函数8随机过程可以根据其统计特性是否随采样时刻而变化来进行分类。统计特性依赖于采样时刻的过程非平稳过程统计特性不依赖于采样时刻的过程平稳过程平稳过程9平稳过程的特点u集合平均值为常数 u相关函数仅仅依赖于时差10二、时间平均.各态历经过程随机过程的每一个样本函数可以在时域内求得:时间平均值 对第k个样本求平均 补充：求函数的平均值 11 自相关函数 在一般情况下，对于不同的样本将得到不同的 和12各态历经过程定义：各个样本的时域统计值都是等同的，而且任一个样本函数在时域的统计值与任一时刻的随机变量的统计值相等。13结论：v集合平均与时刻t1无关，而时间平均与样本标号k无关，是过程为各态历经的充要条件。v从任何一个样本得出的时间特性就等于集合平均特性，将使数据处理容易的多。举例：半主动悬挂特性评价 v各态历经过程一定是平稳的，反之不一定。平稳+各个样本的统计特性相同 各态历经 14例5-1：求正弦函数的相关函数 物理意义：表示样本函数 与其延时 时刻得到的 之间波形的相似程度。，相似程度最高。，相似程度最低。15各个样本函数在时刻 t1的值 构成一个随机变量 ，考察 不大于某个特定的值 x这一随机事件，可以得出发生这一事件的概率 ，它是 x的确定函数，在一般情形下也依赖于采样时刻 t1。三、概率分布、概率密度（1）对于单个随机变量，最完整的统计描述是给出它的概率分布或概率密度。16由概率论公理有：17 随机过程 在时刻 的概率分布函数 18概率密度函数显然有：取值于区间 的概率为：19此外，概率密度函数有下列性质：，极小值，极大值 对于平稳过程来说，其概率分布函数和概率密度函数也不依赖于采样时刻。20（2）多个随机变量的联合概率分布设：采样时刻t1与t2的两个随机变量 定义：随机过程 于时刻 t1与t2的二维联合概率分布函数：不大于x1，同时 不大于x2的联合概率 21性质：22定义：二维概率密度函数 类似地，可以定义随机过程X(t)的n维概率分布与 n维概率密度函数。对平稳过程，其二维概率密度只是时差 的函数 23四、矩 从随机变量的概率分布出发，可以确定其它一系列统计特性。定义：随机变量X(t1)的n次矩 24一次矩：二次矩：25v均值x可视为信号的静态部分vx(t)-x则视为信号围绕其均值波动的动态成分v此动态成分的均方值即为方差。二次中心矩：26对于各态历经过程，可以直接从时间平均求得各次矩。27除某些特殊情况外，确定随机变量的概率分布函数和概率密度函数都比较困难。在随机振动中经常遇见的正态分布过程和某些各态历经过程，却可以用一定的程序来计算。五、确定随机变量的概率分布函数和概率密度 函数28各态历经过程的集合概率与任何样本的时间概率相同。设各态历经过程的一个样本函数如图所示，T表示样本总长，幅值小于 x所对应的各个时间区间为 。概率密度函数为：单个样本的概率分布函数按如下公式计算：样本的概率分布函数：29确定各态历经过程分布函数和密度函数的步骤在样本曲线上，划一根平行于x轴的水平线，其幅值为x。用几何关系求出x(t)的幅值在此水平线下的时间区段 。30例5-2：各态历经过程的一个样本是图示三角波，求其：概率分布函数、概率密度函数及均值、均方值和方差。解：31 水平线x 几何关系：每个周期T：取n个周期，则：32综上33计算均值、均方值和方差：或：34作业：2.设 为常数，X为随机变量，求的平均值与方差。1.用时间平均法求例5-2的均值、均方值和方差355-4 相关函数引言：均值和方差只是描述了随机过程单一时刻（随机变量）的数学特征，要描述两个不同时刻的随机变量之间的联系则要引入相关函数。一、自相关函数二、互相关函数36设随机过程 在两个任意时刻 的随机变量为 是这两个随机变量的二维概率密度函数。的自相关函数，是随机变量 乘积的集合平均。描述随机过程在两个不同时刻之间的线性依赖关系。一、自相关函数定义：37对于平稳过程 38随机过程在时刻t和t+形成的状态（随机变量）X(t)和X(t+)之间的相关程度。物理意义：规范化自协方差（自相关系数）：39自协方差 40自相关函数的性质：自相关函数是偶函数 周期平稳过程的自相关函数也是周期函数，其周期与过程的周期相同。41自相关函数的性质：=0时的自相关函数就是均方值42 如果随机过程不是周期过程，则：43 自相关函数是一个有界函数 一般越大，则两时刻的随机变量X(t1)和X(t1+)之间的相关性愈差。，Rx()。44例5-3：初相位是随机的正弦函数随机过程x=Asin(t+)，其中是随机变量，其均值为零。原随机过程是周期各态历经的，其周期为。而Rx()亦是周期函数,4546例5-4：交流电电路中的电压和电流分别为v=Vsint，i=Isin(t+)瞬时功率P=vi=VIsintsin(t+)相应的平均功率为：而自相关函数 比较P和Rx()，可见Rx()体现了随机过程的平均功率随时差的变化。47二、互相关函数两个随机过程：48定义：的互相关函数。描述：两个随机过程之间的线性依赖关系。一般：49对于平稳过程：对于各态历经过程：其中x(t)、y(t+)分别是随机过程 的 代表性样本函数。50（1）互相关函数一般不是的偶函数，一般也不在 0时取极大值。Rxy()与Ryx()=Rxy(-)之间一般并没有关系。互相关函数的性质：51（2）Rxy()是有界函数 证明：52互协方差：结论：53互相关系数：当X，Y的均值都等于零时，互相关系数就等于规范化互协方差。定义：规范化互协方差54（3）Rxy()与Rx(0)，Ry(0)之间的关系55（4）两个统计独立的随机变量一定是不相关的，但两个不相关的随机变量不一定是统计独立的。证明：565-5 功率谱密度函数一、自功率谱密度函数二、互功率谱密度函数v自相关函数Rx()描述“平均功率”随时差的变化“平均功率”的时间结构。v功率谱密度S x(f)：描述“平均功率”在频域（谱域）的分布频率结构。v二者在不同的域（时域或频域）反映着同一个统计特性。在不同的场合，各有所长，相辅相成。57一、自功率谱密度函数定义：圆频率（角频率）存在上述傅立叶变换的条件：一般地，Rx()Rx()的傅立叶变换一般是存在的。58为什么称为“功率谱功率谱”？设 是作用在R1上的电压信号，则 是瞬时功率信号，而平均功率一方面，此式表示平均功率 的时间结构，即各个瞬时的功 率 对于平均功率的贡献。另一方面，又表示了平均功率的频率结构，即各种频率的功 率成分Sx(f)df对于平均功率的贡献，因此称为功率谱。59为什么称为功率谱“密度密度”量纲：60自谱密度Sx(f)的性质：（1）Sx(f)0（2）61相应地，62（3）随机过程的自谱在整个频域上的积分等于 随机过程的均方值。（4）双边谱工程上，把自谱定义在正半轴上，称为单边谱。63（5）导数过程的自谱64从Parseval 定理角度来定义功率谱密度信号在时域的总能量等与它在频域的总能量65 设 是平稳随机过程的一个样本函数，一般情况下它不一定能满足绝对可积的条件，为此引入辅助函数：根据Parseval 定理6667 对于各态历经过程：68二、互功率谱密度函数定义：对于各态历经过程：对于平稳过程：69互功率谱密度函数性质（自学）（1）互谱一般是复函数（2）（3）互谱是有界函数（4）两个不相关且均值为零的随机过程70作业：5.15.35.55.9715-6 线性系统在随机激励下的响应v 设系统受到的平稳随机激励为 ，它的一个 样本函数 ，引起的响应可根据Duhamel积分得到：v讨论线性振动系统只受一个平稳随机激励时的情况。v所有的样本函数引起的响应的全体 为一个随机过程，讨论：统计性质以及它与 激励 的统计性质的关系。721 响应的均值 对于平稳过程：因此 732.响应的自相关函数 根据Duhamel积分，有 74因此 线性振动系统受到平稳过程激励后的响应也是平稳过程。如果激励是各态遍历过程，则响应也是各态遍历过程。753.响应的自谱以及均方值 在上面的证明过程中用到了傅里叶变换的性质：76响应的均方值 774.激励与响应的互相关函数和互谱激励与响应的互相关函数为 78例5-5 如果激励为白噪声，求激励与响应的互相关 函数。激励的自相关函数为 79激励为白噪声时激励与响应的互相关函数 在白噪声激励下激励与响应的互相关函数与系统的脉冲响应成正比。在实际中可以用这个性质来测定系统的脉冲响应。80激励与响应的互谱81例5-6：求单自由度系统的基础以白噪声运动时 响应的自谱和均方值 解：基础运动的自谱为 基础运动时系统的频响函数为 因此响应的自谱为 82响应的均方值为：83大作业1：汽车的二自由度模型 为悬挂质量（车身质量），为非悬挂质量（车轮质量），为悬挂刚度，悬挂阻尼系数，为车轮刚度84推导：1 运动微分方程 2 固有频率和振型 3 车身位移 对 的频响函数 4 车身加速度 对 的频响函数 5 悬挂动位移 对 的频响函数 6 悬挂动速度 对 的频响函数85大作业2：汽车的单自由度模型路面功率谱密度其中，路面不平度系数 参考空间频率，B级路面：F级路面：861 证明：路面速度谱2 求车身位移响应的自谱和均方值3 求车身加速度响应的自谱和均方值87xiexie!xiexie!谢谢！谢谢！xiexie!xiexie!谢谢！谢谢！