Fonctions

Fonctions de Bessel de 2ème espèce modifiée.

%-------------------------------------------------------------------
% utilisation fonction besselk
% d'aprés 'Apprendre et maitriser Matlab'
% erreur corrigée
%-----------------------------------------------------------------
% f de bessel d'ordre v
% sol particulière de y"+(1/t)y'+(1-[v/t]²)y=0
% on remplace t par l'imaginaire pur jt
% sol générale
%        y = A I (t) + B K (t)
%               v         v
%   Iv : Bessel 1ère espèce mod  besseli
%   Kv : Bessel 2ème espèce mod  besselk
% avec                           2i
%              v   inf      [t/2]
%    I  = [t/2]   Sigma------------------
%     v            i=0   i! gamma(v+i+1)
%
%            I  (t) - I (t)
%            (-v)      v
%    K  = --------------------
%     v      (2/pi)sin(v pi )
%
%----------------------------------------------------------------
echo on;
%----------------------------------------------------------------
%  fonction Bessel 2ème espèce modifiée
% voir contenu fichier pour explications
%----------------------------------------------------------------
echo off;
clear all;
t=20:0.5:24; 		% domaine
nu=0:5;     		% ordre nu varie de 0 à 5 pour comparaison
b2m=besselk(nu,t'); 	% erreur corrigée t'=transposée de t
	% ici on a inversé d'abord t puis nu
	% et il faut transposer b2m
mesh(t,nu,b2m'); 	% maillage3D
ylabel('ordre');
xlabel('t');
zlabel('Inu(t)');
title('Bessel 2ème espèce modifiée');

% vous pouvez aussi essayer ce qui suit à la place de ce qui précède
%mesh(nu,t,b2m); % maillage3D
%xlabel('ordre');
%ylabel('t');
%zlabel('Inu(t)');
%title('Bessel 2ème espèce modifiée');