Sure, https://github.com/SubTachyon/Astronomy-3-Lab/blob/master/ellipsifon_v2.
Right now the code generates its own data, to make things easier. But if you want to plot to actual data comment out line 61 and uncomment lines 58 & 59. I uploaded two sets of data...
I made a small modification to the code so now instead of using deviation of point p_i from p^,_i I use deviation of p_i from p_c where p_c is the closest point to p_i. Is that what you meant?
I've noticed no improved accuracy of this fit after this change, in fact I didn't observe any...
First of all thank you for taking your time and looking at the code. :)
The function I am minimizing is deviation - line 63. The relevant part probably being:
...
for j in range(len(data[:,0])):
deviations[j] = np.sqrt((data[j,0]-dummy_ellipse[j,0])**2 +...
https://github.com/SubTachyon/Astronomy-3-Lab/blob/master/ellipsifonis my code at the moment and it produces something like this. Red ellipse is the fit, blue ellipse is the data.
The minimization seems to do a half-decent job at getting the orbital parameters correct but it fails to move the...
So I looked at the MCMC algorithm and I think I am a little bit over my head. :P I am not very well educated in statistics.
I went back and gave scipy.optimize.minimize another shot, but it's not performing very well with 5 variables. In addition to that my dummy ellipse is not centred with the...
I did look at scipy.optimize but primarily at the leastsq function which was probably the wrong place to focus on. I'll look around some more now that I know it's not the only option are out there, starting with the MCMC. Thank you. :)
Hello again! So I made a function that generates N number of points on an ellipse that is defined by the input variables a, e, i, ω, and Ω. It works fine.
The problem I have now is having to do squared minimization of a function with 5 variables. How do I do that? Pretty much all the examples I...
Hey, sorry if this is not the right section to post in, the topic is a bit ambiguous.
I've generated a set of coordinate points for the orbit of Mercury (3d cartesian) and now I want to fit an ellipse through it so I can get an accurate estimate of the location of the perihelion. I am using...
Thank you, that was a useful link.
''Those who harbor no qualms about pseudo-tensors will say that radiant energy becomes gravitational energy.''
That's what I thought about as well, though I am not entirely satisfied with it as an answer and I guess that until I cover GR in more depth I...
I don't believe you read my question properly. It has nothing to do with velocity; I'm concerned with the energy of the waves.
If a photon is created with the frequency of 100 Hz and after traveling between two distant galaxies due to the expansion of the universe it's frequency drops to 10...
I am assuming the answer to my question is no, but what am I missing?
My reasoning is very basic: E=hf, therefore as the universe expands the wavelength of all the far traveling radiation is increased and due to constant velocity c their frequency must decrease which translates into them losing...
So one guy comes to me and says that he doesn't believe that a supernova 50ly away would be dangerous enough to pose a lethal threat to us. I contradicted him because based on my general reading supernovas (depending on their magnitude) between at least a 100 to 3000 light years away would...
I was pondering our reliance on and the limitations of todays electrical batteries. It seems to me (uneducated member of the public) that while technology continues to develop with an impressive momentum our energy storage methods struggle to fuel these ever evolving technologies effectively...