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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   ' C( Q8 V- t5 R' o+ I* o: p/ E+ r
2. -1 0 -1 1 0 0 0;
   * t$ I  P4 }\" M7 N, O0 H9 U
3. 0 -1 1 0 1 1 0;8 m( S# C% |7 I+ x% {* ~0 B, Z
   
4. 0 0 0 -1 -1 0 -1;
   7 |' @, y2 {4 q+ O\" x3 D
5. 0 0 0 0 0 -1 1 ]+ Q. k. E0 M5 O1 e/ z! a: }0 A
   
6. for i=1:57 [  j- Q9 d; k$ m  c; y) U: H
   
7. a=c;7 n& l6 S+ A% p! x) \  |% Q
   
8. b=a(i,:)* k$ O( X8 O\" @5 y4 A0 |\" T
   
9.      for j=1:42 ^/ {5 S5 i7 ?) k* I\" S7 }- R
   
10.         b(i,j);
    ' r8 @* V* x  t0 }
11.             if b(i,j)~=0
    ) x% t! t' C5 x* x% M1 h2 m4 s
12.                a(:,j)=[]3 B5 c; }. f  v  ]
    
13.             end( \+ f, T5 M% `  b1 `! M: w4 H
    
14.      end  ]5 `' K2 O' q6 f' |7 Y4 _
    
15. a(i,:)=[]
    ; {, h( K4 H5 q; ~/ X
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;
   ' x6 h/ u# I7 u: r, Z8 U4 M, u
2. -1 0 -1 1 0 0 0;
   ! O9 X! M4 y6 |: U% W
3. 0 -1 1 0 1 1 0;
   6 [8 M8 \$ T  |1 j$ D! H
4. 0 0 0 -1 -1 0 -1;0 W5 H! J$ ^1 f\" z% `
   
5. 0 0 0 0 0 -1 1 ]; n  {, B6 @5 s: S
   
6. for i=1:5
   & O  u  X6 I* s3 F7 Z2 i
7. a=c;
   2 k, o! ]2 c3 E8 y- [) Y2 S
8. b(i,:)=a(i,:)! m& g* d; d. ^
   
9.      for j=1:47 N) i( F. l% q6 z
   
10.         b(i,j);
    ) e: m! q( }9 }% \
11.             if b(i,j)~=0
    $ ^2 h( F6 S1 D\" u. S* j' d/ f4 y
12.                a(:,j)=[]% p9 k2 N  v1 N4 Y
    
13.             end
    0 D$ [8 D3 z! R! y+ g: I9 J
14.      end
    ! u1 h0 s# G  F
15. a(i,:)=[]6 l- p, m\" Z, x* n\" c
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
! g- B. {6 G$ T; c! z* M改成这样的就可以运行了

  Y+ D1 s( L; Y1 H1 K不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   \" G  ~: c. c, O+ |$ L
2. -1 0 -1 1 0 0 0;1 u: k: K7 d: ^
   
3. 0 -1 1 0 1 1 0;9 Z1 F5 Z$ P) k7 ^2 ^0 M. E
   
4. 0 0 0 -1 -1 0 -1;8 @# g: A5 d- z
   
5. 0 0 0 0 0 -1 1 ]
   ; O3 y) A, o$ ]9 |\" p% q9 X. s
6. for i=1:57 \5 S& Z  E; w+ I( Y( L
   
7. a=c;
   + J! D% C' T+ k3 u, W
8. b=[];+ C  y3 a4 R/ C, q9 p2 u2 Q
   
9. del_col=[];
   : d; Q# P! J: Z
10. b(i,:)=a(i,:)0 B7 d9 M7 c, O6 k
    
11.      for j=1:7* \) d1 z0 C/ V$ W* P8 f& P  i
    
12.         b(i,j);
    : N1 _+ ^\" n& w. n, }, U! E1 H- \
13.          if b(i,j)~=0
    $ C, y$ W  C( U5 @, m
14.            del_col=[del_col,j];
    \" W. e$ n, c/ F: ]) h
15.          end
    $ G3 x0 }6 L8 [8 _1 o1 \0 _% w# F
16.      end- w( [! ~7 n; ?( {- |+ z: F
    
17.      a(:,del_col)=[]
    ' y8 {3 W' y, a6 \0 h' }1 _8 v
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
1 z" W* L* c. [3 C% T6 M你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

2 m' u' {( ~. ~) k# l9 T/ |2 g可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！