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

勇往直前11

代码如下k=1;
- Y3 Y( d( {# s, f. pfor i=1:788    for j=i+1:788
6 X! U5 h; A9 b        P(k,1)=i;
/ x" R' g! W/ p' Z* C9 ~) ?        P(k,2)=j;
5 U' V8 s( x, m! L; g        k=k+1;
    end
* H2 e+ t' r8 \) m5 Tend

很简单的代码，但是运行非常耗时要运行几天看起来.求向量化怎么改，或者简化就行
% H& B  I- f- S4 A* t2 {/ y

llbnmo

帮你顶一下

IBB-EUT

/ P, O2 G' d8 m+ A* k2 s在循环开始前，预设好矩阵P的维度，就可以大幅节省时间了
% f8 v0 m' J- N, W6 mN = 788;
& V% ~; E( `5 D! d( fP = zeros( N * ( N - 1 ) / 2, 2 );

zzz.dan

N = 788;
M1 = tril( ones( N - 1, 1 ) * [ 1 : N - 1 ] );
1 X; F: u2 H8 J% x$ r& {M2 = tril( [ 2 : N ].' * ones( 1, N - 1 ) );
P = [ M1( M1~=0 ), M2( M2~=0 ) ];