%Scrip que simula el código miller.

%Por: Diego Barragán Guerrero

%diegokillemall@yahoo.com

%h-->vector de bits a codificar

%est_inicial-->estado inicial de la señal (puede ser 1 ó -1)

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

h=[1 1 0 0 1 0 1 1 0 0 0 1 1 1 1 0];%Señal

est_inicial=-1;%Estado inicial

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

clc;%Limpia el command window

close all;%Cierra todas las figuras.

con=est_inicial;%Puede ser 1 ó -1

long=length(h);%Número de bits de la señal

n=1;%Estado inicial para el ciclo "while"

ac=[];%Matriz vacía que contendrá la señal codificada.

bits=[];%Matriz vacía que contendrá la señal original.

h(long+1)=0;%Valor de extensión de la señal

while n<=long%Codifica hasta terminar los bits de la señal.

    if h(n)==1 %Si el bit es 1

        bit=[ones(1,100)];

        s=[con*ones(1,50) -con*ones(1,50)];

        con=con*-1;%Conmuta el estado de la señal.

    else %Si el bit es 0

        bit=[zeros(1,100)];

        s=[con*ones(1,100)];

        if h(n+1)==0%Si el siguiente bit es 0

            con=con*-1;%Conmuta el estado de la señal

        end

    end

    ac=[ac s];%Acumulador de la señal miller.

    bits=[bits bit];%Acumulador de la señal original.

    n=n+1;%Incremento del contador del loop

    s=[];%Reset de la matriz temporal s.

end

subplot(2,1,1);plot(bits,'LineWidth',2);

title('Señal de entrada');

set(gca,'xtick',0:100:100*long)%Grid solo para espacios de 100

axis([0 100*(length(h)-1) -2 2])%Seteo de los ejes

grid on %Activa el grid

%*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*

subplot(2,1,2);plot(ac,'LineWidth',2)%Grafica la señal.

title('Señal Codificada MILLER')

set(gca,'xtick',0:100:100*long)%Grid solo para espacios de 100

axis([0 100*(length(h)-1) -2 2])%Seteo de los ejes

grid on %Activa el grid