Skip to content

rapport labo 5 - #5

Open
bouayoubegrine wants to merge 2 commits into
dpo-mth8408:mainfrom
bouayoubegrine:TP5
Open

bouayoubegrine wants to merge 2 commits into
dpo-mth8408:mainfrom
bouayoubegrine:TP5

Conversation

@bouayoubegrine

Copy link
Copy Markdown

No description provided.

Comment thread rapport.qmd Outdated
@github-actions

Copy link
Copy Markdown

Status:
Success--Here-is-the-PDF

Comment thread rapport.qmd
\begin{align*}
\min_{x, y} \ & \int_0^1 f(x(t), \dot{x}(t), t) \, \mathrm{d}t \\
\text{s.t.} \ & \dot{y}(t) = h(x(t), \dot{x}(t), t) \\
x(0) = x_0, \ & x(1) = x_1, \ y(1) = y(0),

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

La dernière condition est incorrecte. On doit avoir

$$y(0) = \int_0^0 h = 0$$

et

$$y(1) = \int_0^1 h = 0.$$

De façon générale,

$$y(t) = \int_0^t h.$$

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

$y(t) = \int_0^1 h$ est une primitive de $h$, qui s'annule en $0$, mais il n'est pas nécessaire que toutes le primitives s'annulent en $0$ (se different avec une constante qui change rien dans ce problème). Je pense pas que $y(1) = y(0)$ est incorrect dans le sens qu'il change le problème (peut etre incomplète). J'ai expliqué après qu'on peut fixer la constante à $0$.

Comment thread rapport.qmd
& \dot{x}(t) = u(t) \, \dot{y}(t) = \sqrt{1 + u(t)^2} - L\\
& x(0) = a, \quad x(1) = b, \\
& y(1) = y(0) = 0.
\end{align*}

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ok.

Comment thread rapport.qmd

#### Commenter les résultats

- Premièrement, la corde prend une forme parabolique, ce qui paraît intuitif. De même, la fonction $p_1(t)$ est une fonction non constante, comme vu en cours.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

En fait, ce n'est pas une parabole.

Comment thread rapport.qmd

- Premièrement, la corde prend une forme parabolique, ce qui paraît intuitif. De même, la fonction $p_1(t)$ est une fonction non constante, comme vu en cours.

- En ce qui concerne le temps d'exécution, le solveur utilisé a nécessité beaucoup d'itérations pour converger. On observe une convergence superlinéaire lors des dernières itérations.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

👍

Comment thread rapport.qmd

- Le problème suivant est de nature \emph{bang-bang}. Ce que je remarque, c’est que le réservoir 1 dépasse, en un temps inférieur à 1, la moitié du réservoir. À partir de ce moment, le contrôle devient nul afin que le réservoir soit à moitié plein à $t_f = 1$.

- Dans ce deuxième problème, j’ai imposé que le réservoir ne soit pas rempli à plus de la moitié avant $t_f$. On obtient alors un contrôle non lisse, avec de fortes variations, qui tente de se compenser pour que le réservoir atteigne la moitié de sa capacité à l’instant final.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Intéressant.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants