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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   5 i0 S3 ?7 g$ X: G# \# M
2. -1 0 -1 1 0 0 0;
   $ R, M: I- g9 f: p
3. 0 -1 1 0 1 1 0;
   \" ]3 E* V1 [# b0 R& X$ ?; q$ I
4. 0 0 0 -1 -1 0 -1;# g7 O5 {3 ^8 L$ o. S- \: i, P
   
5. 0 0 0 0 0 -1 1 ]9 e8 |( z! B0 a; B/ N
   
6. for i=1:5
   * l0 @8 c: l6 q% Q
7. a=c;; l( N# R1 k9 n2 q' C3 [
   
8. b=a(i,:)+ _4 {, t, Y, s\" t0 b6 g
   
9.      for j=1:4* X, q\" w0 `- t5 X. V5 w- ?; |: J
   
10.         b(i,j);
    : b. G: j\" ?) `- Y0 ]2 ?
11.             if b(i,j)~=0
    9 X1 G$ P& F; N
12.                a(:,j)=[]
    \" `5 c. \- Y( s6 k& A9 `' S
13.             end
    ' W) t2 k; d) N# G8 C- y* c
14.      end& z: T( a: x% S5 Q+ _
    
15. a(i,:)=[], r; G, x; q+ ]: a! X3 A
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;0 \- U5 v3 C0 ~3 @! l. M' e  A
   
2. -1 0 -1 1 0 0 0;
   2 m- J$ W2 K( K  g/ r$ ^3 R% k7 c
3. 0 -1 1 0 1 1 0;
   7 m8 @9 T  X7 U4 T
4. 0 0 0 -1 -1 0 -1;5 a0 v1 W7 V& S; ~- ^  G+ \
   
5. 0 0 0 0 0 -1 1 ]
   5 P/ N7 {1 G8 B5 ]: B7 N2 Q
6. for i=1:5  Y7 s9 F* M* W- I+ `) i& q
   
7. a=c;+ B$ j) R' j9 W8 \1 y& x; \
   
8. b(i,:)=a(i,:)
   - g/ |  _' }( w
9.      for j=1:4
   \" ?8 a( W/ {/ A* W
10.         b(i,j);0 l9 t7 |. b0 X# G
    
11.             if b(i,j)~=0: S0 i+ m\" a4 J, W3 E6 {
    
12.                a(:,j)=[]7 o) k% X7 X  i! h# p
    
13.             end
    5 t. p6 E& b3 W+ f
14.      end
    $ p6 v0 c- H( ~+ l: y# s
15. a(i,:)=[], Y6 G9 d! r2 @6 {7 W\" f  l. w& T
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
8 M  ^: r( ~2 {1 d改成这样的就可以运行了

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

madio

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

1. c=[1 1 0 0 0 0 0;
   ) L' A8 `+ U2 q, S6 r: `! X
2. -1 0 -1 1 0 0 0;3 T/ _8 k7 e7 m' o! ], e
   
3. 0 -1 1 0 1 1 0;7 u6 k7 |0 Y2 }5 h
   
4. 0 0 0 -1 -1 0 -1;
   - I9 R. [% L% c. A1 Q
5. 0 0 0 0 0 -1 1 ]
   # E4 [3 F7 m& I! M3 g2 \
6. for i=1:5
   + A2 F1 X* V  ^- r6 v+ j& e4 {0 t& e
7. a=c;/ p3 i9 Z& g' O7 G% y) p
   
8. b=[];7 @5 k8 T' l! t5 {$ w
   
9. del_col=[];. |6 y2 h0 }0 w
   
10. b(i,:)=a(i,:)
    8 w3 p& _) N  \/ m9 z' L2 I3 [
11.      for j=1:7, i% g& ^; ]+ W* v* h
    
12.         b(i,j);
    : K( H$ x& F) T4 p4 ?
13.          if b(i,j)~=05 M8 o3 V9 U/ E+ U- t) _. p
    
14.            del_col=[del_col,j];
    . r+ S) W9 _5 X4 ?0 K
15.          end' Y1 X% e8 I& @( U/ A( v  z
    
16.      end! O( F\" e! i+ W4 J( X
    
17.      a(:,del_col)=[]
    # o\" t7 {2 s& g$ K
18. end

复制代码

xiangpeng395

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

/ @+ u+ [3 a/ f6 `好的，谢谢了，

xiangpeng395

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

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