Cours de la semaine

This commit is contained in:
Némunaire 2012-03-17 06:16:57 +01:00
parent e74c3a6dc8
commit cc8094acc0
10 changed files with 792 additions and 3 deletions

View file

@ -611,8 +611,410 @@ $$P_0(x)=f(x_1)$$
Pour $m=0,\ldots,n-1$~:
$$P_{m+1}(x)=P_m(x)+(x-x_1)\ldots(x-x_{m-1}).f[x_1,\ldots,x_{m+1}]$$
%% Il cherchait à faire un tableau
\paragraph{Exercice 3} Sur $[-1,1]$ on considère la fonction $f(x)=X^3-X^2$
\begin{enumerate}
\item Décomposer $f(x)$ dans la base de Chebyshev.
\item En déduire $P^*_2(x)$ polynôme de meilleur approximation de degré 2 de
$f(x)$.\\
Donner l'erreur d'approximation $||f-P^*||$
$$||f||=max|f(x)|\qquad -1\leq x\leq 1$$
\end{enumerate}
\[
\begin{array}{l l}
f(x)= & a_0.T_0(x)+a_1.T_1(x)+a_2.T_2(x)+a_3.T_3(x)$$\\
\begin{cases}
T_0(x)=1\\
T_1(x)=x\\
T_2(x)=2x^2-1\\
T_3(x)=4x^3-3x
\end{cases}
& a_0 + a_1x + a_2(2x^2-1) + a_3 (4x^3 - 3x)\\
& a_0 + a_1x + 2a_2x^2 - a_2 + 4a_3 + 4a_3x^3 - 3a_3x\\
& a_0 - a_2 + (a_1 - 3a_3) x + 2a_2x^2 + 4a_3 + 4a_3x^3\\
\end{array}
\]
\[
\Rightarrow
\begin{cases}
a_0 - a_2 = 0\\
a_1 - 3a_3 = 0\\
2a_1 = -1\\
4a_3=1
\end{cases}
\Rightarrow
\begin{cases}
a_0=-\frac{1}{2}\\
a_2=-\frac{1}{2}\text{ et }a_1=\frac{3}{4}\\
a_4=\frac{1}{4}
\end{cases}
\]
Donc
$f(x)=-\frac{1}{2}T_0(x)+\frac{3}{4}T_5(x)-\frac{1}{2}T_2(x)+\frac{1}{4}T_3(x)$
$$\epsilon=f(x)-P_2(x)=\frac{1}{4}=T_3(x)$$
$$M=||\epsilon||=\max|\epsilon(x)|=\frac{1}{4}$$
\[
\begin{cases}
T_3(x)=\cos(3\theta)\\
\theta=\arccos(x)
\end{cases}
\]
\subparagraph{Théorème} Si la fonction erreur $\epsilon(x)$ atteint la veuleur
extrême alternée $M$ et $-M$ en au moins $(n+2)$ points, alors
$P_n(x)=P^*_n(x)$.
$$\epsilon(x)=\frac{1}{4}T_3(x)=
\frac{1}{4}\cos(3\theta)\qquad\theta=\arccos(x)$$
Les $x_k=\cos\frac{k\pi}{n}\quad n=3$.
\[
\begin{array}{l l l}
k=0 & x_0=1 & \varepsilon(x_0)=\frac{1}{4}=M\\
k=1 & x_1=\cos\frac{\pi}{3} & \varepsilon(x_1)=-\frac{1}{4}\\
k=2 & x_2=\cos\frac{2\pi}{3} & \varepsilon(x_2)=\frac{1}{4}\\
k=3 & x_3=\cos\pi & \varepsilon(x)=-\frac{1}{4}
\end{array}
\Rightarrow
\begin{cases}
x_0=1\\
x_1=\cos\frac{\pi}{3}=\frac{1}{2}\\
x_2=\cos\frac{2\pi}{3}=-\frac{1}{2}\\
x_3=\cos\pi=-1
\end{cases}
\Rightarrow
P_2(x)=P_2^*(x)
\]
Donc~:
\[
\begin{array}{l l}
P_2^*(x)= & -\frac{1}{2}T_0(x)+\frac{3}{4}T_1(x)-\frac{1}{2}T_2(x)\\
& = -4+\frac{3}{4}x-\frac{1}{4}(2x^3-1)\\
& = -\frac{1}{2}+\frac{3}{4}X - x^2+\frac{1}{2}\\
& = \frac{3}{4}x = \frac{3}{4}x-x^2
\end{array}
\]
\paragraph{Exercice 4} Soit $f(x)$ une fonction continue donnée aux points
$x_i(i=1,2,3,4)$.
$f(0)=-5, f(1)=17, f(2)=115, f(4)=143$
En utilisant le polynôme d'interpolation de Lagrange, interpoler $f(2)$.
$$f(0)=-5, f(1)=17, f(2)=115. f(4)=143$$
$$P_n(x)=\sum_{i=1}^{n+1}f(x_i)$$
$$L_i(x)=\prod^{n+1}_{j=1}\left(\frac{x-x_i}{x_i-x_j}\right)$$
Interpolons $f(2)$~:
$$L_1(2)=\frac{(2-1)(2-3)(2-4)}{(-1)(-3)(-4)}=\frac{-2}{12}=\frac{-1}{6}$$
$$L_2(2)=\frac{(2-0)(2-3)(2-4)}{1(-2)(-3)}=\frac{4}{6}=\frac{2}{3}$$
$$L_3(2)=\frac{2(2-1)(2-4)}{3(3-1)(3-4)}=\frac{-4}{-6}=\frac{2}{3}$$
$$L_4(2)=\frac{2(2-1)(2-3)}{4(3)(1)}=\frac{-2}{12}=\frac{-1}{6}$$
La valeur interpolée est~: $f(2)\simeq P_3(2)=-\frac{1}{6}(-5)+\frac{2}{3}(17)
+\frac{2}{3}(115)-\frac{1}{6}(143)\simeq 65$
\paragraph{Exercice 5} Soit $f\in C^{n+1}[a,b]$, $P_n(x)$ son polynôme
d'interpolation de $f(x)$ aux points $x_i\in[a,b]$ $(i=1\ldots n+1)$
$$\varepsilon(x)=f(x)-P_n(x)\quad\text{erreur d'interpolation}$$
Montrer qu'il existe $\eta_x\in[a,b]$ tel que $\varepsilon(x)=
\prod^{n+1}_{i=1}(x-x_i)\frac{f(\eta_x)}{(n+1)!}$
\subparagraph{Rappel} Théorème de Rolle
$f$ continue, dérivable sur $[a,b]$. Si $f(a)=f(b)$ alors $\exists C\in]a,b]$
tel que $f'(c)=0$.
$$\varepsilon(x_i)=0\text{ car }P_n(x_i)=f(x_i)$$
$$\varepsilon(x)=\prod^{n+1}_{i=1}(x-x_i)A_x\qquad x.f(x_i)\quad x\neq
x_i\quad A_x\text{ quantité inconnue}$$
Soit la fonction
$\varphi(y)=f(y)-P_n(y)-\prod^{n+1}_{i=1}(y-x_i)A_x\quad\forall y\in[a,b]$
\[
\begin{cases}
\varphi(x_i)=0\qquad\forall i=1\ldots n+1\\
\varphi(x)=0
\end{cases}
\]
$\varphi$ admet $(n+2)$ racines. $\varphi\in C^{n+1}$ car $f\in C^{n+1}$.
On applique le théorème de Rolle $\varphi'$ admet au moins $(n+1)$ racines.
On applique le théorème de Rolle $\varphi''$ admet au moins $n$ racines.
On applique le théorème de Rolle $\varphi^(3)$ admet au moins $n-1$ racines.
\hspace{6cm}\vdots
On applique le théorème de Rolle $\varphi^(p)$ admet au moins $(n+1)-p$
racines.
Si $p=n+1$ $\Rightarrow\varphi^{(n+1)}$ admet au moins 1 racine
$\eta_x\quad\varphi^{(n+1)}(\eta_x)=0$
or~:
$$\varphi^{(n+1)}(y)=f^{(n+1)}(y)-0-(n+1)! A_x$$
$$\varphi^{(n+1)}(\eta_x)=0\Rightarrow f^{(n+1)}(\eta_x)-(n+1)!
A_x=0\Rightarrow A_x=\frac{f^{(n+1)}(\eta_x)}{(n+1)!}$$
Donc~:
$$\varepsilon(x)=\prod^{n+1}_{i=1}(x-x_i)\frac{f^{(n+1)(\eta_x)}}{(n+1)!}$$
$$f\left[x_1, \ldots, x_{n+1}\right]=\frac{f\left[x_1,\ldots,
x_n\right]-f\left[x_2,\ldots, x_{n+1}\right]}{x_1 - x_{n+1}}$$
\section{Dérivation numérique}
Lorsqu'une fonction $f(x)$ dont on veut calculer les dérivées siccessives
$f'(x)$, $f''(x),\ldots$ est données par les valeurs $f(x_i)$, on procède de la
manière suivante~:
\begin{enumerate}[(i)]
\item $f(x)\simeq g(x)\qquad\text{(polynôme d'interpolation)}$.
\item $f^(k)(x)\simeq g^{(k)}\simeq g^{(k)}(x)$
\end{enumerate}
\subsection{Dérivée d'un polynôme par division synthétique}
$$f(x)\simeq P_n(x)=\sum^{n}_{j=0}a_j.x^j$$
Soit à calculer~:
$$f^{(k)}(x=t)\qquad\text{dérivée d'ordre }k$$
$$f^{(k)}\simeq P^{(k)}_n(x)=\sum^{n}_{j=k}j(j-1)\ldots(j-k+1)a_j.x^{j-k}$$
Il suffit alors de calculer $P_n^{(k)}(x=t)$.
$$P_n(x) | x-t\qquad\text{(division par )x-t}$$
$$P_n(x)=(x-t)Q_{n-1}(x)+R_n\quad Q_{n-1}(x)=b_{n-1}X^{n+1}+\ldots+b+0$$
\[
\begin{cases}
b_{n-1}=a_n\\
b_{j-1}=t.b_j+a_j\quad j=n-1,\ldots, 1
\end{cases}
\]
\paragraph{Exemple} $$P_4(x)=2x^4-5x^3+x^2-7x+6\qquad t=2$$
\begin{figure}[h]
\centering
\begin{tabular}{c|cccccc}
$j$ & 4 & 3 & 2 & 1 & 0 &\\\hline
$P_4 a_j$ & 2 & -5 & 1 & -7 & 6 &\\
$Q_3 b_j$ & & 2 & -1 & -1 & -9 & -12\\
\end{tabular}
\end{figure}
On peut poursuivre le processus~:
\begin{figure}[h]
\centering
\begin{tabular}{c|cccccc}
$j$ & 4 & 3 & 2 & 1 & 0 & Reste\\\hline
$P_4\quad a_j$ & 2 & -5 & 1 & -7 & 6 &\\
$Q_3\quad b_j$ & & 2 & -1 & -1 & -9 & $-12 = R_4$\\
$Q_2\quad c_j$ & & & 2 & 3 & 5 & $1 = R_3$\\
$Q_1\quad d_j$ & & & & 2 & 7 & $19 = R_2$\\
$Q_0\quad e_j$ & & & & & 2 & $11 = R_1$\\
\end{tabular}
\end{figure}
$$R_0=a_n=2=Q_0$$
\[
\left.
\begin{array}{l}
P_4(x)=(x-2).Q_3(x)-12\\
Q_3(x)=(x-2).Q_2(x)+1\\
Q_2(x)=(x-2).Q_1(x)+19\\
Q_1(x)=(x-2).Q_0(x)+11\\
\end{array}
\right\}
P_4(x)=-12+(x-2)+19(x-2)^2+11(x-2)+2(x-2)^4
\]
Ainsi les $R_k$ sont les coefficients du développement de $P_n(x)$ en puissance
de $(x-t)\quad(t=2)$.
D'une manière générale~:
$$P_n(x)=\sum^n_{l=0}R_e(x-t)^{n-l}$$
Em dérivant $k$ fois~:
$$P_n^{(k)}(x)=\sum^{n-k}_{l=0}(n-l)\ldots(n-l-k+1)R_e(x-t)^{n-l-k}$$
$$x=t\Rightarrow P_n^{(k)}(t)=k!R_{n-k}$$
$$f^{(k)}(t)\simeq k!R_{n-k}$$
\paragraph{Exercice 6} $f(x)$ continue, donnée aux points~: $f(-2)=69$,
$f(2)=17$, $f(3)=24$, $f(5)=-22$
\begin{enumerate}
\item Construire le tableau des différences divisées.
\item En déduire le polynôme d'interpolation de Newton.
\item En utilisant l'algorithme de division synthétique, interpoler $f(4)$ et
donner les approximations des dérivées $f'(4)$, $f''(4)$, $f'''(4)$.
\end{enumerate}
\begin{figure}[h]
\centering
\begin{tabular}{c|cccc}
$x_i$ & 0 & 1 & 2 & 3 \\\hline
$-2$ & \textbf{69} & \textbf{-13} & & \\
$2$ & 17 & & \textbf{4} & \textbf{-2}\\
$3$ & 24 & 7 & -10 & \\
$5$ & -22 & -23 & & \\
\end{tabular}
\end{figure}
Pour calculer les machins de ce tableau~:
$$f[x_i;x_j]=\frac{f[x_i]-f[x_j]}{x_j-x_i}=\frac{f(x_i)-f(x_j)}{x_i-x_j}$$
$$f[-2,2,3]=\frac{f[2,3]-f[-2,2]}{5}$$
\subparagraph{2. Le polynôme d'interpolation de Newton}
$$P_3(x)=f[-2]+(x+2)f[-2,2]+(x+2)(x-2)f[-2,2,3]+(x+2)(x-2)(x-3)f[-2,2,3,5]$$
On remplace ensuite les coefficients~:
$$P_3(x)=69-13(x+2)+4(x+2)(x-2)-2(x+2)(x-2)(x-3)$$
On développe tout pour avoir le polynôme dans sa forme canonique (utile pour la
question suivante)~:
$$P_3(x)=-2x^3+10x^2-5x+3$$
On peut vérifier que ce polynôme interpole bien les points donnés de $f$.
\subparagraph{3. Algorithme de division synthétique}
Pour $Q_1$,0~: 4 (coeff) * -2 (valeur précédente sur la ligne) + 2 (valeur au dessus) = -6
\begin{figure}[h]
\centering
\begin{tabular}{c|ccccc}
$j$ & 3 & 2 & 1 & 0 & Reste\\\hline
$P_3\quad a_j$ & -2 & 10 & -5 & 3 & \\
$Q_2\quad b_j$ & & -2 & 2 & 3 & $15 = R_3$\\
$Q_1\quad c_j$ & & & -2 & -6 & $-21 = R_2$\\
$Q_0\quad d_j$ & & & & -2 & $-14 = R_1$\\
\end{tabular}
\end{figure}
$$P_3(x)=15-21(x-4)-14(x-4)^2-2(x-4)^3$$
$$f(4)\simeq P_3(4)=15$$
$$f^{(k)}(t)\simeq k!_{n-k}$$
$$f'(4)\simeq 1!R_2=-21$$
$$f''(4)\simeq 2!R_1=-28$$
$$f'''(4)\simeq 3!R_0=-12$$
\paragraph{Exercice 7} Soit $f(x)$ continue, donnée aux points $f(-2)=32$,
$f(-1)=26$, $f(0)=30$, $f(1)=28$
\begin{enumerate}
\item Construire le tableau des différences divisées.
\item En déduire le polynôme d'interpolation de Newton.
\item En utilisant l'algorithme de division synthétique, interpoler $f(2)$ et
donner les approximations des dérivées $f'(2)$, $f''(2)$, $f'''(2)$.
\end{enumerate}
\begin{figure}[h]
\centering
\begin{tabular}{c|cccc}
$x_i$ & 0 & 1 & 2 & 3 \\\hline
$-2$ & \textbf{-32} & \textbf{58} & & \\
$-1$ & 26 & & \textbf{27} & \textbf{8}\\
$0$ & 30 & 4 & -3 & \\
$1$ & 28 & -2 & & \\
\end{tabular}
\end{figure}
\subparagraph{2. Le polynôme d'interpolation de Newton}
$$P_3(x)=f[-2]+(x+2)f[-2,-1]+(x+2)(x+1)f[-2,-1,0]+(x+2)(x+1)(x)f[-2,-1,0]$$
On remplace ensuite les coefficients~:
$$P_3(x)=-32-58(x+2)-27(x+2)(x+1)-8(x+2)(x+1)(x)$$
On développe tout pour avoir le polynôme dans sa forme canonique (utile pour la
question suivante)~:
$$P_3(x)=8x^3-3x^2-7x+30$$
On peut vérifier que ce polynôme interpole bien les points donnés de $f$.
\begin{figure}[h]
\centering
\begin{tabular}{c|ccccc}
$j$ & 3 & 2 & 1 & 0 & Reste\\\hline
$P_3\quad a_j$ & 8 & -3 & -7 & 30 & \\
$Q_2\quad b_j$ & & 8 & 13 & 19 & $68 = R_3$\\
$Q_1\quad c_j$ & & & 8 & 29 & $77 = R_2$\\
$Q_0\quad d_j$ & & & & 8 & $45 = R_1$\\
\end{tabular}
\end{figure}
$$P_3(x)=68+77(x-2)+45(x-2)^2+8(x-2)^3$$
$$f(2)\simeq P_3(2)=68$$
$$f^{(k)}(t)\simeq k!_{n-k}$$
$$f'(2)\simeq 1!R_2=77$$
$$f''(2)\simeq 2!R_1=90$$
$$f'''(2)\simeq 3!R_0=48$$
\section{Intégration numérique}
On se propose d'étudier quelques méthodes numériques permettant d'approcher
$\int^b_af(x)dx$. De telles méthodes s'imposent en particulier lorsque la
primitive $F(x)$ de $f(x)$ n'est pas connue, ou dans le cas où $f(x)$ n'est
donnée que par points.
Par exemple pour~: $\int^b_a e^{-x^2}dx\simeq\sum^n_{i=0}\omega_i f(x_i)$
\subsection{Méthode générale}
$$I=\int^b_a f(x)p(x)dx\quad p(x):\text{fonction poids, }>0$$
L'idée est d'utiliser les valeurs de $f$ aux points $x_i\in[a,b]\quad i=0\ldots
n$.
$$\int^b_a f(x)p(x)dx=\sum^a_{i=0}\omega_i.f(x_i)+E\quad E\text{ l'erreur}$$
Les coefficients $\omega_i$ sont déterminés de telle sorte que $E$ soit nulle
lorsque $f(x)$ appartient à une certaine classe de fonctions (les polynômes de
degré $\le N$).
$$f(x)=\phi(x)+E(x)$$
$$\phi(x_i)=f(x_i)\quad i=0\ldots n$$
$$\int^b_a f(x)p(x)dx=\int^b_a\phi(x)p(x)dx+\int^b_a E(x)p(x)dx$$
$$\phi(x)=\sum^n_{j=0}L_j(x)f(x_j)\quad\text{Lagrange}$$
$$\int^b_a f(x)p(x)dx=\int^b_a\left.\sum^n_{j=0}L_j(x)f(x_j)p(x)dx\right.+
\int^b_a E(x)p(x)dx$$
$$=\sum^n_{j=0}\left(\int^b_a L_j(x)p(x)dx\right)f(x_j)+\int^b_a E(x)p(x)dx$$
$$\int^b_a f(x)p(x)dx=\sum^n_{j=0}\omega_jf(x_j)+E$$
$$\text{avec }\omega_j=\int^b_a L_j(x_j)p(x_j)dx\text{ et
}E=\int^b_aE(x)p(x)dx$$
\subsection{Quelques exemples de méthodes composées}
$$[\alpha,\beta]\qquad I=\int^\beta_\alpha f(x)dx$$
La méthode consiste à décomposer $[\alpha,\beta]$ en $k$ intervalles
$[\alpha_i,\alpha_{i+1}]\quad i=0,1,\ldots k$ et ensuite à approcher chaque
$\int^{\alpha_{i+1}}_{\alpha_i}f(x)dx$ en remplaçant $f(x)$ par son polynôme
d'interpolation.
\paragraph{Exemple} $t_i\in[\alpha_i,\alpha_{i+1}]$
$$f(x)\simeq P_0(x)=f(t_i)\quad\forall x\in\left[\alpha_i,\alpha_{i+1}\right]$$
$$\int^{\alpha_{i+1}}_{\alpha_i}f(x)dx\simeq(\alpha_{i+1}-\alpha_i)f(t_i)$$
$$\int^\beta_\alpha
f(x)dx\simeq\sum^{k-1}_{i=0}\left(\alpha_{i+1}-\alpha_i\right)f(t_i)\quad\text{Somme
de Riemann}$$

View file

@ -1,5 +1,5 @@
\title{Mathématiques}
\author{Le prof}
\title{Approximation, interpolation, optimisation}
\author{Le prof de maths}
\date{ING1}
\maketitle