matlab进行数据的后处理（信号处理）

mt917taylor

1.针对软件（如virtual.lab等等）计算出的位移或加速度或声压响应。假如计算频段宽度：10hz-2000hz，间隔10hz
如何将导出的数据利用matlab编程进行计算：
111.计算出整体频段内的均方差
222.如何将数据转化成1/3倍频程的数据
333.如何将将加速度数据转换成加速度功率谱密度形式
444.如何将数据处理成任意间隔频段的数据
谢谢大家

江树空

你这从VL出来的已经是频域数据了
第一个问题：所有数据平方后求和（即能量求和）开根号所得即为所关心的全频段内的RMS值；
第二个问题：计算出10Hz~2000Hz频段内的三分之一倍频程的每一个带宽的每一对上下限频率，然后在每一对上下限内进行能量求和再开根号，即得到1/3倍频程数据；
第三个问题：因为你的已经是频域数据了，只需将你的数据平方后除以你的频率分辨率10Hz即得到功率谱密度；
第四个问题：与第二个类似，任何一对上下限频率内都可以进行能量求和再开根号，此即为你任意间隔段的RMS值。1/n倍频程只不过是这中问题的某种特殊上下限对。

mt917taylor

江树空 发表于 2014-8-28 15:05
你这从VL出来的已经是频域数据了
第一个问题：所有数据平方后求和（即能量求和）开根号所得即为所关心的全 ...

Virtual里的采样频率是多少呀？
1.那由加速度功率谱密度曲线（采样频率已知，间隔1Hz），如何得到其1/3倍频程的数据呢？
2.如何由实测的加速度时域信号，得到其1/3倍频程的加速度功率谱密度数据呢？

江树空

mt917taylor 发表于 2014-8-29 10:28
Virtual里的采样频率是多少呀？
1.那由加速度功率谱密度曲线（采样频率已知，间隔1Hz），如何得到其1/3 ...

你从VL出来的既然是频域数据了，你还关心采样率作甚！你都说了关心的是10Hz~2000Hz，想必采样率不会低于4000Hz，至于具体是多少，看你在得到频域数据的过程中是怎么设置的了。。。
1、首先按照GB/T4431计算出10~2000Hz带宽内的1/3倍频程的上下限对，将你的功率谱密度谱线落在每一个上下限对之内的谱线进行求和，求和后开根号即可。有多少对上下限你计算多少次求和并开根号。
2、确切的讲，你第2个问题说法不准确，1/n倍频程数据不能叫功率谱密度（你这每一个带宽都不相同，何谓密度？建议弄清楚功率谱密度、1/n倍频程声压级、加速度级等名次的概念）。从时域信号得到1/3倍频程有两种方法，一种是FFT方法、一种是倍频程滤波器。前者即计算出功率谱或功率谱密度之后在带宽内求和得到，后者是计算出1/3的上下限对之后，对每一对上下限设置带通滤波器滤波得到。

江树空

江树空 发表于 2014-8-29 16:30
你从VL出来的既然是频域数据了，你还关心采样率作甚！你都说了关心的是10Hz~2000Hz，想必采样率不会低于4 ...

不好意思国标号搞错了，应该是GB/T3241.

gudanyuzhou

坐等大神指教