Simulation de la loi Gaussienne: Résolution 4/4 [MATLAB]

La démonstration (en Partie IV) effectuée à la première page du topic sur la Simulation probabiliste, voir la démo sur ce lien : Démonstration

nous a permis d’établir un modèle exploitable pour simuler la loi Gaussienne centrée réduite. 

La Résolution matlab qui en découle est assez simple, voici un nuage de points dont chaque évènement suit une loi gaussienne centrée réduite.

SCRIPT

%% Script Loi Gaussienne
%
clear all;
close all;
clc;
%
disp(‘Loi Gaussienne centrée réduite : Méthode de Box-Muller’);
fprintf(‘\n’);
k=input(‘Entrer le nombre d échantillons k a testé\nk= ‘);

• for i = 1 : k

% Extraction de l’évènement E1 suivant la loi R
E1=sqrt(-2*log(rand(1,1)) );
% Extraction de l’évènement E2 suivant la loi O
E2=2*pi*rand(1,1);
% Evènement X et Y suivant la loi gaussienne centrée réduite
X(i,1)=E1*cos(E2);
Y(i,1)=E1*sin(E2);
hold on;
plot(X(i,1)’, Y(i,1)’,’om’);

end

title([‘Nuage de points d une loi Gaussienne pour k=’,num2str(k)]);
xlabel(‘X’); ylabel(‘Y’);
hold off;