Uma dúvida, para a solução de sistemas lineares: como concatenar um array (matriz) A, um array (vetor coluna) b, de forma que se tenha a matriz "aumentada" do sistema, A~ = [A b], usando numpy?
2 Respostas
Sugestão de solução (as matrizes de exemplo A
e b
poderiam ser criadas de forma mais automática, mas assim eu achei que o exemplo ficaria mais didático):
import numpy as np
A = np.array([ [1, 2, 3, 4], [5, 6, 7, 8], [9, 10, 11, 12], [13, 14, 15, 16] ])
b = np.array([17, 18, 19, 20])
# Cria uma nova matriz com 1 coluna a mais
Ab = np.zeros((A.shape[0], A.shape[1]+1))
# Copia a matriz original para a nova matriz
Ab[:,:-1] = A
# Copia o vetor, convertido em uma matriz de uma coluna, para a nova matriz
Ab[:,-1:] = b.reshape(A.shape[0], 1)
print('A:')
print(A)
print('')
print('b:')
print(b)
print('')
print('[A b]:')
print(Ab)
Saída desse código:
A:
[[ 1 2 3 4]
[ 5 6 7 8]
[ 9 10 11 12]
[13 14 15 16]]
b:
[17 18 19 20]
[A b]:
[[ 1. 2. 3. 4. 17.]
[ 5. 6. 7. 8. 18.]
[ 9. 10. 11. 12. 19.]
[ 13. 14. 15. 16. 20.]]
Outra maneira é:
import numpy as np
a = np.array([ [1, 2, 3, 4], [5, 6, 7, 8], [9, 10, 11, 12], [13, 14, 15, 16] ])
b = np.array([[17, 18, 19, 20]])
out = np.concatenate((a,b.T), axis = 1)