急求帮改程序将矩阵原本的常数改为3*3矩阵元素(P只看第一個,P有誤已更正)

真龍虎王

3 A0 R/ N7 O0 Q7 Q; y
% e5 x; R, t4 m& ~原本题目:                                     新增
   NO              GO          TO            P
    1                 0             1             7
    2                 1             2             7
    3                 1             4             1
) L) Z+ A5 S; N- }' v' e    4                 2             3             6
原程序:
% x* Y- I( S) o: UNO=[1;2;3;4];
GO=[0;1;1;2];
+ d1 ~, [+ J; J6 iTO=[1;2;4;3];
k=zeros(length(NO),length(NO));
a=sub2ind(size(k),NO,TO);
# q) n4 R& J$ K, j0 tk(a)=-1;
& e/ c2 v2 Q/ b! Ifor i=1:length(NO)
    s=i;
# K& [0 a" }0 z' @3 c% s4 x    while s>=2
5 H. y( g$ ~; `        s=find(TO==GO(s));
        k(NO(s),(TO(i)))=-1;
3 l! H4 `4 @( s1 C    end
" `7 z  v8 X7 r* c  d$ g0 Zend
( {/ D" H" y. Z7 q) w$ L答:
: s/ ^5 }. J6 q, E3 iK=                TO    =             1    2     3    4           =[-1 -1 -1 -1;0 -1 -1 0;0 0 0 -1;0 0 -1 0]
0 f! S$ h, Q6 \            NO  [     ]         1   [  -1   -1    -1   -1
                                   2       0   -1    -1    0
4 a- \- v9 m. a2 j# y1 k- d                                   3       0    0     0   -1
- _; L1 D6 e- q# L0 k                                   4       0    0    -1    0 ]
现在想改多看一个矩阵P    矩阵大小跟NO相同
' @$ Q' @& d% F4 w0 l+ S$ Q% t如 P=[7;7;1;6]
- v' m  b8 B; Q5 i题目想改为
$ b$ f$ E* O  n2 u  [/ A看NO找TO值再看GO值时顺便看P值
原本K是放入-1想改为当:

$ h% P1 W; L* b. x  l. L0 WPS:只看NO找TO时找到的第一个P值如果GO不是0那下一个的P不使用
只用第一次找到的P

P=7    -1改为放入[-1 0 0;0 -1 0;0 0 -1]
  =6    -1改为放入[-1 0 0;0 0 0;0 0 -1]
( ]! A- `4 ]( E  Z- c' E  =5    -1改为放入[0 0 0;0 -1 0;0 0 -1]
  =4    -1改为放入[-1 0 0;0 -1 0;0 0 0]
+ l2 e- C1 x1 i1 P" u/ k1 D1 r  =3    -1改为放入[0 0 0;0 0 0;0 0 -1]
  =2    -1改为放入[0 0 0;0 -1 0;0 0 0]
. \1 ]2 e8 U1 y0 W3 |, ~  =1    -1改为放入[-1 0 0;0 0 0;0 0 0]
( `  s" g3 Y. R2 I& @4 F
那原本的矩阵的零元素改為零3*3矩阵=[0 0 0;0 0 0;0 0 0]
目前卡关是卡在k矩阵大小是会因为P而变以及-1要改放入矩阵
答案会是
-1   0   0    -1   0   0     -1   0   0     -1   0    0
7 B5 q- h, b. z( [. p% w, z% e4 Q8 c0  -1   0     0  -1   0      0   0   0      0   0    0
0   0  -1     0   0  -1      0   0  -1      0   0    0

& c2 N. y+ L7 z0   0   0    -1   0   0     -1   0   0      0   0    0  
0   0   0     0  -1   0      0   0   0      0   0    0
0   0   0     0   0  -1      0   0  -1      0   0    0
) I/ O, o4 N6 q* D& }5 h3 @
; S0 @8 s; W( N# ?, n0   0   0     0   0   0      0   0   0      -1   0    0  
; ?. I, N, c7 f: y0 ~: n. v' v6 f0   0   0     0   0   0      0   0   0      0   0    0
& X- y. p/ Q" `  f6 z+ _: E0   0   0     0   0   0      0   0   0      0   0    0

0   0   0     0   0   0      -1   0   0      0   0    0  
0   0   0     0   0   0      0   0   0      0   0    0
0   0   0     0   0   0      0   0  -1      0   0    0