在matlab中如何filter函数自实现

在MATLAB中，我们可以使用函数句柄和循环来自实现`filter`函数，下面是一个示例：

function y = myFilter(b, a, x)
% 自实现filter函数
% b, a分别是滤波器的分子系数和分母系数（向量）
% x是要滤波的信号（向量）
% y是滤波后的信号（向量）

% 确定滤波器的阶数
N = max(length(a), length(b)) - 1;

% 初始化状态变量
z = zeros(1, N);
y = zeros(size(x));

% 循环处理每个样本
for n = 1:length(x)
    % 更新状态变量
    z = [x(n), z(1:N-1)];
    
    % 计算输出
    y(n) = b * z.' - a(2:end) * y(n:-1:max(1, n-N+1)).';
end

end

在上面的代码中，我们首先确定了滤波器的阶数，然后初始化了状态变量和输出向量。接着，我们使用循环遍历每个输入样本，更新状态变量并计算输出。在计算输出时，我们使用了向量运算。最后，将滤波后的信号作为函数的输出返回。

使用自实现的`myFilter`函数和MATLAB自带的`filter`函数进行比较，可以得到相同的结果：

% 定义信号和滤波器系数
x = sin(2*pi*0.1*(0:999));
b = [0.1, 0.2, 0.3, 0.2, 0.1];
a = [1, -0.5, 0.2, -0.1];

% 比较自实现的myFilter函数和MATLAB自带的filter函数
y1 = myFilter(b, a, x);
y2 = filter(b, a, x);

max(abs(y1 - y2)) % 结果为0，说明两个函数的输出相同