多重for循环速度奇慢，求指点~

小弟程序中使用了一段多重for循环，运行10小时进度才不到三分之一…
网上都说应该尽量向量化，但水平实在有限，不知该如何优化，恳请不吝赐教！论坛币不多，聊表心意！
该段程序如下：

h=waitbar(0,' ');                 
f_bar_star=zeros(2210,200);
sigma=zeros(2210,200);
bty=zeros(3533,200);
for n=1:200                        
    for k=1:2210
        f=zeros(1,3533);
        for i=2:3533
            bty(i,n)=(btp(i+150,n)-btp(i+149,n))/btp(i+149,n);
            S(k,i,n)=single(btsignal(k,i,n));                                       
                                               %btsignal是个三维数组，之前为了能存储，就设为int8格式，这里转回single以便计算自然对数
            f(k,i,n)=log(1+bty(i,n)*S(k,i-1,n))-log(1+bty(i,n)*S0);
            sigma(k,n)=sigma(k,n)+f(k,i,n);
        end
        clear f;                                                 %只能每次小循环都预设和清除一次f，直接预设三维数组就会报out of memory
        f_bar_star(k,n)=1/3533*sigma(k,n);
    end
    waitbar(n/200);
end

不知道你的btp是个变量还是个方程。

如果是个变量，最里面那个循环的第一步可以写成
(btp[k:end,n]-btp[k-1:end-1,N])./btp[k-1:end-1,N]

如果是个方程（否则好像size不对吧...），well，那就没招了。如果你是多核的，用 parfor吧，能把运算时间缩到大约1/#worker

三层嵌套肯定慢了。200*2210*3532次循环。

12核24线程

sddyl 发表于 2013-4-1 04:58
不知道你的btp是个变量还是个方程。

如果是个变量，最里面那个循环的第一步可以写成

sddyl兄：
btp是一个变量，二维数组哈。
(btp[k:end,n]-btp[k-1:end-1,N])./btp[k-1:end-1,N]  →→→→为啥是用方括号呢？

Davidcharge 发表于 2013-4-1 20:46
sddyl兄：
btp是一个变量，二维数组哈。
(btp[k:end,n]-btp[k-1:end-1,N])./btp[k-1:end-1,N]  →→→→ ...

python写串了...应该是括弧 ： ）

sddyl 发表于 2013-4-6 11:22
python写串了...应该是括弧 ： ）

虽然我后来自己搞定了，还是谢谢！
向量化之后原来跑30个小时的程序一分钟就能运行得到结果

求助楼主啊，我也是for循环很慢，都等了30多个小时了，能不能楼主指导一下如何向量化啊，急用，谢谢

能用矩阵就别用for loop

转角遇到爱981 发表于 2013-10-24 10:29
求助楼主啊，我也是for循环很慢，都等了30多个小时了，能不能楼主指导一下如何向量化啊，急用，谢谢

可以加我QQ 182456270