标签：Python sol 0.1 key Ax np sympy numpy

numpy.linalg.solve 貌似不能解出 Ax = 0 的非零解

sympy 这个库是支持线性代数的一些运算的，但是好像又没那么成熟

所以出此下策，用 numpy 矩阵 + sympy 求解多元方程，下面以归一化条件为例求非零特解

import sympy
import numpy as np

A = np.array([[0.8, 0.1, 0.1],
              [0.7, 0.2, 0.1],
              [0.6, 0.3, 0.1]]) - np.eye(3)

dim = A.shape[0]
x = sympy.symbols([f"x{idx}" for idx in range(dim)])
# 创建未知数列表
gen_sol = sympy.solve(np.array(x) @ A, *x)
# 得到 Ax = 0 的通解
cond = sympy.Eq(sum(x), 1)
# 创建归一化条件
equals = [sympy.Eq(key, value) for key, value in gen_sol.items()] + [cond]
# 以通解创建方程组, 附加条件
part_sol = sympy.solve(equals, x)
# 得到特定条件下的特解
assert part_sol, "Ax = 0 在该条件下无解"
result = np.array([part_sol[key] for key in x])
print(result)

标签：Python,sol,0.1,key,Ax,np,sympy,numpy
来源： https://blog.csdn.net/qq_55745968/article/details/123094985

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。