matlab模拟滤波器频率响应,怎样求FIR滤波器在任意频率上的频率响应？

如题，本人设计一个带通滤波器，想计算得到滤波器在任意频点上的脉冲响应。有三种方法：1.利用freqz计算滤波器的频率响应；2.对滤波器系数进行FFT变换。上述两种方法得到的都是离散频点的脉冲响应，频率是固定的离散值，不能满足满足需求。3.计算滤波器的转移函数。方法3代码如下：

clc;

close all;

clear all;

NFFT             = 2*8192; % FFT点数

fs                   = 1250;    % Sample rate，采样速率

Ts                   =1/fs;

%% 50Hz带通滤波处理

f0       = 50;% Hz，带通中心频率

BW     = 2*4; % 通带宽度

f_p1    = f0 - BW/2; %通带左边缘频率，模拟频率, p--pass

f_p2    = f0 + BW/2; %通带右边缘频率，模拟频率

deta1 = 1;                                       % 通带峰值起伏最大值，dB

f_s1     = f_p1 - 2;  %阻带左边缘频率，模拟频率, s--stop

f_s2     = f_p2 + 2; %阻带右边缘频率，模拟频率

delta2 = 40;                                 % 阻带最小衰减，dB

fcuts                       = [f_s1,f_p