Discussion :

[ninikosim]

Bonjour,
j'ai un petit souci en fait je voudrai estimer une densité d'un vecteur de loi normal en utilisant le noyau gaussien j'ai utiliser se programme mais cela ne marche pas

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 
 
% // nc est le nombre de classes dans l'histogramme 
% // n la taille de l'échantillon
% // h le pas de l'ordre de n^(-1/5)
 
function []= estim_gaussien(nc,n,hn)
  nc=100;
  n=1000;
  hn=1/(1000^(1/5));
 
 % // gaussiennes centrées réduites de variance 1
 
 x=randn(n,1);
 
  hist(x,nc) 
 
%   // estimation de la densité de la loi de E par la méthode des noyaux,
 
  C=[min(x)-1:1/n:max(x)+1];
   %//tracé de la vraie densité
  f =  exp(-C.^2/2)/sqrt(2*pi);
  plot(C,f)
 
  for i=1:length(C)
    B(i)=1/(n*hn)*sum((exp(((C(i)-x)/hn).^2)/2)*(1/sqrt(2*pi)));
  end 
  hold on 
  plot(C,B,'r')
 
 
end

merci pour votre aide

[duf42]

Bonjour,

Pourrais-tu nous dire ce qui ne marche pas? Erreur MATLAB (laquelle), mauvais résultat (qu'attends-tu)?

Bonne journée,
Duf

[Jerome Briot]

Il faut modifier le début de ton code pour gérer les argument optionnels :

Code :

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

1
2
3
4
5
6
7
function estim_gaussien(nc,n,hn)

if nargin == 0
  nc=100;
  n=1000;
  hn=1/(1000^(1/5));
end

NARGIN donne le nombre d'argument passé en entrée de ta fonction.

Ou alors supprimer la déclaration de ces trois variables dans la fonction ESTIM_GAUSSIEN et les mettre dans l'appel de la fonction (par exemple dans le Command Window) :

Code :

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

nc=100; n=1000; hn=1/(1000^(1/5)); estim_gaussien(nc,n,hn)

[sparfait]

Si tu as la toolbox "Statistics", je crois que la fonction ksdensity le fait (Bon d'accord il faut avoir la toolbox).

[ninikosim]

bonjour,

oui j'ai la toolbox statics,
j'ai essayé avec la fonction ksdensity mais je n'ai pas compris comment cela marche car elle me renvoi l'allure de la densité. et moi se que je veut c'est l'approximation de le densité par le noyau gaussien.