sympy库解多元符号方程无法跑出结果

最新推荐文章于 2024-09-05 10:28:02 发布

舶船夜雨

最新推荐文章于 2024-09-05 10:28:02 发布

阅读量228

点赞数 1

文章标签： python

本文链接：https://blog.csdn.net/jjjbbbyyyu/article/details/130252626

版权

这段代码展示了如何利用Python的Sympy库来解决一组复杂的线性代数方程。它涉及多个变量和系数，通过定义符号并应用`solve`函数找到解。

摘要由CSDN通过智能技术生成

from sympy import *

init_printing()

s, a1, a2, a3, a4, a5, a6, = symbols('s, a1, a2, a3, a4, a5, a6')

x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47 = symbols('x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47')

result=solve([(s+2*a1+a2+a5+2*a4+a3)*x1-1,

              (s+2*a1+a2+(2/3)*a3+a4)*x2-(0.5*a1+0.25*a2)*x1,

              (s+2*a1+a2+a5+(4/3)*a3+(2/3)*a4)*x3-(1.5*a1+0.75*a2)*x1,

              (s+a1+a2+a5+1.5*a3+1.5*a4)*x4-(a5+2*a4)*x1,

              (s+2*a1+a5+3*a3)*x5-a3*x1,

              (s+2*a1+a2+2*a3)*x6-(a1+a2)*x2-((2/3)*a1+(1/3)*a2)*x3,

              (s+a1+a2+1.5*a3+1.5*a4)*x7-a3*x2-0.25*(a1+a2)*x4,

              (s+2*a1+3*a3)*x8-a4*x2-0.5*a1*x5,

              (s+2*a1+a2+a5+(2/3)*a3+(1/3)*a4)*x9-((4/3)*a1+(2/3)*a2)*x3,

              (s+a1+a2+a5+a3+a4)*x10-(a5+(4/3)*a3)*x3-0.75*(a1+a2)*x4,

              (s+2*a1+a5+2*a3)*x11-(2/3)*a4*x3-1.5*a1*x5,

              (s+a2+a6+3*a4)*x12-(a5+1.5*a3)*x4,

              (s+a1+a5+3*a3)*x13-1.5*a4*x4-(a5+3*a3)*x5,

              (s+2*a1+a2+(2/3)*a3+(1/3)*a4)*x14-(2*a1+a2)*x6-(a1+0.5*a2)*x9,

              (s+2*a1+2*a2+a3+a4)*x15-(4/3)*a3*x6-(a1+a2)*x7-(1/3)*(a1+a2)*x10,

              (s+2*a1+2*a3)*x16-(2/3)*a4*x6-a1*x8-(2/3)*a1*x11,

              (s+a2+2*a4)*x17-1.5*a3*x7-(1/3)*a2*x12,

              (s+a1+3*a3)*x18-1.5*a4*x7-3*a3*x16-0.75*a1*x13,

              (s+2*a1+2*a2+a5)*x19-0.5*(2*a1+a2)*x9,

              (s+a1+a2+a5+0.5*a3+0.5*a4)*x20-(a5+(2/3)*a3)*x9-(2/3)*(a1+a2)*x10,

              (s+2*a1+a5+a3)*x21-(1/3)*a4*x9-(4/3)*a1*x11,

              (s+a2+a6+2*a4)*x22-(a5+a3)*x10-0.75*a2*x12,

              (s+a1+a5+2*a3)*x23-a3*x10-(a5+2*a3)*x11-0.75*a1*x13,

              s*x24-(2*a1+a2)*x14-a5*x13,

              s*x25-(2*a1+a2)*x14-(2*a1+a2)*x19,

              (s+a1+a2+0.5*a3+0.5*a4)*x26-(2/3)*a3*x14-(a1+a2)*x15-0.5*(a1+a2)*x20,

              (s+a3+2*a1)*x27-(1/3)*a4*x14-2*a1*x16-a1*x21,

              (s+a2+2*a4)*x28-a3*x15-a2*x17+(2/3)*a2*x22,

              (s+a1+2*a3)*x29-a4*x15-2*a3*x16-a1*x18-(1/3)*a1*x23,

              s*x30-3*a4*x17-3*a3*x18,

              (s+a1+a2+a5)*x31-(2/3)*a5*x19-0.5*(a1+a2)*x20,

              (s+2*a1+a5)*x32-a1*x21,

              (s+0.5*a2+a6+a4)*x33-(a5+0.5*a3)*x20-(2/3)*a2*x22,

              (s+a1+a5)*x34-0.5*a4*x20-(a5+a3)*x21+(2/3)*a1*x23,

              s*x35-(a6+2*a4)*x22-(a3+2*a4)*x23,

              s*x36-(a1+a2)*x26-(a1+a2)*x31,

              (s+a2+a4)*x37-0.5*a5-a2*x28-0.5*a2*x33,

              (s+a3+a4)*x38-0.5*a4*x26-a3*x27-0.5*a1*x34,

              s*x39-2*a1*x27-2*a1*x32,

              s*x40-2*a4*x28-2*a3*x29,

              (s+a2+a5)*x41-a5*x32-0.5*a2*x33,

              (s+a1+a5)*x42-a5*x32-0.5*a3*x34,

              s*x43-(a6+a4)*x33-(a5+a3)*x34,

              s*x44-a4*x37,-a4*x41,

              s*x45-a4*x37-a1*x42,

              s*x46-a1*x38-a1*x42,

              s*x47-a5*x41-a5*x42],[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47])

print(result)