#!/usr/bin/env python3
# -*- coding: utf-8 -*-

# =============================================================================
# 1C17 Activité 4                                                      CONFIRME
# Réaliser : Simulation de la propagation d'une onde sinusoïdale
# =============================================================================

import numpy  as np
from matplotlib import pyplot as plt
from matplotlib import animation as anim
from math import pi

# =============================================================================
# Question 2 c
# =============================================================================
"""
Compléter les lignes du code selon les consignes à l'aide de la fiche méthode
disponible sur le site sirius.nathan.fr puis exécuter le programme pour simuler
la propagation des vagues destinées aux surfeurs expérimentés du parc à vagues à
la fréquence maximale.
"""

x = ...A compléter...      # Domaine des abcisses (en m) de 0 à 100m 256 points
v = 4.7                    # Célérité (en m/s)
delta_t = (x[1]-x[0])/v    # Intervalle de temps (en s)

A = ...A compléter...      # Amplitude (en m)
T = ...A compléter...      # Période (en s)

# Initialisation de la figure
figure = plt.figure('Propagation d\'une onde sinsusoïdale')
plt.axis([0,100,-10,10])
plt.xlabel(...A compléter...)
plt.ylabel(...A compléter...)

# Initialisation du graphe tracé mis à jour au fur et à mesure
courbe, = ...A compléter... # Graphe sans point

# Fonction fixant l'arrière de l'animation présent à chaque image
def init():
    ...A compléter...
    return ...A compléter...

# Création de la fonction 'onde' appelée à chaque nouvelle image i
def onde(...): 
    t = i*delta_t
    y = A*np.cos((2*pi/T)*t-(2*pi/(v*T))*x)
    ...A compléter...       # Mise à jour du graphe
    return ...A compléter...

# Animation de la figure affichant onde(i), à chaque image i, pour i
# entier de 0 à 255 avec un délai entre 2 images égal à interval en ms
# blit=True indique que seuls les éléments modifiés sont redessinés
animation = anim.FuncAnimation(...,...,...,frames=...,
                               blit=True,interval=delta_t*1000, repeat=False)

