ecrancondutor.m

计算GTD-UTD绕射源代码

function saida = ecrancondutor( frequencia, angulo, filinha)
%
%  Fun玢o que para um determinado valor da frequ阯cia de trabalh,
%  e para uma ampla gama de angulos, retorna o valor do
%  coeficiente de difra玢o Df=D(teta), usando a formula玢o de FELSEN..
%  Depende da polariza玢o.
%
%
%
%  Fonte bibliogr醘ica: Radiowave propagation for modern communications
%			Henry L. Bertoni
%			Prentice-Hall, 2000
%			[Bertoni 2000]
%
%  saida = ecrancondutor( frequencia, angulo, filinha)
%  saida	       : Vector com valore do coeficiente de difusao, para
%		            uma gama de angulos dada pelo array 'angulo';
%  frequencia   : Valor da frequencia de trabalho expressa em MHz;
%  angulo	    : Vector de entrada contendo valores dos 鈔gulos
%                 'teta entre a ISB e a posi玢o;
%  filinha      :
%

k = 0;

s = 0;

lambda = 300/frequencia;

beta = 2*pi/lambda;

saida = 0*angulo;

aux = 0;

for k = 1 : length(angulo),
   teta = angulo(k);
   aux = -0.5/sqrt(2*pi*beta)*(1/cos((pi-teta)/2)-1/cos((pi-teta+2*filinha)/2));  
   saida(k) = abs(aux);
end