%
% Introduction
%==============
%
% 
% Les questions de TD figurent sur ce listing lisez le bien.
%
% L'evalutation sera base sur un compte rendu simple et precis du pour la
% semaine suivante.
% Il y aura un controle continue sur ces TD. Profitez de la presence d'un
% encadrant pour poser des questions PERTINANTES et ainsi avoir une bonne
% note de controle continue.
%
%=========
% partie 1 
%=========
% 
%
%============================================================
% Quantification 
%============================================================
%
clear            % Nettoie l'espace de travail. Taper "help clear"
close all        % Ferme toutes les figures. Taper "help close" 
%
Np=2048;
n=?(0,2*pi,Np);% Creer un vecteur de longueur Np,
cs=cos(2*n);
var_s=sum(cs.^2)/? % calcule la variance.
for B=1:50
  del(B)=2^(-B)*2;
  cs1=round(cs/del(B))*del(B);
  e=?;              % determine l'erreur de quantification.
  var_e1(B)=sum(e.^2)/?;
  var_e2(B)=del(B)^2/12;
  SNR1(B)=10*log10(var_s/var_e1(B));
  SNR2(B)=6.02*B+10.8-20*log10(2/var_s^0.5);
end
plot(log10((var_e1)))
hold on
plot(log10(var_e2),'r')
figure
plot(SNR1)
?
plot(SNR2,'r')
%
% Ajoutez des legendes a vos figures.
%
% Quels commentaires pouvez vous faire sur les courbes ? 
% Cela confirme t-il les formules du cours ?
% 
%========================
%
B=4;
del(B)=2^(-B)*2;
cs1=round(cs/del(B))*del(B);
e=?;
moy1=sum(e)/?
var1=sum(e.^2)/?
[h,x]=hist(e,?);
moy2=sum(h.*x)/?
var2=sum(h.*(x-moy1).^2)/?
figure
plot(?,?)
%Que constatez vous ?