24小时热门版块排行榜    

查看: 1393  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ywglacier

新虫 (初入文坛)


[交流] 如何改三重for循环-新手

各位,本人新手,如何能把下列三重循环改一下。matlab象这种根本运行不了。
简单来讲,就是三维矩阵每个点都利用fzero求一个等于零的值,然后计算Lup_out。但循环太多了。谢谢。如下:

for i=1:2444
    for j=1:2000
        for k=1:2000
            T_surface =fzero(@(Ts)Sdown_n(i,j,k).*(1-0.5)+Lup(Ts)+Ldown(T_n(i,j,k), RH_n(i,j,k))+G(Ts),0);  % Ts是变量,Sdown_n, T_n, RH_n 都是输入的2000*2000*2444的值;Lup, Ldown, G是函数。
     
           Lup_out(i,j,k) = Lup(T_surface);         
        end
    end
end
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

baby_wolf

新虫 (初入文坛)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
xiegangmai(金币+1): 谢谢分享! 2011-12-19 21:29:47
不夸张的说,你这个问题按这种算法,用一台普通计算机1万年都算不出来。先简化问题、优化算法吧
3楼2011-12-19 18:54:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

hitzhang

木虫 (正式写手)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
ywglacier(金币+1): 2011-12-12 18:12:36
臭水沟(金币+1): 谢谢交流~~ 2011-12-12 18:51:31
矩阵太大了,双精度的话可能要8Gb的RAM才能装得下,所以问题可能不在循环上,试试拆分成小矩阵吧,另外建议你预先赋值会快些。
2楼2011-12-12 15:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见