Я старался быть внимательным в просмотре документации и выхожу пустым.

Я пытаюсь разложить или исключить термины из выражения в матричную форму. Моя проблема, кажется, отличается от полиномиального факторинга (поскольку я планирую реализовать функциюphi(x,y,z) = a_1 + a_2*x + a_3*y + a_4*z)

import sympy
from sympy import symbols, pprint
from sympy.solvers import solve

phi_1, phi_2, x, a_1, a_2, L = symbols("phi_1, phi_2, x, a_1, a_2, L")

#Linear Interpolation function: phi(x)
phi = a_1 + a_2*x
#Solve for coefficients (a_1, a_2) with BC's: phi(x) @ x=0, x=L
shape_coeffs = solve([Eq(phi_1, phi).subs({x:0}), Eq(phi_2, phi).subs({x:L})], (a_1, a_2))
pprint(shape_coeffs)
#Substitute known coefficients
phi = phi.subs(shape_coeffs)
pprint(phi)

Это работает, как и ожидалось, однако я хотел бы перейти к матричной форме, где:

я пыталсяfactor(), cancel(), as_coefficient() без успеха. На бумаге это тривиальная проблема. Чего мне не хватает в решении Simpy? Благодарю.

Действительный метод: взятый как ответ

C_1, C_2 = symbols("C_1, C_2", cls=Wild)
N = Matrix(1,2, [C_1, C_2])
N = N.subs(phi.match(C_1*phi_1 + C_2*phi_2))
phi_i = Matrix([phi_1, phi_2])
display(Math("\phi(x)_{answered} = " + latex(N) + "\ * " + latex(phi_i)))