% !TEX encoding = UTF-8 Unicode
% J.Roussel
% Ce document regroupe les codes TIKZ des figures utilisées pour le cours "Notions de probabilités" situé à la page http://femto-physique.fr/physique_statistique/phystat_complement1.php
\documentclass[11pt]{article}
\input{styles_phystat}
\title{Figures TIKZ du cours "notions de probabilité"}
\begin{document}
	
%------- TIKZ : Loi des grands nombres ----------------
\begin{tikzpicture}[scale=0.65]
\begin{axis}[
	title={Tirages de $S=\sum_{i}x_{i}$ avec $N=10$},
	xlabel={tirage n$^\circ$},
	axis x line=bottom,
	axis y line=left,
	xmin=0,
	xmax=100,
	ymin=0,
	ymax=10,
	extra y ticks={5},
	extra y tick labels={$\overline{S}$},]
 	%--- fonctions à tracer ------
 	\addplot+[mark=none] table[x=tirage,y=S]{simu/loiGrandsNombres10.txt};
\end{axis}
\end{tikzpicture}
\begin{tikzpicture}[scale=0.65]
\begin{axis}[
	title={Tirages de $S=\sum_{i}x_{i}$ avec $N=100$},
	xlabel={tirage n$^\circ$},
	axis x line=bottom,
	axis y line=left,
	xmin=0,
	xmax=100,
	ymin=0,
	ymax=100,
	extra y ticks={50},
	extra y tick labels={$\overline{S}$},]
 	%--- fonctions à tracer ------
 	\addplot+[mark=none] table[x=tirage,y=S]{simu/loiGrandsNombres100.txt};
\end{axis}
\end{tikzpicture}
\begin{tikzpicture}[scale=0.65]
\begin{axis}[
	title={Tirages de $S=\sum_{i}x_{i}$ avec $N=1000$},
	xlabel={tirage n$^\circ$},
	axis x line=bottom,
	axis y line=left,
	xmin=0,
	xmax=100,
	ymin=0,
	ymax=1000,
	extra y ticks={500},
	extra y tick labels={$\overline{S}$},]
 	%--- fonctions à tracer ------
 	\addplot+[mark=none] table[x=tirage,y=S]{simu/loiGrandsNombres1000.txt};
\end{axis}
\end{tikzpicture}
%------- FIN FIGURE ------
\begin{tikzpicture}[scale=0.8] 
%------ courbe de Gauss -----
\def \xMoy{2};%moyenne
\def \ecartType {1};%ecart-type
	\begin{axis}[title={$p(x)=\frac{1}{\sigma \sqrt{2\pi}}\mathrm{e}^{-\frac{(x-\overline{x})^{2}}{2\sigma^{2}}}$},xlabel={$x$},axis x line=bottom,axis y
line=left,xmin=-2,xmax=6,ymin=0,ymax=0.5,grid=major]
\addplot+[mark=none,domain=-2:6,samples=200]{(exp(-(x-\xMoy)^2/(2*\ecartType^2)))/(\ecartType*sqrt(6.28))}; 
\addplot+[mark=none,fill=gray!20,draw=black,opacity=0.5,domain=1:3,samples=150]{(exp(-(x-\xMoy)^2/(2*\ecartType^2)))/(\ecartType*sqrt(6.28))}\closedcycle;
	\draw[<->] (axis cs:1,0.1)--(axis cs: 3,0.1) node[midway,above]{$2\sigma$};
	\end{axis}
\end{tikzpicture}
\begin{tikzpicture}[scale=0.8]
	\def \titre {Distribution binomiale de param\`etre $N=20$ et $p=0.4$.}
	\def \xLabel {$S$};
	\def \xMin{0}; 
	\def \xMax{20};
	 \def \yMin{0};
	 \def \yMax{0.25};
	 \begin{axis}[title=\titre,xlabel={\xLabel},extra y ticks={5}, extra y tick labels={$\overline{S}$}, 
	 grid=major,axis x line=bottom,axis y line=left,xmin=\xMin,xmax=\xMax,ymin=\yMin,ymax=\yMax,]
 	%--- fonctions à tracer ------
 	\addplot+[ycomb] table[x=S,y=bernoulli]{simu/loiBernoulli20.txt};
	\addlegendentry{Distribution binomiale};
	\addplot+[mark=none,smooth] table[x=S,y=gauss]{simu/loiBernoulli20.txt};
	\addlegendentry{Approximation gaussienne};
	 \end{axis}
\end{tikzpicture}
\end{document}