On essaie de trouver l'objet triple sur les mesures de gwori. FITS File: "gwori_2004Dec12.oifits" 1- Premier ajustement sur une binaire (deux objets ponctuels). code yorick utilisé: << w=lit_load_data("gwori_2004Dec12.oifits"); lit_load_model(w, LIT_PONCT_XY, 0) lit_grow_model(w, LIT_PONCT_XY, 0) lit_change_parameter(w, 6, user_value=5); r=lit_trfit_wrapper(w); lit_plot_vis2_vs_freq, w, 0, windownb=0; >> Affichage: *** Warning: TRALMOST stopped because rescaled gradient of residuals is inferior to tolerance *** Warning: TRALMOST found a final reduced Chi square of 8.87823, which means a confidence level deduced from reading a Chi2 table for 71 degrees of freedom and this final Chi square as superior value. Conclusion: On atteind un minimum local! code yorick utilisé: << lit_print_param_index_and_values(w); >> _ 1 _ 0.291745 _ 2 _ -0.000660614 mas _ 3 _ 0.00344917 mas _ 4 _ 0.691652 _ 5 _ 1.34402 mas _ 6 _ 4.32574 mas variances: [0.093141,1.56017e+09,2.61618e+09,0.017421,0.0227859,0.00370632] matrice de correlations: [[0.00867524,171534,1.46175e+07,-0.00142169,-0.000599671,-0.000302321], [171534,2.43413e+18,3.5664e+17,-62098.1,-351557,-8233.99], [1.46175e+07,3.5664e+17,6.84438e+18,-2.6288e+06,-1.51312e+06,-929148], [-0.00142169,-62098.1,-2.6288e+06,0.000303491,0.000141851,6.06165e-05], [-0.000599671,-351557,-1.51312e+06,0.000141851,0.000519199,6.27119e-06], [-0.000302321,-8233.99,-929148,6.06165e-05,6.27119e-06,1.37368e-05]] Ce resultat a été obtenu en imposant comme coordonnees (relatives) de depart pour le deuxieme objet :(0.0,5.0) en mas. Mais si on part de coordonnées beaucoup beaucoup plus grandes, comme (10.0,10.0), on atteint un secon minimum local. 2-Deuxieme ajustement sur le meme modele binaire (deux objets ponctuels). code yorick utilisé: << w2=lit_load_data("gwori_2004Dec12.oifits"); lit_load_model(w2, LIT_PONCT_XY, 0) lit_grow_model(w2, LIT_PONCT_XY, 0) lit_change_parameter(w2, 5, user_value=10.0); lit_change_parameter(w2, 6, user_value=15.0); r2=lit_trfit_wrapper(w2); >> Affichage: *** Warning: TRALMOST stopped because tolerance on parameters evolution is too small. No further improvement in the approximate solution is possible. *** Warning: TRALMOST found a final reduced Chi square of 17.3597, which means a confidence level deduced from reading a Chi2 table for 71 degrees of freedom and this final Chi square as superior value. Conclusion: On a atteint un autre minimum local! code yorick utilisé: << lit_print_param_index_and_values(w2); >> _ 1 _ 0.209497 _ 2 _ 1000 mas _ 3 _ -1000 mas _ 4 _ 0.724821 _ 5 _ 8.11835 mas _ 6 _ 16.0992 mas variances: [0.093141,1.56017e+09,2.61618e+09,0.017421,0.0227859,0.00370632] matrice de correlations: [[0.00867524,171534,1.46175e+07,-0.00142169,-0.000599671,-0.000302321], [171534,2.43413e+18,3.5664e+17,-62098.1,-351557,-8233.99], [1.46175e+07,3.5664e+17,6.84438e+18,-2.6288e+06,-1.51312e+06,-929148], [-0.00142169,-62098.1,-2.6288e+06,0.000303491,0.000141851,6.06165e-05], [-0.000599671,-351557,-1.51312e+06,0.000141851,0.000519199,6.27119e-06], [-0.000302321,-8233.99,-929148,6.06165e-05,6.27119e-06,1.37368e-05]] Ces resultats coincident avec ceux deja obtenus par Jean-Philippe Berger. Remarque: Aucune contrainte n'a été imposée sur la somme des coefficients d'intensité (énergie). Cependant, une barre d'erreur globale a été générée en tenant compte de toutes les barres d'erreur, pour pondérer la valeur de visibilité carrée égale à 1 aux frequences spatiales nulles. Le flux n'est pas égal a 1.0 a la fin de l'ajustement, mais il correspond au meilleur ajustement trouvé pour ce modele. Cela laisse supooser qu'un composante manque a ce modele, pour atteindre une énergie proche de 1. 3- L'existence de ces deux minima nous a conduit a testé un ajustement avec trois objets ponctuels. Code utilisé: << w=lit_load_data("gwori_2004Dec12.oifits"); lit_load_model(w, LIT_PONCT_XY, 0) lit_grow_model(w, LIT_PONCT_XY, 0) lit_change_parameter(w, 5, user_value=12); lit_change_parameter(w, 6, user_value=12); r=lit_trfit_wrapper(w); lit_grow_model(w, LIT_PONCT_XY, 0) lit_change_parameter(w, 8, user_value=8.0); lit_change_parameter(w, 9, user_value=15.0); r=lit_trfit_wrapper(w); >> Affichage: *** Warning: TRALMOST stopped because tolerance on parameters evolution is too small. No further improvement in the approximate solution is possible. *** Warning: TRALMOST found a final reduced Chi square of 17.3597, which means a confidence level deduced from reading a Chi2 table for 71 degrees of freedom and this final Chi square as superior value. *** Warning: TRALMOST stopped because tolerance on parameters evolution is too small. No further improvement in the approximate solution is possible. *** Warning: TRALMOST found a final reduced Chi square of 1.65142, which means a confidence level deduced from reading a Chi2 table for 68 degrees of freedom and this final Chi square as superior value. On a enfin un chi2 reduit de l'ordre de 1!!!!! Code utilisé: << lit_print_param_index_and_values(w); >> Affichage: _ 1 _ 0.138392 coefficient d'intensite du premier objet. _ 2 _ 1000 mas abscisse du premier objet _ 3 _ 1000 mas ordonnee du premier objet _ 4 _ 0.215042 coefficient d'intensite du deuxieme objet _ 5 _ 8.70994 mas abscisse (relative /objet 1) du deuxieme objet _ 6 _ 18.8621 mas ordonnee (relative /objet 1) du deuxieme objet _ 7 _ 0.652282 coefficient d'intensite du troisieme objet _ 8 _ 8.09375 mas abscisse (relative /objet 1) du troisieme objet _ 9 _ 15.4316 mas ordonnee (relative /objet 1) du troisieme objet Primaire: intensité:0.652282 Secondaire: intensité:0.215042 distance:3.4854 mas Autre: intensité:0.138392 distance: 17.4254 mas variances; [0.457772,0.735771,0.46207,1.7971,0.00180045,0.00256456,0.190153,0.00371723, 0.00165107] matrice de correlations: [[0.209556,-0.0248094,0.0325603,-0.74674,-0.00045839,0.00102551,0.0774343, -0.00075358,0.000629954], [-0.0248094,0.541358,-0.107758,0.09328,-0.000238774,-7.42451e-05,-0.00865353, -0.00064022,4.5621e-05], [0.0325603,-0.107758,0.213509,-0.182981,0.000196695,0.000208639,0.0186311, 0.000520377,2.90155e-05], [-0.74674,0.09328,-0.182981,3.22957,0.00202875,-0.00457256,-0.340887, 0.00405966,-0.00283322], [-0.00045839,-0.000238774,0.000196695,0.00202875,3.24163e-06,-3.14185e-06, -0.000217917,5.93198e-06,-2.27255e-06], [0.00102551,-7.42451e-05,0.000208639,-0.00457256,-3.14185e-06,6.57699e-06, 0.000486192,-6.46793e-06,4.14559e-06], [0.0774343,-0.00865353,0.0186311,-0.340887,-0.000217917,0.000486192,0.0361583, -0.000447999,0.00030265], [-0.00075358,-0.00064022,0.000520377,0.00405966,5.93198e-06,-6.46793e-06, -0.000447999,1.38178e-05,-4.8769e-06], [0.000629954,4.5621e-05,2.90155e-05,-0.00283322,-2.27255e-06,4.14559e-06, 0.00030265,-4.8769e-06,2.72602e-06]]