O-C diagram för Miror

Inom Variabelsektionen studeras olika typer av föränderliga stjärnor främst med avseende på ljusstyrka och periodicitet. Sektionsledare är Chris Allen som kan ge råd och tips. Diskussioner om variabler
ThomasK
Inlägg: 1259
Blev medlem: 2010-06-29 12:10:36
spamtest: JA

Re: O-C diagram för Miror

Inläggav ThomasK » 2012-03-22 17:40:19

Arbetet fortgår med att lägga in historiska maxima för de 482 mirorna som det nu finns O-C diagram för, drygt 1500 historiska maxima från "Annals of Harvard College Observatory, vol. 55, pp.95-291" finns nu med i diagrammen.

Samtidigt har jag börjat fundera på hur man bästa kan analysera materialet. Några saker jag vill försöka hitta är dels stjärnor med konstant pågående periodförändring (parapoliskt O-C digram), t.ex R Aql, och dels stjärnor där perioden förändras cykliskt (sågtandsformat O-C diagram), t.ex X Aur och ST Cyg eller vid ett enstaka tillfälle. Problemet är att skilja ut dessa underliggande förändringar från det slumpnässiga brus mirornas period har. I O-C diagrammen adderas allt brus allteftersom och man får "random-walk" mönster som kan vara svåra att skilja från riktiga förändringar.

För detta har jag gjort ett simuleringsprogram som simulerar 1000 cykler för en mira med olika värden för 5 parametrar:
1) Konstant periodförändring
2) Slumpmässig normalfördelad ändring av perioden.
3) Slumpmässig normalförelat fel i bestämning av stjärnans maximum.
4) Period för periodiska periodförändringar mellan 2 värden.
5) Amplitud för dessa periodiska förändringar.
Ja har sen kört 10000 simuleringar för olika uppsättningar av dessa 5 parametrar. I varje simulering plockar jag ut ett slumpmässigt utsnitt på 50, 100 eller 200 cykler (som motsvarar den period vi har data för de riktiga stjärnorna) och beräknar O-C diagram och gör Eddington/Plakidis-testet. På så sätt får jag fram vad ett genomsnittligt diagram har för egenskaper med olika uppsättningar av parametrarna.

Tanken är sen att försöka matcha de riktiga O-C diagrammen mot de simulerade och se hur stor variationen är för de 10000 körningarna för diagrammet som matchar bäst. T.ex om ett rikigt O-C diagram med en periodförändring om +0.25 dagar/cykel matchar ett simulerat med förändringen -1 till +1 dag/cykel har man inget signifikant resultat, men om det istället matchar ett simulerat diagram med förändringen +0.2498 till +0.2502 dagar/cykel bör man rätt säkert kunna säga att här har vi en stjärna med underliggande periodförändring.

Tankar om detta och andra idér om hur man kan göra analysen tas tacksamt emot.
Thomas Karlsson
ThomasK
Inlägg: 1259
Blev medlem: 2010-06-29 12:10:36
spamtest: JA

Re: O-C diagram för Miror

Inläggav ThomasK » 2012-07-01 15:19:23

Det stora arbetet med att samla in historiska maxima för Mirorna är nu klart. Databasen innehåller nu drygt 42000 maxima för 488 Miror, varav 33000 är beräknade från AAVSO/AFOEV/VSOLJ-data och 9000 manuellt avskrivna från gamla Harvard-artiklar.

O-C diagrammen är nu som standard beräknade med stjärnans medelperiod under diagrammets tidsperiod istället för perioden från GCVS-katalogen. Under varje diagram finns nu också en länk (LC) till varje stjärnas medelljuskurva.
Thomas Karlsson
Hans Bengtsson
Inlägg: 3900
Blev medlem: 2011-03-04 18:19:54
spamtest: JA

Re: O-C diagram för Miror

Inläggav Hans Bengtsson » 2012-07-01 15:32:27

Hoppas alla förstår vilket enastående jobb du har gjort! Det finns anledning att tro att detta blir den internationella standardsidan när det handlar om periodvariationer hos miror. Kan inte på rak arm se något tidigare arbete av någon svensk amatörastronom som varit av högre kvalitet.

Edit:
... och standardsida inte bara när det gäller periodvariationer utan också medelljuskurvor. Det klassiska bokverket "Studies of Long Period Variables" (Campbell) har äntligen fått en modern och uppdaterad efterföljare.

M v h
Hans Bengtsson
wilt
Inlägg: 1399
Blev medlem: 2011-10-24 08:34:32
spamtest: JA
Ort: Lund
Kontakt:

Re: O-C diagram för Miror

Inläggav wilt » 2012-07-01 16:29:19

Mycket imponerande! Förtjänar verkligen internationell uppmärksamhet!

/Gustav
8 x 21 monokular / 6x30 Carl Zeiss Silvamar / Celestron SkyMaster 15x70 / 9 cm f/14 Maksutov Cassegrain / 20 cm f/6 Newton.

hemsida / blogg / astronomiblogg / SVO: HGG / AAVSO: HGUA
Hans Bengtsson
Inlägg: 3900
Blev medlem: 2011-03-04 18:19:54
spamtest: JA

Re: O-C diagram för Miror

Inläggav Hans Bengtsson » 2012-07-01 20:06:21

Flera av världens ledande variabelspecialister har redan höjt arbetet till skyarna. Några exempel:
"This is a terrific result."
"This is a fantastic job."
"Thomas, Great work!!!"
"Very interesting and gratifying"

M v h
Hans Bengtsson
Användarvisningsbild
Robin Andersson
Inlägg: 761
Blev medlem: 2011-03-04 20:42:46
spamtest: JA
Ort: Göteborg
Kontakt:

Re: O-C diagram för Miror

Inläggav Robin Andersson » 2012-07-01 20:18:42

Grymt jobbat Thomas! Kul att se också att det uppskattas så enormt av de som ännu har sett det (i alla fall vad man kan säga från deras kommentarer) - finns ingen anledning till varför det inte bör uppskattas.
Tror samtidigt Thomas att du har kommit med ett stort + för hela SVO med att du publicerade detta på AAVSO.

Keep up the good work!

mvh
Robin
if there's no life elsewhere in the universe, it's a terrible waste of space
Hans Bengtsson
Inlägg: 3900
Blev medlem: 2011-03-04 18:19:54
spamtest: JA

Re: O-C diagram för Miror

Inläggav Hans Bengtsson » 2012-07-01 20:40:55

Man kan följa de internationella reaktionerna här:

http://www.aavso.org/o-c-diagrams-miras

M v h
Hans Bengtsson
ThomasK
Inlägg: 1259
Blev medlem: 2010-06-29 12:10:36
spamtest: JA

Re: O-C diagram för Miror

Inläggav ThomasK » 2012-07-02 13:12:58

Tack för de vänliga orden! Det var värst vilken uppståndelse det blev på AAVSO nu, jag har inte fått några reaktioner där på de tidigare inläggen om O-C diagrammen. Antagligen beror det på att AAVSO gjort om så man nu ser de senaste foruminläggen på deras förstasida och därför exponeras mer.
Thomas Karlsson
Hans Bengtsson
Inlägg: 3900
Blev medlem: 2011-03-04 18:19:54
spamtest: JA

Re: O-C diagram för Miror

Inläggav Hans Bengtsson » 2012-07-02 17:31:45

Thomas:

En detalj när det gäller medelljuskurvan för U CrA:
Detta är den enda stjärna i det omfångsrika materialet där nollpunkten i x-axeln inte sammanfaller med den heldragna kurvans max. Vet inte riktigt vad det beror på.

mvh
Hans Bengtsson
ThomasK
Inlägg: 1259
Blev medlem: 2010-06-29 12:10:36
spamtest: JA

Re: O-C diagram för Miror

Inläggav ThomasK » 2012-07-02 20:14:14

Hans Bengtsson skrev:Thomas:

En detalj när det gäller medelljuskurvan för U CrA:
Detta är den enda stjärna i det omfångsrika materialet där nollpunkten i x-axeln inte sammanfaller med den heldragna kurvans max. Vet inte riktigt vad det beror på.

mvh
Hans Bengtsson


Ja, jag har sett detta också och hoppades det skulle passera obemärkt :green: U CrA var en av de sista jag tog med, den har egentligen för litet underlag av observationer för att medelljuskurvan ska bli bra. Men då jag inte hade med någon stjärna alls från Södra Kronan och det fanns rätt många historiska maxima för den fick den vara med ändå. Att maxpunkten hamnar snett beror nog på att noll-punkten först bestäms med ett 8-gradspolynom men medelljuskurvan sedan representeras av ett 12-gradspolynom och att maxpunkten mellan dessa kurvor hamnar olika.
Thomas Karlsson
Användarvisningsbild
Robin Andersson
Inlägg: 761
Blev medlem: 2011-03-04 20:42:46
spamtest: JA
Ort: Göteborg
Kontakt:

Re: O-C diagram för Miror

Inläggav Robin Andersson » 2012-07-02 20:54:13

Har du funktionerna för dina två polynom Thomas? Kanske kan komma på varför genom att tänka rent matematiskt. Jobbade väldigt mycket med kurvkonstruktion första året på Chalmers, var dock givetvis bara mesta grunder men det räcker långt.

Mvh
Robin
if there's no life elsewhere in the universe, it's a terrible waste of space
ThomasK
Inlägg: 1259
Blev medlem: 2010-06-29 12:10:36
spamtest: JA

Re: O-C diagram för Miror

Inläggav ThomasK » 2012-07-02 23:08:34

Robin Andersson skrev:Har du funktionerna för dina två polynom Thomas? Kanske kan komma på varför genom att tänka rent matematiskt. Jobbade väldigt mycket med kurvkonstruktion första året på Chalmers, var dock givetvis bara mesta grunder men det räcker långt.

Mvh
Robin


Hej.
De går att ta fram, men jag är bortrest några dagar framöver så jag kan inte titta på det just nu, men jag kan skicka uppgifter till dig senare om du är intresserad.
Thomas Karlsson
Användarvisningsbild
Robin Andersson
Inlägg: 761
Blev medlem: 2011-03-04 20:42:46
spamtest: JA
Ort: Göteborg
Kontakt:

Re: O-C diagram för Miror

Inläggav Robin Andersson » 2012-07-03 09:14:32

ThomasK skrev:Hej.
De går att ta fram, men jag är bortrest några dagar framöver så jag kan inte titta på det just nu, men jag kan skicka uppgifter till dig senare om du är intresserad.


Ja gärna! Jag lovar inte dock att jag kan få ut något men ska försöka. Vart finns medelvärdeskurvan Hans nämnde? Hittar bara en kurva och det är den du har lagt upp bland alla andra.

Mvh
Robin
if there's no life elsewhere in the universe, it's a terrible waste of space
Hans Bengtsson
Inlägg: 3900
Blev medlem: 2011-03-04 18:19:54
spamtest: JA

Re: O-C diagram för Miror

Inläggav Hans Bengtsson » 2012-07-03 09:50:00

Via SVO:s sida "Ljuskurvor och resultat" kan du gå till sidan med O-C-diagram. Där finns, nästan längst upp i den inledande texten, en klickbar länk till medelljuskurvorna. Alternativt kan du, vid O-C-diagrammet för U CrA, klicka på LC så får du medelljuskurvan för just den stjärnan.

Jag har den bestämda uppfattningen att länkarna till sidorna med medelljuskurvor och O-C-diagram tydligt bör lyftas fram i SVO och kanske även Astronet så att de blir enklare att hitta. Dessa sidor representerar något av det bästa som kommit fram ur svensk amatörastronomi.

M v h
Hans Bengtsson
ThomasK
Inlägg: 1259
Blev medlem: 2010-06-29 12:10:36
spamtest: JA

Re: O-C diagram för Miror

Inläggav ThomasK » 2012-07-12 21:14:57

Jag har lagt till en ny funktion på sidan över O-C diagram. Genom att klicka på länken Avg under varje diagram får man ett diagram där O-C värdena är beräknade som ett glidande medelvärde av 3 punkter före och 3 punkter efter varje punkt. På så sätt får man ett diagram där detaljerna framgår mycket bättre.

De glidande medelvärdena har jag också använt för en helt ny sida, http://var.astronet.se/mirainfoper.php, där man ser stjärnornas period och hur den varierat över tiden istället för O-C. Perioden har jag fått fram genom att beräkna lutningen av O-C kurvan i varje punkt i de glidande medelvärdes O-C diagrammen. Resultatet blev riktigt bra tycker jag.
Thomas Karlsson

Återgå till "SAAF - Variabelsektionen"

Vilka är online

Användare som besöker denna kategori: 1 och 0 gäst