Cuprins
Plotarea datelor 1) rezultate din măsurători sau simulări este probabil cea mai frecventă utilizare a programului gnuplot. Gnuplot are două moduri generale de funcționare - interactiv și script. În modul interactiv, fiecare comandă este scrisa la promptul gnuplot, cu posibilitatea de a folosi la sfârsit comanda save pentru a salva comenzile într-un script. Un script scris manual, sau obtinut cu save poate fi pasat ca parametru programului gnuplot pentru a genera grafice.
Datele discrete de intrare corespund funcțiilor date prin puncte. Avem nevoie de un fișier de date de intrare și de câteva comenzi pentru a manipula datele. Vom începe cu plotarea de bază a datelor simple și apoi vom analiza plotarea datelor cu erori.
La început, vom analiza un fișier de date. Acesta poate fi un fișier text care c onține datele ca coloane. Considerăm fișierul de date numit plotting\_data1.dat ce conține ( gnuplot-ul ignoră liniile care încep cu diez \#):
# comentarii # X Y 1 2 2 3 3 2 4 1
Poate fi plotat scriind:
set style line 1 lc rgb '#0060ad' lt 1 lw 2 pt 7 ps 3.5 plot 'data1' with linespoints ls 1
Aici setăm tipul de punct (pt) și dimensiunea punctului (ps), și culoarea (lc) de utilizat. Pentru stilurile de puncte și linii disponibile, puteți să rulați comanda test la promptul gnuplot. Imaginea rezultată este
Dacă avem colecții de puncte care reprezintă date ne-continue, putem indica aceasta prin introducerea unei linii libere între date
# data2 # XY 1 2 2 3 3 2 4 1
Dacă dorim să folosim altă culoare pentru cel de-al doilea set de date care este totuși în același fișier, putem introduce încă o linie liberă. Apoi trebuie să indexăm blocul de date începând cu comanda 'index'.
# data3 # Primul bloc de date (index 0) # X Y 1 2 2 3 # Al doilea bloc de index (index 1) # X Y 3 2 4 1
Și plotăm cu comenzile:
# albastru: set line style 1 lc rgb '# 0060ad' lt 1 lw 2 pt 7 ps 3.5 # roșu: set line style 2 lc rgb '# dd181f' lt 1 lw 2 pt 5 ps 3.5 plot 'data3' index 0 w lp ls 1 , \ '' index 1 w lp ls 2
Atenție, gnuplot dorește toată comanda pe o singură linie, de aceea poate fi necesar un backslash “\”. După cum putem vedea, am adăugat un alt tip de culoare și punct și am reprezentat cele două seturi de date utilizând indexul și separat parcelele cu o virgulă. Pentru a reutiliza ultimul nume de fișier, putem folosi un nume vid \'\'. Rezultatul este
O necesitate frecventă este reprezentarea datelor cu bare de erori pentru a indica de exemplu comportarea funcției într-un punct. În următorul exemplu avem măsurători ale puterii pentru o rezistență dată, stocate în formatul: r, P, Perror care poate semnifica eroarea de măsurare a puterii:
# X Y stddev 50.000000 0.036990 0.007039 47.000000 0.036990 0.007039 44.000000 0.038360 0.007053 41.000000 0.042160 0.007050 38.000000 0.043200 0.007018 35.000000 0.046900 0.007021 32.000000 0.048840 0.006963 29.000000 0.052000 0.006929 26.000000 0.055470 0.006947 23.000000 0.060000 0.006882 20.000000 0.064660 0.006879 17.000000 0.069600 0.006936 14.000000 0.079800 0.007080 11.000000 0.086920 0.007232 8.000000 0.085500 0.007262 5.000000 0.101260 0.008415 2.000000 0.091000 0.011203 0.000000 0.081480 0.011828
Vom plota astfel:
set xrange [ -2 : 52 ] set yrange [ 0 : 0.12 ] set format y '% .0s' plot 'battery.dat' using 1:2:3 \ w yerrorbars ls 1 , \ '' using 1:2 w lines ls 1
Valorile puterii sunt stocate în Watt în fișierul de date, dar au valori mai mici decât 1. De aceea dorim să folosim mW ca unitate de măsură. Așadar, am setat opțiunea de format pentru a spune gnuplot-ului să folosească “mantisa as base of current logscale”, vezi documentația gnuplot . Apoi, cu indicația using, se specifică gnuplot ce coloane din fișierul de date ar trebui să utilizeze. Deoarece vrem să plotăm eroarea puterii și puterea, avem nevoie de trei coloane - 1,2, și 4. Folosind stilul de plotare yerrorbars nu este posibilă combinarea punctelor cu o linie. Prin urmare, adăugăm o a doua linie la comanda plot pentru a combina punctele cu o linie. Acest lucru va rezulta în:
Putem evita comanda set format din ultimul grafic prin manipularea directă a datelor de intrare:
set yrange [ 0:120 ] plot 'battery.dat' using 1:($2*1000):($3*1000)\ w yerrorbars ls 1
Atenție, sunt necesare paranteze în jurul expresiilor pe coloană, iar numărul coloanei este indicat cu caracterul \2*1000):(2}'</code> afișează coloana a doua a fiecărei linii. De exemplu, pentru acest fișier:
10 2 0.2 11 3 0.3 12 2 0.2 13 3 0.1 14 4 0.05
cat trace.out | awk '{print $1+$2, $2 $3, i++;}'
produce
12 20.2 0 14 30.3 1 14 20.2 2 0 3 16 30.1 4 18 40.05 5
Exercițiul 2: În laboratorul 1, folosiți
$tcp attach [open tcp.tr w] $tcp trace cwnd_ $tcp trace rtt_
pentru a explora relația dintre lungimea cozii la bottleneck (link n2 n3), RTT-ul perceput de TCP, și fracțiunea de debit obținută în concurență cu UDP.
$ x=5 $ echo $(($x / 3)) 1
$ awc() { awk "BEGIN{print $*}"; } funcția poate fi adăugata la .bashrc $ awc "$x / 3" 1.66667
$for i in `seq 10 2 20`; do awc "$i / 10"; done 1 1.2 1.4 1.6 1.8 2 $
$ i=1; while [ $i -le 10 ]; do echo -n "$i "; ns ./script.tcl -rate0 "$(($i*100))Kbps"| grep Throughput |head -n1 | awk '{print $2}'; i=$(($i+1)); done | tee ./out.a 1 505744.0 2 989296.0 3 1495040.0 4 2026480.0 5 2504192.0 6 3004096.0 7 3143088.0 8 3120896.0 9 3118560.0 10 3124400.0 $
echo -n
nu printează un CR la sfârsitul liniei, util pentru a ține un experiment pe o linie ns ./script.tcl
produce un output de tipul Throughput 505744.0 bps
din care reținem cu awk doar câmpul 2head -n1
reține doar prima linie de Throughput tee
lasă outputul pe ecran, dar îl salvează și în fișierul out, gata de plotat Citiți ns2 wireless tutorial
Marc Greis tutorial secțiunea IX .
r -t 0.016905500 -Hs 1 -Hd -2 -Ni 1 -Nx 0.00 -Ny 75.00 -Nz 0.00 -Ne -1.000000 -Nl MAC -Nw --- -Ma 0 -Md 1 -Ms 0 -Mt ACK d -t 1.804824308 -Hs 2 -Hd 2 -Ni 2 -Nx 75.00 -Ny 0.00 -Nz 0.00 -Ne -1.000000 -Nl MAC -Nw COL -Ma 13a -Md 2 -Ms 0 -Mt cbr -Is 0.0 -Id 2.1 -It cbr -Il 1590 -If 0 -Ii 144 -Iv 32 -Pn cbr -Pi 34 -Pf 0 -Po 0
s: Send r: Receive d: Drop f: Forward -t double Time (* For Global Setting) -Ni int Node ID -Nx double Node X Coordinate -Ny double Node Y Coordinate -Nz double Node Z Coordinate -Ne double Node Energy Level -Nl string Network trace Level (AGT, RTR, MAC, etc.) -Nw string Drop Reason -Hs int Hop source node ID -Hd int Hop destination Node ID, -1, -2 -Ma hexadecimal Duration -Ms hexadecimal Source Ethernet Address -Md hexadecimal Destination Ethernet Address -Mt hexadecimal Ethernet Type -P string Packet Type (arp, dsr, imep, tora, etc.) -Pn string Packet Type (cbr, tcp) -Ps sequence number (pentru tcp, coloana 47)
Exercițiul 3: modificați simple-wireless.tcl
din Marc Greis sec IX pentru
$ns_ use-newtrace
)Exercițiu 4 dezactivați RTS/CTS folosind
Mac/802_11 set RTSThreshold_ 3000
și comparați performanța TCP cu cazul precendent. Sugestie: plotați evoluția în timp a numerelor de secvență. Numărul de secvență în timp este de fapt throughput. RTS-CTS introduce un overhead, de fapt timp pierdut, care duce la un throughput redus.
cat simple.tr | grep '^r' | grep AGT | grep tcp | grep -v ack | awk '{print $3, $47}'