Kaiser窗结合filtfilt函数设计零相位滤波器的MATLAB代码

以下是使用Kaiser窗设计零相位滤波器的MATLAB代码：

% 设计滤波器
Fs = 1000;          % 采样率
Fpass = 100;        % 通带截止频率
Fstop = 200;        % 阻带截止频率
Rp = 1;             % 通带最大衰减
Rs = 40;            % 阻带最小衰减
Wpass = Fpass/(Fs/2);
Wstop = Fstop/(Fs/2);
[N,Wn] = kaiserord([Wpass Wstop],[1 0],[Rp Rs]);
b = fir1(N,Wn,kaiser(N+1,beta),'noscale');

% 零相位滤波
y = filtfilt(b,1,x);

其中，`kaiserord`函数用于计算Kaiser窗参数，`fir1`函数用于生成滤波器系数，`kaiser`函数用于生成Kaiser窗，`filtfilt`函数用于进行零相位滤波。需要注意的是，上述代码中的变量`beta`是Kaiser窗的形状参数，可以根据需要调整。