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

回复此楼

» 猜你喜欢

上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告已经有11人回复
上海工程技术大学激光智能制造课题组招收博士研究生已经有10人回复
实名举报：华南理工大学物理与光电学院副院长李志远婚内两次出轨女学生（博士）已经有15人回复
申博发邮件已经有11人回复
现代”学阀”该如何界定已经有7人回复
各位大神，目前国内有哪些比较好用的逆合成软件？已经有12人回复
我的奶奶已经有3人回复
国社科系统bug了，是不是要放榜了？已经有9人回复

天天快乐

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

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

myw052

金虫 (正式写手)

博士

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

引用回帖:

12楼: Originally posted by lurencyj at 2012-04-24 19:27:53:
sudo说的那个是点操作，意思是，所有操作都是作用于对应的矩阵元素。这种操作适合矩阵大小一样的表达式中。

包括.* ./ .^

嗯，对
我还有一个问题，下面的网址是我的问题，您能帮忙看下么？
http://muchong.com/bbs/viewthread.php?tid=4429347

这个问题也是牵涉高维数组的。比较难解决。

赞一下

回复此楼

天天快乐

13楼2012-04-24 20:03:25

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