Antena Koeficientas dėmesio apskaičiavimas pagal lauko duomenų

[winglj]

Sveiki, visi.Aš pagrindinio klausimo apie naudojant MATLAB programa apskaičiuoti antena Koeficientas dėmesio naudojant lauko duomenis.

Srityje duomenys yra saugomi ir N-by-M masyvo dydį E srityje, N yra phi kampas nuo 0 laipsnių iki 360 laipsnių ir M yra theta kampas nuo 0 laipsnių iki 180 laipsnių.

Mano apskaičiavimą iš esmės pertvarkyti srityje dydžio radiacijos intensyvumas, tada susumuojant spinduliuotės intensyvumas viso spinduliuotės galią, dalijant spinduliavimo intensyvumas vidutinė spinduliuotės galia Koeficientas dėmesio.

Tačiau, aš visada yra keletas neatitikimų tarp mano skaičiavimus isotropic šaltinį ir kita programinė įranga rezultatus.Ar kas nors atrodo ne mano programa ir duoti man patarimų?Thanks in advance.

Į MATLAB programos parodyta žemiau:

%%%%%%%%%%%%%%%%%%%%%%%%%%
funkcija [RP_DirDB] = RP_Dir (RP_Mag)
% Ši funkcija konwertuje ryšelyje lauko magitude į valdžią Koeficientas dėmesio

% Input:
% RP_Mag: 1-by-3 ląstelių, kiekviena ląstelė yra spinduliavimo dėsningumo
% (Iš viso, teta, phi-komponentams), atstovaujama srityje magitude.
% Produkcija:
% RP_DirDB: 1-by-3 ląstelių, kiekviena ląstelė yra spinduliavimo dėsningumo
% (Iš viso, teta, phi-komponentams), atstovaujama galia
% Koeficientas dėmesio.[R, S] = TotalP (RP_Mag (1));
RP_DirDB (1) = EMag2DirDB (RP_Mag (1), P);

% Šios komandos naudojamos srityje sudedamosiomis dalimis
% RP_DirDB (2) = EMag2DirDB (RP_Mag (2), P);
% RP_DirDB (3) = EMag2DirDB (RP_Mag (3), P);

funkcija [R, S] = TotalP (Mag)
% Ši funkcija gauna normalizuotą viso spinduliuotės galios iš viso magitute E
% Sritys

% Input:
% Mag: N_phi pagal N_theta matricos yra spinduliavimo dėsningumo į
% Absoliučios vertės.
% Produkcija:
% P Normalizuotas atiduotų spinduliuotės galios
% S: bendras plotas

% Gauti radiacijos intensyvumas
Únor = abs (magija). ^ 2;

% Gaukite phi ir teta duomenų
[Nphi, Ntheta] = dydis (Mag);
Phi = linspace (0,2 * pi, Nphi);
Dphi = 2 * Pi / (Nphi-1);
Theta = linspace (0, Pi, Ntheta);
Dtheta = Pi / (Ntheta-1);

% Inicijuoti bendrą spinduliuotės galios P, bendras plotas S delta galios dP
% Ir detalę bendras plotas dS.
P = 0;
S = 0;
dP = [];
DS = [];

% Paskaičiuoti bendrą spinduliuotės galią ir spinduliuotės intensyvumas
dS (1,1) = 2 * pi * (1-cos (Dtheta / 2));
dP (1,1) = Únor (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) = Únor (1 Ntheta) * dS (1 Ntheta);
P = P dP (1 Ntheta);
S = S dS (1 Ntheta);

už i_theta = 2: (Ntheta-1)
už 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) = Únor (i_phi, i_theta) * dS (i_phi, i_theta);
P = P dP (i_phi, i_theta);
S = S dS (i_phi, i_theta);
pabaiga
pabaiga

funkcija [DirDB] = EMag2DirDB (Mag P)
% Ši funkcija konwertuje absoliutus dydis numeris Decibel numeris

% Input:
% Mag: N_phi pagal N_theta matricos yra spinduliavimo dėsningumo į
% Absoliučios vertės.
% Produkcija:
% DirDB: N_phi pagal N_theta matricos yra spinduliavimo dėsningumo į
% Decibel vertė% Gauti radiacijos intensyvumas
Únor = abs (magija). ^ 2;

% Renormalize spinduliavimo intensyvumas Koeficientas dėmesio
Dir = 4 * pi * Únor / P;

% Apskaičiuokite Koeficientas dėmesio į Decibel
DirDB = 10 * log10 (DIR);