Звездный год (365 дней 369 минут), Тропический год(+ 348.5 минут) и звездные сутки(1436 минут) в радиоактивном распаде

Вот здесь на GitHub.com можно скачать все данные и код

Всего данные по почти 19 дней одного года — days.dat
И 2 дня в мае 2005 года — 0505.txt
И 2 дня в мае 2006 года — 0506.txt

Средняя продолжительность тропического года по весеннему равноденствию — это на Января 1, 2000 — 365.2421897 или 365 дней, 5 часов, 48 минут, 45.19 секунд.
Звездный год — это 365 дней и 369 минут
Звездные сутки — 1436 минут — движение Земли в течении которого планета совершает полный оборот относительно сферы неподвижных звезд.

Загружаем данные и расфасовываем по рядам по 60 значений

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))];harra=[0 for t in range(int(len(matrix)/koeff))];
arra1=[[ 0 for e in range(koeff)] for t in range(int(len(matrix)/koeff))];harra1=[0 for t in range(int(len(matrix)/koeff))];
arrad=[[ 0 for e in range(koeff)] for t in range(int(len(matrixd)/koeff))];harrad=[0 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]);

Сортируем — функцией sorted по возрастанию

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;

for jk in range(int(len(matrix)/koeff)): harra[jk]=sorted(arra[jk]);harra1[jk]=sorted(arra1[jk])

z=0;
for jk in range(int(len(matrixd)/koeff)): 
 for mk in range(koeff): arrad[jk][mk]=float(matrixd[z]);z=z+1;

for jk in range(int(len(matrixd)/koeff)): harrad[jk]=sorted(arrad[jk]);

Делим на среднее значение

for jk in range(int(len(matrix)/koeff)): 
 average=mean(harra[jk])
 for mk in range(koeff): harra[jk][mk]=harra[jk][mk]/average
 average1=mean(harra1[jk])
 for mk in range(koeff): harra1[jk][mk]=harra1[jk][mk]/average1

for jk in range(int(len(matrixd)/koeff)): 
 averaged=mean(harrad[jk])
 for mk in range(koeff): harrad[jk][mk]=harrad[jk][mk]/averaged

Далее привожу функцию сравнения:

import numpy
import numpy as np
from numpy.fft import fft, fftfreq
import matplotlib.pyplot as plt
from numpy import linspace, loadtxt, ones, convolve
from numpy import array, zeros, argmin, inf
from numpy.linalg import norm

def _trackeback(D):
    i, j = array(D.shape) - 1
    p, q = [i], [j]
    while (i > 0 and j > 0):
        tb = argmin((D[i-1, j-1], D[i-1, j], D[i, j-1]))
        if (tb == 0):
            i = i - 1
            j = j - 1
        elif (tb == 1):
            i = i - 1
        elif (tb == 2):
            j = j - 1
        p.insert(0, i)
        q.insert(0, j)
    p.insert(0, 0)
    q.insert(0, 0)
    return (array(p), array(q))

def dtw(dr):
    r, c = len(dr), len(dr[0]);
    D = zeros((r + 1, c + 1));
    D[0, 1:] = inf
    D[1:, 0] = inf
    for i in range(r):
        for j in range(c):
            D[i+1, j+1] = dr[i][j];
    for i in range(r):
        for j in range(c):
            D[i+1, j+1] += min(D[i, j], D[i, j+1], D[i+1, j])
    D = D[1:, 1:];
    dist = D[-1, -1] / sum(D.shape);
    return dist, D, _trackeback(D)

def stat(x,y):
    pikk=[[float("inf") for m in range(len(y))] for t in range(len(x))];
    for e in range(len(x)):
     for l in range(len(y)): pikk[e][l]=(((abs(x[e]-y[l])))*(abs(e-l)+1));
    dt=dtw(pikk)
    p1=dt[2][0];p=dt[2][1];
    res=0;uj=0;rmu=0;
    for e in range(len(p)): 
     for l in range(len(p1)):
      pl=p[l];pe=p[e];ml=p1[l];me=p1[e];  
      #pl=l;pe=e;ml=l;me=e;  
      ret=0;re=1;rea=0;red=0;ab=abs(y[pl]-y[pe]);ac=(pe-pl);rmt=abs(1-x[ml]/y[pl]);rmr=abs(1-x[me]/y[pe]);
      au=abs(x[ml]-x[me]);ao=(me-ml);hip=math.sqrt(ab**2+ac**2);hipe=math.sqrt(ao**2+au**2);
      if (1

def dist(rast):
    statistic=0;statis=0;
    for i in range(1440):
     jk=i+000;x=harra1[jk];y=harra[jk+rast];
     rt=stat(x,y);
     statistic+=rt;
     jk=i+000;x=harra[jk];y=harra1[jk+rast];
     rt=stat(x,y);
     statis+=rt;
    return statis,statistic;


d=[];
for i in range(1440): 
 print(i);
 d.append((dist(i),i));
 print(d[-1]);

d004=[];
for i in range(len(d)):
 t4=int(d[i][0][0]+d[i][0][1])/2
 d004.append(t4)

import numpy
numpy.argsort(d004)