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

xiangpeng395

1. c=[1 1 0 0 0 0 0;! X5 j4 \. V  v
   
2. -1 0 -1 1 0 0 0;
   $ m4 E- W% l; T7 A$ H
3. 0 -1 1 0 1 1 0;2 l0 O4 ?9 S1 L( T4 e; \( I9 o
   
4. 0 0 0 -1 -1 0 -1;
   8 H8 y5 f+ q7 X; q: L
5. 0 0 0 0 0 -1 1 ]4 |) U% w: V, c. h  h1 F
   
6. for i=1:5
   \" \6 ?4 e& I2 X) ^5 D
7. a=c;' a# e$ v; c4 R+ W\" ?
   
8. b=a(i,:)
   8 L. R, N* i3 F
9.      for j=1:42 \4 a. d9 J+ W- U
   
10.         b(i,j);
    - O. m1 W, d/ E7 G2 w: i2 e4 @9 [
11.             if b(i,j)~=0; M+ F, z% @+ L2 H- T5 K
    
12.                a(:,j)=[]  V4 k6 {0 J: P
    
13.             end0 U$ i4 ~5 [  f& L. ^2 l
    
14.      end4 V) E! {: I- @# y3 w3 a+ X- h\" s
    
15. a(i,:)=[]
    6 }$ n\" |% X9 M) N# x( j+ K
16. end

复制代码

为什么总是超出维度啊，
/ l7 O1 V, C+ f# ?5 S2 {: [6 J我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
9 U8 n" L  M4 K* Y
; C; z+ p; U$ |7 I1 w

xiangpeng395

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

5 D7 ^2 U3 B  K! [. r/ H) R4 K, _& Z可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！

xiangpeng395

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

好的，谢谢了，

madio

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

1. c=[1 1 0 0 0 0 0;5 [3 Y' j4 V+ j5 q* B
   
2. -1 0 -1 1 0 0 0;0 W. I# W1 ]) t/ Q+ E7 L; {
   
3. 0 -1 1 0 1 1 0;
   , B  a0 U7 g9 Z# g, d6 X. d& k
4. 0 0 0 -1 -1 0 -1;5 p! ~) e4 `0 q; n9 C
   
5. 0 0 0 0 0 -1 1 ]- E- r/ s2 g/ |. v. D
   
6. for i=1:5: Q( s# F- Z4 z0 ?3 F
   
7. a=c;
   $ Q8 E; d. }: I( p6 I5 V  O' Q\" R
8. b=[];
   1 [- ^) e# E. E) I, C
9. del_col=[];
   / E1 T* d3 \5 @  t\" Z0 C
10. b(i,:)=a(i,:)
    0 z& t1 Y( {7 T2 d# W4 ~9 ^$ v6 G) j
11.      for j=1:77 P$ Z' H+ w0 u) M) I
    
12.         b(i,j);! q& U7 I2 N. Q# y
    
13.          if b(i,j)~=0( M* t6 L1 P$ V
    
14.            del_col=[del_col,j];/ m. ~. x\" v* f) m! d( h; J) R* \  \
    
15.          end
    5 H7 j- V+ Q4 G( i
16.      end
      m* e/ X6 `4 O& _9 V9 L1 D2 B
17.      a(:,del_col)=[]
    # L, w) U* u9 T\" B: o8 h
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 20:26
3 t+ a4 o6 B; J  I0 J  i. o改成这样的就可以运行了

2 c7 Z) u: f# g' E0 g- C2 X不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

1. c=[1 1 0 0 0 0 0;! n, g: [! K; D3 l( C
   
2. -1 0 -1 1 0 0 0;
   ; }\" _2 p# {, B( M! }/ T
3. 0 -1 1 0 1 1 0;
   ! P) U- N( S; {2 }7 Y- o9 g
4. 0 0 0 -1 -1 0 -1;0 H* g5 V\" n( d2 `( }2 d# {
   
5. 0 0 0 0 0 -1 1 ]
   ' O4 J\" [* \& O* D# o3 S
6. for i=1:5
   5 ]2 q2 S# \& @0 ^6 q
7. a=c;% Q; A3 [5 `8 S# }. r
   
8. b(i,:)=a(i,:)% z3 D' ~4 B& i\" Q* i9 `9 r9 G
   
9.      for j=1:4, y! f8 C/ S, Y( K+ f! d, r
   
10.         b(i,j);( m1 }' l) y2 R  f$ h/ }
    
11.             if b(i,j)~=0; Y+ N  X; X9 E/ [  \$ c, b
    
12.                a(:,j)=[]( ]  }/ H- K- v2 q- c5 w9 h
    
13.             end/ P2 o\" t. x' B! x! {5 N- o3 c
    
14.      end
    ! Q& ^4 j: l7 O+ Z; b9 }/ p
15. a(i,:)=[]
    9 W. S5 _4 e0 U) h. x3 c9 C
16. end

复制代码

改成这样的就可以运行了