Estimados amigos del foro
No se apenas como funciona Matlab con lo que espero saber explicar lo que necesito para ver si me podeis ayudar.
Hace varias semenas decidi ponerme manos a la obra para poder estudiar la posibilidad de estudiar la existencia de algun patron en los datos que proporcionan los radiotelescopios de SETI. Consegui llegar a esa fuente de datos aunque el inconveniente que el software de tratamiento opera bajo Linux. La otra alternativa es abrir los archivos en formato “dat” con Matlab. Debido el inconveniente de que el foro esta en ingles no puedo entender mucho, pero para abrir estos archivos DAT comentan que debo de introducir un “codigo octava” para traducir los datos . Para dar mas pistas les dejo este supuesto codigo para traducir el archivo, aunque no se como proceder con ello.
Here is an example of reading the first 8 complex samples of an observation file in both Octave and R (for comparison) preceded by a dump of the 8-bit signed coefficient pairs:
rackrman@dev:~/temp/R$ xxd 2010-01-22-kepler-exo4-1420mhz-1-of-3.dat | head -n 1
0000000: 02f4 01ff 0cd6 f620 ecf5 17d2 010f f9fa ....... ........
rackrman@dev:~/temp/R$ octave -q
octave:1> f=fopen("2010-01-22-kepler-exo4-1420mhz-1-of-3.dat");
octave:2> N=8;
octave:3> coeffs=fread(f,N*2,"int8");
octave:4> cmplx(1:N)=coeffs(1:2:N*2)+coeffs(2:2:N*2)*i;
octave:5> disp(cmplx)
2 - 12i 1 - 1i 12 - 42i -10 + 32i -20 - 11i 23 - 46i 1 + 15i -7 - 6i
rackrman@dev:~/temp/R$ R
R version 2.8.1 (2008-12-22)
> f <- file("2010-01-22-kepler-exo4-1420mhz-1-of-3.dat","rb")
> N <- 8
> x <- readBin(f, integer(), n=N*2, size=1, signed=TRUE, endian="little")
> complex(N, x[seq(1,N*2,2)], x[seq(2,N*2,2)])
[1] 2-12i 1- 1i 12-42i -10+32i -20-11i 23-46i 1+15i -7- 6i
