clc;
clear all;
close all;

%% Variablen

c0      =   physconst('lightspeed'); 
f       =   2.45e9; 
lambda  =   c0/f; 
Rgruppe =   6;                        %entfernung
Ptdbm   = 14.8;                       %Gesendete Leistung mit Berücksichtigter Kabeldaempfung
Gt      = 9.8                         %dbi GewinnHorn 9.5499;   

fid = fopen('GruppenBackplane_0degree.txt'); 


matrix2 = fscanf(fid,'%e %e',[73 73]); 
matrix2 = matrix2';

    Prdb1 =   -(matrix2(2,:)); 
    Theta1 =  (matrix2(1,:)); 
    close   
 
P1   = Prdb1 - Ptdbm;    
Pl1  = 10.^(P1./10);
G1   = (1./Gt).*Pl1.*((4.*pi.*Rgruppe).^2./(lambda.^2));
Gr1  = 10.*log10(G1);

% Gr1 = P1  - Gt + 2.*(log(4)+log(pi)+log(Rgruppe)) - 2.*log(lambda)    %Wolfgangs Formel


figure(1);
theta = Theta1;
rho= Gr1;
dirplot(theta,rho,'r-');
title('Antenna Gain in dB');

figure(2);
hold on
plot(Theta1, Gr1,'r-') %