Intégrales de Fresnel

%----------------------------------------------------------
% Intégrales de Fresnel
% d'aprés 'Apprendre et maitriser Matlab'
%----------------------------------------------------------

clear all;
echo on;
%---------------------------------------
%      Intégrales de Fresnel
% fonction de 2ème espèce
%         ______  x 
%  S(x)= V(2/pi) S sin t² dt      
%                 0   
% fonction de 1ère espèce
%         ______  x 
%  C(x)= V(2/pi) S cos t² dt      
%                  0   
%-------------------------------------
echo off;
pas=0.01;
x=-10:pas:10;
n=length(x);
	%calcul des intégrales par méthode des rectangles
for i=1:n                      
   if ( x(i)>=0) 
      xx=0:pas:x(i);
   else 
      xx=x(i):pas:0;
   end
   Sx(i)=sum(sin(xx.^2)*pas);  	% .^élément par élément
   Cx(i)=sum(cos(xx.^2)*pas);
end
Sx=sqrt(2/pi)*Sx;
Cx=sqrt(2/pi)*Cx;

subplot(2,1,1);
plot(x,Sx,'r',x,Cx,'b');
title('intégrales de Fresnel');
xlabel('x');
text(1,0.2,'Sx');
text(-1,0.7,'Cx');
grid;

subplot(2,1,2);
plot(Cx,Sx);
ylabel('Sx');
xlabel('C(x)');
text(0.1,0.5,'Spirale de Cornu');
grid;