Problème de discrétisation MAtlab

[invitee41990b1]

J'ai la fonction suivante sous Matlab. Ca marche très bien. Mon problème c'est la discrétisation, j'ai besoin qu'il fasse les calculs avec un pas plus petit pour r. En effet, à cause de la boucle "for", j'ai un pas de 1 et ça ne va pas du tout!!!

Pouvez-vous m'aider je vous en supplie.



function intensite=propagation(z)

Rmax=73;
intensite=zeros(1,Rmax);

R=73;
D=2*R;
a=14;
l0=0.351;
n=1.5;
k0=2*pi*n/l0;
zr=k0*D^2/4;

for r=1:Rmax
A=0;
for s=1:50
A=A+1+i*exp(i*(zr/(2*z))*(r^2+1))*(i*(-i)^s*(1/r)^s*bessel(s,zr*r/z))-1/(1+2*a*z/zr)*i*(-i)^s*((1+2*a*z/zr)*1/r)^s*bessel(s,zr*r/z);

end
intensite(r)=A;
end
-----

[erff]

Bonjour,

Dans Scilab, on peut faire la chose suivante :

for i=1:0.1:10 (instructions....) end;


Ce qui veut dire que i va de 1 a 10 par pas de 0.1


Je pense que c'est faisable dans Matlab aussi !!!

[invitee41990b1]

Non ce n'est pas possible. Ça me met un message d'erreur style doit être un "interger"!!!

[erff]

Dans ce cas, si on veut par exemple aller de 0 à 10 par pas de 0.1, il suffit de faire cela

r=0

for i=0:100
(...)
r=i./10 (ainsi on incrémente r de 0.1 à chaque itérations)
(...)

end

[A voir en vidéo sur Futura]

[invite455504f8]

une suggestion avec 100 point de discrétisation (en majuscule)

function intensite=propagation(z)

Rmax=73;
intensite=zeros(1,Rmax);

R=73;
D=2*R;
a=14;
l0=0.351;
n=1.5;
k0=2*pi*n/l0;
zr=k0*D^2/4;


IND=0;
NPOINT=100;
for r=LINSPACE(1,Rmax,NPOINT)
IND=IND+1;
A=0;
for s=1:50
A=A+1+i*exp(i*(zr/(2*z))*(r^2+1))*(i*(-i)^s*(1/r)^s*bessel(s,zr*r/z))-1/(1+2*a*z/zr)*i*(-i)^s*((1+2*a*z/zr)*1/r)^s*bessel(s,zr*r/z);

end
intensite(IND)=A;
end