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')
