fir1函数的具体程序

下面是一个使用fir1函数设计低通滤波器的MATLAB程序示例：

% 设计一个10阶低通滤波器，截止频率为1000Hz，采样频率为8000Hz
n = 10;
Wn = 1000 / (8000 / 2);

% 计算窗函数
w = hamming(n+1);

% 计算理想低通滤波器的系数
b_ideal = fir1(n, Wn);

% 计算实际低通滤波器的系数
b = b_ideal .* w';

% 对滤波器系数进行归一化
b = b / b(1);

% 绘制滤波器的幅频响应曲线
freqz(b, 1, 1024, 8000);

这个程序使用fir1函数设计一个10阶低通滤波器，截止频率为1000Hz，采样频率为8000Hz。程序首先计算出窗函数，然后通过fir1函数计算出理想低通滤波器的系数向量b_ideal。接下来，将理想低通滤波器的系数向量与窗函数相乘，得到实际低通滤波器的系数向量b。最后，通过对滤波器系数向量进行归一化处理，得到归一化后的滤波器系数向量。程序最后使用freqz函数绘制出滤波器的幅频响应曲线。