% TP Signaux aléatoires - SRT4
%% TP5 Propriété ergodique d'un signal aléatoire
%
% Nom - Prénom :
%
% On considère le signal aléatoire en temps discret X(n)=\cos(2\pi f n + \phi) où f=1/128 et \phi est une phase aléatoire (variable aléatoire) uniformément répartie sur [0,2\pi[(fonction rand).\\
% Déterminez si X(n) est stationnaire au premier ordre. Calculer la moyenne statistique d'ordre 1 de X(n). La séquence aléatoire X(n) est-elle ergodique relativement à sa moyenne ?
% Générez 4 réalisations de ce signal sur 1024 points. A partir de ces 4 réalisations, estimez la moyenne statistique à chaque instant n (en utilisant la fonction mean), puis déterminez les moyennes temporelles (toujours avec la fonction mean) pour chacune des 4 réalisations. Que constatez-vous ?
% Générez maintenant 40 puis 400 réalisations de 1024 points de ce signal X(n). Dans chaque cas, estimez à nouveau les moyennes statistiques à chaque instant, et calculez les moyennes temporelles. Que constatez-vous ? Est-ce en accord avec la théorie ?
%
%


close all;
clear all;
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  Ergodisme
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%

% E(X(n))=\int_{-\infty}^\infty \cos(2\pi f n + \phi)f_\phi(\phi)d\phi=\frac{1}{2\pi}\int_{0}^{2\pi} \cos(2\pi f n + \phi)d\phi=0.
% E(X(n)) ne dépend pas de n, la séquence aléatoire est donc stationnaire au premier ordre.
% L'ergodisme relativement à la moyenne pour un processus stationnaire se traduit par l'égalité de la moyenne temporelle et de la moyenne statistique. Il est facile de vérifier que la moyenne temporelle d'une réalisation est bien égale à zéro dans notre cas :
% \left\langle X(n)\right\rangle=\lim_{N\rightarrow\infty}\frac{1}{2N+1}\sum_{-N}^N\cos(2\pi f n + \phi)=0
%
%



nbpoints=1024;
freq=1/128;
Points=[0:nbpoints-1];
FPoints=freq*Points;

% génération de  4 réalisations
nbrealisations=4;

% Tracé de 4 réalisations

% Calcul de la moyenne / estimation de la moyenne statistique à chaque
% instant



% Estimation de la moyenne temporelle d'une réalisation


% On remarque que la moyenne temporelle est proche de la valeur théorique zero tandis que l'estimation de la moyenne d'ensemble est loin d'etre égale à zéros pour tout n lorsque nbrealisations=4.
% L'estimation de cette moyenne d'ensemble n'est pas fiable puisqu'elle ne repose pas sur un nombre suffisant de réalisations.
% Plus ne nombre de réalisations est grand, plus la moyenne est proche de zero quel que soit n.