「Bash」カテゴリーアーカイブ

Bash

paste -s

コメントする

テキストファイルから情報を抜き出してx, yのテキストデータを作りたいのだが、
>> のように行を追加するコマンドはあっても

x

というファイルに y を追加して

x y

というファイルを作るコマンドはない…？

仕方ないので、

x
y

というファイルxy.txtを

paste -s xy.txt

で2列のファイルにする。これを繰り返して

x1 y1
x2 y2
x3 y3
…

とデータを追加していく。中間ファイルができてしまうのが気持ち悪い…

Bash

数値の取り出しと演算

コメントする

ファイル中一番最後の原子のz座標を取り出して足し算したいとき。
awkなら読み取った数値（本来は文字列）を自動的に浮動小数点として扱ってくれる。
ただし、表示桁数が少ないので、OFMTを使って増やす。

tail -1 dum.xsf | awk '{ OFMT = "%.8f"} {print $NF+5.95907033}'
Bash

awkで最大値

コメントする

a 11
b 43
c 89
d 62
e 35

こんなデータがhoge.txtとしてあるとき、最大値やら平均やらを出力するコマンドはすぐに見つかったが、最大値を取るアルファベットを出力する場合に少し手こずったので書いておく。

まずは最大値を出力する場合。
$ cat hoge.txt | awk ‘{if(m<$2) m=$2 } END{print m}' 次に、アルファベットを出力する場合。 $ cat hoge.txt | awk '{if(m<$2) {m=$2; n=$1} } END{print n}' 要するにif文の後の実行文が二つになるので、中括弧を一つ増やせばよかっただけでした。

Bash, Gnuplot

データ解析（bash&gnuplot）

コメントする

磁場を振って測った電流のデータを、素子の劣化を考慮して直線補正するためのスクリプト。

まずは、カンマをスペースに置換し、200行毎に空行を２行挟む（gnuplotのindexでデータを分けるため）。

直線によるフィッティングの結果からデータを補正するのだが、
普通に$1、$2とやろうとするとbashの変数と混同するため、
あらかじめ$1を文字列と解釈してくれる変数aを定義しておく。
最後に、補正したデータをテーブルにダンプして終わり。

#!/bin/bash

if [ $# -ne 1 ]; then
  echo "please add the filename!"
  exit 1
fi

tail -701 $1 | sed -e 's/,/   /g' | sed -e '201~200i \\n' > test.txt

a=\$1
b=\$2

gnuplot -e "
  f(x)=a*x+b ; g(x)=c*x+d; h(x)=e*x+f;
  a=1;b=1;c=1;d=1;e=1;f=1;
  fit f(x) 'test.txt' i 0 u 1:2 via a,b;
  fit g(x) 'test.txt' i 1 u 1:2 via c,d;
  fit h(x) 'test.txt' i 2 u 1:2 via e,f;
  unset key;
  set xlabel 'Magnet Voltage (V)';
  set ylabel 'Current (A)';
  plot 'test.txt' i 0 u 1:($b-$a*a) w l, '' i 1 u 1:($b-$a*c) w l, '' i 2 u 1:($b-$a*e) w l;
  pause -1 ;
  set table 'test.table';
  replot;
  "
Bash

sed

コメントする

ある条件に一致する行を見つけたら、その下の固定行数を抜き出したい場合。
sedを使うとそれなりに簡単にできることがわかった。
例は,i= 7を含む行と、加えて241行を抜き出す。

sed -n -e "/ i =    7/,+241p" pos.xyz >> pos-7.xyz

最後の行を削除する。

sed -i -e '$d' foo.txt

-i　オプションはファイルを上書きする命令。

Bash

ゼロ埋め番号づけ

コメントする

lsで見たときにファイルやディレクトリを順番に並べるには、
ゼロ埋めするとよい（1,2,3…でなはく、01,02,03…とする）。スクリプトで巡回するには下のようにする。

  for i in {1..32} ; do
    num=$(printf %02d $i)
    echo $num
  done