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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   ) {8 {( V0 T\" x. j; t1 O7 j
2. -1 0 -1 1 0 0 0;' J\" e7 e2 m$ n) `+ Z4 t( U7 D
   
3. 0 -1 1 0 1 1 0;7 X- ~4 i& _6 l7 K
   
4. 0 0 0 -1 -1 0 -1;
   ' m; {2 O. ]' |% [- [- C
5. 0 0 0 0 0 -1 1 ]
     W9 S+ p; g/ o: \
6. for i=1:5# i\" x- O% q4 m% t# F
   
7. a=c;
   ' r' {$ M* o( g7 ~* o
8. b=a(i,:)8 a/ e6 d- r3 G1 {, c+ m# r( o; W# G
   
9.      for j=1:4
     ?. t* j: e/ ?
10.         b(i,j);: I$ W' r; R. h4 [\" Q) d
    
11.             if b(i,j)~=01 h9 T6 v/ L! z& s4 J. d
    
12.                a(:,j)=[]1 F1 v& e1 Z9 `
    
13.             end
    & }& n0 g$ k+ s\" p: q: ^0 ?
14.      end* \8 h0 X$ J) ]) ~1 A
    
15. a(i,:)=[]
    + E* R9 }: A8 v! S7 S* D
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;4 `! G- [! O8 g2 k5 s+ c
   
2. -1 0 -1 1 0 0 0;) A' H4 L3 p. q! f
   
3. 0 -1 1 0 1 1 0;
   - x& H7 x$ Q. T8 J% {& I7 p1 Z
4. 0 0 0 -1 -1 0 -1;
   # L+ O& t7 x, v\" d- ^, f
5. 0 0 0 0 0 -1 1 ]' G& O# L* w2 E2 y: b& c
   
6. for i=1:5+ {2 Y\" [/ ~+ s
   
7. a=c;& z0 K1 Q1 n/ `, r. j4 b
   
8. b(i,:)=a(i,:)0 f. W# P2 F. Y+ s' {6 F
   
9.      for j=1:4
   \" U0 V' _9 j, A$ d$ t# u+ Z$ d* M
10.         b(i,j);+ ~\" Z( _/ n( t6 ~\" S+ t: z' V
    
11.             if b(i,j)~=0
    6 }* {; a9 X3 E4 `
12.                a(:,j)=[]
    3 N% n\" ^+ e, L5 E( X. v( Q
13.             end
      A$ R* T. {4 a% G+ W6 y& Z
14.      end
    . C  j1 a5 H3 }9 d% p
15. a(i,:)=[]' S' P; j2 ?& U  t
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

$ e* C5 G6 d& X0 ?9 L$ S+ c不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   : V+ T& S; Q: d9 a. z$ |! B. j
2. -1 0 -1 1 0 0 0;8 h$ z$ Q3 n% i
   
3. 0 -1 1 0 1 1 0;1 i$ x* \3 d* @: i; e: c$ Q7 y+ d
   
4. 0 0 0 -1 -1 0 -1;( M2 G- ?, K\" }7 N0 j\" G$ A
   
5. 0 0 0 0 0 -1 1 ]: o/ ]/ Z& W! N3 q/ C' a+ @
   
6. for i=1:5
   $ s1 Q# Z7 k) f% t7 O
7. a=c;
   & A; j5 U\" _0 m; T
8. b=[];
   $ u: z; l\" s6 T0 e
9. del_col=[];
     H! x4 E' r  h4 c  v+ G+ }1 I\" {
10. b(i,:)=a(i,:)# G; ?- E: L+ ^, H: H
    
11.      for j=1:7
    2 ?( \- r* V, X; M9 e1 u
12.         b(i,j);
    1 t: g9 Q6 F7 t7 y
13.          if b(i,j)~=0
    # A: T3 u' _+ b6 i, k
14.            del_col=[del_col,j];
    . {1 m* k: L4 o# V# q8 i
15.          end
    / z( G% W$ r: O+ M* t. ~5 U1 T# i/ q
16.      end* P2 p- t  R$ t+ B- M  k
    
17.      a(:,del_col)=[]1 Z4 R0 `  E; B1 a  R
    
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
; e  {7 @, i1 L6 R. p  I5 O# |你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
" x8 m+ K6 X: R! ?% d$ x  Q! M你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

9 D0 x  F) ?+ n* c可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！