Grado en Ingeniería de Tecnologías de Telecomunicación
Prácticas de Radiación y Propagación
PRÁCTICA 5 – abril 2023
DISEÑO DE UN ARRAY (900 MHz)
Nombre: Manuel Alejandro Perulero Casal
DNI: 49760206M
N1 N2 N3 N4 N5 N6 N7 N8
4 9 7 6 0 2 0 6
Objetivo:
El objeto de esta práctica es realizar el diseño de una antena para una estación base de
telefonía móvil (900 MHz) colocada sobre un cerro de 3N8N7N6/3= 1200.67 metros de
altura (h). Se pretende obtener una cobertura omnidireccional sobre una gran llanura que
lo rodea.
Cuestiones:
A. Calcule cómo deben alimentarse los dipolos de un array de 4 dipolos de longitud λ/2
separados d=(0.7N5)λ = 23.31 cm [ λ =0.33m] entre centros capaz de conseguir la
máxima ganancia sobre una zona situada a D=10Km del cerro.
ℎ 1200.67
𝑡𝑔(𝜃) = = = 0.120067
𝐷 10000
𝜽 = 𝑎𝑟𝑐𝑡𝑔(0.120067) = 𝟔. 𝟖𝟓º
Angulo de máxima radiación = 6.85º
𝑑𝑒𝑠𝑓𝑎𝑠𝑒 = 𝛽 = −𝑘 ∗ 𝑑 ∗ cos(𝜃)
2𝜋
𝑘=
𝜆
2𝜋
𝜷 = −𝑘 ∗ 𝑑 ∗ cos(𝜃) = − ∗ 0.2331 ∗ cos(6.85º) = −𝟒. 𝟑𝟔𝟐𝟔𝒓𝒂𝒅
0.33
Para conseguir la máxima ganancia sobre una zona situada a 10 km del cerro, los
dipolos de un array de 4 dipolos de longitud = 16.65cm, separados 23.31cm entre
centros, deben ser alimentados con una diferencia de fase de -4.39 radianes entre
ellos
1
B. Adjunte la representación del diagrama de radiación de la antena en coordenadas
polares, en escala lineal, detallando (calculándolos) la posición angular de los nulos y
el nivel de los lóbulos secundarios; emplee para ello Matlab e incluya el código del
programa desarrollado.
El código de Octave utilizado es:
% Parámetros del array
N = 4;
lambda = 0.333;
d = 0.70 * lambda;
theta0 = 6.85; % grados
k = 2 * pi / lambda;
beta = -k * d * cosd(theta0);
% Ángulos
theta = linspace(0, 180, 1000); % grados
psi = k * d * cosd(theta) + beta;
% Función de array
AF = abs(sin(N * psi / 2) ./ sin(psi / 2));
AF = AF / max(AF); % Normalización
% --- Representación en coordenadas polares (lineal) ---
theta_rad = deg2rad(theta);
figure;
polar(theta_rad, AF, 'b');
title('Diagrama de Radiación (Polar, escala lineal)');
grid on;
% --- Representación en coordenadas cartesianas (lineal) ---
figure;
plot(theta, AF, 'r', 'LineWidth', 2);
xlabel('Ángulo (°)');
ylabel('Ganancia normalizada');
title('Diagrama de Radiación (Cartesianas, escala lineal)');
grid on;
2
% --- Cálculo de nulos ---
n = 1:(N - 1);
theta_nulls = [];
for i = 1:length(n)
arg = (-beta + 2 * pi * n(i)) / (k * d);
if abs(arg) <= 1
theta_nulls(end+1) = acosd(arg);
end
end
disp(' Nulos del diagrama (en grados):');
disp(theta_nulls');
% Calcular AF en dB para la gráfica
AF_dB_plot = 20 * log10(AF);
AF_dB_plot(AF_dB_plot < -60) = -60;
% Buscar picos (lóbulos) en AF lineal
[pks, locs] = findpeaks(AF, 'MinPeakHeight', 0.05); % descarta ruido
% Eliminar el lóbulo principal (el mayor pico)
[~, max_idx] = max(pks);
pks_sec = pks;
pks_sec(max_idx) = [];
% Convertir picos secundarios a dB
pks_sec_dB = 20 * log10(pks_sec);
SLL_dB = max(pks_sec_dB);
disp(['Nivel máximo de lóbulo secundario (SLL): ', num2str(SLL_dB), ' dB']);
% Mostrar posición angular de los nulos y nivel de los lóbulos secundarios
for i = 1:length(pks_sec)
disp(['Lóbulo secundario ', num2str(i), ': Nivel = ', num2str(pks_sec_dB(i)), ' dB']);
end
3
Con esto obtenemos que:
Figura 1. Diagrama de radiación
en polares
Figura 2. Diagrama de radiación
en cartesianas
o Nulos del diagrama (en grados): No se encontraron nulos en el rango de
ángulos especificado.
o Nivel máximo de lóbulo secundario (SLL): 56.16 dB
o Lóbulos secundarios:
• Lóbulo secundario 1: Nivel = 31.60 dB
• Lóbulo secundario 2: Nivel = 50.73 dB
• Lóbulo secundario 3: Nivel = 53.48 dB
• Lóbulo secundario 4: Nivel = 56.16 dB
4
C. Estime la directividad del array a partir de los anchos de haz de la antena.
𝜆
El ancho de haz lo aproximamos como: 𝐵𝑊−3𝑑𝐵 = 0.886 ∗ 𝑁∗𝑑∗sin(𝜃)
0.333
𝑩𝑾−𝟑𝒅𝑩 = 0.866 ∗ = 2.66𝑟𝑎𝑑 = 𝟏𝟓𝟐. 𝟎𝟏º
4 ∗ 0.2331 ∗ sin(6.85º)
Una vez calculado el ancho de haz de la antena, podemos hacer una aproximación de
4𝜋 𝐵𝑊−3𝑑𝐵 2 2.66 2
su directividad utilizando: 𝐷 ≅ ; 𝛺𝐴 = ( ) =( ) = 0.00205
𝛺𝐴 57.3 57.3
4𝜋
Por tanto, 𝑫 ≅ ≅ 𝟓𝟖𝟔𝟏. 𝟖𝟑
0.00205
De esta manera obtenemos que el ancho de haz de la antena a -3dB es de 152.01º y
su directividad, aproximadamente, igual a 5861.83.
D. Utilizando el programa 4Nec2 compruebe y compare el diagrama de radiación y su
estimación de la directividad de la antena. Capture y adjunte el diagrama obtenido.
¿El máximo del diagrama de radiación de la antena completa se obtiene para la
dirección calculada? ¿Por qué?
La visualización y parámetros de mi simulación es la siguiente:
Figura 3. Visualización de la simulación
5
Figura 4. Parámetros de la simulación
Los resultados son los siguientes:
Figura 5. Plano E Figura 6. Plano H
Los resultados son muy diferentes de los obtenidos teóricamente y en Octave si los
comparamos con los de la simulación, entiendo que debido a que teóricamente
utilizamos modelos simplificados mientras que 4nec2 utiliza el método de
momentos. Además, nosotros presuponemos elementos continuos y perfectos
mientras que el programa utiliza elementos segmentados.
Sin embargo, los resultados difieren tanto que no me da mucha seguridad de que el
procedimiento en cualquiera de los métodos haya sido el correcto.
