difractionlossgraphfresnel.m

fortran 解复方程的代码

function results = difractionlossgraphfresnel( vinicial, vfinal, pontos)
%
% FUN敲O PARA O CALCULO E DESENHO DO GR罠ICO COM O VALOR DAS
% DAS PERDAS POR DIFRAC敲O, PARA UM A GAMA DO PARAMETRO 'v'
%
% Fonte bibliogr醘ica: The Mobile Radio Propagatio Channel
%                      J. D. Parsons
%                      Pentech Press,1992.                  
%                      P醙ina 43
%
% sa韉a = difractionlossgraphfresnel( vinicial, vfinal, pontos)
%
% sa韉a      : '0' se sucesso, '1' se ocorreu erro.
% vinicial   : Valor inicial para o parametro 'v'
% vfinal     : Valor final para o parametro 'v'
% pontos     : Nr. de pontos que comp鮡 a espiral
%

results = 1;
auxiliar1 = 0;
auxiliar2 = 0;
incremento =(vfinal - vinicial)/(pontos);
x = vinicial:incremento:vfinal;
y = vinicial:incremento:vfinal;
z = vinicial:incremento:vfinal;
for n = 1 : length(x),
  v = vinicial + (n-1) * incremento;
  x(n) = cosenointegral( v);
  y(n) = senointegral( v);
  auxiliar1 = power( 2, 0.5)/2;
  auxiliar2 = power((0.5 - x(n)), 2);
  auxiliar2 = auxiliar2 + power((0.5 - y(n)), 2);
  auxiliar2 = power(auxiliar2, 0.5);
  z(n) = 20*log10(auxiliar1) + 20*log10(auxiliar2);
end
x = vinicial:incremento:vfinal;
plot(x,z);
results = 0;