Как стать автором
Обновить

Элементарный счет звездного года (365 дней и 369 минут [365.2524])в радиоактивном распаде

Python *Алгоритмы *Математика *
Recovery mode
def dt2(x,y):
    uu=0;
    for i in range(len(x)-3):
     uu1=abs((x[i]/x[i+1])-(y[i]/y[i+1]));
     uu2=abs((x[i]/x[i+2])-(y[i]/y[i+2]));
     uu3=abs((x[i]/x[i+3])-(y[i]/y[i+3]));
     uu+=uu1+uu2+uu3;
    return uu;

Есть данные за 2 дня мая 2005 года 2 дня мая 2006 года. Цель найти в сумме 1440 сравнений[60*24] звездный год.

Подготовка ввода данных радиоактивного распада.

import math
import random
from statistics import mean
koeff=60;
matrix = [line.strip() for line in open('/users/andrejeremcuk/downloads/0505.txt')];
matrix1 = [line.strip() for line in open('/users/andrejeremcuk/downloads/0506.txt')];
matrixd = [line.strip() for line in open('/users/andrejeremcuk/downloads/days.dat')];


arra=[[ 0 for e in range(koeff)] for t in range(int(len(matrix)/koeff))];
arra1=[[ 0 for e in range(koeff)] for t in range(int(len(matrix)/koeff))];
arrad=[[ 0 for e in range(koeff)] for t in range(int(len(matrixd)/koeff))];

for i in range(len(matrix)): matrix[i]=int(matrix[i]);matrix1[i]=int(matrix1[i]);

for i in range(len(matrixd)): matrixd[i]=int(matrixd[i]);

z=0;
for jk in range(int(len(matrix)/koeff)): #int(len(matrix)/koeff)
 for mk in range(koeff): arra[jk][mk]=float(matrix[z]);arra1[jk][mk]=float(matrix1[z]);z=z+1;

z=0;
for jk in range(int(len(matrixd)/koeff)): 
 for mk in range(koeff): arrad[jk][mk]=float(matrixd[z]);z=z+1;
def distt(rast):
    statistic=0;
    for jk in range(1440):
     x=arra[jk];y=arra1[jk+rast];
     statistic+=dt2(x,y);
    return statistic;

Результаты выполнения кода - на гитхабе подсчитаны все значения от +1 минуты до + 1440 минут сдвиг 2005 года от 2006 года.

>>> distt(369)
30142.867396646434
>>> distt(368)
29888.11808805372-Наименьшее значение из всех сдвигов от 0 до +1440 минут
>>> distt(367)
30128.959797342242
>>> distt(366)
30139.531482650476
>>> distt(365)
30080.92313005219
>>> distt(364)
30119.64335915869
>>> distt(363)
30192.614068153922
>>> distt(362)
30054.03651908298
>>> distt(361)
30053.5344380305
>>> distt(360)
30019.741597642616
>>> distt(359)
30104.687833185995
>>> distt(358)
30210.36149590906
Теги:
Хабы:
Всего голосов 13: ↑1 и ↓12 -11
Просмотры 1.6K
Комментарии Комментарии 7

Работа

Python разработчик
152 вакансии
Data Scientist
97 вакансий