PLAN
Introduction
Développement
Conclusion
Annexes
INTRODUCTION
La programmation est l’ensemble des activités qui
permettent l’écriture des programmes informatiques. L’écriture
d’un programme se fait dans un langage de
programmation .Ainsi, comme langage, nous utilisons ici
MATLAB (Matrix LABoratory) qui permet de faire du calcul
matriciel, de développer et d’exécuter des algorithmes, de
créer des interfaces utilisateurs (IU) et de visualiser des
données. Grâce à ce langage, nous traiterons donc dans la
suite, différents exercices basés sur les notions de fonctions,
tracé de courbes, factorielles, série de Fourier et autres qui
nous permettront de mettre en exergue les différentes
fonctionnalités de ce langage.
Développement
Exercice 1.1 : Tracé de fonctions (tracefoncs)
Tracer les courbes des fonctions suivantes (quand 2 courbes
sont indiquées, les superposer). On écrira un fichier de
commandes contenant toutes les instructions plot séparées par
l’instruction pause.
– y = √(x 2−x ) et y = x − 1/2 sur [1, 3]
– y = sin(sin x) sur [0, 2π]
– y = J0(x) et y = J1(x) (J0 et J1 fonctions de Bessel de première
espèce) sur [0, 10π]
– y = J0(x) et y = sin (x + π/4)/ √ x sur [π/4, 10π]
– y = erf(x) (fonction erreur) sur [−3, 3]
Exercice 1.2 : Calcul de factorielle (facto)
Ecrire un programme qui calcule n! pour une valeur de n que
l’on définira.
Pour une valeur de n=50, après exécution du programme, on
obtient
>> facto
Entrez une valeur 50
Son factoriel est:
3.0414e+64
Exercice 1.3 : Calcul d’une série (calserie)
Calculer la série:
+∞
S= ∑ n1²
n =1
On tronquera la sommation lorsque le terme général sera
inférieur à 10-8 fois la somme partielle. Comparer le résultat
numérique au résultat exact (= π²/6 comme chacun sait :-).
Suite à l’exécution du programme pour n=1 et S=0 au départ,
on obtient le résultat ci-dessous :
>> calserie
n
7798
Tn
1.6445e-08
S
1.6448
Exercice 1.4 : Fonctions définies par morceaux (tracefoncs2)
Tracer les courbes des fonctions suivantes (il n’est pas
obligatoire d’´ecrire un fichier de fonction).
y = 1 + x, pour x ∈ [−1, 0],
y= y = 1 − x, pour x ∈ [0, 1],
y = 0 sinon
y = −(x + 1)² , pour x < −1
, pour x ∈ [−1, 0],
3/2
y= y = (x + 1)
y = x² + 2 x + 1 pour x ∈ [0, 2],
3
y = 11 x − 3
2
pour x > 2
Exercice 1.5 : Ecriture d’une fonction (g et testg)
– Ecrire une fonction appelée g, qui traduit la fonction
mathématique g(x) = x log(x 2 + 1).
– Utiliser cette fonction pour évaluer g(1).
– Tracer la fonction g sur [0, 10]
Suite à l’écriture de la fonction g dans MATLAB, l’on utilise ce
dernier pour évaluer g(1) et tracer la fonction g sur [0, 10] :
>> testg
g(1)
0.3010
Exercice 1.6 : Ecriture d’une fonction n! (ffact)
Ecrire une fonction ffact qui reçoit un argument en entrée et
calcule sa factorielle.
Pour une valeur de n=50, après exécution du programme, on
obtient :
>> fonctionfact
x
3.0414e+64
Exercice 1.7 : Courbes d´définies en polaires (courbespol)
Avec la fonction de conversion pol2cart et l’instruction plot,
tracer les courbes définies en polaires par :
– r = 1 + e cos θ
1
– r = (1+e cos θ)
–r=1
On prendra pour e une valeur entre 0 et 1. Reprendre
l’exercice avec la commande polar (on regardera l’aide en
ligne).
Avec la fonction de conversion pol2cart et l’instruction plot, on
obtient les courbes suivantes :
Pour e=0, on :
Pour e=1, on a :
Après reprise de
l’exercice avec la
commande
polar on obtient :
Pour e=0, on a :
Pour e=1, on a :
Exercice 1.8 : Illustration des séries de Fourier (fourier)
On considère la fonction sur l’intervalle [0, 2π] :
1 si 0 < x ≤ π
f(x)= (A)
0 sinon
1. Ecrire une fonction MATLAB carre codant cette fonction
mathématique. La fonction doit pouvoir s’appliquer
indifféremment `a un scalaire ou terme à terme à un tableau de
dimension quelconque. Ecrire un programme principal pour
tracer la fonction.
2. On peut montrer que cette fonction f(x) coïncide sur [0, 2π]
avec la série de Fourier :
+∞ n
1−(−1)
g(x) = ½+ ∑ n²
sin (nx) (B)
n =1
Ecrire un programme MATLAB qui construit progressivement la
série, et qui trace la fonction obtenue à chaque ajout d’un
nouveau terme. Vérifiez que l’on retrouve bien, pour un grand
nombre de termes, la fonction f(x).
Suite à l’écrire de la fonction MATLAB carre codant cette
fonction mathématique (A) et un programme principal pour
tracer la fonction on obtient ;
Le programme MATLAB qui construit progressivement la série
(B), et qui trace la fonction obtenue à chaque ajout d’un
nouveau terme permet d’obtenir la courbe suivante :
Exercice 1.9 : Utilisation de global (vdw et testvdw)
On souhaite écrire une fonction représentant l’équation d’état
de Van der Waals d’un gaz donné :
p=
Mais on impose que la fonction n’ait que 2 arguments d’entrée
v et T, et 1 argument de sortie p, les autres grandeurs a, b et R
étant considérées comme des constantes (on donne R = 8.32, a
= 0.5547 et b = 3.0516 × 10−5 en unités S.I.).
1. Trouvez une solution avec l’aide du polycopié et de
l’intervenant, et écrivez la fonction.
2. Ecrire un programme qui trace p en fonction de v pour T =
200K, T = 300K et T = 1000K, sur un même graphique. On
prendra v variant entre 1.1b et 10b.
En remplaçant R, a, et b par leur valeurs respectives et en
prenant T = 200K et v= 3.4x10−4, on obtient :
>> vdw
p
5.8903e+05
Pour des variations de T et v, on obtient :
Exercice 1.10 : Ecriture d’une fonction de plusieurs variables
Pour calculer la capacité calorifique Cp d’un gaz en fonction de
la température T (en Kelvin), on dispose de corrélations du type
:
Il est souvent intéressant de calculer également la primitive de
cette fonction, qui s’exprime par :
1. Ecrire dans un fichier une fonction MATLAB appelée cpip qui
reçoit en entrée la température T et les coefficients A, B, C, D,
E, et calcule Cp et Ip en sortie. On spécifiera le nom du fichier
utilisé.
précédente pour tracer Cp(T) et Ip(T) sur l’intervalle T ∈ [250 K,
2. Ecrire un programme MATLAB qui utilise la fonction
3000 K].
3. Pour des raisons pratiques, on veut que la fonction n’ait
comme argument d’entrée que la température T. Comment
faire ? Réécrivez la fonction et le programme.
En prenant pour valeur T=250, A=10, B=12, C=14, D=15,
E=11, on obtient :
>> resutatcpip
Cp =
36.9778
Ip =
1.7507e+03
Le tracé de Cp(T) et Ip(T) sur l’intervalle T ∈ [250 K, 3000 K]
donne l le graphe ci-dessous :
Exercice 1.11 : Projection de Mercator
Pour tracer des cartes géographiques, on utilise la projection
de Mercator qui à un point d’une surface sphérique fait
correspondre un point du plan. La projection est donnée par :
x=θ
π φ
y = log ([tan( 4 + 2 )]
où θ ∈ [0˚, 360˚] est la longitude et φ ∈ [−90˚, 90˚] la latitude.
1. Ecrire une fonction MATLAB effectuant la projection, qui renvoie les coordonnées (x, y) sur la carte
en fonction des deux coordonnées géographiques θ et φ.
2. On vous fournit un fichier texte de données géographiques [Link], organisé en deux colonnes
(latitude, longitude), exprimées en degrés, d’écrivant les frontières des continents. Ecrire un
programme utilisant la fonction précédente pour tracer la projection de Mercator des continents.
Note : on utilisera la fonction plot pour tracer point par point, sinon ça a une drôle de tête. . .