Python NumPy入门指南:数据处理科学计算的瑞士军刀

Python NumPy入门指南:数据处理科学计算的瑞士军刀
作者:唐叔在学习
专栏:唐叔学python
标签:Python NumPy、数据分析、科学计算、机器学习基础、数组操作、Python数据处理、人工智能基础、Python编程

摘要

NumPy是Python科学计算的基础库,提供了高性能的多维数组对象和工具。本文唐叔将带你从零开始了解NumPy的核心概念、常用操作和实际应用场景,助你在数据分析、机器学习等领域快速上手。无论你是Python初学者还是想提升数据处理能力,这篇文章都将成为你的实用指南。

文章目录


一、NumPy是什么?为什么它如此重要?

各位小伙伴们好,我是唐叔!今天我们要聊的是Python科学计算领域的一个重量级选手——NumPy。

简单来说,NumPy(Numerical Python)是Python中用于科学计算的基础包。它提供了一个强大的N维数组对象,以及用于处理这些数组的各种工具。可以说,NumPy是Python数据科学生态系统的基石,几乎所有的机器学习框架(如TensorFlow、PyTorch)和数据分析库(如Pandas)都建立在NumPy之上。

为什么NumPy这么受欢迎?主要有三大原因:

  1. 性能强大:NumPy的底层是用C语言编写的,运算速度比纯Python快几个数量级
  2. 接口简单:提供了大量简洁高效的数组操作函数
  3. 生态完善:与SciPy、Matplotlib、Pandas等库完美配合

二、NumPy安装与基础使用

2.1 安装NumPy

安装NumPy非常简单,使用pip命令即可:

pip install numpy 

如果你使用的是Anaconda,它已经自带了NumPy,无需额外安装。

2.2 导入NumPy

按照惯例,我们通常这样导入NumPy:

import numpy as np 

这样我们就可以用np作为前缀来调用NumPy的各种功能了。

2.3 创建第一个NumPy数组

让我们创建一个简单的数组来感受一下:

import numpy as np # 创建一维数组 arr1 = np.array([1,2,3,4,5])print(arr1)# 创建二维数组 arr2 = np.array([[1,2,3],[4,5,6]])print(arr2)

输出结果:

[1 2 3 4 5] [[1 2 3] [4 5 6]] 

三、NumPy核心功能详解

3.1 数组属性

每个NumPy数组都有一些重要的属性:

arr = np.array([[1,2,3],[4,5,6]])print("数组维度:", arr.ndim)# 2print("数组形状:", arr.shape)# (2, 3)print("数组元素总数:", arr.size)# 6print("数组元素类型:", arr.dtype)# int32或int64

3.2 创建特殊数组

NumPy提供了多种创建特殊数组的方法:

# 创建全零数组 zeros = np.zeros((2,3))# 创建全1数组 ones = np.ones((3,2))# 创建单位矩阵 eye = np.eye(3)# 创建等差数列 lin = np.linspace(0,10,5)# [ 0. 2.5 5. 7.5 10. ]# 创建随机数组 rand = np.random.rand(2,2)
在这里插入图片描述

3.3 数组索引与切片

NumPy的索引和切片操作是数据处理的核心:

arr = np.array([[1,2,3],[4,5,6],[7,8,9]])# 获取单个元素print(arr[1,2])# 6# 获取一行print(arr[1])# [4 5 6]# 获取一列print(arr[:,1])# [2 5 8]# 切片print(arr[0:2,1:3])""" [[2 3] [5 6]] """

3.4 数组运算

NumPy支持各种数学运算:

a = np.array([1,2,3]) b = np.array([4,5,6])# 加法print(a + b)# [5 7 9]# 乘法(元素级)print(a * b)# [4 10 18]# 点积print(np.dot(a, b))# 32# 平方print(a**2)# [1 4 9]# 平方根print(np.sqrt(4))# 2.0# 求和print(a.sum())# 6.0# 均值print(a.mean())# 2.0# 三角函数print(np.sin(a))

常见的数组运算见下表:

运算符号 / 函数例子
逐元素加+a + b
逐元素乘*a * b
矩阵乘@np.dota @ b.T
平方根np.sqrt(a)逐元素
求和a.sum()a.sum(axis=0)按列求和
均值a.mean()

3.5 广播机制

NumPy的广播机制是其强大功能之一:

a = np.array([[1,2,3],[4,5,6]]) b = np.array([10,20,30])# b会被广播到与a相同的形状print(a + b)""" [[11 22 33] [14 25 36]] """

四、NumPy高级功能

4.1 数组变形

arr = np.arange(6)# [0 1 2 3 4 5]# 改变形状print(arr.reshape(2,3))""" [[0 1 2] [3 4 5]] """# 展平数组print(arr.flatten())# [0 1 2 3 4 5]

4.2 数组拼接与分割

a = np.array([[1,2],[3,4]]) b = np.array([[5,6]])# 垂直拼接print(np.vstack((a, b)))""" [[1 2] [3 4] [5 6]] """# 水平拼接 c = np.array([[7],[8]])print(np.hstack((a, c)))""" [[1 2 7] [3 4 8]] """# 数组分割 arr = np.arange(9).reshape(3,3)print(np.hsplit(arr,3))# 分成3个数组

4.3 统计函数

NumPy提供了丰富的统计函数:

data = np.array([[1,2,3],[4,5,6],[7,8,9]])print("平均值:", np.mean(data))print("中位数:", np.median(data))print("标准差:", np.std(data))print("最大值:", np.max(data))print("最小值:", np.min(data))print("每列的和:", np.sum(data, axis=0))print("每行的和:", np.sum(data, axis=1))
在这里插入图片描述

五、NumPy在实际项目中的应用

NumPy 在实际项目中,由于其内置的大量函数,使其在图像处理、机器学习数据预处理、数据计算处理等场景都有广泛使用。可以不夸张的说:搞AI、搞数据,NumPy是地基!

5.1 图像处理

NumPy数组可以用来表示图像:

from PIL import Image import matplotlib.pyplot as plt import numpy as np # 读取图像为NumPy数组 img = np.array(Image.open('image.png'))# # 显示图像# plt.imshow(img)# plt.show()## # 转换为灰度图# gray = np.mean(img, axis=2)# plt.imshow(gray, cmap='gray')# plt.show() fig,(ax1, ax2)= plt.subplots(1,2) ax1.imshow(img) gray = np.mean(img, axis=2) ax2.imshow(gray, cmap='gray') plt.show()
在这里插入图片描述

5.2 机器学习数据预处理

# 特征标准化 data = np.random.rand(100,5)# 100个样本,5个特征# 计算均值和标准差 mean = np.mean(data, axis=0) std = np.std(data, axis=0)# 标准化 normalized =(data - mean)/ std # 打乱数据顺序 np.random.shuffle(normalized)

5.3 数值计算

# 解线性方程组 A = np.array([[2,1],[1,3]]) b = np.array([4,5]) x = np.linalg.solve(A, b)# [1.4, 1.2]# 计算特征值和特征向量 eigenvalues, eigenvectors = np.linalg.eig(A)''' [[-0.85065081 -0.52573111] [ 0.52573111 -0.85065081]] '''

六、NumPy性能优化技巧

  1. 避免循环:尽量使用NumPy内置函数代替Python循环
  2. 使用视图而非副本:切片操作返回的是视图,不会复制数据
  3. 预分配内存:对于大型数组,先创建好再填充
  4. 使用适当的数据类型:如np.float32np.float64节省内存
import time import numpy as np n = 1_000_000 # 不好的做法 start = time.time() result =[]for i inrange(n): result.append(i**2) result = np.array(result)print("it cost: {}".format(time.time()- start))# 好的做法 start = time.time() result = np.arange(n)**2 time.time() result = np.arange(n)**2print("it cost: {}".format(time.time()- start))
在这里插入图片描述

七、总结

NumPy作为Python科学计算的基石,提供了:

  1. 高效的多维数组对象
  2. 丰富的数学函数库
  3. 强大的广播机制
  4. 便捷的线性代数运算
  5. 与其它科学计算库的良好兼容性

掌握了NumPy,你就打开了Python数据科学的大门。无论是数据分析、机器学习还是科学计算,NumPy都是不可或缺的工具。希望这篇文章能帮助你快速上手NumPy,在数据处理的道路上越走越远!

最后提醒:学习NumPy最好的方式就是多练习!打开你的Python环境,跟着示例代码敲一遍,感受NumPy的强大魅力吧!


觉得有用的话,别忘了点赞收藏加关注!

往期文章推荐:

Read more

【C++经典例题】字符串转整数(atoi)的实现与解析

【C++经典例题】字符串转整数(atoi)的实现与解析

💓 博客主页:倔强的石头的ZEEKLOG主页             📝Gitee主页:倔强的石头的gitee主页             ⏩ 文章专栏:C++经典例题                                   期待您的关注   目录 一、问题描述 二、解题思路 三、代码实现 四、代码逻辑详解 1. 变量初始化 2. 忽略前导空格 3. 处理符号 4. 转换数字 5. 返回结果     一、问题描述 LCR 192. 把字符串转换成整数 (atoi) - 力扣(LeetCode) 在编程中,经常会遇到将字符串转换为整数的需求,就像标准库中的 atoi 函数一样。 本题要求实现一个 myAtoi 函数,将输入的字符串转换为 32 位有符号整数,具体规则如下:   1. 读入字符串并丢弃无用的前导空格。

By Ne0inhk
Effective Modern C++ 条款37:使std::thread在所有路径最后都不可结合

Effective Modern C++ 条款37:使std::thread在所有路径最后都不可结合

Effective Modern C++ 条款37:使std::thread在所有路径最后都不可结合 * 引言:线程生命周期的关键问题 * 线程的两种状态:可结合与不可结合 * 可结合(Joinable)状态的特征 * 不可结合(Unjoinable)状态的四种情况 * 为什么可结合性如此重要? * 两种被拒绝的替代方案 * RAII拯救方案:ThreadRAII类 * ThreadRAII实现详解 * 关键设计决策 * 实际应用案例 * 高级讨论:何时选择join或detach * 性能考量与最佳实践 * 结论:让线程管理无忧 BiliBili上对应的视频为:https://www.bilibili.com/video/BV1iZZgBiE9j 引言:线程生命周期的关键问题 在多线程程序设计中,std::thread的管理是一个看似简单实则暗藏玄机的话题。想象一下,你精心设计的并发程序在大多数情况下运行良好,却在某些边缘情况下突然崩溃——这正是许多开发者在使用原生线程时遇到的噩梦场景。本文将深入探讨std::thread对象

By Ne0inhk
《C++ 动态规划》第001-002题:第N个泰波拉契数,三步问题

《C++ 动态规划》第001-002题:第N个泰波拉契数,三步问题

🔥个人主页:Cx330🌸 ❄️个人专栏:《C语言》《LeetCode刷题集》《数据结构-初阶》《C++知识分享》 《优选算法指南-必刷经典100题》《Linux操作系统》:从入门到入魔 《Git深度解析》:版本管理实战全解 🌟心向往之行必能至 🎥Cx330🌸的简介: 目录 前言: 01.第N个泰波拉契数 算法原理(动态规划): 思路: 解法代码(C++): 博主手记(字体还请见谅哈): 02.三步问题 算法原理(动态规划): 思路: 解法代码(C++): 博主手记(字体还请见谅哈): 结尾: 前言: 聚焦算法题实战,系统讲解三大核心板块:“精准定位最优解”——优选算法,“简化逻辑表达,系统性探索与剪枝优化”——递归与回溯,“以局部最优换全局高效”——贪心算法,讲解思路与代码实现,帮助大家快速提升代码能力 01.

By Ne0inhk
【C++】AVL树的底层以及实现

【C++】AVL树的底层以及实现

个人主页 文章目录 * ⭐一、AVL树的概念 * 🎉二、AVL树的性质 * 🏝️三、AVL树的实现 * 1. 树的基本结构 * 2. 树的插入 * 3. 树的旋转 * • 左单旋 * • 右单旋 * • 左右双旋 * • 右左双旋 * 🎡四、AVL树的其它功能 * 1. 树的查找 * 2. 树的遍历 * 3. 树的高度 * 4. 树的大小 * 🚀五、总结 * 1. AVL树的优缺点 * 2. 完整代码 ⭐一、AVL树的概念 AVL树是一种高度平衡的平衡二叉树,相比于搜索二叉树,它的特点在于左右子树都为AVL树且树的高度差的绝对值不超过1。 这里我们会引入一个新的概念叫做平衡因子。平衡因子也就是左右子树的高度差,我们可以通过平衡因子方便我们后续去观察和控制树是否平衡。 🎉二、AVL树的性质 AVL树主要有三大性质: 1.每棵树的左右子树都是AVL树。 2.左子树和右子树的高度之差的绝对值不超过1。 3.

By Ne0inhk