educhan, TSAS and LABPACK etc.
a2ps 入力ファイル > 出力ファイルという形で使います。
gs ポストスクリプトファイル名
lp -d octopus ポストスクリプトファイル名bass では、従来通り ftp でファイルを送っても印刷できます。
lp -d bass ポストスクリプトファイル名
ftp bass
と入力します。anonymous-FTPとは違い、ログイン名やパスワードは
必要ありません。Connected to bass. 220-JC-CONNECT PowerPrint M-2000 Version 1.00 220-Copyright (C) 1995 ACCESS CO.,LTD. 220 Copyright (C) 1995 Japan Computer Industry Inc. Name (bass:hidaka): <-- 何も入力せずに 230 USER hidaka logged in.リターンのみ。 ftp> cd lp_p1 <-- 必ずこのディレクトリ 250 CWD command successful. に変更。 ftp> put ファイル名 <-- 200 PORT command successful. 150 ASCII data connection established. 226 Transfer complete. local: akay.ps remote: akay.ps 8730 bytes sent in 0.032 seconds (2.7e+02 Kbytes/s) ftp> quit <-- 221 Quit.このようにして、プリンタにファイルを送って印刷します。
Ndata =550
Xmean = 1086.10913
Xs.d. = 47.67794
Freq Sxx S2.0 |Sxx'| S0.5 CGSxx
0.00000 283.92383 0.00003 283.92383 0.00002 0.00000
0.00719 67.29022 161.01050 67.29022 109.46802 0.00000
0.01439 99.25716 119.27703 99.25716 133.95329 0.00000
0.02158 107.26489 113.39518 92.25481 75.05566 15.01009
0.02877 58.15707 117.13137 58.15707 55.42382 0.00000
: : : : : :
最初の4行がデータでなくヘッダであることをgnuplotに教えるには、最初の
4行のみ、行頭に#をつける必要があります。vi で編集してもよいのですが、
大量のファイルを一度に処理するためには、awkを使います。
{
if(NR <= 4){ printf("#%s\n", $0) }
else { print $0 }
}
educhanのプロンプトから、
gawk -f header.awk sample.dat > sample2.datと入力すると、#をつけたものが別名のファイルで保存されます。
<tkweb@educhan:293>$ gnuplot G N U P L O T unix version 3.5 patchlevel 3.50.1.17, 27 Aug 93 last modified Fri Aug 27 05:21:33 GMT 1993 Copyright(C) 1986 - 1993 Thomas Williams, Colin Kelley Send comments and requests for help to info-gnuplot@dartmouth.edu Send bugs, suggestions and mods to bug-gnuplot@dartmouth.edu gnuplot> plot "sample2.dat" using 1:6 with lines gnuplotすると、このように描画されます。
画像をPostscript形式で保存したい場合には、gnuplotにもどって、
gnuplot> set terminal postscript Terminal type set to 'postscript' Options are 'landscape monochrome dashed "Helvetica" 14' gnuplot> set output "sample.ps" gnuplot> plot "sample2.dat" using 1:6 with linesとします。終了するには、quit と入力します。
Ndata = 731
Xmean = 821.12994 Ymean = 120.43776
Xs.d. = 31.05704 Ys.d. = 5.84055
Time Freq Sxy Phase
58.30023 0.00000 772.35907 2.53821
58.30023 0.00951 869.93903 1.72086
58.30023 0.01903 608.97809 0.83454
58.30023 0.02854 293.71045 0.61633
58.30023 0.03806 199.28731 0.36296
: : : : : :
58.30023 0.58989 242.67116 -2.86146
58.30023 0.59940 468.79257 -0.00136
59.12136 0.00000 833.32263 2.57939
59.12136 0.00951 946.70050 1.75705
59.12136 0.01903 684.94086 0.85450
gnuplotで3Dグラフを書くためには、入力ファイルの形式が、
time1 x11 y11 time1 x12 y12 time1 x13 y13 time2 x21 y21 time2 x22 y22 time2 x23 y23 time3 x31 y31 time3 x32 y32 time3 x33 y33のように、time1, time2, time3の間に空行がなければなりません。 この用に整形するためのawkスクリプトは、以下のようになります。
{
if (NR == 5){
time = $1
print $0
}
if (NR > 5){
if ($1 == time){
print $0 }
if ($1 != time){
time = $1
printf("\n")
print $0 }
}
}
次にgnuplotを起動して、さまざまなオプションを設定した後に、
splotコマンドで3D描画しなければならないのですが、これを
ひとつのファイル(test.plt)に記述して済ませることができます。
set hidden set ticslevel 0.5 set view 30, 120 set autoscale set logscale z set parametric set xlabel " Time (sec) " set ylabel " Frequency (Hz) " set data style lines set title "Wigner-Ville Distribution" splot "sample2.wvd" pause -1最後の行の"pause -1"は、「リターンキーを押すまで画像を表示する」 という意味です(このままだと描画後すぐに終了してしまうので)。
ksh: sas:not foundと表示される場合には、環境設定ファイル.kshrcの中で、 環境変数PATHを変更する必要があります。.kshrcの 編集には、viなどのテキストエディタを使います。
PATH=.:/usr/bin:/usr/ucb:/usr/5bin:/opt/NeWSprint/bin:/usr/ccs/binこれを、
PATH=.:/usr/bin:/usr/ucb:/usr/5bin:/opt/NeWSprint/bin:/usr/ccs/bin:/usr/local/bin/sas609のように変更すると、sasを起動することができるようになります。
注：以下のURLは、ユーザー側へ連絡なしに変更されることが多々あります。 その場合には、東京大学や 東京大学付属図書館、 各部局のホームページから新しいURLを検索して下さい。