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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   4 S% Z- j3 d3 K: [' V/ l* p) c
2. -1 0 -1 1 0 0 0;
   6 O; M! p+ j9 Q7 c$ j  C\" z
3. 0 -1 1 0 1 1 0;
   . Y+ e9 i! \0 h# V. p+ O
4. 0 0 0 -1 -1 0 -1;
   . r9 G% X3 |4 C; n/ w\" p
5. 0 0 0 0 0 -1 1 ]: x& h! X\" ]4 D! L1 E
   
6. for i=1:5\" u0 N% R8 i5 i. x, }% A1 I
   
7. a=c;
   9 f9 }6 O# [\" h2 C& s
8. b=a(i,:)
   , ~: h; D7 }# k5 R
9.      for j=1:4
   / x' |* h+ n; \
10.         b(i,j);, i0 {7 O\" ^/ X3 Y\" n8 j' k
    
11.             if b(i,j)~=0
    , y  J5 r6 w4 K) P: Y
12.                a(:,j)=[], @1 a5 w  |% `. S* Z( q
    
13.             end
    4 k% q. m$ S2 P
14.      end
    8 p# a! F# r/ v- l  f
15. a(i,:)=[]
    4 q\" p; Z0 _8 t3 ^( ~
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;
   / i$ X1 k# T8 g/ U; o( k
2. -1 0 -1 1 0 0 0;
   % P0 Y* e5 s. \4 f3 ?' t- t: G\" U
3. 0 -1 1 0 1 1 0;2 m; Q2 \6 s' e9 C
   
4. 0 0 0 -1 -1 0 -1;
   7 Z# I7 b) x0 e0 o: o+ n' U( ?
5. 0 0 0 0 0 -1 1 ]
   1 r3 t7 ]3 l1 q, ~# b  ?8 y
6. for i=1:59 u* e\" e- h: X* @
   
7. a=c;
   0 I5 z/ T! R4 L! O0 M
8. b(i,:)=a(i,:)
   ' S+ J# E: ?; ^6 z
9.      for j=1:4) O+ W; C2 X! I
   
10.         b(i,j);
    * _1 e* F5 m/ A- [! F/ n
11.             if b(i,j)~=0
    8 d$ _+ X+ p% @1 {0 E: g
12.                a(:,j)=[]& `0 c9 z$ _) L# H0 E5 E
    
13.             end9 j9 z\" }5 [\" \# J( v- S; t8 M
    
14.      end
    * C( h7 `# K& z
15. a(i,:)=[]- E( V1 D+ o2 l) h
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

! G$ C7 A$ T( E7 w. u$ N不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;$ R2 Y2 ?; z7 J
   
2. -1 0 -1 1 0 0 0;( c/ W9 p8 `% s/ [, ~1 |+ V
   
3. 0 -1 1 0 1 1 0;! I6 Q4 A6 E4 _( O9 q8 Q
   
4. 0 0 0 -1 -1 0 -1;3 E! Y! k( b$ S& F/ }
   
5. 0 0 0 0 0 -1 1 ]+ w' U2 O\" W- h\" _9 g, j. H: k
   
6. for i=1:5
   3 U- N) R  k5 @: N$ q! l
7. a=c;
   : ]9 y: `2 Z9 x
8. b=[];, L1 I, r0 ]* B: S$ P
   
9. del_col=[];
   8 s' ~4 U- E' G, d$ q8 b
10. b(i,:)=a(i,:)
    ( v0 l8 S1 G2 R; R/ I; D( Q
11.      for j=1:7. @5 R4 A6 ~6 ]( c$ ]
    
12.         b(i,j);
    : Z+ ^1 J) C+ X2 R% W9 k
13.          if b(i,j)~=0
    , ~/ d5 K& y+ c. r, i
14.            del_col=[del_col,j];2 a, T: C* _8 D+ j$ N* Y* I. X) V
    
15.          end: r2 Q1 O1 D6 h; h
    
16.      end* R& W1 _! ?- q
    
17.      a(:,del_col)=[]# e& O) K$ E' s\" y
    
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
/ A5 g" X/ Z, d7 }* a你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

6 d7 r/ d9 v  U# j6 V好的，谢谢了，

xiangpeng395

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

. z9 ~( A# O: w可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！