|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
最近需要用Matlab计算数据核矩阵,数据X_joint大小是40000x183,需要40000x40000的核矩阵,编写了下面的代码,但是for循环感觉速度非常慢,有什么办法能改进这个代码吗?) ~9 T; h0 m$ j' D# F6 @8 ?
K = zeros(40000,40000);
/ F5 s- Y( |2 g5 c6 J. R6 L* mfor i = 1 : 400001 Z8 l7 Y5 F' f! D2 m* [1 ~
for j = 1 : 40000
3 X" ^7 s: p$ P. T9 D K(i,j) = exp(-(norm(X_joint(i, - X_joint(j, ))^2/c);%求核矩阵,采用径向基核函数,参数c=2000; Y6 J2 U/ B C0 ?
end
; R# p& y5 }6 {/ Oend: ?% t' L; z) _# ^
3 T: @' _7 x0 ^/ [: C
! v7 G, q+ H a
$ q" R1 V- g$ V! F4 q
7 v" C! O8 j0 m. d
8 u) \7 h! I% ?! t- H |
|