题目## 题目
向量到直线的正交投影是一种将向量投影到直线上的方法,其计算公式为:
P
=
a
⋅
b
a
⋅
a
a
P = \frac{a \cdot b}{a \cdot a} a
P=a⋅aa⋅ba
其中,
a
a
a 是直线的方向向量,
b
b
b 是向量。
标准代码如下
def dot(v1, v2):
return sum([ax1 * ax2 for ax1, ax2 in zip(v1, v2)])
def scalar_mult(scalar, v):
return [scalar * ax for ax in v]
def orthogonal_projection(v, L):
L_mag_sq = dot(L, L)
proj_scalar = dot(v, L) / L_mag_sq
proj_v = scalar_mult(proj_scalar, L)
return [round(x, 3) for x in proj_v]
除此之外,还可以使用numpy库的函数来实现向量到直线的正交投影。
def orthogonal_projection(v, L):
import numpy as np
v = np.array(v)
L = np.array(L)
return np.round((L * np.dot(v, L) / np.dot(L, L)),3).tolist()
压缩行稀疏矩阵(CSR)格式是一种特殊的矩阵存储格式,其特点是只存储非零元素的值、行号和列指针。本质上是一种三元组表示法。
标准代码如下
def compressed_row_sparse_matrix(dense_matrix):
vals = []
col_idx = []
row_ptr = [0]
for row in dense_matrix:
for j, val in enumerate(row):
if val != 0:
vals.append(val)
col_idx.append(j)
row_ptr.append(len(vals))
return vals, col_idx, row_ptr
除此之外,还可以使用scipy库中的csr_matrix函数来实现压缩行稀疏矩阵的存储。
def compressed_row_sparse_matrix(dense_matrix):
from scipy.sparse import csr_matrix
sparse = csr_matrix(dense_matrix)
return sparse.data.tolist(), sparse.indices.tolist(), sparse.indptr.tolist()
313

被折叠的 条评论
为什么被折叠?



