您现在的位置是：网站首页> 编程资料编程资料

python数学建模(SciPy+ Numpy+Pandas)_python_

2023-05-26 534人已围观

简介 python数学建模(SciPy+ Numpy+Pandas)_python_

前言

SciPy 是一个开源的 Python 算法库和数学工具包, Scipy 是基于 Numpy 的科学计算库，用于数学、科学、工程学等领域。
SciPy库依赖于NumPy，它提供了便捷且快速的N维数组操作。
SciPy库的构建与NumPy数组一起工作，并提供了许多用户友好和高效的数字实践，
并在天文学、生物学、气象学和气候科学，以及材料科学等多个学科得到了广泛应用。
NumPy和SciPy易于使用，强大到足以依靠世界上一些顶尖的科学家和工程师。

SciPy 学习

''' SciPy 包含的模块有最优化、线性代数、积分、插值、特殊函数、快速傅里叶变换、 信号处理和图像处理、常微分方程求解和其他科学与工程中常用的计算。 ''' # 安装scipy库： # SciPy终端安装命令：pip install SciPy # https://www.runoob.com/w3cnote/python-pip-install-usage.html Python pip 安装与使用 # 查看scipy版本： import scipy print(scipy.__version__) # SciPy模块功能表 ''' 模块 功能 scipy.cluster 聚类分析等 scipy.constants 物理和数学函数 scipy.fftpack 傅里叶变换 scipy.integrate 积分 scipy.interpolate 插值 scipy.io 数据输入和输出 scipy.linalg 线性代数 scipy.ndimage n维图像 scipy.odr 正交距离回归 scipy.optimize 优化 scipy.signal 信号处理 scipy.sparse 稀疏矩阵 scipy.spatial 空间数据结构和算法 scipy.special 特殊函数 scipy.stats 统计 ''' # 使用 dir() 函数来查看 constants 模块包含的常量： from scipy import constants print(dir(constants)) ''' 单位类型 常量模块包含以下几种单位： 公制单位 二进制，以字节为单位 质量单位 角度换算 时间单位 长度单位 压强单位 体积单位 速度单位 温度单位 能量单位 功率单位 力学单位 ''' print() # SciPy 常量模块: # constants 是 scipy 的常量模块 from scipy import constants # 查看一英亩等于多少平方米： print(constants.acre) # 输出 4046.8564223999992 # SciPy 常量模块 constants 提供了许多内置的数学常数 # 圆周率： pi # 黄金比例： golden from scipy import constants print(constants.pi) # 输出 3.141592653589793 【圆周率】 print(constants.golden) # 输出 1.618033988749895 【黄金比例】

1-SciPy基本操作

1-1求解非线性方程（组）

1-1

解题代码如下：

# scipy.optimize模块的fsolve和root可求非线性方程（组）的解 # 格式： from scipy.optimize import fsolve from scipy.optimize import root # fsolve或root求解非线性方程组时，先把非线性方程组写成 F(x)=0 这样的形式【x:向量；F(x)：向量函数】 fx = lambda x: x**980-5.01*x**979-3.388*x**977\ +7.398*x**978-x**3+5.01*x**2-7.398*x+3.388 x1 = fsolve(fx, 1.5, maxfev=420) # 函数调用420次【调用小了，会报警告】 x2 = root(fx, 1.5) print(x1) # 相当于答案 print() print(x2) # 相当于解题过程

运行x1、x2结果如下：

1-2

解题代码如下：

from scipy.optimize import fsolve, root fs2 = lambda s: [s[0]**2+s[1]**2-1, s[0]-s[1]] s1 = fsolve(fs2, [1, 1]) print() s2 = root(fs2, [1, 1]) print(s1) # 输出 [0.70710678 0.70710678] print() print(s2)

运行s2效果如下：

1-2积分

scipy.integrate模块提供了多种积分模式。

积分主要分为以下两类：

对给定函数的数值积分
对给定离散点的数值积分，函数有trapz

题目：

''' 函数 说明 quad(func, a, b, args) 计算一重数值积分 dblquad(func, a, b, gfun, hfun, args) 计算二重数值积分 tplquad(func, a, b, gfun, hfun, qfun, rfun) 计算三重数值积分 nquad(func, ranges, args) 计算多变量积分 ''' from scipy.integrate import quad def func(x, a, b): return a*x**2+b*x z1 = quad(func, 0, 1, args=(2, 1)) z2 = quad(func, 0, 1, args=(2, 10)) print(z1) # 输出 (1.1666666666666665, 1.2952601953960159e-14) print(z2) # 输出 (5.666666666666667, 6.291263806209221e-14) # 注：输出的后一个值为积分值的绝对误差

1-3最小二乘解

# 最小二乘解 # scipy.optimize 模块求非线性方程组最小二乘解格式： ''' from scipy.optimize import least_squares least_squares(fun, x0) 注：用到loadtxt需自行准备好文件【准备文件】 ''' from scipy.optimize import least_squares import numpy as np s = np.loadtxt('data.txt') x0 = s[0] y0 = s[1] d = s[2] fs = lambda x: np.sqrt((x0-s[0])**2+(y0-s[1])**2-d) xc = least_squares(fs, np.random.rand(2)) print(xc) print() print(xc.s)

1-4最大模特征值及对应的特征向量

题目：

# 4-最大模特征值及对应的特征向量 # 题目描述：求下列矩阵的最大模特征值及对应的特征向量： from scipy.sparse.linalg import eigs import numpy as np m = np.array([ [1, 2, 3], [2, 1, 3], [3, 3, 6] ], dtype=float) a, b = np.linalg.eig(m) c, d = eigs(m, 1) print('最大模特征值为：', c) # 输出 最大模特征值为： [9.+0.j] print('对应的特征向量：\n', d)

运行结果如下：

Numpy学习（续）

# NumPy 广播(Broadcast) # 广播是 numpy 对不同形状的数组进行数值计算的方式， 对数组的算术运算通常在相应的元素上进行。 # 如果两个数组 a 和 b 形状相同，即满足 a.shape == b.shape，那么 a*b 的结果就是 a 与 b 数组对应位相乘。 # 这要求维数相同，且各维度的长度相同。 ''' 对两个数组，分别比较他们的每一个维度（若其中一个数组没有当前维度则忽略），满足： 数组拥有相同形状。 当前维度的值相等。 当前维度的值有一个是 1。 若条件不满足，抛出 "ValueError: frames are not aligned" 异常 ''' import numpy as np a = np.array([3, 6, 9]) b = np.array([2, 4, 6]) c = a * b print(c) # 输出 [ 6 24 54] # 若形状不同时，numpy 将自动触发广播机制 import numpy as np x = np.array([ [4, 2, 5], [5, 2, 0], [2, 6, 1], [1, 4, 5] ]) y = np.array([3, 1, 2]) print(x+y) yy = np.tile(y, (4, 1)) # 重复b的各个维度 print(x+yy)

1-Numpy 数学函数

1-1三角函数

# NumPy 数学函数 # NumPy 包含大量的各种数学运算的函数，包括三角函数，算术运算的函数，复数处理函数等。 # 1-三角函数 # NumPy 提供了标准的三角函数：sin()、cos()、tan()。 import numpy as np lxw = np.array([0, 30, 45, 60, 90]) # sin() zx = np.sin(lxw*np.pi/180) print(zx) # 计算角度的反正弦【单位：弧度】 fzx = np.arcsin(zx) print(fzx) # 检查结果【通过转化为角度制】 jg = np.degrees(fzx) print(jg) # 输出 [ 0. 30. 45. 60. 90.] # cos() yx = np.cos(lxw*np.pi/180) print(yx) # 反余弦 fyx = np.arccos(yx) print(fyx) # 检查结果： jg2 = np.degrees(fyx) print(jg2) # 输出 [ 0. 30. 45. 60. 90.] # tan() zq = np.tan(lxw*np.pi/180) print(zq) # 反正切 fzq = np.arctan(zq) print(fzq) # 检查结果： jg3 = np.degrees(fzq) print(jg3) # 输出 [ 0. 30. 45. 60. 90.]

2-舍入函数

2-1 numpy.around()

# 2-舍入函数 # 2-1 numpy.around() ''' numpy.around() 函数返回指定数字的四舍五入值。 格式： numpy.around(a,decimals) 参数说明： a: 数组 decimals: 舍入的小数位数。 默认值为0。 如果为负，整数将四舍五入到小数点左侧的位置 ''' import numpy as np bl = np.array([15.222, 22.6555, 13.71111]) print(np.around(bl)) # 输出 [15. 23. 14.] print(np.around(bl, 2)) # 输出 [15.22 22.66 13.71] print(np.around(bl, -1)) # 输出 [20. 20. 10.]

2-2 numpy.floor()

# 2-2 numpy.floor() # numpy.floor() 返回小于或者等于指定表达式的最大整数，即向下取整 import numpy as np xx = np.array([23.3, 13.43, 2.9]) print(np.floor(xx)) # 输出 [23. 13. 2.]

2-3 numpy.ceil()

# 2-3 numpy.ceil() # numpy.ceil() 返回大于或者等于指定表达式的最小整数，即向上取整 import numpy as np xs = np.array([23.1, 23.5, 54.9]) print(np.ceil(xs)) # 输出 [24. 24. 55.]

3-Numpy算术函数

NumPy 算术函数包含简单的加减乘除: add()，subtract()，multiply() 和 divide()

倒数：reciprocal()
幂：power()
余数：mod() | remainder()

注：数组必须具有相同的形状或符合数组广播规则

相关代码如下：

import numpy as np sz = np.arange(9, dtype=np.float_).reshape(3, 3) sz2 = np.array([5, 2, 1]) # 注：如果相除，这里是被除数的话，里面不能有0 # 数组相加 xj = np.add(sz, sz2) print(xj) # 数组相减 xj2 = np.subtract(sz, sz2) print(xj2) # 数组相乘 xc = np.multiply(sz, sz2) print(xc) # 数组相除 xc2 = np.divide(sz, sz2) print(xc2) print() # numpy.power() # numpy.power() 函数将第一个输入数组中的元素作为底数，计算它与第二个输入数组中相应元素的幂 import numpy as np m = np.array([1, 4, 8]) # 数组1 mc = np.power(m, 3) # 数组1所有元素对应的3次方 print(mc) # 输出 [ 1 64 512] m2 = np.array([1, 2, 3]) # 数组2 mc2 = np.power(m, m
                提示：
                    本文由神整理自网络，如有侵权请联系本站删除！
                    

                    本站声明： 

                    1、本站所有资源均来源于互联网，不保证100%完整、不提供任何技术支持； 

                    2、本站所发布的文章以及附件仅限用于学习和研究目的;不得将用于商业或者非法用途；否则由此产生的法律后果，本站概不负责！
                
                
                
                                            上一篇：python数学建模是加深Numpy和Pandas学习_python_
                                                                下一篇：python数学建模之Matplotlib 实现图片绘制_python_