From 15aa6e43e54dc42e9f5a547a35b8c607b74cf20f Mon Sep 17 00:00:00 2001
From: taviera <alois.tavier@student.uclouvain.be>
Date: Sun, 21 May 2023 11:55:03 +0200
Subject: [PATCH] 1.6

---
 Notebook_final.ipynb | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)

diff --git a/Notebook_final.ipynb b/Notebook_final.ipynb
index f0c41aa..ce1e648 100644
--- a/Notebook_final.ipynb
+++ b/Notebook_final.ipynb
@@ -689,6 +689,30 @@
     "probleme.solve(solver=cp.GLPK_MI, warm_start=True, verbose=False)"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "4189d852",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "Duree = 8\n",
+    "x_i = cp.Variable(n//Duree, boolean=True) #Allumage\n",
+    "\n",
+    "# Formulation et résolution du problème\n",
+    "objectif = cp.Minimize(c.T@(p_n_i+p_r_i))\n",
+    "contraintes = [T_i[0] == T_initial, T_i[n] == T_final, T_min <= T_i, T_i <= T_max, T_i[1:n+1] == T_i[0:n]+deltaT_i, p_n_i>=0, p_r_i>=0, \n",
+    "               deltaT_i == -(1-eta)*(T_i[0:n]-T_ext)+(cp.multiply(p_n_i,COP_normal(T_ext)))*0.25/(V*Cx)-p_r_i*COP_reverse*0.25/(V*Cx)]\n",
+    "\n",
+    "for i in range(n):\n",
+    "    contraintes.append((p_n_i[i]+p_r_i[i])<=p_max*x_i[i//Duree])\n",
+    "    contraintes.append((p_n_i[i]+p_r_i[i])>=0.25*p_max*x_i[i//Duree])\n",
+    "#contraintes.append((p_n_i+p_r_i)<=p_max)\n",
+    "\n",
+    "probleme = cp.Problem(objectif, contraintes)\n",
+    "probleme.solve(solver=cp.GLPK_MI, warm_start=True, verbose=False)"
+   ]
+  },
   {
    "attachments": {},
    "cell_type": "markdown",
-- 
GitLab