1、整数版本
#include "corr.h"
int * corr(int a [] ,int b[], int la, int lb)
{
int *result;
result = (int *)malloc ((la+lb-1)*sizeof(int));
int i,j;
for (i = 0; i < la+lb-1; i++){
if(i<=lb){
for(j = 1; j <= i; j++){
result[i] = result[i] + a[j]*b[lb-i+j];
}
}
else if(i<=la){
for(j = i-lb+1; j <= i; j++){
result[i] = result[i] + a[j]*b[j-i+lb];
}
}
else{
for(j = i-lb+1; j <= la; j++){
result[i] = result[i] + a[j]*b[j-i+lb];
}
}
}
return result;
}
2、浮点数版本
#include "corr.h"
float * cor