1.方阵的迹

方阵的迹就是方阵的主对角线元素之和

# -*- coding:utf-8 -*-
# @Author: WanMingZhu
# @Date: 2019/8/12 9:37 import numpy as np arr = np.random.randint(1, 5, size=(4, 4))
print(arr)
# 调用np.trace便可求出矩阵的迹
print(np.trace(arr))
"""
[[1 4 2 3]
[3 2 2 2]
[3 2 3 3]
[1 4 2 4]]
10
"""

2.如何求两个矩阵之间的距离

import numpy as np

arr1 = np.array([[0, 1],
[1, 0]]) arr2 = np.array([[1, 1],
[1, 1]])
# 先计算a和b的差,得出arr3
arr3 = arr1 - arr2
# 然后让arr3和arr3本身进行点乘
arr4 = arr3 @ arr3
# 此时arr4的迹便是距离的平方
print(np.trace(arr4)) # 2

3.如何求出一个矩阵的逆矩阵和伴随矩阵

import numpy as np

a = np.array([[1, -2, 1],
[0, 2, -1],
[1, 1, -2]])
# 直接对a使用np.linalg.inv即可求出逆矩阵
a_inv = np.linalg.inv(a)
print(a_inv)
"""
[[ 1. 1. 0. ]
[ 0.33333333 1. -0.33333333]
[ 0.66666667 1. -0.66666667]]
""" # 使用np.linalg.det(a)求出行列式的值
a_det = np.linalg.det(a)
print(a_det) # -3.0000000000000004 # 然后逆矩阵a_inv和行列式的值a_det相乘即可得出伴随矩阵
print(a_inv * a_det)
"""
[[-3. -3. -0.]
[-1. -3. 1.]
[-2. -3. 2.]]
"""

4.如何解多元一次方程

import numpy as np

"""
x+2y+z=7
2x-y+3z=7
3x+y+2z=18
求这个方程组的解
"""
# 首先将系数写下来,排成一个矩阵
a = np.array([[1, 2, 1],
[2, -1, 3],
[3, 1, 2]])
# 将右边的常数写下来,排成一个矩阵
b = np.array([7, 7, 18])
# 求解,将参数传进去
x = np.linalg.solve(a, b)
print(x) # [ 7. 1. -2.]
# 验证
print(a @ x) # [ 7. 7. 18.]

5.求矩阵的秩

import numpy as np

# 生成4行4列的单位矩阵
i = np.eye(4)
print(i)
"""
[[1. 0. 0. 0.]
[0. 1. 0. 0.]
[0. 0. 1. 0.]
[0. 0. 0. 1.]]
"""
# 求出秩
print(np.linalg.matrix_rank(i)) # 4 # 修改一个值
i[0, 0] = 0
print(np.linalg.matrix_rank(i)) # 3

6.求协方差矩阵

import numpy as np

a = np.array([170, 180, 175])
b = np.array([65, 80, 70])
print(np.cov([a, b]))
"""
[[25. 37.5 ]
[37.5 58.33333333]]
"""

7.求相关矩阵

import numpy as np

a = np.array([170, 180, 175])
b = np.array([65, 80, 70])
print(np.corrcoef([a, b]))
"""
[[1. 0.98198051]
[0.98198051 1. ]]
"""

8.生成一个方程

import numpy as np

x = np.poly1d([1, 1, 1, 1])
print(x)
"""
3 2
1 x + 1 x + 1 x + 1
就是x^3 + x^2 + x + 1
"""
# 怎么生成出来的呢?
# 如果矩阵里面有4个元素,那么未知数的最高次幂就是3
# 然后矩阵里面的元素就是对应项的系数
# 比如我想生成 8 * x^5 + 4 * x^3 + 3 * x^2 + x
"""
首先最高次幂是5,就意味着数组里面有6个值
x^5次幂对应的系数是8,说明第一个值是8
没有x^4,说明x^4对应的系数是0,说明第二个元素是0。虽然是0,但是必须要写,不然numpy不知道你要生成最高几次幂的函数
"""
# 别忘记了最后要加上一个0
x = np.poly1d([8, 0, 4, 3, 1, 0])
print(x)
"""
5 3 2
8 x + 4 x + 3 x + 1 x
""" # 另外我向x里面传值,还可以根据当前的自变量计算出对应的函数
print(x(3)) # 2082
print(8 * 3 ** 5 + 4 * 3 ** 3 + 3 * 3 ** 2 + 3) # 2082

9.求出一个函数的导数

import numpy as np

x = np.poly1d([2, 1, 1])
print(x)
"""
2
2 x + 1 x + 1
"""
# 调用deriv进行求导
x1 = x.deriv()
print(x1) # 4 x + 1
# 显然这个结果是正确的
# 同样的,我们也可以计算相应的导数值
print(x1(3)) # 13

numpy中与高等数学有关的函数的更多相关文章

  1. numpy中与金融有关的函数

    fv函数 计算未来的价值 def fv(rate, nper, pmt, pv, when='end'): ... 参数: rate:存款/贷款每期的利率 nper:存款/贷款期数 pmt:存款/贷款 ...

  2. numpy中的argpartition

    numpy.argpartition(a, kth, axis=-1, kind='introselect', order=None) 在快排算法中,有一个典型的操作:partition.这个操作指: ...

  3. python -- numpy 基本数据类型,算术运算,组合,分割 函数

    0 NumPy数组 NumPy数组:NumPy数组是一个多维数组对象,称为ndarray.其由两部分组成: 实际的数据 描述这些数据的元数据 NumPy数组属性: ndim(纬数,x,y 2),sha ...

  4. numpy中linspace用法 (等差数列创建函数)

    linspace  函数 是创建等差数列的函数, 最好是在 Matlab  语言中见到这个函数的,近期在学习Python 中的 Numpy, 发现也有这个函数,以下给出自己在学习过程中的一些总结. ( ...

  5. Numpy中Meshgrid函数介绍及2种应用场景

    近期在好几个地方都看到meshgrid的使用,虽然之前也注意到meshgrid的用法.但总觉得印象不深刻,不是太了解meshgrid的应用场景.所以,本文将进一步介绍Numpy中meshgrid的用法 ...

  6. numpy中的arg系列函数

    numpy中的arg系列函数 觉得有用的话,欢迎一起讨论相互学习~Follow Me 不定期更新,现学现卖 numpy中arg系列函数被经常使用,通常先进行排序然后返回原数组特定的索引. argmax ...

  7. 内置函数和numpy中的min(),max()函数

    内置min()函数 numpy中的min()函数:

  8. Numpy中扁平化函数ravel()和flatten()的区别

    在Numpy中经常使用到的操作由扁平化操作,Numpy提供了两个函数进行此操作,他们的功能相同,但在内存上有很大的不同. 先来看这两个函数的使用: from numpy import * a = ar ...

  9. numpy中的快速的元素级数组函数

    numpy中的快速的元素级数组函数 一元(unary)ufunc 对于数组中的每一个元素,都将元素代入函数,将得到的结果放回到原来的位置 >>> import numpy as np ...

随机推荐

  1. CSS display的几个常用的属性值,inline , block, inline-block

    1.解释一下display的几个常用的属性值,inline , block, inline-block inline(行内元素): 使元素变成行内元素,拥有行内元素的特性,即可以与其他行内元素共享一行 ...

  2. Cinder AZ 与 Nova AZ 的同步问题

    问题 今天处理了一个 Boot from volume 失败的问题,错误日志给出了明确的原因:The instance and volume are not in the same AZ. Build ...

  3. ServletRequest与ServletResponse

    http://lavasoft.blog.51cto.com/62575/275586/ 请求和相应是Web交互最基本的模式,在Servlet中,分别用HttpServletRequest与HttpS ...

  4. Eclipse使用git发布项目到github

    因为一直都在使用svn,今天尝试了下git,记录下来既是方便自己以后查看,也是分享一些经验! 废话不多说,撸起袖子就是干!!! 1.选中要上传的项目右键 2.选中git 3.在图上打钩,点击所要上传的 ...

  5. pagehelper用法

    spring boot 和 mybatis 中使用 pagehelper: 1. 在pom.xml 中加入pagehelper的引用: <dependency> <groupId&g ...

  6. Java连接Sap系统调并调用RFC函数

    参考博客:https://blog.csdn.net/qq_36026747/article/details/81287462                   https://www.cnblog ...

  7. 安装Redis-cluster-gem install redis报错的解决方案

    错误描述: [root@eshop-cache01 local]# gem install redis ERROR: Loading command: install (LoadError) cann ...

  8. JAVA学习记录<一>

    一: JAVA初体验: 1.JAVA简介: 2.环境搭建: 3:MyEclipse的使用简介: 4:程序的移植:项目的导入,导出. 5:学习JAVA的经验: 多写,多问,总结和复习!!!

  9. Bat:IP切换,内外网切换,路由设置内外网同时连接

    1.IP切换: @echo off ::双冒号表示注释,注释只能单独一行写,@表示不让执行的命令又显示在屏幕上 cls ::清屏 color c ::设置显示字体颜色 set IP=10.10.17. ...

  10. pytorch中的激励函数(详细版)

          初学神经网络和pytorch,这里参考大佬资料来总结一下有哪些激活函数和损失函数(pytorch表示)      首先pytorch初始化:   import torch import t ...