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

xiangpeng395

1. c=[1 1 0 0 0 0 0;' n& ^8 f$ S3 V, ?; {
   
2. -1 0 -1 1 0 0 0;4 B- N9 @' ?) c2 D2 L% C. ^4 J
   
3. 0 -1 1 0 1 1 0;
   8 T\" e) I1 x4 q! A1 ?
4. 0 0 0 -1 -1 0 -1;
   2 ?6 v) ]- v2 q6 V/ w- t
5. 0 0 0 0 0 -1 1 ]2 N6 c\" K6 K4 t\" W, t( y0 R2 A
   
6. for i=1:5$ K2 c  C! f8 ]$ K8 r
   
7. a=c;' T( I1 n5 Z+ o  S
   
8. b=a(i,:)
   3 F' i4 \5 w2 ~' L+ Z! ~, z
9.      for j=1:4* N, w5 j( d\" C* s7 m( Y
   
10.         b(i,j);
    - Q' x' I/ S5 G* `
11.             if b(i,j)~=05 P: \+ E* {; h( N0 B  U' \
    
12.                a(:,j)=[]0 y) A\" ?4 t3 e: E5 p
    
13.             end* u* O+ b: u, P
    
14.      end9 e9 A3 Y% n! c( L
    
15. a(i,:)=[]( U- H; P: P  X/ g
    
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
3 f9 s, c2 w0 j: h
2 @' Q  z( {5 ^6 O7 f

madio

1. c=[1 1 0 0 0 0 0;
   , f6 d/ G$ L% u0 {/ }- B
2. -1 0 -1 1 0 0 0;
   8 ?9 r4 ^5 ^7 @* r* |$ O$ R
3. 0 -1 1 0 1 1 0;$ s  n* Z5 M6 E8 i' y% x
   
4. 0 0 0 -1 -1 0 -1;
   % m/ o( W# v9 \; j9 u\" K4 R
5. 0 0 0 0 0 -1 1 ]* c' b1 \- Y  [0 l\" y4 W$ d
   
6. for i=1:5
   2 Z$ t1 V5 M8 u5 m
7. a=c;; U! ?4 _  z( P$ q. {
   
8. b(i,:)=a(i,:)
   ) I' b5 s- V* A; X2 S\" n
9.      for j=1:4
   : ]/ \+ k6 O* h) |
10.         b(i,j);
    2 C' f1 }2 A3 R, n
11.             if b(i,j)~=0! a) B: m/ W$ s' ^6 F( w& `5 E
    
12.                a(:,j)=[]
    & t\" F: {* u( t. Z$ {8 d- w% I
13.             end$ a; c5 g* N& c% U# O+ }% B
    
14.      end
    5 ?: |1 P6 z5 V# j/ ~
15. a(i,:)=[]0 d; e$ p  r, ?) b, h
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

6 s2 W  q0 D7 }, }: a  g$ y不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;$ u. W6 D% |& b7 A3 U2 w
   
2. -1 0 -1 1 0 0 0;9 K, ^\" D( ^2 B  n3 M
   
3. 0 -1 1 0 1 1 0;# [* T+ C# z+ A! b* u
   
4. 0 0 0 -1 -1 0 -1;8 q2 s+ x) q' b+ P) C
   
5. 0 0 0 0 0 -1 1 ]3 y. U$ o# t- Z9 g9 ^' ]( ~' _  \3 M
   
6. for i=1:5
   ; x+ l: c: a0 ?2 m' m9 v& F# Z& t
7. a=c;
   8 b! K  K) r& N( U9 k
8. b=[];
   5 o5 J6 t' t) F+ s
9. del_col=[];
   ! h) h8 i2 |; l' Y; Z8 C0 Q. Y. ]
10. b(i,:)=a(i,:)# ~0 ~4 M1 B) r2 Y7 O  l
    
11.      for j=1:7
    9 @4 `8 S% h+ i0 z' e- W$ r
12.         b(i,j);
    ; r- B\" w% W: K$ H
13.          if b(i,j)~=0
    0 r\" Y  o/ O& V2 [4 x
14.            del_col=[del_col,j];! x7 d5 }+ {  E' G7 w; N0 O$ E
    
15.          end0 H0 M. |) J0 ~, o
    
16.      end
    ' w2 y6 W) U: @& O: i# n3 U, W; u- h
17.      a(:,del_col)=[]
    7 X  g  L1 p; m, _( z
18. end

复制代码

xiangpeng395

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

3 |% g- p1 V: v" R: J* a( [' i好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
) r* S9 f9 @6 S6 l; H你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

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