Matlab全倍频程和三分之一倍频程分析

我试图在Matlab中实现倍频程和1/3倍频程带滤波器。不幸的是，我没有访问声学工具箱的权限，但我尝试过在Matlab中使用fdesign.octave，但不确定我是否走上了正确的道路，并希望得到任何帮助/建议

以下是我的信号（Fs=10000Hz），持续70秒：

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在倍频程或三分之一倍频程频带中生成的数据？

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