power_law_tv.c

最大似然估计算法

#include "timeseries.h"

void power_law_tv(double *tv, int n, double degree) {

	int i, j;

	double *b;

	b = (double *)calloc((size_t)n, sizeof(double));

        b[0] = 1.0;

        for (j = 0; j < n-1; j++)  b[j+1] = ((double) j - degree) / ((double) (j+1));

        tv[0] = b[0];
        for (j = 1; j < n; j++) tv[j] = tv[j-1] * b[j];

	free(b);
}