sympy解符号方程

于 2024-07-08 11:41:10 发布
阅读量1k 收藏 10
点赞数 21
分类专栏： Python 笔记 sympy 文章标签： python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36980284/article/details/140264237
版权
笔记 同时被 3 个专栏收录
90 篇文章 4 订阅
订阅专栏
Python
19 篇文章 0 订阅
订阅专栏
sympy
1 篇文章 0 订阅
订阅专栏 
y1=np.arange(0.1,1,0.1)#竖直方向的点,包气带
y2=np.arange(-4,0,0.1)#竖直方向的点,饱和带
for i,j in zip(y1,y2):
    print(i,j)

0.1 -4.0
0.2 -3.9
0.30000000000000004 -3.8
0.4 -3.6999999999999997
0.5 -3.5999999999999996
0.6 -3.4999999999999996
0.7000000000000001 -3.3999999999999995
0.8 -3.2999999999999994
0.9 -3.1999999999999993

from sympy import *
x = symbols('x')
a = Integral(cos(x)*exp(x), x)
Eq(a, a.doit())

∫ e x cos ⁡ ( x )   d x = e x sin ⁡ ( x ) 2 + e x cos ⁡ ( x ) 2 \displaystyle \int e^{x} \cos{\left(x \right)}\, dx = \frac{e^{x} \sin{\left(x \right)}}{2} + \frac{e^{x} \cos{\left(x \right)}}{2} ∫excos(x)dx=2exsin(x)​+2excos(x)​

import numpy as np
infiltration=np.random.randn(100)
np.savetxt(r'C:\Users\Administrator\Desktop\infiltrate.txt',infiltration)

解方程组

from sympy import *

# 将变量符号化
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
# 解三元一次方程
f1 = x+y+z-2
f2 = 2*x-y+z+1
f3 = x+2*y+2*z-3
r3 = solve([f1, f2, f3], [x, y, z])
print("r3:", r3)

r3: {x: 1, y: 2, z: -1}

"""
    x+y+z-2=0
    2x-y+z+1=0
    x+2y+2z-3=0
"""
from sympy import *


x, y, z = symbols("x y z")

# 默认等式为0的形式
print("======默认等式为0的形式 =======")
eq = [x+y+z-2, 2*x-y+z+1, x+2*y+2*z-3]
result = linsolve(eq, [x, y, z])
print(result)
print(latex(result))

# 矩阵形式
print("======矩阵形式 =======")
eq = Matrix(([1, 1, 1, 2], [2, -1, 1, -1], [1, 2, 2, 3]))
result = linsolve(eq, [x, y, z])
print(result)
print(latex(result))

# 增广矩阵形式
print("======增广矩阵形式 =======")
A = Matrix([[1, 1, 1], [2, -1, 1], [1, 2, 2]])
b = Matrix([[2], [-1], [3]])
system = A, b
result = linsolve(system, x, y, z)
# print(result)
# print(latex(result))
result

======默认等式为0的形式 =======
{(1, 2, -1)}
\left\{\left( 1, \  2, \  -1\right)\right\}
======矩阵形式 =======
{(1, 2, -1)}
\left\{\left( 1, \  2, \  -1\right)\right\}
======增广矩阵形式 =======

{ ( 1 ,   2 ,   − 1 ) } \displaystyle \left\{\left( 1, \ 2, \ -1\right)\right\} {(1, 2, −1)}

from sympy import *
V, ak, A, B,C,D,Aa,Ba,Ca,Da,E,Δ,Λ,λ,Φ,c,S,F \
= symbols('V a_k A B C D Aa Ba Ca Da E Δ Λ λ Φ c S F')
C1,C2,C3,C4=symbols('C1 C2 C3 C4')
eq1=Eq(Λ*E*((2*ak*A-λ*Φ*B)*C1+(2*ak*C-λ*Φ*D)*C2)-C4)
eq1

C:\ProgramData\Anaconda3\lib\site-packages\sympy\core\relational.py:495: SymPyDeprecationWarning: 

Eq(expr) with rhs default to 0 has been deprecated since SymPy 1.5.
Use Eq(expr, 0) instead. See
https://github.com/sympy/sympy/issues/16587 for more info.

  SymPyDeprecationWarning(

− C 4 + E Λ ( C 1 ( 2 A a k − B Φ λ ) + C 2 ( 2 C a k − D Φ λ ) ) = 0 \displaystyle - C_{4} + E Λ \left(C_{1} \left(2 A a_{k} - B Φ λ\right) + C_{2} \left(2 C a_{k} - D Φ λ\right)\right) = 0 −C4​+EΛ(C1​(2Aak​−BΦλ)+C2​(2Cak​−DΦλ))=0

from sympy import *
V, ak, A, B,C,D,Aa,Ba,Ca,Da,E,Δ,Λ,λ,Φ,c,S,F \
= symbols("V a_k \\bar{A} \\bar{B} \\bar{C} \\bar{D} A' B' C' D' E Δ Λ λ Φ c S F")
C1,C2,C3,C4=symbols('C1 C2 C3 C4')
result = solve([Λ*E*((2*ak*A-λ*Φ*B)*C1+(2*ak*C-λ*Φ*D)*C2)-C4,
                Δ*(C1*(A)+C2*(C))-C4,
                V*((2*ak*Aa-S*Ba)*C1+(2*ak*Ca-S*Da)*C2)-F],[C1,C2,C4])
result

{ C 1 : 2 E F C ˉ a k Λ − E F D ˉ Λ Φ λ − F C ˉ Δ 4 A ′ E V C ˉ a k 2 Λ − 2 A ′ E V D ˉ a k Λ Φ λ − 2 A ′ V C ˉ a k Δ − 2 B ′ E S V C ˉ a k Λ + B ′ E S V D ˉ Λ Φ λ + B ′ S V C ˉ Δ − 4 C ′ E V A ˉ a k 2 Λ + 2 C ′ E V B ˉ a k Λ Φ λ + 2 C ′ V A ˉ a k Δ + 2 D ′ E S V A ˉ a k Λ − D ′ E S V B ˉ Λ Φ λ − D ′ S V A ˉ Δ ,   C 2 : − 2 E F A ˉ a k Λ + E F B ˉ Λ Φ λ + F A ˉ Δ 4 A ′ E V C ˉ a k 2 Λ − 2 A ′ E V D ˉ a k Λ Φ λ − 2 A ′ V C ˉ a k Δ − 2 B ′ E S V C ˉ a k Λ + B ′ E S V D ˉ Λ Φ λ + B ′ S V C ˉ Δ − 4 C ′ E V A ˉ a k 2 Λ + 2 C ′ E V B ˉ a k Λ Φ λ + 2 C ′ V A ˉ a k Δ + 2 D ′ E S V A ˉ a k Λ − D ′ E S V B ˉ Λ Φ λ − D ′ S V A ˉ Δ ,   C 4 : − E F A ˉ D ˉ Δ Λ Φ λ + E F B ˉ C ˉ Δ Λ Φ λ 4 A ′ E V C ˉ a k 2 Λ − 2 A ′ E V D ˉ a k Λ Φ λ − 2 A ′ V C ˉ a k Δ − 2 B ′ E S V C ˉ a k Λ + B ′ E S V D ˉ Λ Φ λ + B ′ S V C ˉ Δ − 4 C ′ E V A ˉ a k 2 Λ + 2 C ′ E V B ˉ a k Λ Φ λ + 2 C ′ V A ˉ a k Δ + 2 D ′ E S V A ˉ a k Λ − D ′ E S V B ˉ Λ Φ λ − D ′ S V A ˉ Δ } \displaystyle \left\{ C_{1} : \frac{2 E F \bar{C} a_{k} Λ - E F \bar{D} Λ Φ λ - F \bar{C} Δ}{4 A' E V \bar{C} a_{k}^{2} Λ - 2 A' E V \bar{D} a_{k} Λ Φ λ - 2 A' V \bar{C} a_{k} Δ - 2 B' E S V \bar{C} a_{k} Λ + B' E S V \bar{D} Λ Φ λ + B' S V \bar{C} Δ - 4 C' E V \bar{A} a_{k}^{2} Λ + 2 C' E V \bar{B} a_{k} Λ Φ λ + 2 C' V \bar{A} a_{k} Δ + 2 D' E S V \bar{A} a_{k} Λ - D' E S V \bar{B} Λ Φ λ - D' S V \bar{A} Δ}, \ C_{2} : \frac{- 2 E F \bar{A} a_{k} Λ + E F \bar{B} Λ Φ λ + F \bar{A} Δ}{4 A' E V \bar{C} a_{k}^{2} Λ - 2 A' E V \bar{D} a_{k} Λ Φ λ - 2 A' V \bar{C} a_{k} Δ - 2 B' E S V \bar{C} a_{k} Λ + B' E S V \bar{D} Λ Φ λ + B' S V \bar{C} Δ - 4 C' E V \bar{A} a_{k}^{2} Λ + 2 C' E V \bar{B} a_{k} Λ Φ λ + 2 C' V \bar{A} a_{k} Δ + 2 D' E S V \bar{A} a_{k} Λ - D' E S V \bar{B} Λ Φ λ - D' S V \bar{A} Δ}, \ C_{4} : \frac{- E F \bar{A} \bar{D} Δ Λ Φ λ + E F \bar{B} \bar{C} Δ Λ Φ λ}{4 A' E V \bar{C} a_{k}^{2} Λ - 2 A' E V \bar{D} a_{k} Λ Φ λ - 2 A' V \bar{C} a_{k} Δ - 2 B' E S V \bar{C} a_{k} Λ + B' E S V \bar{D} Λ Φ λ + B' S V \bar{C} Δ - 4 C' E V \bar{A} a_{k}^{2} Λ + 2 C' E V \bar{B} a_{k} Λ Φ λ + 2 C' V \bar{A} a_{k} Δ + 2 D' E S V \bar{A} a_{k} Λ - D' E S V \bar{B} Λ Φ λ - D' S V \bar{A} Δ}\right\} {C1​:4A′EVCˉak2​Λ−2A′EVDˉak​ΛΦλ−2A′VCˉak​Δ−2B′ESVCˉak​Λ+B′ESVDˉΛΦλ+B′SVCˉΔ−4C′EVAˉak2​Λ+2C′EVBˉak​ΛΦλ+2C′VAˉak​Δ+2D′ESVAˉak​Λ−D′ESVBˉΛΦλ−D′SVAˉΔ2EFCˉak​Λ−EFDˉΛΦλ−FCˉΔ​, C2​:4A′EVCˉak2​Λ−2A′EVDˉak​ΛΦλ−2A′VCˉak​Δ−2B′ESVCˉak​Λ+B′ESVDˉΛΦλ+B′SVCˉΔ−4C′EVAˉak2​Λ+2C′EVBˉak​ΛΦλ+2C′VAˉak​Δ+2D′ESVAˉak​Λ−D′ESVBˉΛΦλ−D′SVAˉΔ−2EFAˉak​Λ+EFBˉΛΦλ+FAˉΔ​, C4​:4A′EVCˉak2​Λ−2A′EVDˉak​ΛΦλ−2A′VCˉak​Δ−2B′ESVCˉak​Λ+B′ESVDˉΛΦλ+B′SVCˉΔ−4C′EVAˉak2​Λ+2C′EVBˉak​ΛΦλ+2C′VAˉak​Δ+2D′ESVAˉak​Λ−D′ESVBˉΛΦλ−D′SVAˉΔ−EFAˉDˉΔΛΦλ+EFBˉCˉΔΛΦλ​}

result2=linsolve([Λ*E*((2*ak*A-λ*Φ*B)*C1+(2*ak*C-λ*Φ*D)*C2)-C4,
                Δ*(C1*(A)+C2*(C))-C4,
                V*((2*ak*Aa-S*Ba)*C1+(2*ak*Ca-S*Da)*C2)-F],[C1,C2,C4])
result2

{ ( − 2 C E F a k Λ + C F Δ + D E F Λ Φ λ 4 A C ′ E V a k 2 Λ − 2 A C ′ V a k Δ − 2 A D ′ E S V a k Λ + A D ′ S V Δ − 4 A ′ C E V a k 2 Λ + 2 A ′ C V a k Δ + 2 A ′ D E V a k Λ Φ λ − 2 B C ′ E V a k Λ Φ λ + B D ′ E S V Λ Φ λ + 2 B ′ C E S V a k Λ − B ′ C S V Δ − B ′ D E S V Λ Φ λ ,   2 A E F a k Λ − A F Δ − B E F Λ Φ λ 4 A C ′ E V a k 2 Λ − 2 A C ′ V a k Δ − 2 A D ′ E S V a k Λ + A D ′ S V Δ − 4 A ′ C E V a k 2 Λ + 2 A ′ C V a k Δ + 2 A ′ D E V a k Λ Φ λ − 2 B C ′ E V a k Λ Φ λ + B D ′ E S V Λ Φ λ + 2 B ′ C E S V a k Λ − B ′ C S V Δ − B ′ D E S V Λ Φ λ ,   A D E F Δ Λ Φ λ − B C E F Δ Λ Φ λ 4 A C ′ E V a k 2 Λ − 2 A C ′ V a k Δ − 2 A D ′ E S V a k Λ + A D ′ S V Δ − 4 A ′ C E V a k 2 Λ + 2 A ′ C V a k Δ + 2 A ′ D E V a k Λ Φ λ − 2 B C ′ E V a k Λ Φ λ + B D ′ E S V Λ Φ λ + 2 B ′ C E S V a k Λ − B ′ C S V Δ − B ′ D E S V Λ Φ λ ) } \displaystyle \left\{\left( \frac{- 2 C E F a_{k} Λ + C F Δ + D E F Λ Φ λ}{4 A C' E V a_{k}^{2} Λ - 2 A C' V a_{k} Δ - 2 A D' E S V a_{k} Λ + A D' S V Δ - 4 A' C E V a_{k}^{2} Λ + 2 A' C V a_{k} Δ + 2 A' D E V a_{k} Λ Φ λ - 2 B C' E V a_{k} Λ Φ λ + B D' E S V Λ Φ λ + 2 B' C E S V a_{k} Λ - B' C S V Δ - B' D E S V Λ Φ λ}, \ \frac{2 A E F a_{k} Λ - A F Δ - B E F Λ Φ λ}{4 A C' E V a_{k}^{2} Λ - 2 A C' V a_{k} Δ - 2 A D' E S V a_{k} Λ + A D' S V Δ - 4 A' C E V a_{k}^{2} Λ + 2 A' C V a_{k} Δ + 2 A' D E V a_{k} Λ Φ λ - 2 B C' E V a_{k} Λ Φ λ + B D' E S V Λ Φ λ + 2 B' C E S V a_{k} Λ - B' C S V Δ - B' D E S V Λ Φ λ}, \ \frac{A D E F Δ Λ Φ λ - B C E F Δ Λ Φ λ}{4 A C' E V a_{k}^{2} Λ - 2 A C' V a_{k} Δ - 2 A D' E S V a_{k} Λ + A D' S V Δ - 4 A' C E V a_{k}^{2} Λ + 2 A' C V a_{k} Δ + 2 A' D E V a_{k} Λ Φ λ - 2 B C' E V a_{k} Λ Φ λ + B D' E S V Λ Φ λ + 2 B' C E S V a_{k} Λ - B' C S V Δ - B' D E S V Λ Φ λ}\right)\right\} {(4AC′EVak2​Λ−2AC′Vak​Δ−2AD′ESVak​Λ+AD′SVΔ−4A′CEVak2​Λ+2A′CVak​Δ+2A′DEVak​ΛΦλ−2BC′EVak​ΛΦλ+BD′ESVΛΦλ+2B′CESVak​Λ−B′CSVΔ−B′DESVΛΦλ−2CEFak​Λ+CFΔ+DEFΛΦλ​, 4AC′EVak2​Λ−2AC′Vak​Δ−2AD′ESVak​Λ+AD′SVΔ−4A′CEVak2​Λ+2A′CVak​Δ+2A′DEVak​ΛΦλ−2BC′EVak​ΛΦλ+BD′ESVΛΦλ+2B′CESVak​Λ−B′CSVΔ−B′DESVΛΦλ2AEFak​Λ−AFΔ−BEFΛΦλ​, 4AC′EVak2​Λ−2AC′Vak​Δ−2AD′ESVak​Λ+AD′SVΔ−4A′CEVak2​Λ+2A′CVak​Δ+2A′DEVak​ΛΦλ−2BC′EVak​ΛΦλ+BD′ESVΛΦλ+2B′CESVak​Λ−B′CSVΔ−B′DESVΛΦλADEFΔΛΦλ−BCEFΔΛΦλ​)}

from sympy import *
V, ak, A, B,C,D,Aa,Ba,Ca,Da,E,Δ,Λ,λ,Φ,c,S,F \
= symbols("V a_k \\bar{A} \\bar{B} \\bar{C} \\bar{D} A' B' C' D' E Δ Λ λ Φ c S F")
C1,C2,C3,C4=symbols('C1 C2 C3 C4')
result = solve([Λ*E*((2*ak*A-λ*Φ*B)*C1+(2*ak*C-λ*Φ*D)*C2)-C4,
                Δ*(C1*(A)+C2*(C))-C4,
                V*((2*ak*Aa-S*Ba)*C1+(2*ak*Ca-S*Da)*C2)-F],[C1,C2,C4])
print(result)

{C1: (2*E*F*\bar{C}*a_k*Λ - E*F*\bar{D}*Λ*Φ*λ - F*\bar{C}*Δ)/(4*A'*E*V*\bar{C}*a_k**2*Λ - 2*A'*E*V*\bar{D}*a_k*Λ*Φ*λ - 2*A'*V*\bar{C}*a_k*Δ - 2*B'*E*S*V*\bar{C}*a_k*Λ + B'*E*S*V*\bar{D}*Λ*Φ*λ + B'*S*V*\bar{C}*Δ - 4*C'*E*V*\bar{A}*a_k**2*Λ + 2*C'*E*V*\bar{B}*a_k*Λ*Φ*λ + 2*C'*V*\bar{A}*a_k*Δ + 2*D'*E*S*V*\bar{A}*a_k*Λ - D'*E*S*V*\bar{B}*Λ*Φ*λ - D'*S*V*\bar{A}*Δ), C2: (-2*E*F*\bar{A}*a_k*Λ + E*F*\bar{B}*Λ*Φ*λ + F*\bar{A}*Δ)/(4*A'*E*V*\bar{C}*a_k**2*Λ - 2*A'*E*V*\bar{D}*a_k*Λ*Φ*λ - 2*A'*V*\bar{C}*a_k*Δ - 2*B'*E*S*V*\bar{C}*a_k*Λ + B'*E*S*V*\bar{D}*Λ*Φ*λ + B'*S*V*\bar{C}*Δ - 4*C'*E*V*\bar{A}*a_k**2*Λ + 2*C'*E*V*\bar{B}*a_k*Λ*Φ*λ + 2*C'*V*\bar{A}*a_k*Δ + 2*D'*E*S*V*\bar{A}*a_k*Λ - D'*E*S*V*\bar{B}*Λ*Φ*λ - D'*S*V*\bar{A}*Δ), C4: (-E*F*\bar{A}*\bar{D}*Δ*Λ*Φ*λ + E*F*\bar{B}*\bar{C}*Δ*Λ*Φ*λ)/(4*A'*E*V*\bar{C}*a_k**2*Λ - 2*A'*E*V*\bar{D}*a_k*Λ*Φ*λ - 2*A'*V*\bar{C}*a_k*Δ - 2*B'*E*S*V*\bar{C}*a_k*Λ + B'*E*S*V*\bar{D}*Λ*Φ*λ + B'*S*V*\bar{C}*Δ - 4*C'*E*V*\bar{A}*a_k**2*Λ + 2*C'*E*V*\bar{B}*a_k*Λ*Φ*λ + 2*C'*V*\bar{A}*a_k*Δ + 2*D'*E*S*V*\bar{A}*a_k*Λ - D'*E*S*V*\bar{B}*Λ*Φ*λ - D'*S*V*\bar{A}*Δ)}

from sympy import collect
from sympy.abc import *
from sympy import symbols
a_k,γ,δ,ε,η,ζ,α,β,θ=symbols('a_k γ δ ε η ζ α β θ')
# V, ak, A, B,C,D,Aa,Ba,Ca,Da,E,Δ,Λ,λ,Φ,c,S,F \
# = symbols("V a_k \\bar{A} \\bar{B} \\bar{C} \\bar{D} A' B' C' D' E Δ Λ λ Φ c S F")
eq1=(2*E*F*γ*a_k*Λ - E*F*δ*Λ*Φ*λ - F*γ*Δ)/(4*ε*E*V*γ*a_k**2*Λ - 2*ε*E*V*δ*a_k*Λ*Φ*λ - 2*ε*V*γ*a_k*Δ - 2*ζ*E*S*V*γ*a_k*Λ + ζ*E*S*V*δ*Λ*Φ*λ + ζ*S*V*γ*Δ - 4*η*E*V*α*a_k**2*Λ + 2*η*E*V*β*a_k*Λ*Φ*λ + 2*η*V*α*a_k*Δ + 2*θ*E*S*V*α*a_k*Λ - θ*E*S*V*β*Λ*Φ*λ - θ*S*V*α*Δ)
collect(eq1,[F,V,E])

F ( 2 E a k Λ γ − E Λ Φ δ λ − Δ γ ) V ( 2 E S a k Λ α θ − 2 E S a k Λ γ ζ − E S Λ Φ β θ λ + E S Λ Φ δ ζ λ − 4 E a k 2 Λ α η + 4 E a k 2 Λ γ ε + 2 E a k Λ Φ β η λ − 2 E a k Λ Φ δ ε λ − S Δ α θ + S Δ γ ζ + 2 a k Δ α η − 2 a k Δ γ ε ) \displaystyle \frac{F \left(2 E a_{k} Λ γ - E Λ Φ δ λ - Δ γ\right)}{V \left(2 E S a_{k} Λ α θ - 2 E S a_{k} Λ γ ζ - E S Λ Φ β θ λ + E S Λ Φ δ ζ λ - 4 E a_{k}^{2} Λ α η + 4 E a_{k}^{2} Λ γ ε + 2 E a_{k} Λ Φ β η λ - 2 E a_{k} Λ Φ δ ε λ - S Δ α θ + S Δ γ ζ + 2 a_{k} Δ α η - 2 a_{k} Δ γ ε\right)} V(2ESak​Λαθ−2ESak​Λγζ−ESΛΦβθλ+ESΛΦδζλ−4Eak2​Λαη+4Eak2​Λγε+2Eak​ΛΦβηλ−2Eak​ΛΦδελ−SΔαθ+SΔγζ+2ak​Δαη−2ak​Δγε)F(2Eak​Λγ−EΛΦδλ−Δγ)​

from sympy import *
init_printing()
c1, c2, c3, c4, c5, c6, D1, D2, D3, D4 = symbols('c1, c2, c3, c4, c5, c6, D1, D2, D3, D4')
x1, x2, x3, x4 = symbols('x1 x2 x3 x4 ')
 
result = solve([c1*x1 + c2*x2 - c5*x4 - D1,
           c4*x1 - c1*x2 - c5*x3 - D2,
           c3*x2 + c1*x3 + c6*x4 - D3,
           c3*x1 + c6*x3 - c1*x4 - D4],[x1, x2, x3, x4])
result

C:\ProgramData\Anaconda3\lib\site-packages\IPython\lib\latextools.py:126: MatplotlibDeprecationWarning: 
The to_png function was deprecated in Matplotlib 3.4 and will be removed two minor releases later. Use mathtext.math_to_image instead.
  mt.to_png(f, s, fontsize=12, dpi=dpi, color=color)
C:\ProgramData\Anaconda3\lib\site-packages\IPython\lib\latextools.py:126: MatplotlibDeprecationWarning: 
The to_rgba function was deprecated in Matplotlib 3.4 and will be removed two minor releases later. Use mathtext.math_to_image instead.
  mt.to_png(f, s, fontsize=12, dpi=dpi, color=color)
C:\ProgramData\Anaconda3\lib\site-packages\IPython\lib\latextools.py:126: MatplotlibDeprecationWarning: 
The to_mask function was deprecated in Matplotlib 3.4 and will be removed two minor releases later. Use mathtext.math_to_image instead.
  mt.to_png(f, s, fontsize=12, dpi=dpi, color=color)
C:\ProgramData\Anaconda3\lib\site-packages\IPython\lib\latextools.py:126: MatplotlibDeprecationWarning: 
The MathtextBackendBitmap class was deprecated in Matplotlib 3.4 and will be removed two minor releases later. Use mathtext.math_to_image instead.
  mt.to_png(f, s, fontsize=12, dpi=dpi, color=color)

{ x 1 : D 1 c 1 3 − D 1 c 1 c 3 c 5 + D 1 c 1 c 6 2 + D 2 c 1 2 c 2 + D 2 c 2 c 6 2 + D 2 c 3 c 5 c 6 + D 3 c 1 c 2 c 5 + D 3 c 1 c 5 c 6 − D 4 c 1 2 c 5 + D 4 c 2 c 5 c 6 + D 4 c 3 c 5 2 c 1 4 + c 1 2 c 2 c 4 − 2 c 1 2 c 3 c 5 + c 1 2 c 6 2 + c 2 c 3 c 5 c 6 + c 2 c 4 c 6 2 + c 3 2 c 5 2 + c 3 c 4 c 5 c 6 ,   x 2 : D 1 c 1 2 c 4 + D 1 c 3 c 5 c 6 + D 1 c 4 c 6 2 − D 2 c 1 3 + D 2 c 1 c 3 c 5 − D 2 c 1 c 6 2 − D 3 c 1 2 c 5 + D 3 c 3 c 5 2 + D 3 c 4 c 5 c 6 − D 4 c 1 c 4 c 5 − D 4 c 1 c 5 c 6 c 1 4 + c 1 2 c 2 c 4 − 2 c 1 2 c 3 c 5 + c 1 2 c 6 2 + c 2 c 3 c 5 c 6 + c 2 c 4 c 6 2 + c 3 2 c 5 2 + c 3 c 4 c 5 c 6 ,   x 3 : − D 1 c 1 c 3 c 4 − D 1 c 1 c 3 c 6 + D 2 c 1 2 c 3 − D 2 c 2 c 3 c 6 − D 2 c 3 2 c 5 + D 3 c 1 3 + D 3 c 1 c 2 c 4 − D 3 c 1 c 3 c 5 + D 4 c 1 2 c 6 + D 4 c 2 c 4 c 6 + D 4 c 3 c 4 c 5 c 1 4 + c 1 2 c 2 c 4 − 2 c 1 2 c 3 c 5 + c 1 2 c 6 2 + c 2 c 3 c 5 c 6 + c 2 c 4 c 6 2 + c 3 2 c 5 2 + c 3 c 4 c 5 c 6 ,   x 4 : D 1 c 1 2 c 3 − D 1 c 3 2 c 5 − D 1 c 3 c 4 c 6 + D 2 c 1 c 2 c 3 + D 2 c 1 c 3 c 6 + D 3 c 1 2 c 6 + D 3 c 2 c 3 c 5 + D 3 c 2 c 4 c 6 − D 4 c 1 3 − D 4 c 1 c 2 c 4 + D 4 c 1 c 3 c 5 c 1 4 + c 1 2 c 2 c 4 − 2 c 1 2 c 3 c 5 + c 1 2 c 6 2 + c 2 c 3 c 5 c 6 + c 2 c 4 c 6 2 + c 3 2 c 5 2 + c 3 c 4 c 5 c 6 } \displaystyle \left\{ x_{1} : \frac{D_{1} c_{1}^{3} - D_{1} c_{1} c_{3} c_{5} + D_{1} c_{1} c_{6}^{2} + D_{2} c_{1}^{2} c_{2} + D_{2} c_{2} c_{6}^{2} + D_{2} c_{3} c_{5} c_{6} + D_{3} c_{1} c_{2} c_{5} + D_{3} c_{1} c_{5} c_{6} - D_{4} c_{1}^{2} c_{5} + D_{4} c_{2} c_{5} c_{6} + D_{4} c_{3} c_{5}^{2}}{c_{1}^{4} + c_{1}^{2} c_{2} c_{4} - 2 c_{1}^{2} c_{3} c_{5} + c_{1}^{2} c_{6}^{2} + c_{2} c_{3} c_{5} c_{6} + c_{2} c_{4} c_{6}^{2} + c_{3}^{2} c_{5}^{2} + c_{3} c_{4} c_{5} c_{6}}, \ x_{2} : \frac{D_{1} c_{1}^{2} c_{4} + D_{1} c_{3} c_{5} c_{6} + D_{1} c_{4} c_{6}^{2} - D_{2} c_{1}^{3} + D_{2} c_{1} c_{3} c_{5} - D_{2} c_{1} c_{6}^{2} - D_{3} c_{1}^{2} c_{5} + D_{3} c_{3} c_{5}^{2} + D_{3} c_{4} c_{5} c_{6} - D_{4} c_{1} c_{4} c_{5} - D_{4} c_{1} c_{5} c_{6}}{c_{1}^{4} + c_{1}^{2} c_{2} c_{4} - 2 c_{1}^{2} c_{3} c_{5} + c_{1}^{2} c_{6}^{2} + c_{2} c_{3} c_{5} c_{6} + c_{2} c_{4} c_{6}^{2} + c_{3}^{2} c_{5}^{2} + c_{3} c_{4} c_{5} c_{6}}, \ x_{3} : \frac{- D_{1} c_{1} c_{3} c_{4} - D_{1} c_{1} c_{3} c_{6} + D_{2} c_{1}^{2} c_{3} - D_{2} c_{2} c_{3} c_{6} - D_{2} c_{3}^{2} c_{5} + D_{3} c_{1}^{3} + D_{3} c_{1} c_{2} c_{4} - D_{3} c_{1} c_{3} c_{5} + D_{4} c_{1}^{2} c_{6} + D_{4} c_{2} c_{4} c_{6} + D_{4} c_{3} c_{4} c_{5}}{c_{1}^{4} + c_{1}^{2} c_{2} c_{4} - 2 c_{1}^{2} c_{3} c_{5} + c_{1}^{2} c_{6}^{2} + c_{2} c_{3} c_{5} c_{6} + c_{2} c_{4} c_{6}^{2} + c_{3}^{2} c_{5}^{2} + c_{3} c_{4} c_{5} c_{6}}, \ x_{4} : \frac{D_{1} c_{1}^{2} c_{3} - D_{1} c_{3}^{2} c_{5} - D_{1} c_{3} c_{4} c_{6} + D_{2} c_{1} c_{2} c_{3} + D_{2} c_{1} c_{3} c_{6} + D_{3} c_{1}^{2} c_{6} + D_{3} c_{2} c_{3} c_{5} + D_{3} c_{2} c_{4} c_{6} - D_{4} c_{1}^{3} - D_{4} c_{1} c_{2} c_{4} + D_{4} c_{1} c_{3} c_{5}}{c_{1}^{4} + c_{1}^{2} c_{2} c_{4} - 2 c_{1}^{2} c_{3} c_{5} + c_{1}^{2} c_{6}^{2} + c_{2} c_{3} c_{5} c_{6} + c_{2} c_{4} c_{6}^{2} + c_{3}^{2} c_{5}^{2} + c_{3} c_{4} c_{5} c_{6}}\right\} {x1​:c14​+c12​c2​c4​−2c12​c3​c5​+c12​c62​+c2​c3​c5​c6​+c2​c4​c62​+c32​c52​+c3​c4​c5​c6​D1​c13​−D1​c1​c3​c5​+D1​c1​c62​+D2​c12​c2​+D2​c2​c62​+D2​c3​c5​c6​+D3​c1​c2​c5​+D3​c1​c5​c6​−D4​c12​c5​+D4​c2​c5​c6​+D4​c3​c52​​, x2​:c14​+c12​c2​c4​−2c12​c3​c5​+c12​c62​+c2​c3​c5​c6​+c2​c4​c62​+c32​c52​+c3​c4​c5​c6​D1​c12​c4​+D1​c3​c5​c6​+D1​c4​c62​−D2​c13​+D2​c1​c3​c5​−D2​c1​c62​−D3​c12​c5​+D3​c3​c52​+D3​c4​c5​c6​−D4​c1​c4​c5​−D4​c1​c5​c6​​, x3​:c14​+c12​c2​c4​−2c12​c3​c5​+c12​c62​+c2​c3​c5​c6​+c2​c4​c62​+c32​c52​+c3​c4​c5​c6​−D1​c1​c3​c4​−D1​c1​c3​c6​+D2​c12​c3​−D2​c2​c3​c6​−D2​c32​c5​+D3​c13​+D3​c1​c2​c4​−D3​c1​c3​c5​+D4​c12​c6​+D4​c2​c4​c6​+D4​c3​c4​c5​​, x4​:c14​+c12​c2​c4​−2c12​c3​c5​+c12​c62​+c2​c3​c5​c6​+c2​c4​c62​+c32​c52​+c3​c4​c5​c6​D1​c12​c3​−D1​c32​c5​−D1​c3​c4​c6​+D2​c1​c2​c3​+D2​c1​c3​c6​+D3​c12​c6​+D3​c2​c3​c5​+D3​c2​c4​c6​−D4​c13​−D4​c1​c2​c4​+D4​c1​c3​c5​​}

from sympy.abc import *
expr = x*y + x - 3 + 2*x**2 - z*x**2 + x**3
expr
x**3 - x**2*z + 2*x**2 + x*y + x - 3
collected_expr = collect(expr, x)
collected_expr
x**3 + x**2*(-z + 2) + x*(y + 1) - 3

x 3 + x 2 ( 2 − z ) + x ( y + 1 ) − 3 \displaystyle x^{3} + x^{2} \left(2 - z\right) + x \left(y + 1\right) - 3 x3+x2(2−z)+x(y+1)−3

# -- coding:utf-8 --
from sympy import *
# init_printing()
#用sympy符号运算解方程
t=symbols('t',real = True) # real 保证全是实数，自变量
# h=symbols('h',function = True) # 全部为函数变量
# r=symbols('r',function = True) # 全部为函数变量
h = Function('h')
r = Function('r')
S,a= symbols('S a')
ode=S*h(t).diff(t)+a*h(t)-r(t)
eq=Eq(ode)
eq

S d d t h ( t ) + a h ( t ) − r ( t ) = 0 \displaystyle S \frac{d}{d t} h{\left(t \right)} + a h{\left(t \right)} - r{\left(t \right)} = 0 Sdtd​h(t)+ah(t)−r(t)=0

# eq
dsolve(eq,h(t))

∫ ( a h ( t ) − r ( t ) ) e a t S   d t = C 1 \displaystyle \int \left(a h{\left(t \right)} - r{\left(t \right)}\right) e^{\frac{a t}{S}}\, dt = C_{1} ∫(ah(t)−r(t))eSat​dt=C1​

from sympy import *
f = symbols('f', cls=Function)
x = symbols('x')
eq = Eq(f(x).diff(x,2) - 2*f(x).diff(x,1) + f(x), sin(x))
eq
# print(dsolve(eq, f(x)))

f ( x ) − 2 d d x f ( x ) + d 2 d x 2 f ( x ) = sin ⁡ ( x ) \displaystyle f{\left(x \right)} - 2 \frac{d}{d x} f{\left(x \right)} + \frac{d^{2}}{d x^{2}} f{\left(x \right)} = \sin{\left(x \right)} f(x)−2dxd​f(x)+dx2d2​f(x)=sin(x)

使用scipy求解数值解

import numpy as np
from scipy import integrate
import matplotlib.pyplot as plt
import sympy

def plot_direction_field(x, y_x, f_xy, x_lim=(-5, 5), y_lim=(-5, 5), ax=None):
    f_np = sympy.lambdify((x, y_x), f_xy, 'numpy')
    x_vec = np.linspace(x_lim[0], x_lim[1], 20)
    y_vec = np.linspace(y_lim[0], y_lim[1], 20)

    if ax is None:
        _, ax = plt.subplots(figsize=(4, 4))

    dx = x_vec[1] - x_vec[0]
    dy = y_vec[1] - y_vec[0]

    for m, xx in enumerate(x_vec):
        for n, yy in enumerate(y_vec):
            Dy = f_np(xx, yy) * dx
            Dx = 0.8 * dx**2 / np.sqrt(dx**2 + Dy**2)
            Dy = 0.8 * Dy*dy / np.sqrt(dx**2 + Dy**2)
            ax.plot([xx - Dx/2, xx + Dx/2], [yy - Dy/2, yy + Dy/2], 'b', lw=0.5)

    ax.axis('tight')
    ax.set_title(r"$%s$" %(sympy.latex(sympy.Eq(y_x.diff(x), f_xy))), fontsize=18)

    return ax

x = sympy.symbols('x')
y = sympy.Function('y')
f = x-y(x)**2

f_np = sympy.lambdify((y(x), x), f)
## put variables (y(x), x) into lambda function f.
y0 = 1
xp = np.linspace(0, 5, 100)
yp = integrate.odeint(f_np, y0, xp)
## solve f_np with initial conditons y0, and x ranges as xp.
xn = np.linspace(0, -5, 100)
yn = integrate.odeint(f_np, y0, xn)

fig, ax = plt.subplots(1, 1, figsize=(4, 4))
plot_direction_field(x, y(x), f, ax=ax)
## plot direction field of function f
ax.plot(xn, yn, 'b', lw=2)
ax.plot(xp, yp, 'r', lw=2)
plt.show()

C:\ProgramData\Anaconda3\lib\site-packages\scipy\integrate\odepack.py:247: ODEintWarning: Excess work done on this call (perhaps wrong Dfun type). Run with full_output = 1 to get quantitative information.
  warnings.warn(warning_msg, ODEintWarning)

在这里插入图片描述

小孟的CDN
关注
专栏目录
sympy解符号方程(2).ipynb
07-08
sympy解符号方程(2)——公式的标准格式显示，及求解结果，符号化求解！
「SymPy」符号运算(3) (非)线性方程(组)求解、数列求和、连乘、求极限
weixin_45549554的博客
01-17 1728
SymPy求解非线性方程、(非)线性方程组、不等式(组)、函数/数列极限、求和式、连乘式
Python中使用sympy库解方程
hou09tian的博客
07-27 3025
sympy库是基于Python语言的符号数学计算系统，它的特点是使用符号而不是数字来处理数学问题。
Python Sympy：解方程利器
最新发布
aobulaien001的博客
10-10 4659
本篇主要介绍Sympy在基本的代数运算方面的应用，Sympy最大的特点是让我们可以用数学专业的思维方式来写代码。它的代码与数学的联系非常紧密，数学专业的朋友如果刚接触编程，用这个库来写代码一定会非常亲切。
python数学建模之用sympy.solve求解方程组的解
qq_51570094的博客
04-09 1万+
在sympy.solve（expression）方法的帮助下，我们可以很容易地求解数学方程，它将返回使用sympy.solve（）方法作为参数提供的方程的根。 参考文档： 参考文档https://www.geeksforgeeks.org/python-sympy-solve-method/在下面这个例子中，我们可以看到通过使用sympy.solve（）方法，我们可以求解数学表达式，这将返回该方程的根。 首先将变量符号化，然后在求解。 例1：求解方程组的解，结果是-2，2. from sy..
python解决经典数学问题-使用Python的SymPy库解决数学运算问题的方法
weixin_39969143的博客
11-11 952
摘要：在学习与科研中，经常会遇到一些数学运算问题，使用计算机完成运算具有速度快和准确性高的优势。Python的Numpy包具有强大的科学运算功能，且具有其他许多主流科学计算语言不具备的免费、开源、轻量级和灵活的特点。本文使用Python语言的NumPy库，解决数学运算问题中的线性方程组问题、积分问题、微分问题及矩阵化简问题，结果准确快捷，具有一定的借鉴意义。1.Sympy库简介SymPy一个用于符...
python用sympy解初等方程
weixin_41855010的博客
02-23 1296
我们在初高中乃至大学里面学到的大多数方程，都可以用计算机来求解。在python中，就提供了sympy方法解方程。流程大致可以分成以下几个步骤： 一、导入sympy包体。这个是常规操作，使用improt语句即可 二、设置未知数。我们使用sympy的.Symbol()方法来设置未知数。 三、列方程并且移项获得表达式。我们需要把一个方程移项成equation = 0的范式，然后把equation作为方程...
python/sympy求解矩阵方程的方法
09-19
Sympy是Python的一个库，用于符号计算，它允许用户进行代数方程的解析求解，微积分，矩阵运算等多种数学运算。在处理线性代数问题，特别是矩阵方程时，Sympy提供了强大的支持。 矩阵方程是形式为AX = B的方程，其中...
python中sympy库求常微分方程的用法
12-20
在Python编程环境中，`sympy`库是一个强大的符号计算工具，它允许用户进行复杂的数学操作，包括求解常微分方程（ODE）。本文将详细介绍如何使用`sympy`来解决常微分方程，并展示两个具体的例子，同时讨论如何通过`...
Sympy 模块解数学方程
abc_xian的博客
04-14 935
Sympy 模块解数学方程
python之sympy库--在线性代数领域的应用
热门推荐
yifengchaoran的专栏
03-06 11万+
sympy作为相对比较全的数学计算库，其也包含针对线性代数的符号运算部分，本文着重介绍sympy在处理线性代数相关问题时的使用方法，且基本严格结合线性代数教材（同济大学版），便于大家回顾，如果想了解sympy在初等代数或微积分方面的应用，可以看文章《python之sympy库--数学符号计算与绘图必备》。 一、矩阵运算 1.1 创建矩阵 创建矩阵是使用sympy处理线性代数问题的起点，以下主要介绍通用创建矩阵的方式以及快速创建特殊矩阵的方式，且一部分主要对应于线性代数教材（同济大学版）的第一章和第二章
python 解方程 sympy_Python数据处理篇之Sympy系列(五)---解方程
weixin_28968285的博客
02-09 5026
前言sympy不仅在符号运算方面强大，在解方程方面也是很强大。本章节学习对应官网的：Solvers官方教程(一)求解多元一次方程-solve()1.说明：解多元一次方程可以使用solve()，在sympy里，等式是用Eq()来表示，例如：2x=42x=4 表示为：Eq(x*2, 4)2.源代码："""解下列二元一次方程2x-y=33x+y=7"""# 导入模块from sympy import *...
gitmodel学习笔记（一）：利用sympy分析高等数学
submarineas的博客
06-17 988
本篇是在复习到中途参加的一个关于gitmodel的学习活动，本系列分为三个部分，分别为高等数学、线性代数以及概率论与数理统计。本篇为第一篇——利用sympy分析高等数学，看完活动文档，查找了相关资料后，汇成笔记在这里记录一下。sympy包相当于让python具备了MATLAB与mathematica相同的解部分数学问题的能力，但是也只是能初步解决一些并不复杂的例子，而且因为还是依赖于python的输出终端，感觉output打印出来的结果没有经过太多美化，不如mathematica足够赏心悦目，但至少还是可以
sympy 求微分方程_快速入门 - SymPy（高数篇）
weixin_28965077的博客
12-30 3360
SymPy 是完全使用 Python 编写的一个符号计算库，也是著名开源计算软件 SageMath 的重要组成部分，是数学学习和教学的的好帮手。其功能主要涵盖以下几个方面：多项式的合并、展开、化简微积分运算微分方程求解线性代数运算系数匹配绘图好了，下面列举一些我整理的范例。声明：SymPy 只是一个库，不要拿来和 Mathematica 这类专业软件比。把它当作科学计算的辅助工具，做数值求解前的公...
高等数学--微分方程的求解（sympy库）
m0_54510474的博客
10-28 3111
sympy求解微分方程的方式，sympy求解实际中的微分方程数学模型
对一个线性方程组的求解（用sympy进行符号运算）
HMMPO的博客
11-15 5041
sympy是python的数学符号计算库，用它可以进行数学表达式的符号推导和演算。与一些成熟的专业符号运算软件相比，sympy功能以及速度比不上，但是对于一些计算工具来说，使用简单，体积小。 举个简单的例子 from sympy import * a, b, c, x = symbols("a, b, c, x") solve(a*x**2 + b*x + c, x) 求解一个一元
数据科学之sympy库
酒心★
05-21 6129
sympy是一个Python的科学计算库，用一套强大的符号计算体系完成诸如多项式求值、求极限、解方程、求积分、微分方程、级数展开、矩阵运算等等计算问题。 Python以其语法简单、易上手、丰富的三方库生态，个人认为可以更优雅地解决日常生活、工作遇到的各种计算问题。 安装： pip install sympy 虚数单位i sympy.I #i平方 sympy.I**2 -1的平方根 sympy....
数学神器！Sympy 模块解数学方程解微积分
静觅
11-23 1067
这是「进击的Coder」的第 758篇技术分享作者：Ckend来源：Python 实用宝典“ 阅读本文大概需要 6 分钟。 ”SymPy 是一个 Python 库，专注于符号数学，它的目标是成为一个全功能的计算机代数系统，同时保持代码简洁、易于理解和扩展。举一个简单的例子，比如说展开二次方程：fromsympy import* x = Symbol('x') y = Symbol('y') ...
python解复杂方程_5.5Python数据处理篇之Sympy系列(五)---解方程
weixin_39957265的博客
11-28 1295
目录前言（一）求解多元一次方程-solve()1.说明：解多元一次方程可以使用solve()，在sympy里，等式是用Eq()来表示，例如：\(2x=4\) 表示为：Eq(x*2, 4)2.源代码："""解下列二元一次方程2x-y=33x+y=7"""# 导入模块from sympy import *# 将变量符号化x = Symbol('x')y = Symbol('y')z = Symbol...
sympy解微分方程
09-21
好的，我可以为您演示如何使用SymPy来解微分方程。首先，您需要安装SymPy包并导入它。您可以使用以下命令来安装SymPy: ``` !pip install sympy ``` 然后，您可以在Python中导入SymPy： ```python import sympy as...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

小孟的CDN

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值