Antenne Directivité Calcul par des données de terrain

[winglj]

Salut, tout le monde.J'ai une question de base sur l'utilisation du programme MATLAB pour calculer la directivité de l'antenne en utilisant les données de terrain.

Les données de terrain sont stockés dans un par-M tableau N de grandeur du champ E, N est l'angle phi de 0 degré à 360 degré et M est l'angle theta de 0 degré à 180 degrés.

Mon calcul est essentiellement la conversion de l'ampleur sur le terrain de l'intensité du rayonnement, puis en additionnant l'intensité du rayonnement de la puissance apparente rayonnée totale, puis en divisant l'intensité du rayonnement de la puissance rayonnée moyenne de directivité.

Cependant, j'ai toujours trouvé un certain décalage entre mes calculs pour la source isotrope et les résultats d'autres logiciels.Ce que quelqu'un regarde mon programme et me donner quelques conseils?Merci d'avance.

Les programmes MATLAB sont affichés ci-dessous:

%%%%%%%%%%%%%%%%%%%%%%%%%%
fonction [RP_DirDB] = RP_Dir (RP_Mag)
% Cette fonction convertit un bouquet de magitude champ à la directivité de puissance

Entrée%:
RP_Mag%: 1-by-3 cellules, chaque cellule contient un diagramme de rayonnement
% (Total, theta, phi-composants), représenté par magitude domaine.
% Sortie:
RP_DirDB%: 1-by-3 cellules, chaque cellule contient un diagramme de rayonnement
% (Total, theta, phi-composants), représenté par le pouvoir
directivité%.P = TotalP (RP_Mag (1));
RP_DirDB (1) = EMag2DirDB (RP_Mag (1), P);

% Les commandes suivantes sont utilisées pour les sous-composants de terrain
(2) = RP_DirDB EMag2DirDB% (RP_Mag (2), P);
(3) = RP_DirDB EMag2DirDB% (RP_Mag (3), P);

fonction [P, S] = TotalP (Mag)
% Cette fonction obtient la puissance totale rayonnée normalisée de la magitute total de E
domaines%

Entrée%:
% Mag: N_phi par N_theta matrice contient le diagramme de rayonnement dans
% En valeur absolue.
% Sortie:
% P: puissance totale rayonnée normalisé
% S: surface intégrale

% Obtenir l'intensité du rayonnement
Unor = abs (Mag). ^ 2;

% Obtenez le phi et theta de données
[NPHI, Ntheta] size = (Mag);
Phi = linspace (0,2 pi *, NPHI);
Dphi = 2 * pi / (NPHI-1);
Theta = linspace (0, pi, Ntheta);
Dtheta = pi / (Ntheta-1);

% Initialisation de la puissance apparente rayonnée totale P, zone intégrante S, dP puissance delta
% Et Detal dS zone intégrante.
P = 0;
S = 0;
dP = [];
dS = [];

% Calculer la puissance d'irradiation totale et l'intensité du rayonnement
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);

pour i_theta = 2: (Ntheta-1)
pour i_phi = 1: (INSP-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);
fin
fin

fonction [DirDB] = EMag2DirDB (Mag, P)
% Cette fonction convertit le nombre de décibels magnitude absolue nombre

Entrée%:
% Mag: N_phi par N_theta matrice contient le diagramme de rayonnement dans
% En valeur absolue.
% Sortie:
% DirDB: N_phi par N_theta matrice contient le diagramme de rayonnement dans
valeur en décibels%% Obtenir l'intensité du rayonnement
Unor = abs (Mag). ^ 2;

Renormaliser% l'intensité du rayonnement à directivité
Dir = 4 * pi * Unor / P;

% Calculer la directivité en décibels
DirDB = 10 * log10 (Dir);