matlab频率采样法设计FIR滤波器,用频率采样法设计FIR滤波器

频率采样法是一种常用的FIR数字滤波器设计方法。下面给出一个MATLAB示例，使用频率采样法设计FIR滤波器：

假设我们要设计一个低通滤波器，其通带截止频率为0.4π，阻带截止频率为0.5π，通带最大衰减为1dB，阻带最小衰减为30dB。

首先，我们需要确定滤波器的采样频率。假设采样频率为1kHz，则对应的角频率为2π×1000=2000π。

接下来，我们需要确定滤波器的频率响应。根据设计要求，我们可以得到以下频率响应：

- 通带增益为1；
- 通带截止频率为0.4π，对应的角频率为0.8π；
- 阻带截止频率为0.5π，对应的角频率为π；
- 通带最大衰减为1dB，对应的增益为0.891；
- 阻带最小衰减为30dB，对应的增益为0.001。

根据上述频率响应，我们可以得到如下代码：

% 采样频率
Fs = 1000;

% 设计要求
wp = 0.4*pi;
ws = 0.5*pi;
Rp = 1;
Rs = 30;

% 计算通带、阻带截止频率对应的角频率
wp = wp*2*pi/Fs;
ws = ws*2*pi/Fs;

% 计算通带最大衰减、阻带最小衰减对应的增益
Rp = 10^(Rp/20);
Rs = 10^(Rs/20);

% 计算插值点数
N = ceil((8*pi)/(ws-wp));

% 计算插值点的频率
w = linspace(0,2*pi,N);

% 计算插值点的频率响应
Hd = zeros(1,N);
Hd(w<=wp) = 1;
Hd(w>=ws) = 0;
Hd(w>wp & w<ws) = (Rs-Rp)*0.5*(1-cos((w(w>wp & w<ws)-wp)*pi/(ws-wp)))+Rp;

% 进行反变换得到滤波器系数
h = ifft(Hd);
h = real(h);

% 绘制滤波器的频率响应
freqz(h);

运行上述代码后，即可得到该FIR滤波器的频率响应。

需要注意的是，频率采样法在设计FIR滤波器时需要进行插值，插值点数的选择会对滤波器的性能产生影响。一般来说，插值点数越多，滤波器的性能越好，但计算量也会增加。因此，插值点数的选择需要在设计要求和性能要求之间进行平衡。