随机振动信号的特征——PSD(功率谱密度)

这篇具有很好参考价值的文章主要介绍了随机振动信号的特征——PSD(功率谱密度)。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

理论知识

1、背景
信号可分为确定性信号和随机信号。确定性信号是每个时间点上的值可以用某个数学表达式或图标唯一地确定的信号;而随机信号(random signal),幅度未可预知但又服从一定统计特性的信号,又称不确定信号(百度百科的解释)。随机信号是普遍存在的,也不能用一个确切的数学公式来描述,因为也不能准确进行预测的信号。正是因为随机信号是随机的,所以只能用统计的方法进行描述,在一定的准确性或可信性范围内用统计学规律去表征随机信号的特性。
如果随机信号的概率特性不随时间变化而变化,则成为平稳随机信号。

2、统计特征量
随机振动功率谱密度,嵌入式知识储备,概率论,傅立叶分析

以上这些统计学特征用于分析信号的特征,并计算出PSD,功能如下:
数学期望值,描述随机信号的平均值。
方差值,描述随机信号幅度变化的强度。
概率密度函数,是描述信号振幅数值的概率。
相关函数,描述随机信号的每两个具有一定时间间隔的幅度值之间的联系程度的数值,它是时间间隔的一个函数。
功率谱密度,描述随机信号在平均意义上的功率谱特性。

3、功率谱密度
由中心极限定理:大量相互独立的随机变量,其平均值正态分布。可知平稳随机信号是趋于正态分布的,所以可以用统计学上特征量来描述随机信号。平稳随机信号的数学期望基本为零,数学期望为0时,方差等于均方值,也就是平均功率。为了描述平均功率,需要使用频谱分析。对于一个随机信号而言,时域信息是杂乱无章的,唯一的确定性信息但是在统计意义下得到的,即幅值呈正态分布,均方值也就是平均功率是固定的。

根据帕塞瓦尔定理,信号在时域的总功率等于在频域的总功率,所以只需要对时域信号进行傅里叶变换即可转换到频域,得到频域分布,只需要将频域分布平方积分就可以得到功率谱密度(PSD)。
但是因为时域的信号是随机的,无法用数学表达式描述,随机信号的不满足傅里叶变换绝对值可积的条件,严格意义傅里叶变换不存在,也就无法通过傅里叶变换将时域转成频域。(傅里叶理论认为任何复杂的波形都可以分解成不同正弦波,而不同的正弦波也能叠加成复杂波形。)

为了解决这个问题,引入自相关函数,它反映了随机信号本身在不同时刻的相互关系,再直白一点:把一个信号平移一段距离,跟原来有多相似。自相关函数将信号的蕴含的周期信号识别出来,并将相位信息去掉(相位不影响平均功率)。

维纳-辛钦定理:功率谱密度函数与自相关是傅里叶变换对。即一个信号的功率密度谱,就是其自相关函数的傅里叶变换。

总结:随机信号的幅值是满足正态分布的,用它的自相关函数求均方值即得到功率谱密度(PSD),一个随机的信号(随机振动)用PSD来描述其特征。

理论实践
传感器采集的样本为离散的点,没必要把离散点还原为原始信号,直接将离散的点用普估计方法求PSD。普估计方法由周期图法和welch方法,周期图法属于有偏估计,误差很大,所以采用welch方法,python或matlab有welch方法的库函数,可直接调用求得PSD。

python代码示例文章来源地址https://www.toymoban.com/news/detail-782084.html

#该加速度传感器的数据是三轴的X/Y/Z,默认Z轴为1g。
#返回四个值分别为频率和三个轴的PSD值,可用于绘图
def get_psd(payload):
	SAMPLES_NUM = 1024
    TicksInSamplingTimePeriod = int(0.5 + (0.0025 / (1 / 32768)))
    Fs_Acc = 32768 / TicksInSamplingTimePeriod  # 400 Hz
    NFFT_Acc = 512  # SAMPLES_NUM_PER_PIN    # 512
    Sensitivity_Acc = 1 / 1024  # LSB = 1/1024 [g] for the +/-2g range
    window = hamming(M=NFFT_Acc)

    samples = np.array(struct.unpack("<{}h".format(SAMPLES_NUM * 3), payload))#从payload中解析出来数据sample
    samples = np.array(samples) * Sensitivity_Acc
    acc_x_units_g, acc_y_units_g, acc_z_units_g = samples[0::3], samples[1::3], samples[2::3]
    (freq_Acc, Pxxf_Acc) = welch(x=acc_x_units_g, fs=Fs_Acc, window=window, noverlap=NFFT_Acc / 2,
                                 nfft=NFFT_Acc, return_onesided=True, detrend=False)
    Pxxf_AccX_dBg = 10. * np.log10(abs(Pxxf_Acc) + epsilon)
    (freq_Acc, Pxxf_Acc) = welch(x=acc_y_units_g, fs=Fs_Acc, window=window, noverlap=NFFT_Acc / 2,
                                 nfft=NFFT_Acc, return_onesided=True, detrend=False)
    Pxxf_AccY_dBg = 10. * np.log10(abs(Pxxf_Acc) + epsilon)
    (freq_Acc, Pxxf_Acc) = welch(x=acc_z_units_g, fs=Fs_Acc, window=window, noverlap=NFFT_Acc / 2,
                                 nfft=NFFT_Acc, return_onesided=True, detrend=False)
    Pxxf_AccZ_dBg = 10. * np.log10(abs(Pxxf_Acc) + epsilon)
    return freq_Acc, Pxxf_AccX_dBg, Pxxf_AccY_dBg, Pxxf_AccZ_dBg

到了这里,关于随机振动信号的特征——PSD(功率谱密度)的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 单片机项目分享 基于stm32的便携用电功率统计系统 -物联网 嵌入式 单片机

    🔥 这两年开始毕业设计和毕业答辩的要求和难度不断提升,传统的毕设题目缺少创新和亮点,往往达不到毕业答辩的要求,这两年不断有学弟学妹告诉学长自己做的项目系统达不到老师的要求。 为了大家能够顺利以及最少的精力通过毕设,学长分享优质毕业设计项目,今天

    2024年02月19日
    浏览(113)
  • stm32毕设分享 基于stm32的便携用电功率统计系统 -物联网 嵌入式 单片机

    🔥 这两年开始毕业设计和毕业答辩的要求和难度不断提升,传统的毕设题目缺少创新和亮点,往往达不到毕业答辩的要求,这两年不断有学弟学妹告诉学长自己做的项目系统达不到老师的要求。 为了大家能够顺利以及最少的精力通过毕设,学长分享优质毕业设计项目,今天

    2024年02月22日
    浏览(62)
  • 嵌入式QT (Qt 信号与槽)

    因为有了信号与槽的编程机制,在 Qt 中处理界面各个组件的交互操作时变得更加直观和简单。 信号 (Signal)就是在特定情况下被发射的事件。 GUI 程序设计的主要内容就是对界面上各组件的信号的响应,只需要知道什么情况下发射哪些信号,合理地去响应和处理这些信号就

    2024年02月02日
    浏览(51)
  • 嵌入式Qt Qt中的信号处理

    Qt消息模型: - Qt封装了具体操作系统的消息机制 - Qt遵循经典的GUI消息驱动事件模型 Qt中定义了与系统消息相关的概念;   Qt中的消息处理机制: Qt的核心 QObject::cinnect函数: Qt中的“新”: 实验1 初探信号与槽   运行结果: - 只有QObject的子类才能自定义槽 - 定义槽的

    2024年02月22日
    浏览(53)
  • 嵌入式实时操作系统的设计与开发(信号量学习)

    除了临界点机制、互斥量机制可实现临界资源的互斥访问外,信号量(Semaphore)是另一选择。 信号量与互斥量的区别 对于互斥量来说,主要应用于临界资源的互斥访问,并且能够有效地避免优先级反转问题。 对于信号量而言,它虽然也能用于临界资源的互斥访问,但是不能

    2024年02月08日
    浏览(59)
  • 嵌入式I2C 信号线为何加上拉电阻(图文并茂)

    IIC 是一个两线串行通信总线,包含一个 SCL 信号和 SDA 信号,SCL 是时钟信号,从主设备发出,SDA 是数据信号,是一个双向的,设备发送数据和接收数据都是通过 SDA 信号。   在设计 IIC 信号电路的时候我们会在 SCL 和 SDA 上加一个上拉电阻。   今天就来分享下,为什么要在

    2024年02月19日
    浏览(50)
  • lv6 嵌入式开发-Flappy bird项目(信号机制、定时器功能实现)

    目录 1 信号(signal) 2 设置信号响应方式 – signal 3 设置定时器 4 示例 问题: getch()阻塞获取键盘按键输入, 怎么操作才能不影响小鸟下落和管道移动? getch如果阻塞,下面的程序都是无法执行。通过信号机制方式实现。 在Linux中,软中断信号(signal,简称为信号)是在软件层

    2024年02月07日
    浏览(91)
  • 嵌入式操作教程_数字信号处理_音频编解码:3-6 AAC音频解码实验

    了解AAC音频格式,掌握AAC音频解码的原理,并实现将AAC格式的音频解码为PCM 音频编解码的主要对象是音乐和语音,音频的编解码格式可分为无压缩的格式、无损压缩格式、有损音乐压缩格式、有损语音压缩格式和合成算法。本实验中使用的AAC格式属于有损音乐压缩格式。音频

    2024年04月15日
    浏览(67)
  • 物联网选题分享 - 便携式STM32的红外信号分析仪 - 单片机 物联网 嵌入式

    🔥 这两年开始毕业设计和毕业答辩的要求和难度不断提升,传统的毕设题目缺少创新和亮点,往往达不到毕业答辩的要求,这两年不断有学弟学妹告诉学长自己做的项目系统达不到老师的要求。 为了大家能够顺利以及最少的精力通过毕设,学长分享优质毕业设计项目,今天

    2024年02月21日
    浏览(53)
  • 嵌入式教学实验箱_数字信号处理实验箱_操作教程:5-16 灰度图像线性变换(LCD显示)

    学习灰度图像线性变换的原理,掌握图像的读取方法,并实现在LCD上显示线性变换前后的图像。 一般成像系统只具有一定的亮度范围,亮度的最大值与最小值之比称为对比度。由于形成图像的系统亮度有限,常出现对比度不足的弊病,使人眼观看图像时视觉效果很差,通过灰

    2024年02月03日
    浏览(68)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包