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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   5 Q* d8 z8 F2 b. x1 c6 s
2. -1 0 -1 1 0 0 0;
   * J  b3 W1 B! L0 d! @3 p- e
3. 0 -1 1 0 1 1 0;' p. H/ O4 V4 V* P7 C
   
4. 0 0 0 -1 -1 0 -1;: g% _5 i/ E) I7 X( g1 I4 K
   
5. 0 0 0 0 0 -1 1 ]
   % F6 G! v4 w6 G
6. for i=1:5\" W# Z0 Z$ @2 u  A
   
7. a=c;
   7 r: D4 o3 S: ~; t) H: B$ z\" E# ^
8. b=a(i,:)4 [& {( F7 n( D8 k
   
9.      for j=1:4, F* v  v* U: L( o8 g9 ]3 J! C7 n. T' i# f
   
10.         b(i,j);
    - W6 P* q; S( \% {# X) b
11.             if b(i,j)~=07 X. q* Y5 M1 g# }5 L! \; F6 d
    
12.                a(:,j)=[]
    ; r- G) I8 z: d. ~7 h  }
13.             end
    - F* E- F! T% o4 a0 T( a& B
14.      end
    . Q+ B! z/ J4 i6 j' I
15. a(i,:)=[]; l1 U\" {: P4 w. Q8 V9 c
    
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;; E) z/ s( N: m* G0 N% k) z\" z& ^( n
   
2. -1 0 -1 1 0 0 0;\" X3 }\" o) b1 |# Q5 ^5 u
   
3. 0 -1 1 0 1 1 0;: a; _\" F2 \* q7 G6 _% M' R
   
4. 0 0 0 -1 -1 0 -1;
   : u6 k( L7 |2 T1 L+ \) w
5. 0 0 0 0 0 -1 1 ]: Z0 ~4 }* Y/ y7 ]/ k& g+ v- p
   
6. for i=1:5$ C9 F* F; P7 a$ `2 ]4 M( d& U
   
7. a=c;+ {8 q5 q4 y9 P6 _: f+ Q& [8 m
   
8. b(i,:)=a(i,:)
   4 O, K0 m2 {  A
9.      for j=1:4
   % t8 `1 O. X# {0 M1 {
10.         b(i,j);3 L/ G; e$ q0 T: P9 r) x% E
    
11.             if b(i,j)~=0( w- e9 `1 y, l/ S
    
12.                a(:,j)=[]
    8 u4 O8 h4 P- @2 z5 G3 t! m
13.             end
    1 R; H$ l3 u1 ^! W
14.      end& p$ \  f3 ^4 X; j; b/ s8 |$ K, }
    
15. a(i,:)=[]: |6 A4 B+ v% O0 V
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

2 C; |- K8 ~5 [不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;- W) ?5 ^6 {6 e' @: [3 F* e
   
2. -1 0 -1 1 0 0 0;
   ; \5 S. E* F- v) A1 V; E8 m
3. 0 -1 1 0 1 1 0;9 ^+ d: ~; x: ^8 b$ S  t% r- s' N7 ^3 R
   
4. 0 0 0 -1 -1 0 -1;
   # H( G8 l* r& k+ v  ]# Q$ Z
5. 0 0 0 0 0 -1 1 ]: K) A, x\" L8 P% \. k
   
6. for i=1:5% `, d- d( O8 \& {8 i, W
   
7. a=c;
   7 A6 o! _: g$ p, U7 [6 A
8. b=[];3 d; C+ T& O% M! }8 i/ a
   
9. del_col=[];- i$ G5 u7 n$ G4 a. u
   
10. b(i,:)=a(i,:)
    4 r% {& @9 R9 l4 p
11.      for j=1:7' v. V. I; T$ X# t3 R; u
    
12.         b(i,j);
    1 a6 D2 s1 `( o7 R
13.          if b(i,j)~=0
    / M: X9 H1 Y: y
14.            del_col=[del_col,j];
    + e8 x3 _$ m9 D+ ~# {
15.          end
    7 R5 ^4 P2 [& c4 j% Z/ S0 Z9 O
16.      end\" C; Y  `0 M& h1 J: B0 B9 e
    
17.      a(:,del_col)=[]' Q3 Z8 a+ K# ~& L
    
18. end

复制代码

xiangpeng395

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

2 C. Q, }: c% a( i好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
# d  x& O- b$ n( [1 Y/ z8 E你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

, n3 P( ~) b' q" f- P可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！