make_fig3.sh

计算高斯各阶导函数的C程序 Computing Gaussian derivative waveforms of any order. Dgwaveform efficiently cre

#! /bin/sh
# Reproduces Fig. 3 in manuscript
# Author: Werner M. Heigl


# Make 425 kHz pulse

   ./dgwaveform fpeak=425 n=2 sign=-1 shift=0.0385 |
   
# and shift vertically:

   suhtmath key=cdp op=add const=2.5 >d2g_425kHz.su
   
# Make 300 kHz pulse

   ./dgwaveform fpeak=300 n=2 sign=-1 shift=0.0375 nt=285 |
   
# and resample to same dt as 425 kHz pulse:

   suresamp dt=0.000294 nt=285 >d2g_300kHz.su
   
# Concatenate into one file

   cat d2g_300kHz.su d2g_425kHz.su |
   
# and set ntr:

   sushw key=ntr a=2 >ricker.su

# make Postscript file
style="normal"
x1beg=0 x1end=0.1 d1num=0.02
x2beg=-1.5 x2end=4 d2num=10 f2num=-2
label1="Time (ms)" label2="Normalized Amplitude" labelsize=16
wbox=6 hbox=4

   supsgraph <ricker.su >figure_3.eps style=$style \
              x1beg=$x1beg x1end=$x1end d1num=$d1num \
	      x2beg=$x2beg x2end=$x2end d2num=$d2num f2num=$f2num \
	      label1="$label1" label2="$label2" \
	      labelsize=$labelsize \
	      wbox=$wbox hbox=$hbox
	   
exit 0