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

xiangpeng395

1. c=[1 1 0 0 0 0 0;\" g% O1 h4 U# D: g$ |
   
2. -1 0 -1 1 0 0 0;2 l5 Q, d& K# R* Q( d
   
3. 0 -1 1 0 1 1 0;# Q/ ?3 X' J4 R* K\" f& y% a) f4 @
   
4. 0 0 0 -1 -1 0 -1;8 e; Z# g0 I+ o$ t3 {7 X
   
5. 0 0 0 0 0 -1 1 ]
   , o0 a5 _8 V3 c
6. for i=1:5% F, Z* C# k3 x2 {\" v  ?: G. L. J; `
   
7. a=c;% B6 h4 U' i8 V4 m- J
   
8. b=a(i,:)3 w) D/ u. j; \! g3 ^8 G& Q  V8 f
   
9.      for j=1:4: p( G8 R$ z% x, y9 v8 y' O\" G/ m% K7 ~
   
10.         b(i,j);( u9 A# j, H: t! D0 U, j
    
11.             if b(i,j)~=0
    1 `2 w% O3 `! L
12.                a(:,j)=[]
    9 [; v& D- g% V1 C! j/ @
13.             end/ j/ \4 z8 l. O% R( z* P: B
    
14.      end# y6 ]5 _+ {/ v; c6 r, D
    
15. a(i,:)=[]+ U9 o# |7 @; N( u
    
16. end

复制代码

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

3 D3 t$ F. [; q% A8 ?

madio

1. c=[1 1 0 0 0 0 0;: k$ Z8 h. G. R
   
2. -1 0 -1 1 0 0 0;& S& M\" }$ P. g
   
3. 0 -1 1 0 1 1 0;' ]1 R0 x) J8 o0 p/ @* `& [- f
   
4. 0 0 0 -1 -1 0 -1;/ ?) w7 q- T- A7 A! a
   
5. 0 0 0 0 0 -1 1 ]\" w5 D! w6 J7 J4 X6 j: E
   
6. for i=1:5; `: t0 z0 V, N, ~- B
   
7. a=c;
   . g# W' M' k8 J8 y1 k' T5 F- e
8. b(i,:)=a(i,:)
   5 d+ @+ }' D7 z5 ?: a$ E8 u\" b+ T; a+ Y
9.      for j=1:4
   + o4 s$ ], ^2 G( x' z; F
10.         b(i,j);
    \" C( Z: V- u# g, \+ J
11.             if b(i,j)~=05 }( G9 v$ i: K
    
12.                a(:,j)=[]
    1 e9 i, H) |( [  f- f, W
13.             end
    ' S3 K: M- r' ~. V
14.      end
    * w\" _+ X9 [* G& F\" @9 N$ R
15. a(i,:)=[]$ P, x\" h) Z6 P, c2 t\" k7 Z$ i
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

  A" F* N4 f  Z/ s不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;% A+ j3 m* B2 X
   
2. -1 0 -1 1 0 0 0;
   # R1 W! t: _4 U8 ~
3. 0 -1 1 0 1 1 0;: x) U4 X( D3 @+ c8 g\" x5 @* K
   
4. 0 0 0 -1 -1 0 -1;
   0 l- u6 g# S, J$ _( ?2 z+ Y3 o
5. 0 0 0 0 0 -1 1 ]
     E$ z: r; c8 q, r( K
6. for i=1:54 o0 ^* d$ k5 ^9 l9 U, G
   
7. a=c;
   * d( Y) G- W+ \, V0 Y0 X
8. b=[];
   * a$ {- j8 |& u# {$ F! @: X: G
9. del_col=[];
   ' ~0 r' Y9 s; K$ \; v# l2 ~. P
10. b(i,:)=a(i,:)
    ! n9 P  S- z0 @+ t4 ]0 Y
11.      for j=1:7
    , x( [4 W+ t2 O8 ~. J, W6 g
12.         b(i,j);, R: h' h, n+ Y4 o
    
13.          if b(i,j)~=0
    3 v+ A\" Z7 x* m
14.            del_col=[del_col,j];
    \" k( n8 f# K# f# c7 q* Y
15.          end' D8 `; i\" p- ?0 ]! K
    
16.      end
    5 q  ~+ }3 I9 c! C- S$ ]  ^/ I
17.      a(:,del_col)=[]6 D, ]! Y/ \9 _1 Z) ]
    
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
5 @- |( H6 t$ z/ r8 L; K& }6 j你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

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

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