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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   8 j3 x  D1 R  O
2. -1 0 -1 1 0 0 0;4 V2 m2 D0 f- M8 ]
   
3. 0 -1 1 0 1 1 0;
   \" `\" V% h: k4 ~: d$ P+ }
4. 0 0 0 -1 -1 0 -1;0 R) v. q0 q$ E
   
5. 0 0 0 0 0 -1 1 ]
   # |( N0 }+ G! {# n5 G
6. for i=1:5
   ( j! Z% ^\" f* G' ~7 V
7. a=c;
     C' t* G: U# I9 _$ N2 J0 e3 j/ h, C
8. b=a(i,:)
   / W# P2 _; J0 I9 G
9.      for j=1:42 {, Y1 m. Y5 D7 ]6 F* e8 T3 i
   
10.         b(i,j);
      D$ }  F3 C/ J+ K% |, c
11.             if b(i,j)~=0; X0 X' o6 e0 u5 S$ {$ z+ e) ?
    
12.                a(:,j)=[]
    3 c& r8 `2 j: D- f4 C
13.             end
    ! ?: C5 _) a& t4 I9 F6 Q( I1 U, |
14.      end\" I/ B# j+ T  e1 L7 `; P( O9 g
    
15. a(i,:)=[]9 J5 E6 k# {) c. W5 c
    
16. end

复制代码

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

xiangpeng395

madio 发表于 2014-8-30 21:15
( {9 L; o- }5 w* _你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

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

xiangpeng395

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

好的，谢谢了，

madio

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

1. c=[1 1 0 0 0 0 0;
   ; S5 P7 |: w1 n- A7 B
2. -1 0 -1 1 0 0 0;
   ( ~1 z3 l1 F: {5 Z8 B6 k
3. 0 -1 1 0 1 1 0;
   % W% W9 D4 I6 q5 f
4. 0 0 0 -1 -1 0 -1;9 c1 K\" c! ?) e+ B
   
5. 0 0 0 0 0 -1 1 ]
   4 r9 I! U, w* S% A3 q' S\" p& l1 r
6. for i=1:5
   - a: e4 _( w- b4 @4 s. `
7. a=c;: m/ y9 ^8 r7 x- `% i
   
8. b=[];, D  X6 \: X: ?
   
9. del_col=[];
   , J8 ?6 f& ^5 U0 {3 a3 ?
10. b(i,:)=a(i,:)) p5 w. M\" Y: p6 I' L
    
11.      for j=1:7/ \: T1 O6 p2 Z& T7 o# L' G
    
12.         b(i,j);
    7 j4 P3 ]+ A) r' N( ]# n
13.          if b(i,j)~=0
    2 l' e) X* g' o8 U7 c. B( P) V$ S
14.            del_col=[del_col,j];
    6 h& _3 r) p' A0 |
15.          end7 v\" h/ o9 A\" @8 E1 ?\" W
    
16.      end
    & _: L( n\" [: C. M! e  {
17.      a(:,del_col)=[]( Q! Z3 G. X\" |7 ?  S1 u# k4 k
    
18. end

复制代码

xiangpeng395

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

/ W! b  b* \% c2 K8 Y  f不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

1. c=[1 1 0 0 0 0 0;
   9 T! t/ U5 T: N/ Z% d
2. -1 0 -1 1 0 0 0;
   ! N0 B5 g- ]4 o# `- S
3. 0 -1 1 0 1 1 0;! p( U$ R  @' d& g! S
   
4. 0 0 0 -1 -1 0 -1;' V/ o9 \; W1 V* Z
   
5. 0 0 0 0 0 -1 1 ]
   , V. c: Y5 r7 W: a
6. for i=1:5! d- e4 r( Y. X( a+ O: [- ^, _3 K
   
7. a=c;\" p. @- q7 j8 N% N! a) |) f# _
   
8. b(i,:)=a(i,:)+ P4 R. x  n\" I
   
9.      for j=1:4
   & @8 b# r. U  w0 G4 F/ v
10.         b(i,j);/ g) I! T$ E$ U% I8 O; Z6 \
    
11.             if b(i,j)~=0
    $ I. x2 b0 l' A% L! v: g
12.                a(:,j)=[]
    . b3 F3 d4 i/ ]% {. N5 M% k
13.             end9 F9 W  H\" ~2 M9 N; p
    
14.      end% [$ V) @( T5 c7 \6 j+ }! |6 K
    
15. a(i,:)=[]
    4 Q+ h$ U4 H0 _5 x; }
16. end

复制代码

改成这样的就可以运行了