matlab 功率谱密度 汉宁窗_如何理解随机振动的功率谱密度?

7489a6a1459486817c0e891da10d112c.png
一、随机信号和正太分布有什么关系?
二、时域、频域之间功率守恒?
三、自相关又是个什么玩意?

作为一个工程师,很多人对随机振动看着熟悉,却又实际陌生。熟悉是因为几乎每个产品在出厂时都要求要做随机振动试验,陌生是因为当面对用户所给的功率密度谱有时会感到很茫然,尤其是这个功率谱的单位居然是

,简直是反人类,为啥整这么麻烦,不能给个加速度
直接振吗?哈哈,这个还真不能,下面咱们就来抽丝剥茧,看看随机振动功率密度谱到底是咋个回事?

随机振动,是因为振动源是随机信号,顾明思议,就是信号的发生带有随机性,无法用一个明确的函数把它表达出来,一个典型的随机信号大概长这个样子:

24490d0d14099a474aa192aff124c4aa.png

可以看出,乍一看完全没有规律可言,高度不规则,无规律的,不可预估也不可重复,每次测量都不一样。那随机信号是不是就是不可描述的呢?首先我们研究一下什么是随机,随机对应不确定性,不确定性在物理学和数学上是一个不受欢迎的词,老板问你图纸什么时候画完,你敢回答“随机吧”,或者更佛系一点“随缘吧”?相信等待你的可不仅仅是白眼哦,这也是为啥当初波恩将概率解释引入到量子力学时受到爱因斯坦强烈的反对,“上帝是不掷骰子”的典故便来源于此。然后诸多证据表明:也许随机是这个世界的本质特性之一,这可以由一个著名的数学定理来证明:中心极限定理。


一、随机信号和正太分布有什么关系

中心极限定理有一组,但基本可以用一句通俗的话来概括它们:大量相互独立的随机变量,其平均值正态分布。对中心极限定理最形象的解释是高尔顿钉板实验:小球在下落过程中碰到很多个钉子,每次碰撞都是一个二项式的随机过程:以同等的概率通过钉子左侧或者右侧,小球最后到达的位置,是这很多个“左右”随机变量相加后的平均位置。不难看出,这个平均值落在中心处的概率最大,小球聚集最多,但也可能向左或向右偏,偏离越大,小球的数目越少,不同位置的不同小球数便形成了一个“分布”,中心极限定理则是从数学上证明了,这个分布的极限是正态分布。

25a786478c3cb6fdba4477979a6154ac.png
图片来自网络

高尔顿钉板实验实质上是一组二项式分布,从数学上还可以证明:中心极限定理的条件可以从二项分布推广到独立同分布随机序列,以及不同分布的随机序列。也就是说:在一定条件下,各种随意形状概率分布生成的随机变量,它们加在一起的总效应,是符合正态分布的。这也是为啥正太分布这么常见的原因,因为实际上的随机生物过程或物理过程,都不是只由一个单独的原因产生的,它们受到各种各样随机因素的影响。比如产品加工免不了有误差,而误差形成的原因五花八门,各种各样。就算我们打开上帝视角,够分别清楚产生误差的每种单一原因,误差的分布曲线可能不是高斯的,但是,所有误差加累计在一起时,通常得到一个正态分布。下图是30组随机变量,和值随数据增加时分布情况。

e7a457863a6bc816938fe270144f7a94.gif
30组随机变量(每组100万数据),和值的随数据增加时分布情况

总之,中心极限定理告诉我们:无论引起过程的各种因素的基本分布是什么样的,当实验次数 充分大时,所有这些随机分量之和近似是一个正态分布的随机变量。也就是说:对于平稳随机过程而言,其分布是趋于正太的。

10bd2b43b68672d46f549be515c62995.png
图片来自:http://blog.sciencenet.cn

我们知道,正太分布的概率密度的表达式为:

4efb2bae56b7583e448fe2b05b37ab73.png

其中

为数学期望,代表着信号的平均值,

为方差,表示这信号偏离平均值的程度,

为均方值,在工程上可以一般看成信号的平均功率,其平方根为有效值(RMS值),

很容易得到:

所以当平均值为零的时候,均方值与方差相等。从图形上看,平均值

决定了正太分布的位置,均方根值,也就是标准差
,决定了分布的幅度。如果正态分布概率密度一旦确定,其数学期望
和方差
也就确定了。对于随机信号而言,数学期望
基本为零(或者去除直流分量后为零),所以唯一确定的量就剩均方值,此时均方值和方差一致,并且

如果仔细观察上式,可以发现,从量纲来看,

代表了能量概念,而方差
代表了平均功率概念。

前面我们说到,平稳随机信号是趋于正太分布的,决定正太分布的两个参数:平均值

基本为零(不为零时,直流分量也很容易处理),唯一能表征随机信号的就剩均方值(
时等于方差)
,也就是平均功率了。再具体一点,对于平稳随机信号而言,描述
是没有意义的,描述平均功率才有意义。如何才能描述平均功率呢?这时我们就要用到
频谱分析的概念了,可能有些人会觉得这个概念很生疏,我要说我们每天一睁眼就在进行谱分析你会信吗?事实还真是这样。我们的眼睛就是一双最精密的谱分析仪:过滤到可见光以外的所有光线,并且可见光按照频谱(频率)进行精确分类。我们知道光本质是一种电磁波,既然是波,就有波长
和波动频率
,光速
一个定值,波长和波动频率成反比。人眼可识别出波长在400纳米至700纳米的光线,并且根据波长不同(频率不同)将可见光分成红、橙、黄、绿、蓝、靛、紫色。当我们说“一个东西的表面是红色”的时候,这是正常人的定性描述,换成喜欢定量描述(装X)工程师会怎么说呢?“这个东西表面特性会选择性反射一种电磁波,这种电磁波波长约700纳米,波动频率 约
赫兹”。

152f774d04ebd5b2ec111eea49dd5472.png
图片来自:https://www.picquery.com

那频谱分析是什么呢?它是一种将复杂号分解为较简单信号的技术,许多物理信号均可以表示为许多不同频率简单信号的和,找出一个信号在不同频率下的信息(可能是幅度、功率、强度或相位等)的作法就是频谱分析。比如通过一个三棱镜就可进行色散试验,将白光分解成频率各不相同的单色光,我们经常在雨后见到的彩虹也是一种色散现象。

85c420a9935e8fb603b526eec49075bd.gif

图片来源:http://www.bolvan.ph.utexas.edu

在傅里叶变换背后到底有什么小秘密 这篇文章中,我们介绍了傅里叶变换概念,它是频谱分析的基础,这篇文章中介绍了复指数函数是如何构成一个完备正交基的,以及如何进行傅里叶变换及其逆变换。但是,说这么多,频谱分析和随机信号有什么关系?


二、时域、频域之间功率守恒?

频谱分析和平稳随机信号可以通过一个叫“帕斯瓦定理”的公式联系起来,一个经常被低估的定理,大概长这个样子:

其中

的傅里叶变换,这个公式看着还挺对称,被积函数都是自变量的平方,貌似蕴含着什么深刻的意义。翻译过来,这个定理说的是什么事呢?那就是
信号的能量(或者平均功率)无论在时域看,还是在频域看,都是一样的。人类习惯于在时域观察一个信号,那会不会存在另一个物种或者空间,他们喜欢在频域观察呢?但无论如何,帕斯瓦定理告诉我们,不论频域还是时域,信号所代表的能量是守恒的,而守恒的本质是对称性,那是不是预示着时域和频域有着某种高度的对称呢?

好了,扯远了,哲学上的事情不归工程师管,我们要所的事情就是如何定量化的的去描述不确定的东西,比如随机信号。前面我们说了,平稳随机信号

的幅值是呈现正太分布的,其平均值接近为零(或去除直流分量后),均方值或方差(也就是平均功率)是固定的,既然平均功率(或能量)在时域和频域是守恒的,而在时域
是随机的,不可描述,那我们可以换到频域去啊:

其中:

表示了信号的平均功率(或能量)在频域上的分布,即单位频带的功率随频率变化的情况,故称之为信号的自功率谱密度函数,简称自功率谱或自谱。
轴包围的面积等于信号
的平均功率,即
的幅值分布的方差或均方值。同时可以看出,当
表示加速度时(单位为
),
的单位就变成了

三、自相关又是个什么玩意?

现在,事情就貌似变得简单了,我们将时域信号

进行傅里叶变换,得到频域分布
,然后平方积分就可以得到功率谱密度了啊。理论上,我们确实可以这么做,但是实际上,这样做有一个小困难,那就是经典傅里叶变换不是一直存在的,前提的:信号
绝对值可积,即
,而实际上平稳随机变量是不满足的。为了解决随机信号的傅里叶变换问题,我还得引入一个概念:自相关,先看定义:随机过程
的自相关函数定义为在时刻
和时刻
的随机变量乘积的平均值,
是时移,当平均时间
时,平均值的极限便是自相关函数,其数学表达式为:


定义有点拗口,可以先不理解,它说明一个什么事呢?它反映了随机信号本身在不同时刻的相互关系,再直白一点:把一个信号平移一段距离,跟原来有多相似。我们都知道,重要的证件上的照片,一般都要求是近半年的,道理显而易见,时间久了跟现在长得可能就不一样了。下图显示了C罗近十年容颜的变化,都说女大十八变(即使不整容),男的又何尝不是。

81aadbd2ec3f3e936c5496eb457c23bb.gif
图片来源于网络

很显然的一点是:随着时间的推移,C罗的长相是发生明显变化的,时间越短,长相越接近,可以想象,当时间足够长后,不知道我们还能不能认出他来,不信,你看:

4b5b53ef1766c079edc4f74d73f8dc97.png
马爸爸不要告我

对于随机信号,特别白噪声随机信号来说,当时移

非常小时,
相差很小的概率很大,这时
值非常大,表示关系密切。特别当
时,
值最大,等于方均值,也就是平均功率,表示完全相关。当时移
较大时,
相差很小的概率很小。作平均计算正负对消,
值很小。并且随着
值的增大,
值很快衰减到零,表示
之间没有依赖关系,对一般的宽带随机振动,时间间隔很远的二个随机量之间不存在任何固定关系。简单点来说:
对于平稳随机信号,自相关函数将信号的平均功率向
这一点集中,
时自相关函数快速衰减为零。

对于周期信号,自相关函数可以把随机信号中的周期成份检测出来,这是因为任何周期信号在所有的时移上都有一定形状的自相关函数图形。对于周期信号来说,经过一个周期后又精确的重复过去的时间历程,因此当时移超过该周期时,其自相关函数必然重复前一段的形状。所以若在自相关函数图上发现时移趋于无穷大,

,而有某种周期性,则说明该随机振动信号混有周期信号成分。
简单点来说:自相关函数能够检测出信号内部蕴藏的周期组分,而过滤掉了周期组分的相位信息。

对于上面两段话不是很理解的,可以看下面两幅图,分别是白噪声随机信号和随机相位正弦信号及其自相关信号。

d24e5be58572822b3e599e64f79403ac.png

caee9155debf78a83b88e14d0efd636a.png

说了这么多,和我们要进行功率谱分析有什么关系?前面我们分析了,信号的功率在时域和频域都是满足守恒定律的

而功率被定义成幅值的平方的时间平均分量,而这个过程,也可以看成是去除频域谐波分量的相位信息的过程,因为本质来说,一个简谐信号的相位是不影响其功率的。而自相关函数,也具有去除信号相位的功能,那自相关函数和功率密度谱是不是有什么深刻的联系呢?答案是肯定的,那就是维纳-辛钦定理(Wiener–Khinchin theorem),这个定理表明:信号的自相关函数与功率密度谱是一对傅里叶变换对:

也就是说:一个信号的功率密度谱,就是其自相关函数的傅里叶变换

总结一下整个逻辑:对于一个随机信号而言,时域信息是杂乱无章的,唯一的确定性信息但是在统计意义下得到的,即幅值呈正太分布,均方值也就是平均功率是固定的。根据帕斯瓦定理,信号的平均功率在时域和频域是守恒的,按道理说直接对时域信号进行傅里叶变换再取平方就可以。但不幸的事,随机信号的不满足傅里叶变换绝对值可积的条件,严格意义傅里叶变换不存在,于是发明了自相关函数的概念,将信号的蕴含的周期信号识别出来,并将相位信息去掉(相位不影响平均功率),于是就出现了我们在教材上见到的最终形式维纳-辛钦定理:一个信号的功率密度谱,就是其自相关函数的傅里叶变换。提炼一下就是:随机信号→幅值正太分布→均方值(平均功率)→帕斯瓦定理(功率守恒)→自相关函数(去除相位信息)→维纳-辛钦定理(最终形式)。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/333768.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

深度学习pytorch--softmax回归(二)

softmax回归的从零开始实现实验前思考获取和读取数据获取数据集查看数据集查看下载后的.pt文件查看mnist_train和mnist_test读取数据集查看数据迭代器内容初始化模型参数定义softmax函数定义模型定义损失函数计算分类准确率模型评价--准确率开始训练可视化总结完整代码实验前思…

HTML块级元素/块标签/块元素

文章目录块元素的特点块元素清单block level element. 块级元素想在同一行显示需浮动或者 display:inline。 块元素的特点 每个块级元素都是独自占一行,其后的元素也只能另起一行,并不能两个元素共用一行; 元素的高度、宽度、行高、顶边距、…

物联卡查询流量_电信物联卡官网是多少?如何快速查询流量信息?

高速率设备的使用场景需要用到5G,中速率LET-Cat1应用范围更广,而低速率则要靠窄带物联网NB-IOT去维护了。这三种网络制式全都与物联网息息相关,这就能知道为什么国家层面对物联网基础设施建设这么重视了。电信物联卡在智能化硬件中有优秀表现…

java8日期转时间戳_Java 8日期和时间

java8日期转时间戳如今,一些应用程序仍在使用java.util.Date和java.util.Calendar API,包括使我们的生活更轻松地使用这些类型的库,例如JodaTime。 但是,Java 8引入了新的API来处理日期和时间,这使我们可以对日期和时间…

HTML行内元素/行级元素/内联元素/行标签/内联标签/行内标签/行元素

文章目录行内元素的特点行内元素清单可变元素列表inline element. 也叫行级元素、内联元素。行内元素默认设置宽度是不起作用,需设置 display:inline-block 或者 block 才行。 行内元素的特点 可以和其他元素处于一行,不用必须另起一行; 元…

深度学习pytorch--softmax回归(三)

softmax回归的简洁实现获取和读取数据定义和初始化模型softmax和交叉熵损失函数定义优化算法模型评价训练模型小结完整代码前两篇链接: 深度学习pytorch–softmax回归(一) 深度学习pytorch–softmax回归(二) 本文使用框架来实现模型。 获取和读取数据 我们仍然使用Fashion-M…

正则表达式的分类

文章目录一、正则表达式引擎二、正则表达式分类三、正则表达式比较四、Linux/OS X 下常用命令与正则表达式的关系一、正则表达式引擎 正则引擎大体上可分为不同的两类:DFA 和 NFA,而 NFA 又基本上可以分为传统型 NFA 和 POSIX NFA。 DFA(Deterministic …

spock测试_使用Spock测试您的代码

spock测试Spock是针对Java和Groovy应用程序的测试和规范框架。 Spock是: 极富表现力 简化测试的“给定/何时/然后” 语法 与大多数IDE和CI服务器兼容。 听起来不错? 通过快速访问Spock Web控制台,您可以非常快速地开始使用Spock。 当您有…

深度学习pytorch--多层感知机(一)

多层感知机隐藏层激活函数ReLU函数sigmoid函数tanh函数多层感知机小结我们已经介绍了包括线性回归和softmax回归在内的单层神经网络。然而深度学习主要关注多层模型。在本节中,我们将以多层感知机(multilayer perceptron,MLP)为例…

太阳能板如何串联_光伏板清洁专用的清洁毛刷

光伏发电是利用半导体界面的光生伏特效应将光能直接转变为电能的一种技术。主要由太阳电池板(组件)、控制器和逆变器三大部分组成。主要部件由电子元器件构成。太阳能电池经过串联后进行封装保护可形成大面积的太阳电池组件,再配合上功率控制…

java 异步等待_Java中的异步等待

java 异步等待编写异步代码很困难。 试图了解异步代码应该做什么的难度更大。 承诺是尝试描述延迟执行流程的一种常见方式:首先做一件事,然后再做另一件事,以防万一出错时再做其他事情。 在许多语言中,承诺已成为协调异步行为的实…

cass生成曲线要素_干货在线 | CASS入门指南——道路断面计算土方

CASS操作指南——道路断面计算土方法小伙伴们赶紧学起来!道路类的土方工程,主要用CASS的断面法土方计算之道路断面来计算。整个计算过程主要分为以下四步:菜单截图第一步:绘制道路中心线道路的中心线,一般由直线段和缓…

正则表达式的捕获性分组/反向引用

文章目录分组捕获性分组和反向引用分组 正则的分组主要通过小括号来实现,括号包裹的子表达式作为一个分组,括号后可以紧跟限定词表示重复次数。如下,小括号内包裹的 abc 便是一个分组: // (abc) 表示匹配一个或多个"abc"&#xf…

深度学习pytorch--多层感知机(二)

多层感知机的从零开始实现获取和读取数据定义模型参数定义激活函数定义模型定义损失函数训练模型小结我们已经从上一节里了解了多层感知机的原理。下面,我们一起来动手实现一个多层感知机。首先导入实现所需的包或模块。 import torch import numpy as np获取和读取…

jwt同一会话_在会话中使用JWT

jwt同一会话这个话题已经在黑客新闻,reddit和博客上讨论了很多次。 共识是–请勿使用JWT(用于用户会话)。 而且我在很大程度上同意对JWT的典型论点 , 典型的“但我可以使其工作……”的解释以及JWT标准的缺陷的批评 。 。 我不会…

表必须要有主键吗_玄关隔断什么材质好?玄关隔断必须要做吗?

为了避免一到门口就能够看到全部室内的东西,为了更好的保护家居的隐私,目前有很多人都会在玄关的位置加一个隔断,而玄关隔断什么材质好?在做玄关隔断的时候,有些人觉得做了隔断会太浪费空间了,而玄关隔断必须要做吗?…

深度学习pytorch--多层感知机(三)

使用pytorch框架实现多层感知机和实现softmax回归唯一的不同在于我们多加了一个全连接层作为隐藏层。它的隐藏单元个数为256,并使用ReLU函数作为激活函数。#模型的核心代码为:nn.Linear(num_inputs, num_hiddens),nn.ReLU(),nn.Linear(num_hiddens, num_outputs),

mysql 日期索引的使用_日期使用

mysql 日期索引的使用时区糟透了。 特别是夏令时。 我不介意像与此行为相关的编程错误那样,不停移动时钟或失去一个小时的睡眠。 更糟糕的是Java的旧日期/时间API。 Java社区通过JSR 310公开承认了这一点,该JSR 310代替了Java Date&Time API…

axure元件库 文件上传_手把手教你打造一套属于产品经理自己的元件库

之前有篇文章聊完了如何打造一套属于自己的原型图设计规范,今天咱们来聊聊如何打造一套属于自己的元件库。毕竟,每个追求效率的PM,总会拥有一个专属自己的Axure元件库,并不断打磨优化。今天就教大家一步一步创建属于自己的Axure元…

KMP算法笔记

1.KMP算法本质上就是对朴素匹配算法(BF)的一个优化,减少朴素匹配算法中不必要匹配的次数,核心代码和朴素匹配算法差不多,BF是移动字串逐个字符匹配,每次模式串(子串)匹配只移动一个字符单位,而KMP算法是每次模式串匹配…