Критерий сериальной корреляции

[ilyakonst]

Помогите найти ошибку. Числа генерируются с помощью ЛКГ. Нужно получить "С" - коэффициент корреляции. Сама формула в прикрепленном файле. U, V - элементы двух последовательностей. Их начальный элемент(объявлены до цикла for - различны), т.е последовательности будут различны. Коэффициент должен получаться от -1 до 1, у меня же выводит -1.#IO

void korel()
{
int i,n=10;
double C=0, koren;
long long c1 = 0, c2 = 0, c3 = 0, z1 = 0, z2 = 0, z3 = 0, z4 = 0, chisl, znam, U=0, V=0;
U = 212;
V = 76;
for (i = 0; i < n; i++) {

U = (1073741829 * U + 17) % 4294967296;
//printf("%d\n", U);
V = (1073741829 * V + 17) % 4294967296;
//printf("%d\n", V);
c1 += U*V;
c2 += U;
c3 += V;
z1 += U*U;
z2 += U;
z3 += V*V;
z4 += V;
}
chisl = n*c1 - c2*c3;
znam = (n*z1 - z2*z2)*(n*z3 - z4*z4);
koren = sqrt(double(znam));
C = chisl / koren;
printf("%.3f",C);

}

[kutani]

Прицеп потеряли, вот критерий сериальной корреляции из Кнута:


всегда в ЦОС делал по нему )

[f.hump]

может znam у тебя отрицательный получается?

кстати, поменяй % 4294967296 на & 0x0FFFFFFFF, а то меня напрягает немного.

kutani, твоя формула выглядит как самокорреляция, а у тс взаимная корреляция двух полследовательностей. и формульной ошибки я там не вижу.

[kutani]

Оставлю это здесь http://computersbooks.net/index.php?...=20012&page=45

[ilyakonst]

Сообщение от kutani

Прицеп потеряли, вот критерий сериальной корреляции из Кнута:


всегда в ЦОС делал по нему )

нет, моя формула также из Кнута, а то, что вы указали - частный случай моей, который мне не подходит. Прочитайте внимательно условие ее использования

[ilyakonst]

Сообщение от f.hump

может znam у тебя отрицательный получается?

кстати, поменяй % 4294967296 на & 0x0FFFFFFFF, а то меня напрягает немного.

kutani, твоя формула выглядит как самокорреляция, а у тс взаимная корреляция двух полследовательностей. и формульной ошибки я там не вижу.

в каком именно месте, а то я нубас) и есть ли принципиальная разница? я же мог 2^32 записать через pow, но это лишняя операция. Я думаю, что у меня где-то происходит переполнение, но найти не могу

[ilyakonst]

Сообщение от f.hump

может znam у тебя отрицательный получается?

не должно быть, в формуле все в квадрат возводится, не может быть отрицательных чисел.
UPD: знаменатель действительно отрицательный, хотя это невозможно.

void korel()
{
int i,n=10;
double C=0, koren;
long long c1 = 0, c2 = 0, c3 = 0, z1 = 0, z2 = 0, z3 = 0, z4 = 0, chisl, znam, U=0, V=0;
U = 2102;
V = 76;
for (i = 0; i < n; i++) {

U = (1073741829 * U + 17) % 4294967296;
//printf("%d\n", U);
V = (1073741829 * V + 17) % 4294967296;
//printf("%d\n", V);
c1 += U*V;
c2 += U;
c3 += V;
z1 += U*U;
z2 += U;
z3 += V*V;
z4 += V;

}
printf("%d\n", z1); printf("%d\n", z2); printf("%d\n", z3); printf("%d\n", z4); //добавил вывод частей знаменателя
chisl = n*c1 - c2*c3;
znam = (n*z1 - z2*z2)*(n*z3 - z4*z4);
//printf("%lld\n", znam);
koren = sqrt(double(znam));
C = chisl / koren;
printf("%.3f",C);

}

Добавил вывод, там при разных входных параметрах числа, возведенные в квадрат почему-то получается отрицательными.
UPD: изменил количество генерируемых эл-ов, теперь оно четное, всегда получается коэффициент. Но он должен быть в пределах от -1..1, а у меня 10ки тыясч, значит где-то в формуле идет ошибка

[Аватар]

c1 += U*V;

Почему бы здесь не быть переполнению, или ниже по аналогичному

[ilyakonst]

Сообщение от Аватар

Почему бы здесь не быть переполнению, или ниже по аналогичному

проверил, переполнений нет нигде.
все высчитывается

[f.hump]

Аватар прав.
у тебя там в каждой операции целочисленные переполнения.

сделай
c1 = 0, c2 = 0, c3 = 0, z1 = 0, z2 = 0, z3 = 0, z4 = 0, chisl, znam
даблом

и в аккумуляторах делай каст (double)U, (double)V