c语言卷积算法

卷积计算的4个步骤：反折，平移，相乘，相加。将两组数，进行卷积后，变成第三组数。

//---------------------------------------------------------------------------
#include <vcl.h>
#pragma hdrstop
#include"stdio.h"
int  table_a[3]={4,2,3};//原数组A-给定
int  table_b[3]={2,5,1};//原数组B-给定
//---------------------------------------------------------------------------
#pragma argsused
void My_Convolution(int *Tab_A,int *Tab_B)
{
	int temp = 0;
	int i = 0;
	int j = 0;
	int k = 0;
	int Length_A,Length_B;//数组A和数组B的长度
	int ALL_Length;//卷积后的总长度
    unsigned int Temp_A[10000]={0};//此空间要足够大，要满足卷积后的总长度
    unsigned int Conv_out[10000]={0};//用于存放输出结果的数组

	Length_A = sizeof(table_a)/sizeof(unsigned int);//求数组A的长度
	Length_B = sizeof(table_b)/sizeof(unsigned int);//求数组B的长度
	printf("数组A的长度: %d\n",Length_A);//打印A的长度
	printf("数在B的长度: %d\n",Length_B);//打印B的长度
	ALL_Length = Length_A + Length_B -1;//卷积后的总长度
	printf("卷积结果长度: %d\n",ALL_Length);//打印卷积结果长度
	for(i = 0; i < Length_A; i++)//将输入参数Tab_A数组后面补0
	{
		Temp_A[i] += Tab_A[i];//处理后放到临时数组Temp_A中
	}
	printf("The convolution result is: \n");//输出提示，"卷积结果如下："
	//------------卷积核心算法部分---------------
	for(i = 0; i < ALL_Length; i ++)
	{
		for(k = i,j = 0; k >= 0; k --,j ++)
		{
			temp += Temp_A[k]*Tab_B[j];	
		}
		Conv_out[i]=temp;
		temp=0;
		printf("%d\n",Conv_out[i]);
	}
    //------------------------------------------
}
int main(int argc, char* argv[])
{
    My_Convolution(table_a,table_b);
    getchar();
    return 0;
}
//---------------------------------------------------------------------------