Séance 1-Calcul Scientifique

Université Moulay Ismaïl
Faculté des Sciences et Techniques; Errachidia

Département de physique
Calcul scientifique
ENERGIES RENOUVELABLES
Option: Technologies solaires et éoliennes
Groupe facebook: Calcul scientifique s5
Pr. SARA TEIDJ
2018/2019
Ce cours constitue une introduction au Calcul
Scientifique. Son objectif est de présenter des
méthodes numériques permettant de résoudre
avec un ordinateur des problèmes mathématiques qui
ne peuvent pas être traités simplement avec une
feuille et un stylo
2018/2019
Déroulement du cours
 Cours magistral
TP ( Matlab)
 TD
2018/2019
Plan du cours
 Initiation à Matlab.
 Les concepts rattachés au calcul

numérique.
 Résolution des équations non linéaires.
 Résolution des systèmes linéaires.
2018/2019
Initiation à Matlab
2018/2019
Sommaire
 Présentation de Matlab
 Vecteurs et Matrices
 La Programmation Matlab
 Les Scripts
 Les Graphiques
2018/2019
Introduction
Ce chapitre est une initiation au logiciel Matlab et à la

programmation dans cet environnement.
On y exposera les bases de Matlab et de la programmation
d’algorithmes mathématiques élémentaires à l’aide de ce
langage.
2018/2019
Présentation de Matlab
La figure 1 présente l’écran Matlab de base: la fenêtre de commandes.

Le symbole >> s’appelle invite de Matlab. Il invite l’utilisateur à
taper une commande.
2018/2019
Dans la fenêtre de commande, l’utilisateur peut affecter des

valeurs à des variables et effectuer des opérations sur
celles-ci. Par exemple :
>> x=4
x=
4
>> y=2
y=
2
>> x+y
ans =
6
>> x*y
ans =
8
>> 2018/2019
Ici, il faut noter que lorsque l’utilisateur ne fixe pas de

variable de sortie, Matlab place par défaut le résultat d’une
opération dans la variable ans. Il est toujours possible de
connaître les variables utilisées et leur type à l’aide de la
commande who ou bien whos. Par exemple, pour les
manipulations précédentes:
>> whos
Name Size Bytes Class
ans 1x1 8 double array
x 1x1 8 double array
y 1x1 8 double array
Grand total is 3 elements using 24 bytes
>>
2018/2019
La solution de x+y a donc été perdue. Il est donc

préférable de toujours donner des noms aux variables de
sortie :
>> x=4;
>> y=2;
>> a=x+y
a=
6
2018/2019
>> b=x*y
b=
8
>> whos
a 1x1 8 double array
b 1x1 8 double array
>>
2018/2019
Notons au passage que le point-virgule permet de ne pas

afficher la valeur à l’écran, ce qui permettra éventuellement
des programmes plus rapides.
Le signe de pourcentage (%) permet de mettre ce qui suit
sur une ligne en commentaire
(Matlab n’en tiendra pas compte à l’exécution).
2018/2019
La fonction clear permet d’effacer des variables. Par exemple :

>> clear x % on efface x de la mémoire
>> whos
>>
2018/2019
La sortie de la fonction whos donne, entre autre, la classe de la

variable. Plusieurs classes de variables sont disponibles à
l’utilisateur de Matlab. Les classes les plus utiles pour
l’utilisateur débutant sont l’entier, le réel simple, le réel
double et les variables caractère ‘char’.
Pour les variables char, la déclaration se fait entre
apostrophes:
>> mot1 = ’bonjour’
mot1 = bonjour
2018/2019
Il est possible de concaténer des mots à l’aide des

parenthèses carrées (crochets) (la fonction strcat de Matlab
permet d’effectuer sensiblement la même tâche) :
>> mot1 = ’bonjour’;
>> mot2 = ’tout le monde’;
>> mot1_2 = [mot1 ’ ’ mot2] % l’emploi de ’ ’ permet
d’introduire un espace
mot1_2 = bonjour tout le monde.
2018/2019
•Les formats d'affichage des réels

format long % format long à 15 chiffres.
format short e % format court à 5 chiffres avec notation
en virgule flottante.
format long e %format long à 15 chiffres avec notation
en virgule flottante.
•Affichage simple
disp % permet d'afficher un tableau de valeurs numériques
ou de caractères.
Num2str% pour convertir une valeur numérique en une
chaîne de caractères
Input % permet de demander à l'utilisateur d'un
programme de fournir des données( La syntaxe est
var = input(' une phrase ').
2018/2019
Opérations mathématiques de base avec Matlab:
scalaires, vecteurs et matrices
L’élément de base de Matlab est la matrice. C’est-à-dire qu’un scalaire

est une matrice de dimension 1x1, un vecteur colonne de dimension
n est une matrice nx1, un vecteur ligne de dimension n, une
matrice 1xn.
Les scalaires se déclarent directement, par exemple :
>> x = 0;
>> a = x;
2018/2019
Vecteurs:
Les vecteurs ligne se déclarent de la manière suivante :
>> V_ligne = [0 1 2]
V_ligne =
012
Pour les vecteurs colonne, on sépare les éléments par des
points-virgules :
>> V_colonne = [0; 1; 2]
V_colonne =
0
1
2
2018/2019
Il est possible de transposer un vecteur à l’aide de la fonction
transpose ou avec le point apostrophe (.’).Ainsi,
>> V_colonne=transpose(V_ligne)
V_colonne =
0
1
2
>> V_colonne=V_ligne.’
V_colonne =
0
1
2
2018/2019
Le double point (:) est l’opérateur d’incrémentation dans

Matlab. Ainsi, pour créer un vecteur ligne des valeurs de 0 à 1
par incrément de 0.2, il suffit d’utiliser:
>> V= [0:0.2:1]
V=
Columns 1 through 6
0 0.2000 0.4000 0.6000 0.8000 1.0000
2018/2019
Par défaut, l’incrément est de 1. Ainsi, pour créer un vecteur
ligne des valeurs de 0 à 5 par incrément de 1, il suffit
d’utiliser :
>> V=[0:5]
V=
012345
On peut accéder à un élément d’un vecteur et même
modifier celui-ci directement il n’y a pas d’indice 0 dans
les vecteurs et matrices en Matlab) :
>> a=V(2);
>> V(3)=3*a
V=
016345
2018/2019
Les opérations usuelles d’addition, de soustraction et de

multiplication par scalaire sur les vecteurs sont définies
dans MATLAB :
>> V1=[1 2];
>> V2=[3 4];
>> V=V1+V2 % addition de vecteurs
V=
4 6
>> V=V2-V1 % soustraction de vecteurs
V=
2 2
2018/2019
>> V=2*V1 % multiplication par un scalaire

V=
2 4
Dans le cas de la multiplication et de la division, il faut faire
attention aux dimensions des vecteurs en cause.
Pour la multiplication et la division élément par élément, on
ajoute un point devant l’opérateur (.* et ./). Par exemple :
>> V=V1.*V2 % multiplication élément par élément
V=
3 8
2018/2019
>> V=V1./V2 % division élément par élément

V=
0.3333 0.5000
Cependant, Matlab lance une erreur lorsque les dimensions ne
concordent pas. Les messages d’erreur sont utiles pour corriger les
programmes (parenthèse oublié par exemple). Il faut cependant
procéder à la vérification systématique de l’instruction ou du
programme avant de lancer l’exécution (reflexe de base d’un
programmeur):
2018/2019
>> V3=[1 2 3]
V3 =
123
>> V=V1.*V3
??? Error using ==> .* Matrix dimensions must agree.
La multiplication de deux vecteurs est donnée par (*). Ici,
l’ordre a de l’importance :
>> V1=[1 2]; % vecteur 1x2
>> V2=V1.’; % vecteur 2x1
2018/2019
>> V=V1*V2
V=
5
>> V=V2*V1
V=
12
24
Il est aussi possible de concaténer des vecteurs. Par exemple :
>> V1=[1 2];
>> V2=[3 4];
>> V=[V1 V2]
V=
1234
2018/2019
De même, pour les vecteurs colonnes :
>> V1=[1;2];
>> V2=[3;4];
>> V=[V1;V2]
V=
1
2
3
4
2018/2019
 Matrices
On peut aussi créer des matrices à partir de vecteurs, par
exemple,
>> V1=[1 2];
>> V2=[3 4];
>> V=[V1;V2]
V=
12
34
qui n’est pas équivalent à :
>> V1=[1; 2];
>> V2=[3; 4];
2018/2019
>> V=[V1 V2]
V=
13
24
Il faut donc être très prudent dans la manipulation des
vecteurs. Par exemple, une mauvaise concaténation :
>> V1=[1 2];
>> V2=[3; 4];
>> V=[V1;V2]
??? Error using ==> vertcat All rows in the bracketed
expression must have the same number of columns.
2018/2019
Les matrices peuvent aussi être construites directement :

>> M=[1 2; 3 4]
M=
12
34
On peut évidemment avoir accès aux éléments de la matrice
par :
>> m21=M(2,1) % 2eme ligne, 1ere colonne
m21 =
3
2018/2019
On peut aussi ”compter” les éléments. Matlab compte alors

tous les éléments d’une colonne (de haut en bas) avant
d’accéder à la colonne suivante.
Ainsi, dans la matrice 3x3 suivante:
>> A=[1 2 3; 8 5 6;7 8 9]
A=
123
856
789
2018/2019
les valeurs des éléments ai,j sont données par leur rang
affecté par Matlab. Le 4e élément est 2 :
>> a4=A(4)
a4 =
2
Il est aussi possible de stocker dans un vecteur une ou
plusieurs lignes (ou colonnes).Ainsi, si l’on veut stocker la
deuxième colonne de la matrice A :
>> V=A(:,2) % ici, (:) signifie toutes les lignes
V=
2
5
8
2018/2019
De la même manière, si l’on veut stocker les lignes 2 et 3
>> M2=A(2:3,:) % (2:3) signifie ligne 2 à 3
% et (:) signifie toutes les colonnes
M2 =
856
789
Il est possible d’inverser inv(), de transposer transpose() ou
avec l’apostrophe (.’) les matrices :
>> invM=inv(M)
invM =
-2.0000 1.0000
1.5000 -0.5000
2018/2019
>> transpM=M.’
transpM =
13
24
Un des intérêts de Matlab est la possibilité d’utiliser
directement les opérations mathématiques prédéfinies pour les
matrices. L’addition et la soustraction sont directes
(attention aux dimensions) ainsi que la multiplication par un
scalaire :
>> A=[1 2; 3 4];
>> B=[4 3; 2 1];
>> C=A+B % addition
2018/2019
C=
5 5
5 5
>> D=A-B % soustraction
D=
-3 -1
1 3
>> C=3*A % multiplication par un scalaire
C=
3 6
9 12
2018/2019
Pour la multiplication et la division, les opérateurs usuels (* et

/) sont définis pour la multiplication et division matricielles :
>> C=A*B % multiplication de matrices
C=
8 5
20 13
2018/2019
>> D=A/B % division de matrices
D=
1.5000 -2.5000
2.5000 -3.5000
Afin de réaliser la multiplication et la division élément par
élément, on précède les opérateurs par un point (.* et ./) :
>> C=A.*B % multiplication élément par élément
C=
4 6
6 4
2018/2019
>> D=A./B % division élément par élément
D=
0.2500 0.6667
1.5000 4.0000
D’autres opérations sur les matrices seront présentées
dans les sections subséquentes.
Il faut noter certaines matrices spéciales qui peuvent être
utilisées, par exemple la matrice identité :
2018/2019
>> I=eye(3) % matrice identité
I=
1 0 0
0 1 0
0 0 1
On peut aussi déclarer des vecteurs (et des matrices) ne
contenant que des zéros ou des 1.
rand(p) % crée une matrice p × p dont les éléments sont
aléatoires, distribués selon une loi de probabilité uniforme dans
l’intervalle [0, 1[.
>> rand(2)
ans =
0.8147 0.1270
0.9058 0.9134 2018/2019
>> V_nul=zeros(1,2) % un vecteur de 1 ligne, 2 colonnes de 0
V_nul=
00
>> V_un=ones(1,2) % un vecteur de 1 ligne, 2 colonnes de 1
V_un=
11
>> M_un=ones(2,2) % une matrice 2x2 de 1
M_un=
1 1
1 1
2018/2019
Dans certaines applications, il est parfois utile de connaître les

dimensions d’une matrice, et la longueur d’un vecteur
(retournés, par exemple, par une fonction).
Dans ce cas, on utilise les fonctions length et size.
>> V=[0:0.1:10]; % utilisation de length - vecteur 1x101
>> n=length(V)
n=
101
2018/2019
>> M=[1 2 3; 4 5 6]; % utilisation de size - matrice 2x3
>> [n,m]=size(M)
n=
2
m=
3
>> dim=length(M) % utilisation de length sur une matrice
dim=
3
Dans ce cas length donne la plus grande dimension, ici le
nombre de colonnes.
2018/2019
A=[1 2 3 ; 2 4 5 ; 6 7 8 ];
det(A) % calcule le déterminant de A

>> det(A)
ans =
-5.0000
abs(x) % la valeur absolue
>> B=[-5 -2; 2 1];
>> abs(B)
ans =
5 2
2 1
2018/2019
 Extraction d’une sous-matrice

On peut utiliser les deux points pour extraire une sous-
matrice d’une matrice A.
 A(:,j) extrait la jème colonne de A. On considère
successivement toutes les lignes de A et on choisit le
jème élément de chaque ligne.
 A(i,:) extrait la ième ligne de A.
 A(:) reforme le matrice A en un seul vecteur colonne en
concaténant toutes les colonnes de A
2018/2019
 A(j:k) extrait les éléments j à k de A et les stocke dans un

vecteur ligne
 A(:,j:k) extrait la sous-matrice de A formée des colonnes j à k.
 A(j:k,:) extrait la sous-matrice de A formée des lignes j à k.
 A(j:k,q:r) extrait la sous-matrice de A formée des éléments
situés dans les lignes j à k et dans les colonnes q à r.
 Ces définitions peuvent s’étendre à des pas d’incrémentation
des lignes et des colonnes différents de 1.
2018/2019
Racines d'un polynôme
Exemple1 :
3x² - 5x + 2 = 0
On commence par définir un " vecteur " qui contient les
coefficients du polynôme :
p = [ 3 -5 2 ]
p=
3 -5 2
roots(p) :trouve les racines d'un polynôme
ans =
1.0000
0.6667
2018/2019
Détermination des coefficients d’un polynôme à partir ses
racines
exemple2 :
a=[2 1]
a=
2 1
>> poly(a): trouve le polynôme à partir des ses racines
ans =
1 -3 2
(c’est-à-dire : x² -3x +2)
2018/2019
Exemple 3
>> p = [ 3 -5 2 ]
p=
3 -5 2
Calcul de f( x = 1) :
>> polyval(p,1) : évalue le polynôme
ans =
0
2018/2019
Fichiers script et function
Jusqu’à présent, l’utilisation que nous avons faite de Matlab

s’apparente beaucoup à celle d’une calculatrice. Pour des
tâches répétitives, il s’avère beaucoup plus pratique et
judicieux d’écrire des programmes pour effectuer les
calculs désirés.
Il existe deux types de fichiers qui peuvent être
programmés avec Matlab: les fichiers script (M-file) et
function. Dans les deux cas, il faut lancer l’éditeur de
fichier et sauvegarder le fichier avec l’extension .m.
2018/2019
Fichiers script
Comme tout langage, Matlab possède aussi un certain nombre
d’instructions syntaxiques (boucles simples, conditionnelles,
etc...) et de commandes élémentaires (lecture, écriture,
etc...). Ces instructions syntaxiques seront vues dans la
deuxième partie du cours. Dès que le calcul à effectuer
implique un enchaînement de commandes un peu
compliqué, il vaut mieux écrire ces dernières dans un fichier.
2018/2019
Par convention un fichier contenant des commandes Matlab

porte un nom avec le suffixe .m et s’appelle pour cette
raison un M-file ou encore script. On utilisera toujours l’éditeur
intégré au logiciel qui se lance à partir de la fenêtre de
commande en cliquant sur les icônes new Script ou open file
dans la barre de menu.
Par exemple, créons à l’aide de l’éditeur intégré de Matlab
dans le répertoire de travail choisi, déjà déclaré par SetPath, le
fichier test.m.
2018/2019
Supposons qu’il contient les instructions suivantes :
clear all
x=4;
y=2;
a=x +y
b=x*y
Ecrivons dans la fenêtre de commandes le nom du fichier
>> test
a=
6
b=
8
2018/2019
En tapant whos ensuite, on produit la sortie suivante:

2018/2019
Habituellement, on utilise les fichiers script afin de :

– Déclarer des variables ;
– Effectuer des opérations mathématiques ;
– Appeler des fonctions ;
– Tracer des figures ;
– Programmer des algorithmes.
2018/2019
 Fichiers function (M-file function)

L’idée de base d’une fonction est d’effectuer des opérations
sur une ou plusieurs entrées ou arguments pour obtenir un
résultat qui sera appelé sortie. Il est important de noter
que l’appel de la fonction se fait en précisant ses variables
entrées si ces dernières ne sont pas disponibles à l’invite
Matlab.
2018/2019
Par exemple, la fonction suivante admet une seule sortie a
qui constitue le résultat de l’addition des deux arguments
d’entrée x et y:
function a = ma_function(x,y)
a=x+y;
Lorsqu’on tape dans la fenêtre de commandes:
 >> a = ma_function(4,2)
 on obtient
 a=6
 Ensuite, on vérifie que
 >> whos
 Name Size Bytes Class
 a 1x1 8 double array
 Grand total is 1 element using 8 bytes
2018/2019
Modifions la fonction pour lui demander de calculer

aussi le produit de x et y de la façon suivate :
function [a,b] = ma_function2(x,y)
a=x+y;
b=x*y;
Dans ce cas, on vérifie que:
>> [a,b] = ma_function2(4,2)
a=
6
b=
8
2018/2019
>> whos
Grand total is 2 elements using 16 bytes.
On peut éviter l’affichage des sorties en utilisant le point-
virgule :
>> [a,b]=ma_function2(4,2);
2018/2019
Remarquons que nous pouvons appeler la fonction

ma_function2 pour calculer uniquement la
somme de x et y. On obtient alors
>> a = ma_function2(4,2)
a=
6
Habituellement, on utilise les fichiers function afin de :
– Programmer des opérations répétitives ;
– Limiter le nombre de variables dans l’invite Matlab ;
– Diviser le programme (problème) de manière claire.
2018/2019
 Définition d’une fonction par la commande «

inline »
Une fonction ne comportant qu’un petit nombre
d’instructions peut être définie directement
dans la fenêtre de commandes de la manière suivante :
>>angle=inline(’atan(y/x)’)
angle =
Inline function:
angle(x,y) = atan(y/x)
>>angle(5,4)
ans =
0.6747
2018/2019
Les arguments de la fonction angle sont normalement fournis à

l’appel dans l’ordre d’apparition dans la définition de la
fonction. On peut aussi spécifier les arguments d’appel
explicitement
>>f =inline(’sin(alpha*(x+y))’,’x’,’y’,’alpha’)
f=
Inline function:
f(x,y,alpha) =sin(alpha*(x+y))
>>f(0.2,0.3,pi)
ans =
1
2018/2019
Fonctions et représentation graphique sous Matlab
 Graphiques simples
Cette section vise une initiation aux

nombreuses facultés graphiques offertes par
Matlab. Dans toutes les représentations
graphiques, le logiciel se base sur des données
discrètes rangées dans des matrices ou des
vecteurs colonnes.
2018/2019
Par exemple,
pour représenter des courbes du type y = f(x) ou des

surfaces z = f(x, y), les données x,y,z doivent être des vecteurs
colonnes (x et y) ou des matrices (z) aux dimensions
compatibles.
L’instruction de dessin correspondante (par exemple plot(x,y)
pour tracer des courbes planes) est alors utilisée et
éventuellement complétée par des arguments optionnels
(couleur, type de trait, échelle sur les axes, etc...).
La visualisation du résultat s’effectue dans une fenêtre
graphique (avec possibilité de zoom, de rotation,
d’impression).
2018/2019
La fonction plot permet de tracer des courbes

en Matlab.
Les arguments de cette fonction sont les vecteurs
des variables indépendantes et dépendantes (en
alternance), comme dans l’exemple qui suit :
>> x=[0:0.01:2*pi];
>> plot(x,cos(x),x,sin(x))
On aurait aussi pu simplifier par :
>> plot(x’,[cos(x)’ sin(x)’])
2018/2019
Ces graphiques manquent cependant de

clarté. Il faut toujours nommer les axes (et
mettre les unités si possible), proposer une
légende, etc. Ceci est évidemment un peu
long à écrire à l’invite Matlab. Un script est
tout indiqué :
2018/2019
% graphique.m
clear all
close all % ferme les anciennes figures
x=[0:0.01:2*pi]; y1=cos(x); y2=sin(x);
figure(1)
plot(x,y1,’.’,x,y2,’+’) % cos(x) en points, sin(x) en
+
title(’sinus et cosinus’)
xlabel(’x’)
ylabel(’f(x)’)
legend(’cos(x)’,’sin(x)’,0) % le 0 place la légende
à côté des courbes
2018/2019
Pour en savoir plus, particulièrement sur les couleurs et

types de courbes, tapez help plot à l’invite Matlab.
•Fixez la limite des axes et la grille
grid on % ajoute une grille sur le graphique de son
choix
axis([xmin xmax ymin ymax]) % On fixe les limites sur
les axes des abscisses et des ordonnées
2018/2019
 Fonctions mathématiques simples
Les opérateurs algébriques (+, -, *, /, .*,./) ont été
définis précédemment pour les scalaires, vecteurs et
matrices.
On montrera ici (sans être exhaustif), les principales
fonctions mathématiques fournies dans Matlab et leur
utilisation. Pour les fonctions non présentées, l’utilisateur
peut toujours utiliser l’aide des fonctions avec la
fonction help qui prend pour argument le nom de la
fonction. Par exemple, la fonction cosinus :
>> help cos
COS Cosine of argument in radians.
COS(X) is the cosine of the elements of X.
2018/2019
 Fonctions mathématiques usuelles
Toutes les fonctions mathématiques de base sont déjà

programmées dans Matlab.
Toutes les fonctions courantes et moins courantes
existent. La plupart d’entre elles fonctionnent en
complexe. On retiendra que pour appliquer une
fonction à une valeur, il faut mettre cette dernière
entre parenthèses. Exemple :
>> sin(pi/12)
ans =
0.16589613269342
2018/2019
Voici une liste non exhaustive :

– fonctions trigonométriques et inverses : sin, cos, tan, asin,
acos, atan
– fonctions hyperboliques (on rajoute «h») : sinh, cosh, tanh,
asinh, acosh, atanh
– racine, logarithmes et exponentielles : sqrt, log, log10,
exp
– fonctions erreur : erf, erfc
– fonctions de Bessel et Hankel: besselj, bessely, besseli,
besselk, besselh. Il faut deux paramètres : l’ordre de la
fonction et l’argument lui-même.
La notion de fonction est plus générale dans Matlab, et
certaines fonctions peuvent avoir plusieurs entrées (comme
besselj par exemple) mais aussi plusieurs sorties.
2018/2019
 Fonctions matricielles
Toutes les fonctions matricielles de base sont

déjà programmées dans Matlab.
Voici quelques exemples :
Size, length, diag, det, norm, rank, trace, sum,
prod, mean, std, var, max, min, rand, null, inv,
pinv, sort, reshape, fliplr, flipud, tril, triu,…
2018/2019
 Fonctions avancées
Ce sont des fonctions qui interviennent en
analyse numérique telles que: lu, chol, qr, cond,
eig,fzero,…
2018/2019
Exemple de représentation graphique
2018/2019
 Remarquer l’emploi d’une double

apostrophe pour en faire figurer une dans
une chaîne de caractères délimitée
justement par deux apostrophes.
 On peut représenter plusieurs courbes sur
la même figure de plusieurs manières:
d’abord par un seul appel à la fonction plot
2018/2019
2018/2019
clear all
close all
x=linspace(0,pi,30);
y1=cos(x);
plot(x,y1,'o-r')
y2=sin(x);
hold on
plot(x,y2,'x-b')
y3=exp(-x);
hold on
plot(x,y3,'*-g')
2018/2019
Fin de séance
2018/2019

Séance 1-Calcul Scientifique

Transféré par

Droits d'auteur :

Formats disponibles

Séance 1-Calcul Scientifique

Transféré par

Informations du document

Description originale:

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Séance 1-Calcul Scientifique

Transféré par

Droits d'auteur :

Formats disponibles

Université Moulay Ismaïl

Faculté des Sciences et Techniques; Errachidia

Groupe facebook: Calcul scientifique s5

Pr. SARA TEIDJ

 Les concepts rattachés au calcul

 Résolution des équations non linéaires.

 Résolution des systèmes linéaires.

Ce chapitre est une initiation au logiciel Matlab et à la

La figure 1 présente l’écran Matlab de base: la fenêtre de commandes.

Dans la fenêtre de commande, l’utilisateur peut affecter des

Ici, il faut noter que lorsque l’utilisateur ne fixe pas de

La solution de x+y a donc été perdue. Il est donc

Notons au passage que le point-virgule permet de ne pas

La fonction clear permet d’effacer des variables. Par exemple :

La sortie de la fonction whos donne, entre autre, la classe de la

Il est possible de concaténer des mots à l’aide des

•Les formats d'affichage des réels

L’élément de base de Matlab est la matrice. C’est-à-dire qu’un scalaire

Le double point (:) est l’opérateur d’incrémentation dans

Les opérations usuelles d’addition, de soustraction et de

>> V=2*V1 % multiplication par un scalaire

>> V=V1./V2 % division élément par élément

Les matrices peuvent aussi être construites directement :

On peut aussi ”compter” les éléments. Matlab compte alors

Pour la multiplication et la division, les opérateurs usuels (* et

Dans certaines applications, il est parfois utile de connaître les

det(A) % calcule le déterminant de A

 Extraction d’une sous-matrice

 A(j:k) extrait les éléments j à k de A et les stocke dans un

>> polyval(p,1) : évalue le polynôme

Jusqu’à présent, l’utilisation que nous avons faite de Matlab

Par convention un fichier contenant des commandes Matlab

En tapant whos ensuite, on produit la sortie suivante:

Name Size Bytes Class

Habituellement, on utilise les fichiers script afin de :

 Fichiers function (M-file function)

Modifions la fonction pour lui demander de calculer

Remarquons que nous pouvons appeler la fonction

 Définition d’une fonction par la commande «

Les arguments de la fonction angle sont normalement fournis à

Cette section vise une initiation aux

pour représenter des courbes du type y = f(x) ou des

La fonction plot permet de tracer des courbes

Ces graphiques manquent cependant de

Pour en savoir plus, particulièrement sur les couleurs et

 Fonctions mathématiques usuelles

Toutes les fonctions mathématiques de base sont déjà

Voici une liste non exhaustive :

Toutes les fonctions matricielles de base sont

 Remarquer l’emploi d’une double

Vous aimerez peut-être aussi

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.