Dotaz: Bash aritmeticky priemer

[Mountains]

Zdravim, potreboval by som urobit skript ktory nacita hodnoty zo suboru a urobi s nich aritmeticky priemer. Dakujem za pomoc.

Odpovědi

[Jan Martinek]

Tady máš tři různé algoritmy včetně generování testovacího souboru a měření doby výpočtu.

#!/usr/bin/python
import random, time

print 'Wait ...'
fw = file('cisla.txt', 'w')
for i in range(1000000):
    fw.write('%f\n' %(random.random()))
fw.close()

start = time.time()
s = 0
lines = 0
for line in file('cisla.txt'):
    s += float(line)
    lines += 1
print s/lines
print time.time() - start

start = time.time()
s = 0
for i, line in enumerate(file('cisla.txt')):
    s += float(line)
print s/(i + 1)
print time.time() - start

start = time.time()
numbers = [float(line) for line in file('cisla.txt')]
print sum(numbers)/len(numbers)
print time.time() - start

[Michal Kubeček]

Algoritmus vidím jen jeden, zato ve třech implementacích (z toho jedné tak zoufale paměťově neefektivní (aniž by to bylo vyváženo časovou náročností), že by se to dalo používat při výuce jako odstrašující příklad).

[Jan Martinek]

Pro případ, že by se někdo špatně vyspal, uvádím ještě tohle (bez ošetřování chyb):

#include<stdio.h>

int main(void){
    FILE* fr;
    double f, s = 0, lines = 0;

    fr = fopen("cisla.txt", "r");
    while (fscanf(fr, "%lf\n", &f) != EOF){
        s += f;
        lines += 1;
    }
    printf("%12.12f\n", s/lines);
    fclose(fr);
    return 0;
}

Je to cca 3-4 krát rychlejší než v pythonu.

[Mountains]

Dakujem za Vasu ochotu a pomoc

[Semo]

No, kedze su prazdniny, tak to domaca asi nebude, hoci to inak vykazuje vsetky jej znaky.

Najvacsi problem je dostat hodnotu z podshellu von a trik je spocitat priemer na desatinne miesta, hoci to bash nevie.

#!/bin/bash
set -- `while read n; do
  sucet=$((sucet + n))
  lines=$((lines + 1))
  echo $sucet $lines
done <subor.in | tail -n -1`

sucet=$1
lines=$2

echo $((sucet / lines))
# # ak to chces na 2 destinne miesta:
# avgfrac=$((sucet * 100 / lines))
# echo ${avgfrac%??}.${avgfrac: -2}

[lada]

pokud to ma byt ciste v shellu tak asi takto: #!/bin/sh

text=`cat soubor_dat |tr '\n' '+'`0
pocet=`cat soubor_dat |wc -l`
prumer=`echo "($text)/$pocet"|bc -l`
echo $prumer


[Michal Kubeček]

Tak tomu říkám řešení čistě v shellu… :-)

[Zdenek]

Zatim je to ovsem nejbliz k zadani. A nepochybne to spravne *nix reseni

[Semo]

Potrebuje to tr a bc. Riesenie CISTO v bashi mam zatial iba ja

[Jan Martinek]

Dá se to napsat třeba i takhle:

#!/bin/awk -f
BEGIN{lines=0;s=0;}
{lines++;s+=$1;}
END{print "prumer je", s/lines;}

no a awk už tak skoro "v shellu" je, že? Nebo teda v posixu. Teda *nixu. Teda ... no prostě je to víc v shellu než třeba bc. A určitě je to mnohem blíž shellu než python. Ale asi je to trochu dál shellu než třeba "tail". To jsou problémy ...

[Michal Kubeček]

Já si hlavně myslím, že zadání typu "udělejte XY v ZW" jsou dobrá akorát tak pro školní úlohy. V praxi je daleko užitečnější umět zvolit vhodný nástroj pro danou úlohu (zde podle okolností např. awk, bc, perl, python nebo C), než se snažit ji za každou cenu udolat nevhodně zvoleným nástrojem.

[pavel]

Já řeším podobný problém, ale data mám ve sloupcích s oddělovačem ; Jak sčítat data ve sloupci?

[Michal Kubeček]

awk(1)

[Jan Martinek]

Takže bude tolik průměrů, kolik je sloupců?

#!/usr/bin/python
# -*- coding:utf-8 -*-
import random, time

radky = 5
sloupce = 10

fw = file('cisla.txt', 'w')
for i in range(radky):
    cisla_str = ['%f' %(i*j) for j in range(sloupce)]
    radek = ';'.join(cisla_str)
    fw.write('%s\n' %radek)
fw.close()

fr = file('cisla.txt')
sums = [float(x) for x in fr.readline().split(';')]
for i, line in enumerate(fr):
    for j, val_str in enumerate(line.split(';')):
        sums[j] += float(val_str)
i += 2
avgs = [x/i for x in sums]
print sums
print avgs

A na začátku se neví, kolik těch sloupců je, že?

[pavel]

Sloupců se ví (je to druhý až patnáctý), ale počet řádků mám pokaždé jiný...

[Jan Martinek]

Aha. Ale snad by mělo jít použít to, co jsem napsal, protože to funguje univerzálně pro libovolný (rozumný) počet sloupců i řádků. Samo si to příslušná čísla zjistí.