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

勇往直前11

代码如下k=1;
/ S8 q8 ^, z, y- Y7 @! ]& @+ ^6 K& wfor i=1:788    for j=i+1:788
8 q+ o9 U$ V& ^, a  J; W) W        P(k,1)=i;
        P(k,2)=j;
0 q; r! z" t9 D/ M4 K        k=k+1;
    end
  |) t: s  I+ i- cend

很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行
7 _: J5 y1 @3 x/ H/ M

llbnmo

帮你顶一下

IBB-EUT

在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
) k2 O/ l. V& S5 j: QN = 788;
* e. n( V% t1 E: y6 L5 k) ?P = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
M1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
; z2 `' Y/ d) V" ZM2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
P = [ M1( M1~=0 ), M2( M2~=0 ) ];
! p8 V' D  z6 O% Z, }