1/3倍频程代码与最大Z振级：批量处理、一键傻瓜式操作环境振动数据保存工具
在处理环境振动数据时，1/3倍频程分析是一个非常常见的需求。今天我们就来聊聊如何用Matlab
实现这个功能，顺便还能一键保存所有数据和图片，真正做到“傻瓜操作”。
首先，我们需要明确什么是1/3倍频程。简单来说，1/3倍频程就是将频率范围分成若干个1/3倍频带
，每个频带中心频率是前一个频带中心频率的2^(1/3)倍。这样做的目的是为了更好地分析振动信号在不
同频带上的分布情况。
接下来，我们来看一段Matlab代码，它可以实现1/3倍频程分析，并且可以批量输出各点的1/3倍频
程和最大Z振级随距离的衰减关系。
```matlab
% 假设我们有一个振动信号x，采样频率为fs
x = randn(1, 10000); % 随机生成一个振动信号
fs = 1000; % 采样频率为1000Hz
% 定义1/3倍频程的中心频率
f_center = [20 25 31.5 40 50 63 80 100 125 160 200 250 315 400 500 630 800 1000];
% 初始化存储1/3倍频程结果的数组
octave_band = zeros(1, length(f_center));
% 计算每个频带的1/3倍频程
for i = 1:length(f_center)
f_low = f_center(i) / (2^(1/6));
f_high = f_center(i) * (2^(1/6));
[b, a] = butter(2, [f_low f_high]/(fs/2));
y = filter(b, a, x);
octave_band(i) = rms(y);
end
% 输出1/3倍频程结果
disp('1/3倍频程结果：');
disp(octave_band);
% 计算最大Z振级
max_Z_level = max(octave_band);
disp('最大Z振级：');
