Dotaz: Gnuplot - součet křivek

[Martin Čížek]

DD, nevíte zda lze v GNUplotu nějak provést součet křivek? Mám několik souborů s daty (x, y); křivka (funkce) z každého vzniká spojením čarami ("with lines"). Takto vzniklé funkce bych potřeboval sečíst. Problém je, že se hodnoty x v jednotlivých souborech nekryjí.

Pokud to Gnuplot neumí, byl bych vděčný za nasměrování na nějaký skriptík, který datové řady umí sečíst (sice bych si jej uměl si jej napsat, ale teď mne trochu tlačí čas).

Díky.

Odpovědi

[Martin Čížek]

Tak jsem si to nakonec musel naprgat. Pokud by se to někomu hodilo, zde to je... Po aplikaci stačí na výstup v gnuplotu pustit plot 1:($2+$3+...)

#!/bin/bash
# USAGE: mergefuncs FILE...
# Martin Cizek <cizek@webland.cz>

[ $# -gt 0 ] || exit 1;

awk '
function complete_last()
{
    for (i = idx; i <= count; i++)
        res[i] = res[i] " 0"
    idx = 1
    have_prev = 0
}
BEGIN {
    xi = 0
    idx = 1
}
ARGIND < ARGC/2 {
    if (length($1))
        x[xi++] = (double)$1
    next
}
!sorted {
    count = asort(x)
    sorted = 1
    lastind = ARGIND
}
/[^ \t]/{
    if (lastind != ARGIND)
        complete_last()
    lastind = ARGIND
    while (idx <= count && x[idx] <= $1) {
        if ((have_prev && x[idx] == prev_x) || x[idx] == $1)
            res[idx] = res[idx] " " $2
        else if (have_prev && x[idx] >= prev_x)
            res[idx] = res[idx] " " (prev_y + ($2-prev_y)*(x[idx]-prev_x)/($1-prev_x))
        else
            res[idx] = res[idx] " 0"
        idx++;
    } 
}
/[^ \t]/&&!have_prev {
    have_prev = 1
    prev_x = (double)$1
    prev_y = (double)$2
}
END {
    complete_last()
    for (i = 1; i <= count; i++)
        print x[i] " " res[i]
}
' "$@" "$@"

[Hynek (Pichi) Vychodil]

a nevyřešilo by to

join -e 0 -j 1 files...

[Martin Čížek]

Bohužel ne, protože hodnoty x jsou v různých souborech opravdu různé (jsou to časy zachycení paketu). Pro každou x-hodnotu z libovolného souboru se musí na výstup dát také lineární interpolace od všech ostatních funkcí.

Nicméně díky za tip...

[Hynek (Pichi) Vychodil]

jojo, fakt to nejde. Jde to kombinací cut, sort -q, join ale to je pak lepsi to udelat nejakym scriptem (ja tedy asi perlem).