1-线性代数-线性方程组(3)

伊伊_f

已于 2022-08-20 15:25:03 修改

阅读量396

点赞数

CC 4.0 BY-SA版权

分类专栏：数理课程文章标签：矩阵 python 算法

于 2022-08-04 11:28:22 首次发布

本文链接：https://blog.csdn.net/fanxianchao_2012/article/details/126155517

数理课程专栏收录该内容

5 篇文章

订阅专栏

该博客展示了如何利用Python的sympy库解决不同阶数的线性方程组。首先，它定义了符号变量和方程，然后通过`solve`函数求解二阶和四阶线性方程组。结果给出了每个未知数的表达式。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

from sympy.matrices import *
from sympy import symbols,Eq,solve

1.多阶方程组
a11,a12,a13,a14=symbols('a11 a12 a13 a14')
a21,a22,a23,a24=symbols('a21 a22 a23 a24')
a31,a32,a33,a34=symbols('a31 a32 a33 a34')
a41,a42,a43,a44=symbols('a41 a42 a43 a44')
x1,x2,x3,x4=symbols('x1 x2 x3 x4')
b1,b2,b3,b4=symbols('b1 b2 b3 b4')
f1=Eq(a11*x1+a12*x2+a13*x3+a14*x4,b1)
f2=Eq(a21*x1+a22*x2+a23*x3+a24*x4,b2)
f3=Eq(a31*x1+a32*x2+a33*x3+a34*x4,b3)
f4=Eq(a41*x1+a42*x2+a43*x3+a44*x4,b4)
solve([f1,f2,f3,f4],[x1,x2,x3,x4])

结果:

2.二阶方程组
f1=Eq(a11*x1+a12*x2,b1)
f2=Eq(a21*x1+a22*x2,b2)
solve([f1,f2],[x1,x2])

2. 求方程组解
from sympy.matrices import *
from sympy import symbols,Eq,solve
x1,x2,x3,x4=symbols('x1 x2 x3 x4')
f1=2*x1+x2-5*x3+x4-8
f2=x1-3*x2-6*x4-9
f3=2*x2-x3+2*x4+5
f4=x1+4*x2-7*x3+6*x4
fm=solve([f1,f2,f3,f4],[x1,x2,x3,x4])
fm

或者

x1,x2,x3,x4=symbols('x1 x2 x3 x4')
f1=Eq(2*x1+x2-5*x3+x4,8)
f2=Eq(x1-3*x2-6*x4,9)
f3=Eq(2*x2-x3+2*x4,-5)
f4=Eq(x1+4*x2-7*x3+6*x4,0)
fm=solve([f1,f2,f3,f4],[x1,x2,x3,x4])
fm