clear all; close all; clc;

%TX
[y] = audioread('audioIn16k.wav');
S = audioinfo('audioIn16k.wav'); %struct
%[y] = audioread('audioIn32k.wav');
%S = audioinfo('audioIn32k.wav'); %struct

fs = S.SampleRate;
nbits = S.BitsPerSample;

player = audioplayer(y,fs);
%play(player)

figure(1)
%X = linspace(0,2*pi,50)';
%m = 0.3*[cos(X)];
%stem(m)
m = y;
plot(m, 'b')
legend('original','DM')

e(length(m))=0;
e1(length(m))=0;
mq = 0;
delta = 0.04;
numavez = 2;

%delta(length(y)) = 0;
correlacao_delta(length(y))=0;
i = 1;

%for delta = 0.01:0.01:0.2
% correlacao_maxima=max(correlacao_delta);

for n=1:length(m)
    e(n) = m(n) - mq(numavez-1);
    e1(n)=delta*sign(e(n)); %e1 = eq
    mq(n) = e1(n) + mq(numavez-1);
    numavez = n+1;
end

hold on;
stairs(mq, 'k')

for n=1:length(e1)
    if e1(n) >= delta
        ad(n) = 1;
    else
        ad(n) = 0;
    end
end

player2 = audioplayer(mq,fs);
%play(player2)

figure(2)
stairs(ad, 'g')
title 'Saida Digital TX'

%RX
for n=1:length(ad)
    if ad(n) >= 1 %saida do TX
        da(n) = delta; %RX
    else
        da(n) = -delta; %RX
    end
end

figure(3)
stairs(da, 'm')

numavez = 2;
mqrx(length(da))=0;
mqrx(numavez-1) = 0;

for n=1:length(da)
    mqrx(n) = da(n) + mqrx(numavez-1);
    numavez=n+1;
end

figure(4);
plot(mqrx, 'r');
title 'Saida Analógica do RX';