dipole90.plt

这是一个 York大学开发的FDTD代码。很好用的工具

#    G N U P L O T
#    MS-Windows 32 bit version 3.5 (pre 3.6)
#    patchlevel beta 340
#    last modified Tue Nov 25 22:57:44 GMT 1997
#
#    Copyright(C) 1986 - 1993, 1997
#    Thomas Williams, Colin Kelley and many others
#
#    Send comments and requests for help to [e-mail]
#    Send bugs, suggestions and mods to [e-mail]
#
# set terminal windows color "Arial" 10
# set output
set noclip points
set clip one
set noclip two
set bar 1.000000
set border 31
set xdata
set ydata
set zdata
set x2data
set y2data
set boxwidth
set dummy x,y
set format x "%g"
set format y "%g"
set format x2 "%g"
set format y2 "%g"
set format z "%g"
set nogrid
set key title ""
set key right top Right noreverse box -2 samplen 4 spacing 1 width 0
set nolabel
set noarrow
set nolinestyle
set nologscale
set offsets 0, 0, 0, 0
set pointsize 1
set encoding default
set nopolar
set angles radians
set noparametric
set view 60, 30, 1, 1
set samples 100, 100
set isosamples 10, 10
set surface
set nocontour
set clabel '%8.3g'
set nohidden3d
set cntrparam order 4
set cntrparam linear
set cntrparam levels auto 5
set cntrparam points 5
set size ratio 0 1,1
set origin 0,0
set data style lines
set function style lines
set xzeroaxis 0
set x2zeroaxis -2
set yzeroaxis 0
set y2zeroaxis -2
set tics in
set ticslevel 0.5
set ticscale 1 0.5
set mxtics default
set mytics default
set mx2tics default
set my2tics default
set xtics border mirror norotate 
set ytics border mirror norotate 
set ztics border nomirror norotate 
set nox2tics
set noy2tics
set title "Ez field at 11,20,11 - 9m from source - Theta = 90deg." 0.000000,0.000000  ""
set timestamp "" bottom norotate 0.000000,0.000000  ""
set rrange [ * : * ] noreverse nowriteback  # (currently [0:10] )
set trange [ * : * ] noreverse nowriteback  # (currently [-5:5] )
set urange [ -5 : 5 ] noreverse nowriteback
set vrange [ -5 : 5 ] noreverse nowriteback
set xlabel "Frequency (MHz)" 0.000000,0.000000  ""
set x2label "" 0.000000,0.000000  ""
set timefmt "%d/%m/%y\n%H:%M"
set xrange [ * : * ] noreverse nowriteback  # (currently [0:49.965] )
set x2range [ * : * ] noreverse nowriteback  # (currently [-10:10] )
set ylabel "Amplitude (dBV/m)" 0.000000,0.000000  ""
set y2label "" 0.000000,0.000000  ""
set yrange [ * : * ] noreverse nowriteback  # (currently [-51.8166:-35.2177] )
set y2range [ * : * ] noreverse nowriteback  # (currently [-10:10] )
set zlabel "Z - Mesh Units" 0.000000,1.000000  ""
set zrange [ * : * ] noreverse nowriteback  # (currently [0:21.1] )
set zero 1e-008
set lmargin -1
set bmargin -1
set rmargin -1
set tmargin -1
set locale "C"
Eth(d,f,th)=dl**2*f*sin(th*pi/180)/(2*c0*d)
dB(x)=20*log10(abs(x))
dl = 1
c0 = 300000000.0
MHz = 1000000.0
plot 'fdez.dat' title "TLM",dB(Eth(9,x*MHz,90)) title "Theory"
#    EOF