Listato

function [ z ] = dithering( in,bit )

%DITHERING Effetto del dithering

 

%Segnale di in

Fs=100;

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

N=length(x);

y(1:N)=in;

 

%Aggiunta rumore

Q=1/2bit;

a=-Q/2;

b=Q/2;

n=a+(b-a)*rand(1,N);

yn=y+n;

 

plot(x,yn);

title('Segnale in ingresso al quantizzatore');

 

%Quantizzazione

Qyn=quant(yn,Q);

 

%Media

Qy=sum(Qyn)/N;

 

%Visualizzazione risultati

stringa=['Valore di ingresso: y = ' num2str(in)];

disp(stringa);

 

stringa=['Media dell''uscita quantizzata: Qy = ' num2str(Qy)];

disp(stringa);

 

pause;

close;