aHemi_fr

%--------------------------------------------------------------------------
%                                aHemi_fr.m                                   
%--------------------------------------------------------------------------
%  Analyse d'un rsonateur hmisphrique (franais)                                              
%
%    A) ENTRES
%    
%    1. 'a' = rayon du rsonateur (en cm)
%    2. 'er' = constante dilectrique relative du rsonateur
%    3. 'VSWR' = taux d'ondes stationnaires tolr pour le calcul    
%       de la bande passante d'impdance
%
%    B) SORTIES
%    
%    - Frquence de rsonance, facteur Q et bande passante d'impdance
%      pour les deux premiers modes : TE111 et TM101.
%--------------------------------------------------------------------------
%  Rfrence :
%
%  - M. Gastine, L. Courtois et J.J. Dormann, "Electromagnetic resonances
%    of free dielectric spheres", IEEE Transactions on microwave theory
%    and techniques, Vol. 15, Num. 12, dcembre 1967, pp. 694-700.
%--------------------------------------------------------------------------
%  Note : Ce fichier fait partie de l'Outil de conception et d'analyse 
%         des rsonateurs dilectriques (DRA.m). Les quations d'estimation
%         sont tires du livre "Dielectric resonator antenna handbook" de
%         Also Petosa.
%
%  Programmation : Alexandre Perron ([e-mail])               
%  Affiliation : Institut national de la recherche scientifique (INRS)  
%  Dernire modification : Le 22 juillet 2008                              
%--------------------------------------------------------------------------

function aHemi_fr

% Entte :
clc
disp(strvcat('===============================================================',...
             ' Outil de conception et d''analyse des rsonateur dilectriques',...
             '==============================================================='));
disp(sprintf('\n----------- Analyse d''un rsonateur hmisphrique -------------'));

% Entre des paramtres :
a = [];
while isempty(a)||(isnumeric(a) == 0),
   a = input('\nEntrez le rayon (a) du rsonateur (en cm) : ');
end

er = [];
while isempty(er)||(isnumeric(er) == 0),
   er = input('\nEntrez la constante dilectrique (er) du rsonateur : ');
end

VSWR = [];
while isempty(VSWR)||(isnumeric(VSWR) == 0),
   VSWR = input('\nEntrez le VSWR pour le calcul de la bande passante (ex.: 2) : ');
end

%============================%
% Calculs pour le mode TE111 %
%============================%

% Frquence de rsonance (en GHz) :
freq_TE = 2.8316*er^-0.47829*4.7713/a;
% Facteur Q :
facteurQ_TE = 0.08+0.796*er+0.01226*er^2-3e-5*er^3;
% Bande passante :
BW_TE = (VSWR-1)/(sqrt(VSWR)*facteurQ_TE)*100;

%============================%
% Calculs pour le mode TM101 %
%============================%

% Frquence de rsonance (en GHz) :
freq_TM = 4.47226*er^-0.505*4.7713/a;
% Facteur Q :
if (er <= 20)
    facteurQ_TM = 0.723+0.9324*er-0.0956*er^2+0.00403*er^3-5e-5*er^4;
else
    facteurQ_TM = 2.621-0.574*er+0.02812*er^2+2.59e-4*er^3;
end
% Bande passante :
BW_TM = (VSWR-1)/(sqrt(VSWR)*facteurQ_TM)*100;

% Entte :
clc
disp(strvcat('===============================================================',...
             ' Outil de conception et d''analyse des rsonateur dilectriques',...
             '==============================================================='));
disp(sprintf('\n----------- Analyse d''un rsonateur hmisphrique -------------\n'));

% Affichage des paramtres d'entre :
disp(sprintf('Rayon (a) du rsonateur (en cm) = %5.4f',a));
disp(sprintf('Constante dilectrique du rsonateur = %5.4f',er));
disp(sprintf('VSWR pour le calcul de la bande passante = %5.4f',VSWR));

% Affichage des rsultats :
disp(sprintf('\n'));
disp(strvcat('                   Mode TE111                    ', ...
             '-------------------------------------------------'));
disp(sprintf('    Frquence de rsonance (en GHz) = %5.4f',freq_TE));
disp(sprintf('                          Facteur-Q = %5.4f',facteurQ_TE));
disp(sprintf('    Bande passante (en pourcentage) = %5.4f',BW_TE));

disp(sprintf('\n'));
disp(strvcat('                   Mode TM101                    ', ...
             '-------------------------------------------------'));
disp(sprintf('    Frquence de rsonance (en GHz) = %5.4f',freq_TM));
disp(sprintf('                          Facteur-Q = %5.4f',facteurQ_TM));
disp(sprintf('    Bande passante (en pourcentage) = %5.4f',BW_TM));

end