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

xiangpeng395

1. c=[1 1 0 0 0 0 0;: q$ L( C% G/ Y5 X; e
   
2. -1 0 -1 1 0 0 0;
   2 z: V2 A. Z\" A& X
3. 0 -1 1 0 1 1 0;* a7 o5 G7 d( n
   
4. 0 0 0 -1 -1 0 -1;- f5 P6 Z0 Y/ k9 m, J
   
5. 0 0 0 0 0 -1 1 ]+ {+ P! ^$ \/ H5 N5 G6 u
   
6. for i=1:5
   # Y8 n5 |, u4 |5 [- i7 f
7. a=c;
   5 U/ ~& s' @+ f- \9 b! w\" H( `
8. b=a(i,:)3 F  M6 V\" h) Q1 Y
   
9.      for j=1:4* k- t* p0 _, E! K/ }2 ?# y* @( b, i
   
10.         b(i,j);: J' p7 J\" f7 {* u6 ~\" J
    
11.             if b(i,j)~=0
    # e, g! Y' I& ^
12.                a(:,j)=[]
    $ Z8 q# ^5 F) _: W3 z. [
13.             end
    ) z, |7 d& ?( m3 ^2 _5 V; j
14.      end( C9 n\" `; E5 r! |1 `: H+ P
    
15. a(i,:)=[]
    ) ^6 _: G2 B0 D; c! O* p
16. end

复制代码

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

3 U  t1 b$ l$ O/ }& d* D9 ?, x

madio

1. c=[1 1 0 0 0 0 0;
   - D( J2 M$ o1 X- \\" _- X7 u: H
2. -1 0 -1 1 0 0 0;
   $ w& |+ }) d; }# Z$ R: c
3. 0 -1 1 0 1 1 0;* {; u: m9 |, F9 f
   
4. 0 0 0 -1 -1 0 -1;
   % s! X5 B4 ^& f% Y* U) \
5. 0 0 0 0 0 -1 1 ]
   + @9 a5 X& j\" D$ t+ w5 x8 t* _. b5 J
6. for i=1:5
   5 f( O5 l3 V, e5 [; {4 K
7. a=c;
   , w) b/ D9 [2 W% X! H4 z# k
8. b(i,:)=a(i,:)' s8 q4 j$ k: H# \+ d
   
9.      for j=1:4
   7 R& p/ Z* j& M, G% u7 Q8 [( }
10.         b(i,j);
    & c- N7 _1 G+ M$ T
11.             if b(i,j)~=0
    ; Z! t2 k7 g. h
12.                a(:,j)=[]) |3 e& p8 T& b
    
13.             end
    \" e* A, R5 H6 Z$ o& X; l
14.      end! b* y3 l$ p7 f8 j/ ?
    
15. a(i,:)=[]! {- }, I& R+ G' N
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;
   5 M3 O# V) _1 k8 I! u. M
2. -1 0 -1 1 0 0 0;* q: N& Q\" z. ?8 {
   
3. 0 -1 1 0 1 1 0;/ d0 o& c; f; H+ |
   
4. 0 0 0 -1 -1 0 -1;
   $ M$ x- z\" |6 j7 A4 _! P
5. 0 0 0 0 0 -1 1 ]
   4 g/ N2 V\" }2 x- e* a4 }
6. for i=1:5% u7 K' r5 T. X' [* H
   
7. a=c;
   + y3 M3 R8 r6 N8 O
8. b=[];
   1 B6 c7 B: O; h5 H
9. del_col=[];) D; V+ F9 p2 R  l- Y
   
10. b(i,:)=a(i,:)
    6 }7 v# ~: C! _- }1 v; Y+ s1 x
11.      for j=1:7
    # r; X, D7 n: b\" R
12.         b(i,j);7 A# V% }0 m( D! _6 _\" U. ?
    
13.          if b(i,j)~=0
    0 L; ^, `& Y' g\" z* [( D* d6 x0 j0 X
14.            del_col=[del_col,j];+ ]( @6 E, M0 p1 I
    
15.          end
    \" i) }6 Y# \9 C7 u; r
16.      end
    . b3 s( o; j& u; X
17.      a(:,del_col)=[]
    - O( M8 ]  B; j/ b  ]
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

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

) J  v3 q4 I0 @5 ~4 r% W) b7 N可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！