trabajo

clear all;
close all;

f = 1; %Frecuencia de onda seno
fs = 100; %periodo de muestreo de onda seno
t = 0:1/fs:1; %Division de tiempo en segmentos 1/fs

%Ajustes de cambio de fase paralas diferentes señales BBPSK
p1 = 0;
p2 = pi;

%obteniendo el numero de bits a ser modulados
N = input('Ingrese el numero de bits a ser modulados: N = ' );

%generando la señal aleatoriaflujo_bits=round(rand(1,N));

%asignacion de variables dinamicas

time = [];
senal_digital = [];
BPSK = [];
BPSK_N = [];
senal_portadora = [];

% GENERANDO LAS SEÑALES

for ii = 1:1:N% GENERACION DE LA SEÑAL DIGITAL -SEÑAL MODULADORA
if flujo_bits(ii) == 0
bit = zeros(1,length(t));
else
bit = ones(1,length(t));
end
senal_digital = [senal_digital bit];

%GENERANDO LASEÑAL BBPSK
if flujo_bits(ii) == 0
bit = sin(2*pi*f*t+p1);
else
bit = sin(2*pi*f*t+p2);
end
BPSK = [BPSK bit];



%GENERANDO SEÑAL PORTADORA
carrier = sin(2*f*t*pi);
senal_portadora =[senal_portadora carrier];

time = [time t];
t = t + 1;

end

subplot(3,1,1);
plot(time,senal_portadora);
grid on;
title('SEÑAL PORTADORA');
axis tight;

subplot(3,1,2);plot(time,senal_digital,'r');
grid on;
title('SEÑAL DIGITAL ALEATORIA');
axis([0 time(end) -0.5 1.5]);

subplot(3,1,3);
plot(time,BPSK);
grid on;
title('SEÑAL BBPSK');
axis tight;

% Desvanecimientosde la Señal

F = 2; %Khz
dl = 1:1:100; %Distancia de 1 a 100 km, intervalo para desv. lento

figure(2);


L = 32.45 + 20*log10(F) + 20.*log10(dl); %Pérdidas de espacio libre
Pr_Pt= -L; %Relación de potencias
plot(dl, Pr_Pt, 'm')
hold on

%Desvanecimiento lento (Distribución normal)

mu = 0;
sigma = 3;
N = mu + sigma.*randn(1, length(dl)); %Se generan valoresaleatorios con distribución normal de media 0 y sigma 4
L_dl = L + N;
Pr_Pt_dl = -L_dl; %Relación de potencias con desvanecimiento lento
plot(dl, Pr_Pt_dl, 'r')
hold on

%Desvanecimiento rápido...