Matrices

[rafaelcrossover]

Bonjour,
Je crée un programme qui me permet de former une matrice globale à partir de plusieurs (2 ou 3) matrices plus petites. Le problème est que je l’obtiens de manière désordonnée et j’aimerais savoir comment je peux le commander sans changer le résultat.

Petites matrices:

Matrice que je trouve:

Ce que je veux, c’est que dans la position 1.1 de la matrice que je reçois, c’est la somme des éléments présents dans les petites matrices ayant le même indice de 1.1. Par exemple, dans ce cas, dans la première matrice, l'élément 1.1 est égal à 0,667 et dans la deuxième matrice, à 1 car dans la grande matrice en position 1, il devrait être de 1,67, et il l'est, mais pas en position 1.1. De la grande matrice, que puis-je faire pour résoudre cela?

[Bisam]

Tu ne donnes pas le contexte et il faut un peu deviner ce que tu veux faire.
J'ai l'impression que la première ligne et la première colonne de tes matrices sont juste des "noms" qu'il faut rechercher et faire un agrégat des résultats.
La structure matrice n'est certainement pas la meilleure pour faire cela... mais tu n'as pas vraiment d'autre moyen sur ta calculette (à moins d'utiliser Ndless et Micropython)
Bref, voici un algorithme qui doit convenir :

Code: Select all

Define combine(mat1,mat2)=Func
Local mat, noms, indices, i, j, c, i1, i2, n1, n2, n, again

n1 := rowdim(mat1)
n2 := rowdim(mat2)
n := n1-1

# on fait la liste des "noms" présents dans la première matrice
# et on note leur place dans la première ligne de la matrice "indices"
noms := mid(mat►list(mat1[1]),2)
indices := {seq(i, i, 2, n1), seq(0, i, 2, n1)}

# on ajoute les "noms" présents dans la deuxième matrice
# et on note leur place dans la deuxième ligne de la matrice "indices"
For i, 2, n2
  j := 1
  nom := mat1[1,i]
  again := true
  While again
    If j > n Then
      again := false
    ElseIf mat2[1,j] = nom
      indices[2,j] := i
      again := false
    EndIf
    j := j+1
  EndWhile
  If j > n Then
    noms[j] := nom
    indices := augment(indices, {{0}, {j}})
    n := n+1
  EndIf
EndFor

# On fabrique la nouvelle matrice avec les informations récoltées dans les deux matrices
mat := newmat(n+1,n+1)
For i, 1, n
  nom := noms[i]
  mat[1,i] := nom
  mat[i,1] := nom
  i1 := indices[1,i]
  i2 := indices[2,i]
  If i1 ≠ 0 Then
    For j, 2, n1
      c := mat1[j,i1]
      mat[j,i] := mat[j,i] + c
      If j ≠ i
        mat[i,j] := mat[i,j] + c
    EndFor
  EndIf
  If i2 ≠ 0 Then
    For j, 2, n2
      c := mat2[j,i2]
      mat[j,i] := mat[j,i] + c
      If j ≠ i
        mat[i,j] := mat[i,j] + c
    EndFor
  EndIf
EndFor
Return mat
EndFunc


Je n'ai pas vérifié que ça marche, mais c'est facile à comprendre et tu pourras corriger sans peine si j'ai raté un truc.