Foro - Perl en Español

por **lis** » 2013-03-13 16:30 @729

Hola, amigos del foro. Tanto tiempo, espero estén todos bien. Les cuento que tengo una matriz de 365x24 (días por hora) y quiero obtener dos promedios por mes, y como saben cada mes tiene diferente cantidad de días (la matriz está en días julianos), lo que he hecho me resulta a medias ya que solo guarda el último resultado, les muestro el script:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl -w
 
my %hash
    = qw (1 15 16 31 32 45 46 59 60 74 75 90 91 105 106 120 121 135 136 151 152 166 167 181 182 196 197 212 213 227 228 243 244 258 259 273 274 288 289 304 305 319 320 334 335 349 350 365);
 
for my $jul ( sort { $hash{$a} <=> $hash{$b} } keys %hash ) {
 
    open( IN, ">graficos.m" );
    print IN "clear all; clc\n";
    print IN "x= load('prueba.txt');\n";
    print IN "a=x(:,1);   m=x(:,2);  d=x(:,3);  h=x(:,4);  vv=x(:,5);  dv=x(:,6);\n";
    print IN "V=reshape(x(:,5)',365,24); D=reshape(x(:,6)',365,24);  \n";
    print IN "mv=mean(V);  md=mean(D);\n";
    print IN "G2R = (pi / 180);\n";
    print IN "gar = G2R * D;\n";
    print IN "u = (-1) .* V .* sin(gar); v = (-1) .*V .* cos(gar);\n";
    print IN "prom=[];\n";
    print IN "for K = 1:365 \n";
    print IN "Samp{K} = mean(V($jul:$hash{$jul},:));\n"; # realiza un promedio, por ejemplo, desde la fila 1 a la 15 para todas las columnas
    print IN "end;\n";
    print IN "prom=[prom; Samp(K)];\n";
    print IN "save Matriz.txt prom;\n";
    system("/usr/bin/octave -q graficos.m");
}
 
Coloreado en 0.005 segundos,  usando GeSHi 1.0.8.4

Si resultara bien debería generar una matriz de 24x24, pero solo genera una matriz de 1x24.

Por favor les pido que me ayuden. Muchas gracias por su tiempo.

my @array
    = qw (1 15 16 31 32 45 46 59 60 74 75 90 91 105 106 120 121 135 136 151 152 166 167 181 182 196
        197 212 213 227 228 243 244 258 259 273 274 288 289 304 305 319 320 334 335 349 350 365);
 
for (my $i = 0; exists($array[$i]) and $array[$i] < 365; $i += 2) {
    say $array[$i], " -> ", $array[$i+1];
 
    open( IN, ">graficos.m" );
    print IN <<EOF;
clear all; clc
x= load('prueba.txt');
a=x(:,1);   m=x(:,2);  d=x(:,3);  h=x(:,4);  vv=x(:,5);  dv=x(:,6);
V=reshape(x(:,5)',365,24); D=reshape(x(:,6)',365,24);
mv=mean(V);  md=mean(D);
G2R = (pi / 180);
gar = G2R * D;
u = (-1) .* V .* sin(gar); v = (-1) .*V .* cos(gar);
prom=[];
for K = 1:365
Samp{K} = mean(V($array[$i]:$array[$i+1],:));
end;
prom=[prom; Samp(K)];
save Matriz.txt prom;
EOF
    close IN;
 
    system("/usr/bin/octave -q graficos.m");
}
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

Lo que veo en el código es que se genera un archivo graficos.m por cada par clave/valor, y ejecuta octave para esa pareja. Pero entonces, la instrucción 'save' sobre escribirá Matrix.txt en cada pareja, por lo que al final, solo quedará el resultado de la pareja final.

Quizás sea eso: debes ir cambiando el nombre del archivo de salida.

por **lis** » 2013-03-15 08:33 @398

Hola, explorer, muchas gracias por responder. Bueno, lo que hice fue lo siguiente:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl -w
 
my %hash
    = qw (1 15 16 31 32 45 46 59 60 74 75 90 91 105 106 120 121 135 136 151 152 166 167 181 182 196 197 212 213 227 228 243 244 258 259 273 274 288 289 304 305 319 320 334 335 349 350 365);
 
for my $jul ( sort { $hash{$a} <=> $hash{$b} } keys %hash ) {
    $file = $hash{$jul} . '.txt';
 
    open( OUT, ">>Matriz.txt" );
 
    open( IN, ">matriz.m" );
    print IN <<EOF;
clear all; clc; 
x= load('prueba.txt');
a=x(:,1);   m=x(:,2);  d=x(:,3);  h=x(:,4);  vv=x(:,5);  dv=x(:,6);
V=reshape(x(:,5)',365,24); D=reshape(x(:,6)',365,24);
mv=mean(V);  md=mean(D);
G2R = (pi / 180); gar = G2R * D;
u = (-1) .* V .* sin(gar); v = (-1) .*V .* cos(gar);
prom=[];
for K = 1:365
Samp{K} = mean(V($jul:$hash{$jul},:));
end;
format bank;
prom=[prom; Samp{K}];
save $file prom '-ascii';
EOF
    close IN;
    system("/usr/bin/octave -q matriz.m");
 
    open( GR, "<$file" );
    while ( $lineas = <GR> ) {
        chop($lineas);
        @datos = split( " ", $lineas );
        print OUT "@datos\n";
    }
}
 
 
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

Mezclé de los dos scripts, ya que

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

say
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

no me funcionó. Al final resulta lo que quería pero quiero ver si al menos la última parte (donde guardo los archivos en uno solo) se pueda simplificar. Así que nuevamente muchas gracias

por **explorer** » 2013-03-15 20:17 @887

say() no es más que un print(), y está disponible a partir de Perl v5.10.

La última parte creo que se puede simplificar bastante:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl
use strict;
use warnings;
use diagnostics;
 
my %hash
    = qw (1 15 16 31 32 45 46 59 60 74 75 90 91 105 106 120 121 135 136 151 152 166 167 181 182 196
            197 212 213 227 228 243 244 258 259 273 274 288 289 304 305 319 320 334 335 349 350 365);
 
open( OUT, ">>Matriz.txt" );
 
for my $jul ( sort { $hash{$a} <=> $hash{$b} } keys %hash ) {
    my $file = "$hash{$jul}.txt";
 
    open( IN, ">matriz.m" );
    print IN <<EOF;
clear all; clc; 
x= load('prueba.txt');
a=x(:,1);   m=x(:,2);  d=x(:,3);  h=x(:,4);  vv=x(:,5);  dv=x(:,6);
V=reshape(x(:,5)',365,24); D=reshape(x(:,6)',365,24);
mv=mean(V);  md=mean(D);
G2R = (pi / 180); gar = G2R * D;
u = (-1) .* V .* sin(gar); v = (-1) .*V .* cos(gar);
prom=[];
for K = 1:365
Samp{K} = mean(V($jul:$hash{$jul},:));
end;
format bank;
prom=[prom; Samp{K}];
save $file prom '-ascii';
EOF
    close IN;
 
    system("/usr/bin/octave -q matriz.m");
 
    open( GR, "<$file" );
    print OUT <GR>;
    close GR;
}
 
close OUT;
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

ya que, realmente, lo que hacemos es copiar las líneas generadas por octave, al archivo Matrix.txt.

El proceso que hacías de, quitar el último carácter, dividir en @campos separando por los espacios en blanco, y sacando el resultado afuera reuniendo los campos con espacios en blanco, y añadiendo un carácter de fin de línea... pues me parece que sobraba.

Se podría simplificar más si supiéramos cómo decirle a octave que la salida la haga hacia la salida estándar en lugar de a un archivo, o que en vez de escribir a un archivo, vaya agregando la salida.

por **lis** » 2013-03-19 14:42 @654

Hola, explorer, muchas gracias por tu respuesta.

He estado buscando y en vez de crear un archivo para cada resultado y luego unirlos en uno solo, hallé la forma de decirle a octave que lo haga pero me genera un error al escribir lo siguiente:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

open( IN, ">matriz.m" );
print IN <<EOF;
clear all; clc;
x= load('$file');
a=x(:,1);   m=x(:,2);  d=x(:,3);  h=x(:,4);  vv=x(:,5);  dv=x(:,6);
V=reshape(x(:,5),24,366)'; D=reshape(x(:,6),24,366)';
mv=mean(V);  md=mean(D);
G2R = (pi / 180); gar = G2R * D;
u = (-1) .* V .* sin(gar); v = (-1) .*V .* cos(gar); t=1:24;
fid = fopen('$name','wt');
"for K = 1:365
Samp_$variable{K} = mean($variable($jul:$hash{$jul},:));
fprintf('%2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f\n',Samp_$variable{K});
end;
EOF
close IN;
system("/usr/bin/octave -q matriz.m");
 
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

El problema radica en que \n genera un salto de línea en Perl, y no en Octave, por lo que el archivo matriz.m no lo reconoce Octave ya que queda así:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

fprintf('%2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f %2.1f
',Samp_);
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

¿Sabes cómo solucionar esto?

De antemano muchas gracias.

por **explorer** » 2013-03-19 16:36 @733

Prueba a escribirlo así: \\n

por **lis** » 2013-03-20 11:07 @505

Muchas gracias, explorer, funcionó a la perfección.

Foro - Perl en Español

Promedio con rangos para una matriz con Perl y Octave

Promedio con rangos para una matriz con Perl y Octave

Publicidad

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

Re: Promedio con rangos para una matriz con Perl y Octave

¿Quién está conectado?