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

勇往直前11

代码如下k=1;
$ c" y% [& H& r4 rfor i=1:788    for j=i+1:788
9 a6 b" t1 y5 ?( [. G        P(k,1)=i;
% D& q; h4 U* c0 V% j0 `: \4 z        P(k,2)=j;
. O% o5 i5 ?& @# U) o/ g        k=k+1;
    end
end

/ _1 G& c( n9 K% c, K1 k1 q. I很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行

llbnmo

帮你顶一下

IBB-EUT

在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
- _2 a( w9 I# o7 n. VN = 788;
2 {4 k) O0 X( H  v+ W2 ?% FP = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
% X! M: O/ H5 D; ]M1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
M2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
P = [ M1( M1~=0 ), M2( M2~=0 ) ];