MINISTRE DES ENSEIGNEMENTS SUPERIEURS DE LA RECHERCHE ET DE L’INNOVATION

UNIVERSITE DAN DICKO DANKOULODO DE MARADI(UDDM)

FACULTE DES SCIENCES ET TECHNIQUES(FST)

DEPARTEMENT DE PHYSIQUE

RAPPORT DES TRAVAUX PRATIQUE SUR MATLAB

Filière : L3 ENERGIE RENOUVELABLE PHE

Rédigé par : Sani Mamadou Safia

Instructeur : Dr HAMIDINE Mahamadou

Année académique : 2025-2026

 TP N°1

1. Fenêtre de commande et calculs de base

Tache : Executer dans Matlab

3+2*4^2

ans =

35

>> ((3+2)*4)^2

ans =

400

>> 1/12-6^2+2/3

ans =

-35.2500

>> 5*(3/4)+9/4

ans =

>> 4^3*[3/4+9/(2)^3]

ans =

120

>> x=2*pi/3

x=

2.0944

2. Gestion des variables

>> y=sin(x)

y=

0.8660

>> A=4/5*7^3

A=

274.4000

>> B=A^2
B=

7.5295e+04

>> T='bonjour'

T=

bonjour

a. Exécuter les commandes who et whos successivement

>> who

Your variables are:

A B T ans x y

>> whos

Name Size Bytes Class Attributes

A 1x1 8 double

B 1x1 8 double

T 1x7 14 char

ans 1x1 8 double

x 1x1 8 double

y 1x1 8 double

b. Exécuter la commande clear pour effacer x,y et T et apres éxecuter whos

>> clear x y T

>> whos

c. Exécuter clear all et ensuite who

>> clear all

>> who

3. Commande d’affichage
Tache : Exécuter dans Matlab

>> format short

>> x=3+11/16+2^1.2

x=

5.9849

>> format long

>> x=3+11/16+2^1.2
x=

5.984896709994070

>> format short e

>> x=3+11/16+2^1.2

x=

5.9849e+00

>> format long e

>> x=3+11/16+2^1.2

x=

5.984896709994070e+00

>> format rat

>> x=3+11/16+2^1.2

x=

1981/331

>> format bank

>> x=3+11/16+2^1.2

x=

5.98

4. Définitions mathématiques de base

Tache :
i. Trouver la racine carr (sqrt), exp(x), le log(x), log10(x)

>> sqrt(x)

ans =

2.45

>> exp(x)

ans =

397.38

>> le log(x)

Error using <=

Not enough input arguments.

>> log(x)

ans =
 1.79

>> log10(x)

ans =

0.78

ii. Trouver le volume d’une sphère qui a un rayon de 2m

>> r=2

r=

2.00

>> v=4*pi/3*r^3

v=

33.51

5. Nombres complexes
Tache :

>> a=5+i*2;

>> b=1-i*3;

>> format short

>> a+b

ans =

6.0000 - 1.0000i

>> a-b

ans =

4.0000 + 5.0000i

>> (conj(a))

ans =

5.0000 - 2.0000i

>> (conj(b))

ans =

1.0000 + 3.0000i

>> (abs(a))
ans =

5.3852

>> (abs(b))

ans =

3.1623

>> (angle(a))

ans =

0.3805

>> (angle(b))

ans =

-1.2490

6. Création d’un fichier dans Matlab

> % fichier exemple1.m pour calculer l'exponentielle d'un groupe de nbombre

>> clear all %eface toutes les variables du workspace

>> x=[1:2:3:4];

>> y=exp(x)

y=

2.7183 7.3891 20.0855 54.5982

>> %calcul des temperatures en degré celsius

>> clear all %eface toutes les variables du workspace

>> temps=[32,50,65,70,85];

>> T=5/9*(temps-32)

T=

0 10.0000 18.3333 21.1111 29.4444

 TP N°2

1. Création

Tâche 1:

i) Créez un vecteur ligne A et un vecteur colonne B avec trois éléments v1=1.1 ; v2=2.2 et

v3=3.3 :

>> v1=1.1;

>> v2=2.2;

>> v3=3.3;

>> A=[v1 v2 v3]

A=

1.1000 2.2000 3.3000

>> B=[v1;v2;v3]

B=

1.1000

2.2000

3.3000

ii) Créez une matrice C avec deux lignes et trois colonnes. Choisissez des valeurs

quelconques :

>> C=[23 35 44;15 20 67]

C=

23 35 44

15 20 67

iii) Créez la matrice suivante dans Matlab D=[1 1 1 ; 1 2].Que dit Matlab ?

>> D=[1 1 1;1 2]

Error using vertcat

Dimensions of matrices being concatenated are not consistent.

Tâche 2 : Créer une matrice unitaire de 4x4 et 10x10 ; Créer une matrice de 5x5 zéros et aussi une
matrice de 5x6 zéros. Créer une matrice de 8x8 un et une matrice de 10x12 un :

>> eye(4,4)

ans =

1 0 0 0
 0 1 0 0

0 0 1 0

0 0 0 1

>> eye(10,10)

ans =

1 0 0 0 0 0 0 0 0 0

0 1 0 0 0 0 0 0 0 0

0 0 1 0 0 0 0 0 0 0

0 0 0 1 0 0 0 0 0 0

0 0 0 0 1 0 0 0 0 0

0 0 0 0 0 1 0 0 0 0

0 0 0 0 0 0 1 0 0 0

0 0 0 0 0 0 0 1 0 0

0 0 0 0 0 0 0 0 1 0

0 0 0 0 0 0 0 0 0 1

>> zeros(5,5)

ans =

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

>> zeros(5,6)

ans =

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

>> ones(8,8)

ans =
 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

>> ones(10,12)

ans =

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1 1 1 1 1

Tâche 3:

i) créer une liste p des nombres pairs de 0 à 10 :

>> p=0:2:10

p=

0 2 4 6 8 10

ii) créer une liste x de nombres entre 0 et 1 avec un incrément d’un dixième :

>> x=0:0.1:1

x=

0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000

Considérant la liste x ci- dessus, trouvez :

 a. y=e^x

>> y=exp(x)

y=

1.0000 1.1052 1.2214 1.3499 1.4918 1.6487 1.8221 2.0138 2.2255 2.4596 2.7183

b. y=x^2

>> Y=x.^2

Y=

0 0.0100 0.0400 0.0900 0.1600 0.2500 0.3600 0.4900 0.6400 0.8100 1.0000

iii) créer une liste u de nombre de 100 à 80 en décroissant de 5 :

>> u=100:-5:80

u=

100 95 90 85 80

iv) créer une liste z de nombre de 10 à -10 avec un espacement de 2 :

>> z=(10:-2:-10)

z=

10 8 6 4 2 0 -2 -4 -6 -8 -10

Tâche 4:

i) Choisit deux nombres quelconques et un espacement de votre choix et exécuter les commandes
linspace et logspace :

>> x1=15;

>> x2=8;

>> n=4;

>> linspace(x1, x2, n)

ans =

15.0000 12.6667 10.3333 8.0000

>> logspace(x1, x2, n)

ans =

1.0e+15 *

1.0000 0.0046 0.0000 0.0000

ii) Soient les matrices données : A=[1,3,5], B=[2,4,1], C=[1,1;1,2]. Trouvez [A,B],[A;B] et [A,C]

>> A=[1 3 5];

>> B=[2 4 1];

>> C=[1 1;1 2];

>> [A B]

ans =

1 3 5 2 4 1

>> [A;B]

ans =

1 3 5

2 4 1

>> [A C]

Error using horzcat

Dimensions of matrices being concatenated are not consistent.

2. Accès et modifications

Tâche 1:

i) Soit V = [12; 17; –2; 0; 4; 4; 11; 19; 27];

 Trouver l’élément de la position 2 et la position 8 :

>> V=[12;17;-2;0;4;4;11;19;27];

>> V(2)

ans =

17

>> V(8)

ans =

19

 Trouver aussi la liste de tous les éléments :

>> V(:,1)

ans =

12

17

-2

0
 4

11

19

27

 Trouver les éléments entre les positions 4 et 6 :

>> V([4:6])

ans =

ii) Soit la matrice A

a)trouver l’élément qui se trouve à l’intersection de la deuxième ligne et la troisième colonne :

>> A=[1 2 3;4 5 6;7 8 9 ];

>> A(2,3)

ans =

b) Repérer les éléments de la colonne 2 :

>> A(:,2)

ans =

c) repérer les éléments de la colonne 2 et 3 :

>> A(:,2:3)

ans =

2 3

5 6

8 9

d) changer l’élément qui se trouve à l’intersection de la ligne 1 et la colonne 1 par -8 :

>> A(1,1)=-8
A=

-8 2 3

4 5 6

7 8 9

e) vider la ligne 2 de A et aussi vider la colonne 2 de A :

>> A(2,:)

ans =

4 5 6

>> A(:,2)

ans =

f) créer une nouvelle matrice 4x4 former des éléments de la première ligne de A :

>> A=[1 2 3;4 5 6;7 8 9];

>> M=A([1,1,1,1],:)

M=

1 2 3

1 2 3

1 2 3

1 2 3

ii) Aussi former une matrice 3x3 dont la première et la dernière ligne sont formées par la première
ligne de A et la deuxième ligne par la ligne 2 de A.

>> T=A([1,1,2],:)

T=

1 2 3

1 2 3

4 5 6

Tâche 2: trouver la taille, le minimum et le maximum des vecteurs suivants :

>> A=[2;3;4;5;6];

>> B=[8 4 4 1 7 11 2 0];

>> length(A)

ans =

>> length(B)

ans =

>> min(A)

ans =

>> min(B)

ans =

>> max(A)

ans =

>> max(B)

ans =

11

3. Opérations avec les matrices

A. Le module d’un vecteur

Tâche : Trouver le module des vecteurs colonnes suivant : J avec pour éléments 0,3 et 4 ; U avec
éléments i, 1+2i et 4 :

>> J=[0.3 4];

>> J.*J

ans =

0.0900 16.0000

>> sum(J.*J)

ans =

16.0900

>> sqrt(sum(J.*J))

ans =
 4.0112

>> U=[i 1+2*i 4];

>> sqrt(sum(U.*conj(U)))

ans =

4.6904

[Link] cammande abs

Cette commande est utilisée pour calculer la valeur absolue d’un vecteur.

Tâche : choisissez un vecteur quelconque avec des éléments négatifs et trouver sa valeur absolue.

>> X=[-2 -5 -13;-4 -19 -65;-9 -34 -56]

X=

-2 -5 -13

-4 -19 -65

-9 -34 -56

>> abs(X)

ans =

2 5 13

4 19 65

9 34 56

[Link] scalaire et produit vectorielle

Tâche :

i) soit U ; Trouver le produit scalaire U.U :

>> U=[-i ;1+i ;4+4i];

>> dot(U,U)

ans =

35

ii) soient A ;B ; calculer le produit vectorielle A∧B :

>> A=[1 2 3];

>> B=[2 3 4];

>> cross(A,B)

ans =

-1 2 -1
4. Operations de bases avec les matrices

Tâche :

i) effectuer A+B et A-B dans Matlab :

>> A=[5 1;0 9];

>> B=[2 -2;1 1];

>> A+B

ans =

7 -1

1 10

>> A-B

ans =

3 3

-1 8

ii) calculer la transposer dans Matlab de A :

>> A=[-1 2 4;0 1 6;2 7 1];

>> A'

ans =

-1 0 2

2 1 7

4 6 1

iii) Multiplier A et B :

>> A=[2 1;1 2];

>> B=[3 4;5 6];

>> A*B

ans =

11 14

13 16

Multiplier C et D

>> C=[1 4;8 0;-1 3];

>> D=[-1 7 4; 2 1 -2];

>> C*D
ans =

7 11 -4

-8 56 32

7 -4 -10

iv) additionné A et b :

>> A=[1 2 3 4];

>> b=2;

>> A+b

ans =

3 4 5 6

v) Effectuer la division de A et B ;en faisant la division à droite et la division à gauche :

>> A=[1 2 3 4];

>> B=[2 2 3 1];

>> A./B

ans =

0.5000 1.0000 1.0000 4.0000

>> B./A

ans =

2.0000 1.0000 1.0000 0.2500

vi) soit B , trouver la matrix carré de B en effectuant la multiplication par élément :

>> B=[2 4;-1 6];

>> B.^2

ans =

4 16

1 36

5. Déterminant et résolution de systèmes d’équations linéaires

A. Déterminant et résolution de systèmes d’équations linéaires

Tache :
i) Trouver le déterminant de chacune des matrices suivantes :

>> A=[1 3;4 5];

>> det(A)

ans =

-7

>> B=[3 -1 2 4;0 2 1 8;-9 17 11 3;1 2 3 -3];

>> det(B)

ans =

-533

ii) Trouver les solutions du système d’équation suivant :

>> E=[5 2 -9;-9 -2 2;6 7 3]

E=

5 2 -9

-9 -2 2

6 7 3

>> F=[44;11;5]

F=

44

11

>> det(E)

ans =

437

>> F\E

ans =
 0.0725 0.0485 -0.1724

B. la commande rank

Tache : Trouver le rang d’une matrice B

>> B=[2 4;-1 6]

B=

2 4

-1 6

>> rank(B)

ans =

C. L’inverse d’une matrice

Tache : Trouver le déterminant de chacune des matrices et vérifier si réellement c’est

l’inverse de la matrice

>> A=[2 3;4 5]

A=

2 3

4 5

>> S=[1 0 -1 2;4 -2 -3 1;0 2 -1 1;0 0 9 8]

S=

1 0 -1 2

4 -2 -3 1

0 2 -1 1

0 0 9 8

>> det(A)
ans =

-2

>> det(S)

ans =

-108

>> inv(A)

ans =

-2.5000 1.5000

2.0000 -1.0000

>> inv(S)

ans =

-0.9259 0.4815 0.4815 0.1111

-0.6296 0.1574 0.6574 0.0556

-0.5926 0.1481 0.1481 0.1111

0.6667 -0.1667 -0.1667 0

>> A*A^-1

ans =

1 0

0 1

>> A^-1*A

ans =

1 0

0 1

>> S*S^-1

ans =

1.0000 0 0 0
 0.0000 1.0000 -0.0000 0

0.0000 -0.0000 1.0000 0

0 0 0 1.0000

>> S^-1*S

ans =

1.0000 0 0 0

0 1.0000 0 0.0000

0 0 1.0000 0

0 0 0.0000 1.0000

D. Les commandes rref et magic

Tache :

i) réduire l’échelon de A :

>> A=[1 2;4 7]

A=

1 2

4 7

>> rref(A)

ans =

1 0

0 1

ii) Trouver une matrice de 5x5 tel que la somme des éléments de la ligne égale à la somme
des éléments de colonne :
>> magic(5)

ans =

17 24 1 8 15

23 5 7 14 16

4 6 13 20 22

10 12 19 21 3

11 18 25 2 9
 TP N3

1. Courbes dans le plan : Graphisme 2D

Tache 1 : Tracer le graph de la fonction y=cos(x) dans l’intervalle 0≤ x ≤4π avec un pas de 0.01.

>> x=0:0.01:4*pi;

>> y=cos(x);

>> plot(x, y)

Tache 2: Exécuter la tache ci-dessus en augmentant la commande xlabel, ylabel et title :

>> xlabel ('x')

>> ylabel ('cosx')

>> title ('Graph de y=cos(x)')

 Tache 3:

i) En utilisant la commande fplot tracer le graph de exp(2t)sint, pour 0≤t≤4 avec les titres t

pour l’axe des x et f(t) pour l’axe des y :

>> fplot (‘exp (2*t).*sin(t)', [0 4]);

>> xlabel ('x');ylabel ('y');title ('Graph de y=exp (2*t).*sin(t)')

 ii) Utiliser la commande plot avec un intervalle de 0.01 :

>> t=0:0.01:4;

>> z=exp (2*t).*sin(t);

>> plot(t, z)

>> xlabel ('x');ylabel ('y');title ('Graph de y=exp (2*t).*sin(t)')

Tache 4: Soit à tracer la fonction y=tanhx sur l’intervalle -6≤x≤6 avec un espacement de 0.01 :

i) Avec la grille sur le graph :

>> x=-6:0.01:6;

>> y=tanh(x);

>> plot(x, y), grid on

ii) Les axes sont carrées :

>> plot(x, y), axis square

iii) Avec un même facteur d’échelle pour les deux axes :

>> plot(x, y), axis equal

 Tracé plusieurs fonctions sur un même graphique

Tache 1: Tracer les deux fonctions f et g suivantes dans l’intervalle 0≤t≤5 :

i) spécifiez la forme de la courbe de la seconde fonction en laissant la pour la première par défaut :

>> t=0:0.01:5;

>> f=exp (-t);

>> g=exp (-2*t);

>> plot(t, f,t,g,'--')

ii) Donner une autre spécificité des formes des courbes des deux fonctions différentes de celle de i) :

>>plot(t,f,'-.',t,g,'-')
 Tache 2: Tracer deux fonctions de l’énergie potentielle qui sont définies en termes des fonctions
trigonométriques hyperboliques sinh(x) et cosh(x) pour 0≤x≤2. Montrez sur votre graph la légende et
aussi les titres pour l’axe des abscisses et l’axe des ordonnées :

>> x=0:0.01:2; h=sinh(x); k=cosh(x);

>> plot(x,h,x,k);xlabel('x');ylabel('y');title('graph de deux fonctions');legend('sinh(x)','cosh(x)')

>> legend('show')

Mettre des couleurs sur les tracés dans le graphique

Tache 1: Tracez les courbes de sinus et cosinus hyperbolique définies sur -5≤x≤5 avec la courbe de la
première en rouge et celle de la seconde en bleue :

>> x=-5:0.01:5; R=sinh(x); S=cosh(x);

>> plot(x, R,'red', x, S,'b')

 Tache 2:

1. Tracer la courbe de la fonction y :

i) Avec 0≤x≤5 et -1≤y≤1

>> x=0:0.01:5;y=-1:0.01:1;y=exp(-3/2*x).*sin(5*x+3);

>> plot(x,y),axis([0 5 -1 1])

 ii)avec -0.7 ≤y ≤0.3

>> axis([0 5 -0.7 0.3])

[Link] la courbe de la fonction y ,sur :

i) 0≤x≤5

>>x=0 :0.01 :5 ;y=(sin(5*x).^2 ;

>>Plot(x,y)
 ii) Sur 0≤x≤1

>> x=0:0.01:1;y=(sin(5*x)).^2;

>> plot(x,y)
Tache 3: Tracer cote à cote les fonctions suivantes y et z avec 0≤x≤5 et -1≤y≤1 :

>> x=0:0.01:5; y = exp(-1.2*x).*sin(20*x);subplot(1,2,1)

>> plot(x,y)

>> z=exp(-2*x).*sin(20*x);subplot(1,2,2)

>> plot(x,z)

>> plot(x,y),axis([0 5 -1 1])

>> plot(x,z),axis([0 5 -1 1])

 Superposition de deux courbes

Tache : Superposer les courbes des fonctions cos(x) et sin(x) avec 0≤ x ≤2π et -1≤y≤1. Utiliser la
commande linspace pour définir l’espacement des x :

>> x=0:0.01:2*pi;u=cos(x);v=sin(x);

>> plot(x,u),hold on,plot(x,v),hold off ;

>> axis([0 2*pi -1 1]);

>> linspace(0.2*pi,0.01) ;
 2. Courbes polaires et logarithmiques

Tache 1:

i) Tracer le contour de la spirale d’Archimède donnée par r . Ou a est une constante avec a=2 et 0≤
θ≤2π . Prendre un espacement ou incrément de π/90 pour θ . Donner aussi le titre Spirale
d’Archimède :

>> a=2;theta=0:pi/90:2*pi;r=a.*theta;

>> polar(theta,r)

>> title('Spirale d’Archimède')

 Spirale d’Archimède
90
15
120 60

10

150 30
5

180 0

210 330

240 300
270

ii) Tracer aussi le contour polaire de r, avec 0≤ θ≤6π, la courbe doit être en point tillées interrompu et
couleur rouge :

>> theta=[0:pi/90:6*pi];

>> r=1+2*cos(theta);

>> polar(theta,r,'r :')

Tache 2: Tracer la courbe de y, ou 0≤ x≤20 avec un incrément de 0.1 sur l’échelle logarithmique après
avoir exécuté la commande plot :

>>x=0:0.1:20;y=exp(-10*x.^2);plot(x,y)

>> loglog(x,y)
 2. Graphisme en 3D

Tache 1: Tracer les contours de z avec −5 ≤ x ≤ 5, −3 ≤ y ≤ 3

>> x=-5:0.01:5;y=-3:0.01:3;

>> [x,y]=meshgrid(-5:0.01:5,-3:0.01:3);

>> z=x.^2+y.^2;

>> contour(x,y,z)
 3

-1

-2

-3
-5 -4 -3 -2 -1 0 1 2 3 4 5

Tache 2: Tracer les contours de la fonction z avec les titres de chaque contour en utilisant les mêmes
intervalles que la tache précédente :

i) Pour la fonction z tracer le contour en 3D avec 10 niveaux :

>> x=-5:0.01:5;y=-3:0.01:3;

>> [x,y]=meshgrid(-5:0.01:5,-3:0.01:3);z=cos(x).*sin(y) ;

>> contour(z,10);
set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2);
ii) Tracer simplement les isolignes de z , avec −2 ≤ x, y ≤ 2 en donnant les titres de l’axe des x par x et
l’axe y par y. Utiliser ensuite la commande contour (x,y,z,n) avec n=30 et comparer les résultats :

>> x=-2:0.01:2;

>> [x,y]=meshgrid(-2:0.01:2,-2:0.01:2);

>> z=y.*exp(-x.^2-y.^2);

>> contour(x,y,z)

>> xlabel('x');ylabel('y');
 >> contour(x,y,z,30)

iii) Pour embellir l’apparence utiliser les commandes ci-après :

>> surface(x,y,z,'EdgeColor',[.8 .8 .8],'FaceColor','none')

>> grid off

>> view(-15,20)
 Tache 4: Donner le graphique en 3D de :

i) z avec −2𝜋 ≤ x, y ≤ 2𝜋

x=-2*pi:0.01:2*pi;

[x,y]=meshgrid(-2*pi:0.01:2*pi,-2*pi:0.01:2*pi);

>> z=cos(x).*sin(x);

>> mesh(x,y,z);xlabel('x');ylabel('y');zlabel('z')
 ii) z avec −2 ≤ x, y ≤ 2

>> x=-2:0.01:2;

>> [x,y]=meshgrid(-2:0.01:2,-2:0.01:2);

>> z=y.*exp(-x.^2-y.^2);

>> mesh(x,y,z);xlabel('x');ylabel('y');zlabel('z');
Tache 5: Exécuter les commandes surf et surfc pour les fonctions de la tache précédente.

Pour tracer une surface éclairée on utilise la commande surfl :

>>surf(x,y,z);xlabel('x');ylabel('y');zlabel('z');
>>surfc(x,y,z);xlabel('x');ylabel('y');zlabel('z');

>> surfl(x,y,z);xlabel('x');ylabel('y');zlabel('z');
 TP N4

Résolutions des Equations algébriques et calcul de limite, dérivé, intégral

Equation du premier et second degré

Equation du premier degré :

Tache 1 : Trouver la solution de l’équation 2x+5=0 :

>> s=solve('2*x+5')

s=

-5/2

Tache 2 : Trouver a dans l’équation ax+5=0 :

>> solve('a*x+5','a')

ans =

-5/x

Equation du second degré :

Tache : Résoudre l’équation de la seconde dégrée : x2 − 6x − 12 = 0 :

>> solve('x^2-6*x-12')

ans =

3 - 21^(1/2)

21^(1/2) + 3

Tracer le graph des équations symboliques :

Tache :

A. Obtenir le graph de :

i) x2 − 6x − 12 = 0 :

>> ezplot('x^2-6*x-12')
 x 2 -6 x-12
70

60

50

40

30

20

10

-10

-20

-6 -4 -2 0 2 4 6
x

ii) x2 − 6x − 12 = 0 avec −2 < x < 8 :

>> ezplot('x^2-6*x-12',[-2 8])

x 2 -6 x-12
5

-5

-10

-15

-20

-2 -1 0 1 2 3 4 5 6 7 8
x
 iii) x+3=0 :

>> ezplot('x+3')

x+3
10

-2

-4
-6 -4 -2 0 2 4 6
x

iv) x + 3 = 0 avec −4 < x < 4, −2 < y < 2 :

>> ezplot('x+3',[-4 4])

>> axis([-4 4 -2 2])

 x+3
2

1.5

0.5

-0.5

-1

-1.5

-2
-4 -3 -2 -1 0 1 2 3 4
x

[Link] les solutions de x2 + x − √2 = 0 et tracer le graph de la fonction :

>> solve('x^2+x-sqrt(2)')

ans =

- (4*2^(1/2) + 1)^(1/2)/2 - 1/2

(4*2^(1/2) + 1)^(1/2)/2 - 1/2

>> ezplot('x^2+x-sqrt(2)')
 x 2 +x-sqrt(2)

45

40

35

30

25

20

15

10

-5
-6 -4 -2 0 2 4 6
x

Equation d’ordre supérieur à 2 :

Tache :

i) Résoudre (x + 1)(x+1) (x − 2) = 0 :

>> s=solve('(x+1)^2*(x-2)');double(s(1))

ans =

-1

>> s=solve('(x+1)^2*(x-2)');double(s(2))

ans =

-1

>> s=solve('(x+1)^2*(x-2)');double(s(3))

ans =

ii) Trouver les racines de x4 − 5x3 + 4x2 − 5x + 6 = 0 et tracer le

Graphe de la fonction sur −10 < x < 10 et rendre les racines en des valeurs simple :

>> s=('x^4-5*x^3+4*x^2-5*x+6');double(s(1))
 ans =

120

>> s=solve('x^4-5*x^3+4*x^2-5*x+6');double(s(2))

ans =

4.2588

>> s=solve('x^4-5*x^3+4*x^2-5*x+6');double(s(3))

ans =

-0.1876 - 1.1076i

>> s=solve('x^4-5*x^3+4*x^2-5*x+6');double(s(4))

ans =

-0.1876 + 1.1076i

>> ezplot('x^4-5*x^3+4*x^2-5*x+6',[-10 10])

x 4 -5 x 3 +4 x 2 -5 x+6

10000

9000

8000

7000

6000

5000

4000

3000

2000

1000

-10 -8 -6 -4 -2 0 2 4 6 8 10
x

iii) Trouver les racines de x3 + 3x2 − 2x – 6 et tracer le graph de la fonction sur −8 < x < 8, −8 < y < 8 :

>> s=solve('x^3+3*x^2-2*x-6')

s=
 -3

2^(1/2)

-2^(1/2)

>> ezplot('x^3+3*x^2-2*x-6',[-8 8]);axis([-8 8 -8 8])

x 3 +3 x 2 -2 x-6
8

-2

-4

-6

-8
-8 -6 -4 -2 0 2 4 6 8
x

Tracer le graph avec la grille :

>> grid on
 x 3 +3 x 2 -2 x-6
8

-2

-4

-6

-8
-8 -6 -4 -2 0 2 4 6 8
x

Système d’équation :

Tache : Résoudre le système d’équation suivant :

>> s=solve('w+x+4*y+3*z=5','2*w+3*x+y-2*z=1','w+2*x-5*y+4*z=3','w-3*z=9');

>> w=s.w

w=

1404/127

>> x=s.x

x=

-818/127

>> y=s.y

y=

-53/127

>> z=s.z

z=

87/127

Développement et réduction des Equations :

Tache 1: Développer pour obtenir des équations :

i) (x + 2) (x − 3) :

>> syms x

>> expand((x+2)*(x-3))

ans =

x^2 - x - 6

ii)cos(x+y)

>> syms x

>> syms y

>> expand(cos(x+y))

ans =

cos(x)*cos(y) - sin(x)*sin(y)

iii) sin(x-y) :

>> syms x

>> syms y

>> expand(sin(x-y))

ans =

cos(y)*sin(x) - cos(x)*sin(y)

iv) (y-2)(y+8) :

>> syms y

>> expand((y-2)*(y+8))

ans =

y^2 + 6*y – 16

Tache 2: Développer :

i) x(x2 − 2) :

>> syms x

>> collect(x*(x^2-2))

ans =

x^3 - 2*x

ii) (t+3)sint :

>> syms t
 >> collect((t+3)*sin(t))

ans =

sin(t)*t + 3*sin(t)

Tache 3: i) Trouver les facteurs de x2-y2 et x3+y3 :

>> f=x^2-y^2;g=x^3+y^3;

>> sym([f,g])

ans =

[x^2 - y^2, x^3 + y^3]

ii) Réduire en expression simplifiée :

a)

>> q=exp(2*log(3*x));

>> simplify(q)

ans =

9*x^2

b)

>> z=(cos(x))^2-(sin(x))^2;

>> simplify(z)

ans =

cos(2*x)

c)

>> w=(cos(x))^2+(sin(x))^2

w=

cos(x)^2 + sin(x)^2

>> simplify(w)

ans =

Equations avec fonctions exponentielles et Log :

Tache : Résoudre :

i)

>> h=solve('log10(x)-log10(x-3)-1')

h=
 10/3

ii) k et l :

>> k=solve('y-3^(2*x)')

k=

log(y)/(2*log(3))

>> l=solve('y-5^x-1')

l=

log(y - 1)/log(5)

iii) exp(x)+x et tracer la courbe de la fonction avec la grille :

>> m=solve('exp(x)+x')

m=

-lambertw(0, 1)

>> ezplot('exp(x)+x')

>> grid on

exp(x)+x

200

150

100

50

-6 -4 -2 0 2 4 6
x

Limite Dérivation et Intégration :

Tache 1:Trouver la limite quand x→3 des fonctions f et g :

>> syms x
 >> f(x)=(2*x+1)/(x-2);g(x)=x^2+1;

>> limit(f,3)

ans(x) =

>> limit(g,3)

ans(x) =

10

Tache 2: a)Montrer que cette limite n’existe pas :

>> limit((x-3)/abs(x-3),3)

ans =

NaN

b)Trouver la limite à gauche et à droite :

>> limit(((x-3)/abs(x-3)),x,3,'left')

ans =

-1

>> limit(((x-3)/abs(x-3)),x,3,'right')

ans =

Tache 3: i) Trouver les dérivées de f= x2 ; g=sint10t :

>> f=x^2;g=sin(t)*10*t ;

>> diff(f)

ans =

2*x

>> diff(g)

ans =

10*sin(t) + 10*t*cos(t)

ii) Trouver le minimum et le maximum de la fonction f (x) dans l’intervalle [0, 2] :

>> f(x)=x^3-3*x^2+3*x;

>> subs(f,0)

ans =

0
 >> subs(f,2)

ans =

Tache 4: calculer :

i)

>> int(x^5*cos(9*x))

ans =

(40*cos(9*x))/177147 + (40*x*sin(9*x))/19683 - (20*x^2*cos(9*x))/2187 + (5*x^4*cos(9*x))/81 -

(20*x^3*sin(9*x))/729 + (x^5*sin(9*x))/9

ii)

>> int(3*y*sec(x))

ans =

3*y*log((sin(x) + 1)/cos(x))

iii)

>> int(exp(-x^2)*sin(x),0,inf)

ans =

-(pi^(1/2)*erf(i/2)*i)/(2*exp(1/4))