Traitement des spectres radio "Würzburg" avec Scilab
Publié : 31 janv. 2013 10:48
Un truc un peu technique ...
Comme discuté lors de la réunion d'hier soir, voici un petit bout de code pour le logiciel (gratuit) Scilab, permettant de traiter les spectres du radiotélescope Würzburg.
Ce code est "brut de fonderie" et améliorable. L'enregistrer dans le dossier où sont les spectres, dans un ficheir texte style "WurzburgPlot.sce". Il faut changer dans le code le nom de base des fichiers (fn = ...). Il s'agit de faire la moyenne d'un ensemble de spectres. Pour le moment la partie graphique est minimaliste, un simple plot. Mais vous pouvez l'améliorer (je verrai ça quand j'aurais le temps), en consultant la doc de Scilab (http://www.scilab.org" onclick="window.open(this.href);return false;).
Comme discuté lors de la réunion d'hier soir, voici un petit bout de code pour le logiciel (gratuit) Scilab, permettant de traiter les spectres du radiotélescope Würzburg.
Code : Tout sélectionner
// Objet
obj = 'M33';
// Nom de base des fichiers
fn = 'data_M33_20121029_17h12mn26s_';
// Nombre de fichiers
num = 57;
// Nombre de points
sz = 2048;
//
y = zeros(sz,1);
//
for i = 1:num
fni = fn+string(i)+'.class';
dni = fn+string(i);
//table = fscanfMat(fni,"%lg');
fid=mopen(fni,'r');
table = mfscanf(sz,fid,'%f %f %f');
y = y+table(:,1);
x = table(:,2);
v = table(:,3);
mclose(fid);
end
//
y = y/sz;
//
subplot(2,1,1)
plot(v,y);
xlabel('km/s')
subplot(2,1,2)
plot(x,y);
xlabel('MHz')
