close all
clc
clear all


%%%%%%%%%% MODULACION ASK  %%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

nbits = input('Digite la longitud de los datos:  ');
muestras = input('Digite el numero de muestras:  ');

data = round (rand(1,nbits));

data_m = muestreo(data,nbits,muestras);

t = 1/muestras : (1/muestras) :length(data);
seno = sin(2*pi*t);
askt = data_m.*seno;

figure(1);
subplot(4,1,1);
plot(data_m(1:100));
title('Datos Digitales');
subplot(4,1,2);
plot(askt(1:100));
title('Se馻l ASK');

b = potencia(askt);
j=1:length(b);
freq=j*1200/128;

figure(2);
subplot(3,2,1);
plot(freq,b),grid on;
title('Potencia ASK');
subplot(3,2,2);
semilogy(freq,b),grid on;
title('Escala Logaritmica');


%%%%%%%%    MODULACION FSK  %%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

seno1 = sin (4*pi*t);
seno2 = sin (2*pi*t);

data_m =  muestreo(data,nbits,muestras);
data_m1 = muestreo1(data,nbits,muestras);

fskt1 = data_m.*seno1;
fskt2 = data_m1.*seno2;

fsk = fskt1 + fskt2;

figure(1);
subplot(4,1,3);
plot(fsk(1:100));
title('Se馻l FSK');

b = potencia(fsk);
j=1:length(b);
freq=j*1200/128;

figure(2);
subplot(3,2,3);
plot(freq,b),grid on;
title('Potencia FSK');
subplot(3,2,4);
semilogy(freq,b),grid on;
title('Escala Logaritmica');

%%%%%%%%    MODULACION PSK  %%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


data_m =  muestreo(data,nbits,muestras);
data_m1 = muestreo1(data,nbits,muestras);

psk1 = data_m.*seno;
psk2 = data_m1.*-seno;

psk = psk1 + psk2;
figure(1);
subplot(4,1,4);
plot(psk(1:100));
title('Se馻l PSK');

b = potencia(psk);
j=1:length(b);
freq=j*1200/128;

figure(2);
subplot(3,2,5);
plot(freq,b);
title('Potencia PSK'),grid on;
subplot(3,2,6);
semilogy(freq,b);
title('Escala Logaritmica'),grid on;