петля по линиям в файле и вычесть предыдущая линия от текущей строки

$ awk 'NR > 1 { print $0 - prev } { prev = $0 }' <file.dat
-0.043962
-0.020674
0.002742
-0.004449
0.011906
-0.026333
0.001461
0.011245
-0.002031
-0.009281
-0.006112
0.011811
-0.004837

Делать это в цикле оболочки с вызовом bcнеудобно. В приведенном выше примере используется простой скрипт awk, который считывает значения из файла одно за другим и для любой строки после первой печатает разницу, как вы описываете.

Первый блок, NR > 1 { print $0 - prev }, условно печатает разницу между этой и предыдущей строкой, если мы достигли второй строки или дальше(NR— это количество прочитанных записей, а «запись» по умолчанию строка ).

Второй блок, { prev = $0 }, безоговорочно устанавливает prevна значение в текущей строке.

Перенаправить вывод на newfile.dat, чтобы сохранить там результат:

$ awk 'NR > 1 { print $0 - prev } { prev = $0 }' <file.dat >newfile.dat

Связанные:

• Почему использование цикла оболочки для обработки текста считается плохой практикой?

Было некоторое упоминание о медленности вызова bcв цикле. Ниже приведен способ использования одного вызова bcдля выполнения арифметических операций при чтении данных в цикле оболочки. (На самом деле я бы не рекомендовал решать эту проблему таким образом, и я только показываю это здесь. для интересующихся со -процессами вbash):

#!/bin/bash

coproc bc

{
    read prev

    while read number; do
        printf '%f - %f\n' "$number" "$prev" >&"${COPROC[1]}"
        prev=$number

        read -u "${COPROC[0]}" result
        printf '%f\n' "$result"
    done
} <file.dat >newfile.dat

kill "$COPROC_PID"

Значение в ${COPROC[1]}является стандартным дескриптором входного файла bc, а ${COPROC[0]}— стандартным дескриптором выходного файла bc.

Если вы хотели попытаться заставить скрипт оболочки работать, вам просто не хватило инициализации:

f=myfile.dat
prev=0
while read line; do
    bc <<< "$line - $prev"
    prev=$line
done < $f > newfile.dat

... где я также переместил перенаправление за пределы цикла, просто чтобы сэкономить часть операций ввода-вывода.

Решение bcне печатает ведущие нули, в то время как решение awk печатает.

Вы можете использовать перенаправление execдля чтения последовательных строк входного файла из нескольких точек в сценарии -один раз перед циклом (для установки начального значения ),затем многократно в течение этого (для каждого нового значения вычитать):

exec 3<file.dat
read prev<&3
while read curr ; do
        bc <<< "$curr - $prev" >> newfile.dat
        prev=$curr
done <&3

Использование некоторых простых утилит GNU и отсутствие циклов оболочки:

paste -d- <(head -n-1 file.dat) <(tail -n+2 file.dat) | bc

Идея состоит в том, чтобы разбить входной файл на два столбца; сместите второй столбец на 1 строку и вставьте столбцы вместе с -в качестве разделителя. headи tailиспользуются для обрезки последней строки 1-го столбца и первой строки 2-го столбца соответственно для достижения необходимого смещения. Результирующий список является обязательным списком арифметических разностей, который передается в bcдля оценки.

Попробуйте онлайн .

В качестве альтернативы, если вам нравится sed, вы можете сделать это:

sed '1{s/$/-\\/;p;d};${p;d};s/.*/&\n&-\\/' file.dat | bc

Это дублирует каждую строку и вставляет -\в конце второй версии каждой строки. Первая и последняя строки обрабатываются по-разному для создания необходимых выражений. Вывод sed выглядит примерно так:

a-\
b
b-\
c
c-\
d

Это снова действительные арифметические разности, которые bcможно вычислить. Не то, чтобы bcпонимал обратную косую черту продолжения строки -на концах каждой второй строки.

Попробуйте онлайн .

Я использую массивы. Я использую их для всего. Я не могу вспомнить, как работают awk и sed, без тщательного изучения справочных страниц. Вот как бы я это сделал.

f=( $(< file.dat) )
for ((num=1;num<=${#f[@]};num++))
do
    echo $(bc <<< ${f[$num]}-${f[(($num-1))]})>>differences.dat
done

Я так понимаю. Он имеет неприятные особенности некоторых других ответов :зацикливание и вызов bc снова и снова. Однако он читает файл только один раз, как и ответы с использованием sed и awk.

Вы можете попробовать это

num <- as.data.frame(num)
num$sub_num <- num[c(2:14, c("0")), ]
num$diff <- num$num - num$sub_num