Discussion :

[serma]

Bonjour, j'essaie de tracer la fonction donnée avec des sommes pour le cas de N = 5

Code :

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

1
2
3
4
5
6
7
syms k phi n
z = 0:0.1:20;
J = zeros(5,201);
for i = 0:4
    J(i+1,:) = besselj(i,z);
end
F1 = symsum(i.^(-n).*exp(i.*n.*phi).*J(5),k,-5,5);

J'obtiens un tas de nombres complexes, tels que:

(82822462903397921*exp(4*n*phi))/(288230376151711744*4^n), -(73243657325038871*exp(4*n*phi))/(144115188075855872*4^n), -(52092056568810841*exp(4*n*phi))/(72057594037927936*4^n), -(4186064380238411*exp(4*n*phi))/(4503599627370496*4^n), -(40571731930933357*exp(4*n*phi))/(36028797018963968*4^n), -


Mais j'aimerais le tracer à la place, mais en intrigue polaire. Comment je fais ça?
Merci

[Invité]

Bonjour,

Code :

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

1
2
3
4
5
6
7
8
9
syms n r phi;
a = j^(-n);
b = besselj(r, phi);
c = exp(j*n*phi);
f = a*b*c;
v = subs(f, n, -5:5);
s = sum(v);
s = simplify(s);
ezplot(s);