数模学习目录

  • SciPy 学习
    • SciPy基本操作
      • 1-求解非线性方程（组）
        • 解题代码如下：
        • 运行x1、x2结果如下：
  • 解题代码如下：
      • 运行s2效果如下：
    • 2-积分
    • 3-最小二乘解
    • 4-最大模特征值及对应的特征向量
      • 运行结果如下：
  • Numpy学习（续）
    • 1 Numpy 数学函数
      • 1-三角函数
      • 2-舍入函数
        • 2-1 numpy.around()
        • 2-2 numpy.floor()
        • 2-3 numpy.ceil()
    • 2 Numpy算术函数
      • 相关代码如下：
  • Pandas学习（续）
    • # pandas的SettingWithCopyWarning
    • Pandas 数据排序
      • Series的排序：
        • 运行结果分别如下：
    • DataFrame的排序
      • 运行部分结果如下：
    • Pandas字符串处理
      • 遇到的问题：
    • 总结：
      • 每日一言：

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

SciPy 学习

# -*- coding = utf-8 -*-
# @Time : 2022/6/30 8:20
# @Author : lxw_pro
# @File : SciPy.py
# @Software : PyCharm'''SciPy 包含的模块有最优化、线性代数、积分、插值、特殊函数、快速傅里叶变换、
信号处理和图像处理、常微分方程求解和其他科学与工程中常用的计算。'''
# 安装scipy库：
# SciPy终端安装命令：pip install SciPy# https://www.runoob.com/w3cnote/python-pip-install-usage.html  Python pip 安装与使用# 查看scipy版本：
import scipyprint(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 constantsprint(dir(constants))'''
单位类型
常量模块包含以下几种单位：公制单位
二进制，以字节为单位
质量单位
角度换算
时间单位
长度单位
压强单位
体积单位
速度单位
温度单位
能量单位
功率单位
力学单位'''print()
# SciPy 常量模块:
# constants 是 scipy 的常量模块
from scipy import constants# 查看一英亩等于多少平方米：
print(constants.acre)   # 输出 4046.8564223999992# SciPy 常量模块 constants 提供了许多内置的数学常数
# 圆周率：  pi
# 黄金比例： goldenfrom scipy import constantsprint(constants.pi)         # 输出 3.141592653589793 【圆周率】print(constants.golden)     # 输出 1.618033988749895 【黄金比例】

SciPy基本操作

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.388x1 = 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, rootfs2 = 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效果如下：

2-积分

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

积分主要分为以下两类：

  1. 对给定函数的数值积分
  2. 对给定离散点的数值积分，函数有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 quaddef func(x, a, b):return a*x**2+b*xz1 = 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)# 注：输出的后一个值为积分值的绝对误差

3-最小二乘解

# 最小二乘解# scipy.optimize 模块求非线性方程组最小二乘解格式：
'''
from scipy.optimize import least_squares
least_squares(fun, x0)注：用到loadtxt需自行准备好文件【准备文件】
'''from scipy.optimize import least_squares
import numpy as nps = 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)

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

题目

# 4-最大模特征值及对应的特征向量# 题目描述：求下列矩阵的最大模特征值及对应的特征向量：from scipy.sparse.linalg import eigs
import numpy as npm = 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学习（续）

# -*- coding = utf-8 -*-
# @Time : 2022/6/30 14:20
# @Author : lxw_pro
# @File : Numpy 广播-7.py
# @Software : PyCharm# NumPy 广播(Broadcast)
# 广播是 numpy 对不同形状的数组进行数值计算的方式， 对数组的算术运算通常在相应的元素上进行。# 如果两个数组 a 和 b 形状相同，即满足 a.shape == b.shape，那么 a*b 的结果就是 a 与 b 数组对应位相乘。
# 这要求维数相同，且各维度的长度相同。'''
对两个数组，分别比较他们的每一个维度（若其中一个数组没有当前维度则忽略），满足：数组拥有相同形状。
当前维度的值相等。
当前维度的值有一个是 1。若条件不满足，抛出 "ValueError: frames are not aligned" 异常'''import numpy as npa = np.array([3, 6, 9])
b = np.array([2, 4, 6])
c = a * b
print(c)    # 输出 [ 6 24 54]# 若形状不同时，numpy 将自动触发广播机制
import numpy as npx = 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-三角函数

# NumPy 数学函数
# NumPy 包含大量的各种数学运算的函数，包括三角函数，算术运算的函数，复数处理函数等。# 1-三角函数
# NumPy 提供了标准的三角函数：sin()、cos()、tan()。
import numpy as nplxw = 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 npbl = 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 npxx = 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 npxs = np.array([23.1, 23.5, 54.9])
print(np.ceil(xs))      # 输出 [24. 24. 55.]

2 Numpy算术函数

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

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

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

相关代码如下：


import numpy as npsz = 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 npm = 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, m2)       # 数组1作为底数，数组2作为幂
print(mc2)  # 输出 [  1  16 512]print()
# numpy.mod()# numpy.mod() 计算输入数组中相应元素的相除后的余数
# 函数 numpy.remainder() 也产生相同的结果
import numpy as npsz1 = np.array([23, 45, 67])
sz2 = np.array([2, 3, 5])print(np.mod(sz1, sz2))         # 输出 [1 0 2]print(np.remainder(sz1, sz2))   # 输出 [1 0 2]

Pandas学习（续）

# pandas的SettingWithCopyWarning

# pandas的SettingWithCopyWarning报警复现、原因、解决方案# 读取数据
import pandas as pddf = pd.read_csv('nba.csv')print(df.head())# 核心解决问题：pandas的dataframe的修改写操作，只允许在源dataframe上进行，一步到位# 解决方法（两种）：
'''
1-将get+set的两步操作，改成set的一步操作
2-若须处理筛选数据做后续的处理分析，使用copy复制dataframe'''# pandas不允许先筛选子dataframe,在进行修改写入

【注意先准备好csv文件】

Pandas 数据排序

Series的排序：

# -*- coding = utf-8 -*-
# @Time : 2022/6/30 15:42
# @Author : lxw_pro
# @File : pandas 数据排序-9.py
# @Software : PyCharm# Pandas 数据排序'''
Series的排序：
Series.sort_values(ascending=True, inplace=False)
参数说明：· ascending: 默认为True升序排序，False为False· inplace:   是否修改原始SeriesDataFrame的排序：
DataFrame.sort_values(by, ascending=True, inplace=False)
参数说明：· by:字符串或者List<字符串>,单列排序或者多列排序· ascending: bool或者List，升序还是降序· inplace:   是否修改原始DataFrame'''# Series的排序：
import pandas as pddf = pd.read_csv('nba.csv')
print(df.head())        # 输出前五行print(df['Weight'].sort_values())                      # 升序排序print(df['Weight'].sort_values(ascending=False))       # 降序排序

运行结果分别如下：

DataFrame的排序

# DataFrame的排序# 单列排序：
print(df.sort_values(by='Weight'))

运行部分结果如下：

print(df.sort_values(by="Weight", ascending=False))     # 降序排序


# 多列排序：
print(df.sort_values(by=['Age', 'Weight']))


# 两个字段都是降序排序
print(df.sort_values(by=['Age', 'Weight'], ascending=False))


# 分别指定升序还是降序
print(df.sort_values(by=['Age', 'Weight'], ascending=[False, True]))

Pandas字符串处理

之前我就在这个字符串处理的题上出了一些问题(不过当天就解决啦)【今天在来看看】，也就是

df['lrl'].str.replace("%", "").astype("int32")

# Pandas字符串处理：
'''
1-使用方法：先获取Series的属性，然后再属性上调用函数
2-只能在字符串列上使用，不能再数字列上使用
3-DataFrame没有str属性和使用
4-Series.str并不是原生Python字符串，它是封装的一套方法'''# 获取Series的属性
# print(df['Salary'].str)   # 报错【示范】# AttributeError: Can only use .str accessor with string values!
# AttributeError：只能使用。带字符串值的str访问器！# 一定得是字符串列
print(df['College'].str)
# 运行结果为： <pandas.core.strings.accessor.StringMethods object at 0x00000204444EBC48># 判断是不是数字列
print(df['College'].str.isnumeric())# print(df['College'].len)      # 报错【示范】# AttributeError: 'Series' object has no attribute 'len'
# AttributeError：“Series”对象没有属性“len”

# 使用str的startswith、contains等得到bool的Series可以做条件查询
tj = df['Height'].str.startswith("6-2")
print(tj)# print(df[tj].head())    # ValueError：无法使用包含NA/NaN值的非布尔数组进行掩码


# 去掉Height中间的“-”
print(df['Height'].str.replace("-", ""))


# 取第一位数
print(df['Height'].str.replace("-", "").str.slice(0, 1))# 同上
print(df['Height'].str.replace("-", "").str[0:1])

遇到的问题：

  • 自己弄的文件和上面的文件不同，做起来还是增加了点难度，不过那都不是事儿，主要能解决就好，这样错多了也可以更加熟练掌握，是吧！

总结：

今天我学习到了新的库SciPy，挺有意思的，但同时要想学好，确实也不是那么容易的，当然，我也不会止步于此的，继续向前，向前走，向前冲！人生，不就是走走停停嘛，这样才有趣，对吧！！

每日一言：

所谓有趣的灵魂，实际上就是这个人的信息密度和知识层面，都远高于你，并愿意俯下身来听你说那毫无营养的废话和你交流，提出一些你没有听过的观点，颠覆了你短浅的想象力及三观。

【数学建模学习笔记【集训十天】之第四天】相关推荐

  1. 数学建模学习笔记（十）——时间序列模型

    文章目录 一.时间序列综述 二.时间序列数据以及基本概念 三.时间序列分解 四.指数平滑模型 五.一元时间序列分析的模型 六.AR(p)模型 七.MA(q)模型 八.ARMA(p, q)模型 九.模型 ...

  2. 数学建模学习笔记（十二）——奇异值分解

    文章目录 一.综述 二.奇异值分解 三.使用SVD进行降维 四.SVD的评价及应用 一.综述 奇异值分解是线性代数中一种重要的矩阵分解方法,它在图形的压缩等方面具有重要的意义及作用. 二.奇异值分解 ...

  3. 数学建模学习笔记（十五）数据归一化

    在神经网络构建之前,需要对数据进行归一化处理,为什么要归一化? 因为不同数据范围不同,比如一个特征的数据范围为(1,5),另一个为(100,1000),会导致特征之间对结果的影响不同,因此需要将它们归 ...

  4. 数学建模学习笔记（十六）lingo求解整数规划

    之前没用过lingo,这次遇到一个线性整数规划的问题,尝试用了lingo求解,发现上手还是很容易的.将题目和求解放置在此,以便查阅~ 问题: 使用lingo求解,输入以下代码 model: min=1 ...

  5. 【数学建模学习笔记【集训十天】之第六天】

    数模学习目录 Matplotlib 学习 Matplotlib简介 Matplotlib 散点图 运行效果如下: Matplotlib Pyplot 运行效果如下: 关于plot() 运行效果如下: ...

  6. 数学建模学习笔记（2.3）lingo软件求解线性规划问题

    数学建模学习笔记(2.3)lingo软件求解线性规划问题 lingo软件的优势在于体积小,专注于解决优化问题 且编程语言通俗易懂,没有门槛 对于刚刚接触数学建模同学比较友善 当然对于已经参与建模很久的 ...

  7. 数学建模学习笔记（1）数学模型的特点和分类

    数学建模学习笔记(1)数学模型的特点和分类 ps:学习的教材为姜启源著的<数学模型(第四版)> 领取数模资料和更多内容请关注公众号:拾壹纪元 传送门: 线性规划(LP)问题 https:/ ...

  8. 数学建模学习笔记之评价问题聚类分析法

    数学建模学习笔记之评价问题聚类分析法 物以类聚.人以群分. 聚类分析是一个很大的概念,显然根据分类的依据不同会出现很多很多聚类的方法.例如K-Means .Sequential Leader.Mode ...

  9. 【OS学习笔记】十六 保护模式四：进入保护模式与在保护模式下访问内存的汇编代码

    本文记录的是之前四篇文章所对应的汇编代码.四篇文章分别是: [OS学习笔记]十二 现代处理器的结构和特点 [OS学习笔记]十三 保护模式一:全局描述符表(GDT) [OS学习笔记]十四 保护模式二:段 ...

  10. 数学建模学习笔记（一）:插值法

    文章目录 前言 一.一维插值问题的描述 二.常用插值方法 1.Lagrange插值法 2.Newton插值法 三.高次插值的Runge现象 四.Matlab插值 1.一维插值 2.二维插值 3.散乱点 ...

最新文章

  1. xslt中的Javascript取得xml中的参数
  2. 请重新认识你作为程序员的价值
  3. python if main_python中if __name__ == '__main__' ：main(()
  4. 使用脚本安装elasticsearch7.3的记录
  5. 【学术相关】是什么让女性在计算机史上“隐身”了？
  6. Go语言反射之反射调用
  7. 程序员专属段子集锦 5/10
  8. 从头学习Drupal--基本架构三
  9. android 中LayoutInflater 的使用
  10. 你关注过浏览器最小字体为多大吗？
  11. 网络社区划分的算法分类
  12. 说说几种常用的前端缓存
  13. html tbody id 赋值,IE6/IE7/IE8/IE9中tbody的innerHTML不能赋值的完美解决方案
  14. 【怀念】清明遥祭外婆
  15. 我以前的痛苦，你也拥有吗?
  16. http://www.horizon3d.com.cn
  17. python爬取B站动态的评论总数(不含用户评论内容详情)
  18. JSP+ssm计算机毕业设计快递代取系统的设计与实现94wcw【源码、数据库、LW、部署】
  19. 体验最近火爆的ChatGPT，真的被震惊到了
  20. 学习笔记——Kaggle_Digit Recognizer (朴素贝叶斯 Python实现）

热门文章

  1. SATA控制器学习资料整理
  2. Rust能力养成之（10）用Cargo进行项目管理：扩展 调用与优化
  3. 人工智能生态环境预测_人工智能在生态领域的应用研究
  4. [WARNING] The POM for com.tenyears:base-common:jar:1.0 is invalid, transitive dependen
  5. abaqus中python版本_ABAQUS中的python语言入门
  6. 实验10-1 圆形体体积计算器
  7. 自然对数e的推导过程
  8. 腾讯招聘日――招聘新体验
  9. 透过三季度财报看3.0用友：向云而升，做企业数智化的云引擎
  10. linux 各个SIG信号含义