#include "typedef.h"
#include "Convolve.h"


void Convolveok(
     double *Input,			//	原始輸入數(shù)據(jù)
     double *Impulse,			//	沖擊響應(yīng)
     double *Output,			//	卷積輸出結(jié)果
     Word16 length 			//  卷積序列長度
)

{
	int i,k,p;
	double r;

	p=0;
	for (k=0; k<=length-1; k++)
	{		
		Output[k]=0;
		r=0;
		for (i=0; i<=p; i++)
		{
			r=Input[k-i]*Impulse[i];
			Output[k] = Output[k]+r;
		}
			p=p+1;
			if (p>length-1) p=length-1;
			else p=p;
	 }
	
	p=length-2;
	for (k=length; k<=length+length-1; k++)
	{
	 	Output[k]=0;
		r=0;
	   	for (i=0; i<=p; i++)
		 {
		 	r = Input[length-1-i]*Impulse[length-1-p+i];
		 	Output[k] = Output[k]+r;
		 }
	    p=p-1;
	}
	return ;
}