function conversion_ad_220621;%conversion_ad
%
close all,

Vcc=5;
% Número de bits
n=5;
niveles=2^n
delta_v=2*Vcc/niveles

Tmax=51;
tiempo=0:.004:Tmax;
ws=2*pi*[1 1.59 2.37 3.61];
amp=[2 1.5 1 0.5];

fsamp=50;
tsamp=0:1/fsamp:Tmax;

 function equis=onda(t);
 %x=amp(1)*sin(ws(1)*tiempo)+amp(2)*sin(ws(2)*tiempo)+amp(3)*sin(ws(3)*tiempo)+amp(4)*sin(ws(4)*tiempo);
 equis=amp(1)*sin(ws(1)*t)+amp(2)*sin(ws(2)*t)+amp(3)*sin(ws(3)*t)+amp(4)*sin(ws(4)*t);
 end% function onda
  
 function discreto=quantice(x);
 discreto=(floor((x+Vcc)/delta_v)*delta_v+delta_v/2)-Vcc;
 end% function quantize(t); 
 
rampa=-Vcc+(2*Vcc/Tmax)*tiempo;
escalera=quantice(rampa);
figure(1),plot(tiempo,rampa,tiempo,escalera) 
 
senial=onda(tiempo);
muestras=onda(tsamp);

figure(2), plot(tiempo,senial,'b',tsamp,muestras,'.k');
   
quantiz=quantice(muestras);
error=muestras-quantiz; 
 figure(3),plot(tiempo,senial,'b',tsamp,quantiz,'r')%,tsamp,muestras,'.k')
 hold on
 plot(tsamp,error,'g')

 std(error)
 delta_v/sqrt(12)
 
figure(4)
hist(error,20)

% FINAL DEL ARCHIVO
end% function conversion_ad_220621;