Listato

function [ z ] = quantizzazione( bit )

%QUANTIZZAZIONE Varianza dell'errore di quantizzazione

 

%Inizializzazione

Fs=10000;

x=[0:1/Fs:1-1/Fs];

N=length(x);

y=sin(2*pi*3*x);

 

%Quantizzazione

q=1/2bit;

Qy=quant(y,q);

 

e=Qy-y;

 

%Varianza

v1=q2/12;

v2=sum(e.2)/N-(sum(e)/N)2;

 

%Visualizzazione risultati

plot(x,y);

title('Segnale di ingresso');

pause;

 

plot(x,Qy);

title('Segnale quantizzato');

pause;

 

plot(x,e);

title('Errore');

 

stringa=['Q2/12 = ' num2str(v1)];

disp(stringa);

 

stringa=['sigma2 = ' num2str(v2)];

disp(stringa);

 

pause;

close;