找回密码
 注册
关于网站域名变更的通知
查看: 597|回复: 3
打印 上一主题 下一主题

双层嵌套循环赋值怎么优化?

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-1-15 13:34 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
代码如下k=1;
: G# j% X! t2 o0 j1 v0 z1 xfor i=1:788    for j=i+1:788+ B2 W2 S6 U+ r# ^  r) S% v4 ~" T! Z* n
        P(k,1)=i;
. E) ?% n9 O. z& t        P(k,2)=j;
1 k& o( G$ \4 s1 E        k=k+1;
. L$ b* L6 z1 G    end# C. z- X1 ~3 u+ f6 |9 e% e
end
+ J9 u1 ]% t$ u. w& v; u! y5 _* G
7 \5 E4 W1 j* E很简单的代码,但是运行非常耗时要运行几天看起来.求向量化怎么改,或者简化就行
# \' }& t2 O- w  @) d* G3 o

该用户从未签到

2#
发表于 2021-1-15 13:46 | 只看该作者
帮你顶一下

该用户从未签到

3#
发表于 2021-1-15 16:03 | 只看该作者
( h6 _0 D0 I+ R' B1 r
在循环开始前,预设好矩阵P的维度,就可以大幅节省时间了
1 \( P$ U7 D5 e4 q* XN = 788;
0 B( |. r7 o- d4 |P = zeros( N * ( N - 1 ) / 2, 2 );

该用户从未签到

4#
发表于 2021-1-15 16:58 | 只看该作者
N = 788;
3 T8 c$ E% x( t. r) yM1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
) ^' e7 ?% T  Z" l# t2 ^M2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );: s2 G' w% [& R2 u- |5 E
P = [ M1( M1~=0 ), M2( M2~=0 ) ];$ h' ^" h& f9 _/ `7 M" Y% t
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-23 20:46 , Processed in 0.156250 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表