卷积函数的C语言实现
void Convert(double BMWL,double NORMAL,double *BMwindowVct,float *RR,int * rR,int lenz,double *RRproc)
{
?double *X=NULL;
?double *Y=NULL;
?if (NULL==X)
?{
??X=new double[BMWL];?
?}
?if (NULL==Y)
?{
??Y=new double[NORMAL];
?}
?for(int i=0;i<BMWL;i++)
?{
??X[i]=BMwindowVct[i];
?}
?for(int i=0;i<NORMAL;i++)
?{
??Y[i]=RR[(int)rR[i]];
?}
?double *zptr=NULL;
?double s;
?double *xp=NULL;
?double *yp=NULL;
?int n,n_lo,n_hi;
?double *Z=NULL;
?if (NULL==Z)
?{
??Z=new double[lenz];
?}
?zptr=Z;
?for (int i=0;i<lenz;i++)
?{
??s=0.0;
??n_lo=0>(i-NORMAL+1)?0:i-NORMAL+1;
??n_hi=BMWL-1<i?BMWL-1:i;
??xp=X+n_lo;
??yp=Y+i-n_lo;
??for (n=n_lo;n<=n_hi;n++)
??{
???s+=*xp * *yp;
???xp++;
???yp--;
??}
??*zptr=s;
??zptr++;
?}
?for(int i=0;i<lenz;i++)
?{
??RRproc[i]=Z[i];
?}
??delete [] X;
??X=NULL;
??delete [] Y;
??Y=NULL;
??delete [] Z;
??Z=NULL;
?int m=0;
?for(int l=(ceil((double)BMWL/2-2));l<(ceil((double)BMWL/2)+NORMAL-2);l++)
?{
??RRproc[m]=RRproc[l];
??m++;
?}
}
