简易FIR滤波器

代码说明：

1、FIR数字滤波器的设计（1）在主函数中调用dsnFIR函数，滤波器参数设为21阶、带阻、矩形窗、采样率1000Hz、上下截止频率分别为100Hz,400Hz；（2）利用drawFIR函数分别画出所设计的FIR滤波器的幅频特性图和相频特性图。2、FIR数字滤波实现（1）编写一个输入信号子程序func()，输入信号为一连续函数，表达式为；（2）对输入信号进行采样，即       ，这里取采样率，得到离散序列x(n)；（3）取x(n)中前40个点的值，将其转为复数存入复数数组x[64]中，并对x[64]做64点的FFT，得到的X(k)存回x[64]中；（4）已知在主函数中调用dsnFIR函数后得到的输出结果b[21]即为FIR滤波器单位冲激响应h(n)的值，将其转为复数存入复数数组h[64]中，并对h[64]做64点的FFT，得到的H(k)存回h[64]中；（5）已知用FIR滤波器对输入信号进行滤波就是将输入信号x(n)与滤波器单位冲激响应h(n)做卷积运算，故我转到频域采用FFT方法实现对输入信号x(n)的滤波，即将X(k)与H(k)相乘，得到滤波后信号的FFT变换Y(k)，然后只需对Y(k)做IFFT变换即可得到滤波后信号y(n)。