espiraldecornu.m

fortran 解复方程的代码

function saida = espiraldecornu( vinicial, vfinal, pontos)
%
% FUN敲O PARA O CALCULO E DESENHO DA ESPIRAL DE CORNU
%
% Fonte bibliogr醘ica: The Mobile Radio Propagatio Channel
%                      J. D. Parsons
%                      Pentech Press,1992.                    
%  		       P醙ina 42
%
% sa韉a = espiraldecornu( 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
%

incremento =(vfinal - vinicial)/(pontos);
x = vinicial:incremento:vfinal;
y = vinicial:incremento:vfinal;
for n = 1 : length(x),
  v = vinicial + (n-1) * incremento;
  x(n) = cosenointegral( v);
  y(n) = senointegral( v);
end
plot(x,y);