Python数值算法片段

集合${\mathbb{R}}^n$是所有$n$元组的实数的集合。 用集合表示法是 R n = { ( x 1 , x 2 , x 3 , ⋯   , x n ) : x 1 , x 2 , x 3 , ⋯   , x n ∈ R } \mathbb{R}^{n}=\left\{\left(x_{1}, x_{2}, x_{3}, \cdots, x_{n}\right): x_{1}, x_{2}, x_{3}, \cdots, x_{n} \in \mathbb{R}\right\} Rn={(x1​,x2​,x3​,⋯,xn​):x1​,x2​,x3​,⋯,xn​∈R}。 例如，集合${\mathbb{R}}^3$表示三维空间中的实三元组$(x,y,z)$坐标。

${\mathbb{R}}^n$中的向量是${\mathbb{R}}^n$中的$n$元组或点。 向量可以在行向量中水平地写（即向量的元素彼此相邻），或在列向量中垂直地写（即向量的元素彼此顶部）。 如果向量的上下文不明确，则通常意味着该向量是列向量。 向量的第$i$个元素$v$用$v_i$表示。 列向量的转置是相同长度的行向量，而行向量的转置是列向量。 在数学中，转置由上标$T$或$v^T$表示。 零向量是${\mathbb{R}}^n$中包含所有零的向量。

向量的范数是其长度的度量。 有许多方法可以根据所使用的度量（即，所选的距离公式）来定义向量的长度。 最常见的称为$L_2$范数，它是根据您可能从小学时熟悉的距离公式计算得出的。 向量$v$的$L_2$范数由$\Vert v{\Vert }_2$表示，$\Vert v{\Vert }_2=\sqrt{{\sum }_i{v}_i^2}$。 有时也称为欧几里得长度，是指向量在一维，二维或三维空间中的“物理”长度。$L_1$范数或“曼哈顿距离”的计算公式为$\Vert v{\Vert }_1={\sum }_i\left|v_i\right|$，以纽约市的网格状道路结构命名。通常，向量的$p$范数$L_p$为：$L_{\infty }$范数是$p$范数，其中$p=\infty$。$L_{\infty }$范数写为$\Vert v{\Vert }_{\infty }$，它等于$v$中的最大绝对值。

片段1

创建行向量和列向量，并显示向量的形状。

import numpy as np
vector_row = np.array([[1, -5, 3, 2, 4]])
vector_column = np.array([[1], 
                          [2], 
                          [3], 
                          [4]])
print(vector_row.shape)
print(vector_column.shape)

输出

(1, 5)
(4, 1)

注意！ 在Python中，行向量和列向量有些棘手。 从上面可以看到，为了获得1行4列或4行1列向量，我们必须使用list列表进行指定。 您可以定义np.array（[1,2,3,4]），但很快您会发现它不包含有关行或列的信息。

片段2

将上面定义的行向量转置为列向量，并计算其$L_1,L_2$和$L_{\infty }$范数。验证向量的$L_{\infty }$范数等于向量中元素的最大值。

from numpy.linalg import norm
new_vector = vector_row.T
print(new_vector)
norm_1 = norm(new_vector, 1)
norm_2 = norm(new_vector, 2)
norm_inf = norm(new_vector, np.inf)
print('L_1 is: %.1f'%norm_1)
print('L_2 is: %.1f'%norm_2)
print('L_inf is: %.1f'%norm_inf)

详情参阅 - 亚图跨际