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

xiangpeng395

1. c=[1 1 0 0 0 0 0;# T& G# O1 U9 H5 U/ ~
   
2. -1 0 -1 1 0 0 0;
   + Y  N& l; A: Y$ K4 f\" g: L& R
3. 0 -1 1 0 1 1 0;
   + g2 x! C( K& n+ B/ o
4. 0 0 0 -1 -1 0 -1;: x- ]% r% T1 i& e
   
5. 0 0 0 0 0 -1 1 ]  h% w5 \: i5 `9 a8 X
   
6. for i=1:5
     c! D  H2 Q$ k, W; \1 e0 ?. e
7. a=c;/ x7 j4 C  |, f4 F
   
8. b=a(i,:)
   ' ?& ?2 s& T6 H$ ~- [& T
9.      for j=1:4* E0 h9 p3 i1 b+ x( {
   
10.         b(i,j);6 |6 [# K- x& u
    
11.             if b(i,j)~=0# O  k1 k7 c% h
    
12.                a(:,j)=[]
    8 o) V: C\" _* Z+ |+ x0 f' t
13.             end
    % _' m. h9 D+ }& K# i0 {8 x6 H
14.      end
    ) l1 l6 {4 Z, u) m7 U
15. a(i,:)=[]6 x6 c9 q4 h( R2 l# M2 i
    
16. end

复制代码

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

, b: p, q/ G* }' b6 R( I* y

madio

1. c=[1 1 0 0 0 0 0;\" \! z* x5 p* z; o
   
2. -1 0 -1 1 0 0 0;6 ?7 V' S! p4 r  B7 s9 g( [/ u
   
3. 0 -1 1 0 1 1 0;1 Y: t7 U\" U% P* l
   
4. 0 0 0 -1 -1 0 -1;
     k; W1 J& Z5 Y, j0 y* \3 j: J
5. 0 0 0 0 0 -1 1 ]
   , C5 ]& f4 U9 M! w7 @0 b\" f: m
6. for i=1:5\" H' V1 m\" v6 f: W9 f4 Y2 Y
   
7. a=c;
   ' p& j7 T; @7 G9 o\" [& _+ x
8. b(i,:)=a(i,:)
   + U( m. j: g/ r. E6 p2 w( r6 d
9.      for j=1:41 \3 S2 L! w/ L& a
   
10.         b(i,j);+ |% g+ j\" R( Z; N+ y: p8 ]: D( o
    
11.             if b(i,j)~=0
    & j/ S8 d) N% g/ y
12.                a(:,j)=[]+ [; Z/ N2 v3 J# f4 c1 ^, Y
    
13.             end
    # n4 Z  [2 F7 e1 g
14.      end
    9 `+ r. y; s& u3 m* Z. ]
15. a(i,:)=[]
    3 a  S  ]  j; g' H  S
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
, S) w9 l8 A  b: E2 ]1 m% W5 T) K改成这样的就可以运行了

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

madio

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

1. c=[1 1 0 0 0 0 0;4 t- }1 z, X1 a) b& h' Y
   
2. -1 0 -1 1 0 0 0;4 U, H3 f0 w4 l
   
3. 0 -1 1 0 1 1 0;: y# \! h9 o8 h) R1 l( y3 R* k! B0 y
   
4. 0 0 0 -1 -1 0 -1;7 u2 t8 B& H5 y
   
5. 0 0 0 0 0 -1 1 ]: ~) Y7 g# u/ W\" g& z+ y' W( O
   
6. for i=1:56 V( d( K2 {+ J4 f6 k- X2 Y$ e' W
   
7. a=c;
   + E3 K7 [' I# g( H- c8 z
8. b=[];
   $ ?& g  z: m9 H! G
9. del_col=[];
   & e# v0 K- N9 X
10. b(i,:)=a(i,:)- ^) ~. ^4 _! O$ w; C
    
11.      for j=1:7
    1 P2 U' ~& x! {1 z
12.         b(i,j);3 e5 d, g& @  U- b6 L
    
13.          if b(i,j)~=0
    7 i5 r* z; ^4 V4 A
14.            del_col=[del_col,j];
    1 K/ ^# s& R0 j8 B# s
15.          end9 J\" c5 W\" W3 F  k$ M, s. G
    
16.      end& g+ q0 T  o( H& E$ }8 r- L
    
17.      a(:,del_col)=[]; T\" |) R0 o' ^8 A& Y. \0 H& P, y. R
    
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

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

$ u7 j- c. o- Y. A  _/ N可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！