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

xiangpeng395

1. c=[1 1 0 0 0 0 0;' q+ e\" G' s6 n: \6 u
   
2. -1 0 -1 1 0 0 0;5 ]7 @( v/ o% X& @! z3 [7 M- w
   
3. 0 -1 1 0 1 1 0;
   ( t! ]( O3 d, f
4. 0 0 0 -1 -1 0 -1;' O: Z7 }* g! n
   
5. 0 0 0 0 0 -1 1 ]
   % z% H$ i4 `5 R/ A
6. for i=1:5$ d6 O9 S9 a  h# ]6 d# e
   
7. a=c;
   ; I+ i0 M1 V; g8 e! U( x
8. b=a(i,:)
   ) H* s2 b. d( D; T& E/ \0 A5 p9 q
9.      for j=1:4) E& C+ f7 s, u  f$ \3 D
   
10.         b(i,j);$ P\" k4 u# y/ _' S2 ~- N4 \8 t
    
11.             if b(i,j)~=0+ b  E. O& _2 Q* u0 J$ }
    
12.                a(:,j)=[]
    # M2 l8 V4 Z2 T
13.             end4 d  A) m\" ]% b1 h+ y
    
14.      end: f5 T$ S% r' H3 C
    
15. a(i,:)=[]
    4 G+ A9 [9 Z* Y  k6 \, ^
16. end

复制代码

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

% [0 x/ W& m; W! o% h5 W. m& ^

madio

1. c=[1 1 0 0 0 0 0;+ L) N& @( w0 E$ P
   
2. -1 0 -1 1 0 0 0;
   : ?. c, d  w! W' }
3. 0 -1 1 0 1 1 0;7 p5 c/ K) s2 B( Y
   
4. 0 0 0 -1 -1 0 -1;
   # ~* m# t  u6 P( y  T8 B7 }+ A
5. 0 0 0 0 0 -1 1 ]/ W$ i( l3 _2 v- t
   
6. for i=1:54 j7 |5 e8 r, m  L6 q2 M
   
7. a=c;
   5 t& \9 l( {/ A
8. b(i,:)=a(i,:). V4 C' W$ J1 G# e* O- {, R
   
9.      for j=1:4! F4 J/ F8 q) @! g; q
   
10.         b(i,j);4 X\" Z# H$ r$ O
    
11.             if b(i,j)~=02 v+ s/ j. H. n9 G
    
12.                a(:,j)=[]5 L' L4 c7 c& r; v' Q
    
13.             end
    ( [# @1 E# Q* a# [5 _6 J
14.      end  E. A8 s) J\" l+ o+ T- q
    
15. a(i,:)=[]6 u# b- u; Y9 T
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;! y5 o1 V+ f0 H6 Q8 k
   
2. -1 0 -1 1 0 0 0;8 C3 L9 O4 \4 g4 Z- W
   
3. 0 -1 1 0 1 1 0;' i! U1 @1 u# [; p; J! ?1 o( j
   
4. 0 0 0 -1 -1 0 -1;
   6 ~$ \6 g: |\" l  V. T( v
5. 0 0 0 0 0 -1 1 ]- L/ S+ j1 ?6 b& U
   
6. for i=1:5% `- H( j; S% f& i3 U, \; G
   
7. a=c;) O9 @, r2 D$ `# O
   
8. b=[];9 p8 }; R  W+ R6 |2 m  k# T4 I
   
9. del_col=[];- l# o0 m9 B5 X
   
10. b(i,:)=a(i,:). j5 B; |+ ]9 t( T  U, U' U; X/ h: v
    
11.      for j=1:7
    ' |; C0 a3 I/ ?+ e* j+ m2 B
12.         b(i,j);
    ' j' E) s  W# i' Q1 G6 |/ |
13.          if b(i,j)~=0; G9 D7 O. @. }/ }8 d
    
14.            del_col=[del_col,j];, ]2 a! }) ?) C! G, j) D
    
15.          end7 b! V. e  C, D+ s$ X
    
16.      end
    4 e- J# t$ V; V3 g% |7 d! F7 F
17.      a(:,del_col)=[]
    5 s' U( ]; o; e' i
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

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

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