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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   0 n+ d, r( b3 B+ ?+ B
2. -1 0 -1 1 0 0 0;
   ; F4 D. H; O- n) p/ K
3. 0 -1 1 0 1 1 0;5 T1 m9 O9 Q+ X- {3 \
   
4. 0 0 0 -1 -1 0 -1;8 O/ W% q3 O( a$ k; P3 a, d
   
5. 0 0 0 0 0 -1 1 ]
   6 u0 \! H) z) \- s* E5 T. h- A7 N
6. for i=1:5
   * N  J7 j$ n5 h, m; e
7. a=c;
   # u# Y; c) a8 I2 Z6 Y- o2 W# ^
8. b=a(i,:)+ W( D4 Z3 r/ n
   
9.      for j=1:4* D0 a) x8 Q  g/ ?6 I
   
10.         b(i,j);
    8 t, {. B: Q% X
11.             if b(i,j)~=08 n: A\" R3 B( p* z
    
12.                a(:,j)=[]) D) [# f. M0 \% \  R
    
13.             end0 g$ X* [9 E  h* {3 {
    
14.      end
    ( i% Z0 c+ D( q4 B& g6 G+ D9 f
15. a(i,:)=[]6 z$ A  A2 y4 L' w% x
    
16. end

复制代码

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

) V  F5 L8 A% x$ B: U6 P

madio

1. c=[1 1 0 0 0 0 0;; n8 a3 I1 _1 J\" I0 w! d4 o
   
2. -1 0 -1 1 0 0 0;$ g; [# q# l! D, [6 ^& m
   
3. 0 -1 1 0 1 1 0;; G8 h$ F4 O- C/ R$ j
   
4. 0 0 0 -1 -1 0 -1;
   ; F8 W( {6 C4 ^( a
5. 0 0 0 0 0 -1 1 ]
   * d; ~! f8 ~- e0 T
6. for i=1:5
   ; K$ M6 ^, b# q2 W7 C0 D
7. a=c;
   \" F. N. K! [, K5 R: _: d+ |
8. b(i,:)=a(i,:)
   8 F1 b4 c7 t: @$ O7 A) O
9.      for j=1:4% ~9 _6 m+ p0 u8 J
   
10.         b(i,j);
    ! F7 {. }$ J; a4 G6 S
11.             if b(i,j)~=0, a2 B7 a6 r% U- W7 Q# O
    
12.                a(:,j)=[]' g; [5 F' J2 l( D) A+ b, U: C  N+ t
    
13.             end
    / z3 N9 e* L9 ~3 P3 _1 }: x2 y
14.      end) \8 {3 o8 d/ C\" u: O
    
15. a(i,:)=[]
    / n; f8 ~2 A# _/ C7 D  n
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

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

madio

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

1. c=[1 1 0 0 0 0 0;/ {2 {% e3 z! q5 P4 V+ q
   
2. -1 0 -1 1 0 0 0;
   \" p2 C/ ~' g$ j\" v/ K, l* Y
3. 0 -1 1 0 1 1 0;; L2 G% b) O0 [( h- U0 [
   
4. 0 0 0 -1 -1 0 -1;
   2 o6 E+ O. a; G! D
5. 0 0 0 0 0 -1 1 ]6 Y1 L3 H& k* \/ v* Z
   
6. for i=1:55 \; _2 f& O0 Y5 e
   
7. a=c;
   # F: b6 n$ R5 f4 O. z3 W
8. b=[];! D; K( w( ^5 H
   
9. del_col=[];+ Y, i, b5 @5 y- ~3 \\" C9 d2 w
   
10. b(i,:)=a(i,:)
    $ w/ `! \+ W! y
11.      for j=1:7; e, o5 W* @6 j$ R- a
    
12.         b(i,j);+ i0 _% C8 n' a  z\" s# q  Q$ @) H
    
13.          if b(i,j)~=0
    , Z7 J' N, U* E- A\" v! P
14.            del_col=[del_col,j];9 p3 I: Q( @8 h$ c
    
15.          end
    $ v4 M* o' G5 ?# R
16.      end+ c9 Y2 m; j: K5 d7 T' m
    
17.      a(:,del_col)=[]
    3 n- t3 _* D% B\" Y9 \6 r/ j
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
8 z# m: }$ W8 C) b/ V$ I; D6 `你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

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

# ^. N6 k! N1 ^# \) K/ }3 n! D可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！