|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
最近需要用Matlab计算数据核矩阵,数据X_joint大小是40000x183,需要40000x40000的核矩阵,编写了下面的代码,但是for循环感觉速度非常慢,有什么办法能改进这个代码吗?2 b* X7 v) J; g1 G. n w$ {
K = zeros(40000,40000);
2 B# t: m9 k+ s9 ~7 d. L: ?for i = 1 : 40000' x1 q) P/ B: H) [; I
for j = 1 : 40000
0 V. q9 ^6 J S K(i,j) = exp(-(norm(X_joint(i, - X_joint(j, ))^2/c);%求核矩阵,采用径向基核函数,参数c=2000
8 X6 p: t6 H- F/ U: h- wend
6 o0 v0 G! J4 [$ l& \ I3 ?( iend d& Y4 ^4 E# y: q, n1 [9 Z) K
; N, L% ]5 A* L5 u( a( }+ C+ f9 g& M: `0 L
/ ~8 S; {( L1 \, Y" U
r n, {+ `. r& r, m3 a' O! K2 N, P% \7 c% C5 ]) {" C
|
|