% HLANT1.M - Entwurf von Hohlleiterschlitzantennen 
%  26.8.1996 Juergen Hasch, hasch@t-online.de
% Nach:
% O.Nell, K. Solbach, J. Dreier: Rundumstrahlende Hohlleiter-Schlitzantenne
%                                für Horizontalpolarisation
% UKW-Berichte 1/91 und 2/91
%
disp('Programm zur Berechnung von Hohlleiterschlitzstrahlern');
disp('--------------------- 8/96 ---------------------------');
Z_F0=120.0*pi;
C0=2.997831e8;
f0 = input('Frequenz in MHz:');
lambda_0= C0/f0/1e6;
 t=sprintf('Freiraumwellenlaenge in mm: %6.2f',lambda_0*1e3);
 disp(t);
 t=sprintf('Cutoff-Bereich bei a: %6.2f',lambda_0*1e3/2);
 disp(t);
 t=sprintf('Optimale Breite von a: %6.2f',lambda_0*sqrt(3)*1e3/2);
 disp(t);
%
% Hohlleiter-Laenge und -Breite eingeben
%
a = input('a in mm:');
b = input('b in mm:');
lambda_c=2.0*a/1e3;
lambda_h= lambda_0 / sqrt(1.0 - (lambda_0/lambda_c)^2);
N = input('Anzahl der Schlitze:');
 t=sprintf('lambda_h in mm: %3.2f',lambda_h*1e3);
 disp(t);
 t=sprintf('lambda_0 in mm: %3.2f',lambda_0*1e3);
 disp(t);
%
% Gewinn und Keulenweite berechnen
%
g= 10*log10(N * lambda_h/lambda_0);
 t=sprintf('Ungefaehrer Antennengewinn in dBi: %3.1f',g);
 disp(t);
dtheta=50.7 * lambda_0 / N / lambda_h * 2.0;
 t=sprintf('Keulenweite in der Elevation: %3.1f',dtheta);
 disp(t);
%
% Mittenversatz berechnen
%
x = 1e-3 * a/pi*asin(sqrt(lambda_0*b/(N*3.5*lambda_h*a*cos(pi*lambda_0/2.0/lambda_h).^2)));
  t = sprintf('Mittenversatz: %6.2f',x*1e3);
  disp(t);
%
% Schlitzbreite
%
 sb = lambda_h / 30;
  t = sprintf('Schlitzbreite: %6.2f',sb*1e3);
  disp(t);
%
% Schlitzlaenge ermitteln
%
% Stuetzpunkte aus Bild 7; gelten für 1,27GHz, deshalb Umrechnung
%
x_23cm = f0/1270 * x;
sx=[ 0    8.6   10    20   29 ];   % x in mm
sl=[ .475 .491  .494  .519  .55];  % L/lambda_0
L = interp1(sx,sl,x_23cm*1e3,'spline');
 t = sprintf('Schlitzlaenge: %6.2f',L*1e3*lambda_0);
 disp(t);
%
% Schlitze berechnen
%
disp('   Schlitz-Nr    Abstand vom Kurzschluss');
disp('----------------------------------------');
for i=0:(N-1)
  t = sprintf('   %d        %6.2f',i+1,1e3*(lambda_h/4+lambda_h/2*i));
  disp(t)
end
%


