Matlab - Filtro de reconstrucción de señal muestreada/cuantificada

%------------------------------------

A = 1; %amplitud

F = 2;  %frecuencia

duracion=1; %duracion de representacion

desfase=0;

t = 0:0.001:duracion;

%señal periodica

x=A*sin(2*pi*F*t+desfase);

subplot(2,1,1),plot(t,x,'g');

axis([0 duracion -2*A 2*A]);

hold on

%muestreo

Fs= 40;  %%MIN DOS VECES MAYOR QUE LA FRECUENCIA MAX DE LA SEÑAL

Tm= 0:1/Fs:duracion;

h= pulstran(t,Tm,'rectpuls',0.01); %funcion de pulsos unitarios

sampling= h.*x;

subplot(2,1,2),plot(t, sampling);

axis([0 duracion -2*A 2*A]);

%-----------------------------------

%------ cuantificacion uniforme---------------

%Partiendo de la señal x de seno del apartado anterior:

[I,J]=max(x); %obtiene en I la amplitud maxima

num_niveles=6; %numero de escalones al cuantificar

delta=2*I/num_niveles; %Divide el rango de la señal en escalones uniformes

%partition divide el rango de la seÒal en L regiones, porlo que tendra l-1 elementos(los diferentes escalone)

partition=[-1+delta:delta:1-delta];

%asigna un valor a cada region definica con partition(punto medio de cada escalon)

codebook=[-1+(delta/2):delta:1-(delta/2)];

[index,quants,distor]=quantiz(x, partition, codebook);

figure,subplot(2,1,1),plot(t,x,'r')

%dibujamos señal original

hold on

plot(t, quants)

%dibujamos los k, correspondiente a lacuantificacion

title('Señal cuantificada');

hold off

%------------------------------------------------------