Matlab全八度和三分之一八度分析

Question

我正在尝试在 Matlab 中实现倍频程和 1/3-倍频程滤波器。不幸的是，我无法访问 Acoustic 工具箱，但我尝试在 Matlab 中使用 fdesign.octave，但不确定我是否走在正确的道路上 - 并希望得到任何帮助/建议。

以下是我在 70 秒内获取的信号 (Fs=10,000Hz)：

Fs = 10000;            % Sampling frequency
T = 1/Fs;             % Sampling period
L = 700000;             % Length of signal
t = (0:L-1)*T;        % Time vector
S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
x = S + 2*randn(size(t));

三分之一倍频程的实现：

d = fdesign.octave(3,'Class 1','N,F0',6,11000,44100);
Hd = design(d);
y= filter(Hd,x);

原始信号和滤波信号图：

在这一点上，我有一个问题 - 如何将 fdesign.octave 生成的数据绘制在八度或三分之一八度频带中？

Answer 1

对于完整的实现，请使用 Matlab 文件交换中的this。无需声学工具箱。包含的演示应该可以解决您的绘图问题。