给我展示一段2路并行多相滤波器的matlab代码

以下是一段2路并行多相滤波器的MATLAB代码：

% 2路并行多相滤波器设计
M = 32; % 每个通道的滤波器长度
L = 2; % 并行通道数
Fs = 1000; % 采样率
Fpass = 100; % 通带边界频率
Fstop = 200; % 阻带边界频率
Rp = 1; % 通带最大衰减
Rs = 30; % 阻带最小衰减
Wp = Fpass / (Fs/2); % 归一化通带边界频率
Ws = Fstop / (Fs/2); % 归一化阻带边界频率
[N,Wn] = buttord(Wp, Ws, Rp, Rs); % 计算巴特沃斯滤波器阶数和截止频率
[b,a] = butter(N,Wn); % 设计巴特沃斯滤波器
h = impz(b,a,M); % 计算滤波器系数

% 多相滤波器分解
h0 = h(1:L:M); % 通道0的滤波器系数
h1 = h(2:L:M); % 通道1的滤波器系数
h2 = h(3:L:M); % 通道2的滤波器系数
h3 = h(4:L:M); % 通道3的滤波器系数

% 生成输入信号
x = randn(1000,1);

% 2路并行多相滤波器滤波
y0 = conv(h0,x); % 通道0输出
y1 = conv(h1,x); % 通道1输出
y2 = conv(h2,x); % 通道2输出
y3 = conv(h3,x); % 通道3输出

% 合并通道输出
y = [y0(1:2:end) + y1(1:2:end), y2(1:2:end) + y3(1:2:end)];

% 绘制输入信号和输出信号图形
t = (0:length(x)-1)/Fs;
figure;
subplot(2,1,1);
plot(t,x);
title('输入信号');
xlabel('时间（秒）');
ylabel('幅值');
subplot(2,1,2);
plot(t,y);
title('输出信号');
xlabel('时间（秒）');
ylabel('幅值');

该代码实现了一个2路并行多相滤波器的设计，并将其应用于一个随机输入信号上，最终输出合并后的两个通道的滤波结果。