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

勇往直前11

代码如下k=1;
for i=1:788    for j=i+1:788
        P(k,1)=i;
; S/ T- X- |% y! n        P(k,2)=j;
        k=k+1;
    end
end

* D# ^/ G& X% K  _% p9 h; y: M很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行
6 `$ b0 \0 E7 a5 o7 G" I

llbnmo

帮你顶一下

IBB-EUT

在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
5 m( n9 y. D, }( [N = 788;
; Q, |4 n/ t; R1 k6 N& }3 ?: pP = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
M1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
9 a  V5 A" \( ^M2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
! _3 a' {2 O, aP = [ M1( M1~=0 ), M2( M2~=0 ) ];
1 }% ~/ S9 a( e2 L% _: a7 y9 r/ v) {