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

xiangpeng395

1. c=[1 1 0 0 0 0 0;  ]# z8 Q% D$ g
   
2. -1 0 -1 1 0 0 0;$ ?6 y) ?4 `, b# V* F9 o! H
   
3. 0 -1 1 0 1 1 0;) {1 ]. w9 d! @; W2 U7 x
   
4. 0 0 0 -1 -1 0 -1;
   \" q- i3 D/ ^1 k) i
5. 0 0 0 0 0 -1 1 ]
   , c5 Z& T2 L6 \) H+ n9 Y
6. for i=1:5
   , ?* w1 i* p6 h4 Z4 y; \' d7 M% i
7. a=c;\" g- h0 l; i$ i- s# j
   
8. b=a(i,:)3 r8 x/ |( S5 ?6 D8 q/ W
   
9.      for j=1:4
   2 h\" J& F( G, X% X. Q
10.         b(i,j);0 \0 \\" K: z) b
    
11.             if b(i,j)~=0
    * Q  l9 O. r8 r/ s6 V- ?1 Z
12.                a(:,j)=[]
    3 n& h0 U5 @* [* F( t
13.             end
    & ~7 g2 m* C. O% E3 x) G
14.      end1 N* p9 c- J% i/ \9 n/ V! C
    
15. a(i,:)=[]9 W, q' }: ]4 f. s; M5 `: O
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;
   ' B9 h! r: Z8 w$ r4 l5 W4 e( m
2. -1 0 -1 1 0 0 0;* o( o. J# W% z; g7 ?+ x7 d  q  ^9 B
   
3. 0 -1 1 0 1 1 0;4 h7 J; T2 C. ]& S, L
   
4. 0 0 0 -1 -1 0 -1;5 s, }+ R\" O: w7 @9 {/ r5 S
   
5. 0 0 0 0 0 -1 1 ]' [  b3 Q0 R: x! c
   
6. for i=1:5
   2 U% \& p  b6 b7 ?* w$ K  }
7. a=c;
   ( Q1 a( p& F4 [* [- l
8. b(i,:)=a(i,:)
   8 A: @\" U. x! |% H
9.      for j=1:4
   \" n- I7 P; Z* f4 U) L0 O
10.         b(i,j);
    7 K+ A, h1 q, i6 k
11.             if b(i,j)~=0  f; g4 g; b/ H% v. Y) I6 n( J
    
12.                a(:,j)=[]( B2 x1 U3 Z7 t. q8 Y. _! Z
    
13.             end+ Q1 i, s( H7 J2 {
    
14.      end. m0 H/ c- r1 Q0 F$ \\" {& }
    
15. a(i,:)=[]- Q3 c% f9 x5 H$ [. ]8 Q
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
改成这样的就可以运行了

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

madio

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

1. c=[1 1 0 0 0 0 0;7 k/ Q$ e3 W1 o4 F
   
2. -1 0 -1 1 0 0 0;7 G( B5 A; y\" G0 G0 K
   
3. 0 -1 1 0 1 1 0;- r* c% q' v9 p( y$ ?- x+ ?\" ]
   
4. 0 0 0 -1 -1 0 -1;
   : L6 b- H9 X1 g) O. t
5. 0 0 0 0 0 -1 1 ]( y2 N! `# |: g8 }% T- x# n3 [\" @' T7 g$ L
   
6. for i=1:5
   7 j4 @$ s2 n5 e, _+ k
7. a=c;
   \" O7 R( s5 Y* T: \; A1 S, a: j
8. b=[];
   ' W% I) U# Z! l. S
9. del_col=[];
   3 D3 L5 L# \: X$ D) y- M. |: K
10. b(i,:)=a(i,:)
    6 P- K1 K$ q9 C\" Y
11.      for j=1:7: K, a\" t' t1 h
    
12.         b(i,j);' M* G+ k; O% A) i
    
13.          if b(i,j)~=0
    9 V% `: q# k2 Z: Q  f5 g& z) d
14.            del_col=[del_col,j];
    # P; R\" w# b$ n' b7 O+ P  L
15.          end! I. u' [- o9 \9 p
    
16.      end
    - {9 W  p. f5 B+ G1 L- t9 G; ~
17.      a(:,del_col)=[]8 I, q( a\" Y5 h) Y+ X
    
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
  M- i& O8 w/ S/ z# D( @你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！