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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   / v! o5 }8 L! j6 B# I1 L. M
2. -1 0 -1 1 0 0 0;. F9 m0 `+ K5 R0 B6 r- o
   
3. 0 -1 1 0 1 1 0;
   # R9 V0 Z  x6 X1 R9 w6 m
4. 0 0 0 -1 -1 0 -1;! R/ E- K* c5 u\" g\" {: c* |
   
5. 0 0 0 0 0 -1 1 ]
   + a7 x: }& ]/ Q1 Z
6. for i=1:5
   ; {& ?: Y4 x5 K; ^$ i; C
7. a=c;3 T1 O2 H) a3 w8 n& I
   
8. b=a(i,:)
   5 M2 m; j0 K; V) r  s. p( y0 y
9.      for j=1:4
   ) m) B/ I\" z7 d$ h! H, v& G
10.         b(i,j);
    , }* f5 c5 h- l0 w\" V* j: n( |8 F3 v
11.             if b(i,j)~=06 |, {/ f% M5 g% C
    
12.                a(:,j)=[]0 c( [+ [$ C* ?
    
13.             end2 t\" q9 x& F: T
    
14.      end
    . H8 S) N$ L9 y6 g3 A$ L+ J0 v
15. a(i,:)=[]5 n' s/ b/ x6 `: E5 O
    
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
3 N% m7 }7 `/ J  z$ @: p+ l1 F
  W# l- x# ~9 ]" @

madio

1. c=[1 1 0 0 0 0 0;; e: ]- T/ H4 `
   
2. -1 0 -1 1 0 0 0;# |) N0 q7 q7 o
   
3. 0 -1 1 0 1 1 0;0 @$ C0 x3 ]* g# ?! ?/ ?
   
4. 0 0 0 -1 -1 0 -1;7 ?# n$ U\" V: p! \( a9 Z
   
5. 0 0 0 0 0 -1 1 ]
   ( K( Q& u8 l7 d1 T
6. for i=1:51 Y9 C/ Y/ H6 B9 l. V- H9 F( ]
   
7. a=c;' W  Z4 S0 R9 E+ \- X' {
   
8. b(i,:)=a(i,:)
   0 z0 ]$ s& p0 K' w; O: u1 i- S
9.      for j=1:4* }7 t! s5 ^\" R- j$ \6 w( y9 C  B7 h
   
10.         b(i,j);
    \" A2 @! b2 L6 }3 K
11.             if b(i,j)~=0
    . ?8 c8 N( T- _- W4 A3 _+ b
12.                a(:,j)=[]
    6 I# x, b6 w6 ?
13.             end, V2 w$ x3 @, g* e
    
14.      end
    + V9 Y( v3 z4 n2 V
15. a(i,:)=[]
    - C4 W3 @9 w: `2 H# i
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;0 Y9 \- |# R4 m% j3 k
   
2. -1 0 -1 1 0 0 0;
   3 ?* D3 ?# B  e; y6 a( a0 E$ c* {
3. 0 -1 1 0 1 1 0;  ~5 i! T9 @6 N/ W9 ]' s/ C) w
   
4. 0 0 0 -1 -1 0 -1;7 Q8 Q% l3 B7 V, M% d
   
5. 0 0 0 0 0 -1 1 ]* {% A1 c% _& q+ N4 Q
   
6. for i=1:5# z: M\" O' t5 N
   
7. a=c;8 w. d3 e, N$ f0 V/ O! ]' C. I
   
8. b=[];
   - g# X1 V# R& v, A1 o- ?
9. del_col=[];
   3 R; k8 L% A: O0 _! @$ ?+ V
10. b(i,:)=a(i,:)5 s. o5 a+ C# ]# k  S
    
11.      for j=1:7
    1 Z8 E& g) @$ K& n) y  \! `
12.         b(i,j);
    7 |; m' j8 z1 j9 w+ \
13.          if b(i,j)~=0
    0 @$ q. s( R. x1 R, S' Y' b
14.            del_col=[del_col,j];& M  Z, S! \' I5 t, `
    
15.          end
    * c  `4 i$ Q& a' G# k/ p1 q/ ]# `
16.      end$ S\" I4 W* q4 Z) `/ L( K/ V  I
    
17.      a(:,del_col)=[]
    ) u  W* W8 f8 I4 k' z
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
; \2 q6 c( r5 x- M& b你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

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

& i. V# h; Z0 ^# Y/ a! w$ q可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！