24小时热门版块排行榜    

查看: 1724  |  回复: 5

hjnwsuaf

金虫 (著名写手)

[求助] 如何提高多层循环速度?见程序 已有1人参与

程序如图片显示内容。
当参数lp较小时,运算可以可以进行,当很大超过500后,就会导致死机。
如何解决这个问题,提高运算速度呢?

具体程序如下:
clc
clear all
lp=50;
a=-lp:1:lp;
xh=length(a);
tic
for i=1:1:xh
    for j=1:1:xh
        for k=1:1:xh
            rt(i,j,k)=a(i)^3+a(j)^3+a(k)^3;
        end
    end
end
toc

如何提高多层循环速度?见程序
1.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duzhaolong

金虫 (小有名气)

123万物

用矩阵乘法看看,或者换用优化过的算法,3层for循环,数据量大的时候确有点慢

发自小木虫Android客户端
努力感悟,但有所得足矣
2楼2015-12-10 09:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
3楼2015-12-10 09:33:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxws.2008

木虫 (正式写手)

matlab循环慢,用C语言或fortran

发自小木虫Android客户端
4楼2015-12-10 09:46:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxws.2008

木虫 (正式写手)

变成3个一维数组,然后合并,根本不用循环

发自小木虫Android客户端
5楼2015-12-10 09:52:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
hjnwsuaf: 金币+20, ★★★很有帮助 2015-12-10 17:14:02
把三次方运算提到外面,并用矩阵操作试试。我的电脑慢,卡住了。
clc
clear all
lp=50;
a=-lp:1:lp;
a = a.^3;
a = transpose(a);
xh=length(a);
tic
A = repmat(a,[1,xh,xh]);
B = permute(A,[2,1,3]);
C = permute(A,[3,2,1]);
R = A + B + C;
toc
知其然,知其所以然。
6楼2015-12-10 13:39:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hjnwsuaf 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 03:50 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:24 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-02 03:18 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-02 03:11 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:45 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 21:26 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
信息提示
请填处理意见