Discussion :

[el_fuche]

Bonjour à tous,

Je cherche à tracer cette fonction dont l'équation (tirée d'une publication) est la suivante : cf pièce jointe pour plus de clarté

(kp ²)/(1-kp ²)=

((1-σ )*J_1 [n1 (1+∆f/f_r ) ] - n1*(1+∆f/f_a )*J_0 [n_1 (1+∆f/f_r ) ])
--------------------------------------------------------------------------------
((1+σ ) J_1 [n_1 (1+∆f/f_a ) ] )

s et n1 sont des constantes : pour s=0.31, n1 =2.05.
Soit x = ∆f/f_r variant de 0 à 0.2,
J0 et J1 sont les fonctions de Bessel de première espèce d'ordre 0 et 1,

il faut tracer kp = f(x).

Je n'ai pas de problème particulier pour la programmation de cette fonction mais je n'arrive pas obtenir les bonnes valeurs, par exemple :

kp(0.18) = 0.60
kp(0.09) = 0.45
kp(0.05) = 0.15

Perso, j'ai programmé ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
x=0:0.001:0.2;
k=ones(1,length(x));
for i=1:length(x)
    J0=besselj(0,n*(1+x(i)));
    J1=besselj(1,n*(1+x(i)));
    t=(((1-s)*J1)-(n*(1+x(i))*J0))/((1+s)*(n*(1+x(i))*J1));
    k(i)=sqrt(t/(1+t));
end

Je ne trouve pas ce qui ne va pas.

Merci pour votre aide !