python符号运输sympy 每日学习一点点

'''符号运算又称计算机代数，通俗地讲就是用计算机推导数学公式，如对表达式进行因式分解、化简、微分、积分、解代数方程、求解常微分方程等。与数值运算相比，符号计算存在以下的特点：
 （1）运算以推理方式进行，因此不受截断误差和累积误差问题的影响；
（2）符号计算的速度比较慢。
'''
import sympy as sp
# 定义符号变量符号函数
x, y, z, a = sp.symbols('x, y, z, a')  # 可以去逗号，保留空格
f, g = sp.symbols('f, g', cls=sp.Function)  # 定义符号函数    y = sp.Function('y')

x0 = sp.solve(a*x**2 + y*x + z, x)
print(x0)

import sympy as sp
 x1, x2 = sp.symbols('x1, x2')
x1, x2 = sp.var('x1,x2')
s = sp.solve([x1**2 + x2**2-1, x1-x2], [x1, x2])
print(s)

# 求特征值特征向量  matrix  矩阵
a = sp.Matrix([[0, 0, 0, 1], [0, 0, 1, 0], [0, 1, 0, 0], [1, 0, 0, 0]])
print("特征值为：", a.eigenvals())
print("特征向量为：", a.eigenvects())


# 绘制三视图
import matplotlib.pyplot as plt
import numpy as np

plt.rc('font', size=16)
plt.rc('font', family='SimHei')
x = np.linspace(-4, 4, 100)
x, y = np.meshgrid(x, x)
z = 50 * np.sin(x + y)
ax = plt.axes(projection='3d')
ax.plot_surface(x, y, z, color='y')
plt.show()