3

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 2

4

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.]]
2

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)

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .