
Для всех кто хочет познакомиться с книгой Космофизические факторы в случайных процессах Симона Эльевича Шноля оставляю ссылку



Далее мой алгоритм по реальным данным 48 дней с 19 декабря 2012 года по 6 февраля 2013 года распад стронция-90, измерение дозиметром ДП-5В.
start='2012/12/19/23/10/35' stop=2013.02.06.22.13.07
>>> 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 февраля! Это Григорианский календарь