v = zerosNM lengththetalengthphi;for ii = 1 lengththeta for jj = 1 lengthphi k = 2 pi lambda sinthetaii pi 180 cosphijj pi 180 + lambda lambdagsinthetaii pi 180 sinphijj pi 1

function [v] = antenna_pattern(theta, phi, lambda, lambdag, N, M, m0) v = zeros(N*M, length(theta)*length(phi)); for ii = 1 : length(theta) for jj = 1 : length(phi) k = 2 * pi / lambda * [sin(theta(ii) * pi / 180) * cos(phi(jj) * pi / 180) + lambda / lambdag,sin(theta(ii) * pi / 180) * sin(phi(jj) * pi / 180) + lambda / lambdag , 0].'; v(:, (ii-1)length(phi) + jj) = m0 * exp(1j * k.' * (1:NM)); end end end