modelo_inti.m

Modelo de Jakes y Rayleigh

%Modelo Jacks con Rayleigh

clc;
clear all;
close all;

k = 12;% Rayos del "Scattering"
fd = 0.02;% N鷐ero m醲imo para el Doppler Shifting
M = 50;% Senoidales
sumatoria = 0;% Sumatioria
N = 1024;% Tama馻no de la secuencia
ds = 0.1;% Duraci髇
T = (N-1)*ds; 
t = 0:ds:T;

for n=1:M
betan = (pi*n)/(M+1);
thetan = betan + (2*pi*(k-1)/(M+1));
alphan = pi * (n-0.5)/(2*M);
fn = fd * cos(betan);
sumatoria = sumatoria + (cos(betan) + j*sin(betan)) * cos(2*pi*fn*t+thetan);    
end

alpha = 0;
R = 2*sqrt(2)*(sumatoria+ 1/sqrt(2)*(cos(alpha)+ j*sin(alpha)) * cos(2*pi*fd*t));

[corr,m] = xcorr(R,'coeff'); 
m1 = find(m>=0); 

%Gr醘icas

scrsz = get(0,'ScreenSize');

figure('Position',[1 scrsz(4)/2 scrsz(3)/3 scrsz(4)/3], 'Name','Distribuci髇 de Rayleigh','NumberTitle','off')
hist(abs(R))
title('Distribuci髇 Rayleigh')
whitebg([0 .3 .6]);

figure('Position',[450 scrsz(4)/2 scrsz(3)/3 scrsz(4)/3], 'Name','Suma de Sinoidales','NumberTitle','off')
%figure('Name','Suma de Sinoidales','NumberTitle','off')
plot (abs(R))
title('Suma de sinoidales')
grid on
whitebg([1 .3 .5]);

figure('Position',[900 scrsz(4)/2 scrsz(3)/3 scrsz(4)/3], 'Name','Respuesta al Impulso','NumberTitle','off')
plot(abs(ifft(fft(R))));
title('Respuesta al impulso')
grid on
whitebg([.1 .5 1]);

figure('Position',[1 scrsz(1)/2 scrsz(3)/3 scrsz(4)/3], 'Name','Densidad de Potencia','NumberTitle','off')
plot(10*log10((abs(fft(R))).^2))
title('Densidad de potencia del espectro')
grid on
whitebg([.61 .61 .61]);

figure('Position',[450 scrsz(1)/2 scrsz(3)/3 scrsz(4)/3], 'Name','Magnitud de Potencia','NumberTitle','off')
plot(10*log10((abs(xcorr(R))).^2))
title('Magnitud de potencia (dB)')
grid on
whitebg([.9 .71 1]);

figure('Position',[900 scrsz(1)/2 scrsz(3)/3 scrsz(4)/3], 'Name','Correlaci髇','NumberTitle','off')
plot(m(m1),corr(m1));
title('Correlaci髇')
grid on
whitebg([.8 .81 .32]);
colormap(hot(256))