Contact Us!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutSign UpSign In

Real-time collaboration for Jupyter Notebooks, Linux Terminals, LaTeX, VS Code, R IDE, and more,
all in one place. Commercial Alternative to JupyterHub.

| Download
Views: 61
Image: ubuntu2004
Kernel: Python 3 (system-wide)

DM3 : Premières applications scientifiques !!!! 😟

Exercice 1: Troisième loi de Kepler

Dans le tableau en fin d'exercice sont donnés la période de révolution TT et le demi grand-axe aa de l'orbite des planètes du système solaire.

Questions :

PlanèteTT [années]aa [UA]
Mercure0.24080.3871
Vénus0.61520.7233
Terre1.0001.000
Mars1.88081.5237
Jupiter11.8625.2026
Saturne29.4579.5547
Uranus84.02019.218
Neptune164.7730.109
import matplotlib.pyplot as plt plt.rcParams.update({'font.size': 22}) # astuce pour avoir des fontes plus grosses dans le TD import numpy as np T=np.array([0.2408,0.6152,1.000,1.8808,11.862,29.457,84.020,164.77]) a=np.array([0.3871,0.7233,1.000,1.5237,5.2026,9.5547,19.218,30.109])
plt.scatter(a,T,100,c='violet')
<matplotlib.collections.PathCollection at 0x7f0e90a42910>
Image in a Jupyter notebook
for i in np.arange(0,len(T),1): x=T[i]**2/a[i]**3 if x==1: print(x,'vérifié') else: print(x)
0.9996395980610487 1.0001779096233403 1.0 vérifié 0.9999696744677578 0.9992048412523753 0.9947772916493796 0.994582440112815 0.994643132447185

Exercice 2: calcul et tracé de la gravité terrestre en fonction de la distance r au centre de la Terre

La gravité créée par la planète Terre supposée sphérique, immobile et de composition homogène peut se formuler en fonction de la distance r au centre de la planète : g(r)=GM(r)r2,g(r)=-\frac{GM(r)}{r^2},

M(r)M(r) est la masse située dans la sphère de rayon rr :

  • Si on est à l'extérieur de la Terre : M(r)=MTerreM(r)=M_{Terre}

  • Si on est à l'intérieur de la Terre : M(r)=ρV(r)=MTerreVTerreV(r)=MTerrer3RTerre3M(r)=\rho V(r)=\frac{M_{Terre}}{V_{Terre}}V(r)=\frac{M_{Terre}r^3}{R_{Terre}^3}

On donne la constante de gravitation universelle G=6.671011G=6.67\cdot10^{-11} (SI) la masse de la Terre MTerre=5.9721024M_{Terre}=5.972\cdot10^{24} kg le rayon moyen de la Terre RTerre=6371R_{Terre}=6371 km

Questions : -Calculez la masse M(r) et la gravité g(r) pour r variant de 0 à 10 000 km par pas de 100 -Tracez la masse M(r) et la gravité g(r)

G=6.67*(10**-11) Mt=5.972*10**24 Rt=6371
listeMr = [] listegr = [] Mr=0 gr=0 for r in np.arange(0,10100,100): if r>6371: Mr=Mt gr=(-G*Mr)/((r*1000)**2) listeMr.append(Mr) listegr.append(gr) else: Mr=(Mt*r**3)/(Rt**3) gr=(-G*Mr)/((r*1000)**2) listeMr.append(Mr) listegr.append(gr)
/tmp/ipykernel_3160/3939560052.py:13: RuntimeWarning: invalid value encountered in double_scalars gr=(-G*Mr)/((r*1000)**2)
[nan, -0.15403620761836861, -0.30807241523673723, -0.46210862285510584, -0.6161448304734745, -0.770181038091843, -0.9242172457102117, -1.0782534533285804, -1.232289660946949, -1.3863258685653175, -1.540362076183686, -1.6943982838020548, -1.8484344914204234, -2.002470699038792, -2.1565069066571607, -2.310543114275529, -2.464579321893898, -2.618615529512266, -2.772651737130635, -2.926687944749004, -3.080724152367372, -3.2347603599857413, -3.3887965676041096, -3.5428327752224784, -3.6968689828408468, -3.8509051904592155, -4.004941398077584, -4.158977605695952, -4.313013813314321, -4.46705002093269, -4.621086228551058, -4.775122436169427, -4.929158643787796, -5.083194851406164, -5.237231059024532, -5.391267266642902, -5.54530347426127, -5.699339681879638, -5.853375889498008, -6.007412097116375, -6.161448304734744, -6.3154845123531125, -6.469520719971483, -6.623556927589851, -6.777593135208219, -6.931629342826588, -7.085665550444957, -7.239701758063325, -7.3937379656816935, -7.547774173300063, -7.701810380918431, -7.855846588536799, -8.009882796155168, -8.163919003773538, -8.317955211391904, -8.471991419010275, -8.626027626628643, -8.780063834247011, -8.93410004186538, -9.08813624948375, -9.242172457102116, -9.396208664720485, -9.550244872338855, -9.704281079957221, -9.724912109375, -9.427985798816568, -9.1444536271809, -8.873521942526175, -8.614455017301038, -8.366570048309178, -8.129232653061225, -7.901852806982742, -7.683881172839506, -7.474805779696003, -7.274149013878744, -7.081464888888889, -6.896336565096953, -6.718374093439029, -6.547212360289284, -6.382509213267104, -6.22394375, -6.071214753848499, -5.924039262343843, -5.7821512556249095, -5.645300453514739, -5.513251211072665, -5.385781503515414, -5.262681992337165, -5.143755165289257, -5.02881454361823, -4.917683950617284, -4.810196836130902, -4.706195652173913, -4.605531275291941, -4.508062471706655, -4.41365540166205, -4.322183159722222, -4.233525348070996, -4.147567680133278, -4.0642016120804, -3.983324]
r=np.arange(0,10100,100) np.shape(r)
(101,)
plt.subplot(221) plt.plot(r,listeMr) plt.title('masse') plt.subplot(222) plt.plot(r,listegr) plt.title('gravité') plt.subplots_adjust(wspace=0.5, hspace=0.5)
Image in a Jupyter notebook

Exercice 3 : Lecture de données existantes \to Insolation et δ18\delta^{18}O

Pour finir le DM dans la joie, abordons des données typiques des études paléoclimatiques. On s'intéresse à l'indice δ18\delta^{18}O relevé en Antarctique dans des carottes de glace dont les âges vont de 1000 à 248000 ans. Le δ18\delta^{18}O caractérise la teneur en isotope 18^{18}O de l'Oxygène relativement à l'atome standard plus léger 16^{16}O. On va le comparer à des valeurs d'insolation. (Les fichiers sont dans le dossier Data). Dans cet exercice, nous lisons un jeu de données existant.

Commençons par lire les données 18^{18}O de l'Oxygène

donnees_dO = np.loadtxt('./donnees_dO18.dat') donnees_dO
array([[ 1.00000000e+03, -3.51213643e+01], [ 2.00000000e+03, -3.42391180e+01], [ 3.00000000e+03, -3.48951069e+01], [ 4.00000000e+03, -3.49413222e+01], [ 5.00000000e+03, -3.45101193e+01], [ 6.00000000e+03, -3.37657110e+01], [ 7.00000000e+03, -3.43370954e+01], [ 8.00000000e+03, -3.45240044e+01], [ 9.00000000e+03, -3.56496240e+01], [ 1.00000000e+04, -3.60200000e+01], [ 1.10000000e+04, -3.74626250e+01], [ 1.20000000e+04, -4.07313081e+01], [ 1.30000000e+04, -3.98893683e+01], [ 1.40000000e+04, -3.69386731e+01], [ 1.50000000e+04, -3.92323404e+01], [ 1.60000000e+04, -4.04464245e+01], [ 1.70000000e+04, -4.02547235e+01], [ 1.80000000e+04, -3.97228160e+01], [ 1.90000000e+04, -3.95545917e+01], [ 2.00000000e+04, -4.14665505e+01], [ 2.10000000e+04, -4.12683093e+01], [ 2.20000000e+04, -4.38326282e+01], [ 2.30000000e+04, -4.17429825e+01], [ 2.40000000e+04, -4.10700000e+01], [ 2.50000000e+04, -3.99067976e+01], [ 2.60000000e+04, -4.16300000e+01], [ 2.70000000e+04, -4.16104436e+01], [ 2.80000000e+04, -4.14263235e+01], [ 2.90000000e+04, -4.21549568e+01], [ 3.00000000e+04, -3.99663368e+01], [ 3.10000000e+04, -3.86037539e+01], [ 3.20000000e+04, -4.07803782e+01], [ 3.30000000e+04, -4.12541349e+01], [ 3.40000000e+04, -3.93716607e+01], [ 3.50000000e+04, -3.82705418e+01], [ 3.60000000e+04, -4.10067739e+01], [ 3.70000000e+04, -4.20842559e+01], [ 3.80000000e+04, -4.09272845e+01], [ 3.90000000e+04, -4.04070686e+01], [ 4.00000000e+04, -4.01253978e+01], [ 4.10000000e+04, -4.24562049e+01], [ 4.20000000e+04, -3.84373934e+01], [ 4.30000000e+04, -3.80975648e+01], [ 4.40000000e+04, -3.61200000e+01], [ 4.50000000e+04, -4.04591051e+01], [ 4.60000000e+04, -3.96484356e+01], [ 4.70000000e+04, -3.97517073e+01], [ 4.80000000e+04, -3.94703617e+01], [ 4.90000000e+04, -3.91473242e+01], [ 5.00000000e+04, -3.81862078e+01], [ 5.10000000e+04, -3.76933835e+01], [ 5.20000000e+04, -4.10445518e+01], [ 5.30000000e+04, -4.07020200e+01], [ 5.40000000e+04, -4.17362080e+01], [ 5.50000000e+04, -3.84423418e+01], [ 5.60000000e+04, -3.98486213e+01], [ 5.70000000e+04, -4.04953853e+01], [ 5.80000000e+04, -4.00778705e+01], [ 5.90000000e+04, -4.08285217e+01], [ 6.00000000e+04, -4.11335815e+01], [ 6.10000000e+04, -4.04768911e+01], [ 6.20000000e+04, -4.18090322e+01], [ 6.30000000e+04, -4.22671833e+01], [ 6.40000000e+04, -4.17727669e+01], [ 6.50000000e+04, -4.20869449e+01], [ 6.60000000e+04, -4.14914853e+01], [ 6.70000000e+04, -3.97167194e+01], [ 6.80000000e+04, -3.93606372e+01], [ 6.90000000e+04, -3.80284771e+01], [ 7.00000000e+04, -4.25464030e+01], [ 7.10000000e+04, -4.27306289e+01], [ 7.20000000e+04, -3.90576044e+01], [ 7.30000000e+04, -3.76922106e+01], [ 7.40000000e+04, -3.93984494e+01], [ 7.50000000e+04, -4.18933780e+01], [ 7.60000000e+04, -3.96634585e+01], [ 7.70000000e+04, -3.79273442e+01], [ 7.80000000e+04, -3.84795624e+01], [ 7.90000000e+04, -3.77839862e+01], [ 8.00000000e+04, -3.72039901e+01], [ 8.10000000e+04, -3.74147322e+01], [ 8.20000000e+04, -3.65000000e+01], [ 8.30000000e+04, -3.59844384e+01], [ 8.40000000e+04, -3.90528622e+01], [ 8.50000000e+04, -3.99004954e+01], [ 8.60000000e+04, -4.03186657e+01], [ 8.70000000e+04, -3.79827842e+01], [ 8.80000000e+04, -3.77403068e+01], [ 8.90000000e+04, -3.89395136e+01], [ 9.00000000e+04, -3.94023460e+01], [ 9.10000000e+04, -3.96989227e+01], [ 9.20000000e+04, -3.94634157e+01], [ 9.30000000e+04, -3.91951444e+01], [ 9.40000000e+04, -3.92337999e+01], [ 9.50000000e+04, -3.88551453e+01], [ 9.60000000e+04, -3.94804313e+01], [ 9.70000000e+04, -3.86773983e+01], [ 9.80000000e+04, -3.80800922e+01], [ 9.90000000e+04, -3.78258667e+01], [ 1.00000000e+05, -3.73915342e+01], [ 1.01000000e+05, -3.67759085e+01], [ 1.02000000e+05, -3.63851271e+01], [ 1.03000000e+05, -3.67638834e+01], [ 1.04000000e+05, -3.90595507e+01], [ 1.05000000e+05, -3.89084973e+01], [ 1.06000000e+05, -3.67934053e+01], [ 1.07000000e+05, -3.77840679e+01], [ 1.08000000e+05, -3.82439274e+01], [ 1.09000000e+05, -3.83753699e+01], [ 1.10000000e+05, -3.87461895e+01], [ 1.11000000e+05, -3.93538433e+01], [ 1.12000000e+05, -3.90595324e+01], [ 1.13000000e+05, -3.76866156e+01], [ 1.14000000e+05, -3.42465242e+01], [ 1.15000000e+05, -3.41520549e+01], [ 1.16000000e+05, -3.47313739e+01], [ 1.17000000e+05, -3.51851646e+01], [ 1.18000000e+05, -3.79641960e+01], [ 1.19000000e+05, -3.35140120e+01], [ 1.20000000e+05, -3.51046799e+01], [ 1.21000000e+05, -3.53312521e+01], [ 1.22000000e+05, -3.66121380e+01], [ 1.23000000e+05, -3.65946406e+01], [ 1.24000000e+05, -3.69258806e+01], [ 1.25000000e+05, -3.70222270e+01], [ 1.26000000e+05, -3.77745110e+01], [ 1.27000000e+05, -3.16888517e+01], [ 1.28000000e+05, -3.18957082e+01], [ 1.29000000e+05, -3.20822895e+01], [ 1.30000000e+05, -3.19960505e+01], [ 1.31000000e+05, -3.35401997e+01], [ 1.32000000e+05, -3.19290821e+01], [ 1.33000000e+05, -3.66826373e+01], [ 1.34000000e+05, -3.66719353e+01], [ 1.35000000e+05, -3.60502196e+01], [ 1.36000000e+05, -3.78121029e+01], [ 1.37000000e+05, -3.97637781e+01], [ 1.38000000e+05, -3.50546847e+01], [ 1.39000000e+05, -3.86091840e+01], [ 1.40000000e+05, -3.90062163e+01], [ 1.41000000e+05, -3.96522426e+01], [ 1.42000000e+05, -4.00507925e+01], [ 1.43000000e+05, -3.92059381e+01], [ 1.44000000e+05, -3.93107795e+01], [ 1.45000000e+05, -4.11842493e+01], [ 1.46000000e+05, -4.01600261e+01], [ 1.47000000e+05, -4.09860373e+01], [ 1.48000000e+05, -4.02863748e+01], [ 1.49000000e+05, -4.04214788e+01], [ 1.50000000e+05, -4.03169497e+01], [ 1.51000000e+05, -3.98602459e+01], [ 1.52000000e+05, -3.95405688e+01], [ 1.53000000e+05, -3.96986612e+01], [ 1.54000000e+05, -3.66489202e+01], [ 1.55000000e+05, -4.02722625e+01], [ 1.56000000e+05, -3.91785652e+01], [ 1.57000000e+05, -4.05689027e+01], [ 1.58000000e+05, -3.99292972e+01], [ 1.59000000e+05, -3.97437608e+01], [ 1.60000000e+05, -4.00387672e+01], [ 1.61000000e+05, -4.00829948e+01], [ 1.62000000e+05, -3.85966730e+01], [ 1.63000000e+05, -3.96557780e+01], [ 1.64000000e+05, -3.90715278e+01], [ 1.65000000e+05, -3.78276767e+01], [ 1.66000000e+05, -3.85220097e+01], [ 1.67000000e+05, -3.70730892e+01], [ 1.68000000e+05, -3.79103636e+01], [ 1.69000000e+05, -3.80786901e+01], [ 1.70000000e+05, -3.79168094e+01], [ 1.71000000e+05, -3.87222721e+01], [ 1.72000000e+05, -3.26128661e+01], [ 1.73000000e+05, -3.67909323e+01], [ 1.74000000e+05, -3.86241511e+01], [ 1.75000000e+05, -3.82900819e+01], [ 1.76000000e+05, -3.93093771e+01], [ 1.77000000e+05, -3.85166564e+01], [ 1.78000000e+05, -3.96541818e+01], [ 1.79000000e+05, -3.89516712e+01], [ 1.80000000e+05, -4.03134791e+01], [ 1.81000000e+05, -3.78449125e+01], [ 1.82000000e+05, -3.88254228e+01], [ 1.83000000e+05, -3.61453777e+01], [ 1.84000000e+05, -3.42864973e+01], [ 1.85000000e+05, -3.73682331e+01], [ 1.86000000e+05, -3.69530647e+01], [ 1.87000000e+05, -3.67253258e+01], [ 1.88000000e+05, -3.75080872e+01], [ 1.89000000e+05, -3.74457534e+01], [ 1.90000000e+05, -3.74185078e+01], [ 1.91000000e+05, -3.83464475e+01], [ 1.92000000e+05, -3.84283995e+01], [ 1.93000000e+05, -3.73464071e+01], [ 1.94000000e+05, -3.72769013e+01], [ 1.95000000e+05, -3.72108380e+01], [ 1.96000000e+05, -3.65831062e+01], [ 1.97000000e+05, -3.66496801e+01], [ 1.98000000e+05, -3.83045981e+01], [ 1.99000000e+05, -3.71950516e+01], [ 2.00000000e+05, -3.68581502e+01], [ 2.01000000e+05, -3.72867825e+01], [ 2.02000000e+05, -3.82469753e+01], [ 2.03000000e+05, -3.70441237e+01], [ 2.04000000e+05, -3.75352357e+01], [ 2.05000000e+05, -3.94008253e+01], [ 2.06000000e+05, -3.78673074e+01], [ 2.07000000e+05, -3.80594857e+01], [ 2.08000000e+05, -3.72407814e+01], [ 2.09000000e+05, -3.74980720e+01], [ 2.10000000e+05, -3.71289129e+01], [ 2.11000000e+05, -3.61585056e+01], [ 2.12000000e+05, -3.66858189e+01], [ 2.13000000e+05, -3.59393873e+01], [ 2.14000000e+05, -3.67645145e+01], [ 2.15000000e+05, -3.52433829e+01], [ 2.16000000e+05, -3.66384321e+01], [ 2.17000000e+05, -3.55579111e+01], [ 2.18000000e+05, -3.61983863e+01], [ 2.19000000e+05, -3.59851599e+01], [ 2.20000000e+05, -3.50658452e+01], [ 2.21000000e+05, -3.64142079e+01], [ 2.22000000e+05, -3.64529976e+01], [ 2.23000000e+05, -3.51683544e+01], [ 2.24000000e+05, -3.61901565e+01], [ 2.25000000e+05, -3.56084996e+01], [ 2.26000000e+05, -3.69936001e+01], [ 2.27000000e+05, -3.71176021e+01], [ 2.28000000e+05, -3.68476934e+01], [ 2.29000000e+05, -3.67283679e+01], [ 2.30000000e+05, -3.69207554e+01], [ 2.31000000e+05, -3.68893214e+01], [ 2.32000000e+05, -3.67302248e+01], [ 2.33000000e+05, -3.66640244e+01], [ 2.34000000e+05, -3.68130456e+01], [ 2.35000000e+05, -3.69782424e+01], [ 2.36000000e+05, -3.63014765e+01], [ 2.37000000e+05, -3.62471597e+01], [ 2.38000000e+05, -3.65987539e+01], [ 2.39000000e+05, -3.68926689e+01], [ 2.40000000e+05, -3.70851344e+01], [ 2.41000000e+05, -3.71978927e+01], [ 2.42000000e+05, -3.77159586e+01], [ 2.43000000e+05, -3.77803696e+01], [ 2.44000000e+05, -3.73578496e+01], [ 2.45000000e+05, -3.77506626e+01], [ 2.46000000e+05, -3.81433636e+01], [ 2.47000000e+05, -3.79070987e+01], [ 2.48000000e+05, -3.73526722e+01]])

La première colonne de donnees_d0, contient les années, la deuxième colonne de donnees_d0, contient les valeurs de δ18O\delta^{18}O

année=np.loadtxt("donnees_dO18.dat",usecols=0) année
array([ 1000., 2000., 3000., 4000., 5000., 6000., 7000., 8000., 9000., 10000., 11000., 12000., 13000., 14000., 15000., 16000., 17000., 18000., 19000., 20000., 21000., 22000., 23000., 24000., 25000., 26000., 27000., 28000., 29000., 30000., 31000., 32000., 33000., 34000., 35000., 36000., 37000., 38000., 39000., 40000., 41000., 42000., 43000., 44000., 45000., 46000., 47000., 48000., 49000., 50000., 51000., 52000., 53000., 54000., 55000., 56000., 57000., 58000., 59000., 60000., 61000., 62000., 63000., 64000., 65000., 66000., 67000., 68000., 69000., 70000., 71000., 72000., 73000., 74000., 75000., 76000., 77000., 78000., 79000., 80000., 81000., 82000., 83000., 84000., 85000., 86000., 87000., 88000., 89000., 90000., 91000., 92000., 93000., 94000., 95000., 96000., 97000., 98000., 99000., 100000., 101000., 102000., 103000., 104000., 105000., 106000., 107000., 108000., 109000., 110000., 111000., 112000., 113000., 114000., 115000., 116000., 117000., 118000., 119000., 120000., 121000., 122000., 123000., 124000., 125000., 126000., 127000., 128000., 129000., 130000., 131000., 132000., 133000., 134000., 135000., 136000., 137000., 138000., 139000., 140000., 141000., 142000., 143000., 144000., 145000., 146000., 147000., 148000., 149000., 150000., 151000., 152000., 153000., 154000., 155000., 156000., 157000., 158000., 159000., 160000., 161000., 162000., 163000., 164000., 165000., 166000., 167000., 168000., 169000., 170000., 171000., 172000., 173000., 174000., 175000., 176000., 177000., 178000., 179000., 180000., 181000., 182000., 183000., 184000., 185000., 186000., 187000., 188000., 189000., 190000., 191000., 192000., 193000., 194000., 195000., 196000., 197000., 198000., 199000., 200000., 201000., 202000., 203000., 204000., 205000., 206000., 207000., 208000., 209000., 210000., 211000., 212000., 213000., 214000., 215000., 216000., 217000., 218000., 219000., 220000., 221000., 222000., 223000., 224000., 225000., 226000., 227000., 228000., 229000., 230000., 231000., 232000., 233000., 234000., 235000., 236000., 237000., 238000., 239000., 240000., 241000., 242000., 243000., 244000., 245000., 246000., 247000., 248000.])
O=np.loadtxt("donnees_dO18.dat",usecols=1) O
array([-35.12136432, -34.23911798, -34.89510686, -34.94132225, -34.51011934, -33.76571103, -34.33709542, -34.52400442, -35.64962404, -36.02 , -37.46262501, -40.73130806, -39.88936835, -36.93867307, -39.23234035, -40.44642449, -40.25472349, -39.72281597, -39.55459171, -41.46655046, -41.26830934, -43.83262819, -41.74298247, -41.07 , -39.9067976 , -41.63 , -41.61044357, -41.42632353, -42.1549568 , -39.96633676, -38.60375389, -40.78037823, -41.25413487, -39.37166066, -38.27054179, -41.00677388, -42.08425589, -40.9272845 , -40.40706859, -40.12539779, -42.45620486, -38.43739342, -38.09756482, -36.12 , -40.45910509, -39.64843556, -39.75170732, -39.47036173, -39.14732423, -38.18620783, -37.69338349, -41.04455184, -40.70201998, -41.73620798, -38.44234182, -39.84862134, -40.49538527, -40.07787046, -40.82852173, -41.13358146, -40.47689109, -41.80903218, -42.26718325, -41.77276692, -42.08694491, -41.49148528, -39.71671941, -39.36063716, -38.02847713, -42.54640297, -42.7306289 , -39.05760437, -37.69221063, -39.39844942, -41.89337802, -39.6634585 , -37.92734417, -38.47956236, -37.78398617, -37.20399005, -37.41473224, -36.5 , -35.98443842, -39.05286219, -39.90049537, -40.31866569, -37.98278418, -37.74030675, -38.93951361, -39.40234605, -39.69892272, -39.4634157 , -39.19514437, -39.23379995, -38.85514527, -39.48043126, -38.67739825, -38.08009225, -37.82586673, -37.39153419, -36.77590847, -36.38512711, -36.76388336, -39.05955068, -38.90849727, -36.79340534, -37.78406789, -38.24392739, -38.3753699 , -38.74618955, -39.35384331, -39.05953238, -37.68661564, -34.2465242 , -34.15205494, -34.73137388, -35.18516462, -37.96419596, -33.51401204, -35.10467993, -35.33125214, -36.61213796, -36.59464056, -36.92588056, -37.02222703, -37.77451104, -31.68885169, -31.89570818, -32.08228948, -31.99605054, -33.54019974, -31.92908205, -36.68263728, -36.67193528, -36.05021961, -37.81210288, -39.76377807, -35.0546847 , -38.60918396, -39.00621634, -39.65224256, -40.05079254, -39.20593808, -39.31077954, -41.18424931, -40.16002606, -40.98603726, -40.28637484, -40.4214788 , -40.31694966, -39.86024592, -39.54056876, -39.69866123, -36.64892017, -40.27226251, -39.17856524, -40.56890273, -39.92929724, -39.74376084, -40.03876724, -40.08299483, -38.59667297, -39.65577803, -39.07152782, -37.82767674, -38.52200969, -37.07308924, -37.91036362, -38.07869008, -37.91680939, -38.72227206, -32.61286607, -36.79093226, -38.62415107, -38.29008192, -39.30937707, -38.51665637, -39.65418177, -38.95167119, -40.31347912, -37.84491249, -38.82542282, -36.14537768, -34.28649735, -37.36823308, -36.9530647 , -36.72532582, -37.50808715, -37.44575344, -37.41850783, -38.34644747, -38.42839954, -37.34640706, -37.27690127, -37.21083802, -36.58310622, -36.64968006, -38.30459805, -37.19505162, -36.85815024, -37.28678251, -38.2469753 , -37.04412374, -37.53523574, -39.40082535, -37.86730741, -38.05948568, -37.24078141, -37.49807195, -37.12891294, -36.15850555, -36.68581894, -35.93938728, -36.76451447, -35.24338291, -36.63843209, -35.55791113, -36.19838634, -35.98515989, -35.06584523, -36.41420786, -36.45299764, -35.16835444, -36.19015655, -35.60849962, -36.9936001 , -37.11760205, -36.84769344, -36.72836792, -36.92075542, -36.88932143, -36.73022481, -36.66402438, -36.81304561, -36.97824239, -36.3014765 , -36.24715965, -36.59875393, -36.89266889, -37.08513439, -37.19789266, -37.71595863, -37.78036962, -37.35784956, -37.75066258, -38.14336365, -37.90709868, -37.35267223])

Question : Quelle est la plage de temps étudiée ? Année minimale, année maximale

print('la plage est entre :',min(année),'ans','à',max(année),'ans')
la plage est entre : 1000.0 ans à 248000.0 ans

Question : Faites un graphique des δ18O\delta^{18}O en fonction des années

plt.plot(année,O) plt.title('graphique $\delta^{18}O$ en fonction des années') plt.ylabel('$\delta^{18}O$') plt.xlabel('Années en ans')
Text(0.5, 0, 'Années en ans')
Image in a Jupyter notebook

Question : Faites de même pour les données d'insolation

donnees_insolation= np.loadtxt('./donnees_insolation.dat') donnees_insolation
array([[ 1.0000000e+03, 9.3281100e-01], [ 2.0000000e+03, 3.2769570e+00], [ 3.0000000e+03, 7.0075330e+00], [ 4.0000000e+03, 1.1962172e+01], [ 5.0000000e+03, 1.7848958e+01], [ 6.0000000e+03, 2.4265751e+01], [ 7.0000000e+03, 3.0730781e+01], [ 8.0000000e+03, 3.6723055e+01], [ 9.0000000e+03, 4.1729583e+01], [ 1.0000000e+04, 4.5295127e+01], [ 1.1000000e+04, 4.7068648e+01], [ 1.2000000e+04, 4.6840655e+01], [ 1.3000000e+04, 4.4566351e+01], [ 1.4000000e+04, 4.0371120e+01], [ 1.5000000e+04, 3.4537838e+01], [ 1.6000000e+04, 2.7477815e+01], [ 1.7000000e+04, 1.9689633e+01], [ 1.8000000e+04, 1.1711510e+01], [ 1.9000000e+04, 4.0730830e+00], [ 2.0000000e+04, -2.7481760e+00], [ 2.1000000e+04, -8.3621520e+00], [ 2.2000000e+04, -1.2490773e+01], [ 2.3000000e+04, -1.4981564e+01], [ 2.4000000e+04, -1.5810417e+01], [ 2.5000000e+04, -1.5074452e+01], [ 2.6000000e+04, -1.2976025e+01], [ 2.7000000e+04, -9.7991600e+00], [ 2.8000000e+04, -5.8800800e+00], [ 2.9000000e+04, -1.5742360e+00], [ 3.0000000e+04, 2.7770250e+00], [ 3.1000000e+04, 6.8765220e+00], [ 3.2000000e+04, 1.0494468e+01], [ 3.3000000e+04, 1.3482176e+01], [ 3.4000000e+04, 1.5774981e+01], [ 3.5000000e+04, 1.7384905e+01], [ 3.6000000e+04, 1.8385185e+01], [ 3.7000000e+04, 1.8889915e+01], [ 3.8000000e+04, 1.9032421e+01], [ 3.9000000e+04, 1.8945647e+01], [ 4.0000000e+04, 1.8747018e+01], [ 4.1000000e+04, 1.8529075e+01], [ 4.2000000e+04, 1.8356124e+01], [ 4.3000000e+04, 1.8266137e+01], [ 4.4000000e+04, 1.8276541e+01], [ 4.5000000e+04, 1.8392081e+01], [ 4.6000000e+04, 1.8613076e+01], [ 4.7000000e+04, 1.8942238e+01], [ 4.8000000e+04, 1.9388852e+01], [ 4.9000000e+04, 1.9969419e+01], [ 5.0000000e+04, 2.0704481e+01], [ 5.1000000e+04, 2.1611976e+01], [ 5.2000000e+04, 2.2698076e+01], [ 5.3000000e+04, 2.3946962e+01], [ 5.4000000e+04, 2.5311289e+01], [ 5.5000000e+04, 2.6705217e+01], [ 5.6000000e+04, 2.8001646e+01], [ 5.7000000e+04, 2.9035290e+01], [ 5.8000000e+04, 2.9612374e+01], [ 5.9000000e+04, 2.9527145e+01], [ 6.0000000e+04, 2.8584286e+01], [ 6.1000000e+04, 2.6625132e+01], [ 6.2000000e+04, 2.3554535e+01], [ 6.3000000e+04, 1.9364473e+01], [ 6.4000000e+04, 1.4150408e+01], [ 6.5000000e+04, 8.1171050e+00], [ 6.6000000e+04, 1.5720770e+00], [ 6.7000000e+04, -5.0932880e+00], [ 6.8000000e+04, -1.1432886e+01], [ 6.9000000e+04, -1.6981845e+01], [ 7.0000000e+04, -2.1295142e+01], [ 7.1000000e+04, -2.3984956e+01], [ 7.2000000e+04, -2.4753789e+01], [ 7.3000000e+04, -2.3421449e+01], [ 7.4000000e+04, -1.9944873e+01], [ 7.5000000e+04, -1.4430034e+01], [ 7.6000000e+04, -7.1349440e+00], [ 7.7000000e+04, 1.5377040e+00], [ 7.8000000e+04, 1.1059873e+01], [ 7.9000000e+04, 2.0810521e+01], [ 8.0000000e+04, 3.0120493e+01], [ 8.1000000e+04, 3.8327281e+01], [ 8.2000000e+04, 4.4834616e+01], [ 8.3000000e+04, 4.9169996e+01], [ 8.4000000e+04, 5.1032574e+01], [ 8.5000000e+04, 5.0324677e+01], [ 8.6000000e+04, 4.7162420e+01], [ 8.7000000e+04, 4.1864095e+01], [ 8.8000000e+04, 3.4918377e+01], [ 8.9000000e+04, 2.6937266e+01], [ 9.0000000e+04, 1.8600310e+01], [ 9.1000000e+04, 1.0597132e+01], [ 9.2000000e+04, 3.5739170e+00], [ 9.3000000e+04, -1.9121320e+00], [ 9.4000000e+04, -5.4287730e+00], [ 9.5000000e+04, -6.6953380e+00], [ 9.6000000e+04, -5.6011570e+00], [ 9.7000000e+04, -2.2167370e+00], [ 9.8000000e+04, 3.2023110e+00], [ 9.9000000e+04, 1.0220567e+01], [ 1.0000000e+05, 1.8240342e+01], [ 1.0100000e+05, 2.6535026e+01], [ 1.0200000e+05, 3.4300476e+01], [ 1.0300000e+05, 4.0723734e+01], [ 1.0400000e+05, 4.5064138e+01], [ 1.0500000e+05, 4.6737592e+01], [ 1.0600000e+05, 4.5391607e+01], [ 1.0700000e+05, 4.0958166e+01], [ 1.0800000e+05, 3.3673972e+01], [ 1.0900000e+05, 2.4063086e+01], [ 1.1000000e+05, 1.2884015e+01], [ 1.1100000e+05, 1.0500520e+00], [ 1.1200000e+05, -1.0463809e+01], [ 1.1300000e+05, -2.0713079e+01], [ 1.1400000e+05, -2.8863359e+01], [ 1.1500000e+05, -3.4250634e+01], [ 1.1600000e+05, -3.6421145e+01], [ 1.1700000e+05, -3.5153753e+01], [ 1.1800000e+05, -3.0469128e+01], [ 1.1900000e+05, -2.2628698e+01], [ 1.2000000e+05, -1.2123412e+01], [ 1.2100000e+05, 3.5041600e-01], [ 1.2200000e+05, 1.3931837e+01], [ 1.2300000e+05, 2.7654687e+01], [ 1.2400000e+05, 4.0524006e+01], [ 1.2500000e+05, 5.1603032e+01], [ 1.2600000e+05, 6.0100429e+01], [ 1.2700000e+05, 6.5445507e+01], [ 1.2800000e+05, 6.7340296e+01], [ 1.2900000e+05, 6.5780866e+01], [ 1.3000000e+05, 6.1045514e+01], [ 1.3100000e+05, 5.3653118e+01], [ 1.3200000e+05, 4.4299097e+01], [ 1.3300000e+05, 3.3779427e+01], [ 1.3400000e+05, 2.2913345e+01], [ 1.3500000e+05, 1.2473923e+01], [ 1.3600000e+05, 3.1325270e+00], [ 1.3700000e+05, -4.5799000e+00], [ 1.3800000e+05, -1.0295168e+01], [ 1.3900000e+05, -1.3815920e+01], [ 1.4000000e+05, -1.5112329e+01], [ 1.4100000e+05, -1.4309865e+01], [ 1.4200000e+05, -1.1669709e+01], [ 1.4300000e+05, -7.5625680e+00], [ 1.4400000e+05, -2.4364580e+00], [ 1.4500000e+05, 3.2204090e+00], [ 1.4600000e+05, 8.9194170e+00], [ 1.4700000e+05, 1.4210571e+01], [ 1.4800000e+05, 1.8716638e+01], [ 1.4900000e+05, 2.2158416e+01], [ 1.5000000e+05, 2.4368362e+01], [ 1.5100000e+05, 2.5291620e+01], [ 1.5200000e+05, 2.4975669e+01], [ 1.5300000e+05, 2.3551518e+01], [ 1.5400000e+05, 2.1210426e+01], [ 1.5500000e+05, 1.8180215e+01], [ 1.5600000e+05, 1.4704516e+01], [ 1.5700000e+05, 1.1027033e+01], [ 1.5800000e+05, 7.3814780e+00], [ 1.5900000e+05, 3.9865790e+00], [ 1.6000000e+05, 1.0446260e+00], [ 1.6100000e+05, -1.2584160e+00], [ 1.6200000e+05, -2.7531910e+00], [ 1.6300000e+05, -3.2892610e+00], [ 1.6400000e+05, -2.7404850e+00], [ 1.6500000e+05, -1.0138400e+00], [ 1.6600000e+05, 1.9386480e+00], [ 1.6700000e+05, 6.1058290e+00], [ 1.6800000e+05, 1.1402794e+01], [ 1.6900000e+05, 1.7657812e+01], [ 1.7000000e+05, 2.4603454e+01], [ 1.7100000e+05, 3.1874853e+01], [ 1.7200000e+05, 3.9018012e+01], [ 1.7300000e+05, 4.5510583e+01], [ 1.7400000e+05, 5.0796187e+01], [ 1.7500000e+05, 5.4331397e+01], [ 1.7600000e+05, 5.5641858e+01], [ 1.7700000e+05, 5.4381317e+01], [ 1.7800000e+05, 5.0385286e+01], [ 1.7900000e+05, 4.3710275e+01], [ 1.8000000e+05, 3.4650704e+01], [ 1.8100000e+05, 2.3728739e+01], [ 1.8200000e+05, 1.1656869e+01], [ 1.8300000e+05, -7.2202200e-01], [ 1.8400000e+05, -1.2508440e+01], [ 1.8500000e+05, -2.2823554e+01], [ 1.8600000e+05, -3.0880622e+01], [ 1.8700000e+05, -3.6044438e+01], [ 1.8800000e+05, -3.7876803e+01], [ 1.8900000e+05, -3.6169264e+01], [ 1.9000000e+05, -3.0965334e+01], [ 1.9100000e+05, -2.2572737e+01], [ 1.9200000e+05, -1.1562738e+01], [ 1.9300000e+05, 1.2496440e+00], [ 1.9400000e+05, 1.4851690e+01], [ 1.9500000e+05, 2.8109618e+01], [ 1.9600000e+05, 3.9873624e+01], [ 1.9700000e+05, 4.9100751e+01], [ 1.9800000e+05, 5.4977682e+01], [ 1.9900000e+05, 5.7022038e+01], [ 2.0000000e+05, 5.5142948e+01], [ 2.0100000e+05, 4.9649209e+01], [ 2.0200000e+05, 4.1204143e+01], [ 2.0300000e+05, 3.0737200e+01], [ 2.0400000e+05, 1.9330049e+01], [ 2.0500000e+05, 8.0975760e+00], [ 2.0600000e+05, -1.9184740e+00], [ 2.0700000e+05, -9.8319450e+00], [ 2.0800000e+05, -1.4967301e+01], [ 2.0900000e+05, -1.6890340e+01], [ 2.1000000e+05, -1.5421014e+01], [ 2.1100000e+05, -1.0635103e+01], [ 2.1200000e+05, -2.8584200e+00], [ 2.1300000e+05, 7.3466160e+00], [ 2.1400000e+05, 1.9205688e+01], [ 2.1500000e+05, 3.1773704e+01], [ 2.1600000e+05, 4.3996254e+01], [ 2.1700000e+05, 5.4791627e+01], [ 2.1800000e+05, 6.3148388e+01], [ 2.1900000e+05, 6.8228156e+01], [ 2.2000000e+05, 6.9459225e+01], [ 2.2100000e+05, 6.6606026e+01], [ 2.2200000e+05, 5.9801643e+01], [ 2.2300000e+05, 4.9536628e+01], [ 2.2400000e+05, 3.6605078e+01], [ 2.2500000e+05, 2.2016783e+01], [ 2.2600000e+05, 6.8901170e+00], [ 2.2700000e+05, -7.6574290e+00], [ 2.2800000e+05, -2.0605945e+01], [ 2.2900000e+05, -3.1107042e+01], [ 2.3000000e+05, -3.8531582e+01], [ 2.3100000e+05, -4.2493604e+01], [ 2.3200000e+05, -4.2856148e+01], [ 2.3300000e+05, -3.9724822e+01], [ 2.3400000e+05, -3.3432518e+01], [ 2.3500000e+05, -2.4515332e+01], [ 2.3600000e+05, -1.3677329e+01], [ 2.3700000e+05, -1.7416230e+00], [ 2.3800000e+05, 1.0412401e+01], [ 2.3900000e+05, 2.1921943e+01], [ 2.4000000e+05, 3.2015976e+01], [ 2.4100000e+05, 4.0082634e+01], [ 2.4200000e+05, 4.5718857e+01], [ 2.4300000e+05, 4.8755129e+01], [ 2.4400000e+05, 4.9252669e+01], [ 2.4500000e+05, 4.7475374e+01], [ 2.4600000e+05, 4.3842796e+01], [ 2.4700000e+05, 3.8872707e+01], [ 2.4800000e+05, 3.3122013e+01]])
année2=np.loadtxt("donnees_insolation.dat",usecols=0) année2
array([ 1000., 2000., 3000., 4000., 5000., 6000., 7000., 8000., 9000., 10000., 11000., 12000., 13000., 14000., 15000., 16000., 17000., 18000., 19000., 20000., 21000., 22000., 23000., 24000., 25000., 26000., 27000., 28000., 29000., 30000., 31000., 32000., 33000., 34000., 35000., 36000., 37000., 38000., 39000., 40000., 41000., 42000., 43000., 44000., 45000., 46000., 47000., 48000., 49000., 50000., 51000., 52000., 53000., 54000., 55000., 56000., 57000., 58000., 59000., 60000., 61000., 62000., 63000., 64000., 65000., 66000., 67000., 68000., 69000., 70000., 71000., 72000., 73000., 74000., 75000., 76000., 77000., 78000., 79000., 80000., 81000., 82000., 83000., 84000., 85000., 86000., 87000., 88000., 89000., 90000., 91000., 92000., 93000., 94000., 95000., 96000., 97000., 98000., 99000., 100000., 101000., 102000., 103000., 104000., 105000., 106000., 107000., 108000., 109000., 110000., 111000., 112000., 113000., 114000., 115000., 116000., 117000., 118000., 119000., 120000., 121000., 122000., 123000., 124000., 125000., 126000., 127000., 128000., 129000., 130000., 131000., 132000., 133000., 134000., 135000., 136000., 137000., 138000., 139000., 140000., 141000., 142000., 143000., 144000., 145000., 146000., 147000., 148000., 149000., 150000., 151000., 152000., 153000., 154000., 155000., 156000., 157000., 158000., 159000., 160000., 161000., 162000., 163000., 164000., 165000., 166000., 167000., 168000., 169000., 170000., 171000., 172000., 173000., 174000., 175000., 176000., 177000., 178000., 179000., 180000., 181000., 182000., 183000., 184000., 185000., 186000., 187000., 188000., 189000., 190000., 191000., 192000., 193000., 194000., 195000., 196000., 197000., 198000., 199000., 200000., 201000., 202000., 203000., 204000., 205000., 206000., 207000., 208000., 209000., 210000., 211000., 212000., 213000., 214000., 215000., 216000., 217000., 218000., 219000., 220000., 221000., 222000., 223000., 224000., 225000., 226000., 227000., 228000., 229000., 230000., 231000., 232000., 233000., 234000., 235000., 236000., 237000., 238000., 239000., 240000., 241000., 242000., 243000., 244000., 245000., 246000., 247000., 248000.])
insolation=np.loadtxt("donnees_insolation.dat",usecols=1) insolation
array([ 0.932811, 3.276957, 7.007533, 11.962172, 17.848958, 24.265751, 30.730781, 36.723055, 41.729583, 45.295127, 47.068648, 46.840655, 44.566351, 40.37112 , 34.537838, 27.477815, 19.689633, 11.71151 , 4.073083, -2.748176, -8.362152, -12.490773, -14.981564, -15.810417, -15.074452, -12.976025, -9.79916 , -5.88008 , -1.574236, 2.777025, 6.876522, 10.494468, 13.482176, 15.774981, 17.384905, 18.385185, 18.889915, 19.032421, 18.945647, 18.747018, 18.529075, 18.356124, 18.266137, 18.276541, 18.392081, 18.613076, 18.942238, 19.388852, 19.969419, 20.704481, 21.611976, 22.698076, 23.946962, 25.311289, 26.705217, 28.001646, 29.03529 , 29.612374, 29.527145, 28.584286, 26.625132, 23.554535, 19.364473, 14.150408, 8.117105, 1.572077, -5.093288, -11.432886, -16.981845, -21.295142, -23.984956, -24.753789, -23.421449, -19.944873, -14.430034, -7.134944, 1.537704, 11.059873, 20.810521, 30.120493, 38.327281, 44.834616, 49.169996, 51.032574, 50.324677, 47.16242 , 41.864095, 34.918377, 26.937266, 18.60031 , 10.597132, 3.573917, -1.912132, -5.428773, -6.695338, -5.601157, -2.216737, 3.202311, 10.220567, 18.240342, 26.535026, 34.300476, 40.723734, 45.064138, 46.737592, 45.391607, 40.958166, 33.673972, 24.063086, 12.884015, 1.050052, -10.463809, -20.713079, -28.863359, -34.250634, -36.421145, -35.153753, -30.469128, -22.628698, -12.123412, 0.350416, 13.931837, 27.654687, 40.524006, 51.603032, 60.100429, 65.445507, 67.340296, 65.780866, 61.045514, 53.653118, 44.299097, 33.779427, 22.913345, 12.473923, 3.132527, -4.5799 , -10.295168, -13.81592 , -15.112329, -14.309865, -11.669709, -7.562568, -2.436458, 3.220409, 8.919417, 14.210571, 18.716638, 22.158416, 24.368362, 25.29162 , 24.975669, 23.551518, 21.210426, 18.180215, 14.704516, 11.027033, 7.381478, 3.986579, 1.044626, -1.258416, -2.753191, -3.289261, -2.740485, -1.01384 , 1.938648, 6.105829, 11.402794, 17.657812, 24.603454, 31.874853, 39.018012, 45.510583, 50.796187, 54.331397, 55.641858, 54.381317, 50.385286, 43.710275, 34.650704, 23.728739, 11.656869, -0.722022, -12.50844 , -22.823554, -30.880622, -36.044438, -37.876803, -36.169264, -30.965334, -22.572737, -11.562738, 1.249644, 14.85169 , 28.109618, 39.873624, 49.100751, 54.977682, 57.022038, 55.142948, 49.649209, 41.204143, 30.7372 , 19.330049, 8.097576, -1.918474, -9.831945, -14.967301, -16.89034 , -15.421014, -10.635103, -2.85842 , 7.346616, 19.205688, 31.773704, 43.996254, 54.791627, 63.148388, 68.228156, 69.459225, 66.606026, 59.801643, 49.536628, 36.605078, 22.016783, 6.890117, -7.657429, -20.605945, -31.107042, -38.531582, -42.493604, -42.856148, -39.724822, -33.432518, -24.515332, -13.677329, -1.741623, 10.412401, 21.921943, 32.015976, 40.082634, 45.718857, 48.755129, 49.252669, 47.475374, 43.842796, 38.872707, 33.122013])
plt.plot(année,insolation) plt.title('graphique insolation en fonction des années') plt.ylabel('Insolation') plt.xlabel('Années en ans')
Text(0.5, 0, 'Années en ans')
Image in a Jupyter notebook

Question : En utilisant la fonction subplots() ou subplot(), faites les deux graphiques côte à côte sur la même figure avec les bons titres d'axes.

plt.subplot(221) plt.plot(année,O) plt.title('graphique $\delta^{18}O$ en fonction des années') plt.subplot(222) plt.plot(année2,insolation) plt.title('graphique insolation en fonction des années') plt.subplots_adjust(wspace=6, hspace=6)
Image in a Jupyter notebook