How to Generate a 3D Radiation Pattern Line-by-Line Similar to Sens... (2024)

% Define your antenna array parameters and call the function

num_ant_X = 1;

num_ant_Y = 2;

num_ant_Z = 2;:

d = 0.5; % Spacing between antennas in wavelength units

c = 3.0e8;

f = 3.5e9;

lambda_0 = c / f;

steering_theta_deg = 90;

steering_phi_deg = 0;

Pattern_array_factor_3d_func(num_ant_X, num_ant_Y, num_ant_Z, lambda_0, d, deg2rad(steering_theta_deg), deg2rad(steering_phi_deg), "omni");

function Pattern_array_factor_3d_func(num_ant_X, num_ant_Y, num_ant_Z, lambda_0, distance_between_antennas_longueur_onde, steering_theta, steering_phi, pattern)

% Calculate the positions of the elements centered around the origin

x_positions = (0:num_ant_X-1) - (num_ant_X - 1) / 2;

y_positions = (0:num_ant_Y-1) - (num_ant_Y - 1) / 2;

z_positions = (0:num_ant_Z-1) - (num_ant_Z - 1) / 2;

x_positions = x_positions * distance_between_antennas_longueur_onde * lambda_0;

y_positions = y_positions * distance_between_antennas_longueur_onde * lambda_0;

z_positions = z_positions * distance_between_antennas_longueur_onde * lambda_0;

% Elevation and azimuth angles

theta = linspace(0, pi, 180);

phi = linspace(-pi, pi, 360);

[THETA, PHI] = meshgrid(theta, phi);

% Array Factor (AF) Initialization

AF = zeros(size(THETA));

% Steering vector components

Us = sin(steering_theta) * cos(steering_phi);

Vs = sin(steering_theta) * sin(steering_phi);

Ws = cos(steering_theta);

% Loop over each antenna element

for x = x_positions

for y = y_positions

for z = z_positions

W = exp(-1j * 2 * pi * (x * Us + y * Vs + z * Ws) / lambda_0);

AF = AF + W .* exp(1j * 2 * pi * (x * sin(THETA) .* cos(PHI) + y * sin(THETA) .* sin(PHI) + z * cos(THETA)) / lambda_0);

end

end

end

% Optional radiation pattern modification (3GPP pattern, etc.)

if strcmp(pattern, "38.901")

radiation_pattern = radiation_pattern_38901(THETA, PHI); % This would need to be defined

AF_global_magnitude_squared = abs(AF).^2 .* radiation_pattern;

else

AF_global_magnitude_squared = abs(AF).^2;

end

AF_global_magnitude_squared_norm = AF_global_magnitude_squared / (num_ant_X * num_ant_Y * num_ant_Z);

% Convert to dBi

AF_magnitude_dBi = 10 * log10(AF_global_magnitude_squared_norm + 0.001);

% Convert to Cartesian coordinates for 3D plot

X = (AF_magnitude_dBi) .* sin(THETA) .* cos(PHI);

Y = (AF_magnitude_dBi) .* sin(THETA) .* sin(PHI);

Z = (AF_magnitude_dBi) .* cos(THETA);

% Plot the 3D radiation pattern

figure;

surf(X, Y, Z, 'FaceAlpha', 0.6, 'EdgeColor', 'none');

colorbar;

title('3D Radiation Pattern');

xlabel('X');

ylabel('Y');

zlabel('Z');

axis equal;

% Extract 2D cut at theta = 90 degrees (elevation pattern)

target_theta = pi / 2;

[~, closest_index_theta] = min(abs(theta - target_theta));

elevation_cut_pattern = AF_magnitude_dBi(:, closest_index_theta);

figure;

polarplot(phi, elevation_cut_pattern);

rlim([min(elevation_cut_pattern) max(elevation_cut_pattern)])

title('Elevation Pattern (dBi)');

% Extract 2D cut at phi = 0 degrees (azimuth pattern)

target_phi = 0;

[~, closest_index_phi] = min(abs(phi - target_phi));

azimuth_cut_pattern = AF_magnitude_dBi(closest_index_phi, :);

figure;

polarplot(theta, azimuth_cut_pattern);

rlim([min(azimuth_cut_pattern) max(azimuth_cut_pattern)])

title('Azimuth Pattern (dBi)');

end