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

xiangpeng395

1. c=[1 1 0 0 0 0 0;' l8 I! A8 T% L3 m! Z# a* i$ G. J6 J
   
2. -1 0 -1 1 0 0 0;- ~2 d2 a* z; T
   
3. 0 -1 1 0 1 1 0;- h1 l  d6 d$ k; s1 C\" V( i
   
4. 0 0 0 -1 -1 0 -1;& F' h+ d, L$ d5 @3 I) ]+ B+ \
   
5. 0 0 0 0 0 -1 1 ]7 O4 k( I\" c5 w
   
6. for i=1:5) l8 @\" V' m) h0 H
   
7. a=c;
   ; a2 m3 h8 V2 i& f, u0 K! V0 n. N. j
8. b=a(i,:)
   \" s- s: A4 |) f) X4 S
9.      for j=1:4. n( t2 m, D\" V( p9 }9 d
   
10.         b(i,j);; Z3 ?( A3 a) n  e2 u. @) V( M
    
11.             if b(i,j)~=0
      e\" v\" Q8 h* c) p7 Z, U& a* \
12.                a(:,j)=[]  R) Y! x* P7 J' h0 V
    
13.             end/ f) F5 v9 \; {
    
14.      end3 Z; Z6 P' R+ `: X: d. c9 @. x
    
15. a(i,:)=[]5 J' g3 y7 z, c9 t. N
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;+ v7 I8 Y% |: \/ ^
   
2. -1 0 -1 1 0 0 0;4 Z* [% b+ p, a# r& ~0 o4 j* g
   
3. 0 -1 1 0 1 1 0;- Z; J/ f; N. v2 G- ]2 G\" T
   
4. 0 0 0 -1 -1 0 -1;8 j: R\" |3 a' w
   
5. 0 0 0 0 0 -1 1 ]0 g7 y, v) z, D( \2 H3 x3 i
   
6. for i=1:50 h7 g' k8 J, s' S/ _! a7 |+ I* Q2 N
   
7. a=c;& S! n/ t% {$ a8 F
   
8. b(i,:)=a(i,:). \3 |! P6 S# [  X\" f8 g
   
9.      for j=1:4( q2 e1 P2 n, F. {/ w\" X9 D: {9 q: |2 l
   
10.         b(i,j);$ t3 X4 W* w+ I
    
11.             if b(i,j)~=06 ^4 B8 J0 s. S
    
12.                a(:,j)=[]& N) s8 a/ m1 v5 g1 k0 L$ s  i8 h
    
13.             end- Y: F+ \3 Y* ^; E# l4 r# m
    
14.      end1 \\" S9 H  k: ^6 F. I5 |
    
15. a(i,:)=[]
    * I, @3 ]& g9 x4 {% M, m
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;
   . p& o& u6 M  c\" V& j, V1 b$ e4 c
2. -1 0 -1 1 0 0 0;% p7 g# n2 r  \9 w' G
   
3. 0 -1 1 0 1 1 0;0 k# S. t( h\" Z) p( p
   
4. 0 0 0 -1 -1 0 -1;/ a\" k, y' G: P# d* Z) A
   
5. 0 0 0 0 0 -1 1 ]
   3 u' O; ^, w5 _! L. u' r6 G) w
6. for i=1:5
   2 e* }7 y+ E: D) P: Z\" h# V
7. a=c;
   ' v* _* d, F2 n9 N+ y
8. b=[];
   % G6 U3 J# @3 _! k+ R
9. del_col=[];( L* \3 x) d3 C$ a+ S
   
10. b(i,:)=a(i,:)
    ) H- ?7 i' |* ~$ N! `
11.      for j=1:7- }. e# p5 N7 N
    
12.         b(i,j);
    ) Z7 R: V; P6 t8 f- j: d
13.          if b(i,j)~=00 b* k) [7 h8 `# r+ M5 D& L! A  ]
    
14.            del_col=[del_col,j];% J\" `& ^0 L$ G5 V5 y9 `2 r$ d
    
15.          end
    * M4 x7 M  [1 n& {
16.      end6 v5 G1 Z\" n7 C/ D6 Q+ ?# U6 V6 `0 C
    
17.      a(:,del_col)=[]
    $ V% \5 h7 @8 Q6 Z; e
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

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

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