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

勇往直前11

代码如下k=1;
3 ^/ J; C, C  K- ofor i=1:788    for j=i+1:788
        P(k,1)=i;
        P(k,2)=j;
        k=k+1;
1 `5 j. G0 v5 L# [! C    end
6 M9 @; }; w2 V* B$ d6 F6 jend
$ \5 T; g$ B7 ^* S
$ V. T) u: q0 i很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行

llbnmo

帮你顶一下

IBB-EUT

在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
4 V: @5 {. O: ~& V2 ]N = 788;
' @- P0 d. m9 ]& C' pP = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
3 e5 g' }# J$ ?6 q7 S% R' aM1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
3 W* d% R4 Y3 Y- x1 pM2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
) t' K* [; `% p( u0 gP = [ M1( M1~=0 ), M2( M2~=0 ) ];
: U/ H5 S) r. I# J/ g