\NeedsTeXFormat{LaTeX2e}
\documentclass[10pt,a4paper]{amsart}
\usepackage{amsmath,amsfonts,amssymb,latexsym}
\usepackage{a4,mathptmx,helvet,courier}
\usepackage[francais]{babel}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}

\newcounter{exonum}[section]
\newenvironment{exo}{\begin{enumerate}\stepcounter{exonum}%
\renewcommand{\theenumi}{\thesection.\arabic{exonum}}%
\renewcommand{\labelenumi}{\bf\theenumi.}\item}{\end{enumerate}}

\newcommand{\N}{\mathbb{N}}
\newcommand{\Z}{\mathbb{Z}}
\newcommand{\Q}{\mathbb{Q}}
\newcommand{\R}{\mathbb{R}}
\newcommand{\C}{\mathbb{C}}
\usepackage{graphicx}

\begin{document} 

\pagestyle{empty}
\thispagestyle{empty}

{\footnotesize\noindent
Universit\'e Joseph Fourier, Grenoble I   \hfill $\bullet$ \hfill
L3 M\'ethodes num\'eriques                       \hfill $\bullet$ \hfill
Ann\'ee 2013/2014           \\[-2mm]\hrule}


\begin{center}
%  {\Large\bf Math\'ematique assist\'ees par ordinateur} \\
  Examen du 18 juin 2014, de 9h \`a 12h. \\
  \it Documents, calculatrices et ordinateurs ultraportables
  d\'econnect\'es du r\'eseau autoris\'es.\\
  Ce sujet comporte 2 pages. 
Bar\^eme donn\'e \`a titre indicatif et non contractuel.
\end{center}


\section{M\'ethode de la puissance (5 points)}
Soit la matrice
$M=\left(\begin{array}{ccc}
7 & 1 & -121 \\
1 & 0 & 5 \\
-121 & 5 & -3
\end{array}\right) $.
\begin{enumerate}
\item
Prendre un vecteur al\'eatoire $v$ \`a coordonn\'ees dans l'intervalle
$[0,1]$, lui appliquer 29 et 30 fois la matrice, 
en d\'eduire une estimation $\lambda$
 de la valeur propre $l$ de $M$ qui est la plus grande en module.
%Donner une majoration de $|l|$.
\item Soit $A$ une matrice de taille $n$ qui admet
une base orthonormale de vecteurs propres.
On suppose qu'on a d\'etermin\'e (par la m\'ethode de la puissance)
un vecteur norm\'e $v$ et une constante $\lambda$ tels que
$\| Av-\lambda v\| < \varepsilon$. Soit $(v_1,...,v_n)$
les coordonn\'ees de $v$ dans la base propre orthonormale
de $A$ associ\'ee aux valeurs propres $(\lambda_1,...,\lambda_n)$.
Calculer $\| Av-\lambda v\|$, en d\'eduire que l'une des valeurs
propres au moins v\'erifie $| \lambda_k-\lambda | < \varepsilon$.
\item 
Peut-on appliquer le r\'esultat du (2) \`a $M$~?
Donner un encadrement de $l$ (justifier).
\item Proposer une méthode permettant de calculer numériquement les
autres valeurs propres de $M$.
\end{enumerate}

\section{Système linéaire (2 points)}
Soit $M$ la matrice de l'exercice 1 et $b=(1,2,3)$ un vecteur de $\R^3$
connu avec une précision relative de \verb|1e-4| pour la norme euclidienne.
Donner la solution $x$ du système linéaire $Ax=b$ en arrondissant
\`a la pr\'ecision ad\'equate (on justifiera).

\section{Méthode itérative: 6 points}
On d\'efinit sur l'intervalle $[2,3]$ la fonction $f$ par
  $f(x)=\cos(x)+2x-\exp(x)$ et on cherche \`a r\'esoudre 
$f(x)=\lambda, \lambda \in \R$
\begin{enumerate}
\item D\'eterminer les valeurs de $\lambda$ telles que 
l'\'equation $f(x)=\lambda$ admette une solution unique sur $[2,3]$.
\item On r\'e\'ecrit l'\'equation sous la forme
  $x=\frac12(\lambda+\exp(x)-\cos(x))$,
peut-on appliquer la m\'ethode du point fixe pour la r\'esoudre~?
Si oui, donner (et justifier!) un encadrement de la solution \`a
\verb|1e-4| pr\`es par cette m\'ethode pour $\lambda=-10$.
\item Donner une suite it\'erative $(u_n)$ convergeant vers la solution 
en utilisant la m\'ethode de Newton, justifier la convergence pour
la valeur de $u_0$ choisie. D\'eterminer la solution de l'\'equation
$f(x)=-10$ avec une pr\'ecision de \verb|1e-8|.
\item \'Ecrire une fonction Xcas prenant en argument $\lambda$
et renvoyant une valeur approch\'ee de la solution \`a \verb|1e-8|
pr\`es de l'\'equation $f(x)=\lambda$
par une m\'ethode it\'erative de votre choix.
\end{enumerate}

\section{M\'ethode des trap\`ezes et polyn\^omes trigonom\'etriques (6 points)}
\'Etant donn\'e une  fonction continue 
$f : [0,2\pi]\rightarrow \mathbb{R}$,  on note $I(f) = \int_0^{2\pi} f(x) {\rm d} x$
son int\'egrale sur $[0,2\pi]$ et ${\rm Trap}_N (f)$ la valeur approch\'ee de $I(f)$ par la m\'ethode des 
trap\`ezes  avec un pas constant $h=2\pi/N$, o\`u  $N>0$ est un entier. 
\begin{enumerate}
\item
Soit $f$ un polyn\^ome trigonom\'etrique 
de degr\'e inf\'erieur ou \'egal \`a $N-1$, de la forme
%
$$
f(x) = \sum_{k=-N+1}^{N-1} c_k e^{{\rm i} k x}
$$
%
avec $c_k \in \mathbb{C}$, $c_{-k} = \overline{c}_k$, $k=-N+1,\ldots, N-1$.
Montrer que ${\rm Trap}_N (f)$ donne le r\'esultat exact pour l'int\'egrale $I(f)$.
\item On consid\`ere la fonction $g(x) =\exp ( \sin x)$.
En utilisant la question 1 et le d\'eveloppement de l'exponentielle en 0, montrer que
l'erreur $E(g) = |I (g) - {\rm Trap}_N (g)|$ est major\'ee par $c/N!$, o\`u $c$ est une constante que vous d\'eterminerez.
\item
D\'eterminez \`a l'aide de Xcas l'augmentation de la pr\'ecision avec $h$ quand on divise $h$ par deux
(vous pourrez par exemple utiliser votre programme pour la m\'ethode des trap\`ezes r\'ealis\'e en TP).
Comparez cette pr\'ecision avec celles obtenues par les m\'ethodes des rectangles \`a gauche et \`a droite. Pouvez-vous expliquer les similitudes/diff\'erences 
des r\'esultats obtenus par ces trois m\'ethodes ?   
\end{enumerate}

\section{Ordres des m\'ethodes de Runge-Kutta (5 points)}
On consid\`ere une m\'ethode   de Runge-Kutta pour r\'esoudre l'\'equation diff\'erentielle 
$y'(t) = f( y(t), t  )$ sur l'intervalle $[0,1]$, donn\'ee par le tableau de coefficients
%
$$
\begin{array}{lllllll}
c_0=0 &  \vline &       0     &                &          &                 & \\
c_1   & \vline & \lambda_{10} &       0        &          &                 & \\
c_2   & \vline & \lambda_{20} &   \lambda_{21}  &    0      &                 & \\
\vdots & \vline & \vdots     &    \vdots      &   \ddots   & \ddots                 & \\
c_q    &  \vline & \lambda_{q0} & \lambda_{q 1} & \cdots   & \lambda_{q\,q-1} &  0 \\
\hline
1      &  \vline &   \mu_0     &  \mu_1         & \cdots   & \mu_{q-1}       & \mu_q 
\end{array}   
$$
%
On a donc $z_{n+1} = z_n + h \Phi (z_n, t_n,h)$ o\`u $h$ est le pas (suppos\'e constant), $z_0 = y(0)$,  $t_n = n h$ et
%
$$
\Phi (z, t,h ) = \sum_{i=0}^q \mu_i f (  z_{n,i}, t + h c_i )
$$
avec
$$z_{n,0} = z \; \text{  et  }\; 
z_{n,i} = z + h \sum_{j=0}^{i-1} \lambda_{i j} f (  z_{n,j}, t + h c_j  )\;,\;i=1,\ldots, q \;.
$$
%   
On note $p$ l'ordre de la m\'ethode. On rappelle que $p>0$.
%
\begin{enumerate}
\item Supposons que $f(y,t )=f(t)$ soit ind\'ependante de $y$. On calcule l'int\'egrale 
$\int_0^1 f(t) {\rm d} t$ en r\'esolvant l'\'equation diff\'erentielle ci-dessus avec la condition initiale $y(0)=0$ par la m\'ethode
de Runge Kutta RK4. Cela correspond-t-il \`a une  m\'ethode d'int\'egration num\'erique connue, si oui laquelle? Vous justifierez votre r\'eponse.
\item En utilisant un r\'esultat g\'en\'eral du cours sur l'ordre des m\'ethodes \`a un pas, appliqu\'e \`a une fonction $f(y,t )=f(t)$ ind\'ependante de $y$,
 montrer que pour $q$ quelconque l'on a 
$$
\sum_{i=0}^q \mu_i c_i^l = \frac{1}{l+1} \;\; \text{ pour } \;\; l=0,\ldots, p-1\;.
$$
\item En d\'eduire que la m\'ethode d'int\'egration num\'erique
%
$$
\int_0^1 g(u) {\rm d} u \simeq \sum_{i=0}^q \mu_i g ( c_i) 
$$
%
est d'ordre au moins $p-1$.
\item On rappelle qu'une m\'ethode d'int\'egration num\'erique \'el\'ementaire \`a $m$ points d'interpolation
est toujours d'ordre strictement inf\'erieur \`a $2m$ (voir l'exercice 2 de la feuille 7 de TD/TP). En d\'eduire
une majoration de $p$ en fonction de $q$. 
\end{enumerate}
%
\end{document}


