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

勇往直前11

代码如下k=1;
- z6 `+ j! a3 ]1 |7 k8 f9 Vfor i=1:788    for j=i+1:788
7 f( @4 u" t, ]% R        P(k,1)=i;
$ q1 C' e* B% [* j; n        P(k,2)=j;
5 N1 M6 k3 W2 J+ i( L! e0 s        k=k+1;
) P" \2 x1 I: L2 T- p( e4 x    end
end
( Z1 V; q! c! w3 o0 |
: v) D( [0 x/ c1 g1 f( u很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行
8 ~: z# f" }4 l

llbnmo

帮你顶一下

IBB-EUT

2 c7 J7 d8 t$ @/ J/ I在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
' b1 J$ \: {* z# mN = 788;
P = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
) t: Z6 d0 P  H2 O2 vM1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
M2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
P = [ M1( M1~=0 ), M2( M2~=0 ) ];
2 y; r( K5 [) ]5 w/ A1 T$ G