【矩阵分析】向量求和转为矩阵

B = [ B 1 , B 2 , ... , B N ] T ∈ R N × 1 \boldsymbol B = [B_1,B_2,\dots,B_N]^T \in \mathbb{R}^{N \times 1} B=[B1,B2,...,BN]T∈RN×1, W = [ W 1 , W 2 , ... , W N ] T ∈ R N × N r \boldsymbol W = [\boldsymbol W_1,\boldsymbol W_2,\dots,\boldsymbol W_N]^T \in \mathbb{R}^{N \times N_r} W=[W1,W2,...,WN]T∈RN×Nr, H = [ H 1 , H 2 , ... , H N ] T ∈ R N × N r \boldsymbol H = [\boldsymbol H_1,\boldsymbol H_2,\dots,\boldsymbol H_N]^T \in \mathbb{R}^{N \times N_r} H=[H1,H2,...,HN]T∈RN×Nr

将 ∑ n ∈ N ( W n H n B n ) \sum_{n \in \mathcal{N}}(\boldsymbol W_n \boldsymbol H_nB_n) ∑n∈N(WnHnBn)写成矩阵形式
∑ n ∈ N ( W n H n B n ) = Tr ( W ⊤ diag ( B ) H ) \sum_{n \in \mathcal{N}} (\boldsymbol{W}_n \boldsymbol{H}_n B_n) =\text{Tr}(\boldsymbol{W}^\top \text{diag}(\boldsymbol{B}) \boldsymbol{H}) ∑n∈N(WnHnBn)=Tr(W⊤diag(B)H)

代码验证

matlab 复制代码
% 定义参数
N = 5; % 样本数
Nr = 3; % 每个向量的列维度

% 随机生成矩阵和向量
W = rand(N, Nr);  % N x Nr 矩阵
H = rand(N, Nr);  % N x Nr 矩阵
B = rand(N, 1);   % N x 1 列向量

% 逐项求和方式计算
sum_result = 0;
for n = 1:N
    sum_result = sum_result + W(n, :) * H(n, :)' * B(n);
end

% 矩阵形式计算
diag_B = diag(B); % 对角矩阵
trace_result = trace(W' * diag_B * H);

% 显示结果
disp('逐项求和结果:');
disp(sum_result);
disp('矩阵形式结果:');
disp(trace_result);

% 验证是否相等
if abs(sum_result - trace_result) < 1e-10
    disp('验证成功:两者相等!');
else
    disp('验证失败:两者不相等!');
end
相关推荐
】余1853816280011 小时前
矩阵系统源码搭建与定制化开发,支持OEM
线性代数·矩阵
YuTaoShao2 天前
【LeetCode 热题 100】73. 矩阵置零——(解法一)空间复杂度 O(M + N)
算法·leetcode·矩阵
Dark__Monarch2 天前
二元一次方程
线性代数
Kaltistss2 天前
240.搜索二维矩阵Ⅱ
线性代数·算法·矩阵
说私域2 天前
视频号账号矩阵运营中定制开发开源 AI 智能名片 S2B2C 商城小程序的赋能研究
人工智能·矩阵·开源
张晓~183399481213 天前
数字人源码部署流程分享--- PC+小程序融合方案
javascript·小程序·矩阵·aigc·文心一言·html5
峙峙峙4 天前
线性代数--AI数学基础复习
人工智能·线性代数
我爱C编程4 天前
基于拓扑结构检测的LDPC稀疏校验矩阵高阶环检测算法matlab仿真
算法·matlab·矩阵·ldpc·环检测
CVer儿4 天前
svd分解求旋转平移矩阵
线性代数·算法·矩阵
张晓~183399481214 天前
数字人分身+矩阵系统聚合+碰一碰发视频: 源码搭建-支持OEM
线性代数·矩阵·音视频