diff --git a/module3/exo3/exercice.ipynb b/module3/exo3/exercice.ipynb index c5de64b39c45910f30cd85ea2bc9003c4dea2e2c..c0b480e343662a4a5ac24b960c0586a7bb9a61a7 100644 --- a/module3/exo3/exercice.ipynb +++ b/module3/exo3/exercice.ipynb @@ -12,7 +12,7 @@ }, { "cell_type": "code", - "execution_count": 81, + "execution_count": 1, "metadata": { "hideCode": true, "hidePrompt": true @@ -44,7 +44,7 @@ }, { "cell_type": "code", - "execution_count": 82, + "execution_count": 2, "metadata": { "hideCode": true, "hidePrompt": true @@ -329,7 +329,7 @@ "[756 rows x 10 columns]" ] }, - "execution_count": 82, + "execution_count": 2, "metadata": {}, "output_type": "execute_result" } @@ -364,7 +364,7 @@ }, { "cell_type": "code", - "execution_count": 83, + "execution_count": 3, "metadata": { "hideCode": true, "hidePrompt": true @@ -651,7 +651,7 @@ "[756 rows x 10 columns]" ] }, - "execution_count": 83, + "execution_count": 3, "metadata": {}, "output_type": "execute_result" } @@ -675,7 +675,7 @@ }, { "cell_type": "code", - "execution_count": 84, + "execution_count": 4, "metadata": { "hideCode": true, "hidePrompt": true @@ -725,7 +725,7 @@ }, { "cell_type": "code", - "execution_count": 71, + "execution_count": 5, "metadata": { "hideCode": true, "hidePrompt": true @@ -747,7 +747,7 @@ "dtype: bool" ] }, - "execution_count": 71, + "execution_count": 5, "metadata": {}, "output_type": "execute_result" } @@ -769,7 +769,7 @@ }, { "cell_type": "code", - "execution_count": 85, + "execution_count": 6, "metadata": { "hideCode": true, "hidePrompt": true, @@ -917,7 +917,7 @@ "1958-07-01 315.22 315.86 315.19 " ] }, - "execution_count": 85, + "execution_count": 6, "metadata": {}, "output_type": "execute_result" } @@ -944,7 +944,7 @@ }, { "cell_type": "code", - "execution_count": 86, + "execution_count": 7, "metadata": { "hideCode": true, "hidePrompt": true @@ -981,7 +981,7 @@ }, { "cell_type": "code", - "execution_count": 88, + "execution_count": 8, "metadata": { "hideCode": true, "hidePrompt": true @@ -1039,7 +1039,7 @@ }, { "cell_type": "code", - "execution_count": 90, + "execution_count": 9, "metadata": { "hideCode": true, "hidePrompt": true, @@ -1114,7 +1114,7 @@ }, { "cell_type": "code", - "execution_count": 92, + "execution_count": 11, "metadata": { "hideCode": true, "hidePrompt": true @@ -1144,7 +1144,7 @@ "s5 = data_cube.co2 - data_cube.co2.mean()\n", "s5_fft = np.fft.fft(s5) # calcul de la transformee de Fourier et des frequences\n", "\n", - "n = s1.size\n", + "n = s5.size\n", "dt = 1/12\n", "pds_s5 = (dt/n) * np.abs(s5_fft)**2 # densité spectral de puissance\n", "freq = np.fft.fftfreq(n, d=dt) # fréquences associées\n", @@ -1171,7 +1171,7 @@ }, { "cell_type": "code", - "execution_count": 93, + "execution_count": 12, "metadata": { "hideCode": true, "hidePrompt": true @@ -1187,9 +1187,9 @@ } ], "source": [ - "indide_f = np.where(pds_s1 > 5)\n", + "indide_f = np.where(pds_s5 > 5)\n", "f_lente = indide_f[0][0] * f1\n", - "indide_f = np.where(pds_s1 > 20)\n", + "indide_f = np.where(pds_s5 > 20)\n", "f_rapide = indide_f[0][-1] * f1\n", " \n", "print(f' Estimation de la frequence lente : {f_lente} \\n Esitmation de la fréquence rapide : {f_rapide}')" @@ -1197,7 +1197,7 @@ }, { "cell_type": "code", - "execution_count": 94, + "execution_count": 13, "metadata": { "hideCode": true, "hidePrompt": true, @@ -1225,7 +1225,7 @@ "\n", "p0 = [f0,ampl0,phase0,offset0]\n", "x = np.array(data_cube['Date 2'])\n", - "popt_l, pcov_l = curve_fit(my_sin, x , s1, p0 = p0)\n", + "popt_l, pcov_l = curve_fit(my_sin, x , s5, p0 = p0)\n", "print(popt_l)\n", "ylent = my_sin(x,popt_l[0],popt_l[1],popt_l[2],popt_l[3])\n", "fl = p0[0]\n", @@ -1234,7 +1234,7 @@ "f0 = f_rapide \n", "ampl0= 3.5\n", "p0 = [f0,ampl0,phase0,offset0]\n", - "popt_r, pcov_r = curve_fit(my_sin, x , s1, p0 = p0)\n", + "popt_r, pcov_r = curve_fit(my_sin, x , s5, p0 = p0)\n", "print(popt_r)\n", "yrapide = my_sin(x,popt_r[0],popt_r[1],popt_r[2],popt_r[3]) \n", "fr = p0[0]" @@ -1242,7 +1242,7 @@ }, { "cell_type": "code", - "execution_count": 95, + "execution_count": 14, "metadata": { "hideCode": true, "hidePrompt": true, @@ -1270,14 +1270,14 @@ "ax2.grid(linestyle='--', linewidth=1)\n", "\n", "ax1.set(title = 'Variation concentration CO2 lente',xlabel='Période',ylabel='Concentration (ppm)')\n", - "ax1.plot(x,ylent,x,s1);\n", + "ax1.plot(x,ylent,x,s5);\n", "ax2.set(title = 'Variation concentration CO2 rapide',xlabel='Période',ylabel='Concentration (ppm)')\n", - "ax2.plot(x,yrapide,x,s1);" + "ax2.plot(x,yrapide,x,s5);" ] }, { "cell_type": "code", - "execution_count": 96, + "execution_count": 15, "metadata": { "hideCode": true, "hidePrompt": true @@ -1309,7 +1309,7 @@ }, { "cell_type": "code", - "execution_count": 102, + "execution_count": 16, "metadata": { "hideCode": true, "hidePrompt": true @@ -1636,7 +1636,7 @@ "[63 rows x 12 columns]" ] }, - "execution_count": 102, + "execution_count": 16, "metadata": {}, "output_type": "execute_result" } @@ -1648,7 +1648,7 @@ }, { "cell_type": "code", - "execution_count": 103, + "execution_count": 17, "metadata": { "hideCode": true, "hidePrompt": true @@ -1836,7 +1836,7 @@ "max 411.850000 410.080000 408.550000 408.430000 410.280000 411.850000 " ] }, - "execution_count": 103, + "execution_count": 17, "metadata": {}, "output_type": "execute_result" } @@ -1847,7 +1847,7 @@ }, { "cell_type": "code", - "execution_count": 163, + "execution_count": 18, "metadata": { "hideCode": true, "hidePrompt": true, @@ -2002,7 +2002,7 @@ "1962 0.7890 0.8740 0.9562 " ] }, - "execution_count": 163, + "execution_count": 18, "metadata": {}, "output_type": "execute_result" } @@ -2041,121 +2041,30 @@ "hidePrompt": true }, "outputs": [], - "source": [] + "source": [ + "d" + ] }, { "cell_type": "code", - "execution_count": 164, + "execution_count": null, "metadata": { "hideCode": true, "hidePrompt": true }, - "outputs": [ - { - "data": { - "text/plain": [ - "0.1246" - ] - }, - "execution_count": 164, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ "0.0416+0.083" ] }, { "cell_type": "code", - "execution_count": 165, + "execution_count": null, "metadata": { "hideCode": true, "hidePrompt": true }, - "outputs": [ - { - "data": { - "text/html": [ - "
\n", - "\n", - "\n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - "
indexYrMnDate 1Date 2s1s2s3s4s5s6
period
2020-01-0174420201438452020.041413.37413.33413.11413.05413.37413.33
\n", - "
" - ], - "text/plain": [ - " index Yr Mn Date 1 Date 2 s1 s2 s3 s4 \\\n", - "period \n", - "2020-01-01 744 2020 1 43845 2020.041 413.37 413.33 413.11 413.05 \n", - "\n", - " s5 s6 \n", - "period \n", - "2020-01-01 413.37 413.33 " - ] - }, - "execution_count": 165, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ "df[-1:]" ]