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

xiangpeng395

1. c=[1 1 0 0 0 0 0;) X( x& l. ?$ n8 F* A3 p: C\" r
   
2. -1 0 -1 1 0 0 0;
   $ {\" F+ e1 S7 q' s/ `. |* h
3. 0 -1 1 0 1 1 0;8 q7 `: P3 K/ p
   
4. 0 0 0 -1 -1 0 -1;
   $ Y7 F, ]9 b2 ~6 ]$ W
5. 0 0 0 0 0 -1 1 ]8 O. U( X  q8 F
   
6. for i=1:5) F9 s& R. b3 n+ a' l\" F  {9 z
   
7. a=c;( l9 z. E- h% D7 T! f1 o; h
   
8. b=a(i,:)
   3 M0 x7 `+ g( k* Q  d
9.      for j=1:4
   ' O; }, r/ V& b4 k2 y
10.         b(i,j);0 q- u2 @2 D3 m9 ?% o
    
11.             if b(i,j)~=0
    ( J5 |4 {: f* D5 ?4 c% Z
12.                a(:,j)=[]5 X# {# j7 z2 U% q3 v
    
13.             end
    / ?( R7 x) L/ u\" g' w8 ^
14.      end- w1 `# }9 Q* X8 a  z8 k
    
15. a(i,:)=[]; H# o) l  ~3 _* d
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;/ \% N3 V8 V8 J& {0 ^. ~7 A9 Q
   
2. -1 0 -1 1 0 0 0;9 [' j0 `$ S) i  v# f1 m$ u  {; m
   
3. 0 -1 1 0 1 1 0;
   0 Q9 J, n) n5 k8 A
4. 0 0 0 -1 -1 0 -1;
   3 V5 V% F* \- q/ U4 R) w* @- P
5. 0 0 0 0 0 -1 1 ]
   9 C# ^/ o; `2 S. n
6. for i=1:59 |0 t# \/ }( h2 W4 _+ x  W- c
   
7. a=c;
   - m( ~! w/ g  b
8. b(i,:)=a(i,:)
   . h& I6 n5 i1 [2 b4 v( B
9.      for j=1:4) ~: T& D+ {2 t  z  t) Z
   
10.         b(i,j);
    * N+ [7 _  O7 b, J
11.             if b(i,j)~=0
    3 i7 J: l# d* t8 W
12.                a(:,j)=[]5 V9 a4 k7 Y( M) J0 Y
    
13.             end- b7 J# v) U6 w% o' q
    
14.      end5 C5 m- u0 e  d  q' c
    
15. a(i,:)=[]5 v% ]5 P1 w' h( ?
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;2 G- q+ k. a9 F/ i
   
2. -1 0 -1 1 0 0 0;8 O8 t! S0 H8 p( }; T- N& v$ x
   
3. 0 -1 1 0 1 1 0;' Y& G. T, w% b\" z
   
4. 0 0 0 -1 -1 0 -1;9 m4 O; f0 X; \5 k5 Q% K( P+ k9 H
   
5. 0 0 0 0 0 -1 1 ]2 o# `& O) P! M6 S9 Z) G
   
6. for i=1:5. p0 u# g% y  Z% e- w
   
7. a=c;/ [\" D# H4 F, v5 I  p
   
8. b=[];  @, d2 `+ y+ D, x% l+ p: U) v
   
9. del_col=[];
   ( k' Y; ~# E- s* x2 a  g
10. b(i,:)=a(i,:)
      Q0 I8 `( D2 L4 ]$ ^) ~
11.      for j=1:7# M8 j) I. j& m( U+ ~0 Y7 U; D+ {$ A
    
12.         b(i,j);
    $ t# _: X1 C, D* k' C% l2 z
13.          if b(i,j)~=03 p, ?# C2 g3 \7 c- G2 m% v$ }  h/ u
    
14.            del_col=[del_col,j];
    \" s( M5 A. @3 q. x1 l/ }% d+ x7 Q6 p
15.          end3 c* R0 }) Y* u# T0 R* m( l4 v
    
16.      end
    9 S3 B, W% v1 Q. O2 R
17.      a(:,del_col)=[]
    ' s' T2 w$ B, K8 m5 S- E; Z
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
8 j/ _& h; U3 Y- Y# u# }  @" d你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

& _% _' L1 y3 Z6 \9 W# E好的，谢谢了，

xiangpeng395

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

  P4 M9 U$ B* M: `可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！