INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA
MECÁNICA Y ELÉCTRICA
UNIDAD CULHUACAN
Ingeniería en Computación
Asignatura: Modulación Digital
Profesor: Dr. Clara Cruz Ramos
Práctica No.10: Modulaciones digitales “16PSK y 32PSK”
Equipo No. 10
Alumnos:
Camacho Suárez Agustín
Del Ángel Ruíz Jesús Roberto
Grupo: 6CM11
Fecha de Entrega: 14 de junio del 2018
Objetivos. –
Objetivo general: Familiarizarse con las técnicas de modulación por
desplazamiento de fase 16PSK y 32PSK.
Objetivos específicos:
Obtener las formas de onda obtenidas mediante las técnicas de modulación
16PSK y 32PSK.
Analizar las formas de ondas, resultados de las modulaciones.
Comprender los principios de las modulaciones 16PSK y 32PSK.
Desarrollo. –
Compilar los programas de las modulaciones 16PSK y 32PSK en diferentes
archivos.
Detallar las características de las formas de ondas obtenidas.
El programa debe ser diseñado de tal forma que se pueda cambiar la señal de
información por el usuario y que al final la mande a imprimir, la señal portadora y
la señal modulada, ya sea 16PSK u 32PSK.
Programa (modulación 16PSK)
clear all;
close all;
str2=('Digite la informacion binaria, usando [] y espacio entre bits');
disp(str2)
bits=input('Ingresar el valor: ');
fprintf('\n\n\n');
%bits=[0,0,1,0,0,1,1,1,0,1,1,1];
%Desfase 16PSK
P1=(0.0625*pi); %11.25 Grados
P2=(0.1875.*pi);%33.75 Grados
P3=(0.3125*pi); %56.25 Grados
P4=(0.4375*pi); %78.75 Grados
P5=(0.5625*pi); %101.25 Grados
P6=(0.6875.*pi);%123.75 Grados
P7=(0.8125*pi); %146.25 Grados
P8=(0.9375*pi); %168.75 Grados
P9=(1.0625*pi); %191.25 Grados
P10=(1.1875.*pi);%213.75 Grados
P11=(1.3125*pi); %236.25 Grados
P12=(1.4375*pi); %258.75 Grados
P13=((281.25*pi)/180); %281.25 Grados
P14=(1.6875.*pi);%303.75 Grados
P15=(1.8125*pi); %326.25 Grados
P16=(1.9375*pi); %348.75 Grados
%Definicion de Frecuencias y periodos
f=1;
fs=100;
acs=5; %amplitud
t=0:1/fs:1; % Periodo para la señal digital y portadora
t1=0:1/fs:2; % Periodo para la señal 8PSK
%Definicion de variable a usar
tiempo=[];
tiempo1=[];
Digital=[];
Portadora=[];
M16PSK=[];
% Secuencia de datos
for i=1:1:length(bits)
if bits(i)==0
z=zeros(1,length(t));
Digital = [Digital z];
end
if bits(i)==1
o=ones(1,length(t));
Digital = [Digital o];
end
% Portadora
Portadora=[Portadora acs*(sin(2*pi*f*t))];
tiempo=[tiempo t];
t= t+1;
end
%Señal 16PSK
for iiii=1:4: length(bits)
% Caso 1 = 0000 Desfase de 11.25
if bits(iiii)==0 && bits(iiii+1)==0 && bits(iiii+2)==0 &&
bits(iiii+3)==0
bits0000= acs*sin(2*pi*f*t1 + P1);
M16PSK=[M16PSK (bits0000)];
end
% Caso 2 = 0001 Desfase de 33.75
if bits(iiii)==0 && bits(iiii+1)==0 && bits(iiii+2)==0 &&
bits(iiii+3)==1
bits0001= acs*sin(2*pi*f*t1 + P2);
M16PSK=[M16PSK (bits0001)];
end
% Caso 3 = 0010 Desfase de 56.25
if bits(iiii)==0 && bits(iiii+1)==0 && bits(iiii+2)==1 &&
bits(iiii+3)==0
bits0010= acs*sin(2*pi*f*t1 + P3);
M16PSK=[M16PSK (bits0010)];
end
% Caso 4 = 0011 Desfase de 78.75
if bits(iiii)==0 && bits(iiii+1)==0 && bits(iiii+2)==1 &&
bits(iiii+3)==1
bits0011= acs*sin(2*pi*f*t1 + P4);
M16PSK=[M16PSK (bits0011)];
end
% Caso 5 = 0100 Desfase de 101.25
if bits(iiii)==0 && bits(iiii+1)==1 && bits(iiii+2)==0 &&
bits(iiii+3)==0
bits0100= acs*sin(2*pi*f*t1 + P5);
M16PSK=[M16PSK (bits0100)];
end
% Caso 6 = 0101 Desfase de 123.75
if bits(iiii)==0 && bits(iiii+1)==1 && bits(iiii+2)==0 &&
bits(iiii+3)==1
bits0101= acs*sin(2*pi*f*t1 + P6);
M16PSK=[M16PSK (bits0101)];
end
% Caso 7 = 0110 Desfase de 146.25
if bits(iiii)==0 && bits(iiii+1)==1 && bits(iiii+2)==1 &&
bits(iiii+3)==0
bits0110= acs*sin(2*pi*f*t1 + P7);
M16PSK=[M16PSK (bits0110)];
end
% Caso 8 = 0111 Desfase de 168.75
if bits(iiii)==0 && bits(iiii+1)==1 && bits(iiii+2)==1 &&
bits(iiii+3)==1
bits0111= acs*sin(2*pi*f*t1 + P8);
M16PSK=[M16PSK (bits0111)];
end
% Caso 9 = 1000 Desfase de 191.25
if bits(iiii)==1 && bits(iiii+1)==0 && bits(iiii+2)==0 &&
bits(iiii+3)==0
bits1000= acs*sin(2*pi*f*t1 + P9);
M16PSK=[M16PSK (bits1000)];
end
% Caso 10 = 1001 Desfase de 213.75
if bits(iiii)==1 && bits(iiii+1)==0 && bits(iiii+2)==0 &&
bits(iiii+3)==1
bits1001= acs*sin(2*pi*f*t1 + P10);
M16PSK=[M16PSK (bits1001)];
end
% Caso 11 = 1010 Desfase de 236.25
if bits(iiii)==1 && bits(iiii+1)==0 && bits(iiii+2)==1 &&
bits(iiii+3)==0
bits1010= acs*sin(2*pi*f*t1 + P11);
M16PSK=[M16PSK (bits1010)];
end
% Caso 12 = 1011 Desfase de 258.75
if bits(iiii)==1 && bits(iiii+1)==0 && bits(iiii+2)==1 &&
bits(iiii+3)==1
bits1011= acs*sin(2*pi*f*t1 + P12);
M16PSK=[M16PSK (bits1011)];
end
% Caso 13 = 1100 Desfase de 281.25
if bits(iiii)==1 && bits(iiii+1)==1 && bits(iiii+2)==0 &&
bits(iiii+3)==0
bits1100= acs*sin(2*pi*f*t1 + P13);
M16PSK=[M16PSK (bits1100)];
end
% Caso 14 = 1101 Desfase de 303.75
if bits(iiii)==1 && bits(iiii+1)==1 && bits(iiii+2)==0 &&
bits(iiii+3)==1
bits1101= acs*sin(2*pi*f*t1 + P14);
M16PSK=[M16PSK (bits1101)];
end
% Caso 15 = 1110 Desfase de 326.25
if bits(iiii)==1 && bits(iiii+1)==1 && bits(iiii+2)==1 &&
bits(iiii+3)==0
bits1110= acs*sin(2*pi*f*t1 + P15);
M16PSK=[M16PSK (bits1110)];
end
% Caso 16 = 1111 Desfase de 348.75
if bits(iiii)==1 && bits(iiii+1)==1 && bits(iiii+2)==1 &&
bits(iiii+3)==1
bits1111= acs*sin(2*pi*f*t1 + P7);
M16PSK=[M16PSK (bits1111)];
end
tiempo1= [tiempo1 t1];
t1=t1+2;
end
% Representacion de las señales
% Señal Digital
subplot(3,1,1);
plot(tiempo, Digital,'lineWidth',2.5);
title('Señal Digital');
axis([0 16 -0.5 1.5]);
xlabel('Tiempo [S]');
ylabel('Amplitud [V]');
grid on;
%Señal Portadora
subplot(3,1,2);
plot(tiempo,Portadora,'lineWidth',2.5);
title('Señal Portadora');
xlabel('Tiempo [S]');
ylabel('Amplitud [V]');
grid on;
% Señal 16PSK
subplot(3,1,3);
plot(tiempo1,M16PSK,'lineWidth',2.5);
title('Señal 16PSK');
xlabel('Tiempo [S]');
ylabel('Amplitud [V]');
grid on;
Programa (modulación 32PSK)
clear all;
close all;
str2=('Digite la informacion binaria, usando [] y espacio entre bits');
disp(str2)
bits=input('Ingresar el valor: ');
fprintf('\n\n\n');
%bits=[0,0,1,0,0,1,1,1,0,1,1,1];
%Desfase 32PSK
P1=(0.03125*pi); %5.625 Grados
P2=(0.09375*pi);%16.875 Grados
P3=(0.15625*pi); %28.125 Grados
P4=(0.21875*pi); %39.375 Grados
P5=(0.28125*pi); %50.625 Grados
P6=(0.34375*pi);%61.875 Grados
P7=(0.40625*pi); %73.125 Grados
P8=(0.46875*pi); %84.375 Grados
P9=(0.53125*pi); %95.625 Grados
P10=(0.59375*pi);%106.875 Grados
P11=(0.65625*pi); %118.125 Grados
P12=(0.71875*pi); %129.375 Grados
P13=(0.78125*pi); %140.625 Grados
P14=(0.84375*pi);%151.875 Grados
P15=(0.90625*pi); %163.125 Grados
P16=(0.96875*pi); %174.375 Grados
P17=(1.03125*pi); %185.625 Grados
P18=(1.09375*pi);%196.875 Grados
P19=(1.15625*pi); %208.125 Grados
P20=(1.21875*pi); %219.375 Grados
P21=(1.28125*pi); %230.625 Grados
P22=(1.34375*pi);%241.875 Grados
P23=(1.40625*pi); %253.125 Grados
P24=(1.46875*pi); %264.375 Grados
P25=(1.53125*pi); %275.625 Grados
P26=(1.59375*pi);%286.875 Grados
P27=(1.65625*pi); %298.125 Grados
P28=(1.71825*pi); %309.375 Grados
P29=(1.78175*pi); %320.625 Grados
P30=(1.84375*pi);%331.875 Grados
P31=(1.90625*pi); %343.125 Grados
P32=(1.96875*pi); %354.375 Grados
%Definicion de Frecuencias y periodos
f=1;
fs=100;
acs=5; %amplitud
t=0:1/fs:1; % Periodo para la señal digital y portadora
t1=0:1/fs:2; % Periodo para la señal 8PSK
%Definicion de variable a usar
tiempo=[];
tiempo1=[];
Digital=[];
Portadora=[];
M32PSK=[];
% Secuencia de datos
for i=1:1:length(bits)
if bits(i)==0
z=zeros(1,length(t));
Digital = [Digital z];
end
if bits(i)==1
o=ones(1,length(t));
Digital = [Digital o];
end
% Portadora
Portadora=[Portadora acs*(sin(2*pi*f*t))];
tiempo=[tiempo t];
t= t+1;
end
%Señal 32PSK
for iiiii=1:5: length(bits)
% Caso 1 = 00000 Desfase de 5.625
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits00000= acs*sin(2*pi*f*t1 + P1);
M32PSK=[M32PSK (bits00000)];
end
% Caso 2 = 00001 Desfase de 16.875
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits00001= acs*sin(2*pi*f*t1 + P2);
M32PSK=[M32PSK (bits00001)];
end
% Caso 3 = 00010 Desfase de 28.125
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits00010= acs*sin(2*pi*f*t1 + P3);
M32PSK=[M32PSK (bits00010)];
end
% Caso 4 = 00011 Desfase de 39.375
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits00011= acs*sin(2*pi*f*t1 + P4);
M32PSK=[M32PSK (bits00011)];
end
% Caso 5 = 00100 Desfase de 50.625
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits00100= acs*sin(2*pi*f*t1 + P5);
M32PSK=[M32PSK (bits00100)];
end
% Caso 6 = 00101 Desfase de 61.875
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits00101= acs*sin(2*pi*f*t1 + P6);
M32PSK=[M32PSK (bits00101)];
end
% Caso 7 = 00110 Desfase de 73.125
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits00110= acs*sin(2*pi*f*t1 + P7);
M32PSK=[M32PSK (bits00110)];
end
% Caso 8 = 00111 Desfase de 84.375
if bits(iiiii)==0 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits00111= acs*sin(2*pi*f*t1 + P8);
M32PSK=[M32PSK (bits00111)];
end
% Caso 9 = 01000 Desfase de 95.625
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits01000= acs*sin(2*pi*f*t1 + P9);
M32PSK=[M32PSK (bits01000)];
end
% Caso 10 = 01001 Desfase de 106.875
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits00000= acs*sin(2*pi*f*t1 + P10);
M32PSK=[M32PSK (bits01001)];
end
% Caso 11= 01010 Desfase de 118.125
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits01010= acs*sin(2*pi*f*t1 + P11);
M32PSK=[M32PSK (bits01010)];
end
% Caso 12 = 01011 Desfase de 129.375
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits01011= acs*sin(2*pi*f*t1 + P12);
M32PSK=[M32PSK (bits01011)];
end
% Caso 13 = 01100 Desfase de 140.625
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits01100= acs*sin(2*pi*f*t1 + P13);
M32PSK=[M32PSK (bits01100)];
end
% Caso 14 = 01101 Desfase de 151.875
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits01101= acs*sin(2*pi*f*t1 + P14);
M32PSK=[M32PSK (bits01101)];
end
% Caso 15 = 01110 Desfase de 163.125
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits01110= acs*sin(2*pi*f*t1 + P15);
M32PSK=[M32PSK (bits01110)];
end
% Caso 16 = 01111 Desfase de 174.375
if bits(iiiii)==0 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits01111= acs*sin(2*pi*f*t1 + P16);
M32PSK=[M32PSK (bits01111)];
end
% Caso 17 = 10000 Desfase de 185.625
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits10000= acs*sin(2*pi*f*t1 + P17);
M32PSK=[M32PSK (bits10000)];
end
% Caso 18 = 10001 Desfase de 196.875
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits10001= acs*sin(2*pi*f*t1 + P18);
M32PSK=[M32PSK (bits10001)];
end
% Caso 19 = 10010 Desfase de 208.125
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits10010= acs*sin(2*pi*f*t1 + P19);
M32PSK=[M32PSK (bits10010)];
end
% Caso 20 = 10011 Desfase de 219.375
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits10011= acs*sin(2*pi*f*t1 + P20);
M32PSK=[M32PSK (bits10011)];
end
% Caso 21 = 10100 Desfase de 230.625
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits10100= acs*sin(2*pi*f*t1 + P21);
M32PSK=[M32PSK (bits10100)];
end
% Caso 22 = 10101 Desfase de 241.875
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits10101= acs*sin(2*pi*f*t1 + P22);
M32PSK=[M32PSK (bits10101)];
end
% Caso 23 = 10110 Desfase de 253.125
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits10110= acs*sin(2*pi*f*t1 + P23);
M32PSK=[M32PSK (bits10110)];
end
% Caso 24 = 10111 Desfase de 264.375
if bits(iiiii)==1 && bits(iiiii+1)==0 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits10111= acs*sin(2*pi*f*t1 + P24);
M32PSK=[M32PSK (bits10111)];
end
% Caso 25 = 11000 Desfase de 275.625
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits11000= acs*sin(2*pi*f*t1 + P25);
M32PSK=[M32PSK (bits11000)];
end
% Caso 26 = 11001 Desfase de 286.875
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits11001= acs*sin(2*pi*f*t1 + P26);
M32PSK=[M32PSK (bits11001)];
end
% Caso 27 = 11010 Desfase de 298.125
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits11010= acs*sin(2*pi*f*t1 + P27);
M32PSK=[M32PSK (bits11010)];
end
% Caso 28 = 11011 Desfase de 309.375
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==0 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits11011= acs*sin(2*pi*f*t1 + P28);
M32PSK=[M32PSK (bits11011)];
end
% Caso 29 = 11100 Desfase de 320.625
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==0
bits11100= acs*sin(2*pi*f*t1 + P29);
M32PSK=[M32PSK (bits11100)];
end
% Caso 30 = 11101 Desfase de 331.875
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==0 && bits(iiiii+4)==1
bits11101= acs*sin(2*pi*f*t1 + P30);
M32PSK=[M32PSK (bits11101)];
end
% Caso 31 = 11110 Desfase de 343.125
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==0
bits11110= acs*sin(2*pi*f*t1 + P31);
M32PSK=[M32PSK (bits11110)];
end
% Caso 32 = 11111 Desfase de 354.375
if bits(iiiii)==1 && bits(iiiii+1)==1 && bits(iiiii+2)==1 &&
bits(iiiii+3)==1 && bits(iiiii+4)==1
bits11111= acs*sin(2*pi*f*t1 + P32);
M32PSK=[M32PSK (bits11111)];
end
tiempo1= [tiempo1 t1];
t1=t1+2;
end
% Representacion de las señales
% Señal Digital
subplot(3,1,1);
plot(tiempo, Digital,'lineWidth',2.5);
title('Señal Digital');
axis([0 20 -0.5 1.5]);
xlabel('Tiempo [S]');
ylabel('Amplitud [V]');
grid on;
%Señal Portadora
subplot(3,1,2);
plot(tiempo,Portadora,'lineWidth',2.5);
title('Señal Portadora');
xlabel('Tiempo [S]');
ylabel('Amplitud [V]');
grid on;
% Señal 32PSK
subplot(3,1,3);
plot(tiempo1,M32PSK,'lineWidth',2.5);
title('Señal 32PSK');
xlabel('Tiempo [S]');
ylabel('Amplitud [V]');
grid on;
Pruebas y Resultados. -
Modulación 16PSK
Datos de código:
Amplitud “acs” = 5 [V]
Frecuencia de la portadora= 1 [Hz]
Frecuencia de muestreo= 100 [Hz]
Según con lo obtenido en la gráfica, 16PSK es otra forma de modulación angular
de amplitud constante, es decir, que con esta modulación son posibles dieciséis
fases de salida diferentes.
Se ve que, conforme a la señal de información, cada que sean cuatro bits, habrá
cambio de fase en los símbolos de la señal modulada.
Se propone una frecuencia de portadora pequeña, para que se pueda observar el
cambio de fase, ya que i fuera mas grande, no se podría ver dicho efecto.
Modulación 32PSK
Datos de código:
Amplitud “acs” = 5 [V]
Frecuencia de la portadora= 1 [Hz]
Frecuencia de muestreo= 100 [Hz]
Según los resultados obtenidos en la gráfica, un PSK de treinta y dos fases es una
técnica de codificación en donde M= 32, es decir, que en el modulador existen
treinta y dos posibles fases de salida.
Se ve que, conforme a la señal de información, cada que sean cinco bits, habrá
cambio de fase en los símbolos de la señal modulada.
Se propone una frecuencia de portadora pequeña, para que se pueda observar el
cambio de fase, ya que i fuera más grande, no se podría ver dicho efecto.
Conclusiones. –
