project24.m

shadowing effects in wireless channel

% INITIALIZE ===========================================================

close all
clear

% INPUT PARAMETERS ==================================================

fMHz=2000;  % frequency in MHz


% SECONDARY PARAMETERS ==============================================

f=fMHz*1e6;
lambdac=300/fMHz;
kc=2*pi/lambdac;


% GEOMETRIC INPUTS =================================================

elev_a=30;     % elevation (deg)
azim_a=10;     % azimuth  (deg)

elev_b=30;
azim_b=0;

elev_c=30;
azim_c=-30;

D=10;        % distance from MS to building
hRx=1.5;      % MS antenna hieght

% buldings

hB1=13;       % bulding height (m)
x1=0;        % beginning of building
x2=20;       % end of building

hB2=10;       % bulding 2 height (m)
x3=30;        % beginning of building 2
x4=45;        % end of building 2

hB3=15;       % bulding 3 height (m)
x5=50;        % beginning of building 3
x6=70;        % end of building 3
 
hB4=12;       % bulding 4 height (m)
x7=80;        % beginning of building 4
x8=90;        % end of building 4

% mobile route ========================================================

stepxRx=lambdac;
xRx=[-20:stepxRx:110];

% Building outline  

BuidlingOutline=[min(xRx),0; x1,0; x1,hB1; x2,hB1; x2,0; ...
    x3,0; x3,hB2; x4,hB2; x4,0; x5,0; x5,hB3; x6,hB3; x6,0; ...
    x7,0; x7,hB4; x8,hB4; x8,0; max(xRx),0];

% Intersection point of ray with screen =============================

x0a=D*tand(azim_a)+xRx;
y0a=D*tand(elev_a)/cosd(azim_a)+hRx;
d2a=(y0a-hRx)/sind(elev_a);

x0b=D*tand(azim_b)+xRx;
y0b=D*tand(elev_b)/cosd(azim_b)+hRx;
d2b=(y0b-hRx)/sind(elev_b);

x0c=D*tand(azim_c)+xRx;
y0c=D*tand(elev_c)/cosd(azim_c)+hRx;
d2c=(y0c-hRx)/sind(elev_c);

% Fresnel radius

R1a=sqrt(lambdac*abs(d2a));
R1b=sqrt(lambdac*abs(d2b)); 
R1c=sqrt(lambdac*abs(d2c));

% Integration limits ==================================================

% vertical 
% link a
v11a=sqrt(2)*(0-y0a)/R1a;
v21a=inf;

v12a=sqrt(2)*(hB1-y0a)/R1a;
v22a=inf;

v13a=sqrt(2)*(0-y0a)/R1a;
v23a=inf;

v14a=sqrt(2)*(hB2-y0a)/R1a;
v24a=inf;

v15a=sqrt(2)*(0-y0a)/R1a;
v25a=inf;

v16a=sqrt(2)*(hB3-y0a)/R1a;
v26a=inf;

v17a=sqrt(2)*(0-y0a)/R1a;
v27a=inf;

v18a=sqrt(2)*(hB4-y0a)/R1a;
v28a=inf;

v19a=sqrt(2)*(0-y0a)/R1a;
v29a=inf;


% link b
v11b=sqrt(2)*(0-y0b)/R1b;
v21b=inf;

v12b=sqrt(2)*(hB1-y0b)/R1b;
v22b=inf;

v13b=sqrt(2)*(0-y0b)/R1b;
v23b=inf;

v14b=sqrt(2)*(hB2-y0b)/R1b;
v24b=inf;

v15b=sqrt(2)*(0-y0b)/R1b;
v25b=inf;

v16b=sqrt(2)*(hB3-y0b)/R1b;
v26b=inf;

v17b=sqrt(2)*(0-y0b)/R1b;
v27b=inf;

v18b=sqrt(2)*(hB4-y0b)/R1b;
v28b=inf;

v19b=sqrt(2)*(0-y0b)/R1b;
v29b=inf;

% link c
v11c=sqrt(2)*(0-y0c)/R1c;
v21c=inf;

v12c=sqrt(2)*(hB1-y0c)/R1c;
v22c=inf;

v13c=sqrt(2)*(0-y0c)/R1c;
v23c=inf;

v14c=sqrt(2)*(hB2-y0c)/R1c;
v24c=inf;

v15c=sqrt(2)*(0-y0c)/R1c;
v25c=inf;

v16c=sqrt(2)*(hB3-y0c)/R1c;
v26c=inf;

v17c=sqrt(2)*(0-y0c)/R1c;
v27c=inf;

v18c=sqrt(2)*(hB4-y0c)/R1c;
v28c=inf;

v19c=sqrt(2)*(0-y0c)/R1c;
v29c=inf;


% horizontal ======================================

% link a
u11a=-inf;
u21a=sqrt(2)*(x1-x0a)/R1a;
% figure,plot(xRx,u21a)

u12a=sqrt(2)*(x1-x0a)/R1a;
u22a=sqrt(2)*(x2-x0a)/R1a;
% figure,plot(xRx,u12a,'r',xRx,u22a,'b')

u13a=sqrt(2)*(x2-x0a)/R1a;
u23a=sqrt(2)*(x3-x0a)/R1a;
% figure,plot(xRx,u13a,'r',xRx,u23a,'b')

u14a=sqrt(2)*(x3-x0a)/R1a;
u24a=sqrt(2)*(x4-x0a)/R1a;
% figure,plot(xRx,u14a,'r',xRx,u24a,'b')

u15a=sqrt(2)*(x4-x0a)/R1a;
u25a=sqrt(2)*(x5-x0a)/R1a;
% figure,plot(xRx,u15a,'r',xRx,u25a,'b')

u16a=sqrt(2)*(x5-x0a)/R1a;
u26a=sqrt(2)*(x6-x0a)/R1a;
% figure,plot(xRx,u16a,'r',xRx,u26a,'b')

u17a=sqrt(2)*(x6-x0a)/R1a;
u27a=sqrt(2)*(x7-x0a)/R1a;
% figure,plot(xRx,u17a,'r',xRx,u27a,'b')

u18a=sqrt(2)*(x7-x0a)/R1a;
u28a=sqrt(2)*(x8-x0a)/R1a;
% figure,plot(xRx,u18a,'r',xRx,u28a,'b')

u19a=sqrt(2)*(x8-x0a)/R1a;
u29a=inf;
% figure,plot(xRx,u19a)

% link b
u11b=-inf;
u21b=sqrt(2)*(x1-x0b)/R1b;
% figure,plot(xRx,u21b)

u12b=sqrt(2)*(x1-x0b)/R1b;
u22b=sqrt(2)*(x2-x0b)/R1b;
% figure,plot(xRx,u12b,'r',xRx,u22b,'b')

u13b=sqrt(2)*(x2-x0b)/R1b;
u23b=sqrt(2)*(x3-x0b)/R1b;
% figure,plot(xRx,u13b,'r',xRx,u23b,'b')

u14b=sqrt(2)*(x3-x0b)/R1b;
u24b=sqrt(2)*(x4-x0b)/R1b;
% figure,plot(xRx,u14b,'r',xRx,u24b,'b')

u15b=sqrt(2)*(x4-x0b)/R1b;
u25b=sqrt(2)*(x5-x0b)/R1b;
% figure,plot(xRx,u15b,'r',xRx,u25b,'b')