<Data> -> <Load Waves> -> <Load General Binary File...>
Input File
Data Type: 16 bit signed integer
Byte to skip at start of file: 0
Number of arrays in the file: 10
Number of points in array: auto
Byte order: Low byte first
Points in file are interleaved: no check
Path: Fileをクリックするとディアログが開くので、ダウンロードしておいたfakeEPSC.datを選択。
Output Waves
Data Type: Double float
Base name: w
Overwrite existing waves: check
Scaling: check
Scaling Offset: 0
Scaling multiplier: 0.0625
これでDo itすると、データが読み込まれ、w0からw9までのwavesが作成されます。
実際のコマンドとしては、
GBLoadWave/O/B/N=w/T={16,4}/Y={0, 0.0625}/W=10 "(ここはファイルパス):fakeEPSC.dat"
が使われています。
<Windows> -> <New Graph...>
Y Wavesの欄で、w0からw9までを選択。
X Wavesの欄は、_calculated_を選択。
そしてDo itすると、グラフが現れます。
使用されたコマンドは、
Display w0,w1,w2,w3,w4,w5,w6,w7,w8,w9
です。
波形はそれらしい形をしています。amplitude（Y軸の大きさ）もよいとしましょう。しかしX軸のスケールがなされていないのがわかります。
<Data> -> <Change Wave Scaling...>
Set X Property: check
Unit Type: Numeric
Unit: ms
Start: 0
Delta: 0.05
SetScale Mode: Start and Delta
Wavesの欄でw0からw9までを選ぶ。
そしてDo Itする。
使用されたコマンドは、
SetScale/P x 0,0.05,"ms", w0,w1,w2,w3,w4,w5,w6,w7,w8,w9
先ほど作成したグラフのX軸の目盛が0 から100 msまでに変わっているはずです。
duplicate/o/d w0, wt
そしてこのトレースを表示します。
display wt
Baselineの値として、X軸の0から10までの範囲のデータの平均を使用することにします。
<Analysis> -> <Wave Stats>
Wave: wt
X Range: 0と10
これでDo It。
使用されているコマンドは、
WaveStats/R=(0,10) wt
です。 いろいろ数値が表示されますが、平均値はV_avgです。
この値を変数baseとして保存しておきます。
variable base = V_arg
この部分は、
variable base
base = V_avg
と2行に分けて書くの同じです。（変数は先ず宣言しておくことが必要です。数値の場合はvariable、文字列の場合はstringで宣言します。）
WaveStats/R=(10,20) wt
で、最小値はV_minです。応答の大きさは、V_min - baseで、この値をamp1として保存します。
variable amp1
amp1 = V_min - base
WaveStats/R=(40,50) wt
variable amp2
amp2 = V_min - base
variable PPR
PPR = amp2/amp1
データの確認には、
print amp1, amp2, PPR
WaveStats/R=(0,10) wt
base = V_avg
WaveStats/R=(10,20) wt
amp1 = V_min - base
WaveStats/R=(40,50) wt
amp2 = V_min - base
PPR = amp2/amp1
print amp1, amp2, PPR
となります。
<Windows> -> <Procedure Window>
を開いて、次の10行を書き入れる。
macro pprCalculate()
WaveStats/R=(0,10) wt
base = V_avg
WaveStats/R=(10,20) wt
amp1 = V_min - base
WaveStats/R=(40,50) wt
amp2 = V_min - base
PPR = amp2/amp1
print amp1, amp2, PPR
end
書き終わったら、下のcompileボタンをクリック。問題がなければ、ボタンが消えるはずです。なお、行の始めを字下げしているのは、見易くするためです。
<Macro> -> <pprCalculate>
もしくはコマンドから、pprCalculate()と入力し、Enter。
duplicate/o/d w1, wt
pprCalculate()
function pprCalculate()
wave w = $("wt")
variable base, amp1, amp2, PPR
WaveStats/Q/R=(0,10) w
base = V_avg
WaveStats/Q/R=(10,20) w
amp1 = V_min - base
WaveStats/Q/R=(40,50) w
amp2 = V_min - base
PPR = amp2/amp1
print amp1, amp2, PPR
end
少し説明を付け加えると、functionaの中では、これまで使用してきたwaveとvariableをそのまま使用することが出来ません。functionの中では、function固有のvariable、waveを宣言しなくてはなりません。そのためbaseやamp1などの変数を宣言しています。またwave w = $("wt")は、wtという名前のwaveをこのfunctionの中では、wと呼びますという宣言です。
WaveStatsのQオプションはQuietオプションで、結果を表示しないようにします。
入力したらCompileボタンをクリック。問題がなければボタンは消えるはずです。
そしてコマンドからpprCalculate()と入力してEnter
function pprCalculate(n)
variable n
string s
wave w = $("wt")
variable base, amp1, amp2, PPR
s = "duplicate/o/d w" + num2istr(n) + ",wt"
execute s
WaveStats/Q/R=(0,10) w
base = V_avg
WaveStats/Q/R=(10,20) w
amp1 = V_min - base
WaveStats/Q/R=(40,50) w
amp2 = V_min - base
PPR = amp2/amp1
print amp1, amp2, PPR
end
pprCalculate(0)、pprCalculate(3)などを試してみる。
説明が必要なのは、s = ....の部分です。duplicateはfunction内ではそのまま使用できません。しかしexecuteをしようすることにより、コマンド入力に相当することが出来ます。たとえば、w7をコピーしたい場合、duplicate/o/d w7, wtというコマンドが必要です。num2istr()は数値を文字に変換する関数で、num2istr(7)で文字"7"を作れますここでは、duplicate...の部分に"7"を加え、さらに", wt"を加えてコマンドの文字列を作成しています。
新たなfunctionを作成します。（pprCalculate()の下に書き加える。）
function allCalc()
variable i
for(i=0;i<10;i+=1)
pprCalculate(i)
endfor
end
そしてallCalc()と入力してEnter.
function pprCalculate(n)
variable n
string s
wave w = $("wt")
wave amp1 = $("vAmp1")
wave amp2 = $("vAmp2")
wave PPR = $("vPPR")
variable base
s = "duplicate/o/d w" + num2istr(n) + ", wt"
execute s
WaveStats/Q/R=(0,10) w
base = V_avg
WaveStats/Q/R=(10,20) w
amp1[n] = V_min - base
WaveStats/Q/R=(40,50) w
amp2[n] = V_min - base
PPR[n] = amp2[n]/amp1[n]
print n, amp1[n], amp2[n], PPR[n]
end
function allCalc()
variable i
make/o/d/n=10 vAmp1
make/o/d/n=10 vAmp2
make/o/d/n=10 vPPR
for(i=0;i<10;i+=1)
pprCalculate(i)
endfor
end
display vPPR
適当にグラフの表示方法を変更
