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

回复此楼

» 猜你喜欢

药化及相关博士的申请已经有3人回复
一篇MDPI论文改变了学习工作和生活已经有4人回复
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇已经有4人回复
太白金星有点烦已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有3人回复
河北省自然基金已经有8人回复
收到国自然专家邀请后几年才会有本子送过来评已经有3人回复
考博已经有3人回复
有没有快的中文核心比较快录用的，纳米材料光催化已经有4人回复
有人投过CCC中国控制会议吗？已经有3人回复

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

天天快乐

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

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

lijie169

铜虫 (著名写手)

应助: 159 (高中生)
金币: 161.5
散金: 2018
红花: 7
帖子: 1338
在线: 491.5小时
虫号: 1180408
注册: 2010-12-30
专业: 分离过程

【答案】应助回帖

★
感谢参与，应助指数 +1
myw052: 金币+1, 嗯，这个道理我还是懂的。只是对于高维数组有时候想要实现比较麻烦。 2012-04-24 18:56:50

matlab上说矩阵运算比for循环高效很多

赞一下

回复此楼

4楼2012-04-24 12:24:47

已阅回复此楼关注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