13426c051f58eb0be12ede8db913f48283a33afc
[Projet_Recherche_Operationnelle.git] / rapport / ProjetOptimRO.tex
1 \documentclass[12pt,oneside,a4paper]{book}
2
3
4 %%%%%Packages
5
6
7 \usepackage{latexsym}
8 \usepackage{amsmath}
9 \usepackage{mathtools}
10 \usepackage{amssymb}
11 \usepackage[utf8]{inputenc}
12 \usepackage[francais]{babel}
13 \usepackage{color}
14 \usepackage{geometry}
15 \usepackage{graphicx}
16 \usepackage{amsfonts}
17 \usepackage[T1]{fontenc}
18 \usepackage{multirow}
19 \usepackage{fancyhdr}
20 \usepackage{tocbibind}
21 \usepackage{lmodern}
22
23
24 %%%%%Marges & en-t\^etes
25
26 \geometry{hmargin=2.3cm, vmargin=3cm}
27 \fancyhf{} % supprime les en-t\^etes et pieds pr\'ed\'efinis
28 \fancyhead[FC]{\bfseries\thepage} % N∞page centre bas
29 \fancyhead[HC]{\footnotesize\leftmark} % chapitre centre haut
30 \renewcommand{\headrulewidth}{0.2pt} % filet en haut
31 \addtolength{\headheight}{0.5pt} % espace pour le filet
32 \renewcommand{\footrulewidth}{0.2pt} % filet en bas
33
34
35 %%%%%Th\'eor\`eme et d\'efinitions
36
37 \newtheorem{Def}{D\'efinition}
38 \newtheorem{Not}[Def]{Notation}
39 \newtheorem{Th}{Th\'eor\`eme}
40 \newtheorem{Prop}[Th]{Proposition}
41 \newtheorem{Cor}[Th]{Corollaire}
42 \newtheorem{Rmq}{Remarque}
43
44 \newcommand{\norme}[1]{\left\Vert #1\right\Vert}
45
46 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
47 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
48
49 \begin{document}
50
51 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
52 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
53
54 %%%%%Page de garde
55
56 \begin{center}
57
58 %\includegraphics[scale=0.5]{logo_sciences_rvb.png}\\
59 \includegraphics[scale=0.5]{polytech.png}\\
60
61 \vspace*{0.5cm}
62
63 \footnotesize{
64 \large \bf D\'epartement d'Informatique, Réseaux et Multimédia\\
65 \large \bf 5ème année\\
66 }
67
68 \vspace*{0.5cm}
69
70 %\large{Master 2 Professionnel\\
71 %Math\'ematiques et Informatique des Nouvelles Technologies\\}
72
73 \large{Projet \\ en \\ Optimisation et Recherche Opérationnelle \\}
74
75 \vspace*{0.7cm}
76
77 \begin{tabular}{c}
78 \hline
79 ~ \\
80 \LARGE\textbf {Programmation Séquentielle Quadratique} \\
81 \LARGE\textbf {en} \\
82 \LARGE\textbf {Optimisation non linéraire sous contraintes} \\
83 ~ \\
84 \hline
85 \end{tabular}
86
87 \vspace*{0.7cm}
88
89 \includegraphics[scale=0.4]{CE.PNG}\\
90
91 \vspace*{0.5cm}
92
93 \large par\\
94
95 %\large \bsc{}\\
96 %\normalsize{M\'emoire encadr\'e par :} \large St\'ephane \bsc{Ballet}\\
97
98 \vspace*{0.2cm}
99 \large {\bf Jérôme \bsc{Benoit} et Sylvain \bsc{Papa}}\\
100
101 %\vspace*{0.1cm}
102
103 % \large sous la direction de \\
104
105 %\vspace*{0.1cm}
106
107 %Eric Audureau et Thierry Masson
108
109 %\vspace*{1cm}
110
111 \vspace*{1cm}
112
113 %\normalsize{Licence de Mathématiques 3ème année}
114 \normalsize{Année 2018-2019}
115
116 \end{center}
117
118 \thispagestyle{empty}
119
120 \newpage
121
122
123 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
124 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
125
126
127 \pagestyle{plain}
128 \frontmatter
129
130
131 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
132 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
133
134
135 %%%%%Table des mati\`eres
136
137 \tableofcontents
138
139 \begin{figure}[!b]
140 \begin{center}
141 %\includegraphics{logo_fac2}
142 \includegraphics[scale=0.04]{amu}
143 \end{center}
144 \end{figure}
145
146 \newpage
147
148
149 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
150 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
151
152
153 \mainmatter
154 \pagestyle{fancy}
155
156
157 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
158 \chapter{Introduction générale}
159
160 \vspace{.5em}
161
162 \section{Qu'est-ce que la recherche opérationnelle?}
163
164 \subsection{Présentation rapide}
165
166 La recherche opérationnelle est une discipline dite "hybride" au confluent de plusieurs disciplines dont principalement l'analyse numérique, les probabilités, la statistique et l'algorithmie.
167 \newline
168 On la considère usuellement comme une sous discipline des mathématiques de la décision.
169
170 \subsection{Définition de la problèmatique}
171
172 Définissons le problème central $ \mathcal{P} $ que ce propose de résoudre la recherche opérationnelle.
173 \begin{Def}
174 Soient $(n, p, q) \in \mathbb{N}^3$, $x \in \mathbb{R}^n$, une fonction $g: \mathbb{R}^n \longrightarrow \mathbb{R}^p$ représentant les contraintes d'inégalités, une fonction $h: \mathbb{R}^n \longrightarrow \mathbb{R}^q$ représentant les contraintes d'égalités et une fonction dite objectif $J: \mathbb{R}^n \longrightarrow \mathbb{R}$.
175 \newline
176 La problèmatique $ \mathcal{P} $ se définit par :
177 $$
178 \mathcal{P} \left \{
179 \begin{array}{r c l}
180 \displaystyle\min_{x \in \mathbb{R}^n} J(x) \\
181 g(x) \leq 0 \\
182 h(x) = 0
183 \end{array}
184 \right .
185 $$
186 \end{Def}
187 \begin{Def}
188 On définit $ \mathcal{C} $ l'ensemble des contraintes par :
189 $$ \mathcal{C} = \left \{ x \in \mathbb{R}^n \ | \ g(x) \leq 0 \land h(x) = 0 \right \} $$
190 \end{Def}
191 Elle se doit de résoudre les problèmes d'existence d'une solution ($ \mathcal{C} \neq \emptyset $) ainsi que de construction d'une solution.
192
193 \section{Qu'est-ce que l'optimisation?}
194
195 La recherche d'un optimum au problème $ \mathcal{P} $ est l'activité principale de l'optimisation.
196 \begin{Def}
197 Soient une fonction $ f: \mathbb{R}^n \longrightarrow \mathbb{R} $
198 et $ x^\ast, h \in \mathbb{R}^n $.
199 On dit que $ f $ est différentiable en $ x^\ast $ si il existe une application linéraire $ d_{x^\ast}f $ de $ \mathbb{R}^n $ dans $ \mathbb{R} $ telle que
200 \[
201 f(x^\ast + h) = f(x^\ast) + d_{x^\ast}f(h) + \underset{h \rightarrow 0}{\mathrm{o}}(\norme{h})
202 \]
203 Autrement dit il existe une application $ \varepsilon_{x^\ast} $ définie sur le voisinage de $ 0 $ dans $ \mathbb{R}^n $ et à valeurs dans $ \mathbb{R} $
204 telle que $ \lim\limits_{h \rightarrow 0} \varepsilon_{x^\ast}(h) = 0 $ et
205 \[
206 f(x^\ast + h) = f(x^\ast) + d_{x^\ast}f(h) + \norme{h}\varepsilon_{x^\ast}(h)
207 \]
208 \end{Def}
209 On peut confondre la somme des dérivées partielles et la fonction linéaire $ d_{x^\ast}f $.
210 \begin{Def}
211 Soit une fonction $ f: \mathbb{R}^n \longrightarrow \mathbb{R} $ différentiable.
212 \newline
213 Le gradient de $ f $, noté $\nabla f$, en $ x^\ast \in \mathbb{R}^n$ se définit par :
214 \[
215 \nabla f(x^\ast) = (\frac{\partial f}{\partial x_1}(x^\ast),\ldots,\frac{\partial f}{\partial x_n}(x^\ast))
216 \]
217 \end{Def}
218 Dans le cas où $ J $ est continûment différentiable et ses dérivées sont continues (ou de classe $ \mathcal{C}^1 $),
219 une condition suffisante et nécessaire pour que $ x^\ast \in \mathbb{R}^n $ soit un de ses extremums local ou global est que $ \nabla f(x^\ast) = 0 $.
220 \newline
221 Dans ce projet, nous nous proposons d'étudier une des méthodes d'optimisation non linéaire avec contraintes nommée programmation quadratique séquentielle.
222
223 % Dans cette section nous prenons appui sur l'ouvrage {\it Optimisation et contrôle des systèmes linéaires} \cite{Berg} de Maïtine Bergounioux \footnote{Maïtine Bergounioux, {\it Optimisation et contrôle des systèmes linéaires}, Dunod, 2001.}.
224 % Nous utiliserons aussi l'ouvrage de Francis Filbet\footnote{Francis Filbet, {\it Analyse numérique - Algorithme et étude mathématique}, Dunod, 2009.}, {\it Analyse numérique - Algorithme et étude mathématique} \cite{Filb}.
225
226 %{\it La relativité}, Que sais-je?, 4ème édition, puf, 2000, \cite{Mavr};
227 %ainsi que Jean Hladik, {\it La relativité selon Einstein}, L'esprit des sciences, Ellipses, 2000, \cite{Hlad}.
228
229
230 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
231
232 \chapter{Méthodes de programmation quadratique séquentielle}
233
234 % \section{Cahier des charges}
235 %
236 % Il s'agit de travailler en binôme ou bien seul sur des sujets complémentaires et d'approfondissement du cours. Le travail en question effectué durant les TDs consistera
237 % à effectuer un dossier sur un thème. Le dossier devra être tapé en Latex ou Tex puisque il peut y avoir des formules de mathématiques ou de physiques. Il pourra aussi comporter une partie "implémentation effective" d'algorithmes (en annexe).
238 %
239 % \vspace{.5em}
240 %
241 % Sur la fond, toutes les sources de connaissance utilisées devront être citées. En particulier, la méthodologie universitaire sera privilégiée
242 % (citations en note de bas de page et dans le corps du document, liste des références en fin de document dans la bibliographie, etc...).
243 % Wikipédia pourra être utilisé mais cela devra être mentionné en tant que référence (note de bas de page ou citation dans le corps du document).
244 % L'accent sera essentiellement mis sur la démarche scientifique utilisée à égal niveau avec le contenu acquis des connaissances.
245 %
246 % \vspace{.5em}
247 %
248 % Plusieurs sources devront être croisées afin de prétendre au maximum de vraisemblance
249 % et d'objectivité scientifique. Le document ne devra pas excéder 10 pages.
250 % On privilégiera les qualités de synthèse, d'organisation ainsi que du contenu du document.
251 %
252 % \section{Proposition de sujets}
253 %
254 % \subsection{Analyse numérique}
255 %
256 % \vspace{.5em}
257 %
258 % 1) Méthode des moindres Carrés (cas général, cas pondéré, cas des équations non linéaires).
259 %
260 % \vspace{.5em}
261 %
262 % 2) Méthode de Newton-Raphson (cas d'une variable, cas de deux variables) - Application: extrema d'une fonction à deux variables.
263 %
264 % \vspace{.5em}
265 %
266 % 3) Autres méthodes: méthode de Jacobi, de Gauss-Seidel, etc....
267 %
268 % \vspace{.5em}
269
270 \section{Optimisation}
271
272 % \vspace{.5em}
273
274 % \subsubsection{Optimisation sans contrainte}
275 %
276 % {\bf A- Algorithmes déterministes}
277 %
278 % \vspace{.5em}
279 %
280 % 1) Régression linéaire sans contrainte (pré-requis: Méthode des moindres carrés).
281 %
282 % \vspace{.5em}
283 %
284 % 2) Méthodes de descente: la méthode du gradient (à pas constant ou à pas variable ou à pas optimal).
285 %
286 % \vspace{.5em}
287 %
288 % 3) Méthode de Newton (ou méthode dite de la tangente) et application à la recherche d'extrema.
289 %
290 % \vspace{.5em}
291 %
292 % 4) Méthodes de descente: méthode du gradient conjugué (cas linéaire et cas général)
293 %
294 % \vspace{.5em}
295 %
296 % 5) Méthode de relaxation
297 %
298 % \vspace{.5em}
299 %
300 % {\bf B- Algorithmes probabilistes ou dit stochastiques}
301 %
302 % \vspace{.5em}
303 %
304 % 1) Dynamique de métropolis (prérequis: chaines de Markov)
305 %
306 % \vspace{.5em}
307 %
308 % 2) Recuit simulé sur un ensemble fini et application au problème du voyageur de commerce (prérequis: dynamique de métropolis)
309 %
310 % \vspace{.5em}
311
312 \subsubsection{Optimisation ou minimisation avec contraintes}
313
314 % \vspace{.5em}
315 %
316 % 1) Régression linéaire avec contraintes (prérequis: méthode des moindres carrés, conditions ou équations dites de Karush-kuhn-Tucker (KKT)) .
317 %
318 % \vspace{.5em}
319 %
320 % 2) Cas de la programmation linéaire (prérequis: Lagrangien et multiplicateurs de Lagrange, conditions de KKT).
321 %
322 % \vspace{.5em}
323 %
324 % 3) Algorithmes: méthode du gradient projeté, méthode de Lagrange-Newton pour des contraintes en égalité,
325 % méthode de Newton projetée pour des contraintes de bornes, méthodes de pénalisation,
326 % méthodes de programmation quadratique successive (SQP Sequential Quadratic Programming),
327 % méthode de dualité (méthode d'Uzawa, prérequis: théorie de la dualité convexe) etc...
328 %
329 % \vspace{.5em}
330 %
331 % \subsection{Recherche opérationnelle}
332 %
333 % \vspace{.5em}
334 %
335 % \subsubsection{La programmation linéaire (cas particulier de l'optimisation avec contraintes)}
336 %
337 % 1) Méthode d'énumération.
338 %
339 % \vspace{.5em}
340 %
341 % 2) Méthode du simplexe.
342 %
343 % \vspace{.5em}
344 %
345 % 3) Application à des problèmes de R.O:
346 %
347 % \vspace{.5em}
348 %
349 % \hspace{.3em} 3.1) Fêtes de Pâques: A l'approche des fêtes de Pâques, un artisan chocolatier décide de confectionner des oeufs en chocolats. En allant inspecter ses réserves, il constate qu'il lui reste 18 kg de cacao, 8 kg de noisettes et 14 litres de lait. Ce chocolatier a deux spécialités: l'oeuf {\it extra} et l'oeuf {\it sublime}. Un oeuf {\it extra} nécessite 1kg de cacao, 1 kg de noisettes et 2 litres de lait tandis qu'un oeuf {\it sublime} nécessite 3 kg de cacao, 1 kg de noisettes et 1 litre de lait. Il fera un bénéfice de 20 euros en vendant un oeuf {\it extra}, et de 30 euros en vendant un oeuf {\it sublime}.
350 %
351 % \vspace{.5em}
352 %
353 % \hspace{.6em} a) \'Ecrire ce problème sous la forme d'un problème de programmation linéaire.
354 %
355 % \vspace{.5em}
356 %
357 % \hspace{.6em} b) Combien d'oeufs extra et sublime doit-il fabriquer pour faire le plus grand bénéfice?
358 %
359 % \vspace{.5em}
360 %
361 % \hspace{.3em} 3.2) Organisation du travail: La fabrication d'une pièce $P_1$ a un prix de revient de 150 euros et celle d'une pièce $P_2$ coûte 100 euros. Chaque pièce est traitée successivement dans trois ateliers. Le nombre d'heures-machines par pièce est indiqué dans le tableau suivant :
362 %
363 % \vspace{.5em}
364 %
365 % \begin{center}
366 % $
367 % \begin{array}{|c|c|c|c|}
368 % \hline
369 % Atelier & A & B & C \\
370 % \hline
371 % Pièce 1 & 3 h & 5 h & 2 h \\
372 % \hline
373 % Pièce 2 & 1 h & 3 h & 3 h \\
374 % \hline
375 % \end{array}
376 % $
377 % \end{center}
378 %
379 % \vspace{.5em}
380 %
381 % Pour éviter le chômage technique, l'atelier A doit obligatoirement fournir 1200 heures machines, l'atelier B doit obligatoirement fournir 3000 heures machines et l'atelier C doit obligatoirement fournir 1800 heures machines.
382 %
383 % \hspace{.6em} a) \'Ecrire ce problème sous la forme d'un problème de programmation linéaire.
384 %
385 % \vspace{.5em}
386 %
387 % \hspace{.6em} b) Combien faut-il fabriquer de pièces $P_1$ et $P_2$ pour minimiser le coût de revient de l'ensemble de la production et pour assurer le fonctionnement des trois ateliers excluant tout chômage technique?
388 %
389 % \vspace{.5em}
390
391 \bibliographystyle{plain}
392 \bibliography{stdlib_sbphilo}
393
394 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
395
396 \end{document}
397
398
399 \begin{thebibliography}{6}\input{MemoireM2Ballet6.synctex.gz(busy)}
400
401 %\bibitem[1]{BL} Jean-Pierre \bsc{Bourguignon} et David \bsc{Langlois}, Cours de M1, Module Relativité Générale,
402 %Ecole Polytechnique, ParisTech, 2011.\\
403
404 %\bibitem[2]{G} Gilles \bsc{Cohen-Tannoudji}, Einstein et la refondation relativiste de la physique, 2005.\\
405
406 %\bibitem[3]{D} Pierre \bsc{Duhem}, La théorie physique, son objet, sa structure, Vrin, 2007.\\
407
408 %\bibitem[4]{E1} Albert \bsc{Einstein}, Die formale grundlage der allgemeinen Relativittstheorie. Kniglich Preussische
409 %Akademie der Wissenschaften (Berlin),Sitzungsberichte: pp 1030-1085. \\
410
411 %\bibitem[5]{G} Christian \bsc{Godin}, Dictionnaire de philosophie, Fayard Edition du temps, 2004.\\
412
413 %\bibitem[6]{H} Jean \bsc{Hladik}, La Relativité selon Einstein, L'Esprit des Sciences, Ellipses.\\
414
415 %\bibitem[7]{IS} \bsc{Iftime} and \bsc{Stachel}, The hole argument for covariant theories, arKiv:gr-qc/0512021v2, 8 avril 2006.\\
416
417 %\bibitem[8]{K} \bsc{Kant}, Critique de la raison pure, Traduction, présentation, notes par Alain Renaut, GF-Flammarion, 2006.\\
418
419 %\bibitem[9]{K2} \bsc{Kant}, Prolégomènes à toute métaphysique future, Traduction de Louis Guilermit, Vrin, 1986.\\
420
421 %\bibitem[10]{KU} Thomas \bsc{Kuhn}, La structure des révolutions scientifiques, Flammarion Champs Sciences, 2008.
422
423 %\bibitem[11]{L} Marc \bsc{Lachièze-Rey}, Initiation à la cosmologie, 3ème édition, Dunod, 2000.\\
424
425 %\bibitem[12]{Mas} Thierry \bsc{Masson}, Cours de géométrie différentielle, groupe et algèbre de Lie, fibrés et connexions, 2010.\\
426
427 %\bibitem[13]{Poi} Henri \bsc{Poincaré}, La Science et L'Hypothèse, Flammarion, Paris, 1968.\\
428
429 %\bibitem[14]{Mav} Stamatia \bsc{Mavridès}, La Relativité, Que sais-je, 4ème édition, PUF, 2000.\\
430
431 %\bibitem[15]{R} Robert \bsc{Rynasiewicz}, The Lessons of the Hole Argument, The British Journal of the Philosophy of Science,
432 %vol; 45 (2), 407-436, Oxford University Press, Oxford Journals, 1994. \\
433
434 %\bibitem[16]{S} Standford Encyclopedia of Philosophy.\\
435
436 %\bibitem[17]{W} Wikipedia.\\
437
438 %\bibitem[1]{Bachtold} {\bf Manuel Bächtold}, L'interprétation de la mécanique quantique, une approche pragmatique, Collection vision des sciences, Hermann, 2008 .\\
439
440 %\bibitem[2]{Aspect} {\bf Alain Aspect}, Présentation naïve des inégalités de Bell, 2004.\\
441
442 % \bibitem[3]{Basda} {\bf Jean-Louis Basdevant et Manuel Joffre}, Mécanique Quantique, Les éditions de l'Ecole Polytechnique, 2006.\\
443
444 %\bibitem[4]{Diu} {\bf Bernard Diu}, Le congrès de Solvay de 1927: petite chronique d'un grand évènement, Bibnum.\\
445
446 %\bibitem[1]{B} \bsc{Aristote}, Métaphysique, traduction J.Tricot, Vrin, 1974.\\
447
448 \end{thebibliography}