Antenne Directivity Beregning av feltdata

[winglj]

Hei, alle.Jeg har et grunnleggende spørsmål om bruk MATLAB program å beregne antenne directivity med feltdata.

Feltet data blir lagret i en N-by-M rekke omfanget av E-feltet, N er phi vinkel fra 0 grader til 360 grader og M er Theta vinkel fra 0 grader til 180 grader.

Mitt regnestykke er utgangspunktet konvertere feltet styrke til stråling intensitet og summing av stråling intensiteten til totalt radiated makt, deretter dele stråling intensiteten til gjennomsnittet radiated strøm til directivity.

Men jeg har alltid funnet enkelte avvik mellom mine beregninger for isotropic kildekode og annen programvare resultater.Kan noen se på programmet og gi meg noen tips?Takk på forhånd.

Den MATLAB programmene er vist nedenfor:

%%%%%%%%%%%%%%%%%%%%%%%%%%
function [RP_DirDB] = RP_Dir (RP_Mag)
% Denne funksjonen konverterer en rekke felt magitude til makten directivity

% Input:
% RP_Mag: 1-av-3 celler, hver celle inneholder en stråling mønster
% (Totalt, Theta, Phi-komponenter) representert ved feltet magitude.
% Output:
% RP_DirDB: 1-av-3 celler, hver celle inneholder en stråling mønster
% (Totalt, Theta, Phi-komponenter) representert ved makten
% Directivity.[P, S] = TotalP (RP_Mag (1));
RP_DirDB (1) = EMag2DirDB (RP_Mag (1), P);

% Følgende kommandoer brukes til feltet sub-komponenter
% RP_DirDB (2) = EMag2DirDB (RP_Mag (2), P);
% RP_DirDB (3) = EMag2DirDB (RP_Mag (3), P);

function [P, S] = TotalP (Mag)
% Denne funksjonen får normalisert totale radiated kraft fra den totale magitute av E
% Felt

% Input:
% Mag: N_phi-by-N_theta matrise, inneholder stråling mønster
% Absoluttverdien.
% Output:
% P: Normaliserte totalt radiated strøm
% S: Integral areal

% Få stråling intensitet
Unor = abs (Mag). ^ 2;

% Få Phi og Theta data
[Nphi, Ntheta] = size (Mag);
Phi = linspace (0,2 * pi, Nphi);
Dphi = 2 * pi / (Nphi-1);
Theta = linspace (0, pi, Ntheta);
Dtheta = pi / (Ntheta-1);

% Initialiser totalt radiated strøm P, integrert område S, Delta strøm dP
% Og detal integral området dS.
P = 0;
S = 0;
dP = [];
dS = [];

% Calculate total stråling makt og stråling intensitet
dS (1,1) = 2 * pi * (1-cos (Dtheta / 2));
dP (1,1) = Unor (1,1) * dS (1,1);
P = P dP (1,1);
S = S dS (1,1);

dS (1, Ntheta) = 2 * pi * (1-cos (Dtheta / 2));
dP (1, Ntheta) = Unor (1, Ntheta) * dS (1, Ntheta);
P = P dP (1, Ntheta);
S = S dS (1, Ntheta);

for i_theta = 2: (Ntheta-1)
for i_phi = 1: (Nphi-1)
dS (i_phi, i_theta) = Dphi * (cos (Theta (i_theta)-Dtheta / 2)-cos (Theta (i_theta) Dtheta / 2));
dP (i_phi, i_theta) = Unor (i_phi, i_theta) * dS (i_phi, i_theta);
P = P dP (i_phi, i_theta);
S = S dS (i_phi, i_theta);
slutt
slutt

function [DirDB] = EMag2DirDB (Mag, P)
% Denne funksjonen konverterer de absolutte styrke for å desibelnivået nummer

% Input:
% Mag: N_phi-by-N_theta matrise, inneholder stråling mønster
% Absoluttverdien.
% Output:
% DirDB: N_phi-by-N_theta matrise, inneholder stråling mønster
% Desibelnivået verdi% Få stråling intensitet
Unor = abs (Mag). ^ 2;

% Renormalize av stråling intensiteten til directivity
Dir = 4 * pi * Unor / P;

% Beregn directivity i desibel
DirDB = 10 * log10 (Dir);