Platon Technologies
not logged in Login Registration
EnglishSlovak
open source software development celebrating 10 years of open source development! Friday, March 29, 2024

File: [Platon] / doc / diplomova-praca-rajo / images / vdrift_T-300K-normal-model.gpi (download)

Revision 1.2, Sun Jan 16 15:01:31 2005 UTC (19 years, 2 months ago) by rajo


Changes since 1.1: +2 -2 lines

Fix: gnuplot 4.0 has different syntax (incompatible with gnuplot 3.7.x)
Change lines from:
    set linestyle 1 linetype 7 linewidth 1
to:
    set line style 1 linetype 7 linewidth 1

#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "vdrift_T-300K-normal-model.ps"
set terminal postscript eps enhanced 16
set output "vdrift_T-300K-normal-model.eps"

set line style 1 linetype 7 linewidth 1
set line style 2 linetype 1 linewidth 3

set multiplot
set size 1.0,1.0
set lmargin 8
set rmargin 1


set mytics 2
set mxtics 2

set format y "%2.2f"

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
set origin 0.0,0.0
set size 0.5,0.5
set key bottom
set ytics 0.2

set xlabel "{/Symbol t}\n a)"
set ylabel "u_d({/Symbol t})"

set label 100000 "{/Symbol a} = 100" at 5.5,0.87
set label 1000   "{/Symbol a} = 1.0"  at 3.5,0.75
set label 500    "{/Symbol a} = 0.5" at 5.5,0.70
set label 200    "{/Symbol a} = 0.2" at 5.5,0.52
set label 100    "{/Symbol a} = 0.1" at 3.5,0.45
set label 10     "{/Symbol a} = 0.01" at 5.5,0.42
set label 1      "{/Symbol a} = 0.001" at 5.5,0.18


plot [:8][0:1] \
    "normal-model-N-5e5_T-300_alpha-100_beta-0_speed-0.dat"        using 2:6 notitle with lines 1, \
    "normal-model-N-5e5_T-300_alpha-1_beta-0_speed-0.dat"        using 2:6 notitle with lines 1, \
    "normal-model-N-5e5_T-300_alpha-0.5_beta-0_speed-0.dat"        using 2:6 notitle with lines 1, \
    "normal-model-N-5e5_T-300_alpha-0.2_beta-0_speed-0.dat"        using 2:6 notitle with lines 1, \
    "normal-model-N-5e5_T-300_alpha-0.1_beta-0_speed-0.dat"        using 2:6 notitle with lines 1, \
    "normal-model-N-5e5_T-300_alpha-0.01_beta-0_speed-0.dat"    using 2:6 notitle with lines 1, \
    "normal-model-N-5e5_T-300_alpha-0.001_beta-0_speed-0.dat"    using 2:6 notitle with lines 1


# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+

set origin 0.5,0.0
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top
#set ytics 0.005
set xlabel "{/Symbol a}\n b)"
set ylabel "v_{drift} [m.s^{-1}]"

set nolabel

set nomultiplot

set logscale xy
set grid
set key bottom samplen 1


set xtics ("10^{-4}" 1e-4, "10^{-3}" 1e-3, "10^{-2}" 1e-2, "10^{-1}" 1e-1, "1" 1, "10^1" 10, "10^2" 1e2 )
set ytics ("10^4" 1e4, "10^3" 1e3, "10^{-2}" 1e-2, "10^{-1}" 1e-1, "10^0" 1, "10^1" 10, "10^2" 1e2 )

# fit [1e-4:1e-2] [] f(x) "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 via a; replot
# a               = 0.667449         +/- 0.002494     (0.3736%)
# 
# fit [1e-4:1e-3] [] f(x) "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 via a; replot
# a               = 0.653661         +/- 0.002711     (0.4147%)
# fit [1e-3:1e-2] [] f(x) "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 via a; replot
# a               = 0.668264         +/- 0.002686     (0.4019%)

k = 891.053540967816
a = 0.668264

plot [][1:1e4] \
    "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 title "MC" with points 1, \
    k * sqrt(x) title "~E^{0.50}" with lines 2, \
    k * x**a title "~E^{0.67}" with lines 1


Platon Group <platon@platon.org> http://platon.org/
Copyright © 2002-2006 Platon Group
Site powered by Metafox CMS
Go to Top