Correction
U n i v e r s i t é M o u l a y I s m aı̈ l Contrôle Continu RNE
E.N.S.A.M.–M e k n è s 3ème Année ENSAM, Semestre 1
Jeudi 07/11/2019 - Durée : 40 mn
Nom : ..................................................... Prénom : ................................. Code ENSAM : ..........................
C ET M P T TR U V • Cochez votre code d’étudiant ENSAM ci-contre caractère par car-
0 1 2 3 4 5 6 7 8 9 actère.
0 1 2 3 4 5 6 7 8 9 • Aucun document n’est autorisé.
0 1 2 3 4 5 6 7 8 9 • Pas de calculatrice.
0 1 2 3 4 5 6 7 8 9
• Les réponses fausses ou incohérentes retirent des points.
0 1 2 3 4 5 6 7 8 9
Partie I : Cocher une seule réponse Question 5 Définir une fonction Matlab rungekutta2 qui résout le
Question 1 Avec un pas du temps h = 12 , quelle approximation de la valeur au problème de Cauchy par le schéma de Runge-Kutta d’ordre 2 (RK2).
temps 1 de l’éventuelle solution x(t) du problème de Cauchy : Cases réservées au correcteur
x(0) = 4, et x0 (t) = −x(t)2 + 8t ∀ t ∈ [0, 1], fournit le schéma de Runge Kutta
d’ordre 2 ? [ RK2 est donné par : xn+1 = xn + h f tn + h2 , yn + h2 f (tn , yn ) ].
function [T,Y]=rungekutta2(F,y0,t0,tf,N)
65/32 113/128 Aucune
255/32 11/64 213/64
Question 2 Avec un pas du temps h = 21 , quelle approximation de la valeur au
temps 1 de l’éventuelle solution x(t) du problème de Cauchy :
x(0) = 2, et x0 (t) = −x(t)3 + t2 x(t) ∀ t ∈ [0, 1], fournit le schéma d’Euler explicite ?
27/128 7/4 1/16
7/32 Aucune 1/4
y 000 (t) + y(t)y 00 (t) − y 0 (t)2 = 0, t ∈ [0, tf ]
Partie II : Soit l’EDO suivante :
y(0) = y 0 (0) = 0, y 0 (0) = 1.
Question 3 Écrire ce problème sous forme d’un système différentiel du 1er ordre
: Y 0 (t) = F (t, Y (t)). Cases réservées au correcteur
Question 6 Compléter le programme Matlab ci-dessous pour réaliser le calcul de
la solution approchée sur l’intervalle [0, tf ] avec tf = 10 en utilisant la méthode de
Runge-Kutta d’ordre 2, et proposer une représentation graphique sur cet intervalle.
Cases réservées au correcteur
%FICHIER : "prg_principal.m"
%Programme principal
clc;
clear all;
close all;
tf=10; % temps final
Question 4 Écrire une fonction Matlab fCauchy qui calcule les valeurs de N=100; % Nombre de pas de temps
F (t, Y ). Cases réservées au correcteur
function Yp=fCauchy(t,Y)