matlab 自己编写filter函数

Matlab 中可以自定义滤波器函数，通常通过设计数字信号处理(DSP)滤波器结构如 FIR (有限 impulse response)、IIR (无限 impulse response) 或者窗函数等来实现。以下是创建简单 FIR 滤波器的基本步骤：

1. **确定滤波器类型和参数**：你需要明确滤波器的类型（低通、高通、带通、带阻等），以及频率响应、截止频率、采样率等参数。

2. **设计滤波器系数**：
- 对于FIR滤波器，你可以使用`fir1`, `fir2`, `firwin`, 或 `firpm` 等函数来直接设计线性相位滤波器，或者使用`designfilt` 函数设计更复杂的滤波器。
- 对于IIR滤波器，可以使用`remez`, `butter`, `cheby1`, `cheby2`, `ellip`等函数来设计 Butterworth, Chebyshev Type I, II, 或 Elliptic 滤波器，然后用`tf2ss`将传递函数转换为状态空间形式。

3. **编写滤波器函数**：

function y = myCustomFilter(x, taps, isIIR)
    % x: 输入信号
    % taps: 如果isIIR为true，taps是零极点对数组；如果是FIR，它是滤波器系数向量
    % isIIR: true表示IIR滤波器，false表示FIR

    if isIIR
        [num, den] = tf(taps);
        sys = ss(num, den, zeros(size(num)), ones(1, size(den)));
        y = lsim(sys, x);  % 使用连续时间模拟
    else
        y = filter(taps, 1, x);  % FIR滤波器，1是单位阶跃响应
    end
end

4. **测试和应用滤波器**：
- 传入输入信号 `x` 和相应的滤波器参数到自定义函数中，即可得到过滤后的信号 `y`。