学习记录——倍频程、中心频率的理解

还是《从这里学NVH——噪声、振动、模态分析的入门与进阶》这本书。

我们由傅里叶变换得到频谱之后,横坐标就是频率,如果我们等间隔的选取区间,由于人的听觉范围是20Hz——20KHz,是很大的数据量,处理起来很麻烦,所以我们不选取等宽频程得到做法,而是选取倍频程。

接下来就是理解什么是倍频程——就是对于相邻的两个频宽来说,后一个的频宽是前一个的倍数,类似于1、2、4、8这种关系。值得一提的是,我们一般说的倍频程,其实是1/1倍频程,那我们常见的1/3倍频程,就是把每一个1/1倍频程又分为了三部分。

下面这个图可以帮助我们理解。

在这里插入图片描述

参考链接:

倍频程的介绍

倍频的概念

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
对振动加速度信号按照1/3倍频程中心频率进行处理,需要进行如下步骤: 1. 读取振动加速度信号数据并选择感兴趣的时间段 2. 对信号进行1/3倍频程滤波 3. 计算有效加速度 下面是具体实现步骤: 1. 读取振动加速度信号数据并选择感兴趣的时间段 假设振动加速度信号数据保存在名为"acceleration_data.csv"的文件中,可以使用MATLAB中的csvread函数读取数据。然后,可以选择感兴趣的时间段,例如从第1000个采样点开始到第2000个采样点结束。代码如下: ``` data = csvread('acceleration_data.csv'); % 读取数据 t = data(:, 1); % 获取时间向量 a = data(:, 2); % 获取加速度向量 % 选择感兴趣的时间段 t_interest = t(1000:2000); a_interest = a(1000:2000); ``` 2. 对信号进行1/3倍频程滤波 可以使用MATLAB中的fir1函数设计一个1/3倍频程滤波器,并使用filtfilt函数对信号进行滤波。代码如下: ``` fs = 1000; % 采样频率 fc = fs/6.3/3; % 1/3倍频程中心频率 [b, a] = fir1(100, fc/(fs/2)); % 设计滤波器 % 对信号进行滤波 a_filtered = filtfilt(b, a, a_interest); ``` 3. 计算有效加速度 可以使用rms函数计算经过滤波后的信号的有效加速度。代码如下: ``` a_rms = rms(a_filtered); ``` 完整的MATLAB代码如下: ``` data = csvread('acceleration_data.csv'); % 读取数据 t = data(:, 1); % 获取时间向量 a = data(:, 2); % 获取加速度向量 % 选择感兴趣的时间段 t_interest = t(1000:2000); a_interest = a(1000:2000); fs = 1000; % 采样频率 fc = fs/6.3/3; % 1/3倍频程中心频率 [b, a] = fir1(100, fc/(fs/2)); % 设计滤波器 % 对信号进行滤波 a_filtered = filtfilt(b, a, a_interest); % 计算有效加速度 a_rms = rms(a_filtered); ``` 这段代码可以读取名为"acceleration_data.csv"的振动加速度信号数据,选择从第1000个采样点到第2000个采样点的时间段,并对该时间段内的信号进行1/3倍频程滤波,最后计算滤波后信号的有效加速度。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值