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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   0 S8 F3 h\" ~( c4 ?* o
2. -1 0 -1 1 0 0 0;  i4 [\" Y, j\" w/ I2 H4 p# j0 D
   
3. 0 -1 1 0 1 1 0;/ A; U$ m4 a: h4 }# A7 D$ Y
   
4. 0 0 0 -1 -1 0 -1;3 E* D/ O9 O* e, g2 S# }) y
   
5. 0 0 0 0 0 -1 1 ]( q/ U/ h( k; g# n' V9 `0 `$ i0 ?
   
6. for i=1:5
   / h3 T* y0 K3 i5 u1 |8 r, `8 O: g% W
7. a=c;
   ! Y  U. M) Q- e
8. b=a(i,:)
   % R( p3 M\" X1 ~
9.      for j=1:4
   / r, L0 I$ ]\" H) i
10.         b(i,j);
    ( d: Y! V\" x! A; T- h' |0 g\" @\" S
11.             if b(i,j)~=0\" P' M& Q' O0 Q7 j2 m) p+ A+ F) S
    
12.                a(:,j)=[]
      A/ g6 f, Z: j/ }7 M
13.             end
    ( O. B* k+ k' {# {: a7 B
14.      end
    6 G% ?% t$ ?6 `8 H5 `6 F
15. a(i,:)=[]: H- B$ Z3 I' x
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;% g) }: z# Z7 }. h- J
   
2. -1 0 -1 1 0 0 0;' Z! G* a2 J% ~, T
   
3. 0 -1 1 0 1 1 0;: A- r8 H+ }) _0 K9 W8 ~
   
4. 0 0 0 -1 -1 0 -1;
   9 l2 C5 Y; ?0 Y& p
5. 0 0 0 0 0 -1 1 ]
   : _- Y; x4 f  k: f
6. for i=1:5
   ( y( a9 z- p0 o. o6 N
7. a=c;
   $ J! U5 T; l- p7 J) @
8. b(i,:)=a(i,:)' w5 s: w- Z7 v% I
   
9.      for j=1:44 K5 _4 m\" r2 E$ n' I. z
   
10.         b(i,j);1 c: F1 O0 \) E5 k8 a/ d
    
11.             if b(i,j)~=0
    * ]% `* ?& A9 D; D9 p5 O
12.                a(:,j)=[]
    2 m- z# O5 {- \/ |9 @  ?/ I2 |+ n$ E
13.             end
    * |. F/ @% ?7 A8 Z
14.      end  S4 m7 d\" U* \* y- ]
    
15. a(i,:)=[]4 u4 s9 ]3 N4 A
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
+ {, f% {4 a3 O' O$ {0 [3 v# e# H改成这样的就可以运行了

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

madio

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

1. c=[1 1 0 0 0 0 0;
   : R) h$ ?! l' ~0 V% d+ C1 ]
2. -1 0 -1 1 0 0 0;
   6 ]  c- X; O) N- ]
3. 0 -1 1 0 1 1 0;  z. Q( y/ B  J8 n1 b
   
4. 0 0 0 -1 -1 0 -1;% i. }; |. ~/ l( B- }6 M% n3 L2 M
   
5. 0 0 0 0 0 -1 1 ]2 _' u  q( P1 M+ A
   
6. for i=1:51 }% g- F9 x8 m1 l. S' O  m) t1 _' l
   
7. a=c;
   0 q2 K2 G4 A  Y$ T( r: q) W
8. b=[];6 p8 K9 }; i; H0 o$ h; s3 D. E; J
   
9. del_col=[];
   6 X) n# x' C# x+ H
10. b(i,:)=a(i,:)
    , b9 {# U4 R2 {( T
11.      for j=1:7
    2 j; A2 {) M% m' W- ?6 H
12.         b(i,j);
    + f, w/ [2 A/ I6 I
13.          if b(i,j)~=0
    + \- e8 D7 @! O$ n: s
14.            del_col=[del_col,j];1 ?$ e( V/ E% @9 U8 Y3 \
    
15.          end9 ?( t, n: I6 M1 {% a
    
16.      end* K5 p4 h$ ~5 B5 K7 ]  X
    
17.      a(:,del_col)=[]
    9 Y1 |6 \5 y1 f3 ?$ X/ h
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

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

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