24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

myw052

金虫 (正式写手)

博士

应助: 9 (幼儿园)
金币: 542
散金: 300
红花: 2
帖子: 409
在线: 164.1小时
虫号: 1320903
注册: 2011-06-11
性别: GG
专业: 固体力学

[求助] matlab for循环速度提高方法

下面的程序是我用matlab写的一个程序，程序中需要用到for循环，由于是高位矩阵，在运行的时候循环速度效率很低，造成计算成本很大。在许多论坛上寻找好的方法来解决matlab运行中循环的速度提高方法，始终没有找到。有人说用矢量化解决问题，但是描述的不具体，操作起来有困难。有人说用“少拿多取”的思想，但是在简单情况下，可以理解，像如下的情况，似乎难以实行。

故在本版发帖，希望斑竹或者各位程序达人能出来指点一二。本人不胜感激。
（注：此贴绝不是无病呻吟）。

function for_cycle()
tic
clc
clear
close all
beta=rand(10^5,10,10);
gama=rand(10^5,10,10);

[Row,Column,Volume]=size(beta);
gama_u=rand(1,1);
m=2;

for k=1:Volume
k
for i=1:Row
      for j=1:Column
         D(i,j,k)=beta(i,j,k)/(beta(i,j,k)+(1-beta(i,j,k))*(gama(i,j,k)-(gama(i,j,k)/gama_u)^m)/(gama(i,j,k)-1));
      end
end
end

D
toc
end

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有8人回复
三无产品还有机会吗已经有6人回复
投稿返修后收到这样的回复，还有希望吗已经有7人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

天天快乐

1楼2012-04-24 09:38:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhfzh

木虫 (正式写手)

应助: 19 (小学生)
金币: 6620.5
散金: 200
红花: 4
帖子: 685
在线: 665小时
虫号: 715595
注册: 2009-03-05
专业: 原子和分子物理

你拭下用矢量化计算
就是把for,end去了，常数乘个同维度的单位阵，所有的乘除号前加点
还有你算beta,gamma的时候不要重复调用，用个临时变量存起来
还慢的话你用C语言写吧

赞一下

回复此楼

5楼2012-04-24 12:31:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 16 个回答

sudo

木虫 (正式写手)

程序强帖: 16
应助: 6 (幼儿园)
金币: 1297.6
散金: 1486
红花: 20
帖子: 588
在线: 641小时
虫号: 1211394
注册: 2011-02-24
性别: GG
专业: 文艺美学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xzhdty: 金币+1, 欢迎常来 2012-04-24 13:54:50
myw052: 金币+7, ★★★★★最佳答案, 您的程序写的很好，矢量化确实很厉害。 2012-04-24 18:52:29

把规模改小了，方便对比测试

CODE:

tic

clc

clear

close all

beta=rand(10,10,10);

gama=rand(10,10,10);

[Row,Column,Volume]=size(beta);

gama_u=rand(1,1);

m=2;

D=zeros(size(beta));

for k=1:Volume                    

    for i=1:Row                    

        for j=1:Column                

            D(i,j,k)=beta(i,j,k)/(beta(i,j,k)+(1-beta(i,j,k))*(gama(i,j,k)-(gama(i,j,k)/gama_u)^m)/(gama(i,j,k)-1));

        end

    end

end

D

toc

tic

D_vectorized = beta ./ ( beta + (ones(size(beta)) - beta) .* (gama-(gama/gama_u).^m) ./ (gama - ones(size(gama))) )

toc