Discussion :

[Suntory]

Bonjour,

Je cherche a plot un cercle ondule (https://www.mypandakitchen.com/1656-...ond-ondule.jpg) dont je controle la periode d'oscillation.
Est-ce que quelqu'un a une idee s'il vous plait ?
Je vous en remercie

code cercle normal :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
from matplotlib import pyplot as plt
import numpy as np
from math import pi
import random
 
x=np.linspace(-10,10,100)
y=x
 
X, Y = np.meshgrid(x,y)
 
circle = (X)**2 + (Y)**2 - 5.
 
plt.contour(X,Y,circle,[0])
plt.show()

[wiztricks]

Salut,

Est-ce que quelqu'un a une idee s'il vous plait ?

Si le rayon du cercle est une fonction périodique, il sera sans doute plus simple de tracer le cercle en calculant x, y en fonction du rayon plutôt qu'en fonction de l'équation du cercle.

- W

[Suntory]

Mon rayon r est une fonction periodique.
x1 et x2 sont calcules en fonction de r et j'obtiens un joli papillon ...

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
from matplotlib import pyplot as plt
import numpy as np
from math import pi
theta = np.linspace(-pi,pi, 100)
 
r = np.sin(theta*2)
 
x1 = r*np.cos(theta)
x2 = r*np.sin(theta)
 
fig, ax = plt.subplots(1)
ax.plot(x1, x2)
plt.show()

[wiztricks]

Regardez le graphique que vous cherchez à reproduire.

Si la valeur du rayon est 1, il varie plutôt de 5 a 10% autour de celle ci, ce qui donne entre 0.9 et 1.1.
Votre rayon varie entre -1 et +1!

Puis vous voyez aussi que pendant que theta va de -pi a +pi, il y a plein d'oscillations autour du rayon.

Rendre compte de çà, c'est juste des maths à poser avant de coder.

- W

[Suntory]

Bonjour,

Merci beacoup.
Voici le resultats avec le rayon r compris entre 0.9< r <1.1.

tan(theta) = y/x du coup theta = arctan(y/x)

2 version : cartesienne (x,y) et (r,theta)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
from matplotlib import pyplot as plt
import numpy as np
from math import pi
 
x=np.linspace(-10,10,1000)
y=x
 
X, Y = np.meshgrid(x,y)
 
circle = (X)**2 + (Y)**2 - (np.sqrt(1+ (np.sin(np.arctan(Y/X)*10))/10))**2
 
plt.contour(X,Y,circle)
plt.show()
 
 
theta = np.linspace(-pi,pi, 1000)
 
r = 1+ (np.sin(theta*10))/10
 
x1 = r*np.cos(theta)
x2 = r*np.sin(theta)
 
fig, ax = plt.subplots(1)
ax.plot(x1, x2)
plt.show()

[Suntory]

Maintenant si je veux faire la meme chose mais pour une sphere (balle de golf).

x**2 + y**2 + z**2 -R**2 avec R=f(X,Y,Z)=f(periodique)

La je ne sais pas comment exprimer le deuxieme angle phi=f(x,y,z) ?

D'habitude :
theta = arctan(y/x)
phi = arcos(z,R)

Mais ici R n'est plus une constante