用matlab实现，采样率为0.96MHZ，将信号经过抽取后，通过多相滤波器，得到20个信道，信道划分方式采用50%交叠，每个子信道带宽48KHz，得到20个信道的频谱

以下是一个 Matlab 示例代码，实现信号抽取、多相滤波器和频谱绘制的功能：

% 假设原始信号为 x，采样率为 fs
fs = 960000;
t = 0:1/fs:1-1/fs;
x = sin(2*pi*1000*t);

% 设置抽取因子和滤波器参数
R = 20; % 抽取因子
M = 101; % 滤波器长度
L = 20; % 信道数
f_cutoff = 24000; % 截止频率
h = fir1(M-1, f_cutoff/(fs/2)); % 低通滤波器系数

% 信号抽取
x_up = zeros(1, length(x)*R);
x_up(1:R:end) = x;

% 多相滤波器
y = zeros(L, length(x_up)/R);
for k = 1:L
    y(k,:) = filter(h, 1, x_up.*exp(-1j*2*pi*(k-1)*(0:length(x_up)-1)/L));
end

% 计算频谱并绘图
f = linspace(0, fs/2, M*100);
H = zeros(L, length(f));
for k = 1:L
    H(k,:) = freqz(y(k,:), 1, f, fs);
end
figure;
plot(f/1000, abs(H));
xlabel('Frequency (kHz)');
ylabel('Magnitude');
title('Filter Bank Frequency Response');

该示例代码假设原始信号为 1 秒钟的 1000 Hz 正弦波，采样率为 960000 Hz。抽取因子为 20，滤波器长度为 101，信道数为 20，信道划分方式采用 50% 交叠，每个子信道带宽为 48 kHz。通过多相滤波器将原始信号分成 20 个子信道，计算每个子信道的频谱并绘制出来。注意，滤波器的截止频率应为实际带宽除以抽取因子，即 24000 Hz。