卷积函数的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;
}
if (NULL==Y)
{
Y=new double;
}
for(int i=0;i<BMWL;i++)
{
X=BMwindowVct;
}
for(int i=0;i<NORMAL;i++)
{
Y=RR[(int)rR];
}
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;
}
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=Z;
}
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=RRproc;
m++;
}
}
页:
[1]