funcaotransicao.m

计算GTD-UTD绕射源代码

function [saida, s_saida] = funcaotransicao( x, y, frequencia)
%
% Gera玢o da fun玢o de transi玢o complexa F(S) 
% que torna o coeficiente de difrac玢o 
% Dt(teta)=D(teta).F(S) cont韓ua e tata = teta(x,y).
%
% function [saida, s] = funcaotransicao( x, y, frequencia)
%
% saida       : Vector de sa韉a com o fun玢o de transi玢o F(S);
% s_saida     : Vector de sa韉a com os valores do parametro 'S';
% x	        : Dist鈔cia horizontal da descontinuidade espacial;
% y	        : Vector com as eleva珲es em rela玢o ao plano horizontal
%		          que cont閙 a descontinuidade espacial (edge);
% frequencia  : Valor da frequencia de trabalho (MHz).
%

saida = y*0;
s_saida = y*0;
s = 0;
tgteta2 = 0;
ro = 0;
k = 0;
fs = 0;
lambda = 300/frequencia;
beta = 2*pi/lambda;
csi = 0;
k1 = 0;
k2 = 0;
k3 = -1/sqrt(2*pi*beta);
f = 0;
g = 0;
num1 = 0;
den1 = 0;
num2 = 0;
den2 = 0;
cdaux = 0;
num3 = 0;
den3 = 0;
costeta = 0;
senteta = 0;
for k = 1 : length(y),
  ro = power( (x*x + y(k)*y(k)), 0.5);
  tgteta2 = power( (y(k)/x), 2);   
  s = beta*ro*tgteta2/2;
  k1 = power( (2*pi*s), 0.5);
  csi = power( (2*s/pi), 0.5);
  num1 = 1 + 0.926*csi;
  den1 = 2 + 1.792*csi + 3.104*power( csi, 2);
  f = num1/den1;
  num2 = 1;
  den2 = 2 + 4.142*csi + 3.492*power( csi, 2) + 6.670*power( csi, 3);
  g = num2/den2;
  k2 = f + j*g;
  fs = k1*k2;
  saida(k) = fs;
  s_saida(k) = s;
end