# 01、Python科学计算和数据分析库简介
Python通过众多出色的第三方库,为科学计算和数据分析提供了支持。主要包括:
(1)NumPy。NumPy是Python数据处理的底层库,是高性能科学计算和数据分析的基础,许多其它科学计算库都基于NumPy库。
(2)SciPy。SciPy是Python解决科学计算中常见问题的工具集团,其不同子模块对应于不同的应用,包括:线性代数、插值、积分、最优化、图像处理、统计等等。
(3)pandas。Pandas是Python的高级数据分析工具库。
(4)Matplotlib。Matplotlib是Python最流行的绘图程序库。Seaborn和Prettyplotlib是在matplotlib的基础上改进功能的绘图程序库。
(5)SciKit。SciKits程序库为SciPy、NumPy和Python提供了专业化的扩展。例如,scikit-learn是用于机器学习和数据挖掘的库。
(6)IPython。IPython是基于Python Shell的交互式解释器,功能更加强大,被广泛应用于交互式科学计算。
# 02、使用IPython Notebook进行科学计算
在案例研究1中简单介绍了Python面向科学计算的环境Anaconda。安装Anaconda时会安装Jupyter Notebook(以前的版本也叫IPython Notebook)。
Jupyter Notebook是一种交互式的数据分析与记录工具,它既是一个交互计算平台,又是一个记录计算过程的“笔记本”,它是数据分析、科学计算以及交互计算的“利器”。
Jupyter Notebook的特点是支持可重复性的互动计算,即可以重复更改并且执行曾经的输入记录。它可以记录演算过程,并保存成其他很多格式,例如Python脚本、HTML、PDF等。很多课程、数据和博客是用Notebook写作的。
Jupyter Notebook使用浏览器作为界面,向后台的IPython服务器发送请求,并显示结果。在浏览器的界面中使用单元格(Cell)输入保存各种信息。Cell主要有两种类型,即Code(输入、编辑和执行Python代码)和MarkDown(输入、编辑和显示Markdown格式的文本)。
Jupyter Notebook由服务端和客户端两部分组成。服务端可以运行在本机,也可以运行在远程服务器上,它主要包含负责运算的iPython kernel,以及一个HTTP/HTTPS服务器(Tornado),主要负责代码的解释和计算;而客户端是浏览器,主要负责与用户进行交互,接受用户的输入,以及渲染输出。
【例1】Jupyter Notebook使用示例。
(1)执行Windows菜单命令“开始”|“所有应用”|Anaconda3(64-bit)|Jupyter Notebook,启动本地服务器,并在默认浏览器中打开主页。如图1所示。
#
■ 图1 启动Jupyter(IPython)Notebook
(2)新建一个Notebook文件。单击如图1中的“New”下拉菜单,选择“Python 3(ipykernel)”,新建一个名为Untitled的Notebook。
(3)输入Markdown文本。在第一个单元格输入框中,输入“使用IPython Notebook进行科学计算”,并设置单元格的类型为Markdown。如图2所示。
■ 图2 输入Markdown文本
(4)输入并执行Python代码。使用菜单命令Insert|Insert Cell Below或者快捷键B,在下方插入一个单元格,输入图3中所示的代码,按快捷键Ctrl+Enter,执行并显示结果。
■ 图3 输入并执行Python代码
03、Python在线性代数中的应用
线性代数(linear algebra)是数学的一个分支,主要研究矩阵理论、向量空间、线性变换和有限维线性方程组等内容。使用Python的科学计算库NumPy和SciPy,可以进行线性代数的学习和研究。
3.1 创建向量和矩阵(数组)
通常使用NumPy的array函数可以创建向量和矩阵(数组),一维数组表示向量、二维数组表示矩阵,也可以使用函数zeros()、zeros_like()、ones()、ones_like()、empty()、empty_like()、arange()、linspace()、numpy.random.rand()、numpy.random.randn()、fromfunction()、fromfile()等方便创建特殊数组。
【例2】创建向量和矩阵(数组)示例如图4所示。
■ 图4 创建向量和矩阵(数组)示例
3.2 矩阵的运算
使用NumPy包,矩阵的加法和减法可以直接使用运算符+和-,乘法则使用函数dot()。
【例3】矩阵的运算。
3.3 线性方程组的求解
已知线性方程组AX=b,则X=A-1b。例如:
x + 3y + 5z = 10
2x + 5y + z = 8
2x + 3y + 8z = 3
则:
【例4】线性方程组求解。