def animate(i):
plt.clf()
ti = t[i]
print(ti)
xi = [xp[:i + 1] for xp in x] # On stocke l'avancement
vi = [vp[:i + 1] for vp in v] # jusqu'à l'instant présent
plt.suptitle('Oscillateur de Landau, $t={}$'.format(round(ti, 2)))
plt.subplot(3, 1, 1) # Première sous-figure
plt.plot([t[0], t[-1]], [np.sqrt(ell0**2 - d**2)] * 2, 'k')
plt.ylabel('$x$')
plt.title('Evolution temporelle de $x$ (en fonction de $t$)')
for j in range(len(xi)):
plt.plot(t[:i + 1], xi[j], color=colors[j])
plt.ylim(max(0, min(x[-1])), xi[-1][0])
if i < decalage: plt.xlim(0, t[decalage])
else: plt.xlim(t[i - decalage], t[i])
plt.subplot(3, 1, 2) # Seconde sous-figure
portrait_de_phase(xi,
vi,
fantome=50,
clearfig=False,
color=colors,
ylim=vlim)
plt.xlabel('')
plt.subplot(3, 1, 3) # Troisième sous-figure
diagramme_energetique(xi, vi, Em, color=colors, clearfig=False, fantome=50)
def animate(i):
plt.clf()
ti = t[i]
print(ti)
xi = [xp[:i+1] for xp in x] # On stocke l'avancement
vi = [vp[:i+1] for vp in v] # jusqu'à l'instant présent
plt.suptitle('Oscillateur de Landau, $t={}$'.format(round(ti,2)))
plt.subplot(3,1,1) # Première sous-figure
plt.plot([t[0],t[-1]],[np.sqrt(ell0**2-d**2)]*2,'k')
plt.ylabel('$x$')
plt.title('Evolution temporelle de $x$ (en fonction de $t$)')
for j in range(len(xi)):
plt.plot(t[:i+1],xi[j],color=colors[j])
plt.ylim(max(0,min(x[-1])),xi[-1][0])
if i < decalage: plt.xlim(0,t[decalage])
else: plt.xlim(t[i-decalage],t[i])
plt.subplot(3,1,2) # Seconde sous-figure
portrait_de_phase(xi,vi,fantome=50,clearfig=False,color=colors,ylim=vlim)
plt.xlabel('')
plt.subplot(3,1,3) # Troisième sous-figure
diagramme_energetique(xi,vi,Em,color=colors,clearfig=False,fantome=50)
def animate(i):
plt.clf()
ti = t[i]
print(ti)
xi = [xp[:i+1] for xp in x] # On stocke l'avancement
vi = [vp[:i+1] for vp in v] # jusqu'à l'instant présent
plt.suptitle('Pendule pesant, $t={}$'.format(round(ti,2)))
plt.subplot(2,2,1) # Première sous-figure
plt.ylabel('$x$')
plt.title('Evolution temporelle de $\\theta$ (en fonction de $t$)')
for j in range(len(xi)):
plt.plot(t[:i+1],xi[j],color=colors[j])
xlim = min(x[-1]),xi[-1][0]
plt.ylim(xlim)
if i < decalage: plt.xlim(0,t[decalage])
else: plt.xlim(t[i-decalage],t[i])
plt.subplot(2,2,3) # Seconde sous-figure
portrait_de_phase(xi,vi,fantome=50,clearfig=False,color=colors,ylim=vlim,xlim=xlim)
plt.xlabel('')
plt.subplot(1,2,2) # Troisième sous-figure
diagramme_energetique(xi,vi,Em,color=colors,clearfig=False,fantome=50,xlim=xlim)
for xi,vi in zip(x0,v0): # Itération sur les conditions initiales
print(xi,vi) # Un peu de feedback
sol = sp.integrate.odeint(landau,[xi,vi],t) # On intègre
x.append(sol[:,0]) # et on stocke à la fois les positions
v.append(sol[:,1]) # et les vitesses
fig = plt.figure(figsize=(10,10)) # Création de la figure
vlim = (np.min(v),np.max(v)) # Limites verticales (horizontales imposées par les CI)
base_name='PNG/M4_oscillateur_de_landau_tilte_portrait_de_phase'
for i,ti in enumerate(t): # On regarde à chaque instant d'intégration
print(ti) # Un peu de feedback
xi = [xp[:i+1] for xp in x] # On stocke l'avancement
vi = [vp[:i+1] for vp in v] # jusqu'à l'instant présent
plt.suptitle('Oscillateur de Landau tilte, $t={}$'.format(round(ti,2)))
plt.subplot(2,1,1) # Première sous-figure
portrait_de_phase(xi,vi,fantome=50,clearfig=False,color=colors,ylim=vlim)
plt.xlabel('')
plt.subplot(2,1,2) # Seconde sous-figure
diagramme_energetique(xi,vi,Em,color=colors,clearfig=False,fantome=50)
plt.savefig('{}_{:04d}.png'.format(base_name,i))
plt.clf()
from film import make_film # Boîte à outil visuelle
make_film(base_name) # et fabrication du film correspondant
base_name = 'PNG/M4_pendule_simple_portrait_de_phase_zoom'
for i, ti in enumerate(t): # Affichage progressif
print(ti) # Un peu de feed back
thi = [th_p[:i + 1] for th_p in th] # On ne prend que jusqu'�
thpi = [thp_p[:i + 1] for thp_p in thp] # l'instant présent
plt.suptitle('Pendule simple, $t={}$'.format(round(ti, 2)))
plt.subplot(2, 1, 1) # Sous-figure du haut
portrait_de_phase(thi,
thpi,
fantome=20,
clearfig=False,
color=colors,
xlim=th_lim,
ylim=thp_lim)
plt.xlabel('')
plt.subplot(2, 1, 2) # Sous-figure du bas
diagramme_energetique(thi,
thpi,
Em,
color=colors,
clearfig=False,
fantome=20,
xlim=th_lim)
plt.savefig('{}_{:04d}.png'.format(base_name, i))
plt.clf() # Nettoyage
from film import make_film # Boîte à outil pour faire un film
make_film(base_name) # et appel effectif
fig = plt.figure(figsize=(10, 10)) # Création de la figure
# Limites verticales (horizontales imposées par les CI)
vlim = (np.min(v), np.max(v))
base_name = 'PNG/M4_oscillateur_de_landau_portrait_de_phase'
for i, ti in enumerate(t): # On regarde à chaque instant d'intégration
print(ti) # Un peu de feedback
xi = [xp[:i + 1] for xp in x] # On stocke l'avancement
vi = [vp[:i + 1] for vp in v] # jusqu'à l'instant présent
plt.suptitle('Oscillateur de Landau, $t={}$'.format(round(ti, 2)))
plt.subplot(2, 1, 1) # Première sous-figure
portrait_de_phase(xi,
vi,
fantome=50,
clearfig=False,
color=colors,
ylim=vlim)
plt.xlabel('')
plt.subplot(2, 1, 2) # Seconde sous-figure
diagramme_energetique(xi, vi, Em, color=colors, clearfig=False, fantome=50)
plt.savefig('{}_{:04d}.png'.format(base_name, i))
plt.clf()
from film import make_film # Boîte à outil visuelle
make_film(base_name) # et fabrication du film correspondant
os.system("pause")
sol = sp.integrate.odeint(pendule, [thi, thpi], t) # Intégration
th.append(sol[:, 0]) # Ajout des positions
thp.append(sol[:, 1]) # et vitesses correspondantes
fig = plt.figure(figsize=(10, 10)) # Création de la figure
th_lim = (np.min(th), np.max(th)) # Limites en theta
thp_lim = (np.min(thp), np.max(thp)) # Limites en theta point
base_name = 'PNG/M4_pendule_simple_portrait_de_phase_zoom'
for i, ti in enumerate(t): # Affichage progressif
print(ti) # Un peu de feed back
thi = [th_p[:i + 1] for th_p in th] # On ne prend que jusqu'à
thpi = [thp_p[:i + 1] for thp_p in thp] # l'instant présent
plt.suptitle('Pendule simple, $t={}$'.format(round(ti, 2)))
plt.subplot(2, 1, 1) # Sous-figure du haut
portrait_de_phase(thi, thpi, fantome=20, clearfig=False,
color=colors, xlim=th_lim, ylim=thp_lim)
plt.xlabel('')
plt.subplot(2, 1, 2) # Sous-figure du bas
diagramme_energetique(thi, thpi, Em, color=colors,
clearfig=False, fantome=20, xlim=th_lim)
plt.savefig('{}_{:04d}.png'.format(base_name, i))
plt.clf() # Nettoyage
from film import make_film # Boîte à outil pour faire un film
make_film(base_name) # et appel effectif
os.system("pause")