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

勇往直前11

代码如下k=1;
for i=1:788    for j=i+1:788
+ @7 v$ U& Y+ |$ T        P(k,1)=i;
        P(k,2)=j;
        k=k+1;
    end
1 U8 J( ~2 D( s2 Y  {9 n" Q0 B0 fend
4 o# ]) _" ]* k# `% ?% [% v
" D+ Y1 {( w- X. M. Y/ s- J很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行
6 G( [, Z0 B4 ^* `

llbnmo

帮你顶一下

IBB-EUT

在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
1 Q/ d' B  U) o, EN = 788;
8 i2 P& v9 h+ D* y( y. u* H: Q2 _P = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
* [# l& Z% W0 C2 A) n& s) IM1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
M2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
P = [ M1( M1~=0 ), M2( M2~=0 ) ];