求助一个滤波问题

小卿

自己利用了原子提供的DAC例程为基础做一个滤波器。由于硬件已经完成了，所以只能利用软件来做滤波。滤波的参数是用MATLAB产生的（IIR数字滤波）。也在晚上查过了很多资料，也试过一个认为网上认为行的代码，可以一敲完代码，编译就出错了。然后自己按着书上的公式写了一个，可是波形是不能实现滤波的。希望有谁也做过软件来滤波的帮帮忙看下是什么问题

main.c如附件所示

adlu

可以参考TI的DSPF_sp_iir( )函数实现方法，见附件P79-80。
TMS320C67x DSP Library Programmers Reference Guide.pdf (400.27 KB, 下载次数: 299)

2016-5-26 14:58 上传

点击文件名下载附件

void DSPF_sp_iir (float* restrict r1,
  const float* x,
  float* restrict r2,
  const float* h2,
  const float* h1,
  int nr
)
{
  int i, j;
  float sum;
  for (i = 0; i < nr; i++)
  {
    sum = h2[0] * x[4+i];
    for (j = 1; j <= 4; j++)
      sum += h2[j] * x[4+i&#8722;j] &#8722; h1[j] * r1[4+i&#8722;j];
    r1[4+i] = sum;
    r2 = r1[4+i];
  }
}
这个是处理数据块的，要实现对数据流的处理要自己改进。
MATLAB辅助设计可参考
http://blog.csdn.net/syrchina/article/details/9053039

xianshasaman

我也不太懂IIR滤波，但是你不把算式变换成差分式的话怎么在程序里用呢？

小卿

xianshasaman 发表于 2016-5-26 13:17
我也不太懂IIR滤波，但是你不把算式变换成差分式的话怎么在程序里用呢？

我的第一个截图里的第二个FOR里面就是做运算呀

xianshasaman

小卿 发表于 2016-5-26 13:49
我的第一个截图里的第二个FOR里面就是做运算呀

所以你就在单片机上对时域的数据直接做卷积了？
且不说你代码写到对不对，这么干肯定没有好结果啊

xianshasaman

xianshasaman 发表于 2016-5-26 14:34
所以你就在单片机上对时域的数据直接做卷积了？
且不说你代码写到对不对，这么干肯定没有好结果啊

啊，我看错啦，忽略我的话吧

小卿

xianshasaman 发表于 2016-5-26 14:41
啊，我看错啦，忽略我的话吧

单片机可以直接处理时域的值吧 我不断地取点出来就可以啦，何况公式也是对每一个点来做卷积运算的

小卿

adlu 发表于 2016-5-26 14:57
可以参考TI的DSPF_sp_iir( )函数实现方法，见附件P79-80。

好的  看看