Python这一编程语言已经问世将近30年了,在这期间,Python不仅完成了自身的进化,还获得了大量的用户。现在Python作为最具人气的编程语言,受到了许多人的喜爱。
接下来我们将使用Python实现深度学习系统。不过在这之前,本文将简单地介绍一下Python,看一看它的使用方法。
Python是一个简单
、易读
、易记
的编程语言,而且是开源
的,可以免费
地自由使用。Python可以用类似英语的语法编写程序,编译起来也不费力,因此我们可以很轻松的使用Python。特别是对首次接触Python的编程人士来说,Python是最合适不过的语言。事实上,很多高校和大专院校的计算机课程均采用Python作为入门语言。
此外,使用Python不仅可以写出可读性高
的代码,还可以写出性能高(处理速度快
?)的代码。在需要处理大规模的数据或者要求快速响应的情况下,使用Python可以稳妥地完成。因此,Python不仅受到初学者的喜爱,同时也受到专业人士的喜爱。实际上,Google、Microsoft、Facebook等战斗在IT行业最前沿的企业也经常使用Python。
综上,Python是最适合数据科学领域的编程语言。而且,Python具有受众广的优秀品质,从初学者到专业人士都在使用。因此,为了完成从零开始实现深度学习的目标,Python可以说是最适合的工具。
下面,我们首先要将Python安装到当前环境(电脑)上。这里说明一下安装时需要注意的一些地方。
Python有Python2和Python3两个版本。如果我们调查一下目前Python的使用情况,会发现除了最新的版本3.x以外,旧的版本2.x仍在被大量使用。因此,在安装Python时,需要慎重选择安装Python的哪个版本。这是因为两个版本之间没有兼容性(严格地讲,是没有“向后兼容性”),也就是说,会发生用Python3.x写的代码不能被Python2.x执行的情况。我们将使用Python3.x,只安装了Python2.x的读者建议另外安装一下Python3.x。
注:在本博客编写的时候,官方已停止对Python2.x的支持。
我们的目标是从零开始实现深度学习。因此,除了Numpy
库和Matplotlib
库之外,极力避免使用外部库。之所以使用这两个库,是因为它们可以有效的促进深度学习的实现。
Numpy
是用于数值计算的库,提供了很多高级的数学算法和便利的数组(矩阵)操作方法。本书中将使用这些便利的方法来有效促进深度学习的实现。
Matplotlib
是用来画图的库。使用Matplotlib能将实验结果可视化,并在视觉上确认深度学习运行期间的数据。
Python的安装方法有很多种,我们推荐使用Anaconda这个发行版。发行版集成了必要的库,使用户可以一次性完成安装。Anaconda是一个侧重于数据分析的发行版,前面说的Numpy
、Matplotlib
等有助于数据分析的库都包含在其中。
如前所述,接下来我们将使用Python3.x版本,因此Anaconda发行版也要安装3.x的版本。请读者从官网下载与自己的操作系统相应的发行版,然后安装。
完成Python的安装后,要先确认一下Python的版本。打开命令行窗口,输入python --version
命令,该命令会输出已经安装的Python版本的信息。
$ python --version
Python 3.8.5
如上所示,显示了Python3.8.5(根据实际安装的版本,版本号可能不同),说明已正确安装了Python3.x。接着输入Python,启动Python解释器。
$ python
Python 3.8.5 (default, Sep 3 2020, 21:29:08) [MSC v.1916 64 bit (AMD64)] :: Anaconda, Inc. on win32
Type "help", "copyright", "credits" or "license" for more information.
>>>
Python解释器也被称为“对话模式”,用户能够以和Python对话的方式进行编程。比如,当用户查询“1+2等于几?”的时候,Python解释器会回答“3”,所谓对话模式,就是指这样的交互
。现在,我们实际输入一下看看。
>>> 1 + 2
3
Python解释器可以像这样进行交互式
的编程。下面,我们使用交互模式,来看几个简单的Python编程的例子。
>>> 1 + 2
3
>>> 1 - 2
-1
>>> 4 * 5
20
>>> 7 / 5
1.4
>>> 3 ** 2
9
*
表示乘法,/
表示除法,**
表示乘方(3 ** 2是3的2次方)。另外,在Python2.x中,整数除以整数的结果是整数,比如7 / 5
的结果是1。但在Python3.x中,整数除以整数的结果是小数(浮点数)。
编程中有数据类型
(data type)这一概念。数据类型表示数据的性质,有整数、小数、字符串等类型。Python中的type()
函数可以用来查看数据类型。
>>> type(10)
<class 'int'>
>>> type(2.7)
<class 'float'>
>>> type("hello")
<class 'str'>
根据上面的结果可知,10是int类型
(整型),2.7是float类型
(浮点型),hello是str类型
(字符串型)。另外,“类型”和“类”这两个词有时用作相同的意思。这里,对于输出结果
可以使用x或y等字母定义变量
(variable)。此外,可以使用变量进行计算,也可以对变量赋值。
>>> x = 10 # 初始化
>>> print(x)
10
>>> x = 100 # 赋值
>>> print(x)
100
>>> y = 3.14
>>> x * y
314.0
>>> type(x * y)
<class 'float'>
Python是属于“动态类型语言
”的编程语言,所谓动态,是指变量的类型是根据情况自动决定的。在上面的例子中,用户并没有明确指出“x的类型是int(整型)”,是Python根据x被初始化为10,从而判断出x的类型为int的。此外,我们也可以看到,整数和小数相乘的结果是小数(数据类型的自动转换)。另外,#
是注释的意思,它后面的文字会被Python解释器忽略。
除了单一的数值,还可以用列表(数组)汇总数据。
>>> a = [1, 2, 3, 4, 5] # 生成列表
>>> print(a) # 输出列表的内容
[1, 2, 3, 4, 5]
>>> len(a) # 获取列表的长度
5
>>> a[0] # 访问第一个元素的值
1
>>> a[4]
5
>>> a[4] = 99 # 赋值
>>> print(a)
[1, 2, 3, 4, 99]
元素的访问是通过a[0]
这样的方式进行的。[]
中的数字称为索引(下标),索引从0开始(索引0对应第一个元素)。此外,Python的列表提供了切片
(slicing)这一便捷的标记法。使用切片不仅可以访问某个值,还可以访问列表的子列表(部分列表)。
>>> print(a)
[1, 2, 3, 4, 99]
>>> print(a)
[1, 2, 3, 4, 99]
>>> a[0:2] # 获取索引为0到2(不包括2!)的元素
[1, 2]
>>> a[1:] # 获取从索引为1的元素到最后一个元素
[2, 3, 4, 99]
>>> a[:3] # 获取从第一个元素到索引为3(不包括3!)的元素
[1, 2, 3]
>>> a[:-1] # 获取从第一个元素到最后一个元素的前一个元素之间的元素
[1, 2, 3, 4]
>>> a[:-2] # 获取从第一个元素到最后一个元素的前二个元素之间的元素
[1, 2, 3]
进行列表切片时,需要写成a[0:2]
这样的形式。a[0:2]
用于取出从索引为0的元素到索引为2的元素的前一个元素之间的元素。另外,索引-1
对应最后一个元素,-2
对应最后一个元素的前一个元素。
列表根据索引,按照0,1,2,···的顺序存储值,而字典则以键值对的形式存储数据。字典就像《新华字典》那样,将单词和它的含义对应着存储起来。
>>> me = {
"height" : 180} # 生成字典
>>> me["height"]
180
>>> me["weight"] = 70
>>> print(me)
{
'height': 180, 'weight': 70}
Python中bool
型。bool型取True
或False
中的一个值。针对bool型的运算符包括and
、or
和not
(针对数值的运算符有+
、-
、*
、/
等,根据不同的数据类型使用不同的运算符)。
>>> hungry = True # 饿了?
>>> sleepy = False # 困了?
>>> type(hungry)
<class 'bool'>
>>> not hungry
False
>>> hungry and sleepy # 饿并且困
False
>>> hungry or sleepy # 饿或者困
True
根据不同的条件选择不同的处理分支时可以使用if/else
语句。
>>> hungry = True
>>> if hungry:
... print("I'm hungry") # 使用空白字符进行缩进
...
I'm hungry
>>> hungry = False
>>> if hungry:
... print("I'm hungry") # 使用空白字符进行缩进
... else:
... print("I'm not hungry")
... print("I'm sleepy")
...
I'm hungry
I'm sleepy
Python中的空白字符具有重要的意义。上面的if语句中,if hungry:
下面的语句开头有4个空白字符。它是缩进的意思,表示当前面的条件成立时,此处的代码会被执行。这个缩进也可以用tab
表示,Python中推荐使用空白字符。(注:Python中使用空白字符表示缩进。一般而言,每缩进一次,使用4个空白字符。)
进行循环处理时,可以使用for
语句。
>>> for i in [1, 2, 3]:
... print(i)
...
1
2
3
这是输出列表[1, 2, 3]中元素的例子。使用for ··· in ··· :
语句结构,可以按顺序访问列表等数据集合中的各个元素。
可以将一连串的处理定义成函数
(function)。
>>> def hello():
... print("Hello World!")
...
>>> hello()
Hello World!
此外,函数可以取参数。
>>> def hello():
... print("Hello World!")
...
>>> hello()
Hello World!
>>> def hello(object):
... print("Hello " + object + "!")
...
>>> hello("cat")
Hello cat!
另外,字符串的拼接可以使用+
。
关闭Python解释器时,在命令行窗口中按Ctrl + Z
键即可。
到目前为止,我们看到的都是基于Python解释器的例子。Python解释器能够以对话模式执行程序,非常便于进行简单的实验。但是,想进行一连串的处理时,因为每次都需要输入程序,所以不太方便。这时,可以将Python程序保存为文件,然后(集中地)运行这个文件。下面,我们来看一个Python脚本文件的例子。
打开文本编辑器,新建一个hungry.py
的文件。hungry.py只包含下面一行语句。
print("I'm hungry!")
接着,打开命令行窗口,移至hungry.py文件所在的位置。然后,将hungry.py文件名为作为参数,运行python命令。
$ python hungry.py
I'm hungry!
这样,使用python hungry.py
命令就可以执行这个Python程序了。
前面我们了解了int
和str
等数据类型(通过type()
函数可以查看对象的类型)。这些数据类型是“内置”的数据类型,是Python中一开始就有的数据类型。现在,我们来定义新的类。如果用户自己定义类的话,就可以自己创建数据类型。此外,也可以定义原创的方法(类的函数)和属性。
Python中使用class
关键字来定义类,类要遵循下述格式(模板)。
class 类名:
def __init__(self, 参数, ...): # 构造函数
...
def 方法名1(self, 参数, ...): # 方法1
...
def 方法名2(self, 参数, ...): # 方法2
...
这里有一个特殊的__init__
方法,这是进行初始化的方法,也称为构造函数
(constructor),只在生成类的实例时被调用一次。此外,在方法的第一个参数中明确地写入表示自身(自身的实例)的self
是Python的一个特点(学过其他编程语言的人可能会觉得这种写self
的方式有一点奇怪)。
下面我们通过一个简单的例子来创建一个类。
class man:
def __init__(self, name):
self.name = name
print("Initialized!")
def hello(self):
print("Hello " + self.name + "!")
def goodbye(self):
print("Good-bye " + self.name + "!")
m = man("David") # Initialized!
m.hello() # Hello David!
m.goodbye() # Good-bye David!
这里我们定义了一个新类man。上面的例子中,类man生成了实例(对象)m。
类man的构造函数(初始化方法)会接收参数name,然后用这个参数初始化实例变量self.name
。实例变量是存储在各个实例中的变量。Python中可以像self.name这样,通过在self后面添加属性名来生成或访问实例变量。
在深度学习的实现中,经常出现数组和矩阵的计算。Numpy的数组类(numpy.ndarray
)提供了很多便捷的方法,在实现深度学习时,我们将使用这些方法。接下来,我们简单介绍一下后面会用到的Numpy。
Numpy是外部库。这里所说的“外部”是指不包含在标准版Python中。因此,我们首先要导入Numpy库。
>>> import numpy as np
Python中使用import
语句来导入库。这里的import numpy as np
,直译的话就是“将Numpy作为np导入”的意思。通过写成这样的形式,之后Numpy相关的方法均可通过np来调用。
要生成Numpy数组,需要使用np.array()
方法。np.array()
接收Python列表作为参数,生成Numpy数组(numpy.ndarray)。
>>> x = np.array([1.0, 2.0, 3.0])
>>> print(x)
[1. 2. 3.]
>>> type(x)
<class 'numpy.ndarray'>
下面是Numpy数组的算术运算的例子。
>>> x = np.array([1.0, 2.0, 3.0])
>>> y = np.array([2.0, 4.0, 6.0])
>>> x + y # 对应元素的加法
array([3., 6., 9.])
>>> x - y
array([-1., -2., -3.])
>>> x * y
array([ 2., 8., 18.])
>>> x / y
array([0.5, 0.5, 0.5])
这里需要注意的是,数组x和数组y的元素个数是相同的(两者均是元素个数为3的一维数组)。当x和y的元素个数相同时,可以对各个元素进行算术运算。如果元素个数不同就会报错,所以元素个数保持一致非常重要。另外,“对应元素的”英文是element-wise,比如“对应元素的乘法”就是element-wise product。
Numpy数组不仅可以进行element-wise,比如“对应元素的惩罚”就是element-wise product。
Numpy数组不仅可以进行element-wise元素,也可以和单一的数值(标量)组合起来进行运算。此时,需要在Numpy数组的各个元素和标量之间进行运算。这个功能也被称为广播
。
>>> x = np.array([1.0, 2.0, 3.0])
>>> x / 2.0
array([0.5, 1. , 1.5])
Numpy不仅可以生成一维数组(排成一列的数组),也可以生成多维数组。比如,可以生成如下的二位数组(矩阵)。
>>> A = np.array([[1, 2], [3, 4]])
>>> print(A)
[[1 2]
[3 4]]
>>> A.shape
(2, 2)
>>> A.dtype
dtype('int32')
这里生成了一个2x2的矩阵A。另外,矩阵A的形状可以通过shape
查看,矩阵元素的数据类型可以通过dtype
查看。下面,我们来看一下矩阵的算术运算。
>>> B = np.array([[3, 0], [0, 6]])
>>> A + B
array([[ 4, 2],
[ 3, 10]])
>>> A.shape
(2, 2)
>>> A.dtype
dtype('int32')
和数组的算术运算一样,矩阵的算术运算也可以在相同形状的矩阵间以对应元素的方式进行。并且,也可以通过标量(单一数值)对矩阵进行算术运算。这也是基于广播的功能。
>>> print(A)
[[1 2]
[3 4]]
>>> A * 10
array([[10, 20],
[30, 40]])
Numpy数组(np.ndarray)可以生成N维数组,既可以生成一维数组、二维数组、三维数组等任意维数的数组。数学上将一维数组称为向量,将二维数组称为矩阵。另外,可以将一般化之后的向量或矩阵等统称为张量(tensor)。本书基本上将二维数组称为“矩阵”,将三维数组及三维以上的数组称为“张量”或“多维数组”。
Numpy中,形状不同的数组之间也可以进行计算。之前的例子中,在2x2的矩阵A和标量10之间进行了乘法运算。在这个过程中,标量10被扩展成了2x2的形状,然后再与矩阵A进行乘法运算。这个巧妙地功能称为广播
(broadcast)。
我们通过下面这个运算再来看一下广播的例子。
>>> A = np.array([[1, 2], [3, 4]])
>>> B = np.array([10, 20])
>>> A * B
array([[10, 40],
[30, 80]])
在这个运算中,一维数组B被“巧妙地”变成了和二维数组A相同的形状,然后再以对应元素的方式进行运算。
综上,因为Numpy有广播功能,所以不同形状的数组之间也可以顺利地进行运算。
元素地索引从0开始。对各个元素的访问可按如下方式进行。
>>> X = np.array([[51, 55], [14, 19], [0, 4]])
>>> print(X)
[[51 55]
[14 19]
[ 0 4]]
>>> X[0] # 第0行
array([51, 55])
>>> X[0][1] # (0, 1)的元素
55
也可以使用for语句访问各个元素。
>>> for row in X:
... print(row)
...
[51 55]
[14 19]
[0 4]
除了前面介绍的索引,Numpy还可以使用数组访问各个元素。
>>> X = X.flatten() # 将X转换为一维数组
>>> print(X)
[51 55 14 19 0 4]
>>> X[np.array([0, 2, 4])] # 获取索引为0、2、4的元素
array([51, 14, 0])
运用这个标记法,可以获取满足一定条件的元素。例如,要从X中抽出大于15的元素,可以写成如下形式。
>>> X > 15
array([ True, True, False, True, False, False])
>>> X[X > 15]
array([51, 55, 19])
对Numpy数组使用使用不等号运算符(上例中是X > 15
),结果会得到一个布尔型的数组。上例中就是使用这个布尔型数组取出了数组的各个元素(取出True对应的元素)。
注:Python等动态类型语言一般比C和C++等静态类型语言(编译型语言)运算速度慢。实际上,如果是运算量大的处理对象,用C/C++写程序更好。为此,当在Python中追求性能时,人们会用C/C++来实现处理的内容。Python则承担“中间人”的角色,负责调用那些用C/C++写的程序。Numpy中,主要的处理也都通过C或C++实现的。因此,我们可以在不损失性能的情况下,使用Python便利的语法。
在深度学习实验中,图形的绘制和数据可视化非常重要。Matplotlib是用于绘制图形的库,使用Matplotlib可以轻松地绘制图形和实现数据的可视化。这里,我们来介绍一下图形的绘制方法和图像的显示方法。
可以使用matplotlib的pyplot模块绘制图形。话不多说,我们来看一个绘制sin
函数曲线的例子。
import numpy as np
import matplotlib.pyplot as plt
# 生成数据
x = np.arange(0, 6, 0.1)
y = np.sin(x)
# 绘制图形
plt.plot(x, y)
plt.show()
这里使用Numpy的arange方法生成了[0, 0.1, 0.2, …, 5.8, 5.9]的数据,将其设为x。对x的各个元素,应用Numpy的sin函数np.sin()
,将x,y的数据传给plt.plot方法,然后绘制图形。最后,通过plt.show()
显示图形。
在刚才的sin函数的图形中,我们尝试追加cos函数的图形,并尝试使用pyplot的添加标题和x轴标签名等其他功能。
import numpy as np
import matplotlib.pyplot as plt
# 生成数据
x = np.arange(0, 6, 0.1) # 以0.1为单位,生成0到6的数据
y1 = np.sin(x)
y2 = np.cos(x)
# 绘制图形
plt.plot(x, y1, label="sin")
plt.plot(x, y2, linestyle="--", label="cos") # 用虚线绘制
plt.xlabel("x") # x轴标签
plt.ylabel("y") # y轴标签
plt.title("sin & cos") # 标题
plt.legend()
plt.show()
pyplot中还提供了用于显示图像的方法imshow()
。另外,可以使用matplotlib.image
模块的imread()
方法读入图像。下面我们来看一个例子。
import matplotlib.pyplot as plt
from matplotlib.image import imread
img = imread("image.jpg") # 读入图像(设定合适的路径!)
plt.imshow(img)
plt.show()
这里,我们假定图像image.jpg
在当前目录下。读者根据自己的环境,可能需要变更文件名或文件路径。
本章重点介绍了实现深度学习(神经网络)所需的编程知识,以为学习深度学习做好准备。从下一章开始,我们将通过使用Python实际运行代码,逐步了解深度学习。
上次更新时间:2021.02.02