> restart:

TP 04

Paramètres de dispersion

et histogrammes

Rappel de consignes

(1) titre en "Times New Roman", corps 36,

NOM Prénom - groupe

TP04 - date

(2) sauvegarder "souvent"

(3) imprimer : paginer (Format/Page_Number)

puis File/Preview.

enfin imprimer en deux colonnes

>

On charge quelques bibliothèques

> with(pldx); # il faut la version 27

> with(simul); # il faut la version 19

>

Paramètres de dispersion

> rra:= rand(8):

répétitions et répétitions : comparer

> rra( ) $ 3 ;

> 'rra( )' $ 3 ;

Création d'une liste de nombres

> n:= 100;

> li:= [ 'rra()' $ n];

> nops(li);

> add(X, X=li)/n ; m:= evalf( %);

Une syntaxe plus stricte, avec affichage avant calculs

> Sum( 'li[k]', k=1..n); value(%);

> stats[describe, mean] (li) ;

> macro(moy= stats[describe, mean] );
moy(li);

> v:= add( (li[k]-m) ^2, k=1..n)/n ;

Koenig

> map( x-> x^2, li) : moy(%) - m^2;

> macro(var= stats[describe, variance] );
var(li); evalf(%);

>

>

Visualisation

> macro(tal= stats[transform, tally] );

> ta:= tal(li);

> macro(histo = xhisto );

> histo(ta, area=1);

> les_interv:= [ seq( k-1/2..k+1/2, k=0..7)] ;

> macro(talto = stats[transform, tallyinto] );

> ta2:= talto(li, les_interv);

> histo(ta2, area=1); pl1:= %:

>

Visualisation générale

> xima:=proc(pl); op(1,pl) : convert(%,list): map(op,%):
map2(op,2,%): max(op(%)): end:

> [m, t, t=0..xima(pl1)]; plot(%, color=red, linestyle=24): pl2:= %:

> display(pl1, pl2);

> [m+sqrt(v), t, t=0..xima(pl1)], [m-sqrt(v), t, t=0..xima(pl1)] ;
plot( {%}, color=red, linestyle=16): pl3:= %:

> display(pl1, pl2, pl3);

>

>