Aide sur programme Pendule

Bonjour, excusez moi si je suis dans la mauvaise section ...

Alors voilà, j'ai à disposition un programme, qui, par Turtle , lance une application d'un pendule harmonique.

Ce que je voudrais savoir faire, c'est pouvoir modifier couleur, et forme du pendule , et y associer un graphique ( y = A Sin (wt - kx), càd une formule traitant d'un pendule "réel" amorti, ainsi que de montrer cette amortissement sur l'animation . Quelqu'un pourrait-il m'aider ?

Voici ce que j'ai déjà : https://mega.co.nz/#!XAx30SJK!Cmvcs0...pjCTDlajwBYULs
(les graphiques sont à modifier car ils n'appartiennent pas au meme programme car j'ai du changer de sujet)
je vous remercie ! :)

Salut,

Le lien que tu donnes conduit à un site qui demande que l'on agrée à des conditions dont on ignore tout. Pas bien ça. Ensuite l'archive rar téléchargée n'est pas lisible.

S'il s'agit d'un script Python, il est préférable de le copier-coller ici entre des balises code. (le # dans la barre d'outils)

Bonjour, désolé pour cette réponse tardive, mais je n'ai pas reçu de mail. Je vous prie de m'excuser pour le lien que j'ai fourni ... Voici les programmes :
a) code lançant une animation de pendule par Turtle

Code:

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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
from math import sin, pi
from time import sleep
from turtle import *
 
GA = 9.80665 # Gravitational Acceleration (meters per second squared)
FORM = 'Time={:6.3f}, Angle={:6.3f}, Speed={:6.3f}'
 
def main():
    length = 9.0            # Of pendulum (meters)
    ngol = - GA / length    # Negative G over L
    total_time = 0.0        # Seconds
    angle = 1.0             # Initial angle of pendulum (radians)
    speed = 0.0             # Initial angular velocity (radians/second)
    time_step = 0.05        # Seconds
    while total_time < 30.0:
        total_time += time_step
        speed += ngol * sin(angle) * time_step
        angle += speed * time_step
        #print(FORM.format(total_time, angle, speed))
        if draw(angle, length): break
        sleep(time_step)
 
def init():
    setup()
    mode('logo')
    radians()
    speed(0)
    hideturtle()
    tracer(False)
    penup()
 
def draw(angle, length):
    if speed() != 0: return True
    clear()
    setheading(angle + pi)
    pensize(max(round(length), 1))
    pendown()
    forward(length * 25)
    penup()
    dot(length * 10)
    home()
    update()
 
if __name__ == '__main__':
    init()
    main()
    bye()

b) Graphiques déjà réalisés à modifier :

Code:

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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
# -*- coding: cp1252 -*-
import matplotlib.pyplot as plt
import pylab
import numpy
 
tabx = []
taby = []
# Avoir le jeu de data
 
while True:
    x = raw_input("Enter X ")
    y = raw_input("Enter Y ")
    tabx.append(float(x))
    taby.append(float(y))
    another = raw_input("Continue ? y/n \n")
    if another != "y":
        break
# le remplir temps que l'utilisateur a pas mis autre chose que y comme rep
 
temp_min = 0.1
temp_max = 20
#axes horizontaux
 
a,b = pylab.polyfit(tabx, taby, 1)
# calcul du moindre carre
 
def droite_moindre_carre(x):
    return a * x + b
 
x_moindre_carre = numpy.linspace(min(tabx), max(tabx), 100)
#linspace(a,b,n) => fournit une liste de valeurs de a à b, en n étapes
y_moindre_carre = map(droite_moindre_carre, x_moindre_carre)
# calcul des valeurs de la courbe de moidnres carrés
# map => Utilise une fonction sur chaque element d'une liste
 
 
plt.axis([temp_min, temp_max,0,100])
plt.plot(tabx, taby, 'or')
# o => point, r => red
plt.plot(x_moindre_carre, y_moindre_carre, 'b')
# b = blue, pas de o => line
plt.show()
# trace le graphique

Un grand merci à vous ! :)

Up !

Avez-vous vu ma réponse ? ^^