sympy 求解矩阵方程为空？

如何解决sympy 求解矩阵方程为空？

在这里尝试解决 A @ B = 0 类型的方程，A 是已知矩阵，B 矩阵包含未知数，但代码给出 [] 输出，subs(b,1) 真解必须为 [0.618,1.0 ] 和 [-1.62,1.0] 谢谢。

from scipy import linalg
import numpy as np
import math
from sympy import *
import matplotlib.pyplot as plt

M = np.diag([2,2])  # ton
K = np.array([[600,-300],[-300,300]
              ])

M_inv = linalg.inv(M)
la,v = linalg.eig(M_inv @ K)


wn = np.sqrt(la)
wn.sort()

def matris(w):
    A = K - w**2*M
    return A

matris_1 = matris(wn[0])
matris_2 = matris(wn[1])

a,b,c,d,e,f = symbols('a,f',commutative=False)

fi_matris = Matrix([a,b])


def mod(matris_n):
    eq1 = matris_n @ fi_matris
    eq2 = eq1.subs(b,1)
    eq3 = Eq(eq2,0)
    solv1 = solve([eq3],(fi_matris))
    return solv1


mod_1 = mod(matris_1)
mod_2 = mod(matris_2)

在放置一个零向量后，我要解决的方程是：

Eq(Matrix([
[485.410196624968*a - 300.0],[185.410196624968 - 300.0*a]]),Matrix([
[0],[0]]))

所以 a 必须是 0.618 而 b 已经设置为 1.0 但给 []

顺便说一句，我编辑了其中一个函数以获得该功能：

zeros_matrice = zeros(2,1)

def mod(matris_n):
    eq1 = matris_n @ fi_matris
    eq2 = eq1.subs(b,zeros_matrice)
    solv1 = solve(eq3)
    return solv1

更准确地说，上面的 eq 在这种情况下是 eq3

我再次编辑该函数，只是为了解决一行：

def mod(matris_n):
    eq1 = matris_n[:1] @ fi_matris
    eq2 = eq1.subs(b,0)
    solv1 = solve(eq3)
    return solv1,eq3


mod_1 = mod(matris_1)
mod_2 = mod(matris_2)

在这种情况下输出是 ([],False)