每个节点的生成树数目程序

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   & u( R: B: q7 v% J- v1 G% O
2. -1 0 -1 1 0 0 0;& b: r# \2 q\" y2 Z
   
3. 0 -1 1 0 1 1 0;
   2 j9 s2 D; q+ l+ _1 D% n& T
4. 0 0 0 -1 -1 0 -1;# B- r1 A# l, A/ [4 O8 O8 E
   
5. 0 0 0 0 0 -1 1 ]2 I' D, s# u1 d6 L7 p% p
   
6. for i=1:5
   9 H  ?8 `# p0 i
7. a=c;
   4 m5 D+ r% Y, O\" f6 [% t6 Q
8. b=a(i,:)
   1 j- r) {8 ~: t! T; L9 {. ^% C
9.      for j=1:46 O\" g\" s5 T. n5 p: Z
   
10.         b(i,j);
    1 c9 w2 O$ q7 A7 C$ T
11.             if b(i,j)~=0
    ' n\" O% O2 e4 [1 h! N
12.                a(:,j)=[]# y- h, h  k2 o! y, V
    
13.             end1 t* B1 ?+ R- U1 g\" Z7 k- d
    
14.      end1 r3 Z1 _# N1 z) B  k$ Q
    
15. a(i,:)=[]  W% D7 k# T% K3 Y: N- v- B' I5 O
    
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
3 p0 n) X/ P8 x/ S) A" {

madio

1. c=[1 1 0 0 0 0 0;
   % u& b  P% m2 g$ I9 N/ z
2. -1 0 -1 1 0 0 0;
   ! q\" J% K% v- z* d5 B* O1 r( {. m
3. 0 -1 1 0 1 1 0;) w* i: ~2 ~/ `8 P
   
4. 0 0 0 -1 -1 0 -1;% `, K+ s' \0 h
   
5. 0 0 0 0 0 -1 1 ]
   . x/ P4 ?; I+ X; j$ m, `4 O
6. for i=1:5\" e& F2 R( i3 E! U9 T  R
   
7. a=c;  `) b' ]; b: S9 W+ L5 l# c( `
   
8. b(i,:)=a(i,:)
   \" \9 H; _$ y( Q9 J! |4 J
9.      for j=1:4
   ' e- j$ U* W\" S\" U
10.         b(i,j);1 K5 _) r+ y- E: M
    
11.             if b(i,j)~=0
    : c4 v' Q: u0 g( Y% ^2 b0 O
12.                a(:,j)=[], F) r. `( U& l
    
13.             end: t3 y3 v) N\" l. G% F/ Q
    
14.      end
    ) a6 O* \9 W% Z6 T
15. a(i,:)=[]: |4 U4 Y( t5 G- ~2 W; j0 t8 x9 Q% X
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
7 A7 o' s6 O# l3 o改成这样的就可以运行了

不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

1. c=[1 1 0 0 0 0 0;: n6 y9 w! }  e$ l
   
2. -1 0 -1 1 0 0 0;
   9 z4 Z5 }' M% U& r+ ~! R
3. 0 -1 1 0 1 1 0;6 j8 B8 N5 C2 d: p2 k# [* E
   
4. 0 0 0 -1 -1 0 -1;
   . n& K% H5 s) W' a! L; V
5. 0 0 0 0 0 -1 1 ]
   1 I  t- A; ^\" Z5 i
6. for i=1:5
   ) N+ P3 {6 V8 n1 \+ R
7. a=c;+ _' U0 x6 a, [
   
8. b=[];1 ?3 V& [1 U' W8 {) G- Q+ g3 T, K
   
9. del_col=[];
   & d& m  @3 N+ E: k& T
10. b(i,:)=a(i,:)7 E. }8 \$ {) @0 u* i2 L8 j
    
11.      for j=1:7
      }  N, ^( c, ~- K- C# X
12.         b(i,j);3 m: R. K) l( z# K! P
    
13.          if b(i,j)~=0
    & ?/ K, O+ ~: O\" v3 M4 r; H- w* Q
14.            del_col=[del_col,j];
    # U& B' m0 P1 [$ \& B2 i
15.          end
    / V1 a: E7 L2 f* z2 c
16.      end
    7 H6 x# O. r5 ~\" G! \\" n
17.      a(:,del_col)=[]\" Z6 u4 L4 z- W( {2 y
    
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
2 m: I3 j3 f3 P/ f' e; o+ k你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

, F- Q- i3 F, E7 l- X可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！