关于数字信号处理课程设计 线性卷积在DSP芯片上的实现

#include <math.h>
#include <stdio.h>

#define Nx a //序列x[k]的长度
#define Nh b //序列h[k]的长度
#define Length Nx+Nh-1 //卷积或相关结果序列的长度

int x[Nx]={};
int h[Nh]={};
int y[Length];

void convolution(int*, int*, int*);

void main()
{
int k;
//调用卷积程序
convolution(x, h, y);
for(k=0;k<=Length-1;k++)
fprintf(stdout,"%i ",y[k]); // 打印卷积计算结果
fprintf(stdout,"\n ");

}

void convolution(int x[Nx] , int h[Nh], int y[Length] )
{
int i,k;
for (k=0;k<=Length-1; k++)
{
int temp=0;
for(i=0; i<Nx; i++)
{ if ((k-i)>=0 && (k-i) <=Nh-1)
temp=temp+x[i]*h[k-i];
}
y[k]=temp;
}
}

这就是卷积的C语言程序，将序列和序列的长度补充进去就好了

算法是很简单的哇，而且很多DSP会针对卷积，相关等操作来做指令优化的，优化的好写出来比C语言还少，跑得也快。