clear all;
clc;
%%%%%%%%%%% Primero simule la señal detectada para generar la señal que se muestreará. En esta simulación, configure la señal detectada para que sea modulación BPSK, fc = 1Mhz
%%%%% % %%%%% La longitud de la señal de banda base es de 100 bits, la tasa de bits de la señal fb = 1 Mbps, la frecuencia de muestreo fs = 100 Msps
signal_base = randint (1,100,2);
signal_base_bpsk = 2 * signal_base-1;
T = 1/1000000;
t0 = 0 : (T / 100) :( 99 * T / 100);
para i = 1: 100,
si signal_base_bpsk (i)> 0
signal_base_pure (((i-1) * 100 + 1): ((i) * 100) ) = sin (2 * pi * 1000000 * t0); de lo
contrario,
signal_base_pure (((i-1) * 100 + 1): ((i) * 100)) = - sin (2 * pi * 1000000 * t0);
end ;
end;
signal_base_awgn = awgn (signal_base_pure, 0);
%%%%%%%%% recopilará los datos bajo fs de 100Msps y realizará una interpolación triple %%%%%%%%%%%%%%
t1 = 0: (T / 100) :( 100 * TT / 100);
s1 = señal_base_awgn;
t2 = 0: (T / 300) :( 100 * TT / 100);
s2 = interp1 (t1, s1, t2) ; %%%% 3 interpolación lineal
%%%%%%%%%% respectivamente para la secuencia de datos s2 después de la interpolación, la secuencia de datos s3 se retrasa 100Ts y la secuencia de datos s4 se retrasa 250Ts.
s3_temp = [zeros (1,100) s2];
s3 = s3_temp (1: 2048);
s4_temp = [zeros (1,250) s2];
s4 = s4_temp (1: 2048);
s6_temp = [zeros (1,120) s2];
s6 = s6_temp (1: 2048);
s7_temp = [zeros (1300) s2];
s7 = s7_temp (1: 2048);
s5 = s2 (1: 2048);
%%%%%%%%%% calcular s5 y s3 por separado , La función de correlación cruzada entre s4
N = 2048;
r35 = ceros (1,500);
para m = 1: 500,
para n = 1: Nm,
r35 (m) = r35 (m) + s5 (n) * s3 ( n + m);
fin;
fin;
r35 = r35 / N;
r45 = ceros (1500);
para m = 1: 500,
para n = 1: Nm,
r45 (m) = r45 (m) + s5 (n) * s4 (n + m);
final;
final;
r45 = r45 / N;
r65 = ceros (1500);
para m = 1: 500,
para n = 1: Nm,
r65 (m) = r65 (m) + s5 (n) * s6 (n + m);
fin;
fin;
r65 = r65 / N;
r75 = ceros (1500);
para m = 1: 500,
para n = 1: Nm,
r75 (m) = r75 (m) + s5 (n) * s7 (n + m);
fin;
fin;
r75 = r75 / N;