Table of Contents

Laboratorul 2

Cuprins

  1. introducere în gnuplot
  2. introducere în awk
  3. comenzi utile în shell
  4. wireless ns-2

1. Introducere în gnuplot

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.

Date simple

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 \#):

data1
 # 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
 # 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
 # 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

Date cu erori

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:

battery.dat
# 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 \$.

În ultimul grafic adăugăm date teoretice și o legendă:
<code gnuplot>
# Legendă
set key at  50,112
set xlabel 'Resistance (ohm)'
set ylabel 'Power (mW)'
# Curba teoretică
P(x) = 1.53**2 * x/(5.67 + x)**2 * 1000
plot 'battery.dat' using 1:($2*1000):($3*1000) \
                    title "Power" w yerr ls 2 , \
      P(x) title 'Theory' w lines ls 1
</code>

În general, legenda este activată cu comanda **set key**, iar poziția sa poate fi specificată prin **set key top left**,
etc. Puteți, de asemenea, să o setați direct la un punct așa cum am
făcut-o aici, pentru a avea suficient spațiu între legendă și
ticuri. Cuvântul cheie din titlu în comanda plot specifică textul
care va fi afișat în legendă.

{{ :isrm:laboratoare:02:gnuplot:battery.png?nolink |}}


  * Alte comenzi utile gnuplot:
<code gnuplot>
set term png       # setează terminalul într-un format
set out 'file.png' # va plota într-un fisier
set xrange [0:4]   # domeniul funcției
set ytics 0.2      # axa y etichetată la 0.2
test               # afișează puncte/linii cu stilul default
replot             # redesenează
help plot with     # manual la prompt
save 'file.plot'   # creează un script pt plotul curent
</code>
 


 **Exercițiul 1** realizați folosind ''gnuplot'' graficul din laboratorul 1. Indicați semnificațiile axelor și legenda. Salvați imaginea în format ''png/svg'', și scriptul pentru restaurarea imaginilor.
       * Exemplu<code gnuplot>
plot 'out.tr' using 1:2 t "TCP" with lp, \
     'out.tr' using 1:3 t "UDP" with lines lw 3
set xlabel 'Time[s]'
set ylabel 'Bandwidth[Mbps]'
set grid ytics
set term png
set out 'lab02bw.png' 
replot
save 'lab02bw.plot'  # se va putea reface graficul cu gnuplot ./lab02bw.plot
</code>  {{:isrm:laboratoare:02:lab02bw.png?300|}}

=== 2. Introducere în awk ===
    * AWK (K vine de la Kernighan) este mic, simplu, și rapid, spre deosebire de perl sau python. Nu poți face tot ce faci în perl/python, dar poți face foarte ușor multe taskuri de procesare de text. Are o sintaxă apropiată de C, dar preferă datele organizate pe coloane, ca foarte multe date în rețelistică: trace-uri de simulare, tcpdump, loguri, etc. Un mare avantaj este ca poate fi rulat direct de pe linia de comandă, fără a mai folosi un script separat - de multe ori apare într-un pipeline cu cat, sed, tr.      
    * În cazul cel mai des întâlnit, se specifică un program care este rulat succesiv pentru fiecare linie de intrare: <code awk>cat trace.out | awk '{print $2}'</code> afișează coloana a doua a fiecărei linii. De exemplu, pentru acest fișier:

trace.out
10 2    0.2
11 3 0.3
12 2 0.2
 
13 3 0.1
14 4 0.05 

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.

3. comenzi utile de shell

4. wireless ns-2

Citiți ns2 wireless tutorial

Marc Greis tutorial secțiunea IX .

Exercițiul 3: modificați simple-wireless.tcl din Marc Greis sec IX pentru

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.

1) Ghid bazat pe Gnuplotting
2) cat simple.tr | grep '^r.*AGT.*ack.*tcp' | awk '{print \$3, \$47}'| head
3) cat simple.tr | grep '^d' | grep MAC | grep 'tcp' | grep -v 'ack' | wc -l
4) cat simple.tr | grep '^d' | grep IFQ | grep 'tcp' | grep -v 'ack' | wc -l
5) atenție, avem cadre ACK, tcp fără ack , tcp și ack. ACK sunt confirmări 802.11, ack sunt confirmări TCP. Se pierd 11 cadre TCP/date (linii cu tcp, dar fără ack) la nivelul 2, și 83 sunt aruncate din coadă. Dacă destinația nu mai răspunde în aer, 802.11 ajunge la maximum retries, apoi trece la următorul pachet. TCP însă nu trece mai departe după pierdere. Când destinația iese din zona de comunicare, TCP continuă să trimită deoarece are fereastră suficientă, și se pierde din coadă.
6) un ack se generează ca răspuns la un tcp. Dacă tcp sunt pierdute la transmisie (coadă sau aer), se generează mai puține ack
7) ACK=0 pierderi, RTS=14 pierderi. Dacă nodurile sunt în apropiere, conversația RTS-CTS-Date-ACK se desfășoară cu bine, nu se pierd ACK-uri. La distanță mare, destinația nu răspunde la RTS, deci sunt pierdute, și nu se mai ajunge la ACK.