Discussion :

[SpacedCb]

Bonjour,

Je n'arrive pas à créer un routine me permettant de créer une matrice de ce type:

[111 000 000 ; 000 111 000; 000 000 111] pour k = 3
[1111 0000 0000 0000 ; 0000 1111 0000 0000 0000 ; 0000 0000 1111 0000; 0000 0000 0000 1111] pour k = 4

Et ainsi de suite pour n'importe quel k.

j'en ai besoin pour une minimisation de fonction. Je n'ai vraiment pas envie de devoir réorganiser les paramètres à optimiser, ça changerai vraiment trop le code.
Avez vous un idée parce que là je sèche. Il doit bien y a voir une fonction qui fait ce genre de truc non?

Merci d'avance,

cdlt,

[Jean Dumoncel]

Bonjour,

regarde du côté de la fonction blkdiag.

[Jerome Briot]

magelan a donné une bonne solution, en voici une autre :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

diag(ones(1,k)*sum(10.^(0:k-1)))

ou encore :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
M = zeros(k);
M(1:k+1:end) = sum(10.^(0:k-1))

Il nous faudrait néanmoins plus de précision sur le type de ce tableau.
Avec des valeurs numériques, tu ne pourras au mieux obtenir que [111 0 0 ; 0 111 0; 0 0 111] pour k = 3

[SpacedCb]

blkdiag c'est parfait! Merci beaucoup.

Pour ta soluce Dut, le problème c'est qu'il me faut le nombre même nombre de 0 que de 1.
Juste une question au passage, comment vous trouvez ces fonctions? Seulement avec l'aide de Matlab?

Merci énormément en tout cas, je vais pouvoir rajouter mes contraintes dans ma minimisation!

Bonne soirée/nuit

[Jean Dumoncel]

Pour ta soluce Dut, le problème c'est qu'il me faut le nombre même nombre de 0 que de 1.

Donc pour k=3, la matrice est :
[1 1 1 0 0 0 0 0 0 ; 0 0 0 1 1 1 0 0 0; 0 0 0 0 0 0 1 1 1]
c'est bien cela?

Juste une question au passage, comment vous trouvez ces fonctions? Seulement avec l'aide de Matlab?

Pour ma part, soit avec matlab :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

>> docsearch block diagonal matrix

soit sur internet en rajoutant le mot clé matlab.

[SpacedCb]

ok ok, j'avais juste pas les bons mots clefs en fait


Merci en tout cas!