Dans le cas des matrices stockées sous forme sparse, des algorithmes particuliers sont mis en oeuvre.Chacune des méthodes précédentes peut être utilisée de manière spécifique grâce aux commandeschol, lu, qr. Il est également possible d'utiliser des méthodes itératives. Les commandes cgs, bicg,bicgstab mettent par exemple en oeuvre des méthodes de type gradient conjugué.>> A=[1 2 ;3 4]; b=[1 1]';>> x= A\bx =-11>> A*xans =11>>3.6 Les polynômesSous Matlab le polynôme de degré n, p(x) = a n x n + a n-1 x n-1 + . . . + a 1 x + a 0 est défini par un vecteur pde dimension n+1 contenant les coefficients {a i } {i=0,...,n} rangés dans l'ordre décroissant des indices.C'est-à-dire que l'on a p(1) = a n , . . . , p(n+1) = a 0 . La commande polyval permet d'évaluer lepolynôme p (la fonction polynomiale) en des points donnés. La syntaxe est polyval(p,x) où x est unevaleur numérique ou un vecteur. Dans le second cas on obtient un vecteur contenant les valeurs de lafonction polynomiale aux différents points spécifiés dans le vecteur x. Utilisée avec la commandefplot, la commande polyval permet de tracer le graphe de la fonction polynomiale sur un intervalle[x min , x max ] donné. La syntaxe de l'instruction est : fplot('polyval([ a_n, ..., a_0] , x)' , [x_min ,x_max]). Voici par exemple comment définir le polynôme p(x) = x 2 -1. Le graphe de la fonctionpolynomiale est présenté à la figure 3.>> p = [ 1, 0, -1];>> polyval(p,0)ans =-1>> polyval(p,[-2,-1,0,1,2])ans =3 0 -1 0 3>> fplot('polyval([ 1, 0, -1] , x)' , [-3,3]), grid>>Figure 3 : Graphe de la fonction polynomiale p(x) = x 2 -1.24
On obtient les racines <strong>du</strong> polynôme p grâce à l'instruction roots(p).L'instruction poly permet d'obtenir la représentation canonique p(x) = a n x n + a n-1 x n-1 + . . . + a 1 x + a 0d'un polynôme de degré n dont on connaît les n+1 racines {x i } {i=0,...,n} . Les coefficients {a i } {i=0,...,n} sontobtenus sous forme d'un vecteur p et sont rangés dans l'ordre décroissant des indices. C'est-à-dire quep(1) = a n , . . . , p(n+1) = a 0 .>> r = roots(p)r =-1.00001.0000>> poly(r)ans =1.0000 0.0000 -1.0000>>4.1 Les formes d’affichage des réels4. Les entrées-sortiesMatlab dispose de plusieurs formats d'affichage des réels. Par défaut le format est le format court àcinq chiffres. Les autres principaux formats sont :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.Matlab dispose également des formats format short g et format long g qui utilise la « meilleure » desdeux écritures à virgule fixe ou à virgule flottante. On obtiendra tous les formats d'affichage possiblesen tapant help format. On impose un format d'affichage en tapant l'instruction de formatcorrespondante dans la fenêtre de contrôle, par exemple format long. Pour revenir au format pardéfaut on utilise la commande format ou format short.>> pians =3.1416>> format long>> pians =3.14159265358979>> format short e>> pi^3ans =3.1006e+01>> format short g>> pi^3ans =31.006>> format short>>25
- Page 1:
T.P. Identification, Commandeet Sta
- Page 4 and 5:
2 2dh Kim = mg−2dt héquation mé
- Page 6 and 7: 4. Etude similaire pour une balle r
- Page 9 and 10: TP 2 : Simulation de trajectoires d
- Page 11 and 12: TP 3 : Application de la méthode d
- Page 13: Remarquons que nous avons besoin de
- Page 16 and 17: 2. Description du procédéOn consi
- Page 18 and 19: En fin d’exécution le fichier vo
- Page 20 and 21: 5. Calcul d’un régulateur némat
- Page 22 and 23: 3. Modélisation du systèmeUn mod
- Page 24 and 25: D :Dépassementω r0.95KsKs:Gaine s
- Page 26 and 27: ( )()( )( )() ( )( )( 0)() 1⎡ y 0
- Page 28 and 29: Figure 7 : Un signal SBPA et sa fon
- Page 30 and 31: Le tableau suivant montre quelques
- Page 33 and 34: Annexe : Généralités sur Matlab1
- Page 35 and 36: y =0.8660z =-0.5000>> save toto y z
- Page 37 and 38: Plusieurs instructions Matlab peuve
- Page 39 and 40: 1>> isreal(rep)ans =1>>2.1.2 Le typ
- Page 41 and 42: y =1.229999999999999e+02tst =0>> fo
- Page 43 and 44: and(m,1): vecteur colonne de longue
- Page 45 and 46: ans =1 11 11 1>> zeros(2)ans =0 00
- Page 47 and 48: On n'est pas obligé de se limiter
- Page 49 and 50: (1,1) 2(2,1) -1(1,2) -1(2,2) 2(3,2)
- Page 51 and 52: infNaN: infini. Est obtenu quand on
- Page 53 and 54: vecteurs. Il s'obtient grâce à l'
- Page 55: Les fonctions matricielles les plus
- Page 59 and 60: 4.4 Impressions dirigées par forma
- Page 61 and 62: 5. Programmer sous Matlab5.1 Script
- Page 63 and 64: nombre de variables d'entrée utili
- Page 65 and 66: • séquence d'instructions est le
- Page 67 and 68: Interprétation :La séquence d'ins
- Page 69 and 70: séquence d'instructions 2...case c
- Page 71 and 72: if nargout == 2rang = rank(A);if na
- Page 73 and 74: %% Sortie :% ly : indice de ligne d
- Page 75 and 76: On obtient alors le graphe de la fo
- Page 77 and 78: ouge + plus -. pointillé mixteg ve
- Page 79 and 80: supposons que la variable numex con
- Page 81 and 82: Figure 11 : Fenêtre graphique déc
- Page 83 and 84: [C,h] = contour(X,Y,Z,n)>> clabel(C
- Page 85 and 86: Dans les deux cas on obtient le ré
- Page 87 and 88: Figure 19 : Exemple de visualisatio