SigAnlz.c

#include <stdio.h>
#include <math.h>
#include <float.h>

#include "SigAnlz.h"
#include "GaussCoef_SigAnlz.h"

static int constrain(int n, int min, int max){
	if(n<min)
		return min;
	if(n>max)
		return max;
	return n;
}

void threshold_SigAnlz(float* sig, int n, BitVec *bvptr, float k, float diff){
	float imean, irms, istddev, istddevPrev;
	
	istddev = FLT_MAX;	
	do{
		istddevPrev = istddev;
		param_SigAnlz(sig, n, bvptr, &imean, &irms, &istddev);

		for(int i=0; i<n; i++)
			if( sig[i]>imean+k*istddev || sig[i]<imean-k*istddev)
				setBit_BitVec(bvptr, i, 1);
	}while(istddevPrev - istddev > diff);
}

void param_SigAnlz(float *sig, int n, BitVec *bvptr, float *mean, float *rms, float *stddev){
	float tmean=0, trms=0, tstddev=0;
	int nSet=0;

	for(int i=0; i<n; i++){
		if(getBit_BitVec(bvptr,i)==0){
			tmean += sig[i];
			trms += sig[i]*sig[i];
			nSet++;
		}
	}

	tmean /= nSet;
	trms = sqrtf(trms/nSet);

	for(int i=0; i<n; i++){
		if(getBit_BitVec(bvptr,i)==0){
			tstddev += (sig[i]-tmean)*(sig[i]-tmean);
		}
	}

	tstddev = sqrtf(tstddev/nSet);

	*mean = tmean;
	*rms = trms;
	*stddev = tstddev;
}



void boxBlur_SigAnlz(float* ip, float *op, int n, int rad){
	float val;
	for(int i=0; i<n; i++){
		val = 0;
		for(int j=i-rad; j<=i+rad; j++){
			val += ip[constrain(j,0,n-1)];
		}
		op[i] = val/(2*rad+1);
	}
}

void gaussBlur_SigAnlz(float* ip, float *op, int n, int rad){
	for(int i=0; i<n; i++){
		op[i] = 0;
		for(int j=i-rad; j<=i+rad; j++){
			op[i] += ip[constrain(j,0,n-1)]*GaussCoeff[rad][j-i+rad];
		}
	}
}



void ignoreDipBitVec_SigAnlz(BitVec* bv, float *sig, int n, float k){
	float mean, rms, stddev;
	param_SigAnlz(sig, n, bv, &mean, &rms, &stddev);

	BitVec *tbv = new_BitVec(getSize_BitVec(bv));

	for(int i=0; i<n; i++){
		if(getBit_BitVec(bv,i)==0)
			setBit_BitVec(tbv,i,0);
		else if(sig[i]<mean-k*stddev)
			setBit_BitVec(tbv,i,0);
		else
			setBit_BitVec(tbv,i,1);
	}

	copy_BitVec(tbv, bv);
	free_BitVec(tbv);
}