arthur/ArtStudies
1
0
Fork 0
mirror of https://github.com/ArthurDanjou/ArtStudies.git synced 2026-01-14 18:59:59 +01:00
Code Issues Packages Projects Releases Wiki Activity

add: TP1/2

Browse Source
This commit is contained in:
Arthur Danjou
2025-03-24 16:07:17 +01:00
parent 4853ad1d64
commit 52c6012197
2 changed files with 47 additions and 56 deletions
Download Patch File Download Diff File Expand all files Collapse all files

75
M1/Numerical Methods/TP1.ipynb
View File

@@ -170,8 +170,8 @@
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:26.900678Z",
"start_time": "2025-03-24T13:04:26.896193Z"
"end_time": "2025-03-24T13:14:56.636016Z",
"start_time": "2025-03-24T13:14:56.632263Z"
}
},
"cell_type": "code",
@@ -185,85 +185,85 @@
],
"id": "61588e94b56ea779",
"outputs": [],
"execution_count": 22
"execution_count": 41
},
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:27.128889Z",
"start_time": "2025-03-24T13:04:27.123995Z"
"end_time": "2025-03-24T13:14:57.036853Z",
"start_time": "2025-03-24T13:14:57.033785Z"
}
},
"cell_type": "code",
"source": [
"def euler_explicit_sir(N, T, u0, f, args):\n",
"def euler_explicit_sir(N, T, y0, f, args):\n",
" t = np.linspace(0, T, N + 1)\n",
" y = np.zeros((N + 1, len(y0)))\n",
" y[0] = y0\n",
" h = T / N\n",
" for i in range(N):\n",
" y[i + 1] = y[i] + h * F(y[i], t[i], *args)\n",
" y[i + 1] = y[i] + h * f(y[i], t[i], *args)\n",
" return t, y"
],
"id": "a32f0facb958a1a0",
"outputs": [],
"execution_count": 23
"execution_count": 42
},
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:27.315627Z",
"start_time": "2025-03-24T13:04:27.312627Z"
"end_time": "2025-03-24T13:14:57.487407Z",
"start_time": "2025-03-24T13:14:57.484369Z"
}
},
"cell_type": "code",
"source": [
"def heun_sir(N, T, u0, f, args):\n",
"def heun_sir(N, T, y0, f, args):\n",
" t = np.linspace(0, T, N + 1)\n",
" y = np.zeros((N + 1, len(y0)))\n",
" y[0] = y0\n",
" h = T / N\n",
" for i in range(N):\n",
" k1 = h * F(y[i], t[i], *args)\n",
" k2 = h * F(y[i] + k1, t[i + 1], *args)\n",
" k1 = h * f(y[i], t[i], *args)\n",
" k2 = h * f(y[i] + k1, t[i + 1], *args)\n",
" y[i + 1] = y[i] + 0.5 * (k1 + k2)\n",
" return t, y"
],
"id": "ca061c36b09f7200",
"outputs": [],
"execution_count": 24
"execution_count": 43
},
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:27.489907Z",
"start_time": "2025-03-24T13:04:27.485432Z"
"end_time": "2025-03-24T13:14:57.869526Z",
"start_time": "2025-03-24T13:14:57.865654Z"
}
},
"cell_type": "code",
"source": [
"def runge_kutta_order_4_sir(N, T, u0, f, args):\n",
"def runge_kutta_order_4_sir(N, T, y0, f, args):\n",
" t = np.linspace(0, T, N + 1)\n",
" y = np.zeros((N + 1, len(y0)))\n",
" y[0] = y0\n",
" h = T / N\n",
" for i in range(N):\n",
" k1 = h * F(y[i], t[i], *args)\n",
" k2 = h * F(y[i] + 0.5 * k1, t[i] + 0.5 * h, *args)\n",
" k3 = h * F(y[i] + 0.5 * k2, t[i] + 0.5 * h, *args)\n",
" k4 = h * F(y[i] + k3, t[i + 1], *args)\n",
" k1 = h * f(y[i], t[i], *args)\n",
" k2 = h * f(y[i] + 0.5 * k1, t[i] + 0.5 * h, *args)\n",
" k3 = h * f(y[i] + 0.5 * k2, t[i] + 0.5 * h, *args)\n",
" k4 = h * f(y[i] + k3, t[i + 1], *args)\n",
" y[i + 1] = y[i] + (k1 + 2 * k2 + 2 * k3 + k4) / 6\n",
" return t, y"
],
"id": "afc8571b14a6b464",
"outputs": [],
"execution_count": 25
"execution_count": 44
},
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:27.725395Z",
"start_time": "2025-03-24T13:04:27.640769Z"
"end_time": "2025-03-24T13:14:58.376791Z",
"start_time": "2025-03-24T13:14:58.220423Z"
}
},
"cell_type": "code",
@@ -298,13 +298,13 @@
"output_type": "display_data"
}
],
"execution_count": 26
"execution_count": 45
},
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:28.172102Z",
"start_time": "2025-03-24T13:04:27.807260Z"
"end_time": "2025-03-24T13:22:36.716443Z",
"start_time": "2025-03-24T13:22:36.310338Z"
}
},
"cell_type": "code",
@@ -356,7 +356,7 @@
"output_type": "display_data"
}
],
"execution_count": 27
"execution_count": 48
},
{
"metadata": {},
@@ -389,8 +389,8 @@
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:04:30.150912Z",
"start_time": "2025-03-24T13:04:29.984979Z"
"end_time": "2025-03-24T13:46:24.490267Z",
"start_time": "2025-03-24T13:46:24.315085Z"
}
},
"cell_type": "code",
@@ -444,20 +444,7 @@
"output_type": "display_data"
}
],
"execution_count": 29
},
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T12:48:05.902735Z",
"start_time": "2025-03-24T12:48:05.890888Z"
}
},
"cell_type": "code",
"source": "",
"id": "aeb67c5158f41a0b",
"outputs": [],
"execution_count": null
"execution_count": 53
}
],
"metadata": {

28
M1/Numerical Methods/TP2.ipynb
View File

@@ -17,13 +17,12 @@
{
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T13:10:01.794956Z",
"start_time": "2025-03-24T13:09:58.777694Z"
"end_time": "2025-03-24T14:45:40.568860Z",
"start_time": "2025-03-24T14:45:37.895550Z"
}
},
"cell_type": "code",
"source": [
"\n",
"import numpy as np\n",
"from sklearn.neural_network import MLPClassifier\n",
"from sklearn.datasets import make_classification\n",
@@ -42,7 +41,7 @@
" y_pred = model.predict(X_test)\n",
" accuracies.append(accuracy_score(y_test, y_pred))\n",
"\n",
"print(f'Accuracy: {np.mean(accuracies) * 100:.0f}%')\n",
"print(f'Mean Accuracy: {np.mean(accuracies) * 100:.0f}%')\n",
"print(f\"Max accuracy: {np.max(accuracies) * 100:.0f}%\")\n",
"print(f\"Min accuracy: {np.min(accuracies) * 100:.0f}%\")"
],
@@ -52,21 +51,26 @@
"name": "stdout",
"output_type": "stream",
"text": [
"Accuracy: 93%\n",
"Max accuracy: 100%\n",
"Min accuracy: 80%\n"
"Mean Accuracy: 92%\n",
"Max accuracy: 96%\n",
"Min accuracy: 86%\n"
]
}
],
"execution_count": 21
"execution_count": 28
},
{
"metadata": {},
"metadata": {
"ExecuteTime": {
"end_time": "2025-03-24T14:37:53.507776Z",
"start_time": "2025-03-24T14:37:53.505376Z"
}
},
"cell_type": "code",
"outputs": [],
"execution_count": null,
"source": "",
"id": "96b6d46883ed5570"
"id": "96b6d46883ed5570",
"outputs": [],
"execution_count": null
}
],
"metadata": {