|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
, r" ]: D, d* Q0 b) Nmatlab在计算大数据内存以及大矩阵运算时,单核运算显然无法满足高速的运算需求。其实matlab提供多核运算的解决方案,这里先介绍最简单的两种5 w O! q. h* ]8 T/ ?( Q% j5 b
A# V8 K. u4 D①paRFor' B6 h$ Y' y; _, [9 T1 P" S
②distributed
! M5 b: ~) Y. B$ s1 ?$ ^ _# F/ ~+ n1 o9 S6 c: @& ~; O
parfor8 W- }, ?3 v% R! ~
for循环一般最为耗时,解决方法就是采用parfor代替for,parfor可以自动将循环任务分配到多个核中,现在matlab的命令行窗口运行parpool(4)查看是否可以并行计算,没有报错则并行已打开。parfor运行时就会有多个matlab进程在计算。计算效率明显提升。3 X" l% G* ^0 |6 ]* ]4 @
, N& T& z" w: |4 H8 r
未使用parfor的代码:6 B8 \# |( Q$ B% x6 G/ L! V' A2 Q& B
8 u4 q5 u5 J# |- answer = 0;
- for i=1:1000000
- answer = answer + i;
- end. C" L8 Z! i# W$ T# e/ r
% l3 G# q- }9 T% U
4 a' R" a) }7 w6 N
时间:. I" Y' ^: b) c: |: [
9 b T8 _0 W/ p+ X
" x O, Q4 r3 V' s
8 p3 |3 P( g) ?* Q3 l使用parfor的代码:
: W5 X$ g3 X6 p) O. V6 j& s1 p" M
7 O7 p5 |+ T; N, K: ]& h) D: Z- answer = 0;
- parfor i=1:1000000
- answer = answer + i;
- end5 C5 Y( R2 P9 ^6 D2 b
7 o. D4 U8 h8 u$ T
0 w7 K X4 \$ ~( m
时间:* C% h% Y, |2 l9 t( d/ e
! H! _# ]2 Z( s" Q8 n0 S
! m! q: s7 S# u! Y5 Y! o+ r& V( i% J g; f. e
显然加了parfor语句速度更快!
}3 L% _. M" U! f
7 Z' N' S3 W2 J5 O注意:
$ `. D! ]% ]( l0 t' u4 {! x4 r- l7 L* t* B) B
①parfor循环中不能使用迭代或者关联性的赋值语句,因为多个核计算时无法交换数据。
& I! j' k8 e' c" l8 p②parfor只支持一层循环,所以下面这种写法是非并行的,并且可能会出错。
+ w$ ]$ D" K+ N/ K5 F' H' p& t③在程序运行之前MATLAB会提示你哪些地方不能用parfor(红色波浪线显示)
; \( ?& I+ q+ K& C1 S
2 H7 x) n; A6 ?distributed9 w8 N7 b: ?3 b' T% W
对于大矩阵的存储问题以及计算的问题,除了自己分块外,可以使用matlab自带的内存分配解决方法,但需要确保你的电脑内存足够多,计算效率就会大幅提高。不然会有内存溢出的问题。 L8 Q3 H5 [6 N7 b. R6 o) O6 V E
具体方式是
- ^7 z" e4 j! A: G$ r+ J
) n/ ]3 A9 a6 j* `- tic;
- X = randn(10000);
- x_dis = distributed(X);
- clear(X);
- whos
- a = svd(x_dis);
- a_out = gather(a);%distrituted数据类型变为double
- toc;- i) {4 W. K2 G# s' }
1 q/ {1 o8 n& [) q& e
+ K5 G7 ]6 _* S5 n3 E之后可以采用x_dis代替X进行计算即可。gather函数的作用是把distributed数据类型变为double。
# ]! N" S2 l' x) {5 s# ~2 q( m7 L$ v$ R
9 g' d1 I- X4 ~0 |; P" U8 G9 e8 x此外,matlab的大多数工具箱都开始支持并行计算,所以如果该问题可以用自带工具箱解决的可以使用工具箱自带的并行计算方法。6 `# C' c( Y- f0 R) h; _
matlab的Parallel Computing Toolbox工具箱就支持并行。 |
|