Scribd
Importer
24 vues4 pages

TP-01 Solution

Le document présente des exemples de calculs élémentaires et matriciels, y compris des opérations sur des matrices, le calcul de rang, de déterminant et d'inverse. Il aborde également la manipulation de polynômes, la fonction de transfert, et la création de graphiques sinusoïdaux et en 3D. Chaque section contient des commandes spécifiques pour effectuer ces calculs et visualisations.

Transféré par

sebahmechri3
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
24 vues4 pages

TP-01 Solution

Le document présente des exemples de calculs élémentaires et matriciels, y compris des opérations sur des matrices, le calcul de rang, de déterminant et d'inverse. Il aborde également la manipulation de polynômes, la fonction de transfert, et la création de graphiques sinusoïdaux et en 3D. Chaque section contient des commandes spécifiques pour effectuer ces calculs et visualisations.

Transféré par

sebahmechri3
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

1 Calculs élémentaires

Tapez la commande :
1 a = 5 + 9;
2 b = a + 2;
3 disp(a);
4 disp(b);

2 Calculs matriciels
Les matrices données sont :
   
1 3 6 0 3 16
A = 5 4 8 , B = 5 10 8
7 5 2 9 7 2

• Calcul du rang de A :

1 A = [1 3 6; 5 4 8; 7 5 2];
2 rank_A = rank(A);
3 disp(rank_A);

• Somme de A et B :

1 B = [0 3 16; 5 10 8; 9 7 2];
2 C = A + B;
3 disp(C);

• Produit matriciel de A et B :

1 D = A * B;
2 disp(D);

• Multiplier chaque élément de A et B :

1 E = A .* B;
2 disp(E);

• Transposée de A :

1
1 A_T = A’;
2 disp(A_T);

• Déterminant et inverse de A :

1 det_A = det(A);
2 inv_A = inv(A);
3 disp(det_A);
4 disp(inv_A);

3 Polynôme
Le polynôme donné est :
p(x) = 3x2 + 2x + 5
Les racines du polynôme sont calculées par :
1 coeffs = [3 2 5];
2 r = roots(coeffs);
3 disp(r);

Tracer le polynôme sur l’intervalle [−10, 10] :


1 x = -10:0.1:10;
2 p = 3*x.^2 + 2*x + 5;
3 plot(x, p);

4 Fonction de transfert
La fonction de transfert donnée est :
s+2
G(s) =
3s3 + 2s2 + 4s − 1
Pour calculer les pôles et zéros :
1 num = [1 2]; % Numrateur : S + 2
2 den = [3 2 4 -1]; % Dnominateur : 3s^3 + 2s^2 + 4s - 1
3 zeros_G = roots(num);
4 poles_G = roots(den);
5 disp(zeros_G);
6 disp(poles_G);

Pour tracer la réponse fréquentielle de G(s) :


1 sys = tf(num, den);
2 bode(sys);

2
5 Tracer trois courbes sinusoı̈dales
Tracer trois courbes sinusoı̈dales avec des styles différents :
1 t = 0:pi/50:10*pi;
2 y1 = sin(t);
3 y2 = cos(t);
4 y3 = sin(t + pi/4);
5 plot(t, y1, ’g-’, t, y2, ’b--o’, t, y3, ’c*’);

6 Tracer en 3D
Tracer plusieurs ensembles de coordonnées sur un même ensemble d’axes :
1 r = 2; xc = 4; yc = 3;
2 theta = linspace(0, 2*pi);
3 x = r*cos(theta) + xc;
4 y = r*sin(theta) + yc;
5 plot(x, y);
6 axis equal;
7 xlabel(’x’);
8 ylabel(’y’);

premier code

deuxième code

3
troisième code

Vous aimerez peut-être aussi