python矩阵、解方程、绘图操作

小电动车

已于 2022-04-30 16:33:40 修改

阅读量1.7k

点赞数 1

分类专栏： python 文章标签： python

于 2022-04-13 11:53:26 首次发布

本文链接：https://blog.csdn.net/weixin_46143152/article/details/124144785

版权

python 专栏收录该内容

7 篇文章 0 订阅

订阅专栏

这篇博客介绍了如何使用numpy进行矩阵操作，包括创建、转置和求逆。接着利用sympy库解了一组线性方程，并提取了方程解的值。最后，展示了使用matplotlib进行图形绘制的过程，画出了两条曲线并显示了图像。

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

矩阵操作

from numpy import * # 矩阵操作库
V = mat([xd,yd,zd]) # 创建13矩阵
V = V.T # 转置
V = V.I # 求逆
print(diag([1,2,2,1])) # 生成对角数组

解方程

from sympy import solve # 解方程
from sympy.core.symbol import symbols # 定义参数

a0,a1,a2,a3,a4... = sympy.symbols('a0 a1 a2 a3 a4...') # 定义参数
# 列写方程（下面的不全
A = a0 - 20
B = a1 - 0
C = a2 - 0
D = a0 + a1 + a2 + a3 + a4 - 60
# ...

result = sympy.solve([A,B,C,D...],[a0,a1,a2,a3,a4...]) # 求解得到一个字典
theta1 = float(result[a0]) # 取值
for value in result.values(): # 作为数组打印
    print(float(value))

绘图

import matplotlib.pyplot as plt # 导入绘图库
import numpy as np
xs = 2
ys = 6
xf = 12
yf = 3

x = np.linspace(xs,xf,11) # 自变量
y = (yf-ys)/(xf-xs)*(x-xf)+yf # 分解路径为10段

J1 = [2.14,-6.79,-14.6,-21.1,-26.1,-29.7,-32.2,-33.6,-34  ,-33.7,-32.6]
J2 = [139,138 ,136 ,133 ,129 ,125 ,120 ,114 ,108 ,101 ,93.2]

plt.plot(x,J1) # 作图
plt.plot(x,J2)
plt.show() # 显示图片窗口