FIR滤波例程中，为什么滤波输出的数据尾部会有一个尖脉...

lushidegreen

创龙FIR例程的函数如下，为了方便观察，我把表达式中的100Hz改为了10Hz
void FIRTest(void)
{
        // 产生待测试信号 100hz信号 450hz噪声
        unsigned int i;
        for (i=0;i<Tn;i++)
                FIR_In=5*sin(2*PI*10*(i/Fs))+15*sin(2*PI*450*(i/Fs));   //100Hz改为了10Hz

        // FIR 滤波
         DSPF_sp_fir_r2(FIR_In,B,FIR_Out,N,Tn);
}

输入FIR_In就是一个10Hz的正弦信号叠加一个450Hz的正弦信号。

lushidegreen

运行完FIRTest( )函数，绘制输入输出信号波形：
输入FIR_In波形：



输出FIR_Out波形：

lushidegreen

问题：
1，为什么输出的波形尾部会带有一个尖峰脉冲？
2，从程序看，DSPF_sp_fir_r2( )函数用作数据块滤波了，是否合理？
3，如果我要进行音频数据滤波，能一块一块数据进行处理吗？

lushidegreen

在创龙《开发例程使用手册》的FIR例程描述中，没有看到上述的尖峰信号，是因为在绘图设置的时候把Display Data Size参数设置位100，即只显示100个数据，所以看到的FIR_Out波形似乎是完美的。

如果把该参数设置成1024，即显示所有的FIR_Out数据，会发现波形尾部是不一样的。

weibo1988821

lushidegreen 发表于 2015-7-10 10:38
在创龙《开发例程使用手册》的FIR例程描述中，没有看到上述的尖峰信号，是因为在绘图设置的时候把Display D ...

你考虑到了滤波器系数延迟吗？和你Matlab仿真一致吗？

lushidegreen

经过摸索，终于知道怎么回事了。
FIR滤波器的延迟就是它的长度N。
经实验验证，DSPF_sp_fir_r2( )函数的输出数据长度为Tn，而输入数据长度其实是Tn+N。
验证方法如下：
1.把输入信号数组长度扩展到Tn+N，即float FIR_In[Tn+N].
2.测试函数改为如下，
void FIRTest(void)
{
        // 产生待测试信号 10hz信号 450hz噪声
        unsigned int i;
        for (i=0;i<Tn+N;i++)
                FIR_In=5*sin(2*PI*10*(i/Fs))+15*sin(2*PI*450*(i/Fs));
        // FIR 滤波
         DSPF_sp_fir_r2(FIR_In,B,FIR_Out,N,Tn);
}

经过上述两步后，重新编译运行，画出的FIR_Out波形就是平滑的了。