Skip to content

M
mooc-rr
Commit dd588c42 authored by cb1e5ba91280d02583d0604166e71c1c's avatar cb1e5ba91280d02583d0604166e71c1c
Browse files
Options

no commit message

parent d35b638c
Show whitespace changes
Inline Side-by-side
Showing with 1 addition and 928 deletions
module3/exo3/exercice.ipynb
View file @ dd588c42
...@@ -1311,934 +1311,7 @@ ...@@ -1311,934 +1311,7 @@
"- Recherche de ces deux fréquences par FFT. \n", "- Recherche de ces deux fréquences par FFT. \n",
"- Nous affinons l'estimation de ces fréquences à l'aide d'une approximation sinusoïdale (fréquence, amplitude, moyenne et phase) par les moindres carrés ordinaires.\n", "- Nous affinons l'estimation de ces fréquences à l'aide d'une approximation sinusoïdale (fréquence, amplitude, moyenne et phase) par les moindres carrés ordinaires.\n",
"\n", "\n",
"Nous considérons pour l'échantillonnage des données un pas constant, correspond à un intervalle de 1 mois $Te = 1/12$. Cet échantillonnage avec le nombre de point de mesure nous donne un échantillonnage en fréquence, ce pas est $1/(N*Te) = 0.016 Hz$ \n" "Nous considérons pour l'échantillonnage des données un pas constant, correspond à un intervalle de 1 mois $Te = \\frac {1}{12}$. Cet échantillonnage avec le nombre de point de mesure nous donne un échantillonnage en fréquence, ce pas est $\\frac {1} {(N*Te)} = 0.016 Hz$ \n"
]
},
{
"cell_type": "code",
"execution_count": 150,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>index</th>\n",
" <th>Yr</th>\n",
" <th>Mn</th>\n",
" <th>Date 1</th>\n",
" <th>Date 2</th>\n",
" <th>s1</th>\n",
" <th>s2</th>\n",
" <th>s3</th>\n",
" <th>s4</th>\n",
" <th>s5</th>\n",
" <th>s6</th>\n",
" </tr>\n",
" <tr>\n",
" <th>period</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1958-06-01</th>\n",
" <td>5</td>\n",
" <td>1958</td>\n",
" <td>6</td>\n",
" <td>21351</td>\n",
" <td>1958.4548</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-10-01</th>\n",
" <td>9</td>\n",
" <td>1958</td>\n",
" <td>10</td>\n",
" <td>21473</td>\n",
" <td>1958.7890</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>312.44</td>\n",
" <td>315.40</td>\n",
" <td>312.44</td>\n",
" <td>315.40</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-02-01</th>\n",
" <td>73</td>\n",
" <td>1964</td>\n",
" <td>2</td>\n",
" <td>23422</td>\n",
" <td>1964.1257</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>320.01</td>\n",
" <td>319.36</td>\n",
" <td>320.01</td>\n",
" <td>319.36</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-03-01</th>\n",
" <td>74</td>\n",
" <td>1964</td>\n",
" <td>3</td>\n",
" <td>23451</td>\n",
" <td>1964.2049</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>320.74</td>\n",
" <td>319.41</td>\n",
" <td>320.74</td>\n",
" <td>319.41</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-04-01</th>\n",
" <td>75</td>\n",
" <td>1964</td>\n",
" <td>4</td>\n",
" <td>23482</td>\n",
" <td>1964.2896</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>321.83</td>\n",
" <td>319.45</td>\n",
" <td>321.83</td>\n",
" <td>319.45</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" index Yr Mn Date 1 Date 2 s1 s2 s3 s4 \\\n",
"period \n",
"1958-06-01 5 1958 6 21351 1958.4548 NaN NaN 317.24 315.14 \n",
"1958-10-01 9 1958 10 21473 1958.7890 NaN NaN 312.44 315.40 \n",
"1964-02-01 73 1964 2 23422 1964.1257 NaN NaN 320.01 319.36 \n",
"1964-03-01 74 1964 3 23451 1964.2049 NaN NaN 320.74 319.41 \n",
"1964-04-01 75 1964 4 23482 1964.2896 NaN NaN 321.83 319.45 \n",
"\n",
" s5 s6 \n",
"period \n",
"1958-06-01 317.24 315.14 \n",
"1958-10-01 312.44 315.40 \n",
"1964-02-01 320.01 319.36 \n",
"1964-03-01 320.74 319.41 \n",
"1964-04-01 321.83 319.45 "
]
},
"execution_count": 150,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df = data.reset_index().copy()\n",
"df['period'] = pd.Series([datetime.date(y,m,1) for y,m in zip(df['Yr'],df['Mn'])])\n",
"df = df.set_index('period') \n",
"#df.head()\n",
"d = df[df.isnull().any(axis=1)]\n",
"d"
]
},
{
"cell_type": "code",
"execution_count": 151,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>index</th>\n",
" <th>Yr</th>\n",
" <th>Mn</th>\n",
" <th>Date 1</th>\n",
" <th>Date 2</th>\n",
" <th>s1</th>\n",
" <th>s2</th>\n",
" <th>s3</th>\n",
" <th>s4</th>\n",
" <th>s5</th>\n",
" <th>s6</th>\n",
" </tr>\n",
" <tr>\n",
" <th>period</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>2019-09-01</th>\n",
" <td>740</td>\n",
" <td>2019</td>\n",
" <td>9</td>\n",
" <td>43723</td>\n",
" <td>2019.7068</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" <td>408.52</td>\n",
" <td>412.05</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-10-01</th>\n",
" <td>741</td>\n",
" <td>2019</td>\n",
" <td>10</td>\n",
" <td>43753</td>\n",
" <td>2019.7890</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" <td>408.67</td>\n",
" <td>412.29</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-11-01</th>\n",
" <td>742</td>\n",
" <td>2019</td>\n",
" <td>11</td>\n",
" <td>43784</td>\n",
" <td>2019.8740</td>\n",
" <td>410.28</td>\n",
" <td>412.56</td>\n",
" <td>410.29</td>\n",
" <td>412.54</td>\n",
" <td>410.28</td>\n",
" <td>412.56</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-12-01</th>\n",
" <td>743</td>\n",
" <td>2019</td>\n",
" <td>12</td>\n",
" <td>43814</td>\n",
" <td>2019.9562</td>\n",
" <td>411.85</td>\n",
" <td>412.78</td>\n",
" <td>411.88</td>\n",
" <td>412.79</td>\n",
" <td>411.85</td>\n",
" <td>412.78</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2020-01-01</th>\n",
" <td>744</td>\n",
" <td>2020</td>\n",
" <td>1</td>\n",
" <td>43845</td>\n",
" <td>2020.0410</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" <td>413.11</td>\n",
" <td>413.05</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" index Yr Mn Date 1 Date 2 s1 s2 s3 \\\n",
"period \n",
"2019-09-01 740 2019 9 43723 2019.7068 408.55 412.06 408.52 \n",
"2019-10-01 741 2019 10 43753 2019.7890 408.43 412.07 408.67 \n",
"2019-11-01 742 2019 11 43784 2019.8740 410.28 412.56 410.29 \n",
"2019-12-01 743 2019 12 43814 2019.9562 411.85 412.78 411.88 \n",
"2020-01-01 744 2020 1 43845 2020.0410 413.37 413.33 413.11 \n",
"\n",
" s4 s5 s6 \n",
"period \n",
"2019-09-01 412.05 408.55 412.06 \n",
"2019-10-01 412.29 408.43 412.07 \n",
"2019-11-01 412.54 410.28 412.56 \n",
"2019-12-01 412.79 411.85 412.78 \n",
"2020-01-01 413.05 413.37 413.33 "
]
},
"execution_count": 151,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df.tail()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Pour les valeur manquantes, nous allons faire une interpolation linéaire afin que la FFT dispose de données prisent à échantillonnage de 1 mois."
]
},
{
"cell_type": "code",
"execution_count": 152,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>index</th>\n",
" <th>Yr</th>\n",
" <th>Mn</th>\n",
" <th>Date 1</th>\n",
" <th>Date 2</th>\n",
" <th>s1</th>\n",
" <th>s2</th>\n",
" <th>s3</th>\n",
" <th>s4</th>\n",
" <th>s5</th>\n",
" <th>s6</th>\n",
" </tr>\n",
" <tr>\n",
" <th>period</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1958-03-01</th>\n",
" <td>2</td>\n",
" <td>1958</td>\n",
" <td>3</td>\n",
" <td>21259</td>\n",
" <td>1958.2027</td>\n",
" <td>315.700</td>\n",
" <td>314.44</td>\n",
" <td>316.19</td>\n",
" <td>314.91</td>\n",
" <td>315.70</td>\n",
" <td>314.44</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-04-01</th>\n",
" <td>3</td>\n",
" <td>1958</td>\n",
" <td>4</td>\n",
" <td>21290</td>\n",
" <td>1958.2877</td>\n",
" <td>317.450</td>\n",
" <td>315.16</td>\n",
" <td>317.30</td>\n",
" <td>314.99</td>\n",
" <td>317.45</td>\n",
" <td>315.16</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-05-01</th>\n",
" <td>4</td>\n",
" <td>1958</td>\n",
" <td>5</td>\n",
" <td>21320</td>\n",
" <td>1958.3699</td>\n",
" <td>317.510</td>\n",
" <td>314.71</td>\n",
" <td>317.86</td>\n",
" <td>315.06</td>\n",
" <td>317.51</td>\n",
" <td>314.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-06-01</th>\n",
" <td>5</td>\n",
" <td>1958</td>\n",
" <td>6</td>\n",
" <td>21351</td>\n",
" <td>1958.4548</td>\n",
" <td>316.685</td>\n",
" <td>314.95</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-07-01</th>\n",
" <td>6</td>\n",
" <td>1958</td>\n",
" <td>7</td>\n",
" <td>21381</td>\n",
" <td>1958.5370</td>\n",
" <td>315.860</td>\n",
" <td>315.19</td>\n",
" <td>315.86</td>\n",
" <td>315.22</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-08-01</th>\n",
" <td>7</td>\n",
" <td>1958</td>\n",
" <td>8</td>\n",
" <td>21412</td>\n",
" <td>1958.6219</td>\n",
" <td>314.930</td>\n",
" <td>316.19</td>\n",
" <td>314.00</td>\n",
" <td>315.29</td>\n",
" <td>314.93</td>\n",
" <td>316.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-09-01</th>\n",
" <td>8</td>\n",
" <td>1958</td>\n",
" <td>9</td>\n",
" <td>21443</td>\n",
" <td>1958.7068</td>\n",
" <td>313.210</td>\n",
" <td>316.08</td>\n",
" <td>312.46</td>\n",
" <td>315.35</td>\n",
" <td>313.21</td>\n",
" <td>316.08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-07-01</th>\n",
" <td>738</td>\n",
" <td>2019</td>\n",
" <td>7</td>\n",
" <td>43661</td>\n",
" <td>2019.5370</td>\n",
" <td>411.850</td>\n",
" <td>411.03</td>\n",
" <td>412.36</td>\n",
" <td>411.57</td>\n",
" <td>411.85</td>\n",
" <td>411.03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-08-01</th>\n",
" <td>739</td>\n",
" <td>2019</td>\n",
" <td>8</td>\n",
" <td>43692</td>\n",
" <td>2019.6219</td>\n",
" <td>410.080</td>\n",
" <td>411.62</td>\n",
" <td>410.23</td>\n",
" <td>411.81</td>\n",
" <td>410.08</td>\n",
" <td>411.62</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-09-01</th>\n",
" <td>740</td>\n",
" <td>2019</td>\n",
" <td>9</td>\n",
" <td>43723</td>\n",
" <td>2019.7068</td>\n",
" <td>408.550</td>\n",
" <td>412.06</td>\n",
" <td>408.52</td>\n",
" <td>412.05</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-10-01</th>\n",
" <td>741</td>\n",
" <td>2019</td>\n",
" <td>10</td>\n",
" <td>43753</td>\n",
" <td>2019.7890</td>\n",
" <td>408.430</td>\n",
" <td>412.07</td>\n",
" <td>408.67</td>\n",
" <td>412.29</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-11-01</th>\n",
" <td>742</td>\n",
" <td>2019</td>\n",
" <td>11</td>\n",
" <td>43784</td>\n",
" <td>2019.8740</td>\n",
" <td>410.280</td>\n",
" <td>412.56</td>\n",
" <td>410.29</td>\n",
" <td>412.54</td>\n",
" <td>410.28</td>\n",
" <td>412.56</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-12-01</th>\n",
" <td>743</td>\n",
" <td>2019</td>\n",
" <td>12</td>\n",
" <td>43814</td>\n",
" <td>2019.9562</td>\n",
" <td>411.850</td>\n",
" <td>412.78</td>\n",
" <td>411.88</td>\n",
" <td>412.79</td>\n",
" <td>411.85</td>\n",
" <td>412.78</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2020-01-01</th>\n",
" <td>744</td>\n",
" <td>2020</td>\n",
" <td>1</td>\n",
" <td>43845</td>\n",
" <td>2020.0410</td>\n",
" <td>413.370</td>\n",
" <td>413.33</td>\n",
" <td>413.11</td>\n",
" <td>413.05</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>743 rows × 11 columns</p>\n",
"</div>"
],
"text/plain": [
" index Yr Mn Date 1 Date 2 s1 s2 s3 \\\n",
"period \n",
"1958-03-01 2 1958 3 21259 1958.2027 315.700 314.44 316.19 \n",
"1958-04-01 3 1958 4 21290 1958.2877 317.450 315.16 317.30 \n",
"1958-05-01 4 1958 5 21320 1958.3699 317.510 314.71 317.86 \n",
"1958-06-01 5 1958 6 21351 1958.4548 316.685 314.95 317.24 \n",
"1958-07-01 6 1958 7 21381 1958.5370 315.860 315.19 315.86 \n",
"1958-08-01 7 1958 8 21412 1958.6219 314.930 316.19 314.00 \n",
"1958-09-01 8 1958 9 21443 1958.7068 313.210 316.08 312.46 \n",
"... ... ... .. ... ... ... ... ... \n",
"2019-07-01 738 2019 7 43661 2019.5370 411.850 411.03 412.36 \n",
"2019-08-01 739 2019 8 43692 2019.6219 410.080 411.62 410.23 \n",
"2019-09-01 740 2019 9 43723 2019.7068 408.550 412.06 408.52 \n",
"2019-10-01 741 2019 10 43753 2019.7890 408.430 412.07 408.67 \n",
"2019-11-01 742 2019 11 43784 2019.8740 410.280 412.56 410.29 \n",
"2019-12-01 743 2019 12 43814 2019.9562 411.850 412.78 411.88 \n",
"2020-01-01 744 2020 1 43845 2020.0410 413.370 413.33 413.11 \n",
"\n",
" s4 s5 s6 \n",
"period \n",
"1958-03-01 314.91 315.70 314.44 \n",
"1958-04-01 314.99 317.45 315.16 \n",
"1958-05-01 315.06 317.51 314.71 \n",
"1958-06-01 315.14 317.24 315.14 \n",
"1958-07-01 315.22 315.86 315.19 \n",
"1958-08-01 315.29 314.93 316.19 \n",
"1958-09-01 315.35 313.21 316.08 \n",
"... ... ... ... \n",
"2019-07-01 411.57 411.85 411.03 \n",
"2019-08-01 411.81 410.08 411.62 \n",
"2019-09-01 412.05 408.55 412.06 \n",
"2019-10-01 412.29 408.43 412.07 \n",
"2019-11-01 412.54 410.28 412.56 \n",
"2019-12-01 412.79 411.85 412.78 \n",
"2020-01-01 413.05 413.37 413.33 \n",
"\n",
"[743 rows x 11 columns]"
]
},
"execution_count": 152,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df = df.interpolate(method='linear', limit_direction='forward',limit=3)\n",
"df"
]
},
{
"cell_type": "code",
"execution_count": 153,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>index</th>\n",
" <th>Yr</th>\n",
" <th>Mn</th>\n",
" <th>Date 1</th>\n",
" <th>Date 2</th>\n",
" <th>s1</th>\n",
" <th>s2</th>\n",
" <th>s3</th>\n",
" <th>s4</th>\n",
" <th>s5</th>\n",
" <th>s6</th>\n",
" </tr>\n",
" <tr>\n",
" <th>period</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1964-01-01</th>\n",
" <td>72</td>\n",
" <td>1964</td>\n",
" <td>1</td>\n",
" <td>23391</td>\n",
" <td>1964.0410</td>\n",
" <td>319.5700</td>\n",
" <td>319.5300</td>\n",
" <td>319.37</td>\n",
" <td>319.32</td>\n",
" <td>319.57</td>\n",
" <td>319.53</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-02-01</th>\n",
" <td>73</td>\n",
" <td>1964</td>\n",
" <td>2</td>\n",
" <td>23422</td>\n",
" <td>1964.1257</td>\n",
" <td>320.2425</td>\n",
" <td>319.4975</td>\n",
" <td>320.01</td>\n",
" <td>319.36</td>\n",
" <td>320.01</td>\n",
" <td>319.36</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-03-01</th>\n",
" <td>74</td>\n",
" <td>1964</td>\n",
" <td>3</td>\n",
" <td>23451</td>\n",
" <td>1964.2049</td>\n",
" <td>320.9150</td>\n",
" <td>319.4650</td>\n",
" <td>320.74</td>\n",
" <td>319.41</td>\n",
" <td>320.74</td>\n",
" <td>319.41</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-04-01</th>\n",
" <td>75</td>\n",
" <td>1964</td>\n",
" <td>4</td>\n",
" <td>23482</td>\n",
" <td>1964.2896</td>\n",
" <td>321.5875</td>\n",
" <td>319.4325</td>\n",
" <td>321.83</td>\n",
" <td>319.45</td>\n",
" <td>321.83</td>\n",
" <td>319.45</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-05-01</th>\n",
" <td>76</td>\n",
" <td>1964</td>\n",
" <td>5</td>\n",
" <td>23512</td>\n",
" <td>1964.3716</td>\n",
" <td>322.2600</td>\n",
" <td>319.4000</td>\n",
" <td>322.34</td>\n",
" <td>319.49</td>\n",
" <td>322.26</td>\n",
" <td>319.40</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-06-01</th>\n",
" <td>77</td>\n",
" <td>1964</td>\n",
" <td>6</td>\n",
" <td>23543</td>\n",
" <td>1964.4563</td>\n",
" <td>321.8900</td>\n",
" <td>319.7500</td>\n",
" <td>321.64</td>\n",
" <td>319.52</td>\n",
" <td>321.89</td>\n",
" <td>319.75</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-07-01</th>\n",
" <td>78</td>\n",
" <td>1964</td>\n",
" <td>7</td>\n",
" <td>23573</td>\n",
" <td>1964.5383</td>\n",
" <td>320.4400</td>\n",
" <td>319.7800</td>\n",
" <td>320.19</td>\n",
" <td>319.55</td>\n",
" <td>320.44</td>\n",
" <td>319.78</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-08-01</th>\n",
" <td>79</td>\n",
" <td>1964</td>\n",
" <td>8</td>\n",
" <td>23604</td>\n",
" <td>1964.6230</td>\n",
" <td>318.6900</td>\n",
" <td>320.0000</td>\n",
" <td>318.24</td>\n",
" <td>319.58</td>\n",
" <td>318.69</td>\n",
" <td>320.00</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-09-01</th>\n",
" <td>80</td>\n",
" <td>1964</td>\n",
" <td>9</td>\n",
" <td>23635</td>\n",
" <td>1964.7077</td>\n",
" <td>316.7100</td>\n",
" <td>319.6500</td>\n",
" <td>316.64</td>\n",
" <td>319.60</td>\n",
" <td>316.71</td>\n",
" <td>319.65</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-10-01</th>\n",
" <td>81</td>\n",
" <td>1964</td>\n",
" <td>10</td>\n",
" <td>23665</td>\n",
" <td>1964.7896</td>\n",
" <td>316.8700</td>\n",
" <td>319.9100</td>\n",
" <td>316.59</td>\n",
" <td>319.62</td>\n",
" <td>316.87</td>\n",
" <td>319.91</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-11-01</th>\n",
" <td>82</td>\n",
" <td>1964</td>\n",
" <td>11</td>\n",
" <td>23696</td>\n",
" <td>1964.8743</td>\n",
" <td>317.6800</td>\n",
" <td>319.5800</td>\n",
" <td>317.76</td>\n",
" <td>319.64</td>\n",
" <td>317.68</td>\n",
" <td>319.58</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-12-01</th>\n",
" <td>83</td>\n",
" <td>1964</td>\n",
" <td>12</td>\n",
" <td>23726</td>\n",
" <td>1964.9563</td>\n",
" <td>318.7100</td>\n",
" <td>319.4900</td>\n",
" <td>318.89</td>\n",
" <td>319.66</td>\n",
" <td>318.71</td>\n",
" <td>319.49</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" index Yr Mn Date 1 Date 2 s1 s2 s3 \\\n",
"period \n",
"1964-01-01 72 1964 1 23391 1964.0410 319.5700 319.5300 319.37 \n",
"1964-02-01 73 1964 2 23422 1964.1257 320.2425 319.4975 320.01 \n",
"1964-03-01 74 1964 3 23451 1964.2049 320.9150 319.4650 320.74 \n",
"1964-04-01 75 1964 4 23482 1964.2896 321.5875 319.4325 321.83 \n",
"1964-05-01 76 1964 5 23512 1964.3716 322.2600 319.4000 322.34 \n",
"1964-06-01 77 1964 6 23543 1964.4563 321.8900 319.7500 321.64 \n",
"1964-07-01 78 1964 7 23573 1964.5383 320.4400 319.7800 320.19 \n",
"1964-08-01 79 1964 8 23604 1964.6230 318.6900 320.0000 318.24 \n",
"1964-09-01 80 1964 9 23635 1964.7077 316.7100 319.6500 316.64 \n",
"1964-10-01 81 1964 10 23665 1964.7896 316.8700 319.9100 316.59 \n",
"1964-11-01 82 1964 11 23696 1964.8743 317.6800 319.5800 317.76 \n",
"1964-12-01 83 1964 12 23726 1964.9563 318.7100 319.4900 318.89 \n",
"\n",
" s4 s5 s6 \n",
"period \n",
"1964-01-01 319.32 319.57 319.53 \n",
"1964-02-01 319.36 320.01 319.36 \n",
"1964-03-01 319.41 320.74 319.41 \n",
"1964-04-01 319.45 321.83 319.45 \n",
"1964-05-01 319.49 322.26 319.40 \n",
"1964-06-01 319.52 321.89 319.75 \n",
"1964-07-01 319.55 320.44 319.78 \n",
"1964-08-01 319.58 318.69 320.00 \n",
"1964-09-01 319.60 316.71 319.65 \n",
"1964-10-01 319.62 316.87 319.91 \n",
"1964-11-01 319.64 317.68 319.58 \n",
"1964-12-01 319.66 318.71 319.49 "
]
},
"execution_count": 153,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d = df[(df.Yr >= 1964) & (df.Yr < 1965)]\n",
"d"
]
},
{
"cell_type": "code",
"execution_count": 154,
"metadata": {},
"outputs": [],
"source": [
"data_cube = df.copy()\n",
"popt, pcov = curve_fit(func_cube,data_cube['Date 2'],data_cube['s1'])\n",
"data_cube['reg_cube'] = fcube(data_cube['Date 2'])\n",
"data_cube['co2'] = data_cube['reg_cube']-data_cube['s1']\n",
"\n",
"nb = len(data_cube['Date 2'])\n",
"ecart = []\n",
"for i in np.arange(0,nb-1):\n",
" ecart.append(data_cube['Date 2'][i+1] - data_cube['Date 2'][i])\n",
"dt = np.mean(ecart)"
] ]
}, },
{ {
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment