Discussion :

[thtghgh]

Bonjour à tous,

J'utilise matlab afin de modéliser des lois de probabilités.
Je voudrais afficher une courbe sur mon histogramme mais je n'y arrive pas.
Voici l'histogramme

J'aimerais afficher la densité de cette loi, la loi laplace gauss (0,1).

Pouvez vouis m'aider?

Merci

[duf42]

Bonjour,

Pourrais-tu nous montrer ton code? Normalement un simple HOLD ON devrait suffir...

Bonne apm,
Duf

[EDIT] Un exemple vite fait à partir de l'aide:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
figure;
 
ax = axes('NextPlot','add');
 
x = -2.9:0.1:2.9;
y = randn(10000,1);
hist(ax,y,x) 
 
plot(ax, [-2.9:0.1:2.9], 200*abs(sin([-2.9:0.1:2.9])), '-r');

[thtghgh]

voici mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
%% Simulation de la loi Laplace-Gauss
 
clear all
clc
clf
 
n=1000;
 
U=rand(n,1);
V=rand(n,1);
X=sqrt(-log(U)).*cos(2*pi*V);
 
hist(X)

J'aimerais afficher la densité de cette loi qui est f(x)=exp(-x²/2)/sqrt(2*pi) pour x appartenant à R.
Lorque je fais hold on et plot(f) ca ne marche pas..

En fait, le souci c'est que sur ma figure vu que je représente des échantillons les ordonnées vont jusqu'a 250 alors que la densité ne dépasse pas 1/sqrt(2pi) donc en fait il y a un problème d'échelle.
Mais la forme de l'histogramme est la meme que celle de la densité.
L'idéal serait d'avoir un autre axe des ordonnées a droite et que la courbe de la densité se réfère à cet axe.

[phryte]

Bonjour.
Un essai :

Code :

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

1
2
3
4
5
6
hist(X)
x=-3:0.1:3;
f=exp(-x.^2/2)/sqrt(2*pi);
f=max(hist(X))*f/max(f);
hold on
plot(x,f,'r')

[duf42]

Sinon en modifiant un peu mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
figure;
 
ax1 = axes();
ax2 = axes();
 
x = -2.9:0.1:2.9;
y = randn(10000,1);
hist(ax1,y,x) 
 
plot(ax2, [-2.9:0.1:2.9], 200*abs(sin([-2.9:0.1:2.9])), '-r');
 
set(ax2,'XTick',[],'XTickLabel',{''},'YTick',[],'YTickLabel',{''},'Color','none');