xameus / mn Goto Github PK

Index invalide.
at line       8 of function descente called by :  
Z = descente(l, m, Y')

J'ai mis un petit disp(i) pour voir quand est-ce que ça plante : au dernier rang.
C'est à cause de : m(i) (m étant en fait la matrice sous diagonale)

l est la diagonale et m la sous-diagonale.

Il manque des étapes dans l'exercice 2, dans le rendu.
Serait-il souhaitable de les rajouter?

Il manque un graphe pour la question 12

Normalement dans LZ = Y, L est une matrice "carrée" (car LL), Y est indiqué comme étant un "vecteur".

Quel type utiliser pour Z et, de façon plus générale, pour les vecteurs ?

Je propose d'utiliser une matrice à 1 colonne n lignes pour représenter les vecteurs, car cela permet des produits avec des matrices nxn (contrairement à des matrices à 1 ligne n colonnes).

N.B. : Il suffirait donc de déclarer plutôt les vecteurs comme : Y = [1; 2; 3; 4; 5] au lieu de Y = [1, 2, 3, 4, 5]

On peut constater qu'en augmentant n et n_t à au moins 1000 et qu'en mettant un nombre d'itération arbitrairement grand dans la boucle, nos résultats deviennent excellents.
A mes yeux, il serait souhaitable de:

• n, n_t à 1000

• Le nombre d'itérations à 10000 (en affichant toutes les 1000 itérations par exemple)

Ne marche pas...

Je pense que c'est en raison des variables globale. Mais étrange qu'elle ne fonctionne pas.

Je pense que l'on pourrait refactoriser le code afin de le clarifier.

• Il y a beaucoup de redondance dans le code
• A priori Scilab ne gère pas les lambda fonctions, qui sont ici réellement utiles.
• Nous avons beaucoup de fonctions qui dépendent de variables non-passées dans les paramètres.
• Nous avons tout dans un seul et même fichier, avec les tests unitaires à l'intérieur, et à chaque exécution nous exécutons tout le fichier.

Devrait-t-on y remédier?