Discussion :

[AC88mm]

Bonsoir à tous, j'ai un problème que je n'arrive pas à résoudre.

J'ai une colonne de valeurs, appelée C, du type:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
C'= 1 2 4 2 5 1 3 6 4 5 5 3

Je veux pouvoir tracer une courbe qui serait une espèce d'histogramme de cette colonne.

Mettons que j'ai y fois le nombre x. Je voudrais pouvoir placer le point de coordonnées x en abscisse et y en ordonnée (nombre d'occurrences de la valeur x), pour avoir la courbe y=f(x).

Avec un histogramme, donc, c'est faisable rapidement (moyennant quelques ajouts à hist pour rendre la chose exploitable avec les bandes), mais un plot

J'ai longuement réfléchi au problème, pour stocker le nombre d'occurrences des différentes valeurs etc., mais je dois dire que là, je ne sais du tout comment faire
Et je n'ai pas trouvé de commande du genre histplot

Pourriez vous m'aider? Merci d'avance!

bonus: et une courbe logarithmique en log(nombre d'occurrences de x)=f(log(valeur x)), y a t-il moyen de le faire?

[rostomus]

Salut,

Essaie:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
x=round(10*rand(1,100));
[n,out]=hist(x,unique(x));
plot(out,n)

pour l'echelle logarithmique, il y a la fonction loglog:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

loglog(out,n)

Bonus: pour les semilogarithmique: semilogx et semilogy.

[AC88mm]

Superbe, ça marche! Merci, merci, et merci!

Et allez, j'ose le demander:

Je tombe donc sur une certaine courbe qui a l'apparence d'une droite. Comment faire pour tracer la droite qui approche au mieux la courbe? Je pense qu'il doit y avoir une commande puissante toute faite, mais je ne la connais pas. Merci encore?

[rostomus]

Tu peux utiliser polyfit.

[EDIT]
ou bien la division a gauche, exemple:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
 
x=1:5;
y=5*x+2+randn(1,5);
c=[x' ones(5,1)]\y';
plot(x,y,'*')
hold all
yy=c(1).*x+c(2);
plot(x,yy)

[LordPeterPan2]

Superbe, ça marche! Merci, merci, et merci!

Et allez, j'ose le demander:

Je tombe donc sur une certaine courbe qui a l'apparence d'une droite. Comment faire pour tracer la droite qui approche au mieux la courbe? Je pense qu'il doit y avoir une commande puissante toute faite, mais je ne la connais pas. Merci encore?

Si x et y sont des vecteurs de taille 1xn et que y=a*x+b+e avec a et b constant et e alors on trouve a et b minimisant e en faisant le code suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
A=[x(:) ones(size(x(:)))];
B=y(:);
X=A\B
a=X(1)
b=X(2)