12.07.2015 Views

Décomposition d'une fraction rationnelle en ... - xavierdupre.fr

Décomposition d'une fraction rationnelle en ... - xavierdupre.fr

Décomposition d'une fraction rationnelle en ... - xavierdupre.fr

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Décomposition <strong>d'une</strong> <strong><strong>fr</strong>action</strong> <strong>rationnelle</strong><strong>en</strong> élém<strong>en</strong>ts simples


REPÈRES 2Repères1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32. Préliminaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42.1 Opérations sur les polynômes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42.2 Opérations sur les matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53. Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103.1 Principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103.2 Construction du système linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103.3 Un exemple de programme informatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124. Le programme informatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144.1 Tableaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144.2 Structures possibles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144.3 Expression de la décomposition à l'aide de ces structures . . . . . . . . . . . . . . . . . . . 15La table des matières est détaillée à la n du docum<strong>en</strong>t.


1. Introduction 3Chapitre 1IntroductionIl existe une méthode systématique qui permet de décomposer une <strong><strong>fr</strong>action</strong> rationelle <strong>en</strong> élém<strong>en</strong>ts simples :la première étape est d'écrire la forme a priori de la décomposition, puis, après mise sous même dénominateur,d'id<strong>en</strong>tier les coeci<strong>en</strong>ts de la décomposition. Cette méthode revi<strong>en</strong>t à résoudre un système linéairede n inconnues à n équations. Les n inconnues sont les n coeci<strong>en</strong>ts de la décomposition de la <strong><strong>fr</strong>action</strong><strong>en</strong> élém<strong>en</strong>ts simples. Un théorème arme que cette décomposition existe et est unique ce qui nous assureque le système linéaire à résoudre admet nécessairem<strong>en</strong>t une solution unique.Les chapitres qui suivront ont pour de prés<strong>en</strong>ter un algorithme de décomposition des <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s<strong>en</strong> élém<strong>en</strong>ts simples.


2. Préliminaires 4Chapitre 2Préliminaires2.1 Opérations sur les polynômesLa multiplication de polynômes sera la seule opération sur les polynômes utilisée par l'algorithme. Néanmoins,l'addition sera égalem<strong>en</strong>t prés<strong>en</strong>tée.Un polynôme P est déni par :D'après ces notations :{ un degré dPun tableau de coeci<strong>en</strong>ts {a 0 , · · · ,a dP }P (x) =Par dénition, un polynôme est constant si et seulem<strong>en</strong>t si d P = 0, par conv<strong>en</strong>tion, le polynôme nul aurapour degré 0.2.1.1 L'additionOn considère deux polynômes :P de degré d P et de coeci<strong>en</strong>ts {a 0 , · · · ,a dP } et Q de degré d Q et de coeci<strong>en</strong>ts { b 0 , · · · ,b dQ}.Soit R = P + Q, on sait que d R max (d P ,d Q ), et es coeci<strong>en</strong>ts de R sont {c 0 , · · · ,c dR }.⎧⎨∀i ∈ {0, · · · ,d R } , c i =⎩La valeur précise de d R est donnée par :d P ∑i=0a i x ia i + b i si i min {d P ,d Q }a i si i d P et i > d Qb i si i > d P et i d Qd R ={ max {dP ,d Q } si d P ≠ d Qinf {i | i ∈ {0, · · · ,d R } et ∀k > i, c k = 0} si d P ≠ d Q


2. Préliminaires 5On pr<strong>en</strong>d cette précaution pour s'assurer que le coeci<strong>en</strong>t du terme de plus haut degré est non nul. Dansle programme, cette addition n'est pas utilisée.2.1.2 La multiplicationOn considère deux polynômes P de degré d Pet de coeci<strong>en</strong>ts {a 0 , · · · ,a dP } et Q de degré d Q et decoeci<strong>en</strong>ts { b 0 , · · · ,b dQ}.P =Q =d P ∑i=0d Q∑i=0a i x ib i x iPar conséqu<strong>en</strong>t, soit R le polynôme vériant R = P Q, de degré d R et de coeci<strong>en</strong>ts {c 0 , · · · ,c dR }, alors :Dans le cas où P ≠ 0 et Q ≠ 0 :d R ={ 0 si P = 0ouQ = 0d P + d Q sinon( dP) ⎛ ∑R = a i x i ⎝R =R =i=0d R ∑∑i=0 k+l=ikd Rd R ∑i=0d Q∑i=0ld Qa k b l x imin{i,d P }∑k=max{0,i−d Q}⎞b i x i ⎠a k b i−k x i2.2 Opérations sur les matricesL'algorithme fait interv<strong>en</strong>ir deux opérations sur les matrices : la multiplication et l'inversion.Une matrice M est représ<strong>en</strong>tée de manière informatique par :2.2.1 La multiplication⎧⎪⎨⎪⎩un nombre de lignes L Mun nombre de colonnes C Mun tableau de coeci<strong>en</strong>ts (a ij ) 1iLM1jC MOn considère deux matrices : M de coeci<strong>en</strong>ts (a ij ) 1iLM1jC Met N de coeci<strong>en</strong>ts (b ij ) 1iLN1jC N.


On peut eectuer le produit K = MN si C M = L N :⎧⎪⎨⎪⎩2. Préliminaires 6K = (c ij ) 1iLM1jC N∀i ∈ {1, · · · ,L M } , ∀j ∈ {1, · · · ,C N } , c ij = C M∑2.2.2 Multiplication et opérations sur les lignes et colonnesk=0a ik b kjLes multiplications <strong>en</strong>tre matrices permett<strong>en</strong>t à l'aide de matrices simples de traduire des opérations surles lignes et les colonnes.Exemple :⎛⎝1 0 00 1 −52 0 1⎞⎠⎛⎝⎞L 1L 2⎠L 3} {{ }matrice quelconquecomposée de 3 lignes⎛= ⎝⎞L 1L 2 − 5L 3⎠L 3 + 2L 1Dans le cas général, multiplier à gauche par une matrice carrée M ∈ M n , M = (a ij ) 1indes opérations sur les lignes sur la matrice de droite :⎛(a ij ) ⎜⎝⎛⎞L 1L 2⎟. ⎠ = L n⎜⎝n∑⎞a 1j L jj=1n∑a 2j L jj=1.⎟n∑ ⎠a nj L jj=11jnrevi<strong>en</strong>t à faireDe même, multiplier à droite par une matrice carrée M ∈ M n ,opérations sur les colonnes sur la matrice de gauche :(( ) n∑C1 C 2 . . . C n (aij ) = a i1 C ii=1n∑a i2 C i . . .i=1M = (a ij ) 1in1jn)n∑a in C ii=1revi<strong>en</strong>t à faire desCes précisions sont importantes <strong>en</strong> terme de coût algorithime . En eet, pour eectuer le produit dematrices suivant :⎛⎝1 0 00 1 −52 0 1⎞⎠⎛⎝⎞L 1L 2⎠L 3} {{ }matrice quelconquecomposée de 3 lignesIl est plus judicieux de faire des opérations sur les lignes plutôt que d'eectuer le produit matriciel complet.Le coût d'un produit matriciel est <strong>en</strong> O ( n 3) ceci signie que pour deux matrices carrées d'ordre n, le


2. Préliminaires 7nombre d'opérations (au s<strong>en</strong>s informatique) est un multiple constant (indép<strong>en</strong>dant des matrices choisies)de n 3 , ou 3 boucles imbriquées de n itérations. Or, eectuer les opérations suivantes sur les lignes :⎛⎝⎞L 1L 2 − 5L 3⎠L 3 + 2L 1est beaucoup moins coûteux. Dans ce cas, chaque ligne de la matrice résultat est une somme coeci<strong>en</strong>téede deux autres lignes. Le coût de ce produit sera donc <strong>en</strong> O ( n 2) . Ce sera justem<strong>en</strong>t le seul cas r<strong>en</strong>contrélors du pivot de Gauss.2.2.3 Matrice de passage du pivot de GaussLors du calcul de l'inverse, il est nécessaire de calculer la matrice de passage qui permet de passer de lamatrice initiale à la matrice triangulaire supérieure obt<strong>en</strong>ue par pivot de Gauss.⎛1 0⎞ ⎛0⎝ −2 1 0 ⎠ ⎝}−3 −5{{1}matrice de passageOn peut écrire égalem<strong>en</strong>t :⎛⎝1 0 0−2 1 0−3 −5 1⎞ ⎛⎠ ⎝Si on note la matrice id<strong>en</strong>tité :⎞L 1L 2⎠ =L 31 0 00 1 00 0 1} {{ }matrice id<strong>en</strong>tité⎛⎝⎞L 1L 2 − 2L 1⎠L 3 − 3L 1 − 5L 2} {{ }obt<strong>en</strong>ue directem<strong>en</strong>t par pivot de Gauss⎞ ⎛⎠ ⎝⎛⎝⎞ ⎛L 1L 2⎠ = ⎝L 31 0 00 1 00 0 1On peut écrire le résultat précéd<strong>en</strong>t sous la forme :⎛⎝⎞ ⎛I 1I 2 − 2I 1⎠ ⎝I 3 − 3I 1 − 5I 2⎞ ⎛L 1L 2⎠ = ⎝L 3⎞⎛= ⎝⎞ ⎛L 1L 2 − 2L 1⎠ = ⎝L 3 − 3L 1 − 5L 2⎛⎠ = ⎝⎞I 1I 2⎠I 3⎞ ⎛L 1L 2 − 2L 1⎠ = ⎝L 3 − 3L 1 − 5L 2⎞a 11 a 12 a 130 a 22 a 23⎠0 0 a 33⎞a 11 a 12 a 130 a 22 a 23⎠0 0 a 33⎞a 11 a 12 a 130 a 22 a 23⎠0 0 a 33Le pivot de Gauss consiste à appliquer une série d'opérations sur les lignes à partir de la matrice initialepour la r<strong>en</strong>dre triangulaire supérieure. La matrice de passage s'obti<strong>en</strong>t donc <strong>en</strong> eectuant exactem<strong>en</strong>t lesmêmes opérations sur les lignes à partir de la matrice id<strong>en</strong>tité.2.2.4 L'inversionL'inversion <strong>d'une</strong> matrice utilise le pivot de Gauss. Cette méthode permet de r<strong>en</strong>dre les matrices triangulairessupérieures. C'est égalem<strong>en</strong>t cette méthode qui est utilisée pour calculer son déterminant.La méthode utilisée se décrit simplem<strong>en</strong>t, nous supposons que nous avons une matrice de la forme :


2. Préliminaires 8⎛ ⎞ ⎛L 1L 2⎜ L 3⎟⎝ L 4⎠ = ⎜⎝L 5a 11 ∗ ∗ ∗ ∗0 a 22 ∗ ∗ ∗0 0 a 33 a 34 ∗0 0 a 43 a 44 ∗0 0 a 53 a 54 ∗⎞⎟⎠les * désign<strong>en</strong>t un coeci<strong>en</strong>t quelconque)L'objectif est de r<strong>en</strong>dre nul les coeci<strong>en</strong>ts a 43 et a 53 , pour ce faire, on eectue les opérations sur leslignes comme si on résolvait un système linéaire, ces opérations correspond<strong>en</strong>t égalem<strong>en</strong>t à un produit dematrice :On suppose que a 33 ≠ 0 :P 3L 4 ← L 4 − L 3a 43a 33L 5 ← L 5 − L 3a 53a 33⎛⎞⎛⎞1 0 0 0 0 a 11 ∗ ∗ ∗ ∗0 1 0 0 00 a 22 ∗ ∗ ∗⎜ 0 0 1 0 0⎟⎜0 0 a 33 a 34 ∗⎟⎝ 0 0 −a 43 1 0 ⎠⎝0 0 a 43 a 44 ∗ ⎠0 0 −a 53 0 1 0 0 a 53 a 54 ∗} {{ }} {{ }M 3=⎛⎜⎝a 11 ∗ ∗ ∗ ∗0 a 22 ∗ ∗ ∗0 0 a 33 a 34 ∗⎞⎟⎠0 0 0 a 44 − a 4333a 34 ∗}0 0 0 a 54 − a 53a 33a 44{{∗}M 4Dans le cas où a 33 = 0 et que a 43 ≠ 0, on fait d'abord : L 3 ← L 4 puis on continue le pivot de Gauss :⎛⎞ ⎛⎞⎛1 0 0 0 0 1 0 0 0 0 a 11 ∗ ∗ ∗ ∗0 1 0 0 00 1 0 0 00 a 22 ∗ ∗ ∗⎜ 0 0 1 0 0⎟ ⎜ 0 0 1 1 0⎟⎜0 0 0 a 34 ∗⎟ =⎝ 0 0 −a 43 1 0 ⎠ ⎝ 0 0 0 1 0 ⎠⎝0 0 a 43 a 44 ∗ ⎠0 0 −a 53 0 1 0 0 0 0 1 0 0 a 53 a 54 ∗} {{ }} {{ }P 3M⎛3⎞⎞a 11 ∗ ∗ ∗ ∗0 a 22 ∗ ∗ ∗⎜ 0 0 a 43 a 43 + a 44 ∗⎟⎝ 0 0 0 −a 34 ∗ ⎠}0 0 0 a 54 − a 53a 43a 34{{∗}M 4Dans le cas général, on suppose que M est une matrice carrée d'ordre n, on suppose que M 1 = M et lasuite (P i ) est construite de telle sorte que :P i M i = M( i+1)M i =∀ (k,l) tel que i > k > l 1, a kli = 0a kli1k,ln


2. Préliminaires 9On obti<strong>en</strong>t nalem<strong>en</strong>t que (P n P n−1 . . . P 1 ) M = N avec N matrice triangulaire supérieure. Comme lamatrice P n P n−1 . . . P 1 est inversibles (c'est un produit de matrice inversibles), la matrice M est inversiblesi la matrice N l'est aussi et N est inversible si tous ses coeci<strong>en</strong>ts sur la diagonale sont non nuls.En eet, on note N = ( b kl) , N est triangulaire supérieure donc :1k,lndet N =n∏i=1b ii= det (P n P n−1 . . . P 1 ) × det M} {{ }=1 par construction= det MOn suppose maint<strong>en</strong>ant que la matrice M est inversible, tous les coeci<strong>en</strong>ts diagonaux de N sont doncnon nuls. Dans cette première phase, nous avons, dans l'ordre croissant des lignes, rempli la matrice Mde 0 <strong>en</strong> dessous de sa diagonale, il est alors possible <strong>en</strong> appliquant le même pivot de Gauss mais cette foisdans l'ordre décroissant des lignes de remplir la matrice de 0 au-dessus de sa diagonale.Il existe donc une suite de matrice (Q i ) tel que (Q n Q n−1 . . . Q 1 ) N = D n où D n est une matrice diagonale.De plus, on peut choisir les matrices Q i de telle sorte que D n = I n où I n est la matrice id<strong>en</strong>tité. Parexemple :Donc :D'où :⎛⎞ ⎛⎞⎛1 0 −b 13 0 0 1 0 0 0 0 b 11 b 12 b 13 0 00 1 −b 23 0 00 1 0 0 00 b 22 b 23 0 0⎜ 0 0 1 0 01⎟ ⎜ 0 0b 330 0⎟⎜0 0 b 33 0 0⎟ =⎝ 0 0 0 1 0 ⎠ ⎝ 0 0 0 1 0 ⎠⎝0 0 0 1 0 ⎠0 0 0 0 1 0 0 0 0 1 0 0 0 0 1} {{ }} {{ }Q 3N⎛3⎞⎜⎝(Q n Q n−1 . . . Q 1 ) (P n P n−1 . . . P 1 ) M = I nM −1 = (Q n Q n−1 . . . Q 1 ) (P n P n−1 . . . P 1 )⎞b 11 b 12 0 0 00 b 22 0 0 00 0 1 0 00 0 0 1 00 0 0 0 1⎟⎠} {{ }N 4


3. Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s 10Chapitre 3Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s3.1 PrincipeTrouver tous les coeci<strong>en</strong>ts de la décomposition a priori <strong>d'une</strong> fonction <strong>rationnelle</strong> revi<strong>en</strong>t à résoudre unsystème linéaire.Soit P et Q deux polynômes à coeci<strong>en</strong>ts réels, Q se décompose <strong>en</strong> :La <strong><strong>fr</strong>action</strong> P Qse décompose donc <strong>en</strong> :s∏Q = α (X − a i ) d ii=1t∏ (X 2 ) d ′+ p i X + q i ii=1PQ = E +s∑ ∑d ioù E est la partie <strong>en</strong>tière.On met tout sous le même dénominateur :P = EQα k d t∑ii=1(X − ak=1 i ) k + ∑′ ii=1 k=1βi kX + γk i(X 2 + p i X + q i ) k( ) (∑+ α s ∏(X + a j ) d ∏ (j X 2 ) d ′ di)+ p i X + q j∑i αi k (X − a i) d i−ki=1 j≠ijk=1+ α t ∑i=1(∏j(X + a j ) d j∏j≠i) ((X 2 ) d ′ d ′+ p i X + q j∑ iik=1(βki X + γik ) (X 2 ) d ′+ p i X + q i −kiIl sut d'id<strong>en</strong>tier ces deux polynômes, ce qui revi<strong>en</strong>t à résoudre un système linéaire. Pour que l'algorithmemarche, il faut néanmoins connaître la factorisation <strong>en</strong> élém<strong>en</strong>ts simples du dénominateur.)3.2 Construction du système linéaireL'illustration par un exemple permet de mieux compr<strong>en</strong>dre l'algorithme qu'une description théoriqueutilisant une multitude d'indice.


3. Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s 11Nous allons donc décortiquer chaque étape de l'algorithme sur la <strong><strong>fr</strong>action</strong> <strong>rationnelle</strong> :3.2.1 Numérotation des coeci<strong>en</strong>tsLa partie <strong>en</strong>tière est de degré 2 (=9-7).La forme a priori de la décomposition de la <strong><strong>fr</strong>action</strong> <strong>rationnelle</strong> est :PQ = X 9(3.1)(X + 1) 3 (X 2 + X + 2) 2PQ = a 1 + a 2 X + a 3 X 2 + a 4X + 1 + a 5(X + 1) 2 + a 6(X + 1) 3 + a 7 + a 8 XX 2 + X + 2 + a 9 + a 10 X(X 2 + X + 2) 2La numérotation suit le s<strong>en</strong>s croissant des puissances <strong>en</strong> comm<strong>en</strong>çant par la partie <strong>en</strong>tière.3.2.2 Construction du système linéaireEn mettant tout sous le même dénominateur, on obti<strong>en</strong>t :P = a 1 Q + a 2 XQ + a 3 X 2 Q +On construit la matrice suivante :⎛⎛⎞a 1 × Qa 2 × XQa 3 × X 2 Qa 4 × Q/ (X + 1)a 5 × Q/ (X + 1) 2=a 6 × Q/ (X + 1) 3a 7 × Q/ ( X 2 + X + 2 )a⎜ 8 × XQ/ ( X 2 + X + 2 )⎝ a 9 × Q/ ( X 2 + X + 2 ) 2 ⎟a 10 × XQ/ ( X 2 + X + 2 ) ⎠ ⎜⎝2a 4QX + 1 + a 5Q(X + 1) 2 + a 6Q(X + 1) 3 + (a 7 + a 8 X) QX 2 + X + 2 + (a 9 + a 10 X) Q(X 2 + X + 2) 2X 0 X 1 X 2 X 3 X 4 X 5 X 6 X 7 X 8 X 9× × × × × × × × × ×a 1 × 4 16 29 33 26 14 5 1a 2 × 4 16 29 33 26 14 5 1a 3 × 4 16 29 33 26 14 5 1a 4 × 4 12 17 16 10 4 1a 5 × 4 8 9 7 3 1a 6 × 4 4 5 2 1a 7 × 1 4 6 4 1a 8 × 1 4 6 4 1a 9 × 1 3 3 1a 10 × 1 3 3 1} {{ }=M (10 lignes et 10 colonnes)On peut écrire ce système sous la forme : (M T est la transposée de M)⎞⎟⎠


3. Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s 12⎛( ) a1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 a 9 a 10 M ⎜⎝Donc par id<strong>en</strong>tication des deux polynômes, on obti<strong>en</strong>t :X 0 ⎞X 1X 2X 3X 4X 5X 6X 7⎟X 8 ⎠X 9=⎛⎜⎝0000000001⎞⎟⎠T ⎛⎜⎝X 0 ⎞X 1X 2X 3X 4X 5X 6X 7⎟X 8 ⎠X 9} {{ }=P numérateur de la <strong><strong>fr</strong>action</strong>(a1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 a 9 a 10)M =(0 0 0 0 0 0 0 0 0 1)D'où :M T ( a 1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 a 9 a 10) T =(0 0 0 0 0 0 0 0 0 1) T} {{ }c'est le système linéaire cherché3.2.3 RésolutionD'après le système précéd<strong>en</strong>t, on déduit que :(a1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 a 9 a 10) T =(MT ) −1 (0 0 0 0 0 0 0 0 0 1) TLa décomposition de la <strong><strong>fr</strong>action</strong> <strong>rationnelle</strong> est donc :PQ = 11 − 5X + X2 − 6,6875X + 1 + 2(X + 1) 2 − 0,25 15 + 4,3125X 5,75 + 5,625X3−(X + 1) X 2 ++ X + 2 (X 2 + X + 2) 23.3 Un exemple de programme informatiqueLa <strong><strong>fr</strong>action</strong> ( 3.1) peut être décrite comme dans la gure 3.1, le résultat est imprimé dans la gure 3.2.Quelques diér<strong>en</strong>ces peuv<strong>en</strong>t être observées <strong>en</strong>tre la solution fournies par le programme et la véritablesolution, ceci est dû aux approximations de calcul des ordinateurs qui code les nombres réels avec <strong>en</strong>viron15 chires de précision. Certains termes nuls peuv<strong>en</strong>t alors être non nuls mais très petits par rapport auxautres coeci<strong>en</strong>ts.


3. Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s 13numerateur degre 9d0 0d1 0d2 0d3 0d4 0d5 0d6 0d7 0d8 0d9 1..................d<strong>en</strong>ominateur produit 2d<strong>en</strong>ominateur puissance 3d<strong>en</strong>ominateur degre 1d0 1d1 1d<strong>en</strong>ominateur puissance 2d<strong>en</strong>ominateur degre 2d0 2d1 1d2 1Fig. 3.1 Description <strong>d'une</strong> <strong><strong>fr</strong>action</strong> <strong>rationnelle</strong><strong><strong>fr</strong>action</strong> a decomposer :numerateur : X ∧ 9nombre de terme au d<strong>en</strong>ominateur : 2terme 1 : (1 + X) ∧ 3terme 2 : (2 + X + X ∧ 2) ∧ 2decomposition <strong>en</strong> elem<strong>en</strong>ts simplespartie <strong>en</strong>tiere : 11 − 5 ∗ X + 1 ∗ X ∧ 2terme <strong>en</strong> (1 + X) : −6.6875terme <strong>en</strong> (1 + X) ∧ 2 : 2terme <strong>en</strong> (1 + X) ∧ 3 : −0.25terme <strong>en</strong>terme <strong>en</strong>Fig. 3.2 Résultat de la décomposition(2 + X + X ∧ 2) : −15 − 4.3125 ∗ X(2 + X + X ∧ 2) ∧ 2 : 5.75 + 5.625 ∗ X


4. Le programme informatique 14Chapitre 4Le programme informatique4.1 TableauxEtant que les tableaux informatique de dim<strong>en</strong>sion n sont indicés de 0 à n − 1 inclus, un polynôme devi<strong>en</strong>t :Et une matrice :⎧⎧⎪⎨⎪⎩a [50] est un tableau appelé a à 50 élém<strong>en</strong>ts numérotés de 0 à 49a[i] représ<strong>en</strong>te l'élém<strong>en</strong>t d'indice iun polynôme s'écrit P = n ∑i=0a [i] X ia [40] [50] est un tableau appelé a de 40 lignes et 50 colonnes⎪⎨⎪⎩a[i][j] représ<strong>en</strong>te l'élém<strong>en</strong>t placé sur la ligne i et la colonne jle dernier élém<strong>en</strong>t est l'élém<strong>en</strong>t a [39] [49]une matrice s'écrit M =⎛⎜⎝a [0] [0] · · · a [0] [n − 1]. · · · .a [m − 1] [0] · · · a [m − 1] [n − 1]⎞⎟⎠ = (a [i] [j]) 0im−10jn−14.2 Structures possiblesDeux types spéciaux peuv<strong>en</strong>t être créés pour cont<strong>en</strong>ir matrices et polynômes.Structure <strong>d'une</strong> matrice :typedef struct {double *C ;int Ligne ;int Colonne ;} Matrice ;


4. Le programme informatique 15Structure d'un polynôme :typedef struct {double *C ;int Degre ;} Polynome ;4.3 Expression de la décomposition à l'aide de ces structuresfacteur−1PQ = partie_<strong>en</strong>tiere+ ∑partie_<strong>en</strong>tiere.DegrePQ = ∑i=0n=0puissance[n]−1∑k=0partie_<strong>en</strong>tiere.C [i] × X i +facteur−1∑n=0puissance[n]−1∑k=0resultat[n][k](elem<strong>en</strong>t[n]) k+1resultat[n][k].Degre∑i=0(elem<strong>en</strong>t[n].Degre∑i=0resultat[n][k].C [i] × X ielem<strong>en</strong>t[n].C [i] × X i ) k+1


BibliographieBibliographie 16


Table des gures 17Table des gures3.1 Description <strong>d'une</strong> <strong><strong>fr</strong>action</strong> <strong>rationnelle</strong> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133.2 Résultat de la décomposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13


Liste des tableauxListe des tableaux 18


Table des matières 19Table des matières1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32. Préliminaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42.1 Opérations sur les polynômes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42.1.1 L'addition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42.1.2 La multiplication . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52.2 Opérations sur les matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52.2.1 La multiplication . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52.2.2 Multiplication et opérations sur les lignes et colonnes . . . . . . . . . . . . . . . . . 62.2.3 Matrice de passage du pivot de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . 72.2.4 L'inversion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73. Application aux <strong><strong>fr</strong>action</strong>s <strong>rationnelle</strong>s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103.1 Principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103.2 Construction du système linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103.2.1 Numérotation des coeci<strong>en</strong>ts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113.2.2 Construction du système linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113.2.3 Résolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123.3 Un exemple de programme informatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124. Le programme informatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144.1 Tableaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144.2 Structures possibles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144.3 Expression de la décomposition à l'aide de ces structures . . . . . . . . . . . . . . . . . . . 15


IndexAaddition. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4Ccoût d'un algorithme . . . . . . . . . . . . . . . . . . . . . . . . . . 6Ddécomposition. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .10Ssystème linéaire . . . . . . . . . . . . . . . . . . . . . . . . . 3, 8, 10Ttableaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14théorème . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3triangulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7Eexemple. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .10expression. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .15F<strong><strong>fr</strong>action</strong> <strong>rationnelle</strong> . . . . . . . . . . . . . . . . . . . . . . . . . . . 10GGauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7Iid<strong>en</strong>tité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7inversion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7Mmatrice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5matrice de passage . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7matrice id<strong>en</strong>tité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7matrice triangulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . 7multiplication. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5, 6Nnumérotation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11Oopérations sur les lignes et les colonnes . . . . . . . . 6Ppivot de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7polynôme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4Rrésolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1220

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!