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

xiangpeng395

1. c=[1 1 0 0 0 0 0;3 P3 u/ }6 n\" H, e
   
2. -1 0 -1 1 0 0 0;4 E8 ?\" O( I6 c6 l. t! Y, S4 v
   
3. 0 -1 1 0 1 1 0;  U2 ]9 b9 M$ z: ^  Z3 Y4 W; m5 w
   
4. 0 0 0 -1 -1 0 -1;
   6 x2 \; {, K- Z  X
5. 0 0 0 0 0 -1 1 ]! x& J\" Y4 j  p4 g. K' v/ B
   
6. for i=1:5
   ! R' {2 G6 N' c* A# @\" `& f1 `3 z4 N
7. a=c;2 z, u! z8 J6 l( n9 j  G4 g4 w: \
   
8. b=a(i,:)
   : T$ |0 m6 U, b5 ?: E3 E* U: ^: L* U1 }
9.      for j=1:4
   9 }) M4 _4 i* G0 T: x
10.         b(i,j);$ h+ C. T  B1 p9 k# V
    
11.             if b(i,j)~=0
    4 U! Z! _! ^0 Y5 o' _( J1 y
12.                a(:,j)=[]
    ! U: d6 {7 l& x3 b
13.             end# z) N; N( l4 \) o( J
    
14.      end- h- B4 C! ]) Y, k4 F- A. Q% q
    
15. a(i,:)=[]
    ( V3 l' p' v3 o6 X' C
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
' k4 g; k: z5 W0 E& |# E
. D3 ]5 {, w7 g2 u6 R7 M% T" o. K

madio

1. c=[1 1 0 0 0 0 0;4 b( C5 I, O4 d5 p( E; h3 z
   
2. -1 0 -1 1 0 0 0;
   + `2 q+ C: q! I0 E0 \
3. 0 -1 1 0 1 1 0;
   , S1 E  P4 Y' H1 ~. t) f! {
4. 0 0 0 -1 -1 0 -1;
   1 v9 q8 Z7 j! D0 U$ |3 X\" F4 z
5. 0 0 0 0 0 -1 1 ]5 X! b0 W) ^8 b# o. b0 B3 Y8 }
   
6. for i=1:5
   ; \0 t/ b6 U. o& ]' L: J
7. a=c;5 r$ n' a, c/ w0 p  z5 Z7 l: O, Q; c; U
   
8. b(i,:)=a(i,:)
   7 f1 r$ Y; C' G$ d! J) i
9.      for j=1:4( p: `* |. \4 J\" q1 B* o
   
10.         b(i,j);
    9 A3 H+ G7 x- A( Q7 w/ a
11.             if b(i,j)~=0
    + ]9 F0 H7 y6 c9 f& X7 d! u
12.                a(:,j)=[]5 q7 B8 L- h: W
    
13.             end
    + U- c9 d, i( g3 u/ ~: O8 \) p
14.      end
    - R7 d( S3 K0 U: ^  j
15. a(i,:)=[]$ `2 f9 R3 h  [
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
8 O, B; ]1 Q* W4 r改成这样的就可以运行了

, f; U' c, e8 [7 l4 A% P不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   & v3 `8 k# @- u# j1 Y( O, R# d\" F7 S
2. -1 0 -1 1 0 0 0;
   1 J\" Y3 ]5 y! A- ?: q9 M, b. l* E
3. 0 -1 1 0 1 1 0;3 F. E- A; G$ p* x/ l0 ~  Z
   
4. 0 0 0 -1 -1 0 -1;
   & k* @! _+ |% W+ i& P2 E7 Q
5. 0 0 0 0 0 -1 1 ]
   6 X/ b% D3 E! L' y
6. for i=1:5  w\" n: e3 S: P& X9 ^
   
7. a=c;
   + W, S7 ^* N, S. @7 K
8. b=[];4 U( s\" r6 L& _/ [. G
   
9. del_col=[];4 _8 f\" T( F, ]7 U4 J; |3 [0 @
   
10. b(i,:)=a(i,:)+ j' n, U8 n% C* o
    
11.      for j=1:7
    ( e4 b( \' R! p  X1 M# D4 x
12.         b(i,j);
      j+ H1 M9 f4 ^( o6 e5 W6 W4 H0 \
13.          if b(i,j)~=08 T3 P; E4 [7 p  ^& s0 R% c
    
14.            del_col=[del_col,j];
    & C( G% H4 B  X2 C, a! Q- I& F
15.          end
    0 a1 \, |8 m8 H! f1 ?/ ?1 U
16.      end# F\" Y* l% T7 {
    
17.      a(:,del_col)=[]
    3 b8 D5 C7 T8 t4 n4 ?' K1 a0 Z; Q6 u' k
18. end

复制代码

xiangpeng395

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

# \5 {0 Q" X! D% V' D  |% v+ z4 g好的，谢谢了，

xiangpeng395

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

4 I6 [. n" w/ x可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！