【NumPy】NumPy线性代数模块详解:掌握numpy.linalg的核心功能

news2024/11/15 8:34:59

🧑 博主简介:阿里巴巴嵌入式技术专家,深耕嵌入式+人工智能领域,具备多年的嵌入式硬件产品研发管理经验。

📒 博客介绍:分享嵌入式开发领域的相关知识、经验、思考和感悟,欢迎关注。提供嵌入式方向的学习指导、简历面试辅导、技术架构设计优化、开发外包等服务,有需要可加文末联系方式联系。

💬 博主粉丝群介绍:① 群内高中生、本科生、研究生、博士生遍布,可互相学习,交流困惑。② 热榜top10的常客也在群里,也有数不清的万粉大佬,可以交流写作技巧,上榜经验,涨粉秘籍。③ 群内也有职场精英,大厂大佬,可交流技术、面试、找工作的经验。④ 进群免费赠送写作秘籍一份,助你由写作小白晋升为创作大佬。⑤ 进群赠送CSDN评论防封脚本,送真活跃粉丝,助你提升文章热度。有兴趣的加文末联系方式,备注自己的CSDN昵称,拉你进群,互相学习共同进步。

NumPy线性代数模块详解:掌握numpy.linalg的核心功能

  • 1. NumPy库介绍
  • 2. linalg 模块介绍
    • 常用的 `numpy.linalg` 函数概述:
    • 2.1 函数定义及参数说明
      • linalg.inv
      • linalg.det
      • linalg.eig
      • linalg.solve
      • linalg.norm
  • 3. 示例代码
    • 3.1 计算矩阵的逆矩阵
    • 3.2 计算矩阵的行列式
    • 3.3 计算矩阵的特征值与特征向量
    • 3.4 解决线性方程组
    • 3.5 计算矩阵或向量的范数
  • 4. 实际应用:主成分分析(PCA)
  • 5. 总结

在这里插入图片描述

1. NumPy库介绍

NumPy(Numerical Python)是Python编程语言的一个核心库,用于大量的科学计算。 NumPy提供了对大型、多维数组和矩阵的支持,并且附带了大量的数学函数库来进行这些数组的操作。它是许多高级数据分析和机器学习库的基础,比如Pandas、SciPy和Scikit-learn。

NumPy的主要优势在于其数组对象(ndarray),这种对象比Python列表更为高效,可以存储同类型的数据元素,并且支持各种复杂的数值运算。对于需要进行大量数值计算和数据处理的应用程序,NumPy是首选工具。

2. linalg 模块介绍

numpy.linalg 模块提供了一组用于线性代数的基础函数。这些函数涵盖了矩阵分解、矩阵特征值与特征向量、求解线性系统等操作。线性代数是科学计算中一个重要的部分,NumPy通过numpy.linalg模块为用户提供高效且功能齐全的线性代数工具。

常用的 numpy.linalg 函数概述:

  • linalg.inv: 计算矩阵的逆矩阵。
  • linalg.det: 计算矩阵的行列式。
  • linalg.eig: 计算矩阵的特征值与特征向量。
  • linalg.solve: 解决线性方程组。
  • linalg.norm: 计算矩阵或向量的范数。

2.1 函数定义及参数说明

linalg.inv

计算逆矩阵。

numpy.linalg.inv(a)

参数:

  • a: 输入方阵。

返回:

  • out: 输入矩阵的逆矩阵。

linalg.det

计算矩阵的行列式。

numpy.linalg.det(a)

参数:

  • a: 输入方阵。

返回:

  • 行列式的值。

linalg.eig

计算矩阵的特征值与特征向量。

numpy.linalg.eig(a)

参数:

  • a: 输入方阵。

返回:

  • w: 特征值数组。
  • v: 特征向量构成的二维数组。

linalg.solve

解决线性方程组。

numpy.linalg.solve(a, b)

参数:

  • a: 系数矩阵。
  • b: 目标矩阵(或向量)。

返回:

  • 解向量或矩阵。

linalg.norm

计算矩阵或向量的范数。

numpy.linalg.norm(x, ord=None, axis=None, keepdims=False)

参数:

  • x: 输入数组。
  • ord: 范数类型(默认为2范数)。
  • axis: 计算范数的维度。
  • keepdims: 布尔值,是否保持原数组的维度。

返回:

  • 范数值。

3. 示例代码

接下来我们通过一些示例代码来展示numpy.linalg模块的具体用法。

3.1 计算矩阵的逆矩阵

在这个示例中,我们将展示如何计算一个方阵的逆矩阵。

import numpy as np

# 创建一个二维数组表示矩阵
A = np.array([[1, 2],
              [3, 4]])

# 计算矩阵的逆矩阵
A_inv = np.linalg.inv(A)
print("Inverse of A:\n", A_inv)

输出如下:

Inverse of A:
 [[-2.   1. ]
  [ 1.5 -0.5]]

3.2 计算矩阵的行列式

行列式是矩阵的重要属性之一,尤其在求解线性方程和矩阵特征值时起重要作用。

import numpy as np

# 创建一个二维数组表示矩阵
A = np.array([[1, 2],
              [3, 4]])

# 计算矩阵的行列式
det_A = np.linalg.det(A)
print("Determinant of A:", det_A)

输出如下:

Determinant of A: -2.0000000000000004

3.3 计算矩阵的特征值与特征向量

特征值与特征向量在很多领域有应用,如振动分析、图像处理和物理学。

import numpy as np

# 创建一个二维数组表示矩阵
A = np.array([[1, 2],
              [2, 1]])

# 计算矩阵的特征值与特征向量
eigenvalues, eigenvectors = np.linalg.eig(A)
print("Eigenvalues of A:", eigenvalues)
print("Eigenvectors of A:\n", eigenvectors)

输出如下:

Eigenvalues of A: [ 3. -1.]
Eigenvectors of A:
 [[ 0.70710678 -0.70710678]
 [ 0.70710678  0.70710678]]

3.4 解决线性方程组

求解形如 (Ax = b) 的线性方程组。

import numpy as np

# 创建系数矩阵A和目标向量b
A = np.array([[3, 1],
              [1, 2]])
b = np.array([9, 8])

# 解决线性方程组
x = np.linalg.solve(A, b)
print("Solution x:", x)

输出如下:

Solution x: [2. 3.]

3.5 计算矩阵或向量的范数

范数是衡量矩阵或向量大小的一种方式。

import numpy as np

# 创建一个二维数组表示矩阵
A = np.array([[1, 2],
              [3, 4]])

# 计算矩阵的Frobenius范数
norm_A = np.linalg.norm(A)
print("Frobenius norm of A:", norm_A)

# 创建一个一维数组表示向量
v = np.array([1, 2, 3])

# 计算向量的2范数(欧几里得范数)
norm_v = np.linalg.norm(v)
print("2-norm of v:", norm_v)

输出如下:

Frobenius norm of A: 5.477225575051661
2-norm of v: 3.7416573867739413

4. 实际应用:主成分分析(PCA)

主成分分析(PCA)是数据降维的经典方法。这里我们展示如何使用 numpy.linalg 进行PCA实现。

import numpy as np

# 创建一个示例数据集
X = np.array([[2.5, 2.4], 
              [0.5, 0.7], 
              [2.2, 2.9], 
              [1.9, 2.2], 
              [3.1, 3.0], 
              [2.3, 2.7], 
              [2, 1.6], 
              [1, 1.1], 
              [1.5, 1.6], 
              [1.1, 0.9]])

# 减去数据的均值
X_mean = X - np.mean(X, axis=0)

# 计算协方差矩阵
cov_matrix = np.cov(X_mean, rowvar=False)

# 计算协方差矩阵的特征值与特征向量
eigenvalues, eigenvectors = np.linalg.eig(cov_matrix)

# 将特征值排序,并获得相应的特征向量
sorted_index = np.argsort(eigenvalues)[::-1]
sorted_eigenvalues = eigenvalues[sorted_index]
sorted_eigenvectors = eigenvectors[:, sorted_index]

# 选择前两个主成分
n_components = 2
eigenvector_subset = sorted_eigenvectors[:, 0:n_components]

# 将数据投影到主成分空间
X_reduced = np.dot(eigenvector_subset.transpose(), X_mean.transpose()).transpose()

print("Reduced data:\n", X_reduced)

输出如下:

Reduced data:
 [[ 0.82797019  0.17511531]
 [-1.77758033  0.14285723]
 [ 0.99219749  0.38437499]
 [ 0.27421042  0.13041721]
 [ 1.67580142 -0.20949846]
 [ 0.9129491   0.17528244]
 [ 0.09910944 -0.3498247 ]
 [-1.14457216  0.04641726]
 [-0.43804614  0.01776463]
 [-1.40196572 -0.384375  ]]

5. 总结

NumPy是进行科学计算的强大工具,该库的numpy.linalg模块则专门提供了各种线性代数运算的支持。在这篇文章中,我们详细介绍了numpy.linalg中的一些常用函数,并通过丰富的示例演示了它们的具体应用。

通过这些示例,我们了解到:

  1. numpy.linalg.inv 用于计算矩阵的逆矩阵。
  2. numpy.linalg.det 用于计算矩阵的行列式。
  3. numpy.linalg.eig 用于计算矩阵的特征值与特征向量。
  4. numpy.linalg.solve 用于解决线性方程组。
  5. numpy.linalg.norm 用于计算矩阵或向量的范数。

此外,我们还展示了numpy.linalg在主成分分析(PCA)中的实际应用,展示了如何使用这些线性代数函数来进行数据降维操作。

掌握NumPy中的linalg模块,不仅能够有效完成各种线性代数计算,还能在实际的数据处理和分析任务中提供重要支持。如果你对更多的NumPy功能感兴趣,建议继续深入学习和探索。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1706175.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

软件测试人员如何规划自己的职业发展路线

在这个飞速发展的时代中,我们每个人都渴望在各自的领域里找到属于自己的一片天空,而对于那些致力于软件测试的朋友们而言,规划好自己的职业发展路线显得尤为重要。 一、明确目标 首先,我们需要确立一个明确的职业发展目标。对于软…

C# yolov8 TensorRT Demo

C# yolov8 TensorRT Demo 目录 效果 说明 项目 代码 下载 效果 说明 环境 NVIDIA GeForce RTX 4060 Laptop GPU cuda12.1cudnn 8.8.1TensorRT-8.6.1.6 版本和我不一致的需要重新编译TensorRtExtern.dll,TensorRtExtern源码地址:https://githu…

《爷爷的信》获短片金奖:电影新星的摇篮与短片的春天

随着上汽大众杯澳涞坞全球青年电影短片大赛金奖的公布,我们迎来了短片创作的一个崭新里程碑。这一赛事不仅为青年电影人搭建了一个展示才华的舞台,更预示着短片艺术即将迈入一个繁荣的新纪元。 澳涞坞集团此次大赛的举办,是对青年创作力量的…

美国西储大学(CRWU)轴承故障诊断——连续小波(CWT)变换

1.数据集介绍 2.代码 import random import matplotlib matplotlib.use(Agg) from scipy.io import loadmat import numpy as npdef split(DATA):step = 400;size = 1024;data = []for i in range(1, len(DATA) - size, step):data1 = DATA[i:i + size]data.append(data1)rand…

【NumPy】深入理解NumPy的cov函数:计算协方差矩阵的完整指南

🧑 博主简介:阿里巴巴嵌入式技术专家,深耕嵌入式人工智能领域,具备多年的嵌入式硬件产品研发管理经验。 📒 博客介绍:分享嵌入式开发领域的相关知识、经验、思考和感悟,欢迎关注。提供嵌入式方向…

全网首发!精选32个最新计算机毕设实战项目(附源码),拿走就用!

Hi 大家好,马上毕业季又要开始了,陆陆续续又要准备毕业设计了,有些学生轻而易举就搞定了,有些学生压根没有思路怎么做,可能是因为技术问题,也可能是因为经验问题。 计算机毕业相关的设计最近几年类型比较多…

基于半花青的近红外荧光材料

一、基于半花青的酯酶响应的光声探针: 参考文献:De Novo Design of Activatable Photoacoustic/Fluorescent Probes for Imaging Acute Lung Injury In Vivo | Analytical Chemistry (acs.org) 1.季铵盐的结构改造: 之前的基于半花青的近红外…

Linux网络_网络基础预备

文章目录 前言一、网络基础知识网络协议协议分层OSI七层模型TCP/IP五层(或四层)模型 认识IP地址认识MAC地址数据包封装和分用 前言 Linux系统编程已经告一段落,但是我们在学习LInux系统编程所积累的知识,将仍然与后面网络知识强相关,学习网络…

Nature plants|做完单细胞还可以做哪些下游验证实验

中国科学院分子植物科学中心与南方科技大学在《Nature Plants》期刊上(IF18.0)发表了关于苜蓿根瘤共生感知和早期反应的文章,该研究首次在单细胞水平解析了结瘤因子处理蒺藜苜蓿(Medicago truncatula)根系24小时内特异细胞类型的基因表达变化…

c++(四)

c(四) 运算符重载可重载的运算符不可重载的运算符运算符重载的格式运算符重载的方式友元函数进行运算符重载成员函数进行运算符重载 模板定义的格式函数模板类模板 标准模板库vector向量容器STL中的listmap向量容器 运算符重载 运算符相似,运…

上周暗网0day售卖情报一览

黑客声称以 1,700,000 美元出售 Outlook RCE 漏洞 0Day 令人担忧的是,一个名为“Cvsp”的威胁参与者宣布出售所谓的 Outlook 远程代码执行 (RCE) 漏洞 0day。这一所谓的漏洞旨在针对跨 x86 和 x64 架构的各种 Microsoft Office 版本,对全球用户构成重大安…

Facebook之魅:数字社交的体验

在当今数字化时代,Facebook作为全球最大的社交平台之一,承载着数十亿用户的社交需求和期待。它不仅仅是一个简单的网站或应用程序,更是一个将世界各地的人们连接在一起的社交网络,为用户提供了丰富多彩、无与伦比的数字社交体验。…

什么是NAND Flash ECC?

在存储芯片行业,数据完整性和可靠性是至关重要的。为了确保数据的准确性和防止数据丢失,ECC(错误校正码)在NAND Flash存储中扮演了关键角色。MK米客方德将为您解答NAND Flash ECC的基本概念、工作原理及其在实际应用中的重要性。 …

RGB 平均值统计

任务:有一一对应的图片多组如下,希望统计灰色部分原有grb平均值,彩色部分rgb平均值。 方法:由下图对各个像素分析,分为3类,并记录坐标,根据坐标统计上图的rgb平均值,结果放在一张Exc…

群晖异地组网-节点小宝搭建使用指南(全平台异地组网)

内网穿透,对于经常传输小文件、远程控制NAS的朋友来说是够用了,但是对经常异地端到端大文件传输需求的朋友来说就差点事,有没有一种免费、速度快、配置简单得方式的呢,答案是有的。节点小宝异地组网是一个非常不错的方式&#xff…

表空间[MAIN]处于脱机状态

达梦数据库还原后&#xff0c;访问数据库报错&#xff1a;表空间[MAIN]处于脱机状态 解决方法&#xff1a; 1&#xff1a;检查备份文件 DMRMAN 中使用 CHECK 命令对备份集进行校验&#xff0c;校验备份集是否存在及合法。 ##语法&#xff1a;CHECK BACKUPSET <备份集目录…

大字体学生出勤记录系统网页HTML源码

源码介绍 上课需要一个个点名记录出勤情况&#xff0c;就借助AI制作了一个网页版学生出勤记录系统&#xff0c; 大字体显示学生姓名和照片&#xff0c;让坐在最后排学生也能看清楚&#xff0c;显示姓名同时会语音播报姓名&#xff0c; 操作很简单&#xff0c;先导入学生姓名…

信息抽取模型TPLinker

1.motivation 早期传统方法首先抽取实体再抽取它们之间的关系&#xff0c;但是忽略了两个任务之间的关联。而后期采取的联合模型都存在着一个严重问题&#xff1a;训练时&#xff0c;真实值作为上下文传入训练&#xff1b;推理时&#xff0c;模型自身生成的值作为上下文传入&a…

代码随想录算法训练营第21天|● 530.二叉搜索树的最小绝对差 ● 501.二叉搜索树中的众数 ● 236. 二叉树的最近公共祖先

二叉搜索树的最小绝对差 题目连接 https://leetcode.cn/problems/minimum-absolute-difference-in-bst/ 思路&#xff1a; 利用二叉搜索树的中序遍历的特性&#xff0c;将二叉树转成有序数组&#xff0c;进而求任意两个数的最小绝对差。 代码 /*** Definition for a bina…

邮箱调用接口的服务有哪些?怎么配置接口?

邮箱调用接口安全性如何保障&#xff1f;使用邮箱服务器的方法&#xff1f; 邮箱调用接口为各种应用和系统提供了便捷的电子邮件发送与接收功能。选择合适的邮箱调用接口服务可以大大提升工作效率和用户体验。本AokSend将探讨一些主要的邮箱调用接口服务。 邮箱调用接口&…