Для всех кто хочет познакомиться с книгой Космофизические факторы в случайных процессах Симона Эльевича Шноля оставляю ссылку
Далее мой алгоритм по реальным данным 48 дней с 19 декабря 2012 года по 6 февраля 2013 года распад стронция-90, измерение дозиметром ДП-5В.
>>> matrix[:60]
[177, 211, 215, 201, 196, 190, 204, 212, 202, 198, 222, 212, 211, 223,
214, 193, 228, 187, 224, 240, 205, 235, 227, 217, 203, 215, 184, 207,
200, 225, 199, 213, 229, 225, 222, 216, 209, 211, 209, 233, 215, 203,
216, 182, 232, 207, 210, 230, 210, 212, 230, 223, 214, 227, 232, 189,
212, 188, 204, 192]Каждая точка данных сколько актов распада было за 1 секунду.
Далее функция сравнения:
def dt2(x,y):
uu=0;
for i in range(len(x)-3):
uu1=0;uu2=0;uu3=0;
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;Это подгрузка данных и почасовая расфасовка:
import pickle
matrix = pickle.load(open("natrix.pkl","rb"));
z=0;koeff=60*60
arra1=[[ 0 for e in range(koeff)] for t in range(int(len(matrix)/koeff))];
for jk in range(int(len(matrix)/koeff)): #int(len(matrix)/koeff)
for mk in range(koeff): arra1[jk][mk]=float(matrix[z]);z=z+1;
Функция попарного сравнения часовых отрезков со сдвигом 'rast' - расстояние.
def distt(rast):
statistic=0;
for jk in range(100):
x=arra1[jk];y=arra1[jk+rast];
statistic+=dt2(x,y);
return statistic;
Далее расчет сравнения 100 раз-часов с отступом почасовых отрезков времени от +300 до +800 часов. Видна Эвекция 111171 (число-дистанция) - +762 часа, Вариация 111321 (число-дистанция) - +328 часа и звездный лунный месяц - + 27.32*24 +656 часов - 111306 (число-дистанция), период между двумя полнолуниями +714 часов 111348 (число-дистанция)
>>> for i in range(300,800):
... print(distt(i),i);
...
111467.51279175737 300
111456.94884998219 301
111560.49898737374 302
111528.18026695744 303
111441.88393507873 304
111592.80776135616 305
111470.9559966939 306
111507.78022934493 307
111637.34079748808 308
111676.25554055597 309
111253.35606458329 310
111624.96022984687 311
111407.84051456978 312
111523.5042529058 313
111389.70065279244 314
111545.27813327349 315
111393.90664361506 316
111652.49042663965 317
111521.86525819649 318
111577.26447673695 319
111515.39756660172 320
111443.3539917646 321
111567.73146520492 322
111376.7129928991 323
111384.5491740358 324
111487.76534411187 325
111513.74776289536 326
111383.04558713648 327
111321.24144567881 328 - Вариация
111747.99312081255 329
111443.02908443383 330
111630.27087613485 331
111747.76721140447 332
111408.39048315817 333
111473.53118176629 334
111451.89327382507 335
111607.40643875868 336
111485.03060886491 337
111501.23435137501 338
111601.95647004734 339
111410.88705609205 340
111650.70924791638 341
111730.29946813371 342
111643.68521665082 343
111580.84176282027 344
111523.86649141887 345
111708.74340456662 346
111546.83299160897 347
111564.68731087203 348
111766.53086589994 349
111568.0886701693 350
111615.46028693976 351
111551.62260653163 352
111484.62983504617 353
111492.76513310458 354
111459.24534684437 355
111548.17695200817 356
111431.8470765421 357
111499.35423415787 358
111570.5551273607 359
111522.82945303114 360
111442.4546211406 361
111474.23166759517 362
111456.69182591826 363
111509.06729575292 364
111533.54900032145 365
111488.16179371731 366
111506.79348614311 367
111392.38802219921 368
111314.84447064178 369
111405.82311124234 370
111555.6924893529 371
111394.66050244718 372
111391.28134166659 373
111536.33411552628 374
111142.5935395788 375
111409.94785461755 376
111422.4326514421 377
111300.00290925261 378
111243.36558475497 379
111354.52943074693 380
111394.58194618016 381
111401.69942294215 382
111445.82953222671 383
111462.88163528839 384
111434.54761647561 385
111289.4079435053 386
111421.68333790186 387
111119.08006356718 388
111584.7780903453 389
111657.57383452768 390
111609.28292854635 391
111366.08273892236 392
111309.67490526177 393
111478.10431404825 394
111532.30967275598 395
111501.37011207212 396
111470.90452555324 397
111503.16372508004 398
111435.35205085679 399
111622.98194374684 400
111564.61467133387 401
111501.83122777389 402
111590.44319354721 403
111639.44527028559 404
111402.99066323852 405
111543.13887331121 406
111384.27870217085 407
111583.51827420754 408
111541.21218902063 409
111526.73584354557 410
111428.92484516394 411
111382.81138668982 412
111459.81440797258 413
111375.50565111767 414
111414.7834907925 415
111333.54693831854 416
111230.87584491842 417
111449.57300995327 418
111616.7090290731 419
111421.46436695165 420
111592.51585545261 421
111509.57049229131 422
111483.37508128239 423
111459.94887361131 424
111507.68710382003 425
111469.26912182695 426
111564.02079979074 427
111460.79267005807 428
111565.74105313339 429
111499.73061479755 430
111572.61502813357 431
111435.31877209507 432
111450.8521073132 433
111418.22985358536 434
111641.24911547151 435
111436.60280910032 436
111363.28570314485 437
111409.22746443994 438
111442.31610779677 439
111124.48511434972 440
111413.10504161498 441
111459.85800722202 442
111474.77979627208 443
111476.62462880297 444
111260.57725222025 445
111350.42947877642 446
111384.07572112922 447
111453.01114534325 448
111245.37322874447 449
111436.01808735648 450
111202.69701845391 451
111277.51973466863 452
111374.32633939214 453
111200.57186932137 454
111386.89937843832 455
111388.34628396771 456
111581.64177064892 457
111351.49655337271 458
111420.9103491089 459
111496.68861312151 460
111559.0797683151 461
111256.37298408982 462
111363.65954451491 463
111318.0567273391 464
111308.22912921835 465
111521.4242317889 466
111375.52582695743 467
111324.11385585566 468
111518.97290071417 469
111272.0094827227 470
111270.61985193074 471
111377.22872581487 472
111390.79279516307 473
111206.98537066634 474
111259.18825290748 475
111171.20070675813 476
111455.60746726347 477
111490.25283821547 478
111470.99392572229 479
111542.01546152157 480
111302.52298411136 481
111498.28985056473 482
111427.62018429863 483
111482.05544116956 484
111245.34184517301 485
111415.9453824231 486
111515.92485115012 487
111394.54296691522 488
111307.8271024343 489
111546.41960655594 490
111387.23469231493 491
111559.64612628415 492
111253.66198135744 493
111414.91552146313 494
111487.79155993623 495
111415.01322135837 496
111413.0182757029 497
111628.76294166897 498
111475.89810856844 499
111575.71863230762 500
111499.90799916143 501
111409.13274648543 502
111310.37743330943 503
111529.28931245313 504
111458.367796574 505
111622.51015956169 506
111635.18750777635 507
111415.6642962675 508
111294.18976000235 509
111409.20867293808 510
111380.17611770755 511
111454.02674934926 512
111682.289035573 513
111495.10554424785 514
111581.98198680482 515
111439.74475835847 516
111582.99950256338 517
111443.70773136121 518
111416.32589558353 519
111583.80036732467 520
111649.26897995496 521
111726.6230230274 522
111491.76749192119 523
111481.89448510073 524
111412.0401168958 525
111512.58563533353 526
111461.29192633531 527
111614.97639075777 528
111441.02000879629 529
111558.97142197369 530
111580.27989782614 531
111625.34453434803 532
111413.50677126601 533
111521.93772730826 534
111769.09673339466 535
111339.48592554551 536
111637.7548091163 537
111506.8043850952 538
111329.39610563427 539
111506.8260320286 540
111498.51179849287 541
111490.27873692919 542
111482.18626765408 543
111595.98388129947 544
111678.65505128623 545
111551.48823362622 546
111656.90005989568 547
111636.22770871657 548
111674.49942720564 549
111484.32789977764 550
111513.34539934693 551
111444.77020043157 552
111618.60570395408 553
111605.54047067594 554
111594.00595341698 555
111497.3872519592 556
111617.23837290473 557
111559.89735976834 558
111607.69628597856 559
111419.41524183573 560
111670.84787266798 561
111641.05981174755 562
111674.00560730195 563
111647.74241840806 564
111633.63051011974 565
111721.25132577177 566
111345.38016119805 567
111724.41486866983 568
111364.4114299918 569
111653.65878095517 570
111465.14912991022 571
111606.46222957493 572
111598.2059385263 573
111583.3643582366 574
111412.975164046 575
111655.39448785366 576
111392.97974349317 577
111500.89759992674 578
111447.78323775239 579
111733.18088175352 580
111406.03835255414 581
111412.95936489775 582
111508.81074111394 583
111663.66854229083 584
111531.37366906407 585
111699.26065105238 586
111632.15328001621 587
111238.91451884362 588
111628.85337583453 589
111629.362809086 590
111595.0096770417 591
111627.63709723114 592
111592.9692767375 593
111329.86772352368 594
111715.19762988444 595
111524.17426665635 596
111738.69738823752 597
111629.86027291436 598
111414.89923723976 599
111498.49645815085 600
111552.45954857652 601
111530.73178188212 602
111590.9378403336 603
111621.95788869068 604
111568.33956571495 605
111565.34931557626 606
111522.95965409772 607
111529.15659256016 608
111793.36764499605 609
111363.30240604545 610
111585.45770287639 611
111308.40607653258 612
111533.97126231804 613
111479.92454616942 614
111293.94631660989 615
111462.24047276816 616
111638.02088438811 617
111305.0813248197 618
111542.76684692693 619
111472.28493376159 620
111381.5632300407 621
111484.41659856569 622
111467.88940238862 623
111641.85737803677 624
111255.6176620698 625
111581.06776223767 626
111492.93405786266 627
111475.51444910164 628
111454.60971598487 629
111652.85429735159 630
111376.1010876821 631
111490.5359707696 632
111505.85695737772 633
111479.70486302527 634
111469.74292701905 635
111496.61588053148 636
111611.11775961649 637
111564.619223655 638
111506.59299069928 639
111368.62470428496 640
111491.23744678978 641
111414.64740661591 642
111529.0400169789 643
111414.61369951541 644
111407.60137571274 645
111310.20736329927 646
111674.64796265004 647
111433.41260867512 648
111428.47039360942 649
111330.28160322034 650
111641.85326982709 651
111286.53170293715 652
111648.57459561471 653
111404.14272340055 654
111375.89324315138 655
111306.287688769 656 - звездный лунный месяц 27.32
111458.86211424357 657
111306.45987539327 658
111437.17622510383 659
111514.15983308444 660
111590.39531771312 661
111499.98923691109 662
111455.73791422772 663
111468.04994900993 664
111517.16793340974 665
111352.95859353546 666
111504.47072828247 667
111223.99164658932 668
111617.75293480075 669
111467.6512476625 670
111505.7992966194 671
111500.63587164156 672
111394.3971065766 673
111372.848475461 674
111422.98829025679 675
111311.45974264054 676
111290.66853472979 677
111338.90673886561 678
111464.30539328471 679
111516.4527185785 680
111343.95635821263 681
111524.8943704452 682
111378.64153399499 683
111446.62774924764 684
111324.34310023636 685
111510.40916661297 686
111252.68491274222 687
111286.60616123048 688
111429.7893178485 689
111424.56970844763 690
111589.81952522197 691
111513.13661581733 692
111490.32714369766 693
111461.25582563298 694
111469.93206554669 695
111439.0897739749 696
111511.77236829432 697
111277.02749721157 698
111534.25979787145 699
111430.03380832344 700
111506.3417029236 701
111393.76605525531 702
111463.86617904271 703
111289.68135536784 704 - средний период между фазами 29.5306
111502.08413640247 705
111543.62598318302 706
111511.0055832941 707
111484.87453800725 708
111327.37846270134 709
111357.14884623184 710
111397.85510440728 711
111518.1853387865 712
111520.00937452086 713
111348.37899050966 714 - период между двумя полнолуниями
111447.04287010482 715
111489.88182361223 716
111438.51221367109 717
111565.90815651458 718
111318.32765322436 719
111453.39901758776 720
111345.81331769945 721
111281.08128626294 722
111357.78956737187 723
111362.8786641281 724
111296.92198031071 725
111226.28372062222 726
111265.01760740408 727
111368.27123096139 728
111394.95994345455 729
111303.20850572128 730
111152.1162904444 731
111292.2090495194 732
111560.61404495954 733
111311.566794089 734
111252.07227318751 735
111359.36852949043 736
111429.63144408909 737
111299.74063544259 738
111357.88922415949 739
111121.63278620502 740
111427.17360124063 741
111420.71157250201 742
111252.54480590341 743
111361.2190104509 744
111329.42545147346 745
111228.60952739263 746
111355.4044499485 747
111407.19683977273 748
111431.19111508688 749
111449.45277587471 750
111487.48344317154 751
111362.74219350381 752
111468.79417999917 753
111509.76063304825 754
111217.8628127819 755
111354.8765878456 756
111417.0579459793 757
111116.8595867451 758
111210.24350904173 759
111241.12933352943 760
111279.21817827334 761
111171.34866525831 762 - Эвекция
111440.34946114497 763
111477.87208356385 764
111398.30951888424 765
111239.28358866072 766
111251.5774729378 767
111370.60752179877 768
111307.28011819102 769
111572.06651634726 770
111531.10499431276 771
111359.13079492023 772
111463.05614494417 773
111424.57994703732 774
111342.21533951438 775
111414.58778554782 776
111344.42981543491 777
111199.31007224412 778
111097.63107741298 779
111226.6913877649 780
111293.62656806431 781
111338.50584269498 782
111200.8392985953 783
111283.33993454164 784
111293.5019388204 785
111394.54518746195 786
111214.69591897709 787
111285.47353167742 788
111615.53936748505 789
111601.65357491955 790
111260.74347223944 791
111217.75403023534 792
111304.05636036192 793
111478.77131610972 794
111422.84292047277 795
111253.22595135917 796
111371.12748646589 797
111345.848028843 798
111322.62873004763 799В период между 19 декабря 2012 по 6 февраля 2013 года - полнолуние
Пт | 28 декабря 2012 г. | 13:21:12 |
Вс | 27 января 2013 г. | 07:38:35 |
Новолуние
Пт | 11 января 2013 г. | 22:43:36 |
Это были лунные периоды, напоминаю что самый сильный период - звездный год или 365 дней 6 часов и 9 минут. Из-за звездного года у нас високосный каждый 4-ый год. Каждый 100 год не-високосный, а каждый 400-ый опять с 29 февраля! Это Григорианский календарь
