matlab 这两种滤波器代码，分别写出

好的，以下是在 MATLAB 中实现巴特沃斯低通滤波器和递归巴特沃斯低通滤波器的代码示例：

巴特沃斯低通滤波器：

% 滤波器参数
fc = 1000;  % 截止频率
fs = 8000;  % 采样率
order = 4;  % 滤波器阶数

% 计算滤波器系数
[b,a] = butter(order, fc/(fs/2), 'low');

% 生成测试信号
t = 0:1/fs:1;
x = sin(2*pi*1000*t);

% 应用滤波器
y = filter(b,a,x);

% 绘制结果
subplot(211);
plot(t,x);
title('原始信号');
xlabel('时间 (s)');
ylabel('幅值');
subplot(212);
plot(t,y);
title('滤波后信号');
xlabel('时间 (s)');
ylabel('幅值');

递归巴特沃斯低通滤波器：

% 滤波器参数
fc = 1000;  % 截止频率
fs = 8000;  % 采样率
order = 2;  % 滤波器阶数
G = 1;  % 放大倍数

% 计算滤波器系数
K = tan(pi * fc / fs);
Q = 1 / sqrt(2);
b0 = 1 / (1 + K/Q + K^2);
b1 = 2 * b0;
b2 = b0;
a1 = 2 * b0 * (K^2 - 1) / (1 + K/Q + K^2);
a2 = (1 - K/Q + K^2) / (1 + K/Q + K^2);

% 生成测试信号
t = 0:1/fs:1;
x = sin(2*pi*1000*t);

% 应用滤波器
y = zeros(size(x));
for n = 3:length(x)
    y(n) = G * (b0*x(n) + b1*x(n-1) + b2*x(n-2) - a1*y(n-1) - a2*y(n-2));
end

% 绘制结果
subplot(211);
plot(t,x);
title('原始信号');
xlabel('时间 (s)');
ylabel('幅值');
subplot(212);
plot(t,y);
title('滤波后信号');
xlabel('时间 (s)');
ylabel('幅值');

以上是简单的巴特沃斯低通滤波器和递归巴特沃斯低通滤波器的实现，希望对你有帮助。需要注意的是，在实际应用中，需要根据实际需求和性能要求来选择合适的滤波器类型和参数。