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

xiangpeng395

1. c=[1 1 0 0 0 0 0;+ D6 ^$ S; S9 a) B* ~
   
2. -1 0 -1 1 0 0 0;
   9 h4 ^% n! {) B0 e& ]' c
3. 0 -1 1 0 1 1 0;
   7 q2 O5 e! w- N. U$ D
4. 0 0 0 -1 -1 0 -1;  k# {4 t; i( }1 g
   
5. 0 0 0 0 0 -1 1 ]$ s' I& U$ v3 |  |7 Y
   
6. for i=1:5
   7 K% K/ U7 r\" p+ \  Q/ |
7. a=c;% B: }- i. c3 o6 l+ n; h
   
8. b=a(i,:)
   0 e5 i6 C% B( b3 g
9.      for j=1:4
   \" z/ Q  P' Q\" d
10.         b(i,j);
    4 J9 k$ w& t) R2 C* ]  ?\" a) j
11.             if b(i,j)~=0, c! b4 J7 n$ ~7 U\" C. t' d% f
    
12.                a(:,j)=[]: Q; J5 i/ M\" c
    
13.             end
    : c5 G( T% C% @) V
14.      end- @! o7 F) H* f8 ^
    
15. a(i,:)=[]& A( {2 c/ f8 \8 V( N2 B* Z; Z! c
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;
   $ p0 _4 q1 t$ Q8 F+ S
2. -1 0 -1 1 0 0 0;! _. A1 M0 p/ s% S/ p+ n0 H2 o, ?! W
   
3. 0 -1 1 0 1 1 0;9 d2 k# B2 U( Q& }7 z5 T( j
   
4. 0 0 0 -1 -1 0 -1;
   . y, I/ a! S5 O7 A- I7 Q8 D  Z
5. 0 0 0 0 0 -1 1 ]
   # x0 S$ C7 B2 ?0 W
6. for i=1:5
   * W9 J! Q) `) x5 f# U( {( ]& T
7. a=c;; T8 O5 H8 t% N4 q+ z- A; ]
   
8. b(i,:)=a(i,:)
   5 ^: ^- E9 t4 Z6 {: E
9.      for j=1:4* ~2 D9 J0 W, Q5 k\" Y7 u
   
10.         b(i,j);) K, ^- S2 ^0 y7 i
    
11.             if b(i,j)~=0
    5 f$ r0 D3 E, K: c/ _
12.                a(:,j)=[]6 h9 L& H: W$ Y7 G! Q# V: h4 H& m
    
13.             end5 G, q% T6 ^! E7 v
    
14.      end& v8 d1 r- g  a, e. J
    
15. a(i,:)=[]
    ' g3 U( e# G9 h# q0 r8 ?* k
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;
   8 [8 x( E* S& b2 v8 e% u  a# ^1 Z
2. -1 0 -1 1 0 0 0;
   $ ?+ _* u6 O- @* h
3. 0 -1 1 0 1 1 0;
   \" I2 K' i+ p- y4 U
4. 0 0 0 -1 -1 0 -1;
   4 G1 K! |8 R0 L) x\" W5 j( k
5. 0 0 0 0 0 -1 1 ]
   ( ^4 ~6 o/ v/ [- @8 X) |8 }7 J
6. for i=1:5
   9 f& P7 m+ c5 I& e8 u( {/ [  e
7. a=c;3 S0 r! h# {; u2 o* d6 X& b$ k, o
   
8. b=[];
   , d, S6 d' H7 Q0 r3 I! w& G  @  E
9. del_col=[];
   5 ~3 d) W6 {+ s7 f
10. b(i,:)=a(i,:)
    6 A- s: x! A+ \: G% X% w1 t, I- w
11.      for j=1:71 [) d\" L5 K7 l- G5 {4 F
    
12.         b(i,j);! L9 _! y: I+ D
    
13.          if b(i,j)~=0
    3 ^( `7 Y& @+ c\" W
14.            del_col=[del_col,j];6 v) x: g3 b: V# ]; p\" ~) \. j- @
    
15.          end6 X% K: E) ?! w7 o8 Y2 R1 W( W
    
16.      end& @, S7 o% s( ]: r, x# F* V
    
17.      a(:,del_col)=[]
    \" z  T\" d\" P+ L+ u+ H0 o$ B4 \
18. end

复制代码

xiangpeng395

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

1 X7 ^$ ~& b8 l/ ]& Z' u* Q) v好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
! T) ]: c- N( C; u4 q* n你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

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