AQE-WMMSE/MATLAB/src/generateLOSchannel.m at main · AlexanderFernandes96/AQE-WMMSE · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
function [H] = generateLOSchannel(M,aoa,N,aod,LOS)
%Generate M x N matrix with LOS paths
%   aoa and aod:
%   1 for linear array, 2 for rectangular array
%   for rectuangular array use M = [Mw, Mh] or N = [Nw, Nh];
%   assume spacing between adjacent antennas is d = lambda/2


if aoa == 1
    Msize = M;
elseif aoa == 2
    Msize = M(1) * M(2);
else
    error("aoa must be 1, or 2.")
end
if aod == 1
    Nsize = N;
elseif aod == 2
    Nsize = N(1) * N(2);
else
    error("aoa must be 1, or 2.")
end
H = zeros(Msize,Nsize);

for l = 1:LOS
    h = sqrt(1/2)*(randn(1,1) + i*randn(1,1)); % path gain
    if aoa == 1
        theta = unifrnd(0,pi/2);
        AOA = aULA(theta,M);
    elseif aoa == 2
        theta = unifrnd(0,pi/2);
        psi = unifrnd(0,pi);
        AOA = aURA(theta, psi, M(1), M(2));
    end

    if aod == 1
        theta = unifrnd(0,pi/2);
        AOD = aULA(theta,N);
    elseif aod == 2
        theta = unifrnd(0,pi/2);
        psi = unifrnd(0,pi);
        AOD = aURA(theta, psi, N(1), N(2));
    end

    H = H + h*AOA*AOD';
end

end