numpy中的 dot, outer and *

每次用到都得先测试一下,整理一下以后会方便一些,先举几个例子吧:

demo1:

import numpy as np
a = [[1,2,3],[4,5,6]]
a = np.array(a)
b = [[1,2],[4,5],[3,6]]
b= np.array(b)
#一个(23), 一个(32),是可以正常的矩阵相乘的
c = np.dot(a, b)
print c
-->result: array([[18, 30],
              [42, 69]])
说明dot是正常的矩阵相乘的方法
c = a * b
-->result: error 不能相乘
注意:这里用的是array不是matrix, 如果定义为matrix则* 表示矩阵乘,是OK的。

c = np.outer(a,b)
-->result:
array([[ 1,  2,  4,  5,  3,  6],
       [ 2,  4,  8, 10,  6, 12],
       [ 3,  6, 12, 15,  9, 18],
       [ 4,  8, 16, 20, 12, 24],
       [ 5, 10, 20, 25, 15, 30],
       [ 6, 12, 24, 30, 18, 36]])
说明outer是a的第一个元素跟b的每一个元素相乘作为第一行,第二个元素跟b的每一个元素相乘作为第二个元素...
一半outer()用的很少

下面试试*的具体作用:

d = np.array([[1,2,1],[2,3,1]])
c = a*d
result:
array([[ 1,  4,  3],
       [ 8, 15,  6]])
*来说只能是shape一样的才能相乘,与其他shape的矩阵相乘(*)都会报错,
结果是对应元素相乘作为当前位置的结果,结果矩阵的形状保持不变
CP分解(Canonical Polyadic Decomposition)也称为PARAFAC分解,是一种张量分解方法。它的目的是将一个三维张量分解成若干个二维矩阵的乘积的和的形式,其每个二维矩阵对应于张量的一个模态(mode)。CP分解在信号处理、图像处理、机器学习等领域都有广泛的应用。 下面是使用Numpy实现CP分解的代码: ```python import numpy as np def cp_decomposition(X, rank, max_iter=1000, tol=1e-6): """ CP分解算法 :param X: 三维张量,shape=(I, J, K) :param rank: 分解的秩,即每个矩阵的列数 :param max_iter: 最大迭代次数 :param tol: 收敛阈值 :return: 分解得到的三个矩阵A, B, C """ I, J, K = X.shape A = np.random.rand(I, rank) B = np.random.rand(J, rank) C = np.random.rand(K, rank) for i in range(max_iter): # 更新A矩阵 tmp = np.zeros((rank, rank)) for k in range(K): tmp += np.dot(C[k].reshape(-1, 1), C[k].reshape(1, -1) * np.sum(X[:, :, k] * np.dot(A, B.T), axis=2)) A = np.dot(np.sum(X, axis=(1, 2)).T, np.linalg.inv(tmp)) # 更新B矩阵 tmp = np.zeros((rank, rank)) for k in range(K): tmp += np.dot(C[k].reshape(-1, 1), C[k].reshape(1, -1) * np.sum(X[:, :, k] * np.dot(A, B.T), axis=0)) B = np.dot(np.sum(X, axis=(0, 2)).T, np.linalg.inv(tmp)) # 更新C矩阵 tmp = np.zeros((rank, rank)) for k in range(K): tmp += np.dot(A.T, A * np.sum(X[:, :, k] * np.dot(B, C[k].reshape(1, -1)), axis=0).reshape(-1, 1) * C[k].reshape(1, -1)) C = np.dot(np.sum(X, axis=(0, 1)).T, np.linalg.inv(tmp)) # 判断收敛条件 err = np.linalg.norm(X - np.sum(np.array([np.outer(A[i], np.outer(B[j], C[k])) for i in range(I) for j in range(J) for k in range(K)]), axis=0)) if err < tol: break return A, B, C ``` 使用方式如下: ```python # 生成一个随机的三维张量 I, J, K = 5, 6, 7 X = np.random.rand(I, J, K) # 进行CP分解 rank = 3 A, B, C = cp_decomposition(X, rank) # 检验分解结果 X_hat = np.sum(np.array([np.outer(A[i], np.outer(B[j], C[k])) for i in range(I) for j in range(J) for k in range(K)]), axis=0) print(np.linalg.norm(X - X_hat)) ``` 输出结果为一个较小的数,说明分解结果比较准确。
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值