|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
最近需要用Matlab计算数据核矩阵,数据X_joint大小是40000x183,需要40000x40000的核矩阵,编写了下面的代码,但是for循环感觉速度非常慢,有什么办法能改进这个代码吗?
! E! }& D# P) V: i( \: D |* xK = zeros(40000,40000);- s# o) I. @$ _8 b0 V9 f
for i = 1 : 40000
- X; G1 S1 }1 V: t6 X: U+ Bfor j = 1 : 40000
]0 Q& |4 t$ o1 `! u6 q K(i,j) = exp(-(norm(X_joint(i, - X_joint(j, ))^2/c);%求核矩阵,采用径向基核函数,参数c=2000
8 I0 G" C3 p4 T3 s8 Send0 K( E$ l, `; ^, ]% x) ^9 l: V
end. M4 b# {1 `) C( c! G
& K1 I. _& F" Y( Y: p7 z. [. g, t. x" m" ~$ {4 O
. o7 Y; l: U# S' _1 `' s/ N5 k4 I! v: d! o$ G8 K) ^ d
! B8 v. O6 @# ^; Y8 [# z! \ |
|