互相关函数实现（c语言完整工程，利用cmake编译）

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 * corr_f(float a [] ,float b[], int la, int lb)

{

		float *result;

		result = (float *)malloc ((la+lb-1)*sizeof(float));

		int i,j;

		for (i = 0; i < la+lb-1; i++){
			if(i&