MATLAB中怎样加快for循环矩阵赋值运行速度

tryeverythingop

最近需要用Matlab计算数据核矩阵，数据X_joint大小是40000x183，需要40000x40000的核矩阵，编写了下面的代码，但是for循环感觉速度非常慢，有什么办法能改进这个代码吗？
, H& @& M3 I% V. y# t& }K = zeros(40000,40000);
+ `& ^3 Z; b9 L- K7 [8 {6 Z; Bfor i = 1 : 40000
% [9 v2 L0 V9 I9 x/ y. `for j = 1 : 40000
- O' S7 u7 C6 {7 K6 C     K(i,j) = exp(-(norm(X_joint(i, - X_joint(j,))^2/c);%求核矩阵，采用径向基核函数，参数c=2000
0 B+ c+ C# ^* D' D. @; L5 x) Send
; \" q& k9 e, Jend


6 \, i# B. \, y0 q/ \

. ]) P+ C2 q, q" \: @. U1 {
) A. A: V! s) T$ x

qq666888qqw

形如 norm(X_joint(i, - X_joint(j,) 的项，实际上对于 i, j 是对称的，这上面可以省掉一半运算步数，对角线是零，你只需要计算上三角或者下三角就够了。
1 N3 ~0 a8 ~" G5 {/ q4 d但本身两层循环的步数实在太多，而且每次都要计算两个183维数据的差再求模，除非你能预先通过其他评估，把183个维度合并或者降低到远小于此的维度，不太容易改进了。